View source for Shock waves, mathematical theory of

A mathematical description of properties, motions and interactions with the surrounding medium of a surface of discontinuity of the parameters of a medium (a so-called shock wave). In a wider and more abstract sense the mathematical theory of shock waves describes the properties of a surface of discontinuity of a solution of a first-order quasi-linear hyperbolic system of partial differential equations (cf. [[Quasi-linear hyperbolic equations and systems|Quasi-linear hyperbolic equations and systems]]). The mathematical theory of shock waves arose in connection with problems of the motion of gases and compressible fluids in the second half of the 19th century. Its foundations were established in the works of S. Earnshaw, B. Riemann, W. Rankine, H. Hugoniot (cf., for example, [[#References|[1]]]–[[#References|[4]]]).

As an idealization of real gases and fluids one considers a medium without dissipation, lacking viscosity and thermal conductivity. In the process of motion in such idealized media there may arise discontinuities in the distribution of all parameters of the flow (the density, the pressure, the temperature, the velocity, etc.). The set of discontinuity points of the parameters of the flow may be very complicated. Only the simplest basic cases have been systematically considered, when this set has piecewise-smooth surfaces of discontinuity consisting of discontinuity points of the parameter of the first kind. In the general case two-dimensional surfaces of discontinuity move in the three-dimensional space <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849401.png" /> with the flow of time. Shock waves are one of the possible types of surfaces of discontinuity.

The appearance of discontinuities very much complicates the mathematical formulation of the problem of the flow of ideal gases and fluids, because discontinuous functions cannot be solutions of the differential equations of gas dynamics (hydrodynamics). Therefore the flows with surfaces of discontinuity are described by generalized solutions of systems of quasi-linear equations of gas dynamics (cf. [[Gas dynamics, equations of|Gas dynamics, equations of]]) and the mathematical theory of shock waves forms part of the theory of generalized solutions of the systems of integral conservation laws of gas dynamics.

==Surfaces of discontinuity.==
On the surfaces of discontinuity the conditions that follow from the integral conservation laws of mass, momentum and energy have to be satisfied. The only exception are the discontinuities at the moment of starting the motion (the so-called initial discontinuities), which may be arbitrary. Let <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849402.png" /> be a smooth surface of discontinuity of the parameters of a flow of a gas (fluid), and let <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849403.png" /> be the normal velocity of the motion of the surface of discontinuity. Only homogeneous media which are characterized by a density <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849404.png" />, a pressure <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849405.png" />, an internal energy of a unit mass of gas <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849406.png" />, and a velocity vector of the motion of the medium <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849407.png" /> will be considered here.

In a point of the surface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849408.png" />, let <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s0849409.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494010.png" /> is the normal and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494011.png" /> the tangential velocity component (relative to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494012.png" />) of the velocity vector <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494013.png" />. The continuity conditions for the flow of mass, momentum and energy on the surface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494014.png" /> are expressed in the form of 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/s/s084/s084940/s08494015.png" /></td> <td valign="top" style="width:5%;text-align:right;">(1)</td></tr></table>

where the square brackets denote the jump in the value of the quantity enclosed in passing from one side of the surface to the other, i.e. <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494016.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494017.png" /> are the limit values of the quantity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494018.png" /> at the point <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494019.png" /> when approaching from the two sides. Conditions (1) are intrinsic boundary conditions which are added to the equations of gas dynamics on surfaces of discontinuity of the parameters of the medium.

There are two types of discontinuity: a tangential discontinuity for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494020.png" />; <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494021.png" />, and shock waves for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494022.png" />. The vector <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494023.png" /> is called the mass flow across unit area of the moving surface of discontinuity (shock wave) <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494024.png" />.

==Tangential and contact discontinuities.==
For a tangential discontinuity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494025.png" /> the pressure is continuous, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494026.png" /> on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494027.png" />, while the quantities <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494028.png" /> may have arbitrary jumps on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494029.png" />. If at least one of the quantities <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494030.png" /> is different from zero, the discontinuity is also called a contact discontinuity. In the case of a contact discontinuity the surface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494031.png" /> is the dividing boundary of the media with different properties, in particular with different equations of state; it is formed by the lines of flow emitted from the surface of the initial contact discontinuity.

A tangential discontinuity (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494032.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494033.png" />) is unstable, because under the action of even a small viscosity intrinsic to the gases and fluids, the initial tangential discontinuity represented by the component <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494034.png" /> of the velocity vector diffuses into a much wider zone of continuous transition. With the exception of singular cases, when this discontinuity is realized in a short interval of time (tangential discontinuity in a nozzle, initial regions of a zone of mixing of different flows, etc.), one does not consider it as an admissible discontinuity. This is a manifestation of the limitations imposed by considering  "ideal"  media or media with a small viscosity, since in an ideal medium the tangential discontinuity is preserved. Therefore the question of the admissibility or non-admissibility of the tangential discontinuity of the flow of an ideal medium has to be analyzed for each separate concrete problem.

==Pure contact discontinuities.==
Now <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494035.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494036.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494037.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494038.png" />. These may collapse by processes of molecular diffusion; as a rule, these processes are so slow that the contact discontinuity may be regarded as stable over finite intervals of time.

However, a contact discontinuity may also become unstable and rapidly collapse if in the course of a large interval of time the acceleration of the medium on the contact boundary is directed towards the side of the medium with lower density. In this case the contact discontinuity collapses and an  "intermixing"  of media that were originally on different sides of the contact discontinuity starts. A simple example of a similar, essentially many-dimensional, instability is the collapse in the field of gravity of a horizontal contact discontinuity between a light fluid, found in the lower part of a vessel, in the presence of a much heavier fluid above it (so-called Rayleigh–Taylor instability).

==Shock waves.==
In this case <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494039.png" /> and (1) implies <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494040.png" />, i.e. the tangential discontinuity of the velocity vector is absent. Condition (1) on the shock wave <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494041.png" /> leads to the three 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/s/s084/s084940/s08494042.png" /></td> <td valign="top" style="width:5%;text-align:right;">(2)</td></tr></table>

which are called the Rankine–Hugoniot conditions. If <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494043.png" />, then the shock wave is called straight, in the opposite case oblique. The elimination of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494044.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494045.png" /> from (2) gives the relation

<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/s/s084/s084940/s08494046.png" /></td> <td valign="top" style="width:5%;text-align:right;">(3)</td></tr></table>

which relates only thermodynamic parameters of the medium on the two sides of the shock wave. This relation is called the adiabatic shock, or the Hugoniot adiabatic.

Let the surface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494047.png" /> be moving in the direction <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494048.png" /> from the left to the right. Then if <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494049.png" />, the shock wave relative to the medium moves to the right, and the matter crosses the shock wave <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494050.png" /> moving relative to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494051.png" /> from the right to the left. Conversely, if <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494052.png" />, then the shock wave moves to the left relative to the medium. The medium to the right of the shock wave for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494053.png" /> and to the left of the shock wave for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494054.png" /> is called the medium in front of the shock wave, the other medium is the medium beyond the shock wave. The parameters of the medium in front of the shock wave are <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494055.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494056.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494057.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494058.png" />; the parameters of the medium beyond the shock wave are <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494059.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494060.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494061.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494062.png" />. In the process of motion the part in front of the shock wave passes across the surface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494063.png" /> of the shock wave and joins the medium beyond the shock wave <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494064.png" />. Moreover, it passes through a zone of high gradients (a shock wave zone) in which one cannot neglect dissipative processes — friction and heat conduction — in spite of small viscosity and thermal conductivity coefficients of the matter. In view of the irreversibility of the dissipative processes, the condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494065.png" /> must hold, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494066.png" /> is the entropy of the gas. This inequality does not allow an exchange of states <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494067.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494068.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494069.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494070.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494071.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494072.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494073.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494074.png" /> (as is allowed by (2), (3)), but defines the position of these states relative to the surface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494075.png" />.

The condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494076.png" /> is called the stability condition of the shock wave. The question of the realizability (admissibility) of a shock wave, on which (2), (3) and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494077.png" /> hold, can be solved comparatively simply only for a so-called normal gas. A gas (fluid) is called a normal gas if its equations of state have 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/s/s084/s084940/s08494078.png" /></td> <td valign="top" style="width:5%;text-align:right;">(4)</td></tr></table>

For a normal gas the condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494079.png" /> guarantees the stability of the shock wave; from the Hugoniot conditions (2), the state <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494080.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494081.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494082.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494083.png" /> beyond the shock wave front is uniquely determined by the given state <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494084.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494085.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494086.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494087.png" />, and the mass flow <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494088.png" /> if

<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/s/s084/s084940/s08494089.png" /></td> </tr></table>

From the stability condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494090.png" /> for the shock wave there follow the properties: a) the shock wave moves with supersonic speed through the medium before the front and with subsonic speed through the medium beyond the front (Tsemplen's theorem), i.e.

<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/s/s084/s084940/s08494091.png" /></td> </tr></table>

b) the shock wave leads to compressed matter and increases the pressure in it, i.e. <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494092.png" /> implies <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494093.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494094.png" />.

==Adiabatic shock.==
For a normal gas the adiabatic shock (3) is represented in the <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494095.png" />-plane by a curve that is downward convex, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494096.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494097.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494098.png" />, passing through the point <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s08494099.png" />, called the centre of the adiabatic shock. The graph of an adiabatic shock for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940100.png" /> lies above the adiabatic Poisson graph <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940101.png" />, and lies below it for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940102.png" />; at the point <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940103.png" /> these two adiabatics have contact of at least the second order. The state <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940104.png" /> beyond the front of the shock wave corresponds to the left branch of the adiabatic shock <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940105.png" />.

For a medium with equations of state

<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/s/s084/s084940/s084940106.png" /></td> </tr></table>

(the so-called ideal gas) conditions (4) hold and equation (3) takes 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/s/s084/s084940/s084940107.png" /></td> <td valign="top" style="width:5%;text-align:right;">(5)</td></tr></table>

Equation (5) defines a hyperbola with asymptotes <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940108.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940109.png" />. The limiting compression of the matter is attained behind the infinitely strong shock wave:

<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/s/s084/s084940/s084940110.png" /></td> </tr></table>

The solution of (2) for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940111.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940112.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940113.png" /> in terms of the given <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940114.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940115.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940116.png" />, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940117.png" /> is provided by the formulas

<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/s/s084/s084940/s084940118.png" /></td> <td valign="top" style="width:5%;text-align:right;">(6)</td></tr></table>

For gases with anomalous thermodynamic properties, when conditions (4) are violated, the requirement <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940119.png" /> does not guarantee stability (admissibility) of the shock transition. The question of the conditions for the admissibility of a shock transition is very complicated and, for gases with arbitrary equations of state satisfying only the necessary thermodynamic conditions, has not been resolved. The case most studied is that when of the conditions (4) only the condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940120.png" /> is violated and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940121.png" /> may have alternating sign. In this case the adiabatic shock (3) contains states <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940122.png" /> for which the shock transition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940123.png" /> is unstable even when the condition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940124.png" /> holds.

==The zone of shock transition and its width.==
This is the zone of continuous change of the parameters of a medium from its values before the front to its values beyond the shock wave front for real media, having finite viscosity and thermal conductivity. It is assumed here that a continuous transition exists only for stable shock waves. Therefore the study of continuous shock transitions gives admissibility criteria for shock waves. The character of change of the parameters of the medium in the transition zone is also studied.

For a narrow transition zone any shock wave may be considered as locally planar. Therefore the shock transition may be described as the solution of the one-dimensional equations of gas dynamics with plane symmetry for a viscous and thermally-conducting gas:

<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/s/s084/s084940/s084940125.png" /></td> <td valign="top" style="width:5%;text-align:right;">(7)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940126.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940127.png" /> are the viscosity and thermal conductivity coefficients (here assumed to be constant) and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940128.png" /> is the bulk viscosity. Equations (7) are supplemented by the equation of state <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940129.png" /> from <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940130.png" />. In a coordinate system moving with the velocity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940131.png" /> of the shock wave, the shock transition may be described as a stationary solution of equations (7). For these one obtains two 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/s/s084/s084940/s084940132.png" /></td> <td valign="top" style="width:5%;text-align:right;">(8)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940133.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940134.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940135.png" /> are constant.

Equations (8) must be solved by a solution <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940136.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940137.png" /> such that <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940138.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940139.png" /> as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940140.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940141.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940142.png" /> as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940143.png" />. Necessary conditions for the existence of such a solution,

<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/s/s084/s084940/s084940144.png" /></td> </tr></table>

reduce to Hugoniot's conditions (2).

A qualitative study of the portrait of the integral curves of the system (8) leads to the conclusion that there exists a unique integral curve for this problem provided conditions (4) are satisfied for the equations of state of the gas; the parameters <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940145.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940146.png" /> satisfy Hugoniot's conditions; the stability condition for shock waves <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940147.png" /> holds. The  "width"  of the zone of shock transition is infinite and the approximation to the limiting value takes place exponentially. In the case of a non-viscous thermally-conducting gas (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940148.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940149.png" />), which may be considered as the limit <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940150.png" /> of the general case, there exists a continuous solution of the problem of shock transition if the curve <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940151.png" /> in the plane of the variables <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940152.png" /> is monotone. In the opposite case (which happens for sufficiently strong shock waves) the limit for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940153.png" /> of the integral curves may have density discontinuity at constant temperature (the so-called isothermal jump). This means that the solutions of the equations of gas dynamics for thermally-conducting media without viscosity may be discontinuous. If the medium has viscosity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940154.png" /> and vanishing thermal conductivity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940155.png" />, then there always exists a continuous transition <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940156.png" />. From this one may conclude that the solutions of the equations of gas dynamics for viscous media cannot have a discontinuity of shock wave type, i.e. for them the Cauchy problem has a global solution (for any <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940157.png" />). However, a rigorous proof of this fact is still lacking.

Similarly one can study shock transitions for an ideal gas: <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940158.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940159.png" /> (where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940160.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940161.png" /> are constants). The equations (8) may in this case be integrated for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940162.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940163.png" />:

<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/s/s084/s084940/s084940164.png" /></td> </tr></table>

The effective width of the zone of shock transition is determined as follows:

<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/s/s084/s084940/s084940165.png" /></td> </tr></table>

A calculation shows that the width <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940166.png" /> of the shock wave has the order of magnitude of the free path of a molecule.

This justifies the point of view according to which the flow of a gas is divided into domains of reversible processes, where the flow can be described by the equations of gas dynamics without taking account of dissipative terms, and domains of irreversible processes, represented by narrow zones which may effectively be described by a moving surface of discontinuity (shock wave). More precisely, the behaviour of the parameters of the medium in the zone of shock transition is described by the [[Boltzmann equation|Boltzmann equation]] for gas-dynamic processes in non-equilibrium.

For gases with anomalous thermodynamic properties the study of the zone of shock transition is more complicated. The case when <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940167.png" /> alternates in sign has been analyzed. Let <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940168.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940169.png" /> in (8). Then <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940170.png" /> and the determination of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940171.png" /> reduces to the solution of a single 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/s/s084/s084940/s084940172.png" /></td> <td valign="top" style="width:5%;text-align:right;">(9)</td></tr></table>

The solution of this equation <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940173.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940174.png" /> as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940175.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940176.png" /> as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940177.png" />, exists if <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940178.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940179.png" /> lie on the adiabatic shock (3) and if

<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/s/s084/s084940/s084940180.png" /></td> <td valign="top" style="width:5%;text-align:right;">(10)</td></tr></table>

at the points of the curve <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940181.png" /> for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940182.png" />.

From (10) it follows that <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940183.png" />, where if <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940184.png" />, then the segment of the Hugoniot adiabatic (3) between the points <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940185.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940186.png" /> in the plane of the variables <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940187.png" /> is a segment of the line <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940188.png" />. The inequality (10) gives a condition of admissibility of the shock waves in case <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940189.png" /> alternates in sign. The basic properties of the shock waves in this case are modified: a) <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940190.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940191.png" />; and b) the shock waves lead to a non-decreasing entropy <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940192.png" />. There exist compressed shock waves <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940193.png" /> and rarefied shock waves <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940194.png" />.

==The origin of shock waves.==
In ideal media shock waves are found both as the result of the decomposition of the initial discontinuities into stable discontinuities (shock waves, contact and tangent discontinuities), as well as in the process of flow continuous up to a certain moment. E.g. in a one-dimensional isentropic flow of an ideal gas <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940195.png" /> shock waves do not arise if for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940196.png" /> 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/s/s084/s084940/s084940197.png" /></td> </tr></table>

hold, and, conversely, shock waves do arise if one of these conditions is violated.

In the theory of discontinuous solutions of quasi-linear hyperbolic equations and systems, a theory of discontinuities has been developed, similar to the mathematical theory of shock waves in gas dynamics. Basic results (far from complete) were obtained only for the case of two independent variables <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940198.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940199.png" />. One quasi-linear 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/s/s084/s084940/s084940200.png" /></td> <td valign="top" style="width:5%;text-align:right;">(11)</td></tr></table>

can be considered as a conservation law, so that <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940201.png" /> is the flow of the quantity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940202.png" />. The continuity of the flow through the line of discontinuity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940203.png" /> leads to the condition

<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/s/s084/s084940/s084940204.png" /></td> <td valign="top" style="width:5%;text-align:right;">(12)</td></tr></table>

which in analogy to (1) can be called the Hugoniot condition. Like for shock waves in gas dynamics, to (12) one adds the requirement of stability of the discontinuity. This guarantees the uniqueness of the solution of the Cauchy problem for (11) in the class of discontinuous functions. If <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940205.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940206.png" />, then the stability condition takes 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/s/s084/s084940/s084940207.png" /></td> </tr></table>

For quasi-linear systems of conservation laws,

<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/s/s084/s084940/s084940208.png" /></td> <td valign="top" style="width:5%;text-align:right;">(13)</td></tr></table>

the Hugoniot conditions have the form (12), where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940209.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940210.png" />. The problem of finding correct conditions of stability at the discontinuity has been solved only for a restricted class of systems of the form (13).

====References====
<table><TR><TD valign="top">[1]</TD> <TD valign="top">  S. Earnshaw,   ''Phil. Trans. Roy. Soc. London'' , '''150'''  (1860)  pp. 133–148  {{MR|}} {{ZBL|}} </TD></TR><TR><TD valign="top">[2]</TD> <TD valign="top">  B. Riemann,   "Ueber die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite" , ''Gesamm. math. Werke'' , Dover, reprint  (1953)  pp. 156–175  {{MR|}} {{ZBL|}} </TD></TR><TR><TD valign="top">[3]</TD> <TD valign="top">  W.J.M. Rankine,   ''Phil. Trans. Roy. Soc. London'' , '''160'''  (1870)  pp. 277–288  {{MR|}} {{ZBL|}} </TD></TR><TR><TD valign="top">[4]</TD> <TD valign="top">  H. Hugoniot,   ''J. École Polytechn.'' , '''58'''  (1889)  pp. 1–125  {{MR|1630342}} {{MR|1630341}} {{ZBL|}} </TD></TR><TR><TD valign="top">[5]</TD> <TD valign="top">  N.E. Kochin,  , ''Collected works'' , '''2''' , Moscow-Leningrad  (1949)  pp. 5–42  (In Russian)  {{MR|}} {{ZBL|}} </TD></TR><TR><TD valign="top">[6]</TD> <TD valign="top">  L.D. Landau,   E.M. Lifshitz,   "Fluid mechanics" , Pergamon  (1959)  (Translated from Russian)  {{MR|0108121}} {{ZBL|0655.76001}} {{ZBL|0146.22405}} </TD></TR><TR><TD valign="top">[7]</TD> <TD valign="top">  H. Weyl,   "Shock waves in arbitrary fluids"  ''Comm. Pure Appl. Math.'' , '''2'''  (1949)  pp. 103–122  {{MR|0034677}} {{ZBL|0035.42004}} </TD></TR><TR><TD valign="top">[8]</TD> <TD valign="top">  D. Gilbarg,   "The existence and limit behavior of the one-dimensional shock layer"  ''Amer. J. Math.'' , '''73'''  (1951)  pp. 256–274  {{MR|0044315}} {{ZBL|0044.21504}} </TD></TR><TR><TD valign="top">[9]</TD> <TD valign="top">  R. Becker,   ''Z. Phys.'' , '''8'''  (1922)  pp. 321–362  {{MR|2501325}} {{MR|1713813}} {{MR|0219505}} {{ZBL|}} </TD></TR><TR><TD valign="top">[10]</TD> <TD valign="top">  R. Courant,   K.O. Friedrichs,   "Supersonic flow and shock waves" , Interscience  (1948)  {{MR|0029615}} {{ZBL|0041.11302}} </TD></TR><TR><TD valign="top">[11]</TD> <TD valign="top">  L.I. Sedov,   "A course in continuum mechanics" , '''1–4''' , Wolters-Noordhoff  (1971–1972)  (Translated from Russian)  {{MR|}} {{ZBL|0308.73006}} {{ZBL|0284.73001}} {{ZBL|0241.76001}} {{ZBL|0308.73005}} </TD></TR><TR><TD valign="top">[12]</TD> <TD valign="top">  Ya.B. Zel'dovich,   "Theory of combustion and detonation of gases" , Moscow-Leningrad  (1944)  (In Russian)  {{MR|}} {{ZBL|}} </TD></TR><TR><TD valign="top">[13]</TD> <TD valign="top">  B.L. Rozhdestvenskii,   N.N. Yanenko,   "Systems of quasilinear equations and their applications to gas dynamics" , Amer. Math. Soc.  (1983)  (Translated from Russian)  {{MR|0694243}} {{ZBL|}} </TD></TR><TR><TD valign="top">[14]</TD> <TD valign="top">  O.A. Oleinik,   "Discontinuous solutions of non-linear differential equations"  ''Uspekhi Mat. Nauk'' , '''12''' :  3  (1957)  pp. 3–73  (In Russian)  {{MR|0094541}} {{ZBL|}} </TD></TR></table>



====Comments====
The short sentence concerning the use of the Boltzmann equation to describe the behaviour of the fluid dynamical variables in the zone of shock transition requires some clarification. The problem of the structure of a shock wave has been the subject of many studies [[#References|[a1]]]; for weak shocks (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940211.png" /> close to 1) the solution of the Boltzmann equation agrees with the solution of the Navier–Stokes equations for compressible fluids, but when <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940212.png" /> increases (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/s/s084/s084940/s084940213.png" />), the two solutions are widely different.

====References====
<table><TR><TD valign="top">[a1]</TD> <TD valign="top">  C. Cercignani,   "The Boltzmann equation and its applications" , Springer  (1988)  {{MR|1313028}} {{ZBL|0646.76001}} </TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top">  J. Smoller,   "Shock waves and reaction-diffusion equations" , Springer  (1983)  {{MR|0688146}} {{ZBL|0508.35002}} </TD></TR></table>