# to how much turbulent kinetic energy is resolved from Openfoam LES results ?

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

 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: 171 Rep Power: 13 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 vivek05 and Gang Wang like this.

 January 27, 2013, 18:41 #2 Senior Member   Julien de Charentenay Join Date: Jun 2009 Location: Australia Posts: 231 Rep Power: 18 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: 171
Rep Power: 13
Quote:
 Originally Posted by julien.decharentenay 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.
Hi Julien,

Thank you so much for your help. For the method 1 you mentioned, I need to obtain the mean and rms quantities on a run-time 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   Julien de Charentenay Join Date: Jun 2009 Location: Australia Posts: 231 Rep Power: 18 Hi H, You may not need to do it during the simulation (i.e runtime), but just as a post-processing. From a practical perspective, I am not able to help... I have done it in the past but with a different solver (star-cd) at post-processing 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: 209 Rep Power: 13 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(U-Umean)? 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: 372 Rep Power: 14 Pls see below: https://tcg.mae.cornell.edu/pubs/Pope_NJP_04.pdf

April 1, 2016, 08:59
#7
Senior Member

Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13
Quote:
 Originally Posted by julien.decharentenay 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. .

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: 372 Rep Power: 14 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
Posts: 273
Rep Power: 18
I know this might be a little late for you, but in an LES solution the unresolved energy is simply:

Quote:
 autoPtr model ( incompressible::turbulenceModel::New(U, phi, laminarTransport) ); volScalarField k(model->k());

November 30, 2016, 10:05
#10
Senior Member

Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13
Quote:
 Originally Posted by anishtain4 I know this might be a little late for you, but in an LES solution the unresolved energy is simply:
Thank you! Yes, I have solved it!

Best regards,
XU

 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? Bashar and Gang Wang like this.

 January 13, 2017, 17:47 #12 Member   Bashar Join Date: Jul 2015 Posts: 74 Rep Power: 11 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: 209 Rep Power: 13 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(); Far, Bashar, SHUBHAM9595 and 2 others like this.