Welcome to AE Resources
Converted document Converted document
ALGEBRAIC MODELS
Algebraic models are the simplest form of turbulence models. The eddy viscosity is described by a basic algebraic relation. This approach was first introduced by Prandtl (1925) with his mixing-length hypothesis. In analogy to the molecular momentum process, Prandtl postulated that for a shear flow moving in the x direction,
(1) μt = ρl2mix||(dU)/(dy)||
where the mixing length lmix depends on the nature of the flow and is typically space dependent. The mixing length model is said to be incomplete. A prior knowledge of the flow is required to obtain reliable results. For example, experimental data suggested that the mixing length of a free shear flow is proportional to the half-width b of the layer.
(2) lm = 0.180b  for a far wake [10];
lm = 0.080b  for a plane jet [2]; lm = 0.071b  for a mixing layer [8].
Prandtl postulated that the mixing length of a wall bounded flow would be proportional to the distance from the surface. This assumption is consistent with the well known law of the wall, which describes the velocity profile into the log layer. In this region, the flow is sufficiently close to the wall so that the inertial terms can be neglected, but still far enough to neglect the viscous stresses compared to the Reynolds stresses. In a two-dimensional incompressible boundary layer, the Navier-Stokes equations can be written
(3) U(U)/(x) + V(U)/(y) = (ν + νT)(2U)/(y2)
Neglecting the convective terms, the sum of viscous and Reynolds stresses must therefore be constant in the flow completely to the surface, where
(4) (μ + μT)(U)/(y) ~ μ(U)/(y)w = τw = ρu2τ
uτ = (τw ⁄ ρ) is called the friction velocity. It is frequently used to non-dimensionalize the velocity and distance from the wall, through the equations
(5) u +  = U ⁄ uτ y +  = uτy ⁄ ν
where y is the distance from the wall and ν is the molecular kinematic viscosity. Introducing Prandtl’s eddy viscosity (lmix = κy, with κ constant) into Eq.(5↑) and neglecting the viscous stresses, the following equation is obtained.
(6) μT(U)/(y) = ρl2mix(U)/(y)2 = ρκ2y2(U)/(y)2 = ρu2τ
which can be integrated to yield
(7) U = (uτ)/(κ)lny + constant or u +  = (1)/(κ)lny +  + C
figure images/lawofthewall.png
Figure 1 Velocity profile in a typical turbulent boundary layer. Source: Kline (1967) [7]
Next Page →
Next Page →