Initial conditionfor turbulent channel flow in LES
hi, i am solving fully developed plane channel flow by LES. I am using periodic boundary condition in steramwise(x-direction),and spanwise(z-direction).A no slip b.c for top and bottom wall(y-direction).
For initial condition i am taking a fully developed mean profile with random perturbation for the u-component (streamwise)of velocity. and v=0,w=0. like as follows...
For the whole domain---- u(x,y,z)=(1-y**8)+0.005*(RAND(x,y,z)-0.5 v(x,y,z)=0.0 w(x,y,z)=0.0
And in fortran programming language---
do 55 i=1,iim
do 55 j=1,jim
do 55 k=1,kim
i,j,k--is spatial grid in three directions
iim=maximum grid no in x-direction
jim=maximum grid no in y-direction
kim=maximum grid no in z-direction
RAND(IND)------is a inbuilt function in fortran. As, the value of IND , change in the loop--then RAND(IND), will give a value , between--0.0-1.0
my question is---------
1. Is my approach correct.
2.Do i need random values for v and w velocity component.
3.Do i need to change the random function.
Please help me. A warm thanks in advance..
|All times are GMT -4. The time now is 10:02.|