A new paradigm is introduced for modeling reactive shock waves in heterogeneous solids at the continuum level. Inspired by the probability density function methods from turbulent reactive flows, it is hypothesized that the unreacted material microstructures lead to a distribution of heat release rates from chemical reaction. Fluctuations in heat release, rather than velocity, are coupled to the reactive Euler equations which are then solved via the Riemann problem. A numerically efficient, one-dimensional hydrocode is used to demonstrate this new approach, and simulation results of a representative impact calculation (inert flyer into explosive target) are discussed.
A novel multiphase shock tube has been constructed to test the interaction of a planar shock wave with a dense gas-solid field of particles. The particle field is generated by a gravity-fed method that results in a spanwise curtain of 100-micron particles producing a volume fraction of about 15%. Interactions with incident shock Mach numbers of 1.67 and 1.95 are reported. High-speed schlieren imaging is used to reveal the complex wave structure associated with the interaction. After the impingement of the incident shock, transmitted and reflected shocks are observed, which lead to differences in flow properties across the streamwise dimension of the curtain. Tens of microseconds after the onset of the interaction, the particle field begins to propagate downstream, and disperse. The spread of the particle field, as a function of its position, is seen to be nearly identical for both Mach numbers. Immediately downstream of the curtain, the peak pressures associated with the Mach 1.67 and 1.95 interactions are about 35% and 45% greater than tests without particles, respectively. For both Mach numbers tested, the energy and momentum fluxes in the induced flow far downstream are reduced by about 30-40% by the presence of the particle field.
A novel multiphase shock tube has recently been developed to study particle dynamics in gas-solid flows having particle volume fractions that reside between the dilute and granular regimes. The method for introducing particles into the tube involves the use of a gravity-fed contoured particle seeder, which is capable of producing dense fields of spatially isotropic particles. The facility is capable of producing planar shocks having a maximum shock Mach number of about 2.1 that propagate into air at initially ambient conditions. The primary purpose of this new facility is to provide high fidelity data of shock-particle interactions in flows having particle volume fractions of about 1 to 50%. To achieve this goal, the facility drives a planar shock into a spatially isotropic field, or curtain, of particles. Experiments are conducted for two configurations where the particle curtain is either parallel to the spanwise, or the streamwise direction. Arrays of high-frequency-response pressure transducers are placed near the particle curtain to measure the attenuation and shape change of the shock owing to its interaction with the dense gas particle field. In addition, simultaneous high-speed imaging is used to visualize the impact of the shock on the particle curtain and to measure the particle motion induced downstream of the shock.
A novel multiphase shock tube to study particle dynamics in gas-solid flows has been constructed and tested. Currently, there is a gap in data for flows having particle volume fractions between the dusty and granular regimes. The primary purpose of this new facility is to fill that gap by providing high quality data of shock-particle interactions in flows having dense gas particle volume fractions. Towards this end, the facility aims to drive a shock into a spatially isotropic field, or curtain, of particles. Through bench-top experimentation, a method emerged for achieving this challenging task that involves the use of a gravity-fed contoured particle seeder. The seeding method is capable of producing fields of spatially isotropic particles having volume fractions of about 1 to 35%. The use of the seeder in combination with the shock tube allows for the testing of the impingement of a planar shock on a dense field of particles. The first experiments in the multiphase shock tube have been conducted and the facility is now operational.
Veloce is a medium-voltage, high-current, compact pulsed power generator developed for isentropic and shock compression experiments. Because of its increased availability and ease of operation, Veloce is well suited for studying isentropic compression experiments (ICE) in much greater detail than previously allowed with larger pulsed power machines such as the Z accelerator. Since the compact pulsed power technology used for dynamic material experiments has not been previously used, it is necessary to examine several key issues to ensure that accurate results are obtained. In the present experiments, issues such as panel and sample preparation, uniformity of loading, and edge effects were extensively examined. In addition, magnetohydrodynamic (MHD) simulations using the ALEGRA code were performed to interpret the experimental results and to design improved sample/panel configurations. Examples of recent ICE studies on aluminum are presented.
A new approach to explosive sample preparation is described in which microelectronics-related processing techniques are utilized. Fused silica and alumina substrates were prepared utilizing laser machining. Films of PETN were deposited into channels within the substrates by physical vapor deposition. Four distinct explosive behaviors were observed with high-speed framing photography by driving the films with a donor explosive. Initiation at hot spots was directly observed, followed by either energy dissipation leading to failure, or growth to a detonation. Unsteady behavior in velocity and structure was observed as reactive waves failed due to decreasing channel width. Mesoscale simulations were performed to assist in experiment development and understanding. We have demonstrated the ability to pattern these films of explosives and preliminary mesoscale simulations of arrays of voids showed effects dependent on void size and that detonation would not develop with voids below a certain size. Future work involves experimentation on deposited films with regular patterned porosity to elucidate mesoscale explosive behavior.
Explosive initiation and energy release have been studied in two sample geometries designed to minimize stochastic behavior in shock-loading experiments. These sample concepts include a design with explosive material occupying the hole locations of a close-packed bed of inert spheres and a design that utilizes infiltration of a liquid explosive into a well-defined inert matrix. Wave profiles transmitted by these samples in gas-gun impact experiments have been characterized by both velocity interferometry diagnostics and three-dimensional numerical simulations. Highly organized wave structures associated with the characteristic length scales of the deterministic samples have been observed. Initiation and reaction growth in an inert matrix filled with sensitized nitromethane (a homogeneous explosive material) result in wave profiles similar to those observed with heterogeneous explosives. Comparison of experimental and numerical results indicates that energetic material studies in deterministic sample geometries can provide an important new tool for validation of models of energy release in numerical simulations of explosive initiation and performance.
New experimental diagnostics and computational modeling provide an unprecedented means for improving the understanding of energetic material behavior at the mesoscale (grain or crystal ensemble levels). This study focuses on the determination of appropriate constitutive and EOS property data of the constituents of an energetic composite at high stress and moderate strain-rate states. The Sandia Z accelerator is used to determine the mechanical response of energetic composites via isentropic ramp wave compression loading. In this paper we describe an energy source method in CTH that models ramp loading for the analysis of ICE experiments. This approach is applied to design experimental configurations to probe the constituent response of PBX 9501 subjected to {approx}40 Kbar ramp load over 300 ns duration. Multiple VISAR are used to determine the averaged response of the composite material in comparison to the individual constituents including the effects of anisotropy of HMX crystals and the interactions of fine crystallites with binder material.
Binders such as Estane, Teflon, Kel F and HTPB are typically used in heterogeneous explosives to bond polycrystalline constituents together as an energetic composite. Combined theoretical and experimental studies are underway to unravel the mechanical response of these materials when subjected to isentropic compression loading. Key to this effort is the determination of appropriate constitutive and EOS property data at extremely high stress-strain states as required for detailed mesoscale modeling. The Sandia Z accelerator and associated diagnostics provides new insights into mechanical response of these nonreactive constituents via isentropic ramp-wave compression loading. Several thicknesses of samples, varied from 0.3 to 1.2 mm, were subjected to a ramp load of {approx}42 Kbar over 500 ns duration using the Sandia Z-machine. Profiles of transmitted ramp waves were measured at window interfaces using conventional VISAR. Shock physics analysis is then used to determine the nonlinear material response of the binder materials. In this presentation we discuss experimental and modeling details of the ramp wave loading ICE experiments designed specifically for binder materials.
The shock response of heterogeneous materials involves highly fluctuating states and localization effects that are produced by mesostructure. Prior studies have examined this shock behavior in randomized inert and reactive media. In this work, we investigate the shock behavior in a porous lattice consisting of hexagonally packed layers of 500 {micro}m tin spheres impacted at 0.5 km/s. This ordered geometry provides a well-defined configuration to validate mesoscale material modeling based on three-dimensional CTH calculations. Detailed wave fields are experimentally probed using a line-imaging interferometer and transmitted particle velocities are compared to numerical mesoscale calculations. Multiple shock fronts traverse the porous layers whereby particle-to-particle interactions cause stress bridging effects and the evolution of organized wave structures.
The specific problem to be addressed in this work is the secondary combustion that arises from shock-induced mixing in volumetric explosives. It has been recognized that the effects of combustion due to secondary mixing can greatly alter the expansion of gases and dispersal of high-energy explosive. Furthermore, this enhanced effect may be a tailored feature for the new energetic material systems. One approach for studying this problem is based on the use of Large Eddy Simulation (LES) techniques. In this approach, the large turbulent length scales of motion are simulated directly while the small scales of turbulent motion are explicitly treated using a subgrid scale (SGS) model. The focus of this effort is to develop a SGS model for combustion that is applicable to shock-induced combustion events using probability density function (PDF) approaches. A simplified presumed PDF combustion model is formulated and implemented in the CTH shock physics code. Two classes of problems are studied using this model. The first is an isolated piece of reactive material burning with the surrounding air. The second problem is the dispersal of highly reactive material due to a shock driven explosion event. The results from these studies show the importance of incorporating a secondary combustion modeling capability and the utility of using a PDF-based description to simulate these events.
In order to provide real-time data for validation of three dimensional numerical simulations of heterogeneous materials subjected to impact loading, an optically recording velocity interferometer system (ORVIS) has been adapted to a line-imaging instrument capable of generating precise mesoscopic scale measurements of spatially resolved velocity variations during dynamic deformation. Combining independently variable target magnification and interferometer fringe spacing, this instrument can probe a velocity field along line segments up to 15 mm in length. In high magnification operation, spatial resolution better than 10 {micro}m can be achieved. For events appropriate to short recording times, streak camera recording can provide temporal resolution better than 0.2 ns. A robust method for extracting spatially resolved velocity-time profiles from streak camera image data has been developed and incorporated into a computer program that utilizes a standard VISAR analysis platform. The use of line-imaging ORVIS to obtain measurements of the mesoscopic scale dynamic response of shocked samples has been demonstrated on several different classes of heterogeneous materials. Studies have focused on pressed, granular sugar as a simulant material for the widely used explosive HMX. For low-density (65% theoretical maximum density) pressings of sugar, material response has been investigated as a function of both impact velocity and changes in particle size distribution. The experimental results provide a consistent picture of the dispersive nature of the wave transmitted through these samples and reveal both transverse and longitudinal wave structures on mesoscopic scales. This observed behavior is consistent with the highly structured mesoscopic response predicted by 3-D simulations. Preliminary line-imaging ORVIS measurements on HMX as well as other heterogeneous materials such as foam and glass-reinforced polyester are also discussed.
The Baer-Nunziato multiphase reactive theory for a granulated bed of energetic material is extended to allow for dynamic damage processes, that generate new surfaces as well as porosity. The Second Law of Thermodynamics is employed to constrain the constitutive forms of the mass, momentum, and energy exchange functions as well as those for the mechanical damage model ensuring that the models will be dissipative. The focus here is on the constitutive forms of the exchange functions. The mechanical constitutive modeling is discussed in a companion paper. The mechanical damage model provides dynamic surface area and porosity information needed by the exchange functions to compute combustion rates and interphase momentum and energy exchange rates. The models are implemented in the CTH shock physics code and used to simulate delayed detonations due to impacts in a bed of granulated energetic material and an undamaged cylindrical sample.
The Greater Confinement Disposal (GCD) facility was established in Area 5 at the Nevada Test Site for containment of waste inappropriate for shallow land burial. Some transuranic (TRU) waste has been disposed of at the GCD facility, and compliance of this disposal system with EPA regulation 40 CFR 191 must be evaluated. We have adopted an iterative approach in which performance assessment results guide site data collection, which in turn influences the parameters and models used in performance assessment. The first iteration was based upon readily available data, and indicated that the GCD facility would likely comply with 40 CFR 191 and that the downward flux of water through the vadose zone (recharge) had a major influence on the results. Very large recharge rates, such as might occur under a cooler, wetter climate, could result in noncompliance. A project was initiated to study recharge in Area 5 by use of three environmental tracers. The recharge rate is so small that the nearest groundwater aquifer will not be contaminated in less than 10,000 years. Thus upward liquid diffusion of radionuclides remained as the sole release pathway. This second assessment iteration refined the upward pathway models and updated the parameter distributions based upon new site information. A new plant uptake model was introduced to the upward diffusion pathway; adsorption and erosion were also incorporated into the model. Several modifications were also made to the gas phase radon transport model. Plutonium solubility and sorption coefficient distributions were changed based upon new information, and on-site measurements were used to update the moisture content distributions. The results of the assessment using these models indicate that the GCD facility is likely to comply with all sections of 40 CFR 191 under undisturbed conditions.
The Greater Confinement Disposal (GCD) facility was established by the Nevada office of the Department of Energy (DOE) in Area 5 at the Nevada Test Site for containment of waste inappropriate for shallow land burial. Some transuranic (TRU) waste has been disposed of at the GCD facility, and compliance of this disposal system with Environmental Protection Agency (EPA) regulations 40 CFR 191 must be evaluated by performance assessment calculations. We have adopted an iterative approach where performance assessment results guide site data collection which in turn influences the parameters and models used in performance assessment. The first iteration was based upon readily available data. The first iteration indicated that the GCD facility would likely comply with 40 CFR 191 and that the downward recharge rate had a major influence on the results. As a result, a site characterization project was initiated to study recharge in Area 5 by use of three environmental tracers. This study resulted in the conclusion that recharge was extremely small, if not negligible. Thus, downward advection to the water table is no longer considered a viable release pathway, leaving upward liquid diffusion as the sole release pathway. This second performance assessment iteration refined the upward pathway models and parameters. The results of the performance assessment using these models still indicate that the GCD site is likely to comply with all sections of 40 CFR 191.
In this work, compaction-induced combustion in packed beds of nitrocellulose-based ball propellants is modeled using a multiphase mixture description. This model is applied to conditions simulating low-velocity impact experiments of Sandusky, et. al. (NSWC). A two- stage combustion model is used whereby compressive reaction begins at the compaction front. Subsequent energy release is delayed following an induction rate law based on time-to-reaction experimental data. Given conditions of sufficient energy release and heat transfer, grain burning is initiated when granular surface temperatures exceed decomposition conditions. Numerical solutions of the one-dimensional multiphase conservation equations are obtained using an adaptive finite element method and calculations are compared to experiments investigating various impact loading conditions on the ball propellants TS3659 and WC140. 11 refs., 12 figs., 1 tab.
