[Jbsl]

Hello friends,
I'm currently working on my own solver of Eulers equations, using cell centered finite volume method with Maccormacks scheme. My domain is channel with circular bump (10% of height), and i have a problem with applying boundary condtition on bottom wall with that bump. I know i have to use formula vn = 0, however, i have no idea how to code it. On top wall it was no problem, since i just set the velocity in cell behind the wall to negative value of the velocity before the wall. So my question is how do i apply the vn = 0 on wall that is not straight? Hope my question makes sense. Sorry for my bad english. Any advice or reference will be welcomed.
Have a nice day.

[FMDenaro]

Originally Posted by Jbsl

Hello friends,
I'm currently working on my own solver of Eulers equations, using cell centered finite volume method with Maccormacks scheme. My domain is channel with circular bump (10% of height), and i have a problem with applying boundary condtition on bottom wall with that bump. I know i have to use formula vn = 0, however, i have no idea how to code it. On top wall it was no problem, since i just set the velocity in cell behind the wall to negative value of the velocity before the wall. So my question is how do i apply the vn = 0 on wall that is not straight? Hope my question makes sense. Sorry for my bad english. Any advice or reference will be welcomed.
Have a nice day.

Just as an example for a 2d problem, if you have a cartesian reference system you just set:

n.v= 0 = nx*u+ny*v -> u/v = -ny/nx

of course if you have also no-slip condition you must associate the tangential condition t.v = 0.