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

 May 16, 2013, 11:05 Question about LES #1 Senior Member   Phoevos Join Date: Mar 2009 Posts: 104 Rep Power: 10 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

 May 16, 2013, 12:10 #2 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 3,751 Rep Power: 41 - 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).

 May 16, 2013, 12:27 #3 Senior Member   Phoevos Join Date: Mar 2009 Posts: 104 Rep Power: 10 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

 May 16, 2013, 12:40 #4 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 3,751 Rep Power: 41 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.

 May 16, 2013, 12:52 #5 Senior Member   Phoevos Join Date: Mar 2009 Posts: 104 Rep Power: 10 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.

 May 16, 2013, 13:12 #6 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 3,751 Rep Power: 41 yes, that's correct

 May 17, 2013, 04:37 #7 Senior Member     Paolo Lampitella Join Date: Mar 2009 Location: Italy Posts: 815 Blog Entries: 17 Rep Power: 23 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.

 May 17, 2013, 05:10 #8 Senior Member   Phoevos Join Date: Mar 2009 Posts: 104 Rep Power: 10 Thank you for your input too Paolo.

 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 martyn88 OpenFOAM Running, Solving & CFD 1 October 18, 2012 03:08 ziooscuro OpenFOAM Running, Solving & CFD 1 March 5, 2011 01:56 usker Siemens 1 June 29, 2008 18:33 Shuo Main CFD Forum 4 July 9, 2007 08:40 toyCFD Main CFD Forum 3 August 23, 2006 15:09

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