SEARCH · Engineering Papers

Results for “shock hydrodynamics”

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

DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗

DOE OSTI2020

Calibration of reactive burn and Jones-Wilkins-Lee parameters for simulations of a detonation-driven flow experiment with uncertainty quantification

Here, uncertainties in the explosive-specific parameters of the Jones-Wilkins-Lee (JWL) equation of state (EOS) are carefully considered in hydrodynamic simulations of an explosive experiment to minimize the error in the flow prediction. Experimental data of the leading shock position in the transverse direction over time serves as the prediction metric for quantifying simulation prediction error. The uncertainty quantification technique, global sensitivity analysis, is utilized to determine the JWL parameters to which the transverse shock propagation is most sensitive. A polynomial response surface (PRS) is constructed in the space of the most influential JWL parameters, and the point of minimum error between the experimental data and the PRS yields calibrated JWL parameters for the experimental flow. The simulation results following the parameter calibration show good agreement with the experimental data. It was found that two significant parameters, the heat release per unit mass of reactant Q and JWL model exponent R 1 are strongly related, which makes it difficult to identify accurate values.

36 MATERIALS SCIENCE↗

NASA NTRS1980

Theoretical interpretation of Traveling Interplanetary Phenomena and their solar origins

A review is made of recent theoretical studies on Traveling Interplanetary Phenomena (TIP) and their possible solar origins in the context of mathematical methods and physical processes. Flare-generated interplanetary shocks are used as examples to illustrate the following set of contrasting approaches to the phenomena: analytical vs numerical, hydrodynamic vs magnetohydrodynamic, dissipative vs non-dissipative, continuum/macroscopic vs kinetic/microscopic. It is suggested in conclusion that future research should emphasize continuum/macroscopic physical processes together with the kinetic approach to deepen understanding of microscopic interactions.

Wu, S. T.↗

NASA NTRS1982

Another self-similar blast wave: Early time asymptote with shock heated electrons and high thermal conductivity

Accurate approximations are presented for the self-similar structures of nonradiating blast waves with adiabatic ions, isothermal electrons, and equation ion and electron temperatures at the shock. The cases considered evolve in cavities with power law ambient densities (including the uniform density case) and have negligible external pressure. The results provide the early time asymptote for systems with shock heating of electrons and strong thermal conduction. In addition, they provide analytical results against which two fluid numerical hydrodynamic codes can be checked.

Cox, D. P.↗

NASA NTRS1988

Cosmic-ray-modified stellar winds. II - A perturbation approach

A perturbation method is developed to describe the modification of a stellar wind flow with termination shock by the Galactic cosmic rays. The perturbation parameter epsilon used in the analysis is the ratio of the galactic cosmic-ray pressure to the thermal gas pressure at a large distance from the star. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, particularly just upstream of the termination shock of the wind. A one-fluid polytropic model is used to describe the thermal stellar wind gas, with the coupling between the cosmic rays and the thermal gas being determined by an average hydrodynamical diffusion coefficient kappa describing the scattering of cosmic rays by magnetohydrodynamic waves traveling in the background flow. The analysis takes into account cosmic-ray modifications of the critical point of the wind, the thermal gas entropy constants on both sides of the shock, and the fluid velocity profile.

Ko, C. M.↗

NASA NTRS1997

A comparative modeling of supernova 1993J

The light curve of Supernova 1993J is calculated using two computational radiation transport approaches. The two approaches are represented by the computer codes STELLA and EDDINGTON. The emphasis is on the shock breakout and the photometry in the U, B and V bands during the first 120 days of the supernova. The STELLA model includes implicit hydrodynamics and is able to model early supernova evolution before the expansion is homologous. The STELLA model employs multi-group photonics and is able to follow the radiation as it decouples from the matter. The EDDINGTON code uses an algorithm for integrating the transport equation which assumes homologous expansion and uses a finer frequency resolution. The agreement between the two codes is considered to be satisfactory only in the case where compatible physical assumptions are made concerning the opacity. The assumptions are justified. The continuum spectrum for SN 1993J is predicted near the shock breakout to be superior to that predicted by standard single energy group hydrocodes. The uncertainties involved in current time dependent models of supernova light curves are discussed.

Blinnikov, Sergei↗

DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗

NASA NTRS1986

Studies on counterstreaming plasma expansion

Recent studies on counterstreaming plasma expansions are summarized. The basic phenomenon of plasma expansion is reviewed, and results from one-dimensional simulations of counterstreaming plasma expansion are discussed. Results from simulations based on an electrostatic particle-in-cell code, in which the dynamics of both the electrons and ions are exactly followed, are discussed. The formation of electrostatic shocks is addressed. Finally, results are presented on the ionospheric plasma expansion along the geomagnetic flux tubes by solving the hydrodynamic equations.

Singh, N.↗

DOE OSTI2022

Production efficiencies of sound waves in the intracluster medium driven by AGN jets

Feedback from active galactic nuclei (AGNs) is believed to be the most promising solution to the cooling flow problem in cool-core clusters. Dissipation of sound waves is considered as one of the possible heating mechanisms; however, its relative contribution to heating remains unclear. To estimate the energy budget for heating, we perform 3D hydrodynamic simulations of AGN jet injections in a Perseus-like cluster and quantify the amount of energy stored in the forms of weak shocks and waves. For this work, e find that, for a single jet injection with typical parameters in cool-core clusters, ${\sim}9{{\ \rm per\ cent}}$ of the total jet energy is stored in compressional waves (including both shocks and waves). However, due to the destructive effects among randomly phased waves as well as the dissipation of shock energies, in our simulations including self-regulated AGN feedback, no more than $3{{\ \rm per\ cent}}$ of the total injected energy goes into compressional waves. We further separate the energy contribution from shocks and waves and find that, for a single outburst, the shocks can only contribute to ${\sim}20{-}30{{\ \rm per\ cent}}$ of the total compressional energy in the inner radii and quickly dissipate away. In the self-regulated case where shocks are repeatedly generated, shocks completely dominate over sound waves in the inner region and can still provide ${\sim}40{-}50{{\ \rm per\ cent}}$ of the total compressional energy at outer radii. Our results suggest that the production of sound waves is not as efficient as what was previously found, and thus sound wave dissipation may be a subdominant source of heating in cool-core clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

Here, the Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2015

Late-Time Evolution of Composite Supernova Remnants: Deep Chandra Observations and Hydrodynamical Modeling of a Crushed Pulsar Wind Nebula in SNR G327.1-1.1

In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 ks Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology; a cometary structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence for mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock(RS) that can occur as a result of a density gradient in the ambient medium and or a moving pulsar that displaces the PWN from the center of the remnant. We present hydrodynamical simulations of G327.1-1.1 that show that its morphology and evolution can be described by a approx. 17,000 yr old composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar's motion. We also show that the RSPWN interaction scenario can reproduce the broadband spectrum of the PWN from radio to gamma-ray wavelengths. The analysis and modeling presented in this work have important implications for our general understanding of the structure and evolution of composite SNRs.

Evolution↗

DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗

DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

DOE OSTI2021

The effect of jet–ejecta interaction on the viewing angle dependence of kilonova light curves

ABSTRACT The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet–ejecta interaction as initial conditions for multidimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at $1.5\, \mathrm{h}$ after merger. We find that on this time-scale, jet shock heating does not affect the kilonova emission for the jet parameters we survey. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near (≲30°) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

Klion, Hannah↗

NASA NTRS1994

On the X-ray spectrum of Kepler's supernova remnant

We have devised a method to do nonequilibrium ionization calculations on the results of two-dimensional hydrodynamical simulations, based on the algorithm of Hughes & Helfand (1985). We have calculated the ionization structure and X-ray emission for a two-dimensional numerical hydrodynamical simulation for the remnant of Kepler's supernova (SN); the hydrodynamical model was presented in a previous paper. In this model, the progenitor of Kepler's SN is assumed to have been a massive runaway star ejected from the Galactic plane. In its red supergiant stage, its dense stellar wind was distorted and compressed into a bow shock by the ram pressure of the tenuous interstellar medium. The subsequent interaction of the supernova ejecta with this asymmetric circumstellar matter produced a strongly asymmetric supernova remnant (SNR). In this paper, we present calculated X-ray spectra for this hydrodynamical model. A comparison with observations implies a moderate overabundance of Fe in Kepler's SNR (only 50% larger than its cosmic value), in contrast to a large (6 to 15) Fe overabundance derived previously. However, we confirm earlier conclusions that Si and S abundances are 2 to 3 times solar. These modest enhancements of Si, S, and Fe may be attributed either to heavy-element enriched SN ejecta or to the initial chemical abundances of the SN progenitor, which originated in the metal-rich inner Galaxy. The comparison of our models with the observed spectra confirm theoretical predictions that moderate electron heating occurs at strong collisionless shock fronts, with the implied electron/mean temperature ratio of approximately 0.5.

Borkowski, Kazimierz J.↗

NASA NTRS2011

A Chandra X-Ray Survey of Ejecta in the Cassiopeia A Supernova Remnant

We present a survey of the X-ray emitting ejecta in the Cassiopeia A supernova remnant based on an extensive analysis of over 6000 spectral regions extracted on 2.5-10" angular scales using the Chandra 1 Ms observation. We interpret these results in the context of hydrodynamical models for the evolution of the remnant. The distributions of fitted temperature and ionization age are highly peaked and suggest that the ejecta were subjected to multiple secondary shocks. Based on the fitted emission measure and element abundances, and an estimate of the emitting volume, we derive masses for the X-ray emitting ejecta as well as showing the distribution of the mass of various elements over the remnant. The total shocked Fe mass appears to be roughly 0.14 Solar Mass, which accounts for nearly all of the mass expected in Fe ejecta. We find two populations of Fe ejecta, that associated with normal Si-burning and that associated with alpha-rich freeze-out, with a mass ratio of approximately 2:1. Surprisingly, essentially all of this Fe (both components) is well outside the central regions of the SNR, presumably having been ejected by hydrodynamic instabilities during the explosion. We discuss this, and its implications for the neutron star kick.

Hwang, Una↗

DOE OSTI2022

Experiments on the single-mode Richtmyer–Meshkov instability with reshock at high energy densities

Here, the hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗