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

Query regarding dissipation term of TKE equation

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By FMDenaro
  • 1 Post By LuckyTran

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 5, 2019, 09:34
Post Query regarding dissipation term of TKE equation
  #1
New Member
 
VL Srinivas
Join Date: Jan 2017
Posts: 4
Rep Power: 9
srinivasvl81 is on a distinguished road
In TKE equation, the dissipation term can be written as:

\epsilon = -2nu*<s_ij s_ij>

where, nu=mu/rho is kinematic viscosity, s_ij= 0.5*{ du_i/dx_j + du_j/dx_i }, mu=dynamic viscosity, rho=density.

here, u= fluctuating velocity component, d=partial derivative operator, <>=time average

Considering a problem of sub-meso scale ocean modelling,

1. Is the density, instantaneous or fluctuating or mean ?

2. After doing a time average<> on the strain rate tensors, s_ij there will be no time dimension left, if the density is instantaneous, will the <s_ij s_ij> be constant over all time steps of density ?

3. If the density is constant, will be dissipation also be constant?
srinivasvl81 is offline   Reply With Quote

Old   June 5, 2019, 10:56
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by srinivasvl81 View Post
In TKE equation, the dissipation term can be written as:

\epsilon = -2nu*<s_ij s_ij>

where, nu=mu/rho is kinematic viscosity, s_ij= 0.5*{ du_i/dx_j + du_j/dx_i }, mu=dynamic viscosity, rho=density.

here, u= fluctuating velocity component, d=partial derivative operator, <>=time average

Considering a problem of sub-meso scale ocean modelling,

1. Is the density, instantaneous or fluctuating or mean ?

2. After doing a time average<> on the strain rate tensors, s_ij there will be no time dimension left, if the density is instantaneous, will the <s_ij s_ij> be constant over all time steps of density ?

3. If the density is constant, will be dissipation also be constant?





The general analysis you find for the turbulence has a hystorical assumption, that is the density is homogeneous. For the ocean it can be introduced the Bousinnesq approximation for the bouyancy to take into account the dependence of density on low temperature gradients (actually you could also add the salinity gradient). However, the kinematic viscosity is assumed to be constant.

Have a look here for the ocean flow model https://www.researchgate.net/publica...physical_flows


Note that for general non-constant density flow problem, the Favre averaging is introduced.
srinivasvl81 likes this.
FMDenaro is offline   Reply With Quote

Old   June 5, 2019, 10:58
Default
  #3
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,762
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
nu also needs to be averaged. Your formula implicitly assumes nu is a constant. Use a formula that doesn't assume this and a lot of your confusion will be cleared away.

Btw, Reynolds averaging works also for unsteady flows, you can extend the <> operation beyond time-averaging and allow for time-dependent flows. Think URANS!

However, the answer is....

1. Density is the reynolds average (which can be the mean, but can also be time-dependent like in unsteady RANS). But it won't be instantaneous, otherwise the TKE equation would not be a TKE equation.
srinivasvl81 likes this.
LuckyTran is offline   Reply With Quote

Old   June 6, 2019, 04:30
Post
  #4
New Member
 
VL Srinivas
Join Date: Jan 2017
Posts: 4
Rep Power: 9
srinivasvl81 is on a distinguished road
Thank you for the replies.

By the way, in the <s_ij s_ij> term, is the multiplication, a dot product or a normal matrix product ?
srinivasvl81 is offline   Reply With Quote

Old   June 6, 2019, 04:35
Default
  #5
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by srinivasvl81 View Post
Thank you for the replies.

By the way, in the <s_ij s_ij> term, is the multiplication, a dot product or a normal matrix product ?



It stands for <S:S>, that means you perform the doct product twice
FMDenaro is offline   Reply With Quote

Reply

Tags
dissipation rate, tke budget


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
Adding diffusion term to interFoam transport equation Gearb0x OpenFOAM Programming & Development 3 February 14, 2023 05:16
Convective term in heat equation for rotating solids vabishek OpenFOAM Programming & Development 8 June 24, 2020 12:10
Discretization of the convective term of the N-S equation Donton Main CFD Forum 1 September 1, 2017 13:01
dissipation term in Navier-Stokes Equation mihirmakwana6 Main CFD Forum 16 July 26, 2015 13:41
Source term energy equation for reactive flows DaIN Main CFD Forum 0 October 6, 2011 16:11


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