# Reynolds stress model (RSM)

### From CFD-Wiki

Littleandy (Talk | contribs) |
(added model variants LRR and SSG with references) |
||

(11 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 [[# | + | 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. |

== Equations == | == 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>\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 exact transport equations for the transport of the Reynolds stresses, <math>\overline{u'_iu'_j}</math> , may be written as follows: | The exact transport equations for the transport of the Reynolds stresses, <math>\overline{u'_iu'_j}</math> , may be written as follows: | ||

+ | <center> | ||

<math> | <math> | ||

- | \frac{\partial}{\partial t}(\rho \overline{u'_iu'_j}) + \frac{\partial}{\partial x_{k}}(\rho u_{k} \overline{u'_iu'_j}) = - \frac{\partial}{\partial x_k}[\rho \overline{u'_iu'_ju'_k} + \overline{p(\delta_{kj}u'_i + \delta_{ik}u'_j)}] </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}[{\mu \frac{\partial}{\partial x_k}(\overline{u'_iu'_j})}] - \rho(\overline{u'_iu'_k}\frac{\partial u_j}{\partial x_k}+\overline{u'_ju'_k}\frac{\partial u_i}{\partial x_k}) - \rho\beta(g_i\overline{u'_j\theta}+g_j\overline{u'_i\theta}) | + | + \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) |

</math> | </math> | ||

- | |||

- | |||

<math> | <math> | ||

- | + \overline{p(\frac{\partial u'_i}{\partial x_j} + \frac{\partial u'_j}{\partial x_i})} - 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(\overline{u'_ju'_m}\epsilon_{ikm} + \overline{u'_iu'_m}\epsilon_{jkm}) + S_{user} | + | -2\rho\Omega_k\left(\overline{u'_ju'_m}\epsilon_{ikm} + \overline{u'_iu'_m}\epsilon_{jkm}\right) + S_{user} |

</math> | </math> | ||

</center> | </center> | ||

Line 34: | Line 29: | ||

<center> | <center> | ||

- | |||

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 | 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> | ||

Line 58: | Line 51: | ||

===Convective Heat and Mass Transfer Modeling=== | ===Convective Heat and Mass Transfer Modeling=== | ||

+ | |||

+ | ===Return-to-isotropy models=== | ||

+ | For an anisotropic turbulence, the Reynolds stress tensor, | ||

+ | <math>\rho\overline{u'_iu'_j}</math> , is usually anisotropic. The | ||

+ | second and third invariances of the Reynolds-stress anisotropic | ||

+ | tensor <math>b_{ij}</math> are nontrivial, where | ||

+ | <math>b_{ij}=\overline{u'_iu'_j}/2k-\delta_{ij} /3 </math> and | ||

+ | <math>k = \overline{u'_iu'_i}/2</math>. 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 <center> | ||

+ | <math> \frac{d b_{ij}}{dt}=-(C_{R}-1) \frac{\varepsilon}{k} | ||

+ | b_{ij}</math> </center>. Here, <math>C_{R}</math> is called the Rotta constant. | ||

== Model constants == | == Model constants == | ||

Line 66: | Line 83: | ||

== 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}} | ||

+ | |||

== Performance, applicability and limitations == | == Performance, applicability and limitations == | ||

== Implementation issues == | == Implementation issues == |

## Revision as of 09:05, 12 August 2010

## 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 [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.

## Equations

The Reynolds stress model involves calculation of the individual Reynolds stresses, , 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, , may be written as follows:

or

Local Time Derivate + = + + + + - + + User-Defined Source Term

where is the Convection-Term, equals the Turbulent Diffusion, stands for the Molecular Diffusion, is the term for Stress Production, equals Buoyancy Production, is for the Pressure Strain, stands for the Dissipation and is the Production by System Rotation.

Of these terms, , , , and do not require modeling. After all, , , , and have to be modeled for closing the equations.

### Modeling Turbulent Diffusive Transport

### 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

For an anisotropic turbulence, the Reynolds stress tensor, , is usually anisotropic. The second and third invariances of the Reynolds-stress anisotropic tensor are nontrivial, where and . 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## Model constants

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

## 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.

## Performance, applicability and limitations

## Implementation issues

## References

**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.