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

help on isotropic turbulence initial condition

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

LinkBack Thread Tools Display Modes
Old   May 19, 2003, 20:25
Default help on isotropic turbulence initial condition
Posts: n/a
Hello All:

I am doing LES on decaying of isotropic turbulence. I get a very strange problem on how to set the initial condition.

I am using Rogallo's code (1981). He defined two random components of the velocity field in the plane orthogonal to k.

Q1 = sqrt( E(k)/4pi k^2 ).exp(i.theta1).cos(alpha) Q2 = sqrt( E(k)/4pi k^2 ).exp(i.theta1).sin(alpha)

Two unit vectors are defined as

e1 = ( k * omega )/ ( |k*omega| ) e2 = ( k * e1 ) / ( |k*e1| )


u = Q1.e1 + Q2.e2

it is supposed to get an isotropic, homogeneous velocity field.

But I can't get it, the initial velocity field I get is not isotropic. The integral length scale L11, L22 and L33 are not equal. Also u_rms = v_rms ~ 0.6, but w_rms is very small ~ 1.e-6. Can't understand why that happens ......

Anyone can give me some suggestions? thank you very much!
  Reply With Quote

Old   May 20, 2003, 14:22
Default Re: help on isotropic turbulence initial condition
Posts: n/a
Do you mean Rogallo's method or did you are actually using the code written by Rogallo?
  Reply With Quote

Old   May 24, 2003, 02:37
Default Re: help on isotropic turbulence initial condition
Posts: n/a
Thank you! both actually. Made some modification of his original code but the method for the isotropical initial condition isnot changed I just found the reason. It becauses of the random function.

But I got another question. It's a free decaying case, u, v, w should have the same behaviour. But I find enrgy of u ( u^2 ) decease fast than v and w. So at the end of run, it becomes anisotropic ......

Anyone has idea why that happens?

I run the same code with another initial condition which got from DNS result. It runs quite well. But we want to try some new initial condition ( pulse initial condition ) for our study.

Thanks a lot!
  Reply With Quote

Old   May 24, 2003, 04:42
Default Re: help on isotropic turbulence initial condition
Posts: n/a

That is not so surprising. In general you will always find some anisotropy at the end of the run, but hopefully it is small. I suppose you have some 'seed' in your random number generator. Changing seed results in statistically independent initial velocity fields but with the same mean energy spectrum etc. If you do a couple of decay runs with different seeds, thus different initial velocity fields, and take the average of the runs then you should find approximately isotropic results.Thus isotropic results are only to be expected in the mean. The anisotropy is probably worse when the integral length scale is large in comparison with the box size.

  Reply With Quote

Old   May 24, 2003, 16:22
Default Re: help on isotropic turbulence initial condition
Posts: n/a
Thank you very much! Tom.

You solved my big problem! Yes, it is true. In my case, the integral length scale is large, about 1/3 of the box size. The reason is my energy peak is at first mode ( k = 1 ), while the integral scale is like uiuj(1)/q^2 in spectral space. Now I understand why other people set the energy peak at high mode like k > 8. I spent a lot of time on this problem because I only change the initial conditon why I can't get the right behavior.

Another consideration is for k = 1, there are only 8 points fall in this shell, it's really hard to get a real random phases. But for k > 8, there are thousands points I think.

I have an another question. If we want to watch the behavior of integral length sacles, I think I have 2 ways to do it now: 1. Generate several initial conditon by choosing different seeds, and average them to get a new isotropic initial condition. Then run with that initial condition once to get the result. 2. Run with diffferent initial conditions, but average the final results instead Which way is better? Thanks a lot! Your help will be highly appreciated!
  Reply With Quote

Old   May 26, 2003, 07:57
Default Re: help on isotropic turbulence initial condition
Posts: n/a

The second method is the best. The first method doesn't help you. You don't get a better initial field by averaging. But you can ofcourse average the integral length scales in the three directions.

By the way, if you want to study the integral length scales it is perhaps a good idea to look at the the two-point correlations. The two-point correlation should approach zero for long correlation lengths. If it doesn't (in other words the box is too small) then you can't say much about the behaviour of the integral length scales.

  Reply With Quote

Old   June 7, 2003, 21:17
Default Re: help on isotropic turbulence initial condition
Posts: n/a
thank you! Tom! spent 2 weeks work on that well, still not quite successful ...
  Reply With Quote

Old   March 12, 2010, 19:03
Senior Member
Join Date: Mar 2010
Posts: 557
Rep Power: 13
cfdnewbie is on a distinguished road
Hello dear colleagues,
I was just wondering whether this forum / topic is still active and if there's someone around who could help me out with Rogallo's procedure?
Thank you all in advance!
A CFDnewbie
cfdnewbie is offline   Reply With Quote

Old   April 28, 2011, 16:55
Question turbulence generator modeling in wind tunnel (orifice perforated plate OPP))
New Member
mehdi farahnak
Join Date: Apr 2011
Posts: 3
Rep Power: 8
farahnak.mehdi is on a distinguished road
hello every one.
i want to model a isotropic turbulent flow with openfoam.
my case is a turbulence generator in a wind tunnel.
i just know the velocity of air at entrance of the tunnel
i want to model it with simplefoam.
what boundary condition must i choose for P and U and ....?
please help me
farahnak.mehdi is offline   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
Compressible Nozzle Flow sebastian OpenFOAM Running, Solving & CFD 14 September 21, 2016 10:47
Initial condition of isotropic turbulence Ye Jianying Main CFD Forum 1 June 12, 2012 12:26
Problems with simulating TurbFOAM barath.ezhilan OpenFOAM 13 July 16, 2009 05:55
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09

All times are GMT -4. The time now is 10:22.