Search NASA⌕ Search

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.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Orbital Modulation of Gamma Rays from PSR J2339-0533

In this work, we report on orbital modulation of the 100–600 MeV gamma-ray emission of the P B = 4.6 hr millisecond pulsar binary PSR J2339-0533 using 11 yr of Fermi Large Area Telescope data. The modulation has high significance (chance probability p ≈ 10 -7 ), is approximately sinusoidal, peaks near pulsar superior conjunction, and is detected only in the low-energy 100–600 MeV band. The modulation is confined to the on-pulse interval, suggesting that the variation is in the 2.9 ms pulsed signal itself. Additionally, this contrasts with the few other known systems exhibiting GeV orbital modulations, as these are unpulsed and generally associated with beamed emission from an intrabinary shock. The origin of the modulated pulsed signal is not yet clear, although we describe several scenarios, including Compton upscattering of photons from the heated companion. This would require high coherence in the striped pulsar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗
DOE OSTI2020

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2025

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗
DOE OSTI2021

Design of a 2 MW ZrC/W-based molten-salt-to-sCO 2 PCHE for concentrated solar power

To increase the power cycle efficiency and lower the levelized cost of electricity (LCOE) of concentrated solar power (CSP) plants, printed circuit heat exchangers (PCHEs) capable of operating above 700 °C with molten chloride salt and a sCO 2 -based fluid are needed. In this paper, the design of a high-pressure, high-temperature, 2 MW PCHE comprised of a thermomechanically-robust, zirconium carbide/tungsten (ZrC/W) composite is conducted for CSP plants. The ZrC/W composite is a material with high thermal conductivity, high stiffness, and high failure strength at high temperatures, along with excellent resistance to thermal cycling and thermal shock. In this work, a thermomechanical design analysis was conducted to select appropriate material thicknesses of the ZrC/W plates, and to determine the geometrical dimensions and the thermal performance of the PCHE. The influences of the plate number and heat exchanger length on power density and pressure drop have also been systematically investigated. Economic analyses were conducted to compare the cost of ZrC/W-based PCHEs to those comprised of IN740H (a state-of-the-art, nickel-based superalloy) and 316 stainless steel (316SS). Finally, at a sufficiently high plate production rate, the manufacturing cost of ZrC/W-based PCHEs can be significantly lower, while achieving a much higher power density, compared with state-of-the-art, nickel alloy-based and stainless steel-based PCHEs.

14 SOLAR ENERGY↗
DOE OSTI2021

High pressure, high flow rate batch mixing apparatus for high throughput experiments

An automated, high pressure, high flow rate batch mixing apparatus has been designed and constructed for rapid, stable, and repeatable mixing of multiple gases and vapors. The apparatus operates as an intermittent batch mixer with cycles of topping off fresh mixture to maintain pressure in an accumulator tank until consumed in an experimental apparatus. At high duty cycles, the apparatus can also function at steady state. This style of mixing is suitable for experiments such as high repetition rate shock tubes and other devices with intermittent flow demands. It is compact and portable, facilitating use in locations such as synchrotron light sources. The entire apparatus is heated to permit the mixing of vapors from species with low volatilities. The apparatus is fully automated and runs for extended periods with the only intervention being to refresh reagent supplies. Here, the accuracy and repeatability of the apparatus were verified by periodic gas sampling and analysis with gas chromatography. Multi-component mixtures spanning a wide range of complexity, dilution, and volatility of constituents have been prepared. The compositions of the majority of the mixture were found to be stable over several filling cycles, repeatable, and with the proper calibration of set-point conditions, accurate. Challenges were encountered preparing a mixture from multi-component liquids, and potential solutions are discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

X-ray Pump-Probe Measurements using a Laser-Plasma Accelerator (Final Report)

This project funded one graduate student, Mario Balcazar, to perform experiments on high power laser facilities to study the application of X-rays generated by laser wakefield accelerated electrons for advanced radiography, in particular X-ray phase contrast imaging of hydrodynamic shocks/instabilities. The work made use of the DOE’s LaserNetUS facilities. We performed dynamic phase-contrast X-ray imaging of the shock generated by a 200 ps duration laser pulse (E = 1 J, I = 10 15 Wcm –2 ) with a 30 µm diameter liquid target, with unprecedented spatio-temporal resolution. This includes multi-shock generation within the liquid jet and plasma instability formation. Innovative electron-beam radiography was used to probe the laser-plasma interplay finding evidence of bilateral heating of the water followed by strong electric field generation. These measurements help explain some of the discrepancies between simulation and experiment and pave the way to better plasma diagnostic systems in HED and ICF physics experiments. This work was performed in collaboration with researchers from Michigan, LLNL, LBL and Sandia National Laboratories, Queen’s University Belfast and Imperial College.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Double sided measurements and associated error bars for planar samples or subassemblies

Many of the high energy density (HED) experiments use targets with planar components usually in a multilayer stack as the physics package which is at the heart of the experiment. The physics package contains the sample part under study as well as other components as shown in Fig. 1, that when exposed to the laser drive applies pressure onto or shock the main component to take it to the desired state in the phase diagram. These other components may include the ablator, pusher, x-ray pre-heat shield, window or other special components. The components are usually glued together. Currently, experiments on NIF and Omega in the Dynamic Materials campaign use several platforms such as Diffraction, EOS and Strength (Ref 1-3) which all use such planar stack physics packages.

42 ENGINEERING↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proto-neutron star convection and the neutrino-driven wind: implications for the r-process

The neutrino-driven wind from proto-neutron stars is a proposed site for r-process nucleosynthesis, although most previous work has found that a wind heated only by neutrinos cannot produce the third r-process peak. However, several groups have noted that introducing a secondary heating source within the wind can change the hydrodynamic conditions sufficiently for a strong r-process to proceed. One possible secondary heating source is gravito-acoustic waves, generated by convection inside the proto-neutron star. As these waves propagate into the wind, they can both accelerate the wind and shock and deposit energy into the wind. Additionally, the acceleration of the wind by these waves can reduce the total number of neutrino captures and thereby reduce the final electron fraction of the wind. In neutron rich conditions, all of these effects can make conditions more favourable for r-process nucleosynthesis. Here, we present a systematic investigation of the impact of these convection-generated gravito-acoustic waves within the wind on potential nucleosynthesis. We find that wave effects in the wind can generate conditions favourable for a strong r-process, even when the energy flux in the waves is a factor of 10 –4 smaller than the total neutrino energy flux and the wind is marginally neutron rich. Nevertheless, this depends strongly on the radius at which the waves become non-linear and form shocks. We also find that both entropy production after shock formation and the acceleration of the wind due to stresses produced by the waves prior to shock formation impact the structure and nucleosynthesis of these winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

CISME report on experiments with HMX powder, PBX 9501 prills, and heated PBX 9501

The ability to predict post-ignition explosive response is relevant to explosives safety, particularly weapons assembly/disassembly operations that must operate under a sufficiently conservative safety basis to meet risk mitigation criteria for inadvertent nuclear detonation (IND). In particular, hazard determinations hinge on whether the PBX 9501 main charge can ultimately attain detonation from an initial, non-shock, ignition stimulus. The CISME experiment uses a spherical explosive charge, which is thermally ignited in the center using a fiber-delivered infrared laser pulse. This ignition mechanism - i.e. a localized thermal hotspot - is representative of the ignition-causing insults typical in accident scenarios, including low-velocity impact, drop, skid, and punch events. In fiscal year 2018 CISME experiments were conducted with pristine PBX 9501 spheres, confined and unconfined. Sub-detonative, low-violence reactions were observed. Here we report experiments conducted since 2018. The explosive charge was varied: experiments were conducted with room temperature pour-density HMX powder and prills of PBX 9501, and heated, thermally damaged PBX 9501. Experiments were conducted in a confined configuration, with two 4 mm-thick stainless-steel hemispherical shells bolted around the sphere. Detonation, via DDT (deflagration-to-detonation transition) was observed in the 3-inch diameter powdered HMX and the 6-inch diameter heated and thermally-damaged PBX 9501 tests. No transition to DDT was observed in smaller diameters for those configurations, and no DDT was observed for pour-density prills of PBX 9501 of any diameter. The original test series demonstrated that the response of pristine, undamaged PBX 9501 to a central ignition is relatively benign. This new test series shows that increased porosity, lack of polymer binders and thermal damage can greatly increase the violence of the post-ignition response, up to and including detonation. These results bridge the explosive response of PBX 9501 from sub-detonative through detonation as a function of diameter and provide valuable test cases for modeling efforts that strive to capture the DDT threshold. By exploring explosive response in a more representative 3D geometry, we can also benchmark these new results to more plentiful DDT datasets generated from testing in 1D tube confinement (the historical test configuration)–the configuration upon which almost all DDT models were developed.

42 ENGINEERING↗
DOE OSTI2020

IHS-60 Stored Energy Analysis

The IHS-60 consists of a General Purpose Heat Source (GPHS) Fueled Clad (FC) that has been hot isostatically pressed to eliminate the gap between the iridium clad and the fuel pellet, thus preventing differential motion between the pellet and the clad in shock and vibration environments. The vented FC is enclosed within a Molybdenum barrier comprising the liner body and the lower support member. The slip- t Molybdenum liner members serve as a mechanical barrier to plutonia contact of the tantalum alloy (T-111) strength member in the event the iridium clad should fail. The slip- t liner is not a barrier to gaseous transport, so a tantalum getter is installed near the vent region of the FC to preclude potential oxygen transport to the T-111 interior surface. The T-111 strength member encloses the Molybdenum liner assembly and provides for both helium gas retention and isotope fuel containment under impact and external pressure environments. It is surrounded by an outer clad of Hastelloy-S which protects the strength member from external oxidation and corrosive sources during nominal operation and in the event of off-nominal operation or accidents. Both the strength member and outer clad parts are closed with seal welds. The strength member and outer clad are the primary containment members of the 60-Watt IHS. The iridium clad contains the fuel pellet during nominal operation and mild accident environments, but it is not considered to be a primary containment member.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗
DOE OSTI2024

High-pressure phase transition of olivine-type Mg 2 GeO 4 to a metastable forsterite-III type structure and their equations of state

Germanates are often used as structural analogs of planetary silicates. We have explored the high-pressure phase relations in Mg 2 GeO 4 using diamond-anvil cell experiments combined with synchrotron X-ray diffraction and computations based on density functional theory. Upon room temperature compression, forsterite-type Mg 2 GeO 4 remains stable up to 30 GPa. At higher pressures, a phase transition to a forsterite-III type (Cmc2 1 ) structure was observed, which remained stable to the peak pressure of 105 GPa. Using a third-order Birch Murnaghan fit to the experimental data, we obtained V 0 = 305.1(3) Å3, K 0 = 124.6(14) GPa, and $K'_0$ = 3.86 (fixed) for forsterite-type Mg 2 GeO 4 and V 0 = 263.5(15) Å 3 , K 0 = 175(7) GPa, and $K'_0$ = 4.2 (fixed) for the forsterite-III type phase. The forsterite-III type structure was found to be metastable when compared to the stable assemblage of perovskite/post-perovskite + MgO, as observed during laser-heating experiments. Understanding the phase relations and physical properties of metastable phases is crucial for studying the mineralogy of impact sites, understanding metastable wedges in subducting slabs, and interpreting the results of shock compression experiments.

58 GEOSCIENCES↗
DOE OSTI2024

Multiphysics Simulations of Thermal Shock Testing of Nanofibrous High Power Targets

Increase of primary beam power for neutrino beam-lines leads to a reduced lifespan for production targets. The field of High Power Targetry (HPT) is generating new concepts to meet the need for robust targets. One idea being investigated by the HPT Research and Development Group at Fermilab is an electrospun nanofiber target. As part of their evaluation, samples with different densities were sent to the HiRadMat facility at CERN for thermal shock tests. The samples with the higher density, irradiated under a high intensity beam pulse, exhibit major damages at the impact site whereas those with the lower density show no apparent damages. The exact cause of this failure was unclear at the time. In this paper, we present the results of multiphysics simulations of the thermal shock experienced by the nanofiber targets that suggest the failure originates from the reduced permeability of the high density sample to airflow. The air present in the porous target expands due to heating from the beam, but is unable to flow freely in the high density sample, resulting in a larger back pressure that blows apart the mat. We close with a discussion on how to further validate this hypothesis.

43 PARTICLE ACCELERATORS↗