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

What does the Boussinesq eddy viscosity hypothesis really means?

Register Blogs Community New Posts Updated Threads Search

Like Tree9Likes
  • 3 Post By FMDenaro
  • 2 Post By sbaffini
  • 1 Post By watermelon
  • 2 Post By LuckyTran
  • 1 Post By sbaffini

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 3, 2022, 23:17
Question What does the Boussinesq eddy viscosity hypothesis really means?
  #1
New Member
 
Songrui LI
Join Date: Jul 2022
Location: shanghai
Posts: 3
Rep Power: 3
Ryan LI is on a distinguished road
Hi CFDers

I just stepped into the mathematics of the RANS model and I got stuck. When I read the Boussinesq hypothesis, also called the eddy viscosity hypothesis, I saw this equation and the name of every term. Such as

\tau_{ij} = 2 \, \mu_t \, S_{ij}^* - \frac{2}{3} \rho k \delta_{ij}

where \mu_t denotes the eddy viscosity, S_{ij}^* the strain rate tensor and k the turbulent kinetic energy.

But I still have concerns,

After some research, here are what I know so far:
the strain rate tensor describes the rate of change of the deformation of the fluid in the neighborhood of a certain point. The turbulent momentum transfer results in the deformation of the fluid, and the deformation gives rise to viscous forces in its interior against the deformation, due to friction between adjacent fluid elements, and due to the internal stress tensor.

And the turbulent kinetic energy defines the energy of motion carried by fluid flow. And it dissipates down the turbulence energy cascade.

Here are what I still don't know:
Why the kinetic energy term in the above equation? I don't know how the turbulent kinetic energy is related to the stress tensor.
Why 2 and - \frac{2}{3} in the above equation? I don't know whether this is pure empirical or if there is a theory behind these factors.

Thank you
Ryan LI is offline   Reply With Quote

Old   July 4, 2022, 11:55
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,776
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Ryan LI View Post
Hi CFDers

I just stepped into the mathematics of the RANS model and I got stuck. When I read the Boussinesq hypothesis, also called the eddy viscosity hypothesis, I saw this equation and the name of every term. Such as

\tau_{ij} = 2 \, \mu_t \, S_{ij}^* - \frac{2}{3} \rho k \delta_{ij}

where \mu_t denotes the eddy viscosity, S_{ij}^* the strain rate tensor and k the turbulent kinetic energy.

But I still have concerns,

After some research, here are what I know so far:
the strain rate tensor describes the rate of change of the deformation of the fluid in the neighborhood of a certain point. The turbulent momentum transfer results in the deformation of the fluid, and the deformation gives rise to viscous forces in its interior against the deformation, due to friction between adjacent fluid elements, and due to the internal stress tensor.

And the turbulent kinetic energy defines the energy of motion carried by fluid flow. And it dissipates down the turbulence energy cascade.

Here are what I still don't know:
Why the kinetic energy term in the above equation? I don't know how the turbulent kinetic energy is related to the stress tensor.
Why 2 and - \frac{2}{3} in the above equation? I don't know whether this is pure empirical or if there is a theory behind these factors.

Thank you



I suggest you to think about the decomposition of the tensor in isotropic and deviatoric parts. You will see the 1/3 appearing.

The model is for the deviatoric part of the stress, thus the isotropic part will be added to the pressure. The 2 appears exactly as in in the molecular diffusion term.
FMDenaro is offline   Reply With Quote

Old   July 7, 2022, 02:53
Default
  #3
New Member
 
Songrui LI
Join Date: Jul 2022
Location: shanghai
Posts: 3
Rep Power: 3
Ryan LI is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I suggest you to think about the decomposition of the tensor in isotropic and deviatoric parts. You will see the 1/3 appearing.

The model is for the deviatoric part of the stress, thus the isotropic part will be added to the pressure. The 2 appears exactly as in in the molecular diffusion term.
Thx Filippo for your kindly response

I followed your lead and I found this section demonstrates what you say and it really helps me.

https://en.wikipedia.org/wiki/Viscou...viscous_stress
Ryan LI is offline   Reply With Quote

Old   July 4, 2022, 15:56
Default
  #4
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,151
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Take the trace of both sides and see what happens
LuckyTran and Ryan LI like this.
sbaffini is offline   Reply With Quote

Old   July 7, 2022, 02:55
Thumbs up
  #5
New Member
 
Songrui LI
Join Date: Jul 2022
Location: shanghai
Posts: 3
Rep Power: 3
Ryan LI is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Take the trace of both sides and see what happens
Yes indeed, I think I should look deeper into the linear algebra
Ryan LI is offline   Reply With Quote

Old   July 7, 2022, 15:40
Default
  #6
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,151
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Quote:
Originally Posted by Ryan LI View Post
Yes indeed, I think I should look deeper into the linear algebra
No, please, don't go into linear algebra, just, literally, take the trace of both sides, which are 3x3 tensors, so just take the sum of their 3 diagonal terms and take note of the definitions of all the terms
Ryan LI likes this.
sbaffini is offline   Reply With Quote

Old   July 4, 2022, 20:55
Default
  #7
New Member
 
Xingguang Zhou
Join Date: May 2022
Location: CHINA
Posts: 13
Rep Power: 4
watermelon is on a distinguished road
Quote:
Originally Posted by Ryan LI View Post
Hi CFDers

I just stepped into the mathematics of the RANS model and I got stuck. When I read the Boussinesq hypothesis, also called the eddy viscosity hypothesis, I saw this equation and the name of every term. Such as

\tau_{ij} = 2 \, \mu_t \, S_{ij}^* - \frac{2}{3} \rho k \delta_{ij}

where \mu_t denotes the eddy viscosity, S_{ij}^* the strain rate tensor and k the turbulent kinetic energy.

But I still have concerns,

After some research, here are what I know so far:
the strain rate tensor describes the rate of change of the deformation of the fluid in the neighborhood of a certain point. The turbulent momentum transfer results in the deformation of the fluid, and the deformation gives rise to viscous forces in its interior against the deformation, due to friction between adjacent fluid elements, and due to the internal stress tensor.

And the turbulent kinetic energy defines the energy of motion carried by fluid flow. And it dissipates down the turbulence energy cascade.

Here are what I still don't know:
Why the kinetic energy term in the above equation? I don't know how the turbulent kinetic energy is related to the stress tensor.
Why 2 and - \frac{2}{3} in the above equation? I don't know whether this is pure empirical or if there is a theory behind these factors.

Thank you
Hello, Songrui, Boussinesq eddy viscosity hypothesis is totally analogized from the N-S constitutive relationship. So, if you were familiar with the N-S relationship, the Boussinesq eddy viscosity will be understood too. This will interpret the coefficient 2.

About the coefficient 2/3, this is from the average of turbulent kinematic energy. Actually here is Pt, a "pressure" due to fluctuating velocity.

Best Regards.
Ryan LI likes this.
watermelon is offline   Reply With Quote

Old   July 5, 2022, 10:53
Default
  #8
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,676
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Really it's about understanding where the 3 in 1/3 comes from. If you're asking why 2, then look at the definition of S and tell me why there is a 1/2 in it.
sbaffini and Ryan LI like this.
LuckyTran is offline   Reply With Quote

Reply

Tags
turbulence modeling


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
Boussinesq eddy viscosity hypothesis in cylindrical coordinates? LuckyTran Main CFD Forum 5 August 12, 2015 09:28
Very low eddy viscosity value in turbulent flow mejahan CFX 14 February 3, 2015 15:55
Eddy viscosity ratio in interDyMFoam? wavemaster OpenFOAM Running, Solving & CFD 1 March 9, 2014 03:57
Calculating Eddy Viscosity in order to visualize at Post-Processing. RTom OpenFOAM Pre-Processing 0 September 4, 2013 14:43
Eddy viscosity hypothesis versus Reynolds stress models JuPa ANSYS 0 August 12, 2013 06:20


All times are GMT -4. The time now is 07:47.