# help on isotropic turbulence initial condition

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

 May 19, 2003, 20:25 help on isotropic turbulence initial condition #1 Hall Guest   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| ) then 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!

 May 20, 2003, 14:22 Re: help on isotropic turbulence initial condition #2 RMG Guest   Posts: n/a Do you mean Rogallo's method or did you are actually using the code written by Rogallo?

 May 24, 2003, 02:37 Re: help on isotropic turbulence initial condition #3 hall Guest   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!

 May 24, 2003, 04:42 Re: help on isotropic turbulence initial condition #4 Tom Guest   Posts: n/a Hi, 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. Tom

 May 24, 2003, 16:22 Re: help on isotropic turbulence initial condition #5 hall Guest   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!

 May 26, 2003, 07:57 Re: help on isotropic turbulence initial condition #6 Tom Guest   Posts: n/a Hi, 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. Tom

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

 March 12, 2010, 19:03 #8 Senior Member   cfdnewbie Join Date: Mar 2010 Posts: 557 Rep Power: 13 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

 April 28, 2011, 16:55 turbulence generator modeling in wind tunnel (orifice perforated plate OPP)) #9 New Member   mehdi farahnak Join Date: Apr 2011 Posts: 3 Rep Power: 8 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post sebastian OpenFOAM Running, Solving & CFD 14 September 21, 2016 10:47 Ye Jianying Main CFD Forum 1 June 12, 2012 12:26 barath.ezhilan OpenFOAM 13 July 16, 2009 05:55 carsten OpenFOAM Bugs 11 September 12, 2008 11:16 lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09

All times are GMT -4. The time now is 20:46.