SEARCH · Engineering Papers

Results for “dynamic shock experiments”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1986

A shock capturing technique for hypersonic, chemically relaxing flows

A fully coupled, shock capturing technique is presented for chemically reacting flows at high Mach numbers. The technique makes use of a total variation diminishing (TVD) dissipation operator which results in sharp, crisp shocks. The eigenvalues and eigenvectors of the fully coupled system, which includes species conversion equations in addition to the gas dynamics equations, are analytically derived for a general reacting gas. Species production terms for a model dissociating gas are introduced and are included in the algorithm. The convective terms are solved using a first-order TVD scheme while the source terms are solved using a fourth-order Runge-Kutta scheme to enhance stability. Results from one-dimensional numerical experiments are shown for a two species and a three species gas.

Eberhardt, S.↗

NASA NTRS2009

Proceedings of the 40th Lunar and Planetary Science Conference

The 40th Lunar and Planetary Science Conference included sessions on: Phoenix: Exploration of the Martian Arctic; Origin and Early Evolution of the Moon; Comet Wild 2: Mineralogy and More; Astrobiology: Meteorites, Microbes, Hydrous Habitats, and Irradiated Ices; Phoenix: Soil, Chemistry, and Habitability; Planetary Differentiation; Presolar Grains: Structures and Origins; SPECIAL SESSION: Venus Atmosphere: Venus Express and Future Missions; Mars Polar Caps: Past and Present; SPECIAL SESSION: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1, Part I; 5 Early Nebula Processes and Models; SPECIAL SESSION: Icy Satellites of Jupiter and Saturn: Cosmic Gymnasts; Mars: Ground Ice and Climate Change; SPECIAL SESSION: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1, Part II; Chondrite Parent-Body Processes; SPECIAL SESSION: Icy Satellites of Jupiter and Saturn: Salubrious Surfaces; SNC Meteorites; Ancient Martian Crust: Primary Mineralogy and Aqueous Alteration; SPECIAL SESSION: Messenger at Mercury: A Global Perspective on the Innermost Planet; CAIs and Chondrules: Records of Early Solar System Processes; Small Bodies: Shapes of Things to Come; Sulfur on Mars: Rocks, Soils, and Cycling Processes; Mercury: Evolution and Tectonics; Venus Geology, Volcanism, Tectonics, and Resurfacing; Asteroid-Meteorite Connections; Impacts I: Models and Experiments; Solar Wind and Genesis: Measurements and Interpretation; Mars: Aqueous Processes; Magmatic Volatiles and Eruptive Conditions of Lunar Basalts; Comparative Planetology; Interstellar Matter: Origins and Relationships; Impacts II: Craters and Ejecta Mars: Tectonics and Dynamics; Mars Analogs I: Geological; Exploring the Diversity of Lunar Lithologies with Sample Analyses and Remote Sensing; Chondrite Accretion and Early History; Science Instruments for the Mars Science Lander; . Martian Gullies: Morphology and Origins; Mars: Dunes, Dust, and Wind; Mars: Volcanism; Early Solar System Chronology; Seek Out and Explore: Upcoming and Future Missions; Mars: Early History and Impact Processes; Mars Analogs II: Chemical and Spectral; Achondrites and their Parent Bodies; and Planning for Future Exploration of the Moon The poster sessions were: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1; LRO and LCROSS; Geophysical Analysis of the Lunar Surface and Interior; Remote Observation and Geologic Mapping of the Lunar Surface; Lunar Spectroscopy; Venus Geology, Geophysics, Mapping, and Sampling; Planetary Differentiation; Bunburra and Buzzard Coulee: Recent Meteorite Falls; Meteorites: Terrestrial History; CAIs and Chondrules: Records of Early Solar System Processes; Volatile and Organic Compounds in Chondrites; Crashing Chondrites: Impact, Shock, and Melting; Ureilite Studies; Petrology and Mineralogy of the SNC Meteorites; Martian Meteorites; Phoenix Landing Site: Perchlorate and Other Tasty Treats; Mars Polar Atmospheres and Climate Modeling; Mars Polar Investigations; Mars Near-Surface Ice; Mars: A Volatile-Rich Planet; Mars: Geochemistry and Alteration Processes; Martian Phyllosilicates: Identification, Formation, and Alteration; Astrobiology; Instrument Concepts, Systems, and Probes for Investigating Rocks and Regolith; Seeing is Believing: UV, VIS, IR, X- and Gamma-Ray Camera and Spectrometer Instruments; Up Close and Personal: In Situ Analysis with Laser-Induced Breakdown Spectroscopy and Mass Spectrometry; Jupiter and Inscrutable Io; Tantalizing Titan; Enigmatic Enceladus and Intriguing Iapetus; Icy Satellites: Cryptic Craters; Icy Satellites: Gelid Geology/Geophysics; Icy Satellites: Cool Chemistry and Spectacular Spectroscopy; Asteroids and Comets; Comet Wild 2: Mineralogy and More; Hypervelocity Impacts: Stardust Models, LDEF, and ISPE; Presolar Grains; Early Nebular Processes: Models and Isotopes; Solar Wind and Genesis: Measurements and Interpretation; Education and Public Outreach; Mercury; Pursuing Lunar Exploration; Sources and Eruptionf Lunar Basalts; Chemical and Physical Properties of the Lunar Regolith; Lunar Dust and Transient Surface Phenomena; Lunar Databases and Data Restoration; Meteoritic Samples of the Moon; Chondrites, Their Clasts, and Alteration; Achondrites: Primitive and Not So Primitive; Iron Meteorites; Meteorite Methodology; Antarctic Micrometeorites; HEDs and Vesta; Dust Formation and Transformation; Interstellar Organic Matter; Early Solar System Chronology; Comparative Planetology; Impacts I: Models and Experiments; Impacts II: Craters and Ejecta; Mars: Volcanism; Mars: Tectonics and Dynamics; Martian Stratigraphy: Understanding the Geologic History of Mars Through the Sedimentary Rock Record; Mars: Valleys and Valley Networks; Mars: Aqueous Processes in Valles Marineris and the Southern Highlands; Mars: Aqueous Geomorphology; Martian Gullies: Morphology and Origins; Mars: Dunes, Dust, and Wind; Mars: Remote Sensing; Mars: Geologic Mapping, Photogrammetry, and Cratering; Martian Mineralogy: Constraints from Missions and Laboratory Investigations; Mars Analogs: Chemical and Physical; Mars Analogs: Sulfates and Sulfides; Missions: Approaches, Architectures, Analogs, and Actualities; Not Just Skin Deep: Electron Microscopy, Heat Flow, Radar, and Seismology Instruments and Planetary Data Systems, Techniques, and Interpretation.

Source record↗

DOE OSTI2023

Performance Characteristics of a Rotating Detonation Combustor Exiting Into a Pressurized Plenum to Simulate Gas Turbine Inlet

Abstract The present study aims to experimentally characterize the performance of a rotating detonation combustion (RDC) system integrated with a pressurized downstream plenum to simulate the high-pressure inlet conditions of power-generating gas turbines. A thorough understanding of the operational behavior including wave mode behavior, static pressure profile along the combustor length, and dynamic features of pressure fluctuations is crucial for successful integration of RDC with the turbine. In this study, two RDC configurations are investigated, RDC with a constant area annulus and RDC with a converging nozzle. In both cases, the RDC flow exited into a plenum chamber kept at pressures varying from 155 kPa to 330 kPa. RDC was operated on methane and oxygen-enriched air to represent reactants used in land-based power generation. Experiments were conducted for the two RDCs configurations operated at three reactant mass flow rates (0.23, 0.32, and 0.46 kg/s). The RDC performance is characterized by time-averaged static pressure measurements, and wave velocity determined by ionization probes. In addition, dynamic pressure measurements were recorded both inside and near the exit of RDC channel to investigate wave interactions between RDC and downstream plenum. Results show that the RDC with the converging nozzle achieved superior performance while minimizing detrimental interactions with the reflected shock and/or acoustic waves from the downstream plenum.

Engineering↗

DOE OSTI2020

A spherical crystal diffraction imager for Sandia’s Z Pulsed Power Facility

Sandia’s Z Pulsed Power Facility is able to dynamically compress matter to extreme states with exceptional uniformity, duration, and size, which are ideal for investigating fundamental material properties of high energy density conditions. X-ray diffraction (XRD) is a key atomic scale probe since it provides direct observation of the compression and strain of the crystal lattice and is used to detect, identify, and quantify phase transitions. Because of the destructive nature of Z-Dynamic Material Property (DMP) experiments and low signal vs background emission levels of XRD, it is very challenging to detect a diffraction signal close to the Z-DMP load and to recover the data. We have developed a new Spherical Crystal Diffraction Imager (SCDI) diagnostic to relay and image the diffracted x-ray pattern away from the load debris field. The SCDI diagnostic utilizes the Z-Beamlet laser to generate 6.2-keV Mn–He α x rays to probe a shock-compressed material on the Z-DMP load. Finally, a spherically bent crystal composed of highly oriented pyrolytic graphite is used to collect and focus the diffracted x rays into a 1-in. thick tungsten housing, where an image plate is used to record the data.

47 OTHER INSTRUMENTATION↗

NASA NTRS1979

Data reduction and analysis for the TRW IMP-7 plasma wave experiment

Highlights from the IMP 7 plasma wave experiment are briefly outlined. The measurements in the tail revealed great complexities involving substorms and fireballs and provided insight into the overall dynamics of the magnetosphere. The low levels of magnetic turbulence may be of significance with respect to the development of tearing mode dissipation for reconnection events. The IMP wave observations in the distant magnetosphere boundary region showed that the boundaries were quite diffuse and frequently purely defined. Finally, nearly simultaneous IMP 7 and 8 comparisons were used to achieve an understanding of shock structure and propagation.

Scarf, F. L.↗

DOE OSTI2020

(U) PRAD0697 & PRAD0698: Complex Loading of CeO₂ Powder

Cerium(IV) oxide (CeO₂) powder is shock compressed using the Precision High Energy-density Liner Implosion eXperiment (PHELIX) platform. Experimental results are compared against several modeling approaches. Compaction behavior is best captured with a P-∝ model, which calculates CeO₂ powder bulk densities within 80-99% of experimental values but overpredicts densi cation at the cylindrical target's outer radius and center by up to 20%. Preliminary calculations suggest that accuracy could be increased with the inclusion of a coupled strength model. Several common computational modeling approaches for the shock compression response of granular materials and the magnetohydrodynamic (MHD) force upon the impactor/liner in pulsed power compression experiments are investigated and analyzed for their validity. The Bi-linear Ramp, P-∝ PACXP, and P-∝ Menikoff-Kober continuum compaction models are calibrated to planar impact Hugoniot data for CeO₂ powder and used to predict the powder's shock compaction response under non-planar shock wave compression. MHD calculations of the PHELIX pulsed power driver are performed using an idealized resistor-inductor-capacitor (RLC) circuit calibrated to previous experiments. All simulations are performed using the LANL code FLAG. Two validation experiments are computationally designed using the calibrated compaction and circuit models, executed using the PHELIX platform on CeO₂ targets with initial porous densities of 3.95 and 4.03 g/cm³, measured with proton radiography, and analyzed against the model predictions. The two P-∝ models more accurately describe CeO₂ powder densi cation than the Bi-linear Ramp model. However, the two P-∝ models overpredict bulk density of the shock compressed CeO₂ powder by up to 20% when the appropriate impact velocities are applied. MHD calculations for both validation experiments underpredict liner impact velocities by 4-11% when using the idealized RLC circuit model calibrated to previous experiments. Compensating underpredictions of impact velocity and overpredictions of powder densication lead to a false accuracy in pre-shot calculations compared to experimental data. To improve correlation between simulations and experiments, the following improvements are suggested: 1. A coupled strength model for CeO₂ powder that updates strength as a function of porosity and applied stress. 2. An improved MHD circuit model that more accurately captures the PHELIX machine.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Multimodal 3D characterization of voids in shock-loaded tantalum: Implications for ductile spallation mechanisms

Predicting the failure of crystalline materials at high strain rates requires knowledge of the underlying failure mechanisms and their dependence on microstructure. However, little experimental consensus exists on the underlying micromechanics of spallation. In this study, a 3D-EBSD characterization experiment is performed on tantalum prior to and after partial spallation by plate impact, which allowed for the statistical assessment of the microstructural neighborhoods surrounding incipient voids. In analyzing the resulting dataset containing 5884 grains and 467 voids, it is observed that the voids were roughly spherical and consistent in size throughout the spalled material. The voids mostly resided at quadruple points, at triple junctions, at grain boundaries, and within grains, in decreasing order of prevalence. Moreover, voids tended to form at grain boundaries with high degrees of plastic incompatibility, growing into the plastically soft grain but orienting primarily with or perpendicular to the loading direction. Here, the statistics from these analyses of 3D microstructural data support dynamic cavitation models for ductile spallation.

36 MATERIALS SCIENCE↗

NASA NTRS1984

Second- and third-order upwind difference schemes for hyperbolic conservation laws

Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.

Yang, J. Y.↗

DOE OSTI2025

Simulation and measurement of the dynamic strain response in a prototypical spallation target

The Second Target Station (STS) at the Spallation Neutron Source (SNS) will address emerging scientific challenges by providing a source of intense cold neutrons to instruments optimized for this source. The STS target will use rotating tungsten blocks and will receive 1.3 GeV proton beam pulses from the SNS accelerator at a repetition rate of 15 Hz. The facility life is planned for 40 years, and each target assembly life is expected to be approximately 10 years. An accurate strain prediction is then critical for fatigue life assessment of STS target blocks because they will be subject to approximately 10 8 beam pulses per lifetime. As an R&D activity, the Los Alamos Neutron Science Center (LANSCE) Weapons Neutron Research (WNR) Target 2 (Blue Room) facility was used to test the strain response of prototypical target blocks to the thermal shock of a proton pulse. The blue room was well suited for a pulsed proton beam impact test of subscale STS target blocks; the 800 MeV proton energy is approximately 60 % of the 1.3 GeV proton energy expected from the SNS accelerator to the STS. The LANSCE Proton Storage Ring (PSR) and SNS are both short-pulse proton beam sources with nominal pulse widths of 250 ns and 661 ns, respectively, so the energy deposition in the target occurs in <1 μs pulse duration. Strain measurements on the outer surface of three target blocks (bare tungsten, tantalum-clad tungsten, niobium-clad tungsten) were recorded for comparison against neutronics and structural simulations. In conclusion, this experiment and the supporting simulations satisfied the following primary research goals for the STS target.

Dynamic strain response↗

DOE OSTI2021

Dynamic compression of TiO 2 to 221 GPa

The high-pressure dynamic response of titanium dioxide (TiO 2 ) is not only of interest because of its numerous industrial applications but also because of its structural similarities to silica (SiO 2 ). We performed plate impact experiments in a two-stage light gas gun, at peak stresses from 64 to 221 GPa to determine the TiO 2 response along the Hugoniot. The lower stress experiment at 64 GPa shows an elastic behavior followed by an elastic–plastic transition, whereas the high stress experiments above 64 GPa show a single wave structure. Previous shock studies have shown the presence of high-pressure phases (HPP) I (26 GPa) and HPP II (100 GPa); however, our data suggest that the HPP I phase is stable up to 150 GPa. Using a combination of data from our current study and our previous Z-data, we determine that TiO 2 likely melts on the Hugoniot at 157 GPa. Furthermore, our data confirm that TiO 2 is not highly incompressible as shown by a previous study.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗

NASA NTRS1947

Pressure Distribution Measurements at High Speed and Oblique Incidence of Flow

The present report contains the results of a series of observations obtained for a wing of symmetrical profile for different angles of yaw. The shock tunnel with 0.4m x 0.4m cross section, of the Institute for Gas Dynamics at L.F.A. - Braunchweig was available for this work. The profile used was a 9 percent, thick hyperbolic profile with maximum thickness 40 percent aft, which was calculated by the method of F. Ringleb (2). Since the conformal transformation of this profile is known, the theoretical pressure distribution could be determined exactly for the case of incompressible flow. Then by the use of the Prandtl rule one may extend the comparison of theory and experiment to the case of higher velocity of incident flow.

Lippisch, A.↗

NASA NTRS2001

Impact and Collisional Processes in the Solar System

In the past year, we have successfully developed the techniques necessary to conduct impact experiments on ice at very low temperatures. We employ the method of embedding gauges within a target to measure the shock wave and material properties. This means that our data are not model dependent; we directly measure the essential parameters needed for numerical simulations of impact cratering. Since then we have developed a new method for temperature control of icy targets that ensures temperature equilibrium throughout a porous target. Graduate student, Sarah Stewart-Mukhopadhyay, is leading the work on ices and porous materials as the main thrust of her thesis research. Our previous work has focused on icy materials with no porosity, and we propose to extend our research to include porous ice and porous ice-silicate mixtures. There is little shockwave data for porous ice, and none of the data was acquired under conditions applicable to the outer solar system. The solid ice Hugoniot is only defined for initial temperatures above -20 C. Our program uniquely measures the properties of ice at temperatures directly applicable to the solar system. Previous experiments were conducted at ambient temperatures soon after removing the target from a cold environment, usually just below freezing, or in a room just below freezing. Since ice has an extremely complicated phase diagram, it is important to conduct experiments at lower temperatures to determine the true outcome of impacts in the outer solar system. This research is complementary to other programs on icy materials. Our work focuses on the inherent material properties by measuring the shock wave directly; this complements the macroscopic observations and immediately provides the parameters necessary to extend this research to the gravity regime. Our numerical simulations of impacts in porous ice under very low gravity conditions, such as found on comets, show that the final crater size and shape is very dependent on the dynamic strength of the material.

Ahrens, Thomas J.↗

NASA NTRS2022

Dayside Transient Phenomena and Their Impact on the Magnetosphere and Ionosphere

Dayside transients, such as hot flow anomalies, foreshock bubbles, magnetosheath jets, flux transfer events, and surface waves, are frequently observed upstream from the bow shock, in the magnetosheath, and at the magnetopause. They play a significant role in the solar wind-magnetosphere-ionosphere coupling. Foreshock transient phenomena, associated with variations in the solar wind dynamic pressure, deform the magnetopause, and in turn generates field-aligned currents (FACs) connected to the auroral ionosphere. Solar wind dynamic pressure variations and transient phenomena at the dayside magnetopause drive magnetospheric ultra low frequency (ULF) waves, which can play an important role in the dynamics of Earth’s radiation belts. These transient phenomena and their geoeffects have been investigated using coordinated in-situ spacecraft observations, spacecraft-borne imagers, ground-based observations, and numerical simulations. Cluster, THEMIS, Geotail, and MMS multi-mission observations allow us to track the motion and time evolution of transient phenomena at different spatial and temporal scales in detail, whereas ground-based experiments can observe the ionospheric projections of transient magnetopause phenomena such as waves on the magnetopause driven by hot flow anomalies or flux transfer events produced by bursty reconnection across their full longitudinal and latitudinal extent. Magnetohydrodynamics (MHD), hybrid, and particle-in-cell (PIC) simulations are powerful tools to simulate the dayside transient phenomena. This paper provides a comprehensive review of the present understanding of dayside transient phenomena at Earth and other planets, their geoeffects, and outstanding questions.

Bow shock↗

DOE OSTI2022

Novel fabrication tools for dynamic compression targets with engineered voids using photolithography methods

Mesoscale imperfections, such as pores and voids, can strongly modify the properties and the mechanical response of materials under extreme conditions. Tracking the material response and microstructure evolution during void collapse is crucial for understanding its performance. In particular, imperfections in the ablator materials, such as voids, can limit the efficiency of the fusion reaction and ultimately hinder ignition. To characterize how voids influence the response of materials during dynamic loading and seed hydrodynamic instabilities, we, in this paper, have developed a tailored fabrication procedure for designer targets with voids at specific locations. Our procedure uses SU-8 as a proxy for the ablator materials and hollow silica microspheres as a proxy for voids and pores. By using photolithography to design the targets’ geometry, we demonstrate precise and highly reproducible placement of a single void within the sample, which is key for a detailed understanding of its behavior under shock compression. This fabrication technique will benefit high-repetition rate experiments at x-ray and laser facilities. Insight from shock compression experiments will provide benchmarks for the next generation of microphysics modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗

NASA NTRS1972

Computation of space shuttle flow fields using noncentered finite-difference schemes.

Second- and third-order, noncentered finite-difference schemes are described for the numerical solution of the hyperbolic equations of fluid dynamics. The advantages of noncentered methods over the more conventional centered schemes are: simpler programming logic, nonhomogeneous terms are easily included, and generalization to multidimensional problems is direct. Second- and third-order methods are compared with regard to dissipative and dispersive errors and shock-capturing ability. These schemes are then used in a shock-capturing technique to determine the inviscid, supersonic flow field surrounding space shuttle vehicles (SSV). Resulting flow fields about typical pointed and blunted, delta-winged SSVs at angle of attack are presented and compared with experiment.

Kutler, P.↗