CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Reynolds stress model (RSM)

Reynolds stress model (RSM)

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
(Return-to-isotropy models)
m (References)
 
(31 intermediate revisions not shown)
Line 1: Line 1:
== Introduction ==
== Introduction ==
-
The Reynold's stress model (RSM) is a higher level, elaborate turbulence model. It is usually called a ''Second Order Closure''. This modelling approach originates from the work by [[#References|[Launder (1975)]]]. In RSM, the eddy viscosity approach has been discarded and the Reynolds stresses are directly computed. The exact Reynolds stress transport equation accounts for the directional effects of the Reynolds stress fields.
+
The Reynolds Stress Models (RSM), also known as the Reynolds Stress Transport Models, are higher level turbulence closures and represent the most complete classical turbulence model. The method of closure employed is usually called a ''Second Order Closure''. This modelling approach originates from the works by [[#References|Chou (1945)]] and [[#References|Rotta (1951)]]. In Reynolds Stress Models, the eddy viscosity approach is avoided and the individual components of the Reynolds stress tensor are directly computed. These models rely on the exact Reynolds stress transport equation. They are able to account for complex interactions in turbulent flow fields, such as the directional effects of the Reynolds stresses.
-
== Equations ==
+
== Shortcomings of Eddy-viscosity based models ==
 +
Eddy-viscosity based models like the  <math> k-\epsilon </math> and the <math>k - \omega</math> models have significant shortcomings in complex, real-life turbulent flows that are often encountered in engineering applications. In flows with high degrees of anisotropy, significant streamline curvature, flow separation, flows with zones of re-circulation or flows influenced by mean rotational effects, the performance of such models is unsatisfactory ([[#References|Craft (1996)]]).
 +
 +
Eddy viscosity based closures cannot account for the limiting states of turbulence. In decaying turbulent flows, it has been observed that the state of turbulence rapidly approaches an isotropic state with equipartition of the turbulent kinetic energy among the components of the Reynolds stresses. Eddy-viscosity based models can never replicate this return to isotropy behaviour. Eddy-viscosity based models cannot replicate the behaviour of turbulent flows in the Rapid Distortion limit, where the turbulent flow essentially behaves as an elastic medium (instead of viscous).
 +
 +
Reynolds Stress Models have an clear advantage over lower-order models for such flows where the transport of terms like the Reynolds stress or heat flux, play important roles. This is because the physics of the transport of second-moment terms can be built into the second-order models.
 +
 +
== Equations ==
The Reynolds stress model involves calculation of the individual Reynolds stresses, <math>\rho\overline{u'_iu'_j}</math> , using differential transport equations. The individual Reynolds stresses are then used to obtain closure of the Reynolds-averaged momentum equation.
The Reynolds stress model involves calculation of the individual Reynolds stresses, <math>\rho\overline{u'_iu'_j}</math> , using differential transport equations. The individual Reynolds stresses are then used to obtain closure of the Reynolds-averaged momentum equation.
Line 12: Line 19:
<center>
<center>
<math>
<math>
-
\frac{\partial}{\partial t}\left(\rho \overline{u'_iu'_j}\right) + \frac{\partial}{\partial x_{k}}\left(\rho u_{k} \overline{u'_iu'_j}\right) = - \frac{\partial}{\partial x_k}\left[\rho \overline{u'_iu'_ju'_k} + \overline{p\left(\delta_{kj}u'_i + \delta_{ik}u'_j\right)}\right] </math>
+
\frac{\partial}{\partial t}\left(\rho \overline{u'_iu'_j}\right) + \frac{\partial}{\partial x_{k}}\left(\rho u_{k} \overline{u'_iu'_j}\right) = - \frac{\partial}{\partial x_k}\left[\rho \overline{u'_iu'_ju'_k} + \overline{p'\left(\delta_{kj}u'_i + \delta_{ik}u'_j\right)}\right] </math>
-
 
+
<math>
<math>
-
  + \frac{\partial}{\partial x_k}\left[{\mu \frac{\partial}{\partial x_k}\left(\overline{u'_iu'_j}\right)}\right] - \rho\left(\overline{u'_iu'_k}\frac{\partial u_j}{\partial x_k}+\overline{u'_ju'_k}\frac{\partial u_i}{\partial x_k}\right) - \rho\beta\left(g_i\overline{u'_j\theta}+g_j\overline{u'_i\theta}\right)
+
  + \frac{\partial}{\partial x_k}\left[{\mu \frac{\partial}{\partial x_k}\left(\overline{u'_iu'_j}\right)}\right] - \rho\left(\overline{u'_iu'_k}\frac{\partial u_j}{\partial x_k}+\overline{u'_ju'_k}\frac{\partial u_i}{\partial x_k}\right)  
</math>
</math>
-
 
-
 
<math>
<math>
-
+ \overline{p\left(\frac{\partial u'_i}{\partial x_j} + \frac{\partial u'_j}{\partial x_i}\right)} - 2\mu\overline{\frac{\partial u'_i}{\partial x_k} \frac{\partial u'_j}{\partial x_k}}
+
+ \overline{p'\left(\frac{\partial u'_i}{\partial x_j} + \frac{\partial u'_j}{\partial x_i}\right)} - 2\mu\overline{\frac{\partial u'_i}{\partial x_k} \frac{\partial u'_j}{\partial x_k}}
</math>
</math>
-
 
-
 
<math>
<math>
-
-2\rho\Omega_k\left(\overline{u'_ju'_m}\epsilon_{ikm} + \overline{u'_iu'_m}\epsilon_{jkm}\right) + S_{user}
+
-2\rho\Omega_k\left(\overline{u'_ju'_m}\epsilon_{ikm} + \overline{u'_iu'_m}\epsilon_{jkm}\right)  
</math>
</math>
</center>
</center>
Line 35: Line 37:
<center>
<center>
-
 
+
Local Time Derivate + <math>C_{ij}</math> = <math>D_{T,ij}</math> + <math>D_{L,ij}</math> + <math>P_{ij}</math> + <math>\phi_{ij}</math> - <math>\epsilon_{ij}</math> + <math>F_{ij}</math>  
-
Local Time Derivate + <math>C_{ij}</math> = <math>D_{T,ij}</math> + <math>D_{L,ij}</math> + <math>P_{ij}</math> + <math>G_{ij}</math> + <math>\phi_{ij}</math> - <math>\epsilon_{ij}</math> + <math>F_{ij}</math> + User-Defined Source Term
+
-
 
+
</center>
</center>
-
where <math>C_{ij}</math> is the [[Convection]]-Term, <math>D_{T,ij}</math> equals the [[Turbulent Diffusion]], <math>D_{L,ij}</math> stands for the [[Molecular Diffusion]], <math>P_{ij}</math> is the term for [[Stress]] Production, <math>G_{ij}</math> equals [[Buoyancy]] Production, <math>\phi_{ij}</math> is for the [[Pressure Strain]], <math>\epsilon_{ij}</math> stands for the Dissipation and <math>F_{ij}</math> is the Production by System Rotation.
+
where <math>C_{ij}</math> is the [[Convection]]-Term, <math>D_{T,ij}</math> equals the [[Turbulent Diffusion]], <math>D_{L,ij}</math> stands for the [[Molecular Diffusion]], <math>P_{ij}</math> is the term for [[Stress]] Production, <math>\phi_{ij}</math> is for the [[Pressure Strain]], <math>\epsilon_{ij}</math> stands for the Dissipation and <math>F_{ij}</math> is the Production by System Rotation.
-
Of these terms, <math>C_{ij}</math>, <math>D_{L,ij}</math>, <math>P_{ij}</math>, and <math>F_{ij}</math> do not require modeling. After all, <math>D_{T,ij}</math>, <math>G_{ij}</math>, <math>\phi_{ij}</math>, and <math>\epsilon_{ij}</math> have to be modeled for closing the equations.
+
Of these terms, <math>C_{ij}</math>, <math>D_{L,ij}</math>, <math>P_{ij}</math>, and <math>F_{ij}</math> do not require modelling. However, <math>D_{T,ij}</math>, <math>\phi_{ij}</math>, and <math>\epsilon_{ij}</math> have to be modelled for closing the equations. The fidelity of the Reynolds stress model depends on the accuracy of the models for the turbulent transport, the pressure-strain correlation and the dissipation terms.
-
===Modeling Turbulent Diffusive Transport===
 
===Modeling the Pressure-Strain Term===
===Modeling the Pressure-Strain Term===
-
 
-
===Effects of Buoyancy on Turbulence===
 
-
 
-
===Modeling the Turbulence Kinetic Energy===
 
-
 
-
===Modeling the Dissipation Rate===
 
-
 
-
===Modeling the Turbulent Viscosity===
 
-
 
-
===Boundary Conditions for the Reynolds Stresses===
 
-
 
-
===Convective Heat and Mass Transfer Modeling===
 
===Return-to-isotropy models===
===Return-to-isotropy models===
Line 81: Line 68:
negligible. The turbulence dissipation and slow pressure-strain-rate
negligible. The turbulence dissipation and slow pressure-strain-rate
are preponderant. Under these cirsumstance, Rotta suggested <center>
are preponderant. Under these cirsumstance, Rotta suggested <center>
-
<math> \frac{d b_{ij}}/{dt}=-(C_{R}-1) \frac{\varepsilon}{k}
+
<math> \frac{d b_{ij}}{dt}=-(C_{R}-1) \frac{\varepsilon}{k}
-
b_{ij}</math> </center>. Here, <math>C_{R}</math> is called Rotta constant.
+
b_{ij}</math> </center>. Here, <math>C_{R}</math> is called the Rotta constant.
 +
 
 +
===Modeling Turbulent Diffusive Transport===
 +
 
 +
===Modeling the Dissipation Rate===
== Model constants ==
== Model constants ==
Line 91: Line 82:
== Model variants ==
== Model variants ==
 +
 +
===LRR, Launder-Reece-Rodi===
 +
 +
{{reference-paper|year=1975|author=Launder, B. E., Reece, G. J. and Rodi, W.|title=Progress in the Development of a Reynolds-Stress Turbulent Closure.|rest=Journal of Fluid Mechanics, Vol. 68(3), pp. 537-566}}
 +
 +
===SSG, Speziale-Sarkar-Gatski===
 +
 +
{{reference-paper|year=1991|author=Speziale, C.G., Sarkar, S., Gatski, T.B.|title=Modeling the Pressure-Strain Correlation of Turbulence: an Invariant Dynamical Systems Approach|rest=Journal of Fluid Mechanics, Vol. 227, pp. 245-272}}
 +
 +
===Johansson & Hallbäck===
 +
 +
{{reference-paper|year=1994|author=Johansson, A. V. and Hallbäck, M. |title=Modelling of rapid pressure—strain in Reynolds-stress closures|rest=Journal of Fluid Mechanics, Vol. 269, pp. 143-168}}
 +
 +
===Mishra & Girimaji===
 +
 +
{{reference-paper|year=2017|author=Mishra, A. A. and Girimaji, S. S.|title=Toward approximating non-local dynamics in single-point pressure–strain correlation closures|rest=Journal of Fluid Mechanics, Vol. 811, pp. 168-188}}
 +
== Performance, applicability and limitations ==
== Performance, applicability and limitations ==
== Implementation issues ==
== Implementation issues ==
== References ==
== References ==
 +
 +
{{reference-book|author=Kemal Hanjalić and Brian Launder|year=2011|title=Modelling Turbulence in Engineering and the Environment: Second-Moment Routes to Closure|rest=ISBN 978-0521883986, Cambridge University Press}}
 +
 +
{{reference-paper|year=1945|author=Chou, P.|title=On velocity correlations and the solutions of the equations of turbulent fluctuation.|rest=Quarterly of Applied Mathematics, Vol. 3(1), pp. 38-54}}
 +
 +
{{reference-paper|year=1951|author=Rotta, J.|title=Statistische theorie nichthomogener turbulenz.|rest=Zeitschrift für Physik A, Vol. 129(6), pp. 547-572}}
{{reference-paper|year=1975|author=Launder, B. E., Reece, G. J. and Rodi, W.|title=Progress in the Development of a Reynolds-Stress Turbulent Closure.|rest=Journal of Fluid Mechanics, Vol. 68(3), pp. 537-566}}
{{reference-paper|year=1975|author=Launder, B. E., Reece, G. J. and Rodi, W.|title=Progress in the Development of a Reynolds-Stress Turbulent Closure.|rest=Journal of Fluid Mechanics, Vol. 68(3), pp. 537-566}}
 +
 +
{{reference-paper|year=1996|author=Craft, T., Launder, B. and K. Suga|title=Development and application of a cubic eddy-viscosity model of turbulence.|rest=International Journal of Heat and Fluid Flow, Vol. 17, pp. 108-115}}
{{Stub}}
{{Stub}}
[[Category:Turbulence models]]
[[Category:Turbulence models]]

Latest revision as of 01:47, 18 April 2017

Contents

Introduction

The Reynolds Stress Models (RSM), also known as the Reynolds Stress Transport Models, are higher level turbulence closures and represent the most complete classical turbulence model. The method of closure employed is usually called a Second Order Closure. This modelling approach originates from the works by Chou (1945) and Rotta (1951). In Reynolds Stress Models, the eddy viscosity approach is avoided and the individual components of the Reynolds stress tensor are directly computed. These models rely on the exact Reynolds stress transport equation. They are able to account for complex interactions in turbulent flow fields, such as the directional effects of the Reynolds stresses.

Shortcomings of Eddy-viscosity based models

Eddy-viscosity based models like the  k-\epsilon and the k - \omega models have significant shortcomings in complex, real-life turbulent flows that are often encountered in engineering applications. In flows with high degrees of anisotropy, significant streamline curvature, flow separation, flows with zones of re-circulation or flows influenced by mean rotational effects, the performance of such models is unsatisfactory (Craft (1996)).

Eddy viscosity based closures cannot account for the limiting states of turbulence. In decaying turbulent flows, it has been observed that the state of turbulence rapidly approaches an isotropic state with equipartition of the turbulent kinetic energy among the components of the Reynolds stresses. Eddy-viscosity based models can never replicate this return to isotropy behaviour. Eddy-viscosity based models cannot replicate the behaviour of turbulent flows in the Rapid Distortion limit, where the turbulent flow essentially behaves as an elastic medium (instead of viscous).

Reynolds Stress Models have an clear advantage over lower-order models for such flows where the transport of terms like the Reynolds stress or heat flux, play important roles. This is because the physics of the transport of second-moment terms can be built into the second-order models.

Equations

The Reynolds stress model involves calculation of the individual Reynolds stresses, \rho\overline{u'_iu'_j} , using differential transport equations. The individual Reynolds stresses are then used to obtain closure of the Reynolds-averaged momentum equation.

The exact transport equations for the transport of the Reynolds stresses, \overline{u'_iu'_j} , may be written as follows:


\frac{\partial}{\partial t}\left(\rho \overline{u'_iu'_j}\right) + \frac{\partial}{\partial x_{k}}\left(\rho u_{k} \overline{u'_iu'_j}\right) = - \frac{\partial}{\partial x_k}\left[\rho \overline{u'_iu'_ju'_k} + \overline{p'\left(\delta_{kj}u'_i + \delta_{ik}u'_j\right)}\right]


 + \frac{\partial}{\partial x_k}\left[{\mu \frac{\partial}{\partial x_k}\left(\overline{u'_iu'_j}\right)}\right] - \rho\left(\overline{u'_iu'_k}\frac{\partial u_j}{\partial x_k}+\overline{u'_ju'_k}\frac{\partial u_i}{\partial x_k}\right)


+ \overline{p'\left(\frac{\partial u'_i}{\partial x_j} + \frac{\partial u'_j}{\partial x_i}\right)} - 2\mu\overline{\frac{\partial u'_i}{\partial x_k} \frac{\partial u'_j}{\partial x_k}}


-2\rho\Omega_k\left(\overline{u'_ju'_m}\epsilon_{ikm} + \overline{u'_iu'_m}\epsilon_{jkm}\right)

or

Local Time Derivate + C_{ij} = D_{T,ij} + D_{L,ij} + P_{ij} + \phi_{ij} - \epsilon_{ij} + F_{ij}

where C_{ij} is the Convection-Term, D_{T,ij} equals the Turbulent Diffusion, D_{L,ij} stands for the Molecular Diffusion, P_{ij} is the term for Stress Production, \phi_{ij} is for the Pressure Strain, \epsilon_{ij} stands for the Dissipation and F_{ij} is the Production by System Rotation.

Of these terms, C_{ij}, D_{L,ij}, P_{ij}, and F_{ij} do not require modelling. However, D_{T,ij}, \phi_{ij}, and \epsilon_{ij} have to be modelled for closing the equations. The fidelity of the Reynolds stress model depends on the accuracy of the models for the turbulent transport, the pressure-strain correlation and the dissipation terms.


Modeling the Pressure-Strain Term

Return-to-isotropy models

For an anisotropic turbulence, the Reynolds stress tensor, \rho\overline{u'_iu'_j} , is usually anisotropic. The second and third invariances of the Reynolds-stress anisotropic tensor b_{ij} are nontrivial, where b_{ij}=\overline{u'_iu'_j}/2k-\delta_{ij} /3 and k = \overline{u'_iu'_i}/2. It is naturally to suppose that the anisotropy of the Reynolds-stress tensor results from the anisotropy of turbulent production, dissipation, transport, pressure-stain-rate, and the viscous diffusive tensors. The Reynolds-stress tensor returns to isotropy when the anisotropy of these turbulent components return to isotropy. Such a correlation is described by the Reynolds stress transport equation. Based on these consideration, a number of turbulent models, such as Rotta's model and Lumley's return-to-isotropy model, have been established.

Rotta's model describes the linear return-to-isotropy behavior of a low Reynolds number homogenous turbulence in which the turbulent production, transport, and rapid pressure-strain-rate are negligible. The turbulence dissipation and slow pressure-strain-rate

are preponderant. Under these cirsumstance, Rotta suggested
 \frac{d b_{ij}}{dt}=-(C_{R}-1) \frac{\varepsilon}{k}
b_{ij}
. Here, C_{R} is called the Rotta constant.

Modeling Turbulent Diffusive Transport

Modeling the Dissipation Rate

Model constants

The constants suggested for use in this model are as follows:

C_s \approx 0.25, C_l \approx 0.25, C_\gamma \approx 0.25

Model variants

LRR, Launder-Reece-Rodi

Launder, B. E., Reece, G. J. and Rodi, W. (1975), "Progress in the Development of a Reynolds-Stress Turbulent Closure.", Journal of Fluid Mechanics, Vol. 68(3), pp. 537-566.

SSG, Speziale-Sarkar-Gatski

Speziale, C.G., Sarkar, S., Gatski, T.B. (1991), "Modeling the Pressure-Strain Correlation of Turbulence: an Invariant Dynamical Systems Approach", Journal of Fluid Mechanics, Vol. 227, pp. 245-272.

Johansson & Hallbäck

Johansson, A. V. and Hallbäck, M. (1994), "Modelling of rapid pressure—strain in Reynolds-stress closures", Journal of Fluid Mechanics, Vol. 269, pp. 143-168.

Mishra & Girimaji

Mishra, A. A. and Girimaji, S. S. (2017), "Toward approximating non-local dynamics in single-point pressure–strain correlation closures", Journal of Fluid Mechanics, Vol. 811, pp. 168-188.

Performance, applicability and limitations

Implementation issues

References

Kemal Hanjalić and Brian Launder (2011), Modelling Turbulence in Engineering and the Environment: Second-Moment Routes to Closure, ISBN 978-0521883986, Cambridge University Press.

Chou, P. (1945), "On velocity correlations and the solutions of the equations of turbulent fluctuation.", Quarterly of Applied Mathematics, Vol. 3(1), pp. 38-54.

Rotta, J. (1951), "Statistische theorie nichthomogener turbulenz.", Zeitschrift für Physik A, Vol. 129(6), pp. 547-572.

Launder, B. E., Reece, G. J. and Rodi, W. (1975), "Progress in the Development of a Reynolds-Stress Turbulent Closure.", Journal of Fluid Mechanics, Vol. 68(3), pp. 537-566.

Craft, T., Launder, B. and K. Suga (1996), "Development and application of a cubic eddy-viscosity model of turbulence.", International Journal of Heat and Fluid Flow, Vol. 17, pp. 108-115.


My wiki