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

Question about LES

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 16, 2013, 12:05
Default Question about LES
  #1
Senior Member
 
Phoevos
Join Date: Mar 2009
Posts: 104
Rep Power: 17
fivos is on a distinguished road
Hi to all,

I am currently trying to implement a LES model (simple Smagorinsky to begin with, later we shall see) in a CFD code. The code has capability of 3D and transient runs, of course and in its current state it is a pressure based RANS code (both compressible & incompressible).

First of all, it is my understanding that the algorithm should be more or less the same as it was in the former RANS form, with the main modification being the calculation of the subgrid turbulent viscosity; in other words it would be like a laminar solver with the viscosity increased by the subgrid viscosity. Is this true?

The main question I have is what is the treatment of the so called trace term of the SGS stress (the one mentioned as tkk in http://www.cfd-online.com/Wiki/Large...imulation_(LES)).

What exactly is this term (I mean what is the equation of this term)? In a work (http://croccolab.umd.edu/publication...FD00_vol13.pdf) it is mentioned that tkk is:

tkk = Ci * 2 * rho * Delta^2 * |S|^2

where Ci is a constant (0.09)
rho the density
Delta is the cell volume ^ (1/3)
and |S| = Sqrt(2*Sij:Sij), with Sij the element of the subgrid stress tensor.

Is the above correct?

I have seen also in several works and books (including Pope, Lesieur) that the tkk term is added to pressure. In others (including the Fluent manual, see also http://www.cfd-online.com/Forums/ope...pressible.html ), it is said that the term is omitted. What is your opinion on the treatment of this term?

Finally, do you have any book/reference to suggest for practical, application oriented LES implementation (I do not want another book pilling a ton pages on how LES is derived and when it comes to the implementation details of the LES equations it presents them in just 10 lines)?

Any comments are welcome
Thank you in advance
fivos is offline   Reply With Quote

Old   May 16, 2013, 13:10
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
- for incompressible flows, the trace is considered into the pressure term, in practice you compute an elliptic equation for a modified pressure field.
- for compressible flows you have to somehow model also the isotropic term.

Note that generally RANS codes are quite low accurate in the time-space discretizations. Often in LES is necessary to build an ad-hoc high accurate code (choosing between implicit or explicit filtering).
FMDenaro is offline   Reply With Quote

Old   May 16, 2013, 13:27
Default
  #3
Senior Member
 
Phoevos
Join Date: Mar 2009
Posts: 104
Rep Power: 17
fivos is on a distinguished road
I see, thank you for your quick answer.

I understand that both spatial and temporal resolutions (and the relevant discretization schemes) should be appropriate in order:
- to properly resolve the larger eddies
- to reduce numerical viscosity/diffusion, which may mask/smooth eddies

What is your opinion about the omission of the term from the pressure equation? Have you tried working with it / without it?

Currently I have implemented the LES model for incompressible flow (I will see later for compressible) and I am testing it on the flow around a square cylinder at Re~21000. I have not implemented the pressure term but up to now subgrid viscosity has reasonable values, similar to the ones reported in literature. However the simulation is still running and it will need quite some time....

The solution of the modified pressure field baffles me ... what will be used to calculate it? The velocity divergence (i.e. the continuity imbalance) ? The imbalance + other terms ?

Thanks
fivos is offline   Reply With Quote

Old   May 16, 2013, 13:40
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
As I wrote before, for incompressible flows (Div V = 0) the isotropic SGS stress is simply disregarded in the filtered equations. As a consequence, the filtered pressure term that you solve from the elliptic equation contains also this term (a simple discussion can be read in the book of Sagaut). This is the reason for your good results without taking it into account.
For compressible flows this is not possible because you have the state equation for the pressure. However, it is somehow debated if the dilatation model improve the results.
FMDenaro is offline   Reply With Quote

Old   May 16, 2013, 13:52
Default
  #5
Senior Member
 
Phoevos
Join Date: Mar 2009
Posts: 104
Rep Power: 17
fivos is on a distinguished road
So (correct me if I am wrong), you imply that in incompressible when I include the subgrid viscosity term, I instantly get instead of pressure, the pressure + the isotropic SGS stress term without having to do anything else (such as include other terms in the SIMPLE or PISO pressure correction equations), right?

Ok, in that case I think I got it,
Thank you very much for your help.
fivos is offline   Reply With Quote

Old   May 16, 2013, 14:12
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
yes, that's correct
FMDenaro is offline   Reply With Quote

Old   May 17, 2013, 05:37
Default
  #7
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,150
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
I just want to add that the reason for this is that for incompressible flows you solve for a scalar field (the pressure) whose gradient is used to correct the velocity and achieve a divergence free velocity field.

Now, in this framework, this scalar field is actually a lagrange multiplier more than just the pressure. So, it does not matter how many unknown gradient terms you have on your right hand side (pressure, SGS stress trace, gravity, etc.), in the exact moment you don't explicitly compute them they build up in the unique lagrange multiplier used to satisfy the divergence free condition for the velocity field.

I am working on a particular LES framework which also requires an explicit modeling term for the pressure itself. So i once tried, for an incompressible flow, to model this term and i also removed the trace free condition for the SGS stress. As a matter of fact i did not get any appreciable difference in the results. Of course, this mostly depends from the specific numerical method too and testing is always required.
sbaffini is offline   Reply With Quote

Old   May 17, 2013, 06:10
Default
  #8
Senior Member
 
Phoevos
Join Date: Mar 2009
Posts: 104
Rep Power: 17
fivos is on a distinguished road
Thank you for your input too Paolo.
fivos is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
LES filtering question martyn88 OpenFOAM Running, Solving & CFD 1 October 18, 2012 04:08
LES supersonic flow ziooscuro OpenFOAM Running, Solving & CFD 1 March 5, 2011 01:56
question regarding LES usker Siemens 1 June 29, 2008 19:33
LES beginner question Shuo Main CFD Forum 4 July 9, 2007 09:40
Basic LES question toyCFD Main CFD Forum 3 August 23, 2006 16:09


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