Hi,
I'd like to define a free slip boundary condition at the wall. I dont want to specify a symmetry condition instead, because I 'd like to use the wall to do some calculations with radiation. How can I achieve it?
Best regards
Thomas
Free slip conditions at solid walls
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
-
David Monfort
Re: Free slip conditions at solid walls
Hi,
Sorry the delay of my answer Thomas... we were kind of busy lately ;)
First of all, and as you wrote, radiative transfer modelling needs that the walls are defined as, well, walls! So you'll have to keep defining the boundary face type as a wall i.e. itypfb(ifac,iphas) = iparoi. Now, in order to define the sleep-condition for the velocities, you will define the boundary condition type f or the variable (here the velocity components) as a symmetry i.e. icodcl(ifac,iu(iphas)) = 4 (and similarly for iv(iphas) and iw(iphas)). Please note that in the previous code, iphas is always equal to 1.
To achieve this, you need to use the Fortran subroutine usclim.f90 for advanced boundary conditions setup. Some explanation is given at the beginning of the subroutine (with comments in English now in the new RC1 version...)
Hope this helps, and works!
Let me know otherwise.
David
Sorry the delay of my answer Thomas... we were kind of busy lately ;)
First of all, and as you wrote, radiative transfer modelling needs that the walls are defined as, well, walls! So you'll have to keep defining the boundary face type as a wall i.e. itypfb(ifac,iphas) = iparoi. Now, in order to define the sleep-condition for the velocities, you will define the boundary condition type f or the variable (here the velocity components) as a symmetry i.e. icodcl(ifac,iu(iphas)) = 4 (and similarly for iv(iphas) and iw(iphas)). Please note that in the previous code, iphas is always equal to 1.
To achieve this, you need to use the Fortran subroutine usclim.f90 for advanced boundary conditions setup. Some explanation is given at the beginning of the subroutine (with comments in English now in the new RC1 version...)
Hope this helps, and works!
Let me know otherwise.
David