CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Issues on modelling of isotropic part of Reynolds (SGS) stresses

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

Like Tree1Likes
  • 1 Post By Likun

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 4, 2015, 04:11
Default Issues on modelling of isotropic part of Reynolds (SGS) stresses
  #1
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 13
Likun is on a distinguished road
Send a message via Skype™ to Likun
Hi Foamers,

Recently I am confused about the modelling of Reynolds or SGS stresses in openFoam. Specifically, I don't understand how the isotropic part the the SGS stresses have been considered. According to the turbulent viscosity assumption, the Reynolds or SGS stresses
\tau_{ij} = -2 \mu_{sgs} \tilde S_{ij}^{D} + \frac{1}{3} \tau_{ij} \delta_{ij} = -2 \mu_{sgs} \tilde S_{ij}^{D} + \frac{2}{3}\bar \rho k \delta_{ij}

where \tilde S_{ij}^{D} = \frac{1}{2}\left( \frac{\partial \tilde u_{i}}{\partial x_{j}} + \frac{\partial \tilde u_{j}}{\partial x_{i}} \right) - \frac{1}{3} \frac{\partial \tilde u_{k}}{\partial x_{k}} \delta_{ij} is deviatoric part of the local strain rate tensor.

As discussed in http://www.cfd-online.com/Forums/ope...sonicfoam.html and http://www.cfd-online.com/Forums/openfoam-solving/58214-calculating-divdevreff.html#post388688, , the isotropic part of Reynolds (SGS) stresses, \frac{1}{3} \tau_{ij} \delta_{ij} or \frac{2}{3}\bar \rho k \delta_{ij}, are not included in the turbulence->divDevRhoReff(U) term of UEqn
Code:
fvVectorMatrix UEqn 
(    
      fvm::ddt(rho, U)     
    + fvm::div(phi, U)   
    + turbulence->divDevRhoReff(U) 
);
And I can not find the source code, where this term has been considered. So I suppose that this term has been simply ignored in openFoam! Someone please correct me if I am wrong.

In many literature, including the book of Pope, it says that this term is absorbed into a modified pressure which leads to the following modelled momentum equation:

\frac{\partial \bar \rho \tilde u_{i}}{\partial t} + \frac{\partial \left(\bar \rho \tilde u_{i} \tilde u_{j} \right)}{\partial x_{j}} = -\frac{\partial \bar p}{\partial x_{i}} \left( \bar p + \frac{2}{3} k \right) + \frac{\partial}{\partial j} \left(2 \left(\mu + \mu_{sgs} \right) \tilde S_{ij}^{D} \right) + \bar S_{M,i}

It is clear that if the turbulent (SGS) kinetic energy k is much smaller than the thermodynamic pressure p, this term is indeed negligible.

So my questions are:

1. Is it true that the \frac{1}{3} \tau_{ij} \delta_{ij}term has been neglected in OpenFoam? If this is true, the pressure (in compressible part) is the real pressure, right?

2. How this term can be included, since we already have the value for k?

Any idea will be appreciated.

Best,
Likun
ztnuaa likes this.
Likun is offline   Reply With Quote

Old   June 25, 2015, 05:57
Default
  #2
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Stockholm, Sweden
Posts: 119
Rep Power: 14
tiam is on a distinguished road
Realized that this issue is not really clear to me either . Any answers since 2013?
tiam is offline   Reply With Quote

Old   January 25, 2016, 02:51
Default
  #3
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 13
Likun is on a distinguished road
Send a message via Skype™ to Likun
New discussions on this issue can find here: http://www.cfd-online.com/Forums/ope...tml#post582044
Likun 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
[snappyHexMesh] crash sHM H25E OpenFOAM Meshing & Mesh Conversion 11 November 10, 2014 11:27
LES modelling issues with waveFoam Agni OpenFOAM 0 September 17, 2012 10:33
Istropic part of Reynolds stress tensor AlmostSurelyRob Main CFD Forum 6 September 11, 2012 09:58
Sampling Reynolds Stresses sven OpenFOAM 4 March 1, 2012 16:15
Deviatoric Part of the reynolds tress tensor make OpenFOAM Running, Solving & CFD 2 September 22, 2008 14:02


All times are GMT -4. The time now is 10:24.