Publications

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.

Abstract not provided.

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.

Abstract not provided.

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.