Reynolds Stress Model (RSM)
Reynolds Stress Model (RSM): Theoretical Foundations
Overview
Professor, what's the difference between the Reynolds Stress Model (RSM) and other RANS models?
The biggest difference is that it does not use the eddy viscosity hypothesis (Boussinesq assumption). Models like k-epsilon or k-omega relate the Reynolds stress $\overline{u_i'u_j'}$ to the mean strain rate via the eddy viscosity $\mu_t$. RSM does not employ this assumption and instead directly solves transport equations for each of the six components of the Reynolds stress tensor.
Six components? Does that mean six more equations?
The Reynolds stress tensor is symmetric, so it has six independent components. Adding the equation for $\varepsilon$ (Dissipation Rate) gives a total of seven additional equations to solve. Compared to two-equation models (like k-omega), the computational cost increases by 2 to 3 times.
Governing Equations
Could you show me the specific equations?
The transport equation for the Reynolds stress $R_{ij} = \overline{u_i'u_j'}$ is as follows.
Let's summarize the meaning of each term.
| Term | Expression | Physical Meaning |
|---|---|---|
| Production $P_{ij}$ | $-R_{ik}\frac{\partial U_j}{\partial x_k} - R_{jk}\frac{\partial U_i}{\partial x_k}$ | Production by mean velocity gradient (can be calculated exactly) |
| Pressure-Strain $\Pi_{ij}$ | $\overline{p'\left(\frac{\partial u_i'}{\partial x_j}+\frac{\partial u_j'}{\partial x_i}\right)}$ | Redistribution of energy among components (requires modeling) |
| Turbulent Diffusion $D_{ij}^T$ | $-\frac{\partial}{\partial x_k}\overline{u_i'u_j'u_k'}$ | Transport by turbulence (requires modeling) |
| Viscous Diffusion $D_{ij}^\nu$ | $\nu\nabla^2 R_{ij}$ | Diffusion by molecular viscosity (exact) |
| Dissipation $\varepsilon_{ij}$ | $2\nu\overline{\frac{\partial u_i'}{\partial x_k}\frac{\partial u_j'}{\partial x_k}}$ | Viscous dissipation (requires modeling) |
So not everything can be calculated exactly.
Exactly. The production term $P_{ij}$ and viscous diffusion can be calculated in closed form, but the pressure-strain term, turbulent diffusion term, and dissipation tensor require modeling. The accuracy of RSM heavily depends on the quality of these models.
Pressure-Strain Term Models
What kinds of models are there for the pressure-strain term?
Let me list some representative models.
| Model | Proposer | Features |
|---|---|---|
| LRR (Linear Return to Isotropy) | Launder, Reece, Rodi (1975) | Linear model. Standard for industrial applications. |
| SSG (Speziale-Sarkar-Gatski) | Speziale et al. (1991) | Quadratic nonlinear model. Improved accuracy. |
| GL (Gibson-Launder) | Gibson, Launder (1978) | Accounts for wall reflection effects. |
In the LRR model, $\Pi_{ij}$ is decomposed as follows.
$\Pi_{ij,1}$ is the slow return term (effect of returning anisotropy to isotropy), $\Pi_{ij,2}$ is the rapid term (redistribution by mean strain), and $\Pi_{ij,w}$ is the wall reflection term.
The Luxury of Solving 7 Equations β Why RSM is the "Pinnacle of Theory"
The Reynolds Stress Model (RSM) solves individual transport equations for the six stress components $\overline{u_i u_j}$ and the dissipation rate Ξ΅, making it the most rigorous model within RANS. By not using the eddy viscosity hypothesis, it can naturally represent turbulence anisotropy caused by swirling flows, curvature, and buoyancy. The trade-off is computational cost: it requires 2-3 times more memory and computation time than two-equation models, and convergence is difficult. While there are real-world scenarios like secondary flows in turbomachinery or cyclone separator design where "RSM is absolutely necessary," it also serves as a prime example that "a more accurate model is not necessarily the model you should always use."
Computational Methods for Reynolds Stress Model (RSM)
Discretization and Solution in FVM
How do you solve the seven equations of RSM?
In a segregated solver, the six Reynolds stress equations and one $\varepsilon$ equation are solved iteratively in sequence. An RSM update step is inserted between the momentum and pressure equations.
1. Solve the momentum equations (using the stress tensor directly from $R_{ij}$, not via eddy viscosity)
2. Solve the pressure correction equation (SIMPLE/PISO)
3. Solve the six Reynolds stress equations
4. Solve the $\varepsilon$ equation
5. Convergence Check, iterate
What do you mean by using Reynolds stress directly in the momentum equations?
In the k-epsilon model, the eddy viscosity hypothesis $-\overline{u_i'u_j'} = \mu_t(\partial U_i/\partial x_j + \partial U_j/\partial x_i) - (2/3)k\delta_{ij}$ treats it as a diffusion-like term. In RSM, $\overline{u_i'u_j'}$ is inserted directly as a source term into the momentum equations. This allows anisotropic effects to be accurately reflected.
Numerical Stability Challenges
I've heard RSM is difficult to converge.
Because it solves seven strongly coupled nonlinear equations, convergence is significantly more difficult compared to k-epsilon. Let me list the main problems and countermeasures.
| Problem | Cause | Countermeasure |
|---|---|---|
| Reynolds Stress Realizability Violation | $R_{ij}$ ceases to be positive definite symmetric | Apply Realizability constraint |
| Overestimation of $\varepsilon$ | Excessive dissipation near walls | Use $\omega$-based RSM (BSL-RSM) |
| Slow Convergence | Strong coupling between equations | Set Under-relaxation to 0.3β0.5 |
| Sensitivity to Initial Values | Initial $R_{ij}$ is non-physical | Initialize with k-epsilon results |
As a practical TIP, when starting an RSM calculation, the golden rule is to first obtain a converged solution with k-epsilon or SST k-omega, then use that result to initialize the RSM.
Omega-based RSM
Are there RSM variants that use $\omega$ instead of $\varepsilon$?
Yes. The BSL-RSM (Baseline RSM) proposed by Menter uses the $\omega$ equation from SST k-omega instead of the $\varepsilon$ equation. This improves near-wall behavior and works well with wall functions.
| RSM Variant | Dissipation Equation | Wall Treatment | Solver Support |
|---|---|---|---|
| LRR-RSM | $\varepsilon$ | Low-Re or Wall Function | Fluent, CFX, OpenFOAM |
| SSG-RSM | $\varepsilon$ | Low-Re | Fluent, CFX, STAR-CCM+ |
| BSL-RSM ($\omega$-based) | $\omega$ | Automatic WT | CFX, Fluent |
| LRR-RSM-w | $\omega$ | Wall Function Compatible | OpenFOAM |
Is BSL-RSM recommended in CFX?
In CFX, BSL-RSM is the default RSM. It can use the same wall treatment as SST k-omega (Automatic Wall Treatment), making it robust with respect to mesh $y^+$.
Scenes Where RSM Achieves a Come-from-Behind Victory in Cyclone Separator Design
In the analysis of cyclone separators (devices that separate particles or droplets via centrifugal force), it is known that standard k-Ξ΅ or SST models can significantly mispredict separation efficiency. This is because strong swirling flows cause severe anisotropy in turbulent stresses, breaking down the eddy viscosity hypothesis. When RSM is applied to this problem, many cases report that swirl velocity distribution and pressure drop align much more closely with experimental values. This is why some companies in chemical plant cyclone design use RSM as standard. "For specific flows, only RSM works"βthis reality justifies the high cost.
Practical CAE quality notes for Reynolds Stress Model (RSM)
Reynolds Stress Model (RSM) should be treated as an engineering model, not as an isolated formula. In fluid simulation, reliable results come from a clear chain of assumptions: governing physics, material data, boundary conditions, numerical discretization, solver settings, and post-processing criteria. Before using this note in a design review, identify which quantities are prescribed, which are solved, and which are only diagnostic indicators.
Model setup checklist
- Define the scope: decide whether Reynolds Stress Model (RSM) is being used for screening, detailed design, failure investigation, or verification of another simulation.
- Check dimensions and units: keep SI units internally and document every conversion applied to loads, geometry, material constants, and time or frequency scales.
- State assumptions explicitly: record linearity, steady-state or transient behavior, small-deformation limits, continuum assumptions, and any symmetry or ideal boundary conditions.
- Compare with a baseline: use a hand calculation, limiting case, mesh refinement trend, or independent solver result before accepting the final value.
Validation signals
| Review item | What to verify | Typical warning sign |
|---|---|---|
| Inputs | Geometry, material data, loads, and constraints match the intended fluid simulation problem. | Correct-looking plots with unrealistic magnitudes or units. |
| Numerics | Mesh, time step, convergence tolerance, and solver options are adequate for Reynolds Stress Model. | Large changes after small mesh or tolerance adjustments. |
| Physics | The selected theory remains valid in the expected stress, temperature, velocity, or frequency range. | Results are used outside the assumptions stated in the model. |
For production use, keep the model file, input table, result plots, and review comments together. This makes Reynolds Stress Model (RSM) traceable and prevents the page from being used as a black-box answer without engineering judgment.