On the Decay of Plane Shock Waves
Laboratory Report No. 423
Aberdeen Proving Ground, Md.
8 November 1943
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. . 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
(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 Rankine-Hugoniot equation
where y denotes the ratio of pressures on either side of the shock front, is compared with that given by
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 = 5c-u.
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. . 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)
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
The foregoing equations are equivalent to
The particular significance of the case
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.
and equation (9) reduces to
Aof the foregoing equation can be written down at once. We have
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))
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
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 Umax behind the shock front by the Rankine-Hugoniot equation
Thus at any given instant the positive phase of the pulse extends from
(24)x = t to .
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
Accordingly the value of u at the new position B′C′ of the shock front is given by
The increment in umax during the time Δt is therefore given by
Thus the differential equation governing the dependence of U′ on time is
This equation can be integrated to give
where U0 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)
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 U0 ≤ 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 U0 = 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 xmax, umax, and U as functions of time also for the case U0 = 1.24 and γ = 1.4.
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
It is possible that the laws
are more general than their derivations from equation (31) and (32) would suggest.
3. . 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, . pp. 375-380) 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))
where a1 and a2 are two arbitrary constants, the general integral of equation (16) can be expressed as the eliminant between the equations
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.
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
where p denotes the pressure, ρ the density, and u the mass velocity at any point. Introduce the two functions
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
(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))
Remembering that θ remains constant during the motion apd using the equations (1) we find that
Now for a perfect gas (with a ratio of specific heats γ) we have
where A is a constant. From this equation we readily derive the formulae
Thus for the case under consideration
(12) and .
Moreover we readily verify from the foregoing that
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.