CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

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

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

Like Tree11Likes
  • 2 Post By hz283
  • 2 Post By wenxu
  • 2 Post By arone
  • 5 Post By wenxu

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 26, 2013, 20:24
Question to how much turbulent kinetic energy is resolved from Openfoam LES results ?
  #1
Senior Member
 
Join Date: Nov 2012
Posts: 171
Rep Power: 13
hz283 is on a distinguished road
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.
hz283 is offline   Reply With Quote

Old   January 27, 2013, 19:41
Default
  #2
Senior Member
 
Julien de Charentenay
Join Date: Jun 2009
Location: Australia
Posts: 231
Rep Power: 17
julien.decharentenay is on a distinguished road
Send a message via Skype™ to 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.
__________________
---
Julien de Charentenay
julien.decharentenay is offline   Reply With Quote

Old   February 5, 2013, 08:04
Default
  #3
Senior Member
 
Join Date: Nov 2012
Posts: 171
Rep Power: 13
hz283 is on a distinguished road
Quote:
Originally Posted by julien.decharentenay View Post
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.
hz283 is offline   Reply With Quote

Old   February 6, 2013, 19:10
Default
  #4
Senior Member
 
Julien de Charentenay
Join Date: Jun 2009
Location: Australia
Posts: 231
Rep Power: 17
julien.decharentenay is on a distinguished road
Send a message via Skype™ to julien.decharentenay
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
julien.decharentenay is offline   Reply With Quote

Old   April 1, 2016, 07:11
Default
  #5
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 12
wenxu is on a distinguished road
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
wenxu is offline   Reply With Quote

Old   April 1, 2016, 07:18
Default
  #6
Senior Member
 
Join Date: Jan 2013
Posts: 372
Rep Power: 14
openfoammaofnepo is on a distinguished road
Pls see below:

https://tcg.mae.cornell.edu/pubs/Pope_NJP_04.pdf
openfoammaofnepo is offline   Reply With Quote

Old   April 1, 2016, 09:59
Default
  #7
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 12
wenxu is on a distinguished road
Quote:
Originally Posted by julien.decharentenay View Post
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
wenxu is offline   Reply With Quote

Old   April 1, 2016, 10:10
Default
  #8
Senior Member
 
Join Date: Jan 2013
Posts: 372
Rep Power: 14
openfoammaofnepo is on a distinguished road
probably you can think about how the k is calculated for the LES model?
openfoammaofnepo is offline   Reply With Quote

Old   November 30, 2016, 11:03
Default
  #9
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 273
Rep Power: 17
anishtain4 is on a distinguished road
I know this might be a little late for you, but in an LES solution the unresolved energy is simply:

Quote:
autoPtr<incompressible::turbulenceModel> model
(
incompressible::turbulenceModel::New(U, phi, laminarTransport)
);
volScalarField k(model->k());
anishtain4 is offline   Reply With Quote

Old   November 30, 2016, 11:05
Default
  #10
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 12
wenxu is on a distinguished road
Quote:
Originally Posted by anishtain4 View Post
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
Bashar and Gang Wang like this.
wenxu is offline   Reply With Quote

Old   December 31, 2016, 08:55
Default
  #11
New Member
 
Arone Brake
Join Date: Dec 2015
Posts: 1
Rep Power: 0
arone is on a distinguished road
How to calculate resolved and unresolved kinetic energy in paraview from openfoam data?
Bashar and Gang Wang like this.
arone is offline   Reply With Quote

Old   January 13, 2017, 18:47
Default
  #12
Member
 
Bashar
Join Date: Jul 2015
Posts: 74
Rep Power: 10
Bashar is on a distinguished road
Hi Wenxu

Can you please share how to practically solve it? i.e. the step to do that? if you have time.
Thank you
Bashar is offline   Reply With Quote

Old   January 18, 2017, 13:01
Default
  #13
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 12
wenxu is on a distinguished road
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.
wenxu is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 05:04.