CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

The initial condition for turbulence channel flow

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   April 19, 2005, 02:48
Default The initial condition for turbulence channel flow
Ma Jianfeng
Posts: n/a
I'm a freshman in DNS and want to do the direct simulation of plane turbulence channel flow. The plane turbulence channel flow has been studied widely and have many available results, So I choose this case as my start.(the flow is 3-D,the periodic boundary condtions are imposed at the streamwise and spanwise, the no-slip boundary condition is imposed at the upper and lower wall) In the simulation, I found I can't given the proper initial condition. From the literature, some authors says they start their simulation from the laminar velocity profile with random noise. So I write the following code as the intial condition of computation.

ran=1.0e-8; srand(2); //create a seed of random-number for (i=2;i<L;i++) {

for (j=1;j<=M;j++)


for (k=2;k<N;k++)


u[i][j][k]=ran*rand()/double(RAND_MAX); v[i][j][k]=ran*rand()/double(RAND_MAX); w[i][j][k]=1.0-uy[j]*uy[j]+ran*rand()/double

(RAND_MAX); ///the main flow profile




} } In the computation, if I large the variable ran, then the pressure may not convergence. If I choose ran=1.0e-8 or even smaller, as the advance of time, these random noise may decrease step by step and the main flow w[i][j][k] didn't changeat all.

So who can help walk out this quandary?? Thanks a lot!!
  Reply With Quote

Old   April 20, 2005, 03:50
Default Re: The initial condition for turbulence channel f
Posts: n/a

It is better to introduce some large scale disturbances. Random noise with a very short wave length damps out rather quickly. Therefore, I just take some random sine/cosine motions with different wave lngths. Details do not matter so much and it does not have to be diverge free, at least if your code is working correctly. The initial fluctuations should be about 5 or 10% of the mean flow. If your code can't handle that then you should have a critical look at your algorithm.


  Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Forces in OF15 richard OpenFOAM Running, Solving & CFD 168 November 14, 2017 12:51
Low Mach number Compressible jet flow using LES ankgupta8um OpenFOAM Running, Solving & CFD 7 January 15, 2011 14:38
Problems with simulating TurbFOAM barath.ezhilan OpenFOAM 13 July 16, 2009 05:55
Computational time sunnysun OpenFOAM Running, Solving & CFD 5 March 16, 2009 04:32
MRFSimpleFoam amp cyclic patches david OpenFOAM Running, Solving & CFD 36 October 21, 2008 21:55

All times are GMT -4. The time now is 07:00.