Turbulence modeling

From Wikipedia, the free encyclopedia

Turbulence modeling is the area of physical modeling where a simpler mathematical model than the full time dependent Navier-Stokes Equations is used to predict the effects of turbulence. There are various mathematical models used in flow modelling to understand turbulence.

Joseph Boussinesq was the first practitioner of this, introducing the concept of eddy viscosity. In this model, the additional turbulent stresses are given by augmenting the molecular viscosity with an eddy viscosity. This can be a simple constant eddy viscosity (which works well for some free shear flows such as axisymmetric jets, 2-D jets, and mixing layers). Later, Ludwig Prandtl introduced the additional concept of the mixing length, along with the idea of a boundary layer. For wall-bounded turbulent flows, the eddy viscosity must vary with distance from the wall, hence the addition of the concept of a 'mixing length'. In the simplest wall-bounded flow model, the eddy viscosity is given by the equation:

\nu_t = \kappa \left|\frac{\partial u}{\partial y}\right|l^2
where:
κ is the von Karman constant (0.41);
\frac{\partial u}{\partial y} is the partial derivative of the streamwise velocity (u) with respect to the wall normal direction (y);
l is the distance from the wall.

This simple model is the basis for the "Law of the Wall", which is a surprisingly accurate model for wall-bounded, attached (not separated) flow fields with small pressure gradients.

More general turbulence models have evolved over time, with most modern turbulence models given by field equations similar to the Navier-Stokes equations.

Among many others, Joseph Smagorinsky (1964) proposed a useful formula for the eddy viscosity in numerical models, based on the local derivatives of the velocity field and the local grid size:

\nu_t = \Delta x \Delta y \sqrt{\left(\frac{\partial u}{\partial x}\right)^2 + \left(\frac{\partial v}{\partial y}\right)^2 + \frac{1}{2}\left(\frac{\partial u}{\partial y} + \frac{\partial v}{\partial x}\right)^2}

[edit] References