Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Detonators”

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 523 records · Page 29

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

Measurements of Argon-39 from locations near historic underground nuclear explosions

Measurement of radioactive gas seepage from an underground nuclear explosion is one of the primary methods to confirm whether an event was nuclear in nature. Radioactive noble gas indicators that are commonly targeted by such measurements (e.g. 133Xe, 37Ar) have half-lives of 35 days or less. Argon-39, an activation product similar to 37Ar, is produced by the interaction between neutrons and potassium in the surrounding geology and has a half-life of 269 years. Measurements taken at three sites near three historic underground nuclear test locations at the Nevada National Security Site have all shown highly elevated levels of 39Ar in subsurface air decades after the test events. Elevated levels of 39Ar (30-50 times background) were also detected in atmospheric air collected roughly 30 cm above ground level near two of these sites, and outside the entrance of the one tunnel site. These measurements demonstrate that 39Ar has the potential to be a long-term signature of an underground nuclear event which can be reliably detected at the surface or in the shallow subsurface. This radionuclide detection of an underground nuclear event decades after the event takes place is in contrast to the commonly held assumption that detecting underground nuclear events via radionuclides at the surface needs to be done in a matter of months. Depending upon what further studies show about the robustness of this signature in a variety of geological settings, it may in fact be easy to detect underground nuclear events at the surface for a very long time post-detonation.

Johnson, Christine M.↗

Analysis of 127 Xe tracer measurements using a net counts method

A suite of measurement systems were deployed as part of the Physical Experiment 1 series of experiments, which involved detonating chemical explosives along with radionuclide tracers in an underground cavity, at the Nevada National Security Site (NNSS) in the United States. One of the radionuclide tracers, 127 Xe was released from the containment following the explosion and detected on a SAUNA Q B sampler situated approximately 3.5 km away. The system uses a beta-gamma coincidence detector system to measure fission product radioisotopes of xenon relevant to nuclear explosion monitoring. In this work we use the coincidence measurement data to analyse and interpret the results from the SAUNA Q B system, to calculate the measured 127 Xe activity concentration(s).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Noninvasive acoustic time-of-flight measurements in heated, hermetically-sealed high explosives using a convolutional neural network

In this work, we present a data-driven technique for measuring the time-of-flight through material sealed within a container. Time-of-flight measurement provides a noninvasive means of quantifying the sound speed profile within a material by transmitting an acoustic burst and then measuring the time required for the burst to arrive at an opposing receiver. In a hermetically-sealed cylindrical container, a portion of the acoustic energy propagates through the material as a bulk wave, while the remainder of the acoustic energy propagates around the container walls as guided waves. As a result, interference from the guided waves obscures the bulk arrival, inhibiting measurement of the sound speed. The technique uses a Convolutional Neural Network (CNN) to identify critical features in the measured waveforms and identify bulk wave arrivals. We demonstrate this time-of-flight measurement technique on high explosive-filled containers as they are heated from room temperature to detonation. This is a particularly challenging application for acoustic time-of-flight measurements as the high explosives have significant sound speed gradients as they undergo heating, and they lead to significant attenuation of the bulk wave, as opposed to the guided waves, which do not suffer significant attenuation. We characterize the performance of the CNN as a function of the high explosive temperature and as a function of the CNN hyperparameters. We then provide physical insight into the error trends.

47 OTHER INSTRUMENTATION↗

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗

Production of desensitized, ultrafine PETN powder

Pentaerythritol tetranitrate (PETN) is a widely studied high explosive (HE), most commonly used in detonator applications. In this work, we use a spray drying technique to manufacture a new “nano-powder” morphology of PETN. The spray dried PETN is several orders of magnitude smaller in particle size than traditionally prepared PETN powders. The spray dried PETN has a mostly spherical and smooth morphology, while traditional crystalline PETN typically has sharp edges and crystal facets. Small-scale sensitivity tests including drop-weight impact, friction, and electrostatic discharge (ESD) indicate the spray dried PETN is less sensitive than traditionally used forms of PETN powders. Furthermore, we also observed no changes in chemical properties (melt temperature and onset of decomposition) after spray drying, and the material remains in the tetragonal, crystalline phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Commercial aluminum powders, Part I: Particle size characterization and slow heating rate thermal analysis.

We characterized nine commercial aluminum (Al) powders using several methods to measure particle characteristics and thermal analysis, with the goal to understand how these parameters influence energy release. Although it is well-known that lot-to-lot variations in commercial nanoparticles are common, the Al powders were more heterogeneous than anticipated – both with regards to particle size distributions and impurities. Manufacturer specifications – often quoted in the literature without confirmation – were not always accurate for the specific sample lots we investigated. In several cases, different conclusions could be drawn from individual particle size techniques; a combination of multiple techniques provides a clearer picture of the powder properties. Thorough characterization of Al powders is required prior to interpretation of experimental results from a variety of applications. In particular, previous studies have shown contradictory results on the influence of Al on detonation performance, perhaps partially due to insufficient characterization of the Al powders themselves.

36 MATERIALS SCIENCE↗

Programming approaches for scalability, performance, and portability of combustion physics codes

Here, this paper presents the process, strategy, and results associated with porting a typical combustion physics flow solver to current state-of-the-art and future massively-parallel computer architectures. Major focus is placed on the distinct algorithmic structure of these types of codes and how it can be integrated with modern programming paradigms for heterogeneous platforms (i.e., distributed many-core systems with accelerators). An end-to-end case study is presented that exemplifies the process in a generic manner, which then serves as a clear guide with respect to the strategy and best practices leading to a robust and adaptable framework that performs well, is durable over time, is portable, and requires minimal human-effort. This end is accomplished beginning with the use of a mature, validated, structured, multiblock code framework optimized for application of both Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS). This code has been ported to a variety of platforms over the past decade, including most recently the Oak Ridge Leadership Computing Facility’s “Summit” Platform. The experience gained on these multiple platforms provides general insights and thus the results presented are not specific to any one code or platform other than the overarching trend toward distributed many-core systems with accelerators in order to move toward exascale performance. The resultant performance and scalability of the ported code is demonstrated on a real-world application; a state-of-the-art rotating detonation rocket engine simulation that matches the complex geometry and boundary conditions imposed as part of a companion experimental campaign.

97 MATHEMATICS AND COMPUTING↗

Optical spectroscopy and modeling of uranium gas-phase oxidation: Progress and perspectives

We report that studies related to U gas-phase oxidation through plasma- and thermo-chemistry are important for many fields, including environmental monitoring, forensic analysis, debris analysis in a weapon detonation event, and nucleation physics. Recently, significant efforts have been made to understand the chemical pathways involved in the progression from U atoms to diatoms (UO) and polyatomic molecules (UxOy), employing optical spectroscopy tools and computational modeling. In many studies, laser ablation of U or a U-containing flow reactor are used as a highly resource-efficient, repeatable, tunable, and lab-scale testbed for studying gas-phase oxidation in U plasmas. The spectroscopic analysis of high-temperature gas-phase oxidation of U is challenging due to the congested U spectra, resolution limitations of instrumentation, and the numerous chemical reaction pathways possible. This article focuses on the current understanding and challenges related to studying U plasma chemistry, specifically U gas-phase oxidation and molecular formation, via optical spectroscopy of plasmas and associated computational and spectral modeling. The physical and chemical processes involved in the evolution from U atoms to U oxide molecules to nanoparticles and agglomerates (i.e., debris) are discussed in the context of optical spectroscopic studies. The article concludes by highlighting opportunities for future research efforts based on existing knowledge published in the literature.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

U-Pu and Ba-Cs isotopic measurements on Trinitite by laser ablation sampling on the Neoma MC-ICP-MS

In this study we present the results of combined U-Pu and Ba-Cs isotope measurements obtained by laser ablation (LA) sampling of two glassy debris fragments (‘Trinitite’) from the world's first atomic bomb detonation conducted in New Mexico on July 16, 1945. Our primary goal in conducting these measurements was to understand whether examination of the U-Pu and Ba-Cs systematics by direct sampling (e.g. without any chemical separation or purification prior to isotope ratio measurement) could yield meaningful information that would differentiate the Trinitite fragments from glassy material lacking a nuclear fission signature. These measurements were conducted on a ThermoFisher Scientific Neoma multi collector – inductively coupled plasma – mass spectrometer (MC-ICP-MS), which is a relatively new MC-ICP-MS platform, so we also examine the behavior of these isotope systems in standards sampled in solution and via LA. Unsurprisingly, the measurements made on purified solutions of the U, Pu, and Ba isotopic standards produce high precision isotope ratios. Furthermore, this extends to the U-Pu measurements made by LA sampling, with the expected degradation in precision and accuracy related to matrix effects and signal intensity fluctuation. However, the Ba-Cs data acquired by LA is of low precision across all of the matrices examined and bears evidence of complex mass fractionation that will require further investigation to resolve. In total, our results indicate that the observed U-Pu isotope data are of sufficient quality to accurately constrain the U and Pu isotopic composition of glass containing sub-ppm levels of these elements which in turn could be used to differentiate glass containing anthropogenic fission products from natural glass whereas the Ba-Cs LA data cannot be used for this purpose until further methodological refinement is performed.

Ba-Cs↗

A new method for quantifying 64 Cu in nuclear debris samples

Quantifying 64 Cu in post-detonation nuclear debris samples can provide important diagnostic information regarding the structural materials used within a nuclear device. However, this task is challenging due to the weak gamma emissions associated with the decay of 64 Cu, its short half-life (12.701 h), and the presence of interfering fission product radioisotopes. Large quantities of debris sample are generally needed to accurately quantify 64 Cu, which can be problematic in sample-limited scenarios where other radiometric analyses are required. Herein, we present a new method for the separation of 64 Cu from solutions of mixed fission products and demonstrate the quantification of its activity through use of gas-flow proportional beta counting. Here, the new method was validated through a series of rigorous tests and was shown to improve the detection limit of 64 Cu by over two orders of magnitude, from 2.5 × 10 6 to 1.3 × 10 4 atoms/sample for 100 min measurements.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Rapid analysis of 237 Np and Pu isotopes in unseparated sample matrices using ICP-MS/MS

Inductively coupled plasma tandem mass spectrometry (ICP-MS/MS) is an emerging technique for measuring actinide isotopes when assessing pre- and post- detonation nuclear material. In this study, ICP-MS/MS was investigated for direct Np and Pu quantitation in unseparated, dissolved bulk soil matrices. To achieve this, purified nitric oxide (NO) was investigated for the reactivity of Th, Np, U, Pu, Am, and Cm. Purifying NO prior to the collision reaction cell (CRC) results in increased sensitivity and allows for higher gas flows to be utilized for the removal of interferences. Here, the interference from uranium hydrides was mitigated to less than 3.85 x 10 -11 . This method was demonstrated on standard reference materials which were measured for 237 Np and 238,239,240 Pu in dilute sample digestions. The 238 Pu measurement was validated by spiking into a standard reference material and was accurately measured with an excess of 85000 of 238 U.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗