CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

LES initial condition

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

Like Tree1Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   April 12, 2004, 17:37
Default LES initial condition
  #1
Vishnu K Vijayaraghavan
Guest
 
Posts: n/a
Hi all,

I am starting LES simulations in Fluent. I am simulating flow through a straight pipe with periodic conditions at the boundaries (Re = 100,000). What initial conditions have to be used? Is there any way to set the initial conditions at different nodes individually?

Any help will be greatly appreciated, Thanks, Vishnu
  Reply With Quote

Old   April 13, 2004, 06:27
Default Re: LES initial condition
  #2
Dmitry
Guest
 
Posts: n/a
hi, it's a good idea to start your simulation using a simple RANS turbulence model. In general, this step is optional and you can set any other initial conditions (IC), for example, from some inlet boundary. If you want to set IC at different nodes individually you can use patch possibility in FLUENT (solve->initialize->patch).

  Reply With Quote

Old   April 14, 2004, 12:07
Default Re: LES initial condition
  #3
Vishnu K Vijayaraghavan
Guest
 
Posts: n/a
Thanks for the reply!!

Is it necessary to initially perturb the flow, or will Fluent add random perturbations automatically?

Also, patch seems to be limited in the sense that only initial conditons in different zones can be specified (is this correct?). I am looking into DEFINE_INIT udf for specifying the initial condition at each node.

Thanks, Vishnu
  Reply With Quote

Old   April 15, 2004, 06:02
Default Re: LES initial condition
  #4
Dmitry
Guest
 
Posts: n/a
In my opinion it's not necessary to initially perturb the flow. In state of the art solutions some experts in this area set up perturbations on inlet boundaries for velocity and turbulent kinetic energy. But I did not perform such simulations and don't know is it really important.

Best regards.
  Reply With Quote

Old   July 30, 2009, 16:33
Default Imitate Periodic Boundary Without Periodic Boundary
  #5
New Member
 
Nicolas Castro
Join Date: Jul 2009
Posts: 4
Rep Power: 7
Nicolas_Castro is on a distinguished road
I have ran and gotten accurate results for flow in a pipe using LES with periodic boundary conditions. I now wish to change/remove my periodic boundaries to an inlet and an outlet. I would like to set the inlet instantaneous values (i.e. u, v, w, p) equal to the values at the outlet (or values at a plane some distance downstream) so that turbulence continues to feed itself. The reason I would like to do this is to be able to introduce particles into the pipe since fluent does not allow discrete particle modeling with periodic boundary conditions. Any help or suggestions would be much appreciated!
Nicolas_Castro is offline   Reply With Quote

Old   August 2, 2009, 18:30
Default
  #6
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Fluent has its own "turbulence" generation methods. Not the best ones but...
sbaffini is offline   Reply With Quote

Old   August 3, 2009, 09:19
Default
  #7
New Member
 
Nicolas Castro
Join Date: Jul 2009
Posts: 4
Rep Power: 7
Nicolas_Castro is on a distinguished road
I know Fluent has it's own turbulence generation that are not very good, which is why I would like to be able to set my inlet equal to instantaneous values I have already calculated that I know have the correct profile and turbulent characteristics. Again I have ran and attained a solution for turbulent pipe flow using periodic boundary conditions. All I want to do now is figure out what I need to come up with a UDF that will copy the values at a plane calculated some distance downstream and set the upstream inlet boundary equal to them so that turbulence is feeding itself without the need for periodic boundaries.
Nicolas_Castro is offline   Reply With Quote

Old   August 4, 2009, 09:18
Default
  #8
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Well...sorry for the stupid response. However what you want to set as a boundary condition is exactly the same as the periodic boundary condition so if it is not possible to use it there should be a proper reason (i guess...). You could use the previously calculated statistics (in the periodic case) in the turbulence generator and this should give you an inlet condition which is not that bad.


In any case, to force the boundary condition you mentioned you will need an UDF...i think that you could try to make a first cycle on the cells of your domain in a define_on_demand macro and store the cell index of those cells which are on the same plane you want to use as a mirror for the inlet.


After that, in the subroutine for the inlet velocity profile you should put an additional cycle on the previously memorized cells to find the cell whose coordinates (except the streamwise one) is the nearest to the point you want to apply the boundary condition. (i'm not sure if this kind of nesting is possible).

Maybe, if the cycle on the inlet cells does not change its order, after the first time step you could memorize the copied cells index in the correct order and you will not need any further cycle to find the nearest point but could access the right cell directly (again, the practical feasibility of this approach should be checked).
sbaffini is offline   Reply With Quote

Old   August 6, 2009, 09:24
Default
  #9
Member
 
vishwanath somashekar
Join Date: Apr 2009
Posts: 41
Rep Power: 8
vishwa is on a distinguished road
Hello All,
I am running a LES simulation in a confined impinging jet reactor. I am using the smagorinsky-lilly model with the default settings available in fluent.

From the experiments we are sure that there is a turbulence in the reactor. However, for the same conditions, in the simulation, it shows that it is totally laminar after I have ran for 1/3rd residence time.

Right now I have plug flow entering the two inlets of the reactor. The inlets are long enough for the flow to develop. Is there a better way to initialize the flow so that the turbulence kicks in.

Regards,
Vishwa
vishwa is offline   Reply With Quote

Old   August 6, 2009, 09:55
Default
  #10
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
There could be a lot of reasons. The first one is surely the model, you should use the dynamic smagorinsky model because the original smagorinsky model is usually found to be too dissipative (especially in the near wall region).

The second one is probably the plug flow inlet. It is usually very difficult to get the transition even with a very fine grid and a real inflow condition thus, even if the Fluent inflow condition is not such a good method, the plug flow is surely not going to help. It will no cost you that much to activate one of the two method available in Fluent.

Another one could depend on the numerics. The convective scheme should be the central one or, at least, the bounded central scheme (in this case you could also perform a simple laminar simulation according to the MILES approach).

Of course the grid should be fine enough to properly represents the turbulent scales. For wall bounded and channel/pipe flows the grid should have the following (somehow relaxed) spacings:

dx+ = 30 - 100 - streamwise
dz+ = 10 - 50 - spanwise
dy+ = 1 - near wall
dy+ = 30 - 50 away from the wall

but these values are not exhaustive and, of course, the jet deserves attention because it introduces its own scales as well as any other geometry dependent energy injection

Also, according to my experience with Fluent, the time step should be fine enough to have a courant number of 0.1 in the ideal case or courant=1-2 in the worst case. Of course you have to use the 2nd order time integration method. Also, to speed-up the simulation, i suggest the fractional step method for the pressure-velocity coupling.

these are general guidelines to be checked with your simulation results
sbaffini is offline   Reply With Quote

Old   August 20, 2009, 11:33
Default
  #11
Member
 
vishwanath somashekar
Join Date: Apr 2009
Posts: 41
Rep Power: 8
vishwa is on a distinguished road
Thanks very much for the reply sbaffini. could you please tell me what do these dx+, dy+ and dz+ that you have mentioned mean?

When using the FSM, I am guessing you need to have a really small time step and how will fluent in this case take care of the convergence ?

Thanks in advance
Vishwa
vishwa is offline   Reply With Quote

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

Old   August 27, 2012, 12:06
Default
  #13
New Member
 
hava
Join Date: Aug 2012
Posts: 11
Rep Power: 4
shmhava is on a distinguished road
Hello,

I hope someone would answer my post.
I'm trying for a very long time to simulate a simple fully developed flow in a square duct (4 walls). The flow is periodic in the streamwise direction. The Re number is 10504. I'm using the LES approach with smagorinsky constant equal to 0.1. I don't know why, but the u_rms, or more precisely, the location of the peak of u_rms is in a very wrong place. According to the literature, it should be somewhere between y+=15-17, but to me it's in y+=27.
The dynamic model helps to correct this result, but at the same time, my friction coefficient becomes too small, and the mean velocity doesn't behave according to the log law.
What am I doing wrong?
Please help me,
Hava
shmhava is offline   Reply With Quote

Old   August 27, 2012, 13:01
Default
  #14
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Dear hava,

are you dealing with the same test case of the following DNS?

http://www.cfmbyexample.com/resources/JFM_Duct.pdf

I have a personal opinion on this database (see the discussion here DNS of square duct at Re_tau = 600).

However, as a matter of fact, i performed some LES computation on this case at the same Re number (also for aspect ratios AR = 2 and 4) with Fluent and the Dynamic Model with a bounded central scheme, PRESTO! for pressure and the fractional step for pressure velocity-coupling. The outcome is the best i ever had with LES in Fluent and it is practically perfect, which is why i started questioning the database above.

To give you some help you should describe your computational setting with more details (number of grid points in the 3 directions, the law of their spatial distribution, maximum and minimum grid spacings in the 3 directions, time step, domain size, convective scheme etc. etc.)
sbaffini is offline   Reply With Quote

Old   August 27, 2012, 13:27
Default
  #15
New Member
 
hava
Join Date: Aug 2012
Posts: 11
Rep Power: 4
shmhava is on a distinguished road
Dear Sbaffini,

Thank you for your reply. I'm dealing exactly with the same test case. As for the information you asked: I'm using a condition of mass flow rate in the periodic condition and not pressure gradient and for that the fractional step is not available. Instead I'm using second order implicit for transient. I also use BCD. For pressure I'm using simplec and standard, because I read that PRESTO! can be dissipative. However, I used once the PRESTO! and it didn't change a thing.
My number of grid points are:120x120x64 which gives: dx+min=1.2, dx+max=13 and dz+=64. Time step is 0.0001s, which according to fluent, results in CFL_max~0.2-0.3. Domain size is HxHx6.4H.

What was your skin friction coefficient? it should be approximately 0.0075.
What was the ratio between the velocity in the center and the bulk velocity? Did you get the secondary flows in the corners?

Thanks a lot,
Hava
shmhava is offline   Reply With Quote

Old   August 27, 2012, 14:34
Default
  #16
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Dear hava,

my friction coefficient is:

cf = 2 * (u_tau / u_bulk)^2 = 2 * tau_wall / (rho * u_bulk^2) = 0.00713

and

u_centerline / u_bulk = 1.228

where:

u_bulk = 33.538762 m/s (volume average of the mean-streamwise-velocity over the whole domain)

tau_wall = 4.0049109 kg/m*s^2 (area average of the mean wall shear stress over the four walls)

rho = 1 kg/m^3

My other settings are:

Domain: 6.4H x H x H with H = 1m and 96 x 64 x 64 points (sin law in the wall normal directions)

Fixed pressure gradient: dp/dx = -16 kg/m^2*s^2

Dynamic Smagorinsky Model

mu = 1/300 kg/m*s

dt = 0.0001 s (second order implicit in time)

NITA - Fractional step with PRESTO!, Least-squares gradient reconstruction and bounded central scheme for the convective term.

Statistics gathered over 30 s (300k time steps)


There are some things i would say:

1) I made some tests in the past between the mass-flow/pressure-gradient option (channel flow at Re_tau = 590) and non sensible differences were present in the nondimensionalaized results. However, you should be aware that fixing the mass flow, all the problems will show up at the wall stress level. In contrast, fixing the pressure gradient, all the problems will show up at the bulk velocity level. Of course, friction coefficient problems will show up in both, but coming from different sources. My dimensional results with fixed mass flow were clearly different from the expected logarithmic profiles. Also, consider that the reference DNS is nondimensionalized by the LOCAL friction velocity.

2) I made some test with several numerical parameters (for this and the channel flow test case) and, while very few differences were found in this case (much more for the channel where i made them with the unbounded scheme), the PRESTO! was found as the best option (less dissipative) for LES. I'm not saying it is always the best option (several times i can't simply use it to maintain the solution stable), just that it is the best for this case with the bounded scheme.
I also tested the gradient options in some different tests (pipe flow as i remember) and, while still getting very few differences, i noticed the least-squares option as the best.

3) By comparison with your settings, i guess that the main difference is in the grid. Consider that one of the main concerns i have on the reference DNS is the obviously large grid step in the stream-wise direction (dx+ = 40 is LES!) together with the 5th order upwind scheme used therein. That's why i used their same grid step in the stream-wise direction.

4) As a matter of fact, i realized that the solution is strongly sensitive to the stram-wise resolution (we also worked on several spatial resolutions). Also, by a sin law wall normal point distribution, i got y+ lower than one at the wall with only 64 points while still having the maximum dy+ lower than 15.

I can't provide additional results or details as this is still a work in progress and the results are not definitive yet but, roughly speaking, we got a full superposition with mean and rms DNS results. We only got a light discrepancy in the secondary flows, mostly on the position, but they are certainly present. Sadly, the available DNS spectra are not very reliable in my opinion so we can't get a full meaningful comparison.

Hope this helps
Far likes this.

Last edited by sbaffini; August 27, 2012 at 14:38. Reason: added some details at the point (1)
sbaffini is offline   Reply With Quote

Old   August 27, 2012, 14:54
Default
  #17
New Member
 
hava
Join Date: Aug 2012
Posts: 11
Rep Power: 4
shmhava is on a distinguished road
Dear Paolo,

Is it possible to create sin law in the wall normal directions in gambit?
I used the successive ratio method, and with 60 nodes in the wall normal direction I received dx+max=60. How did you make it possible?

One more thing. What was your subgrid dynamic viscosity const? The maximum I got was 0.045.

Thank you,
Hava
shmhava is offline   Reply With Quote

Old   August 27, 2012, 15:36
Default
  #18
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Dear hava,

i made my own grid generator in matlab. You can find it attached at some post in here:

Writing mesh in FORTRAN or C

It is a simple one-block structured grid generator but it did the job. By now i can't explain much how it works, but it should be quite simple to understand.

For the dynamic constant i did not pay too much attention and i didn't take any average of it. My experience is that it usually covers the full spectrum of the possible values (0<Cd<0.0529) as it isn't so well behaved and requires clipping. However, if the grid is fine enough you can probably attain lower values (as in your case) but, if you check the values of the constant in time (during the simulation), i'm pretty sure at some point you will also get the maximum.
sbaffini is offline   Reply With Quote

Old   August 28, 2012, 10:58
Default
  #19
New Member
 
hava
Join Date: Aug 2012
Posts: 11
Rep Power: 4
shmhava is on a distinguished road
Dear Paolo,

I hope it's not too rude, but, could you send me your case file? I have no intention to use it for publication purposes or something like that. I just want to see what I'm doing wrong.

Thank you,
Hava
shmhava is offline   Reply With Quote

Old   August 28, 2012, 16:47
Default
  #20
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 518
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Dear hava,

of course you are not rude for asking, but i hope you understand why i can't give you my case files. This is not only because the work is still in progress but also because i am not the only one involved in this work; finally, and this is more trivial, i have no direct access to the data at the moment. However i'm here to give you any useful information and probably for the next week i will manage to also provide some graphical results.

Also, consider that the informations i gave you are pretty much all the ones we shared in our group at the moment the different cases had to be set up (different cases were set up by different people in different places at different moments).
sbaffini is offline   Reply With Quote

Reply

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
SLTS+rhoPisoFoam: what is rDeltaT??? nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13
Forces in OF15 richard OpenFOAM Running, Solving & CFD 166 September 17, 2011 03:24
MRFSimpleFOAM goes divergenced! renyun0511 OpenFOAM Running, Solving & CFD 0 November 19, 2009 03:11
ForcesCoeffs ronaldo OpenFOAM 4 September 14, 2009 07:11
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16


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