The ReynoldsAveraged NavierStokes (RANS) Equations and Models
Key Takeaways

RANS is a numerical method to model a turbulent flow wherein the flow quantities are decomposed into their timeaveraged and fluctuating components (Reynolds decomposition).

The averaging of NavierStokes equations yields a nonlinear Reynolds stress term that requires additional modeling to fully resolve the system > Turbulence model.

Most (if not all) RANS turbulence models are based on empirical observations.
Turbulent fluid flow can be described with a Reynoldsaveraged NavierStokes (RANS) model.
The NavierStokes equations are arguably among the most complex to use and solve in mathematical physics. This set of coupled timedependent nonlinear equations requires many approximations in order to reach analytical solutions, and not all systems are solvable. For more general systems or complex geometries, numerical techniques are needed to solve these equations and derive useful insights into flow behavior. Some mathematical approximations can be used to reduce the complexity of these equations, which will help speed up simulations without losing too much accuracy.
The Reynoldsaveraged NavierStokes (RANS) equations are a reduced form of the general NavierStokes equations. In the RANS equations, the steadystate solution is decoupled from the timevarying fluctuations in the system, the latter of which will account for turbulence in different flow regimes. Further manipulations of the RANS results and the application of empirical turbulence models yield many other CFD models, some of which have made their way into opensource and commercial CFD applications.
An Overview of the ReynoldsAveraged NavierStokes (RANS) Equations
The technique used to derive the RANS equations is called Reynolds decomposition. Reynolds averaging and Reynolds decomposition do not refer directly to a manipulation of the Reynolds number, but rather to an application of time averaging to the NavierStokes equations. Time averaging is often used to reduce complex systems of differential equations into simpler forms by partially or fully eliminating the time variable. In this case, the RANS equations use a solution that is split into a timeindependent mean flow velocity and timevarying fluctuations about the mean, or .
With this decomposition, one can apply some specialized operators and a time averaging operation to derive the following nonlinear equation describing flow (in tensor notation):
Momentum portion of the incompressible Reynoldsaveraged NavierStokes (RANS) equations
The mean of the strain rate tensor is given by:
Mean strain rate in the Reynoldsaveraged NavierStokes (RANS) equations
From the momentum portion of the RANS results, we can see there is a nonlinear term corresponding to the Reynolds stress, this model is called the Reynolds stress model.
This particular term can account for anisotropic turbulence, as shown in the section on eddy viscosity below. The time dependence of this term generally involves averages of a set of nonlinear derivative terms. The general form of this term is:
Timedependence of the Reynolds stress tensor in RANS
While this model can potentially be very complicated, it also gives some freedom to apply the model to very specific situations where certain terms in the strain and/or stress tensors can be ignored. This requires an additional turbulence model to account for these nonlinear terms in the stress tensor. Most of these turbulence models used in the RANS equations are determined from empirical observations rather than being derived from first principles.
Eddy Viscosity in the ReynoldsAveraged NavierStokes Equations
In order to solve the RANS equations shortly after their development, Boussinesq introduced the concept of eddy viscosity, which is used to reduce the Reynolds stress tensor to a simpler form. Physically, eddy viscosity describes turbulent energy transfer via moving eddy currents in a turbulent flow. Mathematically, this links the Reynolds stress term to the turbulent kinetic energy as defined below:
The constant for the flow gradient terms is the turbulent eddy viscosity (defined above as a kinematic viscosity), which could be a simple constant for viscous turbulent flow. The proportionality constant for the Kronecker delta term is the flow velocity squared term:
Once multiplied by the density in the Reynolds stress equation, these terms define the turbulent kinetic energy as well as the anisotropic turbulence. This treatment is applicable in many systems with free shear flow, such as modeling mixing layers, jets, turbulent boundary layers, channel flows, and many other problems.
OneEquation and TwoEquation Models
There are other reduced models that are a subset of the Reynolds stress model shown above, which may have either one or two equations. These different models apply different definitions relating the dynamic eddy viscosity with the bulk fluid’s viscosity. These other models will also lump some terms in the Reynolds stress equation into compact definitions, as described below:
 SpalartAllmaras model:
 k–ε models:
 k–ω models:
In the above definitions, we have a turbulence dissipation rate defined as , and a specific dissipation rate defined as . Systems designers who need to perform CFD simulations will need to choose the appropriate model that matches their system, something which takes significant initial analysis and observation.
Reduce Computational Complexity With Simulation Tools
Although the resulting equations have a new nonlinear stress term that requires additional modeling, splitting the equations can reduce the computational complexity involved in accounting for turbulent flows. As long as a turbulence model can be developed for the system under investigation, the behavior of a given system can be extended to other fluids or geometries. When faster computation is needed, a oneequation or twoequation model could be used to evaluate fluid flow in a system.
Working with the Reynoldsaveraged NavierStokes equations, or with related models, is much easier when you have an advanced set of CFD simulation tools. The mesh generation applications in Pointwise can help you create accurate simulation grids with high mesh order in complex geometries, helping you reduce computational complexity without sacrificing accuracy. Once a mesh is created and ready for further analysis, the set of CFD simulation tools in Omnis will determine the flow solution in your system accurately and efficiently.
Subscribe to our newsletter for the latest CFD updates or browse Cadence’s suite of CFD software to learn more about how Cadence has the solution for you.