Triple deck theory is a theory that describes a three-layered boundary-layer structure when sufficiently large disturbances are present in the boundary layer. This theory is able to successfully explain the phenomenon of boundary layer separation, but it has found applications in many other flow setups as well, including the scaling of the lower-branch instability (T-S) of the Blasius flow, etc. James Lighthill, Lev Landau and others were the first to realize that to explain boundary layer separation, different scales other than the classical boundary-layer scales need to be introduced. These scales were first introduced independently by James Lighthill and E. A. Müller in 1953. The triple-layer structure itself was independently discovered by Keith Stewartson (1969) and V. Y. Neiland (1969) and by A. F. Messiter (1970). Stewartson and Messiter considered the separated flow near the trailing edge of a flat plate, whereas Neiland studied the case of a shock impinging on a boundary layer.
Suppose
and
are the streamwise and transverse coordinate with respect to the wall and
be the Reynolds number, the boundary layer thickness is then
. The boundary layer coordinate is
. Then the thickness of each deck is

The lower deck is characterized by viscous, rotational disturbances, whereas the middle deck (same thickness as the boundary-layer thickness) is characterized by inviscid, rotational disturbances. The upper deck, which extends into the potential flow region, is characterized by inviscid, irrotational disturbances.
The interaction zone identified by Lighthill in the streamwise direction is

The most important aspect of the triple-deck formulation is that pressure is not prescribed, and so it has to be solved as part of the boundary-layer problem. This coupling between velocity and pressure reintroduces ellipticity to the problem, which is in contrast to the parabolic nature of the classical boundary layer of Prandtl.
Flow near the trailing edge of a flat plate
Let the length scales be normalized with the plate length
and the velocity scale by the free-stream velocity
; then the only parameter in the problem is the Reynolds number
. Let the origin
of the coordinate system be located at the trailing edge of the plate. Further let
be the non-dimensional velocity components,
be the non-dimensional pressure field and
be the non-dimensional stream function such that
and
. For shortness of notation, let us introduce the small parameter
. The coordinate for horizontal interaction and for the three decks can then be defined by

As
(or
), the solution should approach the asymptotic behaviour of the Blasius solution, which is given by
![{\displaystyle {\frac {\psi }{\varepsilon ^{4}}}={\sqrt {2}}f_{B}\left({\frac {\eta }{\sqrt {2}}}\right)+{\frac {x}{\sqrt {2}}}\left[f_{B}\left({\frac {\eta }{\sqrt {2}}}\right)-{\frac {\eta }{\sqrt {2}}}f_{B}'\left({\frac {\eta }{\sqrt {2}}}\right)\right]+\cdots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/bccc59e463f7721f8b207ff6c3e2fcbd3cb59bbc)
where
is the Blasisus function which satisfies
subjected to
. As
(or
), the solution should approach the asymptotic behaviour of the Goldstein's near wake, which is given by

where
and
. The Goldstein's inner wake solution is not needed here.
Upper deck
In the upper deck, the solution is found to given by
![{\displaystyle {\frac {\psi }{\varepsilon ^{4}}}={\frac {\zeta }{\varepsilon }}-{\sqrt {2}}\beta -{\frac {\varepsilon }{\pi }}\int _{-\infty }^{+\infty }A'(\chi )\left[\tan ^{-1}\left({\frac {\chi -{\hat {\chi }}}{\zeta }}\right)+{\frac {\pi }{2}}\right]d{\hat {\chi }}+\cdots ,\quad p=\varepsilon ^{2}P(\chi )+\cdots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f1a92c28a7244439ec2dbf61c2cf1a9d54f65fe)
where
. Furthermore, the upper deck problem also provides the relation between the displacement and the pressure function as

in which
stands for Cauchy principal value. One may notice that the pressure function and the derivative of the displacement function (aka transverse velocity) forms a Hilbert transform pair (see Kramers–Kronig relations).
Lower deck
In the lower deck, the solution is given by

where
will satisfy a boundary-layer type equations driven by the pressure gradient
and the slip-velocity of order
generated by the middle deck. It is convenient to introduce
and
, where
and
must satisfy

These equations are subjected to the conditions
- :{\hat {u}}\to {\frac {\xi +A(\chi )}{\sqrt {2\alpha ^{3}}}},}

- :\,\,\,A\to 0,\quad {\text{and}}\quad \chi \to +\infty :\,\,\,A\to {\frac {\mu }{\lambda }}\chi ^{1/3}}

where
. The displacement function
and therefore
must be obtained as part of the solution. The above set of equations may resemble normal boundary-layer equations, however it has an elliptic character since the pressure gradient term now is non-local, i.e., the pressure gradient at a location
depends on other locations as well. Because of this, these equations are sometimes referred to as the interactive boundary-layer equations. The numerical solution of these equations were obtained by Jobe and Burggraf in 1974.