Boundary layer theory

From Thermal-FluidsPedia

Revision as of 07:26, 21 March 2010 by Founding Editor (Talk | contribs)
Jump to: navigation, search


Concepts of the Boundary Layer Theory

Unless the geometry is very simple, it is difficult to solve for the complete viscous fluid flow around a body. A full domain numerical solution is time consuming and impractical, because one needs to solve the full Navier-Stokes equations in the full domain, which are nonlinear, elliptic, and complex.

In 1904, Prandtl discovered that for most practical applications, the influence of viscosity is observed in a very thin domain close to the object, as shown in Fig.1. Outside this region one can assume the flow is inviscid (µ = 0) .

Viscous or momentum boundary layer
Figure 1: Viscous or momentum boundary layer.

The thin region where the effect of viscosity is dominant is called the momentum or viscous boundary layer. The solution of boundary layer analysis can be simplified due to the fact that its thickness is much smaller than the characteristic dimension of the object. The fluid adjacent to the surface of the body has zero relative velocity; ufluid – usurface = 0. This is also called the no slip boundary condition.

One of the assumptions of boundary layer approximations is that the fluid velocity at the wall is at rest relative to the surface. This is true except when the fluid pressure is very low, and therefore, the Kundsen number, Kn= λ/L, of the fluid molecules is much larger than 1. For external flow, the flow next to an object can be divided into two parts. The larger part is related to a free stream of fluid, in which the effect of viscosity is negligible (potential flow theory). The smaller region is a thin layer next to the surface of the body, in which the effects of molecular transport (such as viscosity, thermal conductivity and mass diffusivity) are very important.

Laminar and turbulent boundary layer flow over a flat plate
Figure 2: Laminar and turbulent boundary layer flow over a flat plate.

Potential flow theory neglects the effect of viscosity, and therefore, significantly simplifies the Navier-Stokes equations, which provides the solution for the velocity distribution. A disadvantage of the potential theory is that since second order terms are neglected, the effect of viscosity, no slip, and impermeability boundary conditions at the surface cannot be accounted for. In general, potential flow theory predicts the free stream field accurately, despite its simplicity. Boundary layer thickness is defined as the distance within the fluid in which most of the velocity change occurs. This thickness is usually defined as the thickness in which the velocity reaches 99% of the free stream velocity, u = 0.99U.

Figure 2 illustrates how the momentum boundary layer thickness changes along the plate. Flow is laminar at relatively small values of x, where it is shown as the laminar boundary layer region. As x increases, the fluid motion begins to fluctuate. This is called the transition region. The boundary layer may be either laminar or turbulent in this region. Finally, above a certain value of x, the flow will always be turbulent. There is a very thin region next to the wall in the turbulent region where the flow is still laminar, called the laminar sublayer.

For flow over a flat plate, experimental data indicates

Rex ≤ 2×105 the flow is laminar
2×105 <Rex <3×106 the flow is in transition
Rex ≥ 3×106 the flow is turbulent


\operatorname{Re} _x  = \frac{{\rho U_\infty  x}}
{\mu }

Boundary layer approximation

If one assumes that the boundary layer thickness, δ, is very small compared to the characteristic dimension of the object, one can make the assumption that δ is significantly less than L (δ L), where L is the characteristic dimension of the object. Using scale analysis discussed in Chapter 1, for flow over a flat plate with constant free stream velocity, one can show, in a steady two-dimensional laminar boundary layer representation, that the following conditions are met within the boundary layer region, assuming there are no body forces:

\frac{{\partial ^2 u}}{{\partial y^2 }} \gg \frac{{\partial ^2 u}}{{\partial x^2 }}
u \gg v
\frac{{\partial p}}{{\partial x}} \approx \frac{{dp}}{{dx}}
\frac{{\partial p}}{{\partial y}} \approx 0
Boundary layer concept over a flat plate
Figure 3: Boundary layer concept over a flat plate.

A similar concept exists when there is heat and/or mass transfer between a fluid and the surface of an object. Again, the region in which the effect of temperature or concentration is dominant is, in the most practical case, a region very close to the surface, as shown in Fig. 3. In general, for the case of flow over a flat plate, one can expect three different boundary layer regions with thicknesses δ, δT, and δC, corresponding to momentum, thermal, and concentration boundary layers, respectively. δ, δT, and δC are not necessarily the same thickness and their values, as will be shown later, depend on the properties of the fluid such as kinematic viscosity ( ), thermal diffusivity (α = k/ρcp), specific heat, c, and mass diffusion coefficient, D. Using boundary layer approximation, scale analysis, and order of magnitude, one can show that similar approximations exist for thermal and mass concentration boundary layer analysis.

\frac{{\partial T}}{{\partial y}} \gg \frac{{\partial T}}{{\partial x}},\quad \frac{{\partial ^2 T}}{{\partial y^2 }} \gg \frac{{\partial ^2 T}}{{\partial x^2 }}
\frac{{\partial \omega }}{{\partial y}} \gg \frac{{\partial \omega }}{{\partial x}},\quad \frac{{\partial ^2 \omega }}{{\partial y^2 }} \gg \frac{{\partial ^2 \omega }}{{\partial x^2 }}

where ω is the mass fraction.

Governing Equations for Boundary Layer Approximation

As noted before, one can obtain all pertinent information related to momentum, heat, and mass transfer between a surface and a fluid by focusing on the thin region (boundary layer) next to the surface, and solving the governing equations including the boundary conditions. This provides significant simplification irrespective of whether an analytical or numerical approach to solve the physical problem is used. For most practical applications, the effect of molecular transport due to mass, momentum, energy, and species is dominant in this thin region. The purpose of this section is to develop the transport phenomena equations for boundary layer approximation, starting from the original differential conservation equations developed in Chapter 2. Consider flow over a flat plate as shown in Figure 4.4 with constant free stream velocity, temperature, and mass fraction of U∞, T∞, and ω∞, respectively. The surface wall is kept at constant temperature and concentration. Starting from the differential conservation equations for mass, momentum, energy, and species (as presented in Chapter 2), and assuming two-dimensional, steady, laminar flow, and constant properties in a Cartesian coordinate system, the following results are obtained:

Continuity equation

 \frac{{\partial u}}
{{\partial x}} + \frac{{\partial v}}
{{\partial y}} = 0

Momentum equation in the x-direction

u\frac{{\partial u}}
{{\partial x}} + v\frac{{\partial u}}
{{\partial y}} =  - \frac{1}
{\rho }\frac{{\partial p}}
{{\partial x}} + \nu \left( {\frac{{\partial ^2 u}}
{{\partial x^2 }} + \frac{{\partial ^2 u}}
{{\partial y^2 }}} \right)

Momentum equation in the y-direction

u\frac{{\partial v}}
{{\partial x}} + v\frac{{\partial v}}
{{\partial y}} =  - \frac{1}
{\rho }\frac{{\partial p}}
{{\partial y}} + \nu \left( {\frac{{\partial ^2 v}}
{{\partial x^2 }} + \frac{{\partial ^2 v}}
{{\partial y^2 }}} \right)

Energy equation

u\frac{{\partial T}}
{{\partial x}} + v\frac{{\partial T}}
{{\partial y}} = \alpha \left( {\frac{{\partial ^2 T}}
{{\partial x^2 }} + \frac{{\partial ^2 T}}
{{\partial y^2 }}} \right) + \frac{\nu }
{{cp}}\left( {\frac{{\partial u}}
{{\partial y}}} \right)^2