
Thank you for your fast reply. I don't want to use noslip, Symmetry/Slip, outflow, neutral outflow and inlet boundary conditions. How can I describe a new boundary condition in the GUI?
I tried to describe a new boundary condition in the mfiles. Firstly, I set fea.phys.ns.bdr.sel as follows,
fea.phys.ns.bdr.sel = [ 1, 5, 1, 1, 2, 1, 4, 1, 1, 1 ];
and then for the second boundary I used symmetry/slip boundary condition, so that I get now for second boundary,
'bs_ns', 'n&.u = 0, n.(pI + mu(grad u+grad uT)).t = 0', 'Symmetry/slip', [], { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }, [], { 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip', 'solve_hook_bdrslip';
0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 0, 0, 0, 0, 0, 0, 0, 0, 0 } };
The boundary condition I want to describe is similar to the symmetry boundary condition. When I changed the " n.(pI + mu(grad u+grad uT)).t = 0", the solutions did not change. What can I do?
Thanks..
