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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessing the vulnerability of solar inverters to EMPs: Port testing, PCI modeling, and protection strategies

Renewable energy sources are becoming an ever-larger contributor to the power grid. These renewable energy sources depend upon the power electronic devices, specifically inverters, being essential for connecting Photovoltaic (PV) generation to the grid. However, the Electromagnetic Pulses (EMPs) caused by the high-altitude nuclear explosions can generate fast broad-band pulses with nanosecond rise time, potentially causing damage or destruction to electronic components. To assess the vulnerability of PV inverters to high-altitude EMPs, the port testing and Pulsed Current Injection (PCI) modeling schemes are proposed based on the port impedance analysis. Wide-band frequency measurements are achieved by fusing impedance results from three vector network analyzers. Then, a PCI model is used to simulate the induced response to EMP, with two typical immunity levels of EC5 and EC8 tested. Here, the experiment successfully excites the induced voltage and current under EMP, where the voltage and current can reach 1500V/40A and 8000V/150Aunder EC5 andEC8, respectively. The port vulnerability analysis results demonstrate that only some ports can survive under EC5. To defend against the impact of EMP, three protection strategies are discussed.

42 ENGINEERING↗
DOE OSTI2020

Simulations of Sweeping Wave Propagation in a Boron Carbide Plate

In February 2020, experiments were conducted at Los Alamos National Laboratory measuring the propagation of a sweeping wave through a plate of boron carbide which was induced by detonating an adjacent charge of high explosive. Six such shots were red and in each case progress of the wave along the exposed surface of the plate was tracked using photon doppler velocimetry (PDV). In this report we discuss calculations performed using the Lagrangian hydrocode FLAG to model data taken from these experiments. We discuss the models chosen for the different experimental components, the choice of mesh (focusing on, both, two- and three-dimensional geometries), and compare the results from these calculations to experimentally-extracted PDV traces. Notably, we find that the PDV traces for this experiment are most faithfully reproduced using the full three-dimensional geometry, even when the mesh resolution for the three-dimensional calculations was significantly lower than the resolution for the corresponding two-dimensional calculations. We attribute this to the geometry of the experimental setup. Moreover, we find only negligible differences for the PDV traces obtained using the sesame table equation of state for boron carbide and the Gruneisen equation of state.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetic Energy Release, Plasma Dynamics, and Particle Acceleration in Relativistic Turbulent Magnetic Reconnection

In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be the primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue in studies of relativistic reconnection is the consequence of three-dimensional dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of three-dimensional, fully kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron–electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long-wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. Additionally, we find that plasma dynamics in RTMR is vastly different from its 2D counterpart in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted by the secondary kink instability. This turbulent evolution leads to superdiffusive behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and the timescale for energy release can be very fast and do not depend strongly on the turbulence amplitude. Additionally, the main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the nonideal electric field acceleration plays a subdominant role. We also discuss possible observational implications of three-dimensional RTMR in high-energy astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Rock Valley Dense Gravity Acquisition

Characterizing the shallow structure of the Rock Valley region of the Nevada National Security Site is a critical component of the Rock Valley Direct Comparison project. Geophysical data of the region is needed for operational decisions, to constrain geologic models used for simulation, and to facilitate the analysis of future explosive source data. Local measurements of gravity are a key piece of geophysical information that helps to resolve the underlying geologic composition, fault structure, and density characteristics, yet, in the Rock Valley region these measurements are sparse on the scale of the testbed. In this report, we present the details of a recent gravity data acquisition survey designed to collect a dense dataset in the region of interest that complements the existing gravity work but greatly enhances our resolution. This dataset will be integrated with a complementary Los Alamos National Laboratory gravity collection and combined with the existing seismic data in a joint inversion. These measurements were conducted over two weeks with a portable gravimeter and high-resolution GPS and include repeat measurements at a USGS base station as well as reoccupation of gravity sites in the regional dataset. This collection of over 100 new dense gravity measurements will facilitate refinement of the existing Geologic Framework Model and directly complement newly acquired dense seismic data, ultimately improving the project’s ability to investigate the direct comparison of shallow earthquake and explosive sources.

58 GEOSCIENCES↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Singularity-EOS XCAP Report

Solving the Euler equations is a fundamental component of simulating many physical phenomena ranging from high explosives to astrophysics. An equation of state (EOS) is a required piece that relates any two thermodynamic quantities to all other thermodynamic values. The presence of multiple materials within a control volume further complicates the solution requiring additional equations to describe the interaction of materials at a sub-grid level. Equations of state themselves can also come in many forms ranging from simple algebraic relations to more complicated differential equation models that describe material interactions over a broad range of physical conditions. In the latter case, the EOS is often pre-computed at a given set of grid points and provided in a tabular form where additional properties can be derived from the interpolation functions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

Asymmetries in Core-Collapse Supernovae from Maps of Radioactiver 44Ti in Cassiopeia A

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surroundingmedium1, directly probes the explosion asymmetries. Cassiopeia A is a young2, nearby3, core-collapse4 remnant from which 44Ti emission has previously been detected5-8 but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission9, from optical light echoes10, and from jet-like features seen in the X-ray11 and optical12 ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae.

stellarcore-collapse explosions↗
DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutrino heating in 1D, 2D, and 3D core-collapse supernovae: characterizing the explosion of high-compactness stars

Massive stars can end their lives with a successful supernova explosion (leaving behind a neutron star or, more rarely, a black hole), or a failed explosion that leaves behind a black hole. The density structure of the pre-collapse progenitor star already encodes much of the information regarding the outcome and properties of the explosion. However, the complexity of the collapse and subsequent shock expansion phases prevents drawing a straightforward connection between the pre-collapse and post-explosion properties. In order to derive such a connection several explodability studies have been performed in recent years. However, different studies can predict different explosion outcomes. In this article, we show how compactness, which is related to the average density of the star’s core, has an important role in determining the efficiency of neutrino heating, and therefore the outcome of the explosion. Commonly, high-compactness progenitors are assumed to yield failed explosions, due to their large mass accretion rates, preventing the shock from expanding. We show by analysing ~150 2D flash and F ornax simulations and 20 3D F ornax simulations that this is not the case. Instead, due to the rapid increase of neutrino heating with compactness, high-compactness progenitors lead to successful shock revival. We also show that 1D+ simulations that include v-driven convection using a mixing-length theory approach correctly reproduce this trend. Finally, we compare 1D+ models, which we show can reproduce some aspects of multi-D simulations with reasonable accuracy, with other widely used 1D models in the literature.

neutrinos↗
NASA NTRSDate not supplied

Digital Twins and Living Models at NASA

The idea of a “digital twin” was born at NASA in the 1960s as a “living model” of the Apollo mission. In response to Apollo 13’s oxygen tank explosion and subsequent damage to the main engine, NASA employed multiple simulators to evaluate the failure and extended a physical model of the vehicle to include digital components. This “digital twin” was the first of its kind, allowing for a continuous ingestion of data to model the events leading to up to the accident for forensic analysis and exploration of next steps. Fast forward half a century and NASA, along with others in the aerospace community, continues to develop and utilize high-fidelity digital models of physical systems and components as well as the extreme environments in which they operate. NASA aims to travel further and stay longer in space as we realize the Artemis program, taking us from the moon to Mars by establishing a sustainable presence on the Moon to prepare for missions to Mars. We will no longer be able to rely on constant connectivity with an asset nor be in-the-loop for on-demand human intervention in the event of an anomaly. Further, the importance of digital twins is increasing as we seek alternatives for certification of structures so large that they cannot be fully evaluated in existing test facilities and autonomous systems that are not deterministic. The idea of a digital twin is not a new one. What is new is the scale, ordinality, and non-deterministic nature of the models that are critical to achieving NASA’s goals. Their number and autonomy from each other as well as the reference system is suggestive of a changing ecosystem, returning us the idea of a living model.

Digital Twin↗
DOE OSTI2025

The Hunga Volcanic Eruption Atmospheric Impacts Report

On 15 January 2022 a highly explosive eruption of the Hunga volcano occurred in the Kingdom of Tonga in the South Pacific Ocean (175°24’ W, 20°33’ S). The Volcanic Explosivity Index (VEI) 6 eruption originated from a shallow submarine vent, making it distinct from large subaerial eruptions of recent decades (e.g., 1982 El Chichón, 1991 Mt. Pinatubo). In particular, seawater enhanced explosivity and dampened sulfur dioxide (SO 2 ) emissions. The eruption was the culmination of ~1 month of precursory activity; however, the timing and size of the eruption were unexpected, partly due to the challenges of monitoring submarine volcanoes. The stratospheric hydration caused by the eruption was unprecedented in magnitude, altitude, and duration in the satellite record. This Executive Summary reflects the current assessment of the Hunga eruption and its impact on the climate system. We report key observations of the eruption and its aftermath, as well as simulations of its impact by global chemistry-climate models. The Hunga eruption had an unprecedented impact on the stratosphere and mesosphere due to the plume height and large water content, which increased the global stratospheric water vapour burden by 10%. Most of this water has remained in the atmosphere into 2025. However, Hunga’s net impact on surface climate was small compared to that of earlier large-magnitude volcanic eruptions, due to limited sulfate aerosol loading in the stratosphere and the high altitude of the water vapour injection.

58 GEOSCIENCES↗
DOE OSTI2023

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Towards Resilient Near Real-Time Analysis Workflows in Fusion Energy Science

Nuclear fusion holds the promise of an endless source of energy. Several research experiments across the world and joint modeling and simulation efforts between the nuclear physics and high performance computing communities are actively preparing the operation of the International Thermonuclear Experimental Reactor (ITER). Both experimental reactors and their simulated counterparts generate data that must be analyzed quickly and in a resilient way to support decision making for the configuration of subsequent runs or prevent a catastrophic failure. However, the cost if the traditional techniques used to improve the resilience of analysis workflows, i.e., replicating datasets and computational tasks, becomes prohibitive with explosion of the volume of data produced by modern instruments and simulations. Therefore, we advocate in this paper for an alternate approach based on data reduction and data streaming. The rationale is that by allowing for a reasonable, controlled, and guaranteed loss of accuracy it becomes possible to transfer smaller amounts of data, shorten the execution time of analysis workflows, and lower the cost of replication to increase resilience. We develop our research and development roadmap towards resilient near real-time analysis workflows in fusion energy science and present early results showing that data streaming and data reduction is a promising way to speed up the execution and improve the resilience of analysis workflows.

Suter, Fred↗
NASA NTRS2010

Reconnection Onset in the Breakout Model for CME Initiation

Fast coronal mass ejections (CMEs) are the most massive explosions in the heliosphere, and the primary drivers of geoeffective space weather. Although it is generally agreed that magnetic reconnection is the key to fast CME initiation, different models incorporate reconnection in different ways. One promising model --- the breakout scenario --- involves reconnection in two distinct yet interconnected locations: breakout reconnection ahead of the CME, and flare reconnection behind it. We will discuss what we have learned about the early evolution of breakout and flare reconnection from recent high-resolution 2.5D adaptively refined MHD simulations of CME initiation, including the evolving properties of the breakout and flare current sheets, the conditions that trigger reconnection onset in each sheet, the ensuing positive feedback between breakout and flare reconnections, and implications for electron acceleration in flares.

Karpen, Judy T.↗
DOE OSTI2021

Empirical Acoustic Source Model for Chemical Explosions in Air

Chemical explosions generate pressure disturbances in air that radiate as nonlinear shock waves near the source and transition into acoustic waves with distance. Because low-frequency acoustic waves generally travel large distances without significant loss of energy, they are often used for explosion monitoring and yield estimation. However, quantitative relationships between acoustic energy and explosion yields are required for accurate yield estimation. Here, we develop an empirical acoustic source model for chemical explosions from experimental data. The empirical model returns the acoustic pressure waveform for the detonation of 1 kg of trinitrotoluene, which is conventionally used to represent the explosive release of 4.184 MJ of explosion energy. The full-waveform model can be used to predict acoustic signals for an arbitrary yield of a high-explosive detonation based on the standard scaling law and to estimate acoustic energies in a specific frequency range. We evaluate the accuracy of the acoustic source model independently by estimating the yield of other explosive events that are not included in the model development. Finally, statistical characteristics of the model and their implications for the uncertainty quantification of estimated yields are discussed.

58 GEOSCIENCES↗