CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Missing term in generation of TKE for the k equation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 19, 2021, 11:04
Default Missing term in generation of TKE for the k equation
  #1
Member
 
BM
Join Date: Sep 2021
Posts: 35
Rep Power: 5
blyatman is on a distinguished road
In the equation for G_k in the following link: https://www.afs.enea.it/project/nept...de61.htm#eq-gk, G_k is equal to the product of the Reynolds stress and the velocity gradient. However, it then says G_k=2\mu_tS_{ij}S_{ij}, where S_{ij} is the mean rate of strain tensor, which I'm assuming is just 0.5(\partial_iu_j+\partial_ju_i).

Now, the full Reynolds stress term contains a term for the TKE, as shown here:
https://www.cfd-online.com/Wiki/Bous...ity_assumption

In the Fluent theory, G_k clearly does not contain the TKE term. As far as I can tell, it does not seem to absorbed into another term in the TKE equation. Why is the TKE term not present in the definition of G_k? Is it neglected because it is insignificant?
blyatman is offline   Reply With Quote

Old   September 19, 2021, 14:02
Default
  #2
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
A proper tensorial closure model should give you 4.4-21.

However, the Boussinesq hypothesis deals only with the traceless part of the Reynolds stress tensor and that's why it is written that way (with the -2/3 k magically showing up out of nowhere).

The (2/3)k is often absorbed into the pressure term explicitly or implicitly since there's no way for a solver to differentiate between pressure and 2/3 k unless you couple k into the momentum equations. It usually is very small. So yes, it magically disappears as well.

k is proportional to UI (mean velocity x turbulence intensity) whereas the dynamic pressure is like U^2/2

To fully track across different sources where is the two=thirds k is a bit tedious because you have to write down the full NS equations and all the constitutive relations on the same page and their definitions without vague references on relying on people to know what is a traceless strain rate tensor (because that could be any of 4 different tensors).
LuckyTran is offline   Reply With Quote

Old   September 19, 2021, 21:59
Default
  #3
Member
 
BM
Join Date: Sep 2021
Posts: 35
Rep Power: 5
blyatman is on a distinguished road
Thanks, this mostly makes sense. I'm a bit confused about the part where you said another reason it's done is to decouple the momentum equation from the k equation. From the Boussinesq hypothesis, the Reynolds stress term is contains the eddy viscosity, which depends on k. Thus, aren't the momentum and k equations coupled? Intuitively, this is what I would expect since the turbulence affects the mean flow field. However, I understand the part about how k is small, so it is neglected anyway (regardless of whether it's coupled or not).
blyatman is offline   Reply With Quote

Old   September 19, 2021, 22:26
Default
  #4
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
Quote:
Originally Posted by blyatman View Post
I'm a bit confused about the part where you said another reason it's done is to decouple the momentum equation from the k equation.

Coupling between independent transport equations is different. Of course everything is coupled otherwise we would not need to solve them at the same time (temperature/energy also affects the flow field even though temperature never appears anywhere in the momentum equations). I'm talking about k appearing explicitly in the momentum equations.


If i have a system that I think is:
Code:
ax+by+(p1)=0
but in reality the real system is:
Code:
ax+by+(p2+2/3 k)=0
where the a, x, b, y are the same...

If I solve the first system for p1, I will still end up solving for p2+2/3 k whether I like it or not. Even if I use a magic wand to make the two-thirds k disappear analytically, numerically it shows ups because the a,x,b,y are the same. So the pressure that you solve for from your CFD code isn't actually pressure unless you explicitly put the k in it and solve the second system.
LuckyTran is offline   Reply With Quote

Reply


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
Reading a list of scalar data as a constant source term in the governing equation Lewis Liang OpenFOAM Running, Solving & CFD 2 December 30, 2019 00:08
Add different source term in diffusion equation at each time step Lewis Liang OpenFOAM Programming & Development 1 June 7, 2018 11:10
Discretization of the convective term of the N-S equation Donton Main CFD Forum 1 September 1, 2017 13:01
TKE budget equation in LES doctorWho Main CFD Forum 3 March 12, 2014 12:27
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45


All times are GMT -4. The time now is 03:43.