
[Sponsors] 
to how much turbulent kinetic energy is resolved from Openfoam LES results ? 

LinkBack  Thread Tools  Search this Thread  Display Modes 
January 26, 2013, 19:24 
to how much turbulent kinetic energy is resolved from Openfoam LES results ?

#1 
Senior Member
Join Date: Nov 2012
Posts: 168
Rep Power: 9 
Dear All,
After finishing a LES case using OPenfoam, I have the thinking to check how much turbulent kinetic energy (for instance 70% or 80%) is resolved based on this case's mesh and results. This can give me some feedback about the mesh or some numerical schemes. Does anybody know how to do this for the Openfoam LES results? best regards, H 

January 27, 2013, 18:41 

#2 
Senior Member

Hi H,
In my opinion, you can do this in two ways: 1) From the LES analysis you will have access to the mean and rms of the resolved velocity  you will want to check that the rms of velocity are converged. You can then calculate the resolved tke based on the rms of velocity. The unresolved can be calculated from the contribution of the SGS turbulence model. 2) Alternatively you could also look at the spectral distribution of the resolved turbulent structure using spatial FFT. This may have issues when coherent structures.
__________________
 Julien de Charentenay 

February 5, 2013, 07:04 

#3  
Senior Member
Join Date: Nov 2012
Posts: 168
Rep Power: 9 
Quote:
Thank you so much for your help. For the method 1 you mentioned, I need to obtain the mean and rms quantities on a runtime basis. I think I can these quantities from UPrime2Mean. But I found that even if I defined the fieldAverage1 in controlDict, I can not use these quantities in, e.g. createFields.H files. So how can I access to the these mean and variance during the simulation? Thank you. 

February 6, 2013, 18:10 

#4 
Senior Member

Hi H,
You may not need to do it during the simulation (i.e runtime), but just as a postprocessing. From a practical perspective, I am not able to help... I have done it in the past but with a different solver (starcd) at postprocessing stage only. I can't help with openFoam side... Julien
__________________
 Julien de Charentenay 

April 1, 2016, 06:11 

#5 
Senior Member
Freedom
Join Date: May 2014
Posts: 208
Rep Power: 9 
Dear Julien,
Could you give some more details or some references how to calculate the resolved kinetic energy and unresolved kinetic energy? The resolved kinetic energy is equal to k? K=0.5*sqrt(UUmean)? which variable can be characterize the unresolved kinetic energy? Could you give the formulation of them? Thank you in advance! Best, Xu 

April 1, 2016, 06:18 

#6 
Senior Member
Join Date: Jan 2013
Posts: 368
Rep Power: 10 

April 1, 2016, 08:59 

#7  
Senior Member
Freedom
Join Date: May 2014
Posts: 208
Rep Power: 9 
Quote:
I have read the paper by Pope, thank you very much. But could you give some more details how to calculate the unresolved part? e.g. Could you give the formulation of it? What do you mean " calculated from the contribution of the SGS turbulence model."?? Thank you again! Best, Xu 

April 1, 2016, 09:10 

#8 
Senior Member
Join Date: Jan 2013
Posts: 368
Rep Power: 10 
probably you can think about how the k is calculated for the LES model?


November 30, 2016, 10:03 

#9  
Senior Member
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 272
Rep Power: 14 
I know this might be a little late for you, but in an LES solution the unresolved energy is simply:
Quote:


November 30, 2016, 10:05 

#10 
Senior Member
Freedom
Join Date: May 2014
Posts: 208
Rep Power: 9 

December 31, 2016, 07:55 

#11 
New Member
Arone Brake
Join Date: Dec 2015
Posts: 1
Rep Power: 0 
How to calculate resolved and unresolved kinetic energy in paraview from openfoam data?


January 13, 2017, 17:47 

#12 
Member
Bashar
Join Date: Jul 2015
Posts: 74
Rep Power: 6 
Hi Wenxu
Can you please share how to practically solve it? i.e. the step to do that? if you have time. Thank you 

January 18, 2017, 12:01 

#13 
Senior Member
Freedom
Join Date: May 2014
Posts: 208
Rep Power: 9 
It is not very difficult. You can output this field, volScalarField& kToBeUsed = turbulence>k(); Then do what you want to deal with the kToBeUsed field. If you do not want to ruin the original field of k, use: volScalarField kToBeUsed = turbulence>k();


Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Turbulent kinetic energy and energy dissipation at boundary  mshotm2000  FLOW3D  4  December 12, 2016 09:44 
Contours of turbulent kinetic energy in LES using ANSYS FLUENT 13  khosrow  FLUENT  5  July 29, 2013 05:28 
Subgrid Turbulent Kinetic Energy in LES  saqure  ANSYS  3  September 7, 2012 16:44 
LES: mean velocity and turbulent kinetic energy  MET  FLUENT  8  December 8, 2006 05:08 
Why FVM for highRe flows?  Zhong Lei  Main CFD Forum  23  May 14, 1999 13:22 