in my thesis i use stream function vorticity method to solve the flow in square cavity,the grid is stretched and the coordinates is generalized coordinates, the code run very well in laminar flow,but now ito write the code for turbulent flow using fortran 90,i need to know how to: 1-but boundary conditions for different variables 2- how to get y+ at any point during run 3-what is the limitation on time step

Hi Khairy,

there is a now out of date textbook by Spalding, Runchal, Gosman, and several others. It is however pretty hard to find, infact I searched several University libraries now and could not find it. If you look hard enough you will!

However, a more important question is why you are using stream-function vorticity methods? They have pretty much all been replaced by primitive variable formulations. This is because they do not extend well to 3-D and (as you have nodoubt noticed) are a challenge to formulate appropriate boundary conditions.

This is probably something you should discuss with your advisor?

Regards,

Bak_Flow

Hi Bak_Flow, thanks for reply stream function vorticity is prefered in the case of 2 D flow,incompressible from storage point of view the boundary condition of omga is the main problem but i handle this problem completely using image point method. finally are you write down a turbulence code using finite difference.

Regards,

Khairy

Hi Khairy,

there is a Fortran code in the text that I mentioned. It has been a number of years since I read it so do not recall the limitations ie, laminar, steady, any additional scalar equations solved, etc. It may be a useful start for you, however.

Good luck,

Bak_Flow