Skip to main content

RANS equations and the closure problem

Short answer

Flow variables are split into mean and fluctuating parts; time-averaging the incompressible N-S equations yields RANS equations with Reynolds stresses -ρ⟨u'_i u'_j⟩. More unknowns than equations cause a closure problem. The Boussinesq hypothesis models the stresses via eddy viscosity: -ρ⟨u'_i u'_j⟩ = μ_t(∂ū_i/∂x_j+∂ū_j/∂x_i) - (2/3)ρk δ_{ij}.

Explain Reynolds decomposition and time averaging. Derive the Reynolds-averaged continuity and momentum equations for incompressible flow, and explain the closure problem and the Boussinesq hypothesis.

Reynolds decomposition and averaging

Each flow variable is split into a mean and a fluctuation:

u=uˉ+u′,p=pˉ+p′,uˉ=1T∫0Tu dt,u′‾=0u=\bar u+u',\qquad p=\bar p+p',\qquad \bar u=\frac1T\int_0^Tu\,dt,\qquad \overline{u'}=0

Rules: uˉ‾=uˉ\overline{\bar u}=\bar u, u′‾=0\overline{u'}=0, but u′v′‾≠0\overline{u'v'}\neq0 in general. For unsteady mean flow, ensemble averaging is used.

RANS equations

Substitute into the incompressible continuity and Navier-Stokes equations and average:

Continuity:

∂uˉi∂xi=0\frac{\partial\bar u_i}{\partial x_i}=0

Momentum:

ρ(∂uˉi∂t+uˉj∂uˉi∂xj)=−∂pˉ∂xi+∂∂xj(μ∂uˉi∂xj−ρui′uj′‾)\rho\left(\frac{\partial\bar u_i}{\partial t}+\bar u_j\frac{\partial\bar u_i}{\partial x_j}\right)=-\frac{\partial\bar p}{\partial x_i}+\frac{\partial}{\partial x_j}\left(\mu\frac{\partial\bar u_i}{\partial x_j}-\rho\overline{u'_iu'_j}\right)

The new term −ρui′uj′‾-\rho\overline{u'_iu'_j} is the Reynolds stress tensor: 6 independent components (symmetric) which represent the extra momentum transport by turbulence. The term comes from the averaging of the non-linear convective term.

Closure problem

The averaged equations have 4 equations (continuity + 3 momentum) but 10 unknowns (4 mean variables + 6 Reynolds stresses). Writing equations for the stresses brings in triple correlations u′u′u′‾\overline{u'u'u'}, and so on: every level introduces more unknowns than equations. The set cannot be closed without turbulence models that relate the Reynolds stresses to mean quantities.

Boussinesq hypothesis

Reynolds stresses are assumed to act like viscous stresses, proportional to the mean strain rate, through an eddy viscosity μt\mu_t:

−ρui′uj′‾=μt(∂uˉi∂xj+∂uˉj∂xi)−23ρkδij,k=12ui′ui′‾-\rho\overline{u'_iu'_j}=\mu_t\left(\frac{\partial\bar u_i}{\partial x_j}+\frac{\partial\bar u_j}{\partial x_i}\right)-\frac23\rho k\delta_{ij},\qquad k=\tfrac12\overline{u'_iu'_i}

The problem then becomes one of finding μt\mu_t: zero-equation (mixing length), one-equation (Spalart-Allmaras) and two-equation (k-ε\varepsilon, k-ω\omega) models. Alternatively, Reynolds stress models (RSM) solve transport equations for each stress, without the isotropy assumption.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Applied Computational Fluid Dynamics chapter 5 questions →