Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high explosives, simulations, modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

An integrated EOS, pore‐crush, strength and damage model framework for near‐field ground‐shock

Abstract An integrated Equation of State (EOS) and strength/pore‐crush/damage model framework is provided for modeling near to source (near‐field) ground‐shock response, where large deformations and pressures necessitate coupling EOS with pressure‐dependent plastic yield and damage. Nonlinear pressure‐dependence of strength up to high‐pressures is combined with a Modified Cam‐Clay‐like cap‐plasticity model in a way to allow degradation of strength from pore‐crush damage, what we call the “Yp‐Cap” model. Nonlinear hardening under compaction allows modeling the crush‐out of pores in combination with a fully saturated EOS, that is, for modeling partially saturated ground‐shock response, where air‐filled voids crush. Attention is given to algorithmic clarity and efficiency of the provided model, and the model is employed in example numerical simulations, including finite element simulations of underground explosions to exemplify its robustness and utility.

54 ENVIRONMENTAL SCIENCES↗

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗

Optical and Near-Infrared Observations of SN 2013DX Associated with GRB 130702A

We present optical and near-infrared (NIR) light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be E(sub gamma, iso) = 6.4(+1.3/-1.0) x 10(exp 50) erg (1 keV to 10 MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed g'r'i'z' light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves approx. 20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined SNe Ic, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of approx. 21,000 km/s. We construct a quasi-bolometric (g'r'z'yJ) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond (Delta)t > 40 days. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a Ni-56 mass of M(sub Ni) = 0.37+/- 0.01 Stellar Mass, an ejecta mass of M(sub ej) = 3.1+/- 0.1 Stellar Mass, and a kinetic energy of E(sub K) = (8.2+/- 0.43) x 10(exp 51) erg (statistical uncertainties only), consistent with previous GRB-associated supernovae. When considering the ensemble population of GRB-associated supernovae, we find no correlation between the mass of synthesized Ni-56 and high-energy properties, despite clear predictions from numerical simulations that M(sub Ni) should correlate with the degree of asymmetry. On the other hand, M(sub Ni) clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

gamma-ray burst: individual (GRB 130702A) – supe↗

Experimental investigation of nozzle/plume aerodynamics at hypersonic speeds

Much of the work involved the Ames 16-Inch Shock Tunnel facility. The facility was reactivated and upgraded, a data acquisition system was configured and upgraded several times, several facility calibrations were performed and test entries with a wedge model with hydrogen injection and a full scramjet combustor model, with hydrogen injection, were performed. Extensive CFD modeling of the flow in the facility was done. This includes modeling of the unsteady flow in the driver and driven tubes and steady flow modeling of the nozzle flow. Other modeling efforts include simulations of non-equilibrium flows and turbulence, plasmas, light gas guns and the use of non-ideal gas equations of state. New experimental techniques to improve the performance of gas guns, shock tubes and tunnels and scramjet combustors were conceived and studied computationally. Ways to improve scramjet engine performance using steady and pulsed detonation waves were also studied computationally. A number of studies were performed on the operation of the ram accelerator, including investigations of in-tube gasdynamic heating and the use of high explosives to raise the velocity capability of the device.

Bogdanoff, David W.↗

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗

Simulation Assisted Risk Assessment: Blast Overpressure Modeling

A probabilistic risk assessment (PRA) approach has been developed and applied to the risk analysis of capsule abort during ascent. The PRA is used to assist in the identification of modeling and simulation applications that can significantly impact the understanding of crew risk during this potentially dangerous maneuver. The PRA approach is also being used to identify the appropriate level of fidelity for the modeling of those critical failure modes. The Apollo launch escape system (LES) was chosen as a test problem for application of this approach. Failure modes that have been modeled and/or simulated to date include explosive overpressure-based failure, explosive fragment-based failure, land landing failures (range limits exceeded either near launch or Mode III trajectories ending on the African continent), capsule-booster re-contact during separation, and failure due to plume-induced instability. These failure modes have been investigated using analysis tools in a variety of technical disciplines at various levels of fidelity. The current paper focuses on the development and application of a blast overpressure model for the prediction of structural failure due to overpressure, including the application of high-fidelity analysis to predict near-field and headwinds effects.

Lawrence, Scott L.↗

Solid Rocket Launch Vehicle Explosion Environments

Empirical explosion data from full scale solid rocket launch vehicle accidents and tests were collected from all available literature from the 1950s to the present. In general data included peak blast overpressure, blast impulse, fragment size, fragment speed, and fragment dispersion. Most propellants were 1.1 explosives but a few were 1.3. Oftentimes the data from a single accident was disjointed and/or missing key aspects. Despite this fact, once the data as a whole was digitized, categorized, and plotted clear trends appeared. Particular emphasis was placed on tests or accidents that would be applicable to scenarios from which a crew might need to escape. Therefore, such tests where a large quantity of high explosive was used to initiate the solid rocket explosion were differentiated. Also, high speed ground impacts or tests used to simulate such were also culled. It was found that the explosions from all accidents and applicable tests could be described using only the pressurized gas energy stored in the chamber at the time of failure. Additionally, fragmentation trends were produced. Only one accident mentioned the elusive "small" propellant fragments, but upon further analysis it was found that these were most likely produced as secondary fragments when larger primary fragments impacted the ground. Finally, a brief discussion of how this data is used in a new launch vehicle explosion model for improving crew/payload survival is presented.

Richardson, E. H.↗

Mixed-precision numerics in scientific applications: survey and perspectives

The explosive demand for artificial intelligence (AI) workloads has led to a significant increase in silicon area dedicated to lower-precision computations on recent high-performance computing hardware designs. However, mixed-precision capabilities, which can achieve performance improvements of up to 8x compared to double-precision in extreme compute-intensive workloads, remain largely untapped in most scientific applications. A growing number of efforts have shown that mixed-precision algorithmic innovations can deliver superior performance without sacrificing accuracy. These developments should prompt computational scientists to seriously consider whether their scientific modeling and simulation applications could benefit from the acceleration offered by new hardware and mixed-precision algorithms. In this survey, we (1) review progress across diverse scientific domains—fluid dynamics, weather and climate, quantum chemistry, and computational genomics—that have begun adopting mixed-precision strategies; (2) examine state-of-the-art algorithmic techniques such as iterative refinement, splitting and emulation schemes, and adaptive precision solvers; (3) assess their implications for accuracy, performance, and resource utilization; and (4) survey the emerging software ecosystem that enables mixed-precision methods at scale. We conclude with perspectives and recommendations on cross-cutting opportunities, domain-specific challenges, and the role of co-design between application scientists, numerical analysts, and computer scientists. Collectively, this survey underscores that mixed-precision numerics can reshape computational science by aligning algorithms with the evolving landscape of hardware capabilities.

Graphics processing units↗

Detectability of collective neutrino oscillation signatures in the supernova explosion of a 8.8 M ⊙ star

In order to investigate the impact of collective neutrino oscillations (CNOs) on the neutrino signal from a nearby supernova, we perform three-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 M ⊙ progenitor star. In this work, we find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\bar{ν}_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to ~ 10 kpc. On the other hand, for the normal mass hierarchy, the spectrum of ν e becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to ~ 10 kpc . More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\bar{ν}_e$ in HK becomes softer due to CNO, the spectrum of $ν_e$ in DUNE becomes harder, and vice versa. These synergistic observations in $\bar{ν}_e$ and $ν_e$, by HK and DUNE, respectively, will be an intriguing opportunity to test the occurrence of CNO.

79 ASTRONOMY AND ASTROPHYSICS↗

Bolt cutter functional evaluation

The Aerospace Corporation has been implementing finite difference and finite element codes for the analysis of a variety of explosive ordnance devices. Both MESA-2D and DYNA3D have been used to evaluate the role of several design parameters on the performance of a satellite separation system bolt cutter. Due to a lack of high strain rate response data for the materials involved, the properties for the bolt cutter and the bolt were selected to achieve agreement between computer simulation and observed characteristics of the recovered test hardware. The calculations provided insight into design parameters such as the cutter blade kinetic energy, the preload on the bolt, the relative position of the anvil, and the anvil shape. Modeling of the cutting process clarifies metallographic observation of both cut and uncut bolts obtained from several tests. Understanding the physical processes involved in bolt cutter operation may suggest certain design modifications that could improve performance margin without increasing environmental shock response levels.

S. Goldstein↗

Risk-Based Explosion Assessment Internal Explosion of Exploration Upper Stage in Vehicle Assembly Building

This paper presents risk-based explosion analyses used to determine if an explosion of NASA's Space Launch System (SLS) Exploration Upper Stage (EUS) during fueling in Vehicle Assembly Building (VAB) High Bay 4 can hazard the Space Launch System (SLS) being processed in Bay 3. This paper focuses on blast fragmentation as it was the controlling hazard. An estimate of the amount of Net Explosion Weight (NEW) was used to develop a fragment list (mass, shape, takeoff velocity/angle) given an explosion occurs at ground level in Bay 4. Intervening superstructure that exists between Bay 4 and Bay 3 (that could potentially block EUS fragments) was modeled based on detailed visual images (VISSIM) and AUTOCAD/structural drawings. The fragment list was input to a 3D trajectory application that performs multiple random explosion simulations; each simulation varying fragment ballistic coefficient and takeoff angle/velocity. Individual fragments are tracked until stopped by the intervening superstructure or reach Bay 3, potentially impacting the SLS. The probability of impacting the SLS is then determined from the random simulations; the results are used to help decision-makers determine the efficacy of co-processing within the VAB and decide what mitigations may be appropriate.

Explosion Assessment↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Learning Physically Interpretable Atmospheric Models From Data With WSINDy

The multiscale and turbulent nature of Earth's atmosphere has historically rendered accurate weather modeling a hard problem. Recently, there has been an explosion of interest surrounding data-driven approaches to weather modeling, which in many cases show improved forecasting accuracy and computational efficiency when compared to traditional methods. However, many of the current data-driven approaches employ highly parameterized neural networks, often resulting in uninterpretable models and limited gains in scientific understanding. In this work, we address the interpretability problem by explicitly discovering partial differential equations governing atmospheric phenomena, identifying symbolic mathematical models with direct physical interpretations. The purpose of this paper is to demonstrate that, in particular, the weak-form sparse identification of nonlinear dynamics (WSINDy) algorithm can learn effective atmospheric models from both simulated and assimilated data. Our approach adapts the standard WSINDy algorithm to work with high-dimensional fluid data of arbitrary spatial dimension.

58 GEOSCIENCES↗

The Coronae of Moderate-Mass Giants in the Hertzsprung Gap and the Clump

We have used the Roentgensatellit (ROSAT), the Extreme Ultraviolet Explorer (EUVE), and the Hubble Space Telescope (HST) to measure X-ray and ultraviolet emissions of moderate-mass (APPROX. 2-3 solar mass) giants in the Hertzsprung gap (spectral types early F to mid-G) and the post-helium flash "clump" (approx. G8-K0). Our motivation was to document the evolution of hot coronae (T greater than 10(exp 6)K) along the post-main-sequence trajectories traveled by such stars in order to gain insight concerning the "X-ray deficiency" of the F-GO giants and the strong braking of stellar rotation at the red edge of the Hertzsprung gap. With few exceptions, Hertzsprung gap and clump giants observed by ROSAT PSPC show hot (T approx. 10(exp 7)K) coronal energy distributions, regardless of any X-ray deficiency, EUVE spectra of gap star 31 Com (G0 111) indicate a broad coronal emission measure hump at approx. 10(exp 7.2)K, while the active clump giant beta Ceti (K0 III) displays a sharp peak at approx. 10(exp 6.8)K, as seen previously in the mixed clump/gap binary Capella (alpha Aur: G8 III + G0 III). The gap giants upsilon Peg (F8 III) and 24 UMa (G4 III) have EUV emissions of intermediate temperature (approx. 10(exp 7.0)K). The stars 31 Com, psi(sup 3) Psc (G0 III), and beta Cet exhibit redshifted transition zone (TZ: approx. 10(exp 5)K) lines in HST GHRS spectra, as reported earlier in Procyon (alpha CMi: F5 IV-V) and Capella G0. Such redshifts on the Sun are thought to signify flows in magnetic loops. beta Cas (F2 III)-a rare soft coronal source among the gap stars-displays blueshifts of C iv and 0 iv], although emissions at cooler and hotter temperatures are near the photospheric velocity. The remarkably broad line profiles of the fastest rotating gap giants suggest that the 10(exp 5) K "subcoronal " emission zones extend to h approx. R(sub *) above the photosphere, about 50 scale heights. In contrast to the TZ line redshifts, the upper chromospheric emissions (e.g., Mg II and Si III) of 31 Com and upsilon(sup 3) Psc have blueshifted cores. Blue-asymmetric peaks in the solar Mg a lines are thought to indicate dynamical heating in the chromosphere. Observations of the H(sub I) Ly(alpha) feature of 31 Com taken 9 months apart reveal striking profile changes, reminiscent of those noted previously in the Ly(alpha) blue peak of the Capella G0 star. We used the far-ultraviolet diagnostics, in combination with ROSAT X-ray photometry and EUVE high-excitation line strengths, to constrain physical models of the stellar outer atmospheres. Quasi-static magnetic loops can simulate the empirical coronal emission measures of the giant stars, but the inferred pressures for sensible loop lengths conflict with direct measurements of subcoronal densities. Furthermore, the high rate of emission at approx. 10(exp 5) K cannot be explained by thermal conduction down the legs of hot quasi-static loops. On the other hand, the possible existence of elongated (l - R(sub *) emission structures on the gap giants leads to a speculative scenario to explain the X-ray deficiency. It is based on the increased importance of the dynamical filling phase ("explosive evaporation") of the loop life cycle; conductive cooling, yielding TZ emissions at the footpoints, when the heating is interrupted; and the possibility for transitions between " hot " and " cool " energy balance solutions owing to dynamical suspension and centrifugal trapping of the cooling gas. The long loops might represent a vestigial global " magnetosphere " inherited from the main-sequence phase, which ultimately is disrupted near approx. G0 by the deepening convective envelope and growth of a more solar-like dynamo. Coronal emissions might be boosted temporarily as the X-ray deficiency is removed but soon would be quenched by wind braking previously inhibited by the magnetospheric "dead zone."

Ayres, Thomas R.↗

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

HEAO 3 limits on the Ti-44 yield in Galactic supernovae

Data fron the high-resolution gamma-ray spectroscopy experiment on HEAO 3 have been searched for line emission from the decay of Ti-44 created in recent, as yet unobserved, Galactic supernova explosions, where the ages and locations are unknown. Because the 78 yr mean life of Ti-44 is comparable to the average time between Galactic supernovae, the gamma-ray line emission from its decay should appear as Galactic point sources. No evidence was found for such emission from a point source anywhere in the Galactic plane, with a 1-sigma limit of 8.3 x 10 exp -5 photons/sq cm per sec. Detailed models were developed to simulate the Galactic gamma-ray emission from the decay of Ti-44 produced in both type I and type II supernovae. These models were used with the measured gamma-ray line limits to constrain the supernova yields and recurrence periods.

Mahoney, W. A.↗

FY2021 Improvements to the New CTH Code Verification & Validation Test Suite

Over the past few years, the CTH multiphysics hydrocode has overhauled its software quality and testing processes, implementing current best practices in software quality and building a robust V&V test suite comprised of traditional hydrocode verification problems, including ASC Tri-Lab Test Suite and Enhanced Tri-Lab Test Suite problems, as well as validation problems for some of CTH’s most frequently used equations of state, materials models, and other key capabilities. Substantial progress towards building this new test suite was made in FY19 and FY20. In FY21, the test suite has been expanded to include verification and validation tests of the Steinberg-Guinan-Lund (ST) viscoplastic model and the Johnson Cook (JFRAC) fracture model. Additionally, two new verification tests were added, covering hydrodynamics and high explosive (HE) modeling capabilities: the Kidder Gaussian density problem and the Escape of HE Products (EHEP) problem from the Tri-Lab Test Suite. This report discusses each of these test problems in detail. Verification test results are compared to analytic solutions. Validation test results are compared to experimental data. Wherever possible, convergence or mesh refinement studies are included. Additionally, while implementing the Kidder verification problem, a bug was identified that affects the use of tables to initialize pressure or density in 1D or 2D calculations. A brief discussion of the bug and its fix is included. CTH demonstrates good performance overall on the new test suite problems. Simulation results showed good agreement with analytic solutions for the Kidder problem, with convergence rates ranging between 1.8 and sub-linear, and relatively good agreement for the EHEP problem, though convergence rates for pressure and density were nearly 0. The ST and JFRAC strain rate loading verification tests show good agreement with analytic solutions. Likewise, CTH simulation results show good agreement with experimental validation data, including Taylor rod impact testing, for the materials tested. Future V&V work will focus on adding 2D and 3D versions of existing verification tests as well as adding validation tests of other frequently used capabilities such as other fracture models.

42 ENGINEERING↗

Experimental observations of shock-wave-induced bubble collapse and hot-spot formation in nitromethane liquid explosive

We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗