Nomenclature
-
a Time at which the applied pressure impulse ends.
-
ai Times that define the constant pressure squarewave segments.
-
Ai [V i–1 2 (Xai ) – γbi ]Xai 2
-
bi Weighting coefficients determined from f(t).
-
Bi γbi .
-
D Shock-wave propagation speed.
-
d LC Distance from a line charge.
-
f(t) Function describing the shape of a variable pressure impulse.
-
Hp Total momentum per unit area applied to the snow.
-
Hpi Momentum per unit area applied to the snow in the th square-wave segment.
-
Hs, dHs Momentum per unit area in the snow and its derivative.
-
I0 Instantaneously applied pressure impulse.
-
M(Xf ) Total mass of the compacted material.
-
m, dm Mass per unit area and its derivative.
-
Ρ Shock-wave pressure.
-
Ρ0 Maximum pressure amplitude of the shock wave.
-
P(t) Shock-wave pressure amplitude as a function of time.
-
r Current position of the inner radius of a cavity surface.
-
R Shock-wave propagation radius.
-
R0 Initial radius of a cavity surface in the snow.
-
t, dt Time and its derivative.
-
tf Time required for the shock wave to reach Xf.
-
t0 Time at which the pressure impulse is applied.
-
U0 Distance snow surface has moved.
-
V Snow-particle velocity.
-
Vi Snow-particle velocity in the ith square-wave segment.
-
Xa Shock-wave front position at t = a.
-
Χai Shock-wave front position at t =ai .
-
Xf Shock-wave front position; Xf = Xf−X0 for X0 =0.
-
Xi Current snow-surface position.
-
X0 Initial snow-surface position.
-
α Decay coefficient.
-
β Pressure-dependent relative snow compaction (1 – ρ 0/ρ f).
-
γ P 0 β/ρ f(1 – β).
-
ε R-R0 .
-
ρf Pressure-dependent compacted snow density.
-
ρ0 Initial snow density.
Introduction
An understanding of shock-wave attenuation (that is, the pressure-amplitude reduction of a shock wave over a given propagation distance) in snow is needed to solve a number of problems in planetary science, cold regions and military engineering, shock-wave isolation, and seismicmonitoring. For example, astronomers are interested in determining the stresses associated with impacts between low-density frost clouds and planetary bodies, avalanche-control professionals are interested in knowing the effectiveness of explosives in initiating avalanches, and military engineers arc primarily interested in the extent to which their explosive methods are degraded in snow-covered terrain. Shock-wave isolation, using snow, is a technique for protecting an object from shock damage by causing shock waves to attenuate in a protective layer of snow around the object.
The mechanical properties of snow for high strain-rate, large-amplitude shocks are not now well understood. Consequently, attempts to estimate shock-wave attenuation in snow have reliedon field measurements of pressure attenuation from explosions in snow and on theoretical constitutive descriptions.Reference Mellor,Mellor (1977) used the results from field measurements of explosive detonations in snow to estimate the attenuating properties of snow for spherical shock waves. He estimated that the attenuation of shock-wave pressure, combining geometric spreading and internal dissipation, near an explosion attenuates as R-4 (R being the propagation radius), and close to the elastic limit the attenuation is approximately an R−3 attenuation. Reference Brown,Brown (1980, Reference Brown,1981, Reference Brown,1983) estimated shock-wave attenuation for plane waves in snow. He used two theoretical volumetric constitutive laws, one for medium-to high-density snow and one for low-density snow, to estimate shock-wave attenuation. His results indicate that plane shock waves can attenuate by more than 80-90% after propagating through only 0.06–0.1 m of snow. These are large attenuations that arc, at their maximum, proportional to (X f – X 0)–1.2 ([X f – X 0]) is the propagation distance).
Some of the difficulties in determining shock-wave attenuation in snow are that the attenuation is affected by the pressure, geometry and duration of a shock wave, in addition to the mechanical properties of the snow. The objective of this paper is to use a simple analytical model to show the influence of a pressure impulse (the integral of pressure over the time interval during which the pressure acts) and its geometry on shock-wave pressure attenuation, in addition to the effects caused by the mechanical behavior of snow. The effects of shock-wave geometry are accounted for by calculating the attenuation of plane, cylindrical and spherical shock waves in snow. The effects of shock-wave duration and initial pressure amplitude are analyzed by applying a variable pressure-time profile shock impulse to the snow and calculating pressure as a function of shock-wave propagation distance.
Momentum Model
Shock-wave attenuation, in this model, is the result of a mechanical process in which attenuation occurs through momentum transfer from an applied pressure impulse to the snow. The momentum model, also known as the “snow-plow” model, was used in early attempts to develop constitutive relations to analyze the dynamic behavior of porous materials (Reference Herrmann,Herrmann, 1971). The porous material is assumed to compact to its final density at a negligible stress and then be incompressible (an ideal single-stage locking material). The snow-plow model can be extended to include porous materials that compact along a prescribed pressure-density (P–ρ) curve. Pressure changes in the porous material will cause it to collapse to a density consistent with the (P–ρ) curve. After compaction, the material is assumed to be incompressible until pressure changes cause a collapse to a new density state (an ideal multistage locking material (Reference Hanagud, and ŌuraHanagud, 1907)). The change in momentum caused by a pressure impulse is spread uniformly over the mass of the compressed material producing a constant pressure and particle velocity in the mass behind an advancing shock wave (Reference Herrmann,Herrmann, 1971). The shock wave is lengthened in time and reduced in amplitude as more of the material is compacted by the propagating shock. Hence, attenuation is caused by momentum spreading. Attenuation attributable to plastic deformation, fracturing and release waves are not considered. Although plastic deformation and fracturing may significantly affect the behavior of snow, the compaction is treated phenomenologically since little data exist from which to determine the mechanisms that produce a given compaction curve. Release waves (tensile waves generated at boundaries or following a pressure impulse) are also an important source of attenuation. The momentum model, presented here, is not capable of treating attenuation caused by release waves. A following paper using experimentally derived compaction data for snow will treat the more complicated problem of release-wave attenuation. In this paper, attenuation solutions are first developed for a single-stage locking material and then extended to include a multi-stage locking material.
Instantaneously Applied Impulse
Consider a pressure impulse (I0 ) applied normal to the plane surface of snow, and assume that snow is an ideal single-stage locking material. Snow next to the plane surface will be immediately compacted to its defined final density stage. Since the compacted snow is rigid, it will move at a uniform pressure and particle velocity after the pressure impulse is applied. The stress wave will propagate into the snow as a compaction shock wave, moving with a velocity D at a pressure P, and particle velocity V. At the shock front, these parameters are related by the Rankine Hugoniot jump conditions for the conservation of mass and momentum across the shock front (Reference Kolsky,Kolsky, 1963)
and
where ρ 0 is the initial density of the snow and ρ f is the compacted snow density.
Figure 1 shows the deformation geometry in one dimension for a pressure impulse applied at the initial snow surface X0. The location of the surface of the snow and the shock front at some time after the application of the pressure impulse are Xi and Xf respectively. During the time that the shock front has traveled to Xf, the snow surface has moved a distance U0. U0 is a function of time or,alternatively, a function of shock-wave propagation position Xf . At any time, the location of the snow surface is
The displacement of the snow surface can be calculated by integrating the particle velocity over time using
where t0 is the time that the pressure impulse is applied and tf denotes the time required for the shock wave to reach Xf . The equation describing shock-wave attenuation can be formulated as a function of time or shock-wave propagation distance. Knowing the shock-wave pressure attenuation as a function of propagation distance is of more practical interest than following the pressure as a function of time. Shock-wave attenuation as a function of propagation distance can be developed by reformulating the snow-surface displacement in terms of shock-wave position rather than time by using
and rewriting Equation (4) as
From Equation (1)
where β describes the relative snow compaction (the mechanical behavior of the snow). Substituting Equation (7) into Equation (6) and solving for the position of the surface at time tf as a function of X0 and Xf gives
The mass and the momentum per unit area of the snow contained in a length dϰ; between the current snow-surface position Xi and shock-front position Xf are
and
Since the compacted snow is assumed to be rigid, V and ρ f are constant through the region Xi to Xf where V is equal to the particle velocity at Xf , i.e. V = V(Xf). Consequently, the momentum per unit area in the snow at the time the shock wave has reached Xf is
Substituting Equation (8) into Equation (11) gives
The momentum per unit area in the snow Hs must be equal to the momentum per unit area applied to the snow from the instantaneously applied pressure impulse that is given by
Equating Hs and Hp and solving for V gives
Using Equations (7) and (14) in Equation (2) to determine the pressure at the shock front gives
Figure 2 shows the shock-front pressure for snow subjected to an instantaneous pressure impulse and compaction to a final density occurring in one step (single stage). The pressure attenuates as (Xf – X0 )−2, with proportionality constant determined by pressure impulse, initial snow density and relative snow compaction. The proportionality constant determines the pressure amplitude that the snow is capable of supporting at a given (Xf – X0 ), but not the attenuation with distance. Increasing the pressure-impulse amplitude produces a greater momentum transfer to snow that results in higher shock-wave pressures at a given propagation distance, ρ 0 and β remaining constant. A larger β with constant I0 and ρ 0 causes a lower shock-wave propagation velocity at a given (Xf – X0), as required by mass-conservation Equation (7), which results in reduced pressures in the snow.
When snow compaction is pressure-dependent (Fig. 3), the momentum equation (Equation (11)) cannot be directly integrated as ρ f is not constant through the region Xi to Xf . Since the material behind the shock wave is rigid, V is constant from Xi to Xf and the momentum (Equation (11)) is
where
is the mass of the compacted material behind the shock wave. Although the density varies over the limits of integration, it is only the total mass behind the shock wave that is needed to determine the momentum. The mass can be determined from
where ρf(Xf ) is the final density at the shock front and Xi is determined from Equation (8)). The pressure and ρf(Xf ) at the shock front are determined by solving Equation (15) iteratively until Pf(Xf ) and ρf(Xf ) lie on the specified (Ρ-ρ) compaction curve (Fig. 3). Shock-pressure attenuation is approximately (X f – X 0)−1.5 for the prescribed multi-stage compaction (P–ρ) curve used in this study as compared to (X f – X 0)−2 for single-stage compaction (impulse Β compared to impulse A; (Fig. 4). The pressure attenuation with distance is a direct result of the form of the compaction curve and may be different for various snow types. A single-stage compaction will always have a larger attenuation with distance than a multi-stage compaction. For single-stage compaction, ρ f will remain constant as the shock pressure decreases whereas for multi-stage compaction ρ f decreases with a decrease in pressure.
Constant-Pressure Impulse of Finite Duration
The solution given by Equation (15) is not a very satisfying way of evaluating shock-wave attenuation, since it assumes instantaneous application of the pressure impulse. This is unrealistic, as it implies infinite pressure at x = X0 and does not show how pressure impulses with different amplitudes and durations, but each having the same total momentum, are attenuated. An analytical solution can be derived for a constant-pressure impulse of finite duration, given by
where a is the duration of the applied pressure pulse. The momentum per unit area caused by the pressure impulse is
and
Equations (17) and (18) show that the momentum during application of the pressure impulse will vary with time, while the momentum after the impulse has been applied is constant. Therefore, separate solutions for shock-wave pressure attenuation are needed during the time period of pulse application and for the time period after the pulse has been applied. Equations (17) and (18) can be transformed into spatial coordinates by use of the results of Equations (5) and (7), giving
The limits of integration are 0 to X a = aV/β. Equating momentum in the snow (Equation (12)) to the pressure-impulse momentum (Equation (19)), letting X0 = 0 and differentiating gives
Integrating Equation (20) gives the particle velocity
Hence, during the period of pressure impulse application, the density, relative snow compaction and particle velocity are constant and there is no attenuation of the shock-wave pressure. The distance that the shock wave has traveled at t = a is given by
After the pressure-impulse momentum has been applied to the snow, the relationship between the pressure-impulse momentum and snow momentum is
Solving for the particle velocity and using the definition of Xa in Equation (22) gives
The shock-wave pressure as a function of distance can now be determined from Equations (2), (7), (21) and (24), giving
Since Ρ is constant while the pressure impulse is being applied, ρ f and β also constant and Xa can be directly calculated for either single-stage or multi-stage compaction. After the pressure impulse has ended, solving Equation (25b) for multi-stage compaction requires that ρf(Xf ) be iterated until ρf(Xf ) and P(Xf ) lie on the (P-p) curve (Fig. 3), as was the case for the instantaneously applied pressure impulse. Figure 4 shows the shock-wave attenuation for single-stage compaction (square A) and multi-stage compaction (square B). The shock wave does not begin attenuating until after all of the pressure-impulse momentum has been applied to the snow surface, and then it attenuates as for for the single-stage compaction and approximately as for the multistage compaction curve. The propagation distance at a constant shock-wave pressure differs between squares A and B. This difference is controled by the amount of relative snow compaction. With a large amount of compaction, the shock wave propagates a shorter distance before it attenuates than it does for a smaller amount of compaction. Comparison of the square-wave pressure-impulse (squares A and B) to an instantaneous pressure impulse (impulses A and B) shows that the pressure is lower for the square wave during its duration of application than for an instantaneously applied pressure impulse. After the square-wave momentum has been completely transferred to the snow, both the square wave and instantaneously applied impulse pressures are the same. Changing the duration of a square-wave pressure impulse with constant total momentum will also affect shock-wave pressure attenuation. Decreasing the duration of an impulse causes the initial pressure of a shock wave to increase. The shock wave will also being attenuating at , for single-stage compaction, and , for multi-stage compaction, earlier.
Variable Pressure Impulse of Finite Duration
Simulating a realistic applied pressure impulse requires a function that can represent variable pressure impulses of finite duration. These variable-pressure impulses must be modeled by non-linear differential equations that do not have closed-form solutions and whose numerical solutions can be unstable (Reference Hornbeck,Hornbeck, 1971). In addition, the form of the differential equations depends on the form of the pressure impulse applied to the snow surface. To reduce these computational difficulties, here the variable pressure impulses were approximated by a sequence of square waves. The accuracy of such an approximation depends on the duration specified for each square-wave segment comprising the total applied pressure impulse. Instead of exactly formulating the problem and obtaining a numerical solution, the pressure impulse is approximated and solved analytically within the sequential time steps. The variable-pressure impulse of finite duration beginning at t = 0 is approximated by
A smooth curve, f(t), fitted to the sequence of square waves passes through a series of points, bi (weighting coefficients), at the midpoint between the beginning time and ending time for each square-wave segment f[(a i+1 + ai )/2], where each square-wave segment ends at time ai where i = 1,2,3, …, n – 1 (Fig. 5)).
Solving for the pressure in snow as a function of distance, using Equation (26), requires that the problem be formulated for three conditions. The first is for the time, or spatial increment, of the first square-wave segment. The next solution must account for the momentum that has been applied to the snow at the time, or spatial position, of interest. Finally, a solution must be found after all of the pressure-impulse momentum has been applied to the snow. The conservation of momentum conditions, for X0 = 0, can be expressed as
There are n + 1 equations to solve before the complete particle-velocity solution can be found. Each equation depends on the solutions of the previous segments.
Equation (27a) is a constant-pressure impulse with a particle velocity of
where
Equation (27b) can be solved sequentially to give a set of recursion relations for position, and particle velocity
where
and
After all of the pressure impulse has been applied to the snow, the particle velocity is
where
is the total momentum per unit area applied to the snow. Equations (28) and (29) are used to start the solution and the recursion relations (Equations (30) and (31)) are used to follow the progression of the shock wave. Equation (32) is used to calculate particle velocities after the pressure impulse has ended.
Figure 4 (exponentials A and B) shows the results of using Equation (28) through (32), (2) and (7) to calculate the shock-wave attenuation for an exponentially decaying pressure impulse given by
where α is a decay constant and P0 is the initial pressure, and f (t) = e–αt . Exponential A is the shock-wave pressure attenuation resulting from a single-stage compaction to ρ f and exponential Β is the attenuation that results from a multi-stage compaction. The exponential A and Β shock waves each begin attenuating immediately upon application of the exponential pressure impulse. The attenuation for exponential A gradually increases to as the pressure impulse ends. Attenuation of exponential Β is about after the pressure impulse ends, which is consistent with the results for instantaneous and finite-duration impulse and a multistage compaction. As with the constant-pressure impulse of finite duration, the attenuation of shock waves generated by finite-duration exponential pressure impulses depends on impulse duration. A short-duration impulse will produce a higher initial shock-wave pressure than an impulse with the same total momentum content but of longer duration.
Attenuation of Cylindrical and Spherical Geometry Shock Waves
Geometric spreading can greatly increase shock-wave attenuation in snow. Reference Torvik,Torvik (1971) developed equations describing pressure attenuation for cylindrical and spherical geometry shock waves using a single-stage compaction snow-plow model and assuming instantaneous application of pressure impulse. For a cylindrically spreading shock, the pressure is given by
The pressure for a spherically spreading shock wave is
where
I0 is the instantaneously applied pressure impulse, R0 is the initial radius of the cavity surface in the snow on which the pressure impulse is applied, R is the shock-wave propagation radius and r is the location of the inner radius of the cavity at some time after the application of the pressure impulse (Fig. 6).
When R = R 0 + ϵ, where ϵ ≪ R 0, a cylindrical geometry shock wave attenuates as ϵ –2 for single-stage compaction and about ϵ –1.5 for the multi-stage compaction curve shown in Figure 3, increasing to an attenuation of R –4 for ϵ ≫ R 0 (Fig. 7). A spherical geometry shock wave attenuates as ϵ –2 for single-stage compaction and about ϵ –1.5 for multi-stage compaction when ϵ ≪ R 0, increasing to R –6 for ϵ ≫ R 0 (Fig. 8). Pressure attenuation, for cylindrical and spherical shock waves, varies significantly depending on the initial radius over which the pressure impulse is applied and the distance from the initial radius.
Pressure Attenuation of a Line Charge on Snow
Line charges are used to clear minefields and to remove snow cornices from mountain ridge tops. Therefore, it is of practical interest to know the extent that snow attenuates shock waves produced by line-charge detonations. For a line charge resting on or above a snow cover, the cylindrical pressure wave that hits the snow surface propagates primarily through the air.
Reference Ford,Ford (1986) has measured air-blast pressure, positive-phase duration and total pressure impulse at the ground as a function of lateral distance from the axis of a line-charge explosion (Table 1). The positive-phase duration is defined as the duration of the compressive (positive) shock-wave pressure. It is assumed, for this example, that the air-blast pressure given in Table 1 would also occur over snow. Furthermore, the snow is assumed to be deep and shock-wave pressures are determined only for small propagation depths into the snow so that ϵ ≪ R 0, where R0 is the distance from the line charge to the point of interest on the snow surface. With these conditions, the snow-ground boundary can be neglected and shock-wave pressure in the snow can be estimated from the equations for a one-dimensional variable-pressure impulse of finite duration (Equations (28)–(32), (2) and (7)), and an exponentially attenuating pressure pulse (Equation (33)). Figure 9) shows the shock-wave pressure (air blast) as a function of distance from the line charge (dLC ) along the snow surface, using the pressure-impulse data from Table 1. Figure 9 also shows the calculated pressures as a function of (dLC ) for shock-wave propagation depths into the snow of 0.06, 0.2, 0.4 and 0.6 m calculated using the compaction curve given in Figure 3.
Air-blast shock-wave pressures do not attenuate the same amount for a given depth in the snow for different values of (dLC ). Shock-wave attenuation is greatest near the line charge (dLC ) 0.96 m) and least at the farthest distance (dLC ) 23.97 m). This occurs because the duration of the applied air-blast pressure impulse (positive-phase duration) increases with (dLC ) (Table 1). The increasing positive-phase duration allows the shock wave to propagate through a greater depth of snow before starting the ϵ –1.5 attenuation. These findings suggest that, for shallow shock-wave penetration depth into a snow cover, line-charge pressure reduction depends on the positive-phase duration in addition to the total pressure impulse. Also, shock-wave pressure attenuation in a shallow snow cover can be reduced or eliminated by in-creasing the positive-phase duration of a given pressure impulse.
Discussion and Conclusions
The snow-plow model has a limited ability to describe shock-wave propagation and attenuation in snow because of the simplifying underlying assumptions. The model is, however, capable of illustrating some of the important features and giving conservative estimates of shock-wave attenuation in snow. Shock-wave pressure attenuation predicted by the snow-plow model will be less for a multi-stage snow compaction than for a single-stage compaction to the final snow density.
The snow-plow model predicts that maximum shock-wave pressure attenuations will proceed as a function of for plane waves (where X0 = 0), R−4 for cylindrical waves, and R−6 for spherical waves with a single-stage snow compaction. Attenuation for multi-stage compaction depends on the form of the (P–ρ) curve used to describe snow compaction and ranges from a function of , for the compaction curve used in this study, to the maximum for plane, cyclindrical and spherical shock waves.
Comparison of instantaneous, finite-duration square-wave and finite-duration exponential pressure impulses (each with the same total momentum) shows that shock-wave attenuation can be markedly different during the period of pressure-impulse application. The instantaneously applied pressure impulse applies all of its mom-entum at once, resulting in infinite initial pressure that immediately begins attenuating as for single-stage compaction and as for multi-stage compaction. The square-wave presure impulse maintains a constant pressure equal to the applied pressure until all of its momentum has been transferred to the snow and then begins to attenuate. Finally, the exponential pressure impulse gradually attenuates while its momentum is being transferred to the snow, asymptotically approaching an attenuation for single-stage compaction and an attenuation for multi-stage compaction. Once a pressure impulse has been applied to the snow, pressure attenuation is controled by the mechanical properties of the snow and the magnitude of total pressure impulse.
Reference Mellor,Mellor (1977) estimated that spherical geometry shocks, where ϵ ≫ R 0, attenuate as R−3 to R-4 as compared to the snow-plow model prediction of R-6. It may be that the snow-plow model predictions are grossly in error because of the simple snow-compaction path used in the model or that the field measurements used by Mellor to make his estimates of attenuation were not the result of shock-wave propagation in snow. The pressure sensors used in many of the field tests could not survive near-source pressures and were, consequently, located well outside the zone of shock-wave-induced snow compaction. Outside this zone of snow compaction, pressure attenuation is due to viscous dissipation and geometric spreading that produce much less pressure attenuation than the pore-collapse mechanisms described by the snow-plow model. If, in fact, this was the case, Mellor’s estimates of pressure attenuation in snow for spherical shocks may be too low, but a final conclusion is not possible until better snow-compaction and release paths are used in the model calculations, or more field tests are conducted.
Estimating shock-wave attenuation requires some thoughtful application of the results of this study. Explosive-induced shocks in snow often result from an explosive-charge detonation on the snow or in the air above the snow. In this situation the shock wave will be transmitted into the snow primarily by an air-blast wave propagating over the snow surface. If the initial radius for the air-blast pressure wave is large compared to the snow depth, then the shock wave may be approximated by a plane wave rather than a spherical wave or cylindrical wave. Attenuation will then proceed as ϵ –1.5 rather than as R−6 for spherical waves or R−4 for cylindrical waves.
Pressure attenutation of shock waves in snow can be reduced, delayed or even eliminated by using a sufficiently long-duration pressure impulse. Shock-wave propagation and attenuation in snow are affected by shock-wave geometry, pressure-impulse time variation and the mechanical properties of snow. Consequently, determining the mechanical response of snow to shock loading requires an accurate knowledge of the shock-loading conditions in addition to a knowledge of stresses and particle velocities in the snow.
Acknowledgements
I thank Dr.J. A. Brown (Los Alamos National Laboratory), Dr E. S. Gaffncy (Ktech Corporation) and K. Jones (CRREL) for technically reviewing this paper. The research described in this paper was funded by DA Project 4A7G2784AT42, “Cold regions engineering technology ” Work unit CS/012, “Attenuation of shock waves by snow”.
The accuracy of the references in the text and in this list is the responsibility of the authors, to whom queries should be addressed.