SEARCH · Search NASA

Results for “SHOCK HEATING”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2020

An analytical expression for temperature in a thermodynamically consistent model with a Mie-Gruneisen equation for pressure

The objective of this paper is to develop an analytical expression for temperature in a thermodynamically consistent model with a Mie-Gruneisen equation for pressure. This expression for the temperature is obtained for both constant and variable specific heat models and can be used in common models for the material yield strength that include thermal softening due to increase in temperature. Additionally, the material parameters in the model can easily be calibrated to match known experimental data and the model has been used to determine the value of the shock pressure at the onset of melting during shock loading.

42 ENGINEERING↗

DOE OSTI2024

Particle Image Velocimetry and Total Temperature Characterization of Arc-Heated Mach 4.5 Free Jet

Particle image velocimetry (PIV) measurements were conducted in the Arc-Heated Combustion Tunnel (ACT-II) on a high-enthalpy Mach 4.5 nozzle flow. The primary motivation for the current work is to establish the feasibility of PIV as a high-quality flow diagnostic technique in hypersonic arc-heated flows and produce a high-resolution hypersonic nozzle flow characterization, including nozzle free jet shear layer velocity profiles. Due to the uncertain nature of the gas composition in the flowfield of interest, PIV measurements of a supersonic wedge flow were used to couple the flow velocities to the oblique shock relations. This enabled a PIV-derived total temperature calculation that does not require a priori assumptions of the exact gas composition. A PIV expansion module was developed for this study, and proper seed sizing and material were selected for the present flow. An analysis of the seed particle response from the PIV measurements provided an initial validation of the seed performance while revealing potential issues with consistent seed sizing. These measurements and derived conditions represent a means of acquiring high-resolution velocimetry and temperature datasets that consider uncertainties in the experimental generation of hypersonic high-enthalpy flows.

Engineering↗

DOE OSTI2020

Ion-scale Current Structures in Short Large-amplitude Magnetic Structures

Magnetic reconnection is an explosive phenomenon widely happen in plasma environments. Recent studies show that reconnection could occur in the shock transition region. We study thin current layers in the region in front of the shock and find evidence of reconnection. The analysis also shows the possible generation mechanism of the thin current sheet that further enables reconnection. A large-scale wave is further generated, and the thin current sheets are generated through compression or secondary waves. The thin current structures also lead to electron heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Laboratory Resolved Structure of Supercritical Perpendicular Shocks

Supermagnetosonic perpendicular flows are magnetically driven by a large radius theta-pinch experiment. Fine spatial resolution and macroscopic coverage allow the full structure of the plasma-piston coupling to be resolved in laboratory experiment for the first time. Here, a moving ambipolar potential is observed to reflect unmagnetized ions to twice the piston speed. Magnetized electrons balance the radial potential via Hall currents and generate signature quadrupolar magnetic fields. Electron heating in the reflected ion foot is adiabatic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

DOE OSTI2022

Electron energization dynamics in interaction of self-generated magnetic vortices in upstream of collisionless electron/ion shocks

Abstract Relativistic collisionless shocks are considered responsible for particle energization mechanisms leading to particle acceleration. While electron energization in shock front region of electron/ion collisionless shocks are the most studied, the mechanism of electron energization in interaction with self-generated magnetic vortices (MVs) in the upstream region is still unclear. We investigate electron energization mechanism in the upstream region of electron/ion relativistic collisionless shocks, using two dimensional particle-in-cell (PIC) simulations. We discuss mechanism of electron energization which takes place in the upstream region of the shock, where the counter stream particles interact with incoming flow. The energy gain of electrons happens during their interaction with evolving fields of self-generated magnetic vortices in this region. Three Fermi-like electron energization scenarios are discussed. Stochastic acceleration of electrons in interaction with fields of MV leads to anisotropic heating of fast electrons due to diffusion in the momentum space of electrons and, finally, synergetic effect of evolving fields of MVs leads to the formation of a power-law tail of supra-thermal particles.

Naseri, N.↗

DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

FENIX: Towards a Fully Integrated Multiphysics Framework for Plasma Facing Component Modeling

Computational tools have a crucial role to play in accelerating the deployment of fusion as a clean, reliable, abundant, and sustainable energy source. Multiphysics, high-fidelity simulation capabilities can help model, study, and predict intricate interactions between materials performance, plasma exposure, neutron irradiation, and engineering processes. As such, they can assist in the resolution of scientific and engineering challenges underpinning design, construction, and commission of fusion power plants. To address these needs, ongoing efforts are leveraging the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework and delivering new computational tools for the fusion community. These tools inherit crucial attributes from MOOSE. They are open-source, modular, integrated with nuclear industry-standard software quality assurance processes, and enable multiphysics, multi-fidelity, fully integrated, zero- to three-dimensional, and massively parallel simulations. After a short overview of these capabilities, we will present the development of Fusion ENergy Integrated multiphys-X (FENIX), a MOOSE-based application designed to enable plasma facing component design and performance evaluation. Throughout their lifetime, plasma facing components are exposed to extreme thermal loads, repeated thermal shocks, and irradiation by plasma ions, neutral particles, and high-energy neutrons. Consequently, designing a plasma facing component with acceptable lifetime degradation is extremely challenging. FENIX aims to model the multiphysics environment in which plasma facing components evolve to accelerate their design studies. To that end, FENIX couples existing MOOSE capabilities such as heat transfer, thermomechanics, and thermal hydraulics, with tritium transport via the MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8), with neutronics via the MOOSE-based high-fidelity neutron-photon transport and fluid dynamics code Cardinal, and finally with Particle-in-Cell plasma simulation capabilities being developed in this project. In this study, we present the current FENIX capabilities and preliminary results of its application to model the Tritium Plasma Experiment set up at Idaho National Laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The complex intracluster medium of Abell 1569 and its interaction with central radio galaxies

ABSTRACT We present the first in-depth study of X-ray emission from a nearby (z ∼ 0.0784) galaxy cluster Abell 1569 using an archival Chandra observation. A1569 consists of two unbound subclusters – a northern subcluster (A1569N) hosting a double-lobed radio galaxy 1233+169 at its centre, and a southern subcluster (A1569S) harbouring a wide-angle-tailed (WAT) radio source 1233+168. X-ray emission from A1569N and A1569S extends to a radius r ∼248 kpc and r ∼370 kpc, respectively, indicating that the two gas clumps are group-scale systems. The two subclusters have low X-ray luminosities (∼1042–43 erg s−1), average elemental abundances ∼1/4 Z⊙, low average temperatures (∼2 keV), and lack large (r ≳ 40–50 kpc) cool cores associated with the intracluster gas. We detect a pair of cavities coincident with the radio lobes of 1233+169 in A1569N. The total mechanical power associated with the cavity pair is an order of magnitude larger than the X-ray radiative loss in the cavity-occupied region, providing corroborating evidence for cavity-induced heating of the intragroup gas in A1569N. A1569S exhibits possible evidence for a small-scale cluster-subcluster merger, as indicated by its high central entropy, and the presence of local gas elongation and a density discontinuity in between the bent radio tails of 1233+168. The discontinuity is indicative of a weak merger shock with Mach Number, M ∼ 1.7. The most plausible geometry for the ongoing interaction is a head-on merger occurring between A1569S and a subcluster falling in from the west along the line bisecting the WAT tails.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A Simple Hydride Model for Cerium Ejecta Particles

Cerium ejecta particles form after a shock wave impacts a cerium plate with a known roughness or prescribed surface perturbation. When these particles are ejected into hydrogen or deuterium gas the material reacts exothermically and raises the temperature of the particle relative to the gas. A simple model is developed, assuming that the particle remains intact and spherical, to capture the heat and mass transfer that occurs during this process. Model performance is evaluated by comparing with cerium ejecta experiments where the particles are ejected into deuterium gas at initial pressures of 4 and 8 atm. Overall, the model is able to capture the approximately 400 K increase in particle temperature above the surrounding gas temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A simple hydride model for cerium ejecta particles

36 MATERIALS SCIENCE↗

DOE OSTI2022

A simple hydride model for cerium ejecta particles

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3

In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Comparison of the Core-collapse Evolution of Two Nearly Equal-mass Progenitors

Abstract We compare the core-collapse evolution of a pair of 15.8 M ☉ stars with significantly different internal structures, a consequence of the bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79 M ☉ progenitors have core masses (masses interior to an entropy of 4 k B baryon −1 ) of 1.47 and 1.78 M ☉ and compactness parameters ξ 1.75 of 0.302 and 0.604, respectively. The core-collapse simulations are carried out in 2D to nearly 3 s postbounce and show substantial differences in the times of shock revival and explosion energies. The 15.78 M ☉ model begins exploding promptly at 120 ms postbounce when a strong density decrement at the Si–Si/O shell interface, not present in the 15.79 M ☉ progenitor, encounters the stalled shock. The 15.79 M ☉ model takes 100 ms longer to explode but ultimately produces a more powerful explosion. Both the larger mass accretion rate and the more massive core of the 15.79 M ☉ model during the first 0.8 s postbounce time result in larger ν e / ν ¯ e luminosities and RMS energies along with a flatter and higher-density heating region. The more-energetic explosion of the 15.79 M ☉ model resulted in the ejection of twice as much 56 Ni. Most of the ejecta in both models are moderately proton rich, though counterintuitively the highest electron fraction ( Y e = 0.61) ejecta in either model are in the less-energetic 15.78 M ☉ model, while the lowest electron fraction ( Y e = 0.45) ejecta in either model are in the 15.79 M ☉ model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Thermal-shock experiments for separate-effects validation of UO 2 fuel fracture models

Because of the important role that fracture plays in the behavior of ceramic UO2 fuel in a nuclear reactor environment, fracture models are a major component of fuel performance codes. As with any aspect of fuel performance, it is crucial to validate these fracture models against experimental data; however, obtaining well-controlled data for conditions representative of a reactor environment is difficult. Quenching is proposed here as a relatively simple approach for using a laboratory environment to achieve conditions that approximate those of a reactor environment. In this paper, an experimental apparatus containing a single instrumented fuel pellet in a sealed section of copper tubing is developed. It is then applied to a series of seven experiments in which the apparatus is first heated to a high temperature (580-680°C) by immersion in a molten salt bath, then quenched in a cold bath (-10-4 °C). Development of these experiments was guided by numerical simulations, and post-test simulations were performed to aid in understanding the experimental behavior and assessing the accuracy of the predictions of fracture initiation and propagation. In addition, similar experiments were performed on solid copper rods to provide temperature-dependent heat transfer coefficients for use in simulations of these experiments. Here, this study demonstrates that quenching is a viable approach for generating thermal gradients representative of those in prototypical light-water reactor conditions at powers of about 5–10 kW/m. Fracture is expected to begin at these power levels, and moderate amounts of radial and axial cracking was observed in these quenching tests.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗