CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   to how much turbulent kinetic energy is resolved from Openfoam LES results ? (https://www.cfd-online.com/Forums/openfoam-solving/112358-how-much-turbulent-kinetic-energy-resolved-openfoam-les-results.html)

 hz283 January 26, 2013 20:24

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

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:confused:?

best regards,
H

 julien.decharentenay January 27, 2013 19:41

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.

 hz283 February 5, 2013 08:04

Quote:
 Originally Posted by julien.decharentenay (Post 404401) 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.

 julien.decharentenay February 6, 2013 19:10

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

 wenxu April 1, 2016 06:11

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?

Best,
Xu

 openfoammaofnepo April 1, 2016 06:18

Pls see below:

https://tcg.mae.cornell.edu/pubs/Pope_NJP_04.pdf

 wenxu April 1, 2016 08:59

Quote:
 Originally Posted by julien.decharentenay (Post 404401) 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

 openfoammaofnepo April 1, 2016 09:10

probably you can think about how the k is calculated for the LES model?

 anishtain4 November 30, 2016 11:03

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());

 wenxu November 30, 2016 11:05

Quote:
 Originally Posted by anishtain4 (Post 627605) 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

 arone December 31, 2016 08:55

How to calculate resolved and unresolved kinetic energy in paraview from openfoam data?

 Bashar January 13, 2017 18:47

Hi Wenxu

Can you please share how to practically solve it? i.e. the step to do that? if you have time.
Thank you

 wenxu January 18, 2017 13:01

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();

 All times are GMT -4. The time now is 22:58.