On the Decay of Plane Shock Waves
Ballistic Research
Laboratory Report No. 423
Chandrasekhar/emh
Aberdeen Proving Ground, Md.
8 November 1943
ON THE DECAY OF PLANE SHOCK WAVES
Abstract
In this import the problem of the decay of plane shock waves is considered and it is shown that there exists a special class of shock pulses for which the problem of decay admits of an explicit solution. The case in question arises when shocks of moderate intensities (Mach numbers less than 1.5 if an accuracy of one per cent is demanded) are quantity 5 c — u (where c denotes the local velocity of sound and u the mass velocity) as the shock front is crossed can be neglected. The case of linear shock pulses in which both u and c are linear functions of x behind the shock front is considered particularly in some detail.
1. Introduction. The only case in the theory of the decay of plane shock waves which appears to have been studied in any detail is the one which has recently been investigated by W. G. Penney and H. K. Dasgupta (R.C. 301) by numerical methods. And it appears not to have been recognized that there exists a special class of shock pulses for which the problem of decay admits of an explicit solution. The case in question arises when shock waves of moderate intensities (i.e., shocks with Mach numbers less than 1.5) are considered. For, under these circumstances (when the ratio of pressures on either side of the shock front is less than 2.5) the increase in entropy of an element of gas as it crosses the shock front can be ignored. And moreover the change in the quantity
(1)
(where c denotes the local velocity of sound and u the mass velocity) as we cross the shock can also be neglected. That this is so is apparent from Table 1 where the velocity of sound immediately behind the shock front (in units of the velocity of sound, a, in front of the shock) as determined by the RankineHugoniot equation
(2)

where y denotes the ratio of pressures on either side of the shock front, is compared with that given by
(3)

TABLE 1
In this table y denotes the ratio of the pressures on the two sides of the shock front, U the shock velocity, u the mass velocity behind the front (both in a frame of reference in which the mass velocity in the undisturbed region is zero), c the velocity of sound immediately behind the shock front, P = 5c+u and Q = 5cu.
y  U/a  u/a  c/a  y^{1/7}  P/a  Q/a 
1.0  1.000  0.000  1.000  1.000  5.000  5.000 
1.1  1.042  0.069  1.014  1.014  5.137  5.000 
1.2  1.082  0.132  1.026  1.026  5.265  5.001 
1.3  1.121  0.191  1.038  1.038  5.384  5.001 
1.4  1.159  0.247  1.050  1.049  5.496  5.002 
1.5  1.195  0.299  1.061  1.059  5.602  5.005 
1.6  1.231  0.348  1.071  1.069  5.704  5.007 
1.7  1.265  0.395  1.081  1.078  5.801  5.010 
1.8  1.298  0.440  1.091  1.087  5.894  5.014 
1.9  1.331  0.483  1.100  1.095  5.984  5.018 
2.0  1.363  0.524  1.109  1.104  6.071  5.023 
2.1  1.394  0.564  1.118  1.112  6.155  5.027 
2.2  1.424  0.602  1.127  1.119  6.237  5.034 
2.3  1.454  0.639  1.136  1.127  6.318  5.040 
2.4  1.483  0.674  1.144  1.133  6.395  5.047 
2.5  1.512  0.709  1.152  1.140  6.470  5.053 
It would accordingly appear that a significant case of shock pulses is provided by neglecting entropy changes and considering Q as a constant throughout. As we shall see this leads to an interesting class of shock solutions which does not appear to have been isolated so far.
2. The Solutions for a Special Class of Shock Pulses for which Q = Constant. Measuring the velocities in units of the velocity of sound in the undisturbed air in front of the shock, the equations of motion in Riemann's form are (cf. Penney, R. C., 260; also, the Appendix to this paper)
,(4)
and
,(5)
where
(6) ;
and Q denotes the potential temperature (i.e., the temperature which the element of gas under consideration would have if reduced adiabatically to a certain standard pressure). For shocks of moderate Intensities the term in 0 which, incorporates the changes in entropy can be ignored and we have
(7) ,
and
(8) ,
The foregoing equations are equivalent to
(9) ,
and
(10) ,
The particular significance of the case
(11) constant
(12) for
It appears that after a sufficient length of time the shock pulse must be characterized by Q = 5 throughout. In any event it is clear that in considering the case Q = constant = 5 we are not limiting ourselves to too 'trivial' a class of shock pulses. Assuming then the validity of equation (11) we have.
(13) ;;,
and equation (9) reduces to
(14) .
Letting
(15)
we have
(16) .
A complete integral of the foregoing equation can be written down at once. We have
(17)
where b and q are two arbitrary constants. In terms of this complete integral we can readily write down the general solution of equation (16). But postponing this discussion to §3 and limiting ourselves for the present to the solution (17), we have (cf. eqs. (13) and (15))
(18)
and
(19)
According to equations (13) and (19), at any given instant both c and u are linear functions of x behind the shock front. By a translation of the time axis we can rewrite equations (13) and (19) more conveniently in the forms
(20) ,
and
(21) .
From the foregoing equations it follows that
(22) u = 0, c = 1 for x = t
In other words the point at which u = 0 in the pulse moves fort ward with a uniform velocity equal to the velocity of sound in the undisturbed regions. Starting at this point, x = t, u and c increase linearly with x till they attain their maximum values immediately behind the shock front. And moreover the shock velocity U is related to the mass velocity U_{max} behind the shock front by the RankineHugoniot equation
(23) .
Thus at any given instant the positive phase of the pulse extends from^{[1]}
(24) x = t to .
Fig. 1
Let ABC denote the pulse at time t. At time t + Δt, A moves to A′ and C to C′ where
(25) and ,
where U denotes the shock velocity at time t. During the same time the slope of AB also changes. For at t + Δt the velocity field will he governed by
(26)
Accordingly the value of u at the new position B′C′ of the shock front is given by
(27) .
The increment in u_{max} during the time Δt is therefore given by
(28)

(29) .
Thus the differential equation governing the dependence of U′ on time is
(30) .
This equation can be integrated to give
(31)
where U_{0} denotes the shock velocity at time t = 0. The dependence of the ratio of pressures on the two sides of the shock front on time can, be readily written down from equation (31). We have (for γ = 1.4)
(32) .
Equations (20), (21), (23), (24), (31) and (32) together describe completely the behavior of a linear shock pulse. The only limitation on this solution is that U_{0} ≤ 1.5 for γ = 1.4 if an accuracy of the order of 1% is demanded.
In Fig. 2 we have illustrated the dependence of U on t for the case U_{0} = 1.24 and γ = 1.4. Similarly in Fig. 3 the velocity field in the positive phase of the shock pulse at various instants is illustrated for the same case. And finally in Table II we have tabulated x_{max}, u_{max}, and U as functions of time also for the case U_{0} = 1.24 and γ = 1.4.
FIG. 2
FIG. 3
shock pulse with a Mach number in the neighborhood of 3, it is seen that their curves illustrating the velocity in the pulse at various instants are very similar to those illustrated in Fig. 3.
We may further note that, according to equations (31) and (32), for γ = 1.4
(33) 
It is possible that the laws
(34) 
are more general than their derivations from equation (31) and (32) would suggest.
Table II
qt  qx_{max}  U  u_{max}  qt  qx_{max}  U  u_{max} 
0  0.434  1.24  0.361  004.878  006.000  1.10  0.159 
0.195  0.673  1.22  0.334  008.199  009.616  1.08  0.128 
0.450  0.982  1.20  0.306  015.375  017.28  1.06  0.097 
0.796  1.394  1.18  0.277  035.88  038.77  1.04  0.065 
1.280  1.960  1.16  0.248  146.6  152.5  1.02  0.033 
1.986  2.771  1.14  .219  ∞  ∞  1.00  0.000 
3.074  3.999  1.12  0.189 
3. On the General Solution of ShockPulses with Q = Constant. In §2 we have discussed the special case of linear shock pulses which are characterized by Q = Constant throughout. In this section we shall briefly indicate how the most general shock pulses under the circumstances Q = constant can be constructed. For this purpose we require the general integral of equation (16). As is well known (cf. A.R. Forsyth, Differential Equations. pp. 375380) the general integral can be readily written down in terms of a complete integral, i.e., an integral which contains as many constants as independent variables. Writing the complete integral of equation (16) in the form (cf. eq. (17))
(38)
where a_{1} and a_{2} are two arbitrary constants, the general integral of equation (16) can be expressed as the eliminant between the equations
(36)

where χ is any arbitrary function of a. It is now evident that with the solution in the form (36) we can make φ satisfy any arbitrary distribution at time t = 0. Alternatively we may say that the distribution of c (or equivalently u) at time t = 0 will determine χ(a) thus making the solution determinate. In this fashion the most general form of shock pulses under the assumptions made in §1 can be constructed. In a later report we propose to give examples of shock pulses belonging to this more general class.
APPENDIX
THE EQUATIONS OF MOTION IN RIEMANN'S FORM ALLOWING FOR CHANGES IN ENTROPY
In the text equations (4) and (5) were quoted from Penney (R.C., 260). Since Penney's report may not be generally accessible, it has been thought worth while to include here a brief outline of his derivation of these equations.
The equations of motion of a linear pulse are
,
(1)
,
where p denotes the pressure, ρ the density, and u the mass velocity at any point. Introduce the two functions
,
(2)
,
where
(3) ;,
and c the local sound velocity. In equations (3), θ denotes the temperature which the element of gas under consideration would have when it is reduced to a standard pressure adiabatically. It is evident that θ remains constant during the motion of any element of gas (except when it crosses a discontinuity).
Consider the total differential
(4) (i=1,2).
(4′)
(the + or  sign going respectively with i = 1 or 2) we shall consider the second term in square brackets occurring as the coefficient of dt. We have (cf. eq. (3))
(5)
.
Remembering that θ remains constant during the motion apd using the equations (1) we find that
(6) .
Thus,
(7) .
and
(7) .
Now for a perfect gas (with a ratio of specific heats γ) we have
(9) .
where A is a constant. From this equation we readily derive the formulae
(10) ,
and
(11) .
Thus for the case under consideration
(12) and .
Moreover we readily verify from the foregoing that
(13) .
Combining equations (7), (12) and (13) we obtain the equations quoted in the text.
 ↑ In this paper we do not explicitly discuss the negative (or the suction) phase of the pulse. It is however clear that the discussion of the suction phase will proceed on lines exactly similar to that given for the positive phase.

This work is in the public domain in the United States because it is a work of the United States federal government (see 17 U.S.C. 105).