# Reg. LES in OpenFOAM

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

 May 28, 2010, 11:13 #21 Senior Member   Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 13 The "average" value that you specify is not what is applied to the boundary. If you specify an average value, then the instantaneous mapped data is scaled such that the average of the mapped data matches the specified average. If you want to use mapping for a developing boundary layer, then you need something different. I suggest you try and contact Gavin Tabor at Exeter University. They published some work on this type of stuff.

 July 28, 2010, 14:03 SGS stress #22 Member   Kim Yusik Join Date: Dec 2009 Posts: 39 Rep Power: 8 Hi, Eugene. I am simulating the plane channel flow using LES and I look into your PhD thesis as a reference. It helps alot, thanks. BTW, I also want to see SGS stress on channel as you presented (in your case SGS shrear stress, Sigma_xy). So I am wondering if there is any option to calculate SGS stress. I think I need to calculate SGS stress and U simultaneously because SGS stress = nu_SGS * dU_i/dx_j (U_i is resolved velocity vector, not mean velocity). Correct me If I am wrong. Thanks in advance. Yusik

 August 27, 2010, 18:10 #23 Senior Member   Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 13 Unfortunately, to calculate the mean sgs stress, will require some coding. You need to add something like this to your solver: Before the time loop in pisoFoam, e.g. in createFields.H add: "tmp R(turbulence->R());" for the oodles solver replace with sgsModel->() . Also add the following inside the time loop: R = turbulence->R(); Now, use the fieldAverage function object to get RMean. Hope this helps.

September 8, 2010, 03:26
#24
Member

Kim Yusik
Join Date: Dec 2009
Posts: 39
Rep Power: 8
Quote:
 Originally Posted by eugene Unfortunately, to calculate the mean sgs stress, will require some coding. You need to add something like this to your solver: Before the time loop in pisoFoam, e.g. in createFields.H add: "tmp R(turbulence->R());" for the oodles solver replace with sgsModel->() . Also add the following inside the time loop: R = turbulence->R(); Now, use the fieldAverage function object to get RMean. Hope this helps.

Thanks, Eugene.

To generate R field, I added this code in 'createField.H' :
Info<< "Reading field U\n" << endl;
volSymmTensorField R
(
IOobject
(
"R",
runTime.timeName(),
mesh,
IOobject::AUTO_WRITE
),
mesh
);

rather than
tmp<volSymmTensorField> R(turbulence->R());

because if I use only one line above, I had this error message ' turbulence was not declared in this scope'.

The others are fine. Thanks again.
Yusik

 September 8, 2010, 09:37 #25 Senior Member     Daniel WEI (老魏) Join Date: Mar 2009 Location: Beijing, China Posts: 689 Blog Entries: 9 Rep Power: 13 He was not wrong, I guess you have misunderstood him, "turbulence" in his post is more a general type than a specific variable. Bye __________________ ~ Daniel WEI ------------- Boeing Research & Technology - China Beijing, China Email

 July 24, 2012, 17:25 #26 Senior Member   ehsan Join Date: Mar 2009 Posts: 106 Rep Power: 9 I had a basic question on LES. In OF, we use the same equations for momentum while we selected every type of turbulence models including LES, K-e and so on. So, where is the effect of volume filtering in LES in momentum equation? Thanks

 July 24, 2012, 17:31 #27 Senior Member     Daniel WEI (老魏) Join Date: Mar 2009 Location: Beijing, China Posts: 689 Blog Entries: 9 Rep Power: 13 Generally speaking, filter operation is done "implicitly", which is a common practice in many finite volume code. And you are right, the idea of "implicit" filtering or averaging, makes the gap between RANS and LES no longer big from the programming point of view, and then hybrid methods arise! Go hybrid! __________________ ~ Daniel WEI ------------- Boeing Research & Technology - China Beijing, China Email

 July 24, 2012, 17:39 #28 Senior Member   ehsan Join Date: Mar 2009 Posts: 106 Rep Power: 9 Hi Daniel O.K, But does it mean all LES models in Openfoam are implicit? When we select some models which permit us selecting filter (like Laplace), I think the LES in explicit in these cases, am I wrong? Ehsan

 July 24, 2012, 17:45 #29 Senior Member   ehsan Join Date: Mar 2009 Posts: 106 Rep Power: 9 And I had another key question: If our Momentum equation is the same for RANS and LES, and all the effects of small scales eddies are given by: B=-2v_sgs*D, the final momentum equation in LES reads the same as any other turbulence modelling: v_effective=v+v_sgs (in RANS it is v+v_turb), assume we use OneEqEddy model and solve an equation for k to get v_sgs. This OneEqEddy is typically close to the one in RANS. So, what is the different LES and RANS finally? Why LES gives much better results? Only because this v_sgs is limited to sub grid sclae unlike RANS where v_t is for whole domain? Ehsan

 July 24, 2012, 18:01 #30 Senior Member     Daniel WEI (老魏) Join Date: Mar 2009 Location: Beijing, China Posts: 689 Blog Entries: 9 Rep Power: 13 Are LES and RANS big different? In short: Yes and No. 1. The stress scaling with delta or scaling with wall distance makes a huge difference, the former we call LES, the latter we call RANS. 2. Of cause, there are more differences in numerics so it is not easy to convert a RANS-based code to LES, but if you just consider the turbulence modeling part, it is not that difficult, as you already noticed, their formula are quite similar 3. Both ways, both formulations are correct, so once you know how to use them, their answers should be bother correct and accurate, so saying LES is always better is not correct. 4. Since LES only try to model a small portion of energy in the spectrum, and these small eddies tends to be isotropic away from the wall, so the model accuracy is not so sensitive as a rans model, so even with a smagrinsky model, you can have a pretty good results, no needs always go dynamics, however, as you may experience, building a powerful RANS model is such a great challenge because it is not easy to find an universal one, but to be sure, we still need a new and powerful generation of RANS model. __________________ ~ Daniel WEI ------------- Boeing Research & Technology - China Beijing, China Email

 July 24, 2012, 18:12 #31 Senior Member   ehsan Join Date: Mar 2009 Posts: 106 Rep Power: 9 Thanks a lot, could you please give a clue on my former question, explicit LES in OpenFOAM as well. Regards

 July 24, 2012, 18:21 #32 Senior Member     Daniel WEI (老魏) Join Date: Mar 2009 Location: Beijing, China Posts: 689 Blog Entries: 9 Rep Power: 13 That I dont know very well, hope someone else could shed some lights on it and share some comparison results! My understanding is that through some mathematical operations, the "effect" of of filter is actually only working on the SGS stress term, that is the viscosity formulation in terms of EVM theory. At least for some filters. __________________ ~ Daniel WEI ------------- Boeing Research & Technology - China Beijing, China Email

 July 25, 2012, 06:00 #33 Senior Member   Vesselin Krastev Join Date: Jan 2010 Location: University of Tor Vergata, Rome Posts: 364 Rep Power: 11 To my knowledge an explicit LES model means that the sgs viscosity is modeled explicitly in some way (e. g. with an algebraic expression which uses flow variables, such as the Smagorinsky formulation). On the other side an implicit LES does not have a separately modeled sgs viscosity, but it simply introduces the necessary viscosity level into the flow domain in form of numerical viscosity generated by the numerical schemes. So, generally speaking OpenFOAM treats LES modeling in an explicit way. Regards V.

 July 25, 2012, 08:34 #34 Senior Member   ehsan Join Date: Mar 2009 Posts: 106 Rep Power: 9 If you see in the implementations of LES models in OF, some of theme does not require any filter (like OneEqEddy) but some other categories require it. I think that implicit models does not need filter while explicit ones need. Regards

 July 25, 2012, 08:56 #35 Senior Member   Vesselin Krastev Join Date: Jan 2010 Location: University of Tor Vergata, Rome Posts: 364 Rep Power: 11 In the OneEqEddy model there is still an explicitly calculated sgs viscosity, which depends on the filtering spatial cutoff value "delta", so to me this cannot be defined as an implicit LES model. However, if someone else wants to enter the discussion it would be interesting to further clarify this point. Regards V.

July 25, 2012, 09:05
#36
Senior Member

ehsan
Join Date: Mar 2009
Posts: 106
Rep Power: 9
1- Please see this nice paper which comes from Furbey,

Numerical Simulation of an Oscillating Cylinder Using Large Eddy Simulation and Implicit Large Eddy Simulation

DOI: 10.1115/1.4005766

Journal of Fluids Engineering, 2012, MARCH 2012, Vol. 134

See Eq. 2 and 5 about implicit and explicit LES. This work is done with OpenFOAM.

2- And as I see the 2008 paper (attached some pages), oneEqEddy is explicit, as you told. So where are the implicit models that Bensow used?

Attached Files
 Pages from bensow 2012.pdf (60.7 KB, 89 views) Pages from Bensow 2008.pdf (94.6 KB, 68 views)

 July 25, 2012, 09:35 #37 Senior Member   Vesselin Krastev Join Date: Jan 2010 Location: University of Tor Vergata, Rome Posts: 364 Rep Power: 11 Honestly I don't know how the implicit LES modeling has been implemented in these papers (I don't have access to the full copies), but remember that OpenFOAM is an open code, so the fact that apparently there isn't a "standardized" implicit LES modeling approach in the official releases does not mean that such modeling cannot be done (probably you should ask one of the authors of the papers for more information about it). Anyway, I think that an interesting passage extracted from the second of your attachments is the following: "The discretization scheme used is a central differencing scheme with different TVD limiters applied to the convective term to ensure stability and introduce the necessary dissipation for the implicit modeling." Regards V.

July 25, 2012, 10:04
#38
Senior Member

Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
I was saying "implicit filtering", not Implicit LES (ILES). ILES is a different matter.

Quote:
 "So, generally speaking OpenFOAM treats LES modeling in an explicit way." -vkrastev
I Agree!
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email

August 29, 2012, 06:37
solvers for LES in steady flow
#39
Member

Paula
Join Date: Aug 2012
Posts: 30
Rep Power: 6
Quote:
 Originally Posted by eugene For LES I would suggest pimpleFoam with the appropriate turbulenceProproperties (LES) settings and boundary conditions. pimpleFoam is like pisoFoam but with incomplete convergence of the corrector steps. For external aero LES just use the SpalartAllmaras DES model as SGS model.
Hi Eugene,

You told about the solver pisoFoam to solve a flow for LES, but Iīve read itīs for transient flow... do you know which solver can be used to solve a steady one?

Paula

February 26, 2014, 21:07
#40
Member

Peter
Join Date: Nov 2011
Posts: 45
Rep Power: 6
Quote:
 Originally Posted by alberto ... Depending on the information you have from experiments, you can, at least roughly, estimate the residual kinetic energy at the inlet, and use that value....
Hi, Alberto!

Could you please give me any suggestion on how should I set up k at wall boundary in my LES simulation? The y+ is smaller than 2 in my case.

Best regards,
Peter

 Tags les, pisofoam

 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 cedric_duprat OpenFOAM Running, Solving & CFD 154 September 21, 2016 15:31 fs82 OpenFOAM Bugs 21 November 16, 2009 09:15 fs82 OpenFOAM 6 October 13, 2009 09:58 pete Site News & Announcements 0 June 29, 2009 05:56 grtabor OpenFOAM Running, Solving & CFD 0 March 25, 2008 05:36

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