View source for Boundary-layer theory

An asymptotic approximation of the solution of boundary value problems for differential equations containing a small parameter in front of the highest derivative (singular problems) in subregions where there is a substantial effect from the terms containing the highest derivatives on the solution. The boundary-layer phenomenon arises in narrow zones near the parts of the boundary on which there is a difference in the numbers of boundary conditions for the initial problem and the degenerate one (with the small parameter taking the value zero), as well as near the surfaces of discontinuity of the solution of the degenerate problem.

The solution of the singular problem may be represented as the sum of two expansions. The outer expansion is determined by the method of the small parameter with the part <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173101.png" /> of the boundary conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173102.png" />. The inner expansion decreases rapidly outside the boundary layer and is usually sought for as polynomials in powers of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173103.png" />. To determine these, the differential equations are transformed by variables that depend on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173104.png" /> and that stretch the subregions of the boundary layer. The equations of the boundary layer are derived by equating to zero the coefficients at various powers of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173105.png" /> after substituting the polynomials into the transformed equations. To these one adds conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173106.png" />. The error in the outer expansion, after it has been found, indicates the necessary change of variables. In solving complicated applied problems, the necessary change of variables can be clarified on the basis of physical estimates for the terms in the initial equations and the simplifications corresponding to them. This change should eliminate the highest derivatives in front of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173107.png" />. To solve the problem it is necessary to determine where the boundary layers lie and how the conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173108.png" /> are to be separated into <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b0173109.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731010.png" />. A notable feature of this solution is that hyperbolic equations for the outer expansions and parabolic equations for the inner expansions may correspond to initially elliptic equations.

In the theory of systems of ordinary differential equations

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731011.png" /></td> <td valign="top" style="width:5%;text-align:right;">(1)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731012.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731013.png" /> are <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731014.png" />-dimensional and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731015.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731016.png" /> are <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731017.png" />-dimensional vector functions, an existence and uniqueness theorem for the solution of the Cauchy problem under the conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731018.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731019.png" /> and under certain properties of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731020.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731021.png" /> has been proved. One has also derived properties of the solution as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731022.png" /> (see [[#References|[1]]]). In the case of a boundary value problem for (1) under the conditions

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731023.png" /></td> </tr></table>

where the sum of the numbers of components of the vectors <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731024.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731025.png" /> is <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731026.png" />, there exist, in general, boundary layers in neighbourhoods of the ends of the segment <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731027.png" />. An algorithm has been constructed for finding asymptotics of the solution for this problem. Under certain properties of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731028.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731029.png" /> it has been shown that a solution exists and is unique, and it has been estimated (see [[#References|[3]]]). If a solution to the limit equation <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731030.png" /> is not unique with respect to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731031.png" />, one can construct an internal boundary layer (in a neighbourhood of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731032.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731033.png" />) that separates regions with different solutions of the limit equation. An algorithm has been constructed for a particular type of integro-differential equation giving an asymptotic expansion in <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731034.png" /> for the problem with initial conditions, and some features of the behaviour of the solution have been examined.

In the case of linear ordinary differential equations <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731035.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731036.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731037.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731038.png" /> are differential operators, with boundary conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731039.png" />, one can distinguish the class of problems the solutions of which contain boundary layers, and the concept of regular degeneracy (the solution to the limit equation enables one to satisfy the conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731040.png" />, while the asymptotic solution for the boundary layer enables one to satisfy <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731041.png" />) can been introduced (see ). An iterative process has been constructed for the asymptotic representation of the solution, and estimates have been given for the residual terms in the expansions.

It has been shown (see [[#References|[5]]]) in the theory of boundary layers of a general non-linear second-order ordinary differential equation, subject to certain assumptions, that the solution of the first boundary value problem is made up of an external solution, a boundary layer and a residual term that, together with its first-order derivative, is of order <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731042.png" /> on the segment.

Studies have been made on the behaviour of solutions of boundary value problems of basic types for a linear partial differential equation of the form

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731043.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731044.png" /> is the Laplace operator, in a region <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731045.png" /> with boundary <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731046.png" />. Conditions on the functions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731047.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731048.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731049.png" />, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731050.png" />, the boundary <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731051.png" />, and the functions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731052.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731053.png" /> of the points <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731054.png" /> on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731055.png" /> appearing in the boundary condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731056.png" /> have been given such that <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731057.png" /> in <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731058.png" /> tends uniformly to the solution of the limit equation with this boundary condition on a certain part of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731059.png" /> (absence of a boundary layer) [[#References|[6]]].

For an elliptic second-order equation in a region <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731060.png" /> with boundary <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731061.png" />, using the example of two independent variables

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731062.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731063.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731064.png" /></td> </tr></table>

iterative processes have been constructed solving the problem with the condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731065.png" /> on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731066.png" />, theorems have been proved on the structure of the expansion of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731067.png" /> with respect to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731068.png" /> and estimates have been made of the residual term in this expansion . Similar results have been obtained for equations of higher orders.

A method has been devised [[#References|[7]]] for combining asymptotic expansions for the equation

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731069.png" /></td> </tr></table>

in a rectangle with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731070.png" /> given on the boundary.

Research on boundary-layer theory for non-linear partial differential equations is related mainly to aerohydrodynamics and is based on the Navier–Stokes equations or generalizations of them. Practical requirements have led to the development of the mathematical theory and to methods of handling various problems. Below only laminar flows are considered (see [[#References|[8]]]–[[#References|[10]]]).

The Navier–Stokes equations in the case of the hydrodynamics of planar (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731071.png" />) and axi-symmetric (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731072.png" />) flows of an incompressible fluid with constant viscosity coefficient <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731073.png" /> are:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731074.png" /></td> <td valign="top" style="width:5%;text-align:right;">(2)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731075.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731076.png" /> is the Reynolds number, which is represented in terms of the characteristic values of the velocity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731077.png" /> and the linear dimension <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731078.png" />. The solutions are defined by the boundary conditions on the closed boundary <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731079.png" /> of a region <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731080.png" />, where on the solid contour <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731081.png" /> one has the conditions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731082.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731083.png" /> for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731084.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731085.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731086.png" /> are the tangential and normal components to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731087.png" /> of the vector <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731088.png" />. The initial values of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731089.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731090.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731091.png" /> are given on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731092.png" />.

For <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731093.png" /> small, in a first approximation the asymptotic solution is composed of the solution of (2) for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731094.png" /> with some of the conditions on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731095.png" /> (only the condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731096.png" /> is imposed on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731097.png" />) and the solution of the boundary-layer equations. The equations for the dynamic boundary layer are derived on the assumption that the conditional thickness <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731098.png" /> of the boundary layer and the value of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b01731099.png" /> have orders <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310100.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310101.png" />, and that the terms on the left-hand sides of the latter equations in (2) are of the order of the terms containing <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310102.png" />. Introduction of the variables <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310103.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310104.png" /> leads, as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310105.png" />, to the Prandtl equations (cf. [[Prandtl equation|Prandtl equation]]):

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310106.png" /></td> <td valign="top" style="width:5%;text-align:right;">(3)</td></tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310107.png" /></td> </tr></table>

with the conditions

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310108.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310109.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310110.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310111.png" /> is the distance from the symmetry axis for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310112.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310113.png" /> is a known function. These equations and conditions apply for any curvilinear contour with radius of curvature much larger than <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310114.png" />. In the latter case, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310115.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310116.png" /> are the coordinates along the contour and along the normal to it.

If <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310117.png" /> is constant, the problem reduces to a boundary value problem for an ordinary differential equation. There are also other classes of analogous solutions.

Conditions under which solutions of boundary-layer problems exist are known; one has investigated the problem of uniqueness and stability of the solutions as well as how they result from the solutions for stationary cases [[#References|[11]]]. Solutions have been constructed by the method of straight lines, and it has been shown that they converge.

Boundary-layer equations for a compressible liquid can be derived from the equations for the flow of a viscous and heat-conducting gas; they are much more complicated than (3). Their number is also larger. There is an integral transformation that simplifies these equations in the general case and reduces them to (3) when the Prandtl number <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310118.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310119.png" /> is the heat capacity of the gas at constant pressure and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310120.png" /> is the coefficient of heat conductivity [[#References|[12]]]. Several modifications of the transformation exist. In the general case, the boundary-layer equations describe so-called natural convective flows. If <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310121.png" /> is independent of the temperature and if the Archimedean force is negligible, then the energy equation splits off from the system of boundary layer equations and one speaks of forced convective flow. The energy equation determines the thermal boundary layer, whose thickness differs from <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310122.png" />.

Boundary layers also arise in zones separating flows with different characteristic velocities. Shock waves are also boundary layers.

A distinct class of two-dimensional boundary-layer problems is associated with flows in rotating axi-symmetric plates and bodies.

Not only have methods been developed for solving non-stationary problems, but also problems in which <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310123.png" /> is periodic, when there is stepwise motion from a state of rest, problems with accelerated motion, for the boundary layer behind a shock wave, and problems with a variable temperature over the surface of the solid around which the flow takes place have been solved.

In the boundary-layer theory for three-dimensional flows, methods for obtaining a solution have been developed and cases in which the equations simplify have been studied. The boundary-layer equations for a sliding cylindrical wing of infinite span are analogous to the equations for a two-dimensional boundary layer. Approximate solutions have been obtained for the problem of a boundary layer on a rotating cylindrical propeller blade and on a rotating cylinder in a skew flow, as well as for the boundary-layer problem near the line of intersection between two planes.

These researches on boundary layers in aerohydrodynamics relate to a first approximation in boundary-layer theory. Higher approximations enable one to examine the interactions of boundary layers with the external flow, and to make calculations for moderate values of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310124.png" /> [[#References|[13]]].

The stability of boundary layers enables one to determine the limits of applicability of the theory. There are studies [[#References|[14]]] based on methods of small perturbations with periodic and local initial perturbations. In the case of two-dimensional flows, the analysis of three-dimensional perturbations reduces in linear approximation, on the basis of Squire's theorem, to a two-dimensional analysis with an altered value of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310125.png" />. Non-linear analysis applied to stability loss shows that longitudinal vortices occur.

Physical generalizations of problems in the theory of boundary layers (see [[#References|[15]]]–[[#References|[17]]]) are related to research on multi-phase flows, to the use of real equations of state and transfer coefficients (a complication of the equations), to the study of non-equilibrium flows with diffusion (an extension of the system of equations, which then becomes of parabolic-hyperbolic type), to taking into account ablation of the surface around which the flow takes place (a complication of the boundary conditions and it becomes necessary to consider the thermal conduction in the solid), and also to considering radiation transport (integro-differential equations).

A further development of boundary-layer theory in aerohydrodynamics has been obtained in the study of flows that do not satisfy Prandtl's assumptions, where solutions are obtained as multi-layer asymptotic expansions. The causes of the complexity in the structure of the solution are that there are additional small parameters in the boundary conditions (for example, because the radius of curvature of the contour around which the flow takes place is small), there may be singular points, lines or surfaces in the first approximation of the theory, and bifurcation of the solution is possible.

This class also contains flows around points of separation or connections of the boundary layers to the contour around which the flow takes place and around the points of incidence of shock waves on the boundary layers. The characteristic solutions have ([[#References|[18]]]–[[#References|[20]]]) a three-layer structure. The external solution defines the potential flow perturbed by the boundary layer and is described by equations in perturbations. The middle layer of order <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310126.png" /> can be described by the equations for non-viscous vortex flows with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310127.png" />, and it receives gas from the main part of the preceding boundary layer.

The equations for the third (thinnest) layer nearest to the wall can be derived on the assumption that its length and thickness are, correspondingly, of order <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310128.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310129.png" />. The following variables are introduced in the stationary case:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310130.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310131.png" /></td> </tr></table>

This leads, as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310132.png" />, again to (3) in which <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310133.png" /> is replaced by <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310134.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310135.png" /> by <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310136.png" />. This structure for the solution of a problem is characteristic for a wide class of flows with small perturbations.

Many studies have been made on the dynamics of flows for small <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310137.png" />, where the pressure in the outer supersonic flow varies considerably over short distances (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310138.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310139.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310140.png" /> is the velocity of the gas and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310141.png" /> is the velocity of sound). This includes problems on calculating flows around contours of large local curvature and flow attachment to the surface of the body. In these cases <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310142.png" /> in the middle layer of the three-layer scheme for the solution.

The class of problems in which the perturbed solution occupies a finite region has also been studied. These solutions occur when there is moderate or strong interaction of the boundary layer with an external hypersonic flow (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310143.png" />) or in supersonic flow around a body of finite length in which there are intense gas injections through the surface (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310144.png" />). In these cases, the pressure on the outer boundary of the boundary layer is determined from the solution to the complete problem. The perturbation from the trailing edge of the body propagates upstream, which is caused by the non-uniqueness of the solution in a neighbourhood of the leading edge.

====References====
<table><TR><TD valign="top">[1]</TD> <TD valign="top">  A.N. Tikhonov,   "Systems of differential equations containing small parameters in front of the derivatives"  ''Mat. Sb.'' , '''31''' :  3  (1952)  pp. 575–586  (In Russian)</TD></TR><TR><TD valign="top">[2]</TD> <TD valign="top">  W. Wasov,   "Asymptotic expansions for ordinary differential equations" , Interscience  (1965)</TD></TR><TR><TD valign="top">[3]</TD> <TD valign="top">  A.B. Vasil'eva,   V.F. Butuzov,   "Asymptotic expansions of solutions of singularly perturbed equations" , Moscow  (1973)  (In Russian)</TD></TR><TR><TD valign="top">[4a]</TD> <TD valign="top">  M.I. Vishik,   L.A. Lyusternik,   "Regular degeneracy and boundary layers for linear differential equations with a small parameter"  ''Uspekhi Mat. Nauk'' , '''12''' :  5  (1957)  pp. 3–122  (In Russian)</TD></TR><TR><TD valign="top">[4b]</TD> <TD valign="top">  M.I. Vishik,   L.A. Lyusternik,   "Solutions to some problems on perturbations in the case of matrices and non-self-adjoint differential equations"  ''Uspekhi Mat. Nauk'' , '''15''' :  3  (1960)  pp. 3–80  (In Russian)</TD></TR><TR><TD valign="top">[5]</TD> <TD valign="top">  J.D. Cole,   "Perturbation methods in applied mathematics" , Blaisdell  (1968)</TD></TR><TR><TD valign="top">[6]</TD> <TD valign="top">  O.A. Oleinik,   "on equations of elliptic type with a small parameter in front of the highest derivatives"  ''Mat. Sb.'' , '''31''' :  1  (1952)  pp. 104–117  (In Russian)</TD></TR><TR><TD valign="top">[7]</TD> <TD valign="top">  A.M. Il'in,   E.F. Lelikova,   "The method of matched asymptotic expansions for the equation <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b017/b017310/b017310145.png" /> in a rectangle"  ''Mat. Sb.'' , '''96''' :  4  (1975)  pp. 568–583  (In Russian)</TD></TR><TR><TD valign="top">[8]</TD> <TD valign="top">  N.E. Kochin,   I.A. Kibel',   N.V. Roze,   "Theoretical hydrodynamics" , '''2''' , Interscience  (1964)  (Translated from Russian)</TD></TR><TR><TD valign="top">[9]</TD> <TD valign="top">  L.G. Loitsyanskii,   "Laminar boundary layers" , Moscow  (1962)  (In Russian)</TD></TR><TR><TD valign="top">[10]</TD> <TD valign="top">  H. von Schlichting,   "Grenzschicht-Theorie" , Karlsruhe  (1951)</TD></TR><TR><TD valign="top">[11]</TD> <TD valign="top">  O.A. Oleinik,   "Mathematical problems of boundary-layer theory"  ''Uspekhi Mat. Nauk'' , '''23''' :  3  (1968)  pp. 3–65  (In Russian)</TD></TR><TR><TD valign="top">[12]</TD> <TD valign="top">  A.A. Dorodnitsyn,   "Boundary layers in an incompressible gas"  ''Prikl. Mat. i Mekh.'' , '''6''' :  6  (1942)  pp. 449–486  (In Russian)</TD></TR><TR><TD valign="top">[13]</TD> <TD valign="top">  M. van Dyke,   "Perturbation methods in fluid mechanics" , Parabolic Press  (1975)</TD></TR><TR><TD valign="top">[14]</TD> <TD valign="top">  R. Betchov,   V.O. Criminale,   "Stability of parallel flows" , Acad. Press  (1967)</TD></TR><TR><TD valign="top">[15]</TD> <TD valign="top">  S. Soo,   "Fluid dynamics of multiphase systems" , Blaisdell  (1971)</TD></TR><TR><TD valign="top">[16]</TD> <TD valign="top">  W.H. Dorrance,   "Viscous hypersonic flow" , McGraw-Hill  (1962)</TD></TR><TR><TD valign="top">[17]</TD> <TD valign="top">  W.M. Kays,   "Convective heat and mass transfer" , McGraw-Hill  (1971)</TD></TR><TR><TD valign="top">[18]</TD> <TD valign="top">  V.Ya. Neiland,   "On the theory of separation of laminar boundary layers in supersonic flow"  ''Izv. Akad. Nauk SSSR Mekh. Zhidk. i Gaza'' :  4  (1969)  pp. 53–57  (In Russian)</TD></TR><TR><TD valign="top">[19]</TD> <TD valign="top">  V.Ya. Neiland,   "Asymptotic problems in the theory of viscous flows"  ''Trudy TsAGI'' :  1529  (1974)  (In Russian)</TD></TR><TR><TD valign="top">[20]</TD> <TD valign="top">  K. Stewartson,   "Multistructured boundary layers on flat plates and related bodies"  ''Adv. App. Mech.'' , '''14'''  (1974)  pp. 145–239</TD></TR></table>



====Comments====
The article above presents a one-sided view of the subject and is in many aspects not up-to-date. Furthermore, very substantial contributions in the Western literature are neglected. Mathematically boundary layers arise in problems of singular perturbations and are treated by methods of asymptotic analysis. In linear and non-linear, ordinary and partial, differential equations much more has been accomplished than described above in the first part. Overviews that are partially complementary can be found in [[#References|[a1]]]–[[#References|[a3]]].

From the point of view of applications boundary-layer theory is described above as a branch of fluid- and gasdynamics, where it originated. However, boundary layers and transition layers also occur in many other fields, such as for example combustion, geophysics, free boundary problems, or epidemology. Samples and further literature can be found in [[#References|[a4]]], [[#References|[a5]]].

Systems of (ordinary) differential equations of the type (1) above are often referred to as singular perturbations (cf. also the editorial comments to [[Perturbation theory|Perturbation theory]]).

For an idea of the role of singular perturbations, boundary layers and multiple time scales in control and system theory, cf. [[#References|[a6]]].

====References====
<table><TR><TD valign="top">[a1]</TD> <TD valign="top">  J.L. Lions,   "Perturbation singulières dans les problèmes aux limits et en contrôle optimal" , ''Lect. notes in math.'' , '''323''' , Springer  (1973)</TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top">  W. Eckhaus,   "Asymptotic analysis of singular perturbations" , North-Holland  (1979)</TD></TR><TR><TD valign="top">[a3]</TD> <TD valign="top">  K.W. Chang,   F.A. Howes,   "Nonlinear singular perturbation phenomena: theory and application" , Springer  (1984)</TD></TR><TR><TD valign="top">[a4]</TD> <TD valign="top">  C.M. Brauner (ed.)  B. Gay (ed.)  J. Mathieu (ed.) , ''Singular perturbations and boundary layer theory'' , ''Lect. notes in math.'' , '''594''' , Springer  (1977)</TD></TR><TR><TD valign="top">[a5]</TD> <TD valign="top">  F. Verhulst (ed.) , ''Asymptotic analysis II'' , ''Lect. notes in math.'' , '''985''' , Springer  (1983)</TD></TR><TR><TD valign="top">[a6]</TD> <TD valign="top">  P. Kokotovic (ed.) , ''Singular perturbations and time scales in modelling and control of dynamical systems'' , Univ. Illinois at Urbana-Champaign  (1980)</TD></TR></table>