CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Resolved TKE budget terms calculation from LES Simulations

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 29, 2021, 13:33
Default Resolved TKE budget terms calculation from LES Simulations
  #1
New Member
 
Join Date: Aug 2020
Posts: 19
Rep Power: 6
Agavi is on a distinguished road
Hello everyone,

I am running LES and need to plot the turbulent kinetic energy budget of my resolved flow field. I therefore need the equation of the TKE that comes from the square of the u'_{resolved} = \tilde{u}-\bar{u} components.

TKE Budget equations that I have found online are in the physical sense (ie the one in Pope's book) where the filtering operation is obviously ignored. A DNS simulation could use those equations for the calculation of the TKE budget terms.

I'm looking for the derivation of the resolved TKE equation. This probably should come from a similar derivation to the TKE equation but also including the SGS tensor.

-- This old thread discusses the matter but didn't reach a conclusion - That's why I started this new thread.

TKE budget equation in LES


-- I am probably looking for the derivation of equation (1) in this paper. It accounts for the effect of \tau_{ij}^{SGS} to the resolved TKE balance.

https://www.osti.gov/servlets/purl/1500073

Please let me know if I'm understanding something wrong here!

Thanks

Agave

Last edited by Agavi; May 29, 2021 at 14:49.
Agavi is offline   Reply With Quote

Old   May 29, 2021, 14:33
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Agavi View Post
Hello everyone,

I am running LES and need to plot the turbulent kinetic energy budget of my resolved flow field. I therefore need the equation of the TKE that comes from the square of the u'_{resolved} = \tilde{u}-\bar{u} components.

TKE Budget equations that I have found online are in the physical sense (ie the one in Pope's book) where the filtering operation is obviously ignored. A DNS simulation could use those equations for the calculation of the TKE budget terms.

I'm looking for the derivation of the resolved TKE equation. This probably should come from a similar derivation to the TKE equation but also including the SGS tensor.

-- This old thread discusses the matter but didn't reach a conclusion - That's why I started this new thread.

TKE budget equation in LES

-- Also in this OpenFOAM code the SGS stress tensor is not taken into account for the calculation of the resolved TKE budget and I am not sure why

https://github.com/syavash20/TurbLab.../LESpimpleFoam

-- I am probably looking for the derivation of equation (1) in this paper. It accounts for the effect of \tau_{ij}^{SGS} to the resolved TKE balance.

https://www.osti.gov/servlets/purl/1500073

Please let me know if I'm understanding something wrong here!

Thanks

Agave



What are you searching for, the equation for vbar.vbar or the equation for (v'.v')bar ?
FMDenaro is offline   Reply With Quote

Old   May 29, 2021, 14:38
Default
  #3
New Member
 
Join Date: Aug 2020
Posts: 19
Rep Power: 6
Agavi is on a distinguished road
Hi FMDenaro,

Thanks for your reply.

I'm searching for the equation for the resolved by the LES turbulent kinetic energy, so (u_{resloved}'u_{resolved}')_{bar}.

Best

Agave
Agavi is offline   Reply With Quote

Old   May 29, 2021, 14:45
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Well, if I understand you, you need first to evaluate the statistical averaging of the resolved field, that is <vbar> (<*> is a time averaging) which is a steady field. Then you want to compute the fluctuation field vbar-<vbar> which is an unsteady field representing the part of the fluctuations you can resolve (of course these are not the DNS-based one).

In practice you need to sample the LES velocity for several time-unit and then compute the statistical mean.
FMDenaro is offline   Reply With Quote

Old   May 29, 2021, 15:05
Default
  #5
New Member
 
Join Date: Aug 2020
Posts: 19
Rep Power: 6
Agavi is on a distinguished road
Thanks for your reply!

I do not only need the components (u'u')_bar. I need the whole transport equation of the resolved turbulent kinetic energy (so production, dissipation and diffusion terms - what is explained in post #1 in the thread I referenced TKE budget equation in LES).
Agavi is offline   Reply With Quote

Old   May 29, 2021, 15:40
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Agavi View Post
Thanks for your reply!

I do not only need the components (u'u')_bar. I need the whole transport equation of the resolved turbulent kinetic energy (so production, dissipation and diffusion terms - what is explained in post #1 in the thread I referenced TKE budget equation in LES).



Have you read also the textbook of Sagaut? You will find the equations
FMDenaro is offline   Reply With Quote

Old   May 29, 2021, 16:04
Default
  #7
New Member
 
Join Date: Aug 2020
Posts: 19
Rep Power: 6
Agavi is on a distinguished road
I looked at the equations and from my understanding eqn. 3.30 is the total resolvet kinetic energy (so it comes from UTilde). I'm not sure I understand equation 3.33. What is filtered subgrid kinetic energy?

If our real physical flow field is

U = \bar{U}+ U'_{resolved}+U'_{modelled}

Then it is my understanding that subgrid kinetic energy has to come from (U'_{modelled}U'_{modelled})_{bar}. However I need the equation for the transport of (U'_{resolved}U'_{resolved})_{bar}.

What am I understanding wrong?
Agavi is offline   Reply With Quote

Old   May 29, 2021, 16:19
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Agavi View Post
I looked at the equations and from my understanding eqn. 3.30 is the total resolvet kinetic energy (so it comes from UTilde). I'm not sure I understand equation 3.33. What is filtered subgrid kinetic energy?

If our real physical flow field is

U = \bar{U}+ U'_{resolved}+U'_{modelled}

Then it is my understanding that subgrid kinetic energy has to come from (U'_{modelled}U'_{modelled})_{bar}. However I need the equation for the transport of (U'_{resolved}U'_{resolved})_{bar}.

What am I understanding wrong?



Have a look to this slide. Remember that the fluctuations in LES are always modelled. What you can evaluate are the fluctuations with respect to the statistical mean.
Attached Images
File Type: jpg Image1.jpg (142.7 KB, 31 views)
FMDenaro is offline   Reply With Quote

Old   June 2, 2021, 14:19
Default
  #9
New Member
 
Join Date: Aug 2020
Posts: 19
Rep Power: 6
Agavi is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Remember that the fluctuations in LES are always modelled.
Hi Fillipo,

I'm not sure I understand that. What I know is that in LES some of the turbulent fluctuations are resolved and the rest of them are modelled (subgrid scale fluctuations).

Best

Agave
Agavi is offline   Reply With Quote

Old   June 2, 2021, 14:26
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Agavi View Post
Hi Fillipo,

I'm not sure I understand that. What I know is that in LES some of the turbulent fluctuations are resolved and the rest of them are modelled (subgrid scale fluctuations).

Best

Agave



In LES you assume a decomposion in resolved (filtered) and unresolved components. By definition, the fluctuations of the velocity (with respect to the filtered velocity) are only in the unresolved components.
The resolved components are unsteady and 3d as same as the fluctuations.

Could you address the reference you have read?
FMDenaro is offline   Reply With Quote

Old   June 2, 2021, 14:41
Default
  #11
New Member
 
Join Date: Aug 2020
Posts: 19
Rep Power: 6
Agavi is on a distinguished road
Oh OK. I thought that by 'fluctuations' you meant wrt to the mean flow. Now I get you mean the fluctuations wrt to the filtered field that is only modelled in LES.

Now that we have established the definitions, I'm looking for the tke transport equation for the resolved LES fluctuations with respect to the mean flow.
( basically in your post here TKE budget equation in LES looking for the transport equation of 0.5 <(u'_{RANS}-u'_{LES})^2>. ) Therefore I think that the equations in Sagaut and in the slides you attached can't work for what I need because they are for the kinetic energy in different regions of the spectrum.

Maybe Im wrong and such an equation cannot be derived. However I'm pretty sure this is what is used in the paper I attached in #1 of this thread and I'm looking for the derivation of it
Agavi is offline   Reply With Quote

Old   June 2, 2021, 16:13
Default
  #12
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
u'RANS(x,t) = u(x,t) -<u>(x)


u'LES(x,t) = u(x,t) -u_bar(x,t)


hence



u'RANS(x,t) -u'LES(x,t) =u_bar(x,t)-<u>(x)


If you want to write an equation for the LHS it would simply be equal do du_bar/dt.

The term you asked for is 0.5*<(u_bar-<u>)^2> which is, by definition of <>, not depending on time.
Is that what you asked for?

Last edited by FMDenaro; June 2, 2021 at 17:11. Reason: I edited the equation that I've written before to correct them.
FMDenaro 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 budget terms with openfoam Thangam OpenFOAM Programming & Development 6 July 15, 2016 09:21
Budget terms in LES shahzeb irfan FLUENT 0 August 10, 2011 00:00
Budget terms in LES shahzeb irfan Main CFD Forum 0 August 7, 2011 23:57
What are BUDGET TERMS in Large Eddy Simulation??? mariachi FLUENT 1 January 23, 2010 03:36
The terms that should be treated implicitly in LES ben Main CFD Forum 3 January 28, 2005 03:32


All times are GMT -4. The time now is 12:42.