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

How do we calculate the normal Reynolds stresses? (Linear eddy viscosity models)

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By userid42
  • 1 Post By LuckyTran
  • 1 Post By sbaffini

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   August 18, 2017, 12:19
Default How do we calculate the normal Reynolds stresses? (Linear eddy viscosity models)
  #1
New Member
 
Daniel
Join Date: Aug 2017
Posts: 7
Rep Power: 8
userid42 is on a distinguished road
I have read in multiple places that the normal Reynolds stresses are equal when the linear eddy viscosity model is used.
If they are equal, per the definition of turbulent kinetic energy:

k = \frac{1}{2} (\overline{u'u'} + \overline{v'v'} + \overline{w'w'})

\overline{u'u'} = \frac{2}{3}k.

From what I understand, we can calculate the Reynolds stresses from the Boussinesq approximation. In this approximation, the eddy viscosity is assumed to be isotropic (not turbulence).

The Boussinesq approximation is:

-\rho\overline{u^{'}_{i}u^{'}_{j}}=\mu_t\left(\frac{\partial\overline{u_i}}{\partial{x_j}}+\frac{\partial\overline{u_j}}{\partial{x_i}}-\frac{2}{3}\frac{\partial\overline{u_k}}{\partial{x_k}}\delta_{ij}\right)-\frac{2}{3}\rho k\delta_{ij}.

From this equation, the normal stresses (for incompressible flows) can be calculated by:

\overline{u^{'}u^{'}}=\frac{2}{3}k - \frac{\mu_t}{\rho}\left(\frac{\partial\overline{u}}{\partial{x}}+\frac{\partial\overline{u}}{\partial{x}}\right),

\overline{v^{'}v^{'}}=\frac{2}{3}k - \frac{\mu_t}{\rho}\left(\frac{\partial\overline{v}}{\partial{y}}+\frac{\partial\overline{v}}{\partial{y}}\right),

and

\overline{w^{'}w^{'}}=\frac{2}{3}k - \frac{\mu_t}{\rho}\left(\frac{\partial\overline{w}}{\partial{z}}+\frac{\partial\overline{w}}{\partial{z}}\right).

From the expressions above, I see that the normal stresses are not equal because the velocity gradients will be different. Also, this approach introduces an extra term "-\frac{\mu_t}{\rho}(...)".

Which approach is correct? Am I missing something? Or is this a misconception of what is isotropic?
Feel free to correct me if I stated something wrong.

Thanks!
Ashkan Kashani likes this.

Last edited by userid42; August 18, 2017 at 21:07.
userid42 is offline   Reply With Quote

 


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
Calculate and display eddy viscosity ceyrows OpenFOAM Running, Solving & CFD 4 February 16, 2009 12:57
non-linear eddy viscosity models George Main CFD Forum 0 December 11, 2006 19:14
Why Turbulence models are not universal. Senthil Main CFD Forum 4 July 5, 2000 04:34
Reynolds Stress Models Roued Main CFD Forum 20 February 8, 2000 02:58
Reynolds Stress Models Roued CFX 2 February 7, 2000 10:33


All times are GMT -4. The time now is 00:45.