|November 30, 2007, 13:04||
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..
|Thread||Thread Starter||Forum||Replies||Last Post|
|Extrusion with OpenFoam problem No. Iterations 0||Lord Kelvin||OpenFOAM Running, Solving & CFD||8||March 28, 2016 11:08|
|Floating point exception error||Alan||OpenFOAM Running, Solving & CFD||10||April 6, 2012 14:02|
|Orifice Plate with a fully developed flow - Problems with convergence||jonmec||OpenFOAM Running, Solving & CFD||3||July 28, 2011 05:24|
|MRFSimpleFOAM goes divergenced!||renyun0511||OpenFOAM Running, Solving & CFD||0||November 19, 2009 03:11|
|Negative value of k causing simulation to stop||velan||OpenFOAM Running, Solving & CFD||1||October 17, 2008 05:36|