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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Exploration of a Novel Technique for Waste Heat Recovery Through Molecular Dynamics: Influence of Wettability and Electric Field on Water and Water-Based Nanofluids

Most of the energy produced globally comes by way of a heat engine. The Carnot principle places a limit as to how thermodynamically efficient a heat engine can be. There is no heat engine that can be 100% thermodynamically efficient and as such a substantial proportion of all heat supplied to a heat engine is lost as waste heat. Waste heat therefore is a large energy source ready to be properly utilized. Herein, a novel approach for converting waste heat to electricity is discussed. It involves the use of the liquid to vapor phase change of a material dielectric (water) or electrolyte (nanofluid) in the embodiment of a capacitor for direct thermal to electrostatic energy conversion. While this method of waste heat recovery could potentially be added to the ever expanding portfolio of energy conversion techniques, a number of aspects must be addressed before it can be brought into practice. Water was seen as an ideal dielectric phase change material given its high relative permittivity ratio when in the liquid form as compared to its vapor form. However, given its short voltage holdoff time the phase change of water would need to occur rapidly. This brings up concerns of explosive boiling. Herein, molecular dynamics analysis into the explosive boiling behavior of thin water films gave more insight into how the interaction between the surface and liquid affected explosive boiling onset time. A Lennard-Jones potential with one interaction site and a Morse potential with three interaction sites between water and solid substrate were used. It was found generally that a stronger interaction between water film and substrate led to faster explosive boiling onset times but an increase in the number of interaction sites delayed explosive boiling, even at the same wettability (contact angle). Understanding changes in the density and enthalpy of vaporization of a liquid dielectric such as water in the presence of an electric field is of importance due to the electrostatic nature of the waste heat conversion method under consideration. Specifically, if both density and enthalpy of vaporization are increased, the thermodynamic efficiency of the waste heat conversion method under consideration is decreased. Electric field effects are explored herein via molecular dynamics using two water models, the TIP4P-Ew and SWM4-NDP. The SWM4-NDP model is polarizable while the TIP4P-Ew model is not, which allows for a determination of the importance of model polarizability (i.e. variation in water model dipole moment) on these two properties of water when subjected to an electric field. Herein it was found that both water models respond similarly in terms of density and vaporization enthalpy variance upon the introduction of an electric field. Comparison was also made to the pressure induced by the electric field (electrostriction pressure) by way of a density comparison and it was found that the predicted electrostriction pressure overestimates the pressure experienced by water. Water by itself has a high enthalpy of vaporization, which limits the efficiency of the newly proposed conversion method. Research both experimental and through simulation has shown that the vaporization enthalpy of nanofluids can be engineered via nanoparticle size and material selection. An avenue less explored is manipulating the enthalpy of vaporization by altering the interaction strength between the nanoparticles and the base fluid. In practice this could be achieved through the addition of coatings to the nanoparticles to alter their wettability to the base fluid. This was explored by using a Lennard-Jones potential and Morse potential to model the interaction between base fluid (water) and the nanoparticle. For nanoparticles 2nm in diameter and at weight percentages up to 6%, the change in vaporization enthalpy due to alterations of the interaction strength between nanoparticle and base fluid was not significant (less than a 1% difference) when compared to the effect of altering the weight percentage of nanoparticles in the nanofluid or introducing an electric field. However, the effect of wettability may still become important at other nanoparticle concentrations and sizes. In all, the studies presented here further the understanding of phase change and thermodynamic properties of water and water based nanofluids under an electrostatic field which will help inform the development of a novel approach to waste heat conversion. The reduction of waste heat will improve energy sustainability outlooks.

30 DIRECT ENERGY CONVERSION↗

NASA NTRS1999

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

NASA NTRS2021

Particle Acceleration by Dense Impulsive Structures Moving in Ambient Magnetospheric Plasma. 3-D Hybrid Kinetic Modeling and MMS Observations

High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.

A. S. Lipatov↗

NASA NTRS2012

Probing Pre-Galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature approximately greater than10(exp 4) K.We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1 + z) with values of about a few cm(exp −3). In more massive halos, corresponding to the first galaxies, the density may be larger, n approximately greater than100 cm(exp −3). The resulting afterglow fluxes are weakly dependent on redshift at a fixed observed time, and may be detectable with the James Webb Space Telescope and Very Large Array in the near-IR and radio wavebands, respectively, out to redshift z approximately greater than 20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. The metal absorption line signature is expected to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages (C II, O I, Si II and Fe II). We show that, if the afterglow can be observed sufficiently early, analysis of the metal lines may distinguish whether the first heavy elements were produced in a pair-instability supernova or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

Metal↗

NASA NTRS1996

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗

DOE OSTI2020

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2025

Assessing the Potential Impact of Fugitive Methane Emissions on Offshore Platform Safety

One of the biggest risks to safety on offshore platform safety is the ignition of high-pressure natural gas streams. Currently, the size and number of fugitive emissions on offshore platforms is unknown and methods used to detect fugitives have significant shortcomings. To investigate the frequency, size, and potential impact of fugitives, a data collection exercise was conducted using incidents reported, leak survey data, and independent measurements. The size and number of fugitives on offshore facilities were simulated to investigate likely areas of safety concern. Incident reports indicate in 2021 there were 113 reports of gas leaks on 1119 offshore facilities, suggesting 0.02 fugitives per Type 1 facility (older, shallow-water platforms) and 0.31 fugitives per Type 2 facility (larger deeper-water facilities). Leak survey data report 12 fugitives per Type 1 facility (average emission 0.6 kg CH 4 h −1 leak −1 ) and 15 fugitives per Type 2 facility (average emission 1.5 kg CH 4 h −1 leak −1 ). Reconciliation of direct measurements with a bottom-up model suggests that the number of fugitive emissions generated from the leak report data is an underestimate for Type 1 platforms (44 fugitives facility −1 ; average emission 0.6 kg CH 4 h −1 leak −1 ) and in general agreement for the Type 2 platforms (15 fugitives facility −1 ; average emission 1.5 kg CH 4 h −1 leak −1 ). Analysis of the fugitive emission rates on an offshore platform suggests that gas will not collect to explosive concentration if any air movement is present (>0.36 mph); however, large volumes of air (~600 m 3 ) near representative leaks on the working deck could become explosive in hour-long zero-wind conditions. We suggest that wearable technology could be employed to indicate gas build up, safety regulations amended to consider low-wind conditions and real-world experiments are conducted to test assumptions of air mixing on the working deck.

explosion↗

DOE OSTI2024

A Theory for Neutron Star and Black Hole Kicks and Induced Spins

Abstract Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s −1 range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s −1 range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10 M ⊙ , is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0 M ⊙ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

Astronomy & Astrophysics↗

DOE OSTI2022

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗

NASA NTRS2012

The Effect of Magnetic Topology on the Escape of Flare Particles

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗

NASA NTRS2012

A Model for SEP Escape

Antiochos, S. K.↗

NASA NTRS2010

Hidden Statistics Approach to Quantum Simulations

Recent advances in quantum information theory have inspired an explosion of interest in new quantum algorithms for solving hard computational (quantum and non-quantum) problems. The basic principle of quantum computation is that the quantum properties can be used to represent structure data, and that quantum mechanisms can be devised and built to perform operations with this data. Three basic non-classical properties of quantum mechanics superposition, entanglement, and direct-product decomposability were main reasons for optimism about capabilities of quantum computers that promised simultaneous processing of large massifs of highly correlated data. Unfortunately, these advantages of quantum mechanics came with a high price. One major problem is keeping the components of the computer in a coherent state, as the slightest interaction with the external world would cause the system to decohere. That is why the hardware implementation of a quantum computer is still unsolved. The basic idea of this work is to create a new kind of dynamical system that would preserve the main three properties of quantum physics superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. In other words, such a system would reinforce the advantages and minimize limitations of both quantum and classical aspects. Based upon a concept of hidden statistics, a new kind of dynamical system for simulation of Schroedinger equation is proposed. The system represents a modified Madelung version of Schroedinger equation. It preserves superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. Such an optimal combination of characteristics is a perfect match for simulating quantum systems. The model includes a transitional component of quantum potential (that has been overlooked in previous treatment of the Madelung equation). The role of the transitional potential is to provide a jump from a deterministic state to a random state with prescribed probability density. This jump is triggered by blowup instability due to violation of Lipschitz condition generated by the quantum potential. As a result, the dynamics attains quantum properties on a classical scale. The model can be implemented physically as an analog VLSI-based (very-large-scale integration-based) computer, or numerically on a digital computer. This work opens a way of developing fundamentally new algorithms for quantum simulations of exponentially complex problems that expand NASA capabilities in conducting space activities. It has been illustrated that the complexity of simulations of particle interaction can be reduced from an exponential one to a polynomial one.

Zak, Michail↗

DOE OSTI2022

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗

NASA NTRS2023

Orbital Debris Shape Effect Investigations for Mitigating Risk

NASA’s Orbital Debris Program Office (ODPO) and Hypervelocity Impact Technology (HVIT) team have coordinated to better understand the risks to upper stages and spacecraft from non-spherical orbital debris. It is well understood that fragmentation (collision or explosion) events in orbit produce fragments of various materials, sizes, and shapes. To further characterize these parameters, the ODPO is developing the next-generation Orbital Debris Engineering Model (ORDEM) version 4.0 to include orbital debris shape distributions. Ground-based assets, such as radar and optical sensors, can provide size estimates and some insight into material based on radar return or optical filter photometry/spectroscopy, respectively. Characterizing an object’s shape requires more laboratory analyses to infer how shape affects these measurements. More importantly, in addition to size and material/density, the shape of fragments in orbit will alter the ballistic limit equations used in orbital debris risk assessments with NASA’s Bumper Code. The ODPO plans to release ORDEM 4.0 in the coming years. Performing ground-based laboratory impact tests on high-fidelity spacecraft mockups provides the means to directly measure size, mass, material/density, and shape of fragments, all key parameters needed to characterize real-world break up events. The DebriSat test, the results of which are provided, showcases the details of this type of experiment. The goal of this collaborative research between the ODPO and the HVIT team is to include a shape parameter in the environmental and breakup models used to assess risk for various space structures. This paper examines ground-based laboratory impact tests and the associated fragment shape categories. Provided these defined shapes, the approach is simplified by assuming a right circular cylinder (RCC) approximation with varying length-to-diameter ratios. Highlights of impact tests conducted by the HVIT team using non-spherical projectiles based on the RCC approximation are presented. Hydrocode simulations have also been performed to expand on the complexity of variations with non-spherical projectiles. Lastly, ray-tracing simulations of various RCCs of known material are provided to support the ongoing research on optical reflectance distributions with known shapes and to highlight how this may modify the current optical size estimation model. The status and plan forward are outlined for NASA's orbital debris shape effect investigation using a multidisciplinary approach by the ODPO and the HVIT team.

Heather Cowardin↗

NASA NTRS2023

Orbital Debris Shape Effect Investigations for Mitigating Risk

Heather Cowardin↗

DOE OSTI2024

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗