Turbulence modeling

<h2 id="closure-problem">Closure problem</h2>
The <a href="/facts/Navier%E2%80%93Stokes_equations/48pQugEa">Navier–Stokes equations</a> govern the velocity and pressure of a fluid flow. In a turbulent flow, each of these quantities may be decomposed into a mean part and a fluctuating part. Averaging the equations gives the <a href="/facts/Reynolds-averaged_Navier%E2%80%93Stokes_equations/WpGIsv7x">Reynolds-averaged Navier–Stokes (RANS) equations</a>, which govern the mean flow. However, the nonlinearity of the Navier–Stokes equations means that the velocity fluctuations still appear in the RANS equations, in the nonlinear term 
 
 
 
 −
 ρ
 
 
 
 
 v
 
 i
 
 
 ′
 
 
 
 v
 
 j
 
 
 ′
 
 
 
 ¯
 
 
 
 
 {\displaystyle -\rho {\overline {v_{i}^{\prime }v_{j}^{\prime }}}}
 
 from the convective acceleration. This term is known as the <a href="/facts/Reynolds_stress/PY1P33IL">Reynolds stress</a>, 
 
 
 
 
 R
 
 i
 j
 
 
 
 
 {\displaystyle R_{ij}}
 
.<a class="footnote-ref" id="fnref:2" href="#fn:2">2</a> Its effect on the mean flow is like that of a stress term, such as from pressure or viscosity.
To obtain equations containing only the mean velocity and pressure, we need to close the RANS equations by modelling the Reynolds stress term 
 
 
 
 
 R
 
 i
 j
 
 
 
 
 {\displaystyle R_{ij}}
 
 as a function of the mean flow, removing any reference to the fluctuating part of the velocity. This is the closure problem.

<h2 id="eddy-viscosity">Eddy viscosity</h2>
<a href="/facts/Joseph_Valentin_Boussinesq/m4N2S7FQ">Joseph Valentin Boussinesq</a> was the first to attack the closure problem,<a class="footnote-ref" id="fnref:3" href="#fn:3">3</a> by introducing the concept of <a href="/facts/Viscosity/QLVfrdCz">eddy viscosity</a>. In 1877 Boussinesq proposed relating the turbulence stresses to the mean flow to close the system of equations. Here the Boussinesq hypothesis is applied to model the Reynolds stress term. Note that a new proportionality constant 
 
 
 
 
 ν
 
 t
 
 
 >
 0
 
 
 {\displaystyle \nu _{t}>0}
 
, the (kinematic) turbulence eddy viscosity, has been introduced. Models of this type are known as eddy viscosity models (EVMs).

 
 
 
 −
 
 
 
 
 v
 
 i
 
 
 ′
 
 
 
 v
 
 j
 
 
 ′
 
 
 
 ¯
 
 
 =
 
 ν
 
 t
 
 
 
 (
 
 
 
 
 ∂
 
 
 
 v
 
 i
 
 
 ¯
 
 
 
 
 ∂
 
 x
 
 j
 
 
 
 
 
 +
 
 
 
 ∂
 
 
 
 v
 
 j
 
 
 ¯
 
 
 
 
 ∂
 
 x
 
 i
 
 
 
 
 
 
 )
 
 −
 
 
 2
 3
 
 
 k
 
 δ
 
 i
 j
 
 
 
 
 {\displaystyle -{\overline {v_{i}^{\prime }v_{j}^{\prime }}}=\nu _{t}\left({\frac {\partial {\overline {v_{i}}}}{\partial x_{j}}}+{\frac {\partial {\overline {v_{j}}}}{\partial x_{i}}}\right)-{\frac {2}{3}}k\delta _{ij}}

which can be written in shorthand as

−
        
          
            
              
                v
                
                  i
                
                
                  ′
                
              
              
                v
                
                  j
                
                
                  ′
                
              
            
            ¯
          
        
        =
        2
        
          ν
          
            t
          
        
        
          S
          
            i
            j
          
        
        −
        
          
            
              2
              3
            
          
        
        k
        
          δ
          
            i
            j
          
        
      
    
    {\displaystyle -{\overline {v_{i}^{\prime }v_{j}^{\prime }}}=2\nu _{t}S_{ij}-{\tfrac {2}{3}}k\delta _{ij}}

where

<ul><li>
 
 
 
 
 S
 
 i
 j
 
 
 
 
 {\displaystyle S_{ij}}
 
 is the <a href="/facts/Reynolds-averaged_Navier%E2%80%93Stokes_equations/WpGIsv7x">mean rate of strain tensor</a></li>
<li>
 
 
 
 
 ν
 
 t
 
 
 
 
 {\displaystyle \nu _{t}}
 
 is the (kinematic) turbulence eddy viscosity</li>
<li>
 
 
 
 k
 =
 
 
 
 1
 2
 
 
 
 
 
 
 
 v
 
 i
 
 ′
 
 
 v
 
 i
 
 ′
 
 
 ¯
 
 
 
 
 {\displaystyle k={\tfrac {1}{2}}{\overline {v_{i}'v_{i}'}}}
 
 is the <a href="/facts/Turbulence_kinetic_energy/X5jdtHQg">turbulence kinetic energy</a></li>
<li>and 
 
 
 
 
 δ
 
 i
 j
 
 
 
 
 {\displaystyle \delta _{ij}}
 
 is the <a href="/facts/Kronecker_delta/gcGAy0bm">Kronecker delta</a>.</li></ul>
In this model, the additional turbulence stresses are given by augmenting the <a href="/facts/Molecule/N9ljSfFq">molecular</a> viscosity with an eddy viscosity.<a class="footnote-ref" id="fnref:4" href="#fn:4">4</a> This can be a simple constant eddy viscosity (which works well for some free <a href="/facts/Shear_(fluid)/aWWnzIhE">shear</a> flows such as axisymmetric jets, 2-D jets, and mixing layers).
The Boussinesq hypothesis – although not explicitly stated by Boussinesq at the time – effectively consists of the assumption that the Reynolds stress tensor is aligned with the strain tensor of the mean flow (i.e.: that the <a href="/facts/Shear_stress/aWWnzIhE">shear stresses</a> due to turbulence act in the same direction as the shear stresses produced by the averaged flow). It has since been found to be significantly less accurate than most practitioners would assume.<a class="footnote-ref" id="fnref:5" href="#fn:5">5</a> Still, turbulence models which employ the Boussinesq hypothesis have demonstrated significant practical value. In cases with well-defined shear layers, this is likely due the dominance of streamwise shear components, so that considerable relative errors in flow-normal components are still negligible in absolute terms. Beyond this, most eddy viscosity turbulence models contain coefficients which are calibrated against measurements, and thus produce reasonably accurate overall outcomes for flow fields of similar type as used for calibration.

<h2 id="prandtls-mixing-length-concept">Prandtl's mixing-length concept</h2>
Later, <a href="/facts/Ludwig_Prandtl/I4b4E7sq">Ludwig Prandtl</a> introduced the additional concept of the mixing length,<a class="footnote-ref" id="fnref:6" href="#fn:6">6</a> along with the idea of a <a href="/facts/Boundary_layer/go8gqWM1">boundary layer</a>. For wall-bounded turbulent flows, the eddy viscosity must vary with distance from the wall, hence the addition of the concept of a 'mixing length'. In the simplest wall-bounded flow model, the eddy viscosity is given by the equation:

ν
          
            t
          
        
        =
        
          |
          
            
              
                ∂
                u
              
              
                ∂
                y
              
            
          
          |
        
        
          l
          
            m
          
          
            2
          
        
      
    
    {\displaystyle \nu _{t}=\left|{\frac {\partial u}{\partial y}}\right|l_{m}^{2}}

where

<ul><li>
 
 
 
 
 
 
 ∂
 u
 
 
 ∂
 y
 
 
 
 
 
 {\displaystyle {\frac {\partial u}{\partial y}}}
 
 is the <a href="/facts/Partial_derivative/JWHsBkAu">partial derivative</a> of the streamwise velocity (u) with respect to the wall normal direction (y)</li>
<li>
 
 
 
 
 l
 
 m
 
 
 
 
 {\displaystyle l_{m}}
 
 is the mixing length.</li></ul>
This simple model is the basis for the "<a href="/facts/Law_of_the_wall/2aj4NCEh">law of the wall</a>", which is a surprisingly accurate model for wall-bounded, attached (not separated) flow fields with small <a href="/facts/Pressure_gradient/v7BmlAAS">pressure gradients</a>.
More general turbulence models have evolved over time, with most modern turbulence models given by <a href="/facts/Field_equation/b2LS7sdN">field equations</a> similar to the <a href="/facts/Navier%E2%80%93Stokes_equations/48pQugEa">Navier–Stokes equations</a>.

<h2 id="smagorinsky-model-for-the-sub-grid-scale-eddy-viscosity">Smagorinsky model for the sub-grid scale eddy viscosity</h2>
<a href="/facts/Joseph_Smagorinsky/shUzeX9U">Joseph Smagorinsky</a> was the first who proposed a formula for the eddy viscosity in <a href="/facts/Large_eddy_simulation/ULBNg1SY">Large Eddy Simulation</a> models,<a class="footnote-ref" id="fnref:7" href="#fn:7">7</a> based on the local derivatives of the velocity field and the local grid size:

ν
          
            t
          
        
        =
        Δ
        x
        Δ
        y
        
          
            
              
                (
                
                  
                    
                      ∂
                      u
                    
                    
                      ∂
                      x
                    
                  
                
                )
              
              
                2
              
            
            +
            
              
                (
                
                  
                    
                      ∂
                      v
                    
                    
                      ∂
                      y
                    
                  
                
                )
              
              
                2
              
            
            +
            
              
                1
                2
              
            
            
              
                (
                
                  
                    
                      
                        ∂
                        u
                      
                      
                        ∂
                        y
                      
                    
                  
                  +
                  
                    
                      
                        ∂
                        v
                      
                      
                        ∂
                        x
                      
                    
                  
                
                )
              
              
                2
              
            
          
        
      
    
    {\displaystyle \nu _{t}=\Delta x\Delta y{\sqrt {\left({\frac {\partial u}{\partial x}}\right)^{2}+\left({\frac {\partial v}{\partial y}}\right)^{2}+{\frac {1}{2}}\left({\frac {\partial u}{\partial y}}+{\frac {\partial v}{\partial x}}\right)^{2}}}}

In the context of Large Eddy Simulation, turbulence modeling refers to the need to parameterize the subgrid scale stress in terms of features of the filtered velocity field. This field is called <a href="http://www.scholarpedia.org/article/Turbulence:_Subgrid-Scale_Modeling">subgrid-scale modeling</a>.

<h2 id="spalartallmaras-k-and-k-models">Spalart–Allmaras, k–ε and k–ω models</h2>
The Boussinesq hypothesis is employed in the <a href="/facts/Spalart%E2%80%93Allmaras_turbulence_model/1RlnXcVD">Spalart–Allmaras</a> (S–A), <a href="/facts/K-epsilon_turbulence_model/vzwSeyq0">k–ε</a> (k–epsilon), and <a href="/facts/K-omega_turbulence_model/Pda5KmM1">k–ω</a> (k–omega) models and offers a relatively low cost computation for the turbulence viscosity 
 
 
 
 
 ν
 
 t
 
 
 
 
 {\displaystyle \nu _{t}}
 
. The S–A model uses only one additional equation to model turbulence viscosity transport, while the k–ε and k–ω models use two.

<h2 id="common-models">Common models</h2>
The following is a brief overview of commonly employed models in modern engineering applications.

<ul><li><a href="/facts/Spalart%E2%80%93Allmaras_turbulence_model/1RlnXcVD">Spalart–Allmaras (S–A)</a>

The Spalart–Allmaras model<a class="footnote-ref" id="fnref:8" href="#fn:8">8</a> is a one-equation model that solves a modelled transport equation for the kinematic eddy turbulent viscosity. The Spalart–Allmaras model was designed specifically for aerospace applications involving wall-bounded flows and has been shown to give good results for boundary layers subjected to adverse pressure gradients. It is also gaining popularity in turbomachinery applications.</li><li><a href="/facts/K-epsilon_turbulence_model/vzwSeyq0">k–ε (k–epsilon)</a>

K-epsilon (k-ε) turbulence model<a class="footnote-ref" id="fnref:9" href="#fn:9">9</a> is the most common model used in computational fluid dynamics (CFD) to simulate mean flow characteristics for turbulent flow conditions. It is a two-equation model which gives a general description of turbulence by means of two transport equations (PDEs). The original impetus for the K-epsilon model was to improve the mixing-length model, as well as to find an alternative to algebraically prescribing turbulent length scales in moderate to high complexity flows.</li><li><a href="/facts/K-omega_turbulence_model/Pda5KmM1">k–ω (k–omega)</a>

In computational fluid dynamics, the k–omega (k–ω) turbulence model<a class="footnote-ref" id="fnref:10" href="#fn:10">10</a> is a common two-equation turbulence model that is used as a closure for the Reynolds-averaged Navier–Stokes equations (RANS equations). The model attempts to predict turbulence by two partial differential equations for two variables, k and ω, with the first variable being the turbulence kinetic energy (k) while the second (ω) is the specific rate of dissipation (of the turbulence kinetic energy k into internal thermal energy).</li><li><a href="/facts/SST_(Menter%E2%80%99s_Shear_Stress_Transport)/f1RIeltf">SST (Menter’s Shear Stress Transport)</a>

SST (Menter's shear stress transport) turbulence model<a class="footnote-ref" id="fnref:11" href="#fn:11">11</a> is a widely used and robust two-equation eddy-viscosity turbulence model used in computational fluid dynamics. The model combines the k-omega turbulence model and K-epsilon turbulence model such that the k-omega is used in the inner region of the boundary layer and switches to the k-epsilon in the free shear flow.</li><li><a href="/facts/Reynolds_stress_equation_model/uNg6BR8q">Reynolds stress equation model</a>

The Reynolds stress equation model (RSM), also referred to as second moment closure model,<a class="footnote-ref" id="fnref:12" href="#fn:12">12</a> is the most complete classical turbulence modelling approach. Popular eddy-viscosity based models like the <a href="/facts/K-epsilon_turbulence_model/vzwSeyq0">k–ε (k–epsilon)</a> model and the <a href="/facts/K-omega_turbulence_model/Pda5KmM1">k–ω (k–omega)</a> models have significant shortcomings in complex engineering flows. This arises due to the use of the eddy-viscosity hypothesis in their formulation. For instance, in flows with high degrees of anisotropy, significant streamline curvature, flow separation, zones of recirculating flow or flows influenced by rotational effects, the performance of such models is unsatisfactory.<a class="footnote-ref" id="fnref:13" href="#fn:13">13</a> In such flows, Reynolds stress equation models offer much better accuracy.<a class="footnote-ref" id="fnref:14" href="#fn:14">14</a>

Eddy viscosity based closures cannot account for the return to isotropy of turbulence,<a class="footnote-ref" id="fnref:15" href="#fn:15">15</a> observed in decaying turbulent flows. Eddy-viscosity based models cannot replicate the behaviour of turbulent flows in the Rapid Distortion limit,<a class="footnote-ref" id="fnref:16" href="#fn:16">16</a> where the turbulent flow essentially behaves like an elastic medium.<a class="footnote-ref" id="fnref:17" href="#fn:17">17</a></li></ul>

<h3>Notes</h3>

<h3>Other</h3>
<ul><li>Absi, R. (2019) "Eddy Viscosity and Velocity Profiles in Fully-Developed Turbulent Channel Flows" Fluid Dyn (2019) 54: 137. <a href="https://doi.org/10.1134/S0015462819010014">https://doi.org/10.1134/S0015462819010014</a></li>
<li>Absi, R. (2021) "Reinvestigating the Parabolic-Shaped Eddy Viscosity Profile for Free Surface Flows" Hydrology 2021, 8(3), 126. <a href="https://doi.org/10.3390/hydrology8030126">https://doi.org/10.3390/hydrology8030126</a></li>
<li>Townsend, A. A. (1980) "The Structure of Turbulent Shear Flow" 2nd Edition (Cambridge Monographs on Mechanics), <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 0521298199</li>
<li>Bradshaw, P. (1971) "An introduction to turbulence and its measurement" (Pergamon Press), <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 0080166210</li>
<li>Wilcox, C. D. (1998), "Turbulence Modeling for CFD" 2nd Ed., (DCW Industries, La Cañada), <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 0963605100</li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1">Pope, Stephen (2000). Turbulent Flows. <a href="#fnref:1" class="footnote-back-ref">↩</a></li>
<li id="fn:2">Andersson, Bengt; et al. (2012). Computational fluid dynamics for engineers. Cambridge: Cambridge University Press. p. 83. ISBN 978-1-107-01895-2. <a href="978-1-107-01895-2" target="_blank">978-1-107-01895-2</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></li>
<li id="fn:3">Boussinesq, Joseph (1903). Boussinesq, J. (1903). Thōrie analytique de la chaleur mise en harmonie avec la thermodynamique et avec la thōrie mc̄anique de la lumi_re: Refroidissement et c̄hauffement par rayonnement, conductibilit ̄des tiges, lames et masses cristallines, courants de convection, thōrie mc̄anique de la lumi_re. Gauthier-Villars. <a href="#fnref:3" class="footnote-back-ref">↩</a></li>
<li id="fn:4">John J. Bertin; Jacques Periaux; Josef Ballmann (1992), Advances in Hypersonics: Modeling hypersonic flows, Springer, ISBN 9780817636630 <a href="9780817636630" target="_blank">9780817636630</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></li>
<li id="fn:5">François G. Schmitt (2007), "About Boussinesq's turbulent viscosity hypothesis: historical remarks and a direct evaluation of its validity", Comptes Rendus Mécanique, 335 (9–10): 617–627, doi:10.1016/j.crme.2007.08.004, hdl:20.500.12210/73178, S2CID 32637068 <a href="https://hal.science/hal-00264386/en/" target="_blank">https://hal.science/hal-00264386/en/</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></li>
<li id="fn:6">Prandtl, Ludwig (1925). "Bericht uber Untersuchungen zur ausgebildeten Turbulenz". Z. Angew. Math. Mech. 5 (2): 136. Bibcode:1925ZaMM....5..136P. doi:10.1002/zamm.19250050212. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></li>
<li id="fn:7">Smagorinsky, Joseph (1963). "Smagorinsky, Joseph. "General circulation experiments with the primitive equations: I. The basic experiment". Monthly Weather Review. 91 (3): 99–164. Bibcode:1963MWRv...91...99S. doi:10.1175/1520-0493(1963)091<0099:GCEWTP>2.3.CO;2. <a href="https://doi.org/10.1175%2F1520-0493%281963%29091%3C0099%3AGCEWTP%3E2.3.CO%3B2" target="_blank">https://doi.org/10.1175%2F1520-0493%281963%29091%3C0099%3AGCEWTP%3E2.3.CO%3B2</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></li>
<li id="fn:8">Spalart, Philippe R.; Allmaras, Steven R. (1992). "A one-equation turbulence model for aerodynamic flows". 30th Aerospace Sciences Meeting and Exhibit, AIAA. doi:10.2514/6.1992-439. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></li>
<li id="fn:9">Hanjalic, K.; Launder, B. (1972). "A Reynolds stress model of turbulence and its application to thin shear flows". Journal of Fluid Mechanics. 52 (4): 609–638. Bibcode:1972JFM....52..609H. doi:10.1017/S002211207200268X. S2CID 122631170. <a href="https://www.researchgate.net/publication/233832584" target="_blank">https://www.researchgate.net/publication/233832584</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></li>
<li id="fn:10">Wilcox, D. C. (2008). "Formulation of the k-omega Turbulence Model Revisited". AIAA Journal. 46 (11): 2823–2838. Bibcode:2008AIAAJ..46.2823W. doi:10.2514/1.36541. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></li>
<li id="fn:11">Menter, F. R. (1994). "Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications" (PDF). AIAA Journal. 32 (8): 1598–1605. Bibcode:1994AIAAJ..32.1598M. doi:10.2514/3.12149. S2CID 120712103.[dead link‍] <a href="http://www.academia.edu/download/54560714/1._Two_Equation_Eddy_Viscosity_Turbulence_Models_for_Engineering.pdf" target="_blank">http://www.academia.edu/download/54560714/1._Two_Equation_Eddy_Viscosity_Turbulence_Models_for_Engineering.pdf</a> <a href="#fnref:11" class="footnote-back-ref">↩</a></li>
<li id="fn:12">Hanjalić, Hanjalić; Launder, Brian (2011). Modelling Turbulence in Engineering and the Environment: Second-Moment Routes to Closure. <a href="#fnref:12" class="footnote-back-ref">↩</a></li>
<li id="fn:13">Mishra, Aashwin; Girimaji, Sharath (2013). "Intercomponent energy transfer in incompressible homogeneous turbulence: multi-point physics and amenability to one-point closures". Journal of Fluid Mechanics. 731: 639–681. Bibcode:2013JFM...731..639M. doi:10.1017/jfm.2013.343. S2CID 122537381. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:13" class="footnote-back-ref">↩</a></li>
<li id="fn:14">Pope, Stephen. "Turbulent Flows". Cambridge University Press, 2000. <a href="#fnref:14" class="footnote-back-ref">↩</a></li>
<li id="fn:15">Lumley, John; Newman, Gary (1977). "The return to isotropy of homogeneous turbulence". Journal of Fluid Mechanics. 82: 161–178. Bibcode:1977JFM....82..161L. doi:10.1017/s0022112077000585. S2CID 39228898. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:15" class="footnote-back-ref">↩</a></li>
<li id="fn:16">Mishra, Aashwin; Girimaji, Sharath (2013). "Intercomponent energy transfer in incompressible homogeneous turbulence: multi-point physics and amenability to one-point closures". Journal of Fluid Mechanics. 731: 639–681. Bibcode:2013JFM...731..639M. doi:10.1017/jfm.2013.343. S2CID 122537381. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:16" class="footnote-back-ref">↩</a></li>
<li id="fn:17">Sagaut, Pierre; Cambon, Claude (2008). Homogeneous Turbulence Dynamics. <a href="#fnref:17" class="footnote-back-ref">↩</a></li>
</ol>

Turbulence modeling open-in-new

Turbulence modeling