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 19 records

On the Determination of Ejecta Structure and Explosion Asymmetry from the X-ray Knots of Cassiopeia A

We present a detailed analysis of Chandra X-ray spectra from individual ejecta knots in the supernova remnant Cassiopeia A. The spectra are fitted to give the electron temperature T(sub e), and (single) ionization age n(sub e)t. These quantities are compared with the predictions of self similar hydrodynamic models incorporating time dependent ionization and radiation losses, and Coulomb electron-ion equilibration behind the reverse shock, for a variety of different ejecta density profiles described by a uniform density core and a power law envelope. We find that the ejecta close to the 'jet' region in the NE, but not actually in the jet itself, have a systematically shallower outer envelope than ejecta elsewhere in the remnant, and we interpret this as being due to more energy of the initial explosion being directed in this polar direction as opposed to equatorially. The degree of asymmetry we infer is at the low end of that generally modelled in asymmetric core-collapse simulations, and may be used to rule out highly asymmetric explosion models.

Laming, J. Martin↗

Megafires in a Warming World: What Wildfire Risk Factors Led to California’s Largest Recorded Wildfire

Massive wildfires and extreme fire behavior are becoming more frequent across the westernUnited States, creating a need to better understand how megafire behavior will evolve in our warmingworld. Here, the fire spread model Prometheus is used to simulate the initial explosive growth ofthe 2020 August Complex, which occurred in northern California (CA) mixed conifer forests. Hightemperatures, low relative humidity, and daytime southerly winds were all highly correlated withextreme rates of modeled spread. Fine fuels reached very dry levels, which accelerated simulationgrowth and heightened fire heat release (HR). Model sensitivity tests indicate that fire growth andHR are most sensitive to aridity and fuel moisture content. Despite the impressive early observedgrowth of the fire, shifting the simulation ignition to a very dry September 2020 heatwave predicted a>50% increase in growth and HR, as well as increased nighttime fire activity. Detailed model analysesof how extreme fire behavior develops can help fire personnel prepare for problematic ignitions.

Kevin Varga↗

First Detection of Plasmoids from Breakout Reconnection on the Sun

Transient collimated plasma ejections (jets) occur frequently throughout the solar corona, in active regions, quiet Sun, and coronal holes. Although magnetic reconnection is generally agreed to be the mechanism of energy release in jets, the factors that dictate the location and rate of reconnection remain unclear. Our previous studies demonstrated that the magnetic breakout model explains the triggering and evolution of most jets over a wide range of scales, through detailed comparisons between our numerical simulations and high-resolution observations. An alternative explanation, the resistive-kink model, invokes breakout reconnection without forming and explosively expelling a flux rope. Here we report direct observations of breakout reconnection and plasmoid formation during two jets in the fan-spine topology of an embedded bipole. For the first time, we observed the formation and evolution of multiple small plasmoids with bidirectional flows associated with fast reconnection in 3D breakout current sheets (BCSs) in the solar corona. The first narrow jet was launched by reconnection at the BCS originating at the deformed 3D null, without significant flare reconnection or a filament eruption. In contrast, the second jet and release of cool filament plasma were triggered by explosive breakout reconnection when the leading edge of the rising flux rope formed by flare reconnection beneath the filament encountered the preexisting BCS. These observations solidly support both reconnection-driven jet models: the resistive kink for the first jet, and the breakout model for the second explosive jet with a filament eruption.

Kumar, Pankaj↗

Lunar and Planetary Science XXXV: Mars Impact Cratering

This document discusses the following topics: High Resolution Digital Elevation Models of Pristine Explosion Craters; Crater Degradation in the Martian Highlands: Morphometric Analysis of the Sinus Sabaeus Region and Simulation Modeling Suggest Fluvial Processes; Analysis of Impact Crater Preservation on Mars Using THEMIS Data; Atmospheric Entry Studies and the Smallest Impact Craters on Mars; Updating the Crater Count Chronology System for Mars; Control of Impact Crater-related Fracture Systems on the Subsurface Hydrology and Ground Collapse; Quantitative Analyses of Terrestrial Crater Deposits: Constraining Formation and Sediment Transport Processes on Mars; and Predicted Effects of Surface Processes on Martian Impact Crater Depth/Diameter Relationships

Stephen J Mackwell↗

Simulation of Blast Waves with Headwind

The blast wave resulting from an explosion was simulated to provide guidance for models estimating risks for human spacecraft flight. Simulations included effects of headwind on blast propagation, Blasts were modelled as an initial value problem with a uniform high energy sphere expanding into an ambient field. Both still air and cases with headwind were calculated.

Olsen, Michael E.↗

Postexplosion hydrodynamics of SN 1987A

A comprehensive study is presented of the postexplosion hydrodynamics of the ejecta of SN 1987A encompassing the Rayleigh-Taylor instabilities at the shell interfaces and the dynamical effects of the energy input due to the radioactive decay of Ni-56 and Co-56. The dependence of the hydrodynamical behavior on progenitor structure, initial perturbations, explosion energy, and dimensionality of the numerical simulation is investigated in order to explain the observed high-velocity wings of the iron line of SN 1987A. It is shown that none of these factors leads to the required high-velocity tail of iron. When nickel is artificially premixed in the initial model, it is found that this early mixing has to affect the inner 1.5 solar mass of the ejecta to reproduce the high-velocity wings at a time of 90 d. It is argued that the high-entropy bubble associated with the delayed mechanism is responsible for this premixing by driving convective motions early in the explosion.

Herant, Marc↗

Physics-Based Fragment Acceleration Modeling for Pressurized Tank Burst Risk Assessments

As part of comprehensive efforts to develop physics-based risk assessment techniques for space systems at NASA, coupled computational fluid and rigid body dynamic simulations were carried out to investigate the flow mechanisms that accelerate tank fragments in bursting pressurized vessels. Simulations of several configurations were compared to analyses based on the industry-standard Baker explosion model, and were used to formulate an improved version of the model. The standard model, which neglects an external fluid, was found to agree best with simulation results only in configurations where the internal-to-external pressure ratio is very high and fragment curvature is small. The improved model introduces terms that accommodate an external fluid and better account for variations based on circumferential fragment count. Physics-based analysis was critical in increasing the model's range of applicability. The improved tank burst model can be used to produce more accurate risk assessments of space vehicle failure modes that involve high-speed debris, such as exploding propellant tanks and bursting rocket engines.

Acceleration↗

The Mechanisms for the Onset and Explosive Eruption of Coronal Mass Ejections and Eruptive Flares

We have investigated the onset and acceleration of coronal mass ejections (CMEs) and eruptive flares. To isolate the eruption physics, our study uses the breakout model, which is insensitive to the energy buildup process leading to the eruption. We performed 2.5D simulations with adaptive mesh refinement that achieved the highest overall spatial resolution to date in a CME/eruptive flare simulation. The ultra-high resolution allows us to separate clearly the timing of the various phases of the eruption. Using new computational tools, we have determined the number and evolution of all X- and O-type nulls in the system, thereby tracking both the progress and the products of reconnection throughout the computational domain. Our results show definitively that CME onset is due to the start of fast reconnection at the breakout current sheet. Once this reconnection begins, eruption is inevitable; if this is the only reconnection in the system, however, the eruption will be slow. The explosive CME acceleration is triggered by fast reconnection at the flare current sheet. Our results indicate that the explosive eruption is caused by a resistive instability, not an ideal process. Moreover, both breakout and flare reconnections begin first as a form of weak tearing characterized by a slowly evolving plasmoids, but eventually transition to a fast form with well-defined Alfvenic reconnection jets and rapid flux transfer. This transition to fast reconnection is required for both CME onset and explosive acceleration. We discuss the key implications of our results for CME/flare observations and for theories of magnetic reconnection.

Karpen, Judith T.↗

Investigation of heat transfer in zirconium potassium perchlorate at low temperature: A study of the failure mechanism of the NASA standard initiator

The objective of this work was to study the reasons for the failure of pyrotechnic initiators at very low temperatures (10 to 100 K). A two-dimensional model of the NASA standard initiator was constructed to model heat transfer from the electrically heated stainless steel bridgewire to the zirconium potassium perchlorate explosive charge and the alumina charge cup. Temperature dependent properties were used in the model to simulate initiator performance over a wide range of initial temperatures (10 to 500 K). A search of the thermophysical property data base showed that pure alumina has a very high thermal conductivity at low temperatures. It had been assumed to act as a thermal insulator in all previous analyses. Rapid heat transfer from the bridgewire to the alumina at low initial temperatures was shown to cause failure of the initiators if the wire did not also make good contact with the zirconium potassium perchlorate charge. The mode is able to reproduce the results of the tests that had been conducted to investigate the cause for failure. It also provides an explanation for previously puzzling results and suggests simple design changes that will increase reliability at very low initial temperatures.

Varghese, Philip L.↗

Current loop coalescence model of solar flares

A computer simulation and theoretical study of the physical characteristics of the explosive coalescence of current-carrying loops is presented. Characteristics of the explosive coalescence include a large impulsive increase of the kinetic energies of electrons and ions, the simultaneous heating and acceleration of electrons and ions in high and low energy ranges, and a break in the energy spectra of electrons and ions. A characteristic double subpeak structure is found in the quasi-periodic oscillations found in the time profiles of the solar flares of June 7, 1980 and November 26, 1982 which can be explained in terms of the coalescence instability of two current loops.

Tajima, T.↗

Differential Muon Tomography to Continuously Monitor Changes in the Composition of Subsurface Fluids

Muon tomography has been used to seek hidden chambers in Egyptian pyramids and image subsurface features in volcanoes. It seemed likely that it could be used to image injected, supercritical carbon dioxide as it is emplaced in porous geological structures being used for carbon sequestration, and also to check on subsequent leakage. It should work equally well in any other application where there are two fluids of different densities, such as water and oil, or carbon dioxide and heavy oil in oil reservoirs. Continuous monitoring of movement of oil and/or flood fluid during enhanced oil recovery activities for managing injection is important for economic reasons. Checking on leakage for geological carbon storage is essential both for safety and for economic purposes. Current technology (for example, repeat 3D seismic surveys) is expensive and episodic. Muons are generated by high- energy cosmic rays resulting from supernova explosions, and interact with gas molecules in the atmosphere. This innovation has produced a theoretical model of muon attenuation in the thickness of rock above and within a typical sandstone reservoir at a depth of between 1.00 and 1.25 km. Because this first simulation was focused on carbon sequestration, the innovators chose depths sufficient for the pressure there to ensure that the carbon dioxide would be supercritical. This innovation demonstrates for the first time the feasibility of using the natural cosmic-ray muon flux to generate continuous tomographic images of carbon dioxide in a storage site. The muon flux is attenuated to an extent dependent on, amongst other things, the density of the materials through which it passes. The density of supercritical carbon dioxide is only three quarters that of the brine in the reservoir that it displaces. The first realistic simulations indicate that changes as small as 0.4% in the storage site bulk density could be detected (equivalent to 7% of the porosity, in this specific case). The initial muon flux is effectively constant at the surface of the Earth. Sensitivity of the method would be decreased with increasing depth. However, sensitivity can be improved by emplacing a greater array of particle detectors at the base of the reservoir.

Coleman, Max↗

The delayed-detonation model of Type Ia supernovae. 2: The detonation phase

The investigation, by use of two-dimensional numerical hydrodynamics simulations, of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is continued. Previously we found that the deflagration is insufficient to unbind the star. Expansion shuts off the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. The burning was allowed to develop into a detonation in these nonspherical models. The detonation grows toward spherical symmetry at late times. At these densities (rho approx. 10(exp 7) to 10(exp 8) g cm(exp -3)), either Ni-56 or nuclei of the Si-Ca group are the dominant products of the burning. The bulk yields are sensitive to the density of the star when the transition to detonation occurs. The relevance of the abundances, velocities, mixing, and total energy release to the theory and interpretation of Type Ia supernovae is discussed.

Arnett, David↗

INPIStron switched pulsed power for dense plasma pinches

The inverse plasma switch INPIStron was employed for 10kJ/40kV capacitor bank discharge system to produce focused dense plasmas in hypocycloidal-pinch (HCP) devices. A single unit and an array of multiple HCP's were coupled as the load of the pulsed power circuit. The geometry and switching plasma dynamics were found advantageous and convenient for commutating the large current pulse from the low impedance transmission line to the low impedance plasma load. The pulse power system with a single unit HCP, the system A, was used for production of high temperature plasma focus and its diagnostics. The radially running down plasma dynamics, revealed in image converter photographs, could be simulated by a simple snow-plow model with a correction for plasma resistivity. The system B with an array of 8-HCP units which forms a long coaxial discharge chamber was used for pumping a Ti-sapphire laser. The intense UV emission from the plasma was frequency shifted with dye-solution jacket to match the absorption band of the Ti crystal laser near 500 nm. An untuned laser pulse energy of 0.6 J/pulse was obtained for 6.4 kJ/40 kV discharge, or near 103 times of the explosion limit of conventional flash lamps. For both systems the advantages of the INPIStron were well demonstrated: a single unit is sufficient for a large current (greater than 50 kA) without increasing the system impedance, highly reliable and long life operation and implied scalability for the high power ranges above I(sub peak) = 1 MA and V(sub hold) = 100 kV.

Han, Kwang S.↗

Climatic Effects of Regional Nuclear War

We use a modern climate model and new estimates of smoke generated by fires in contemporary cities to calculate the response of the climate system to a regional nuclear war between emerging third world nuclear powers using 100 Hiroshima-size bombs (less than 0.03% of the explosive yield of the current global nuclear arsenal) on cities in the subtropics. We find significant cooling and reductions of precipitation lasting years, which would impact the global food supply. The climate changes are large and longlasting because the fuel loadings in modern cities are quite high and the subtropical solar insolation heats the resulting smoke cloud and lofts it into the high stratosphere, where removal mechanisms are slow. While the climate changes are less dramatic than found in previous "nuclear winter" simulations of a massive nuclear exchange between the superpowers, because less smoke is emitted, the changes seem to be more persistent because of improvements in representing aerosol processes and microphysical/dynamical interactions, including radiative heating effects, in newer global climate system models. The assumptions and calculations that go into these conclusions will be described.

Oman, Luke D.↗

Star formation in proto dwarf galaxies

The effects of the onset of star formation on the residual gas in primordial low-mass Local-Group dwarf spheroidal galaxies is studied by a series of hydrodynamical simulations. The models have concentrated on the effect of photoionization. The results indicate that photoionization in the presence of a moderate gas density gradient can eject most of the residual gas on a time scale of a few 10 to the 7th power years. High central gas density combined with inefficient star formation, however, may prevent mass ejection. The effect of supernova explosions is discussed briefly.

Noriega-Crespo, A.↗

Reply to Comment by Cole-dai Et Al. on "Climatic Impact of the Long-lasting Laki Eruption: Inapplicability of Mass-independent Sulfur Isotope Composition Measurements"

Here we respond to the comments by Cole-Dai et al. [2014] on our article Schmidt et al. [2012]. Specifically, in response to section 2 of their reply, we argued in Schmidt et al. [2012] that based on previously published estimates of the volatile release height during the 1783-1784 C.E. Laki eruption, the lack of a sulfur massindependent fractionation (MIF) anomaly is expected. In other words, no previous study on Laki ever argued that this eruption emitted SO2 into altitudes >13-15 km. In section 2.3, Cole-Dai et al. [2014] argue that the nonzero Δ33S value of their Laki sample 1 may be explained by a short-lived explosive phase at Laki during which volatiles reached the stratosphere. In Schmidt et al. [2012] in section 2, we argued in agreement with Cole-Dai et al. [2014] (section 3.1) that for a MIF anomaly to be preserved, the Laki volatiles would have had to be emitted in >20 km altitude. Our main point is that eruption column heights >20 km are unlikely based on the historical accounts and plume-rise modeling for the Laki eruption [Stothers et al., 1986; Woods, 1993; Thordarson and Self, 2003]. In Schmidt et al. [2012], we argued that to deduce a short-lived climatic impact of the Laki eruption based on the lack of a MIF anomaly and the length of the sulfate deposition in Greenland ice cores may be misleading because the climatic impact will outlast the radiative forcing of the Laki aerosol cloud. Cole-Dai et al. [2014] acknowledge the latter in their reply in section 4.2. We agreewith Cole-Dai et al. [2014] in that themagnitude and length of the climatic impact during the winter of 1783-1784 depends on the altitude of the volatile release during the eruption (sections 2.3 and 4.2). However, even if we assumed that during Laki all sulfur dioxide (SO2) would have been released in the troposphere, then the aerosol cloud would still be present in the upper troposphere during March 1784, as is evident from independent model simulations of this "tropospheric-only" scenario [Stevenson et al., 2003]. We acknowledge that there is uncertainty on the volatile release height for Laki; however, it is worth considering that those climate model simulations that used an injection altitude between 9 km and 13 km for the Laki SO2 [Highwood and Stevenson, 2003; Oman et al., 2006a, 2006b; Schmidt et al., 2012] best match the observed temperature changes during summer of 1783 [Angell and Korshover, 1985; Brázdil et al., 2010; Briffa et al., 1998; D'Arrigo and Jacoby, 1999; Jacoby et al., 1999; Kington, 1988; Manley, 1974; Parker et al., 1992; Thordarson and Self, 2003]. Based on these model simulations, a climatic impact during the winter of 1783-1784, albeit weaker than during the climactic phases of Laki, is expected (and our argument here does not exclude the role of natural variability in contributing to the cold winter of 1783-1784 as discussed in Schmidt et al. [2012]). Therefore, we continue to argue that for high-latitude eruptions such as Laki, the applicability of sulfur isotopic measurements to interpret the climatic relevance has yet to be demonstrated. Itmay transpire that the interpretation of MIF signals for the climate-relevance of an eruption is valid and unambiguous only for short-lived explosive eruptions in the tropics. In terms of the processes producing a MIF anomaly (section 3.3 in Cole-Dai et al. [2014]), the works by Hattori et al. [2013] and Ono et al. [2013] suggest that there are remaining issues not discussed by Cole-Dai et al. [2014], for instance, self-shielding of SO2 due to high column densities typical for eruptions of Pinatubo-scale and greater, and the preservation of the MIF signature in general.

general circulation model↗

Earth Satellite Population Instability: Underscoring the Need for Debris Mitigation

A recent study by NASA indicates that the implementation of international orbital debris mitigation measures alone will not prevent a significant increase in the artificial Earth satellite population, beginning in the second half of this century. Whereas the focus of the aerospace community for the past 25 years has been on the curtailment of the generation of long-lived orbital debris, active remediation of the current orbital debris population should now be reconsidered to help preserve near-Earth space for future generations. In particular, we show in this paper that even if launch operations were to cease today, the population of space debris would continue to grow. Further, proposed remediation techniques do not appear to offer a viable solution. We therefore recommend that, while the aerospace community maintains the current debris-limiting mission regulations and postmission disposal procedures, future emphasis should be placed on finding new remediation technologies for solving this growing problem. Since the launch of Sputnik 1, space activities have created an orbital debris environment that poses increasing impact risks to existing space systems, including human space flight and robotic missions (1, 2). Currently, more than 9,000 Earth orbiting man-made objects (including many breakup fragments), with a combined mass exceeding 5 million kilograms, are tracked by the US Space Surveillance Network and maintained in the US satellite catalog (3-5). Three accidental collisions between cataloged satellites during the period from late 1991 to early 2005 have already been documented (6), although fortunately none resulted in the creation of large, trackable debris clouds. Several studies conducted during 1991-2001 demonstrated, with assumed future launch rates, the unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects (7-13). In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new satellites (i.e., debris) due to collisions exceeds the loss of objects due to orbital decay. NASA s evolutionary satellite population model LEGEND (LEO-to-GEO Environment Debris model), developed by the Orbital Debris Program Office at the NASA Lyndon B. Johnson Space Center, is a high fidelity three-dimensional physical model that is capable of simulating the historical satellite environment, as well as the evolution of future debris populations (14, 15). The subject study assumed no rocket bodies and spacecraft were launched after December 2004, and no future disposal maneuvers were allowed for existing spacecraft, few of which currently have such a capability. The rate of satellite explosions would naturally decrease to zero within a few decades as the current satellite population ages. The LEGEND future projection adopts a Monte Carlo approach to simulate future on-orbit explosions and collisions. Within a given projection time step, once the explosion probability is estimated for an intact object, a random number is drawn and compared with the probability to determine if an explosion would occur. A similar procedure is applied to collisions for each pair of target and projectile involved within the same time step. Due to the nature of the Monte Carlo process, multiple projection runs must be performed and analyzed before one can draw reliable and meaningful conclusions from the outcome. A total of fifty, 200-year future projection Monte Carlo simulations were executed and evaluated (16).

Liou, Jer-chyi↗