Engineering Papers⌕ Search

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.

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

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1996

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗
DOE OSTI2020

Bipolar planetary nebulae from outflow collimation by common envelope evolution

ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE winds may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they become core-collapse supernovae (SNe). The resulting dense circumstellar medium may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomenon were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward-moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2–4 km s -1 convective motions seed the intrinsic convective instability associated with the high-entropy material deep in the envelope, enabling gas from deep within the envelope to escape and increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by ≈10× along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Reynolds stress closure in jet flows using wave models

A collection of papers is presented. The outline of this report is as follows. Chapter three contains a description of a weakly nonlinear turbulence model that was developed. An essential part of the application of such a closure scheme to general geometry jets is the solution of the local hydrodynamic stability equation for a given jet cross-section. Chapter four describes the conformal mapping schemes used to map such geometries onto a simple computational domain. Chapter five describes a solution of a stability problem for circular, elliptic, and rectangular geometries. In chapter six linear models for the shock shell structure in non-circular jets is given. The appendices contain reprints of papers also published during this study including the following topics: (1) instability of elliptic jets; (2) a technique for predicting the shock cell structure in non-circular jets using a vortex sheet model; and (3) the resonant interaction between twin supersonic jets.

Morris, Philip J.↗
NASA NTRS1980

Physicochemical and dynamical processes in cometary ionospheres. I - The basic flow profile

The paper deals with the basic aspects of solar wind interaction with the ionosphere of a large comet and with detailed hydrodynamic models developed to describe the flow exterior to the tangential discontinuity surface that separates the outflowing purely cometary ions from the inflowing contaminated solar wind. The equation of the basic flow model and the results of a preliminary analysis of these equations are given, and the shape and structure of the shock layer are analyzed. The shock layer structure associated with fan-shaped comet emanating from the cometary nucleus in the general direction to the sun is identified.

Houpis, H. L. F.↗
DOE OSTI2020

Modeling for pRad Explosive Experiments

The hydrodynamics code CTH was used to simulate high-explosive detonations in various geometries to design a planar, high-pressure drive with Taylor wave release for experiments in the proton radiography containment vessel. The most successful geometry is based on a P25 lens with some of its outer diameter and output end removed. Shock planarity was significantly improved relative to an unmodified P25 lens for this case, where its full diameter was not needed, and the calculated timing was adequate to trigger the experiments successfully.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Herbig-Haro objects

The IUE (International Ultraviolet Explorer) spectra of Herbig-Haro (HH) objects and their relation to the very detailed optical spectra available for these objects are studied. Useful information about the physics of HH objects in general and especially about their hydrodynamic models which are contained in these IUE observations is outlined. The merged spectra of high excitation, low excitation HH objects (like HH43 and HH47), and peculiar HH objects are discussed. Results of the spatial variation of lines and continua in HH objects are discussed. Their compatibility with predictions for aerodynamic and especially bow shock models is discussed. The problems arising in the interpretation of the continuous energy distribution in the short wavelength range are discussed. They require the presence of other emission mechanisms in addition to the collisionally enhanced two photon continuum. In HH43 the fluorescent H2 line emission comes from a surprisingly small spatial region. Implications of these results are discussed.

Boehm, K. H.↗
NASA NTRSDate not supplied

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗
DOE OSTI2025

Learning robust parameter inference and density reconstruction in flyer plate impact experiments

Estimating physical parameters or material properties from experimental observations is a common objective in many areas of physics and material science. In many experiments, especially in shock physics, radiography is the primary means of observing the system of interest. However, radiography does not provide direct access to key state variables, such as density, which prevents the application of traditional parameter estimation approaches. Here we focus on flyer plate impact experiments on porous materials, and resolving the underlying parameterized equation of state (EoS) and crush porosity model parameters given radiographic observation(s). We use machine learning as a tool to demonstrate with high confidence that using only high impact velocity data does not provide sufficient information to accurately infer both EoS and crush model parameters, even with fully resolved density fields or a dynamic sequence of images. We thus propose an observable data set consisting of low and high impact velocity experiments/simulations that capture different regimes of compaction and shock propagation, and proceed to introduce a generative machine learning approach which produces a posterior distribution of physical parameters directly from radiographs. We demonstrate the effectiveness of the approach in estimating parameters from simulated flyer plate impact experiments, and show that the obtained estimates of EoS and crush model parameters can then be used in hydrodynamic simulations to obtain accurate and physically admissible density reconstructions. Finally, we examine the robustness of the approach to model mismatches, and find that the learned approach can provide useful parameter estimates in the presence of out-of-distribution radiographic noise and previously unseen physics, thereby promoting a potential breakthrough in estimating material properties from experimental radiographic images.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

The Advanced Simulation and Computing (ASC) Transport project’s Jayenne Implicit Monte Carlo (IMC) transport library now includes an unstructured mesh capability and is available in a Lagrangian Applications Project (LAP) code. In this presentation, we discuss recent work by the LAP and Transport projects that provides an IMC transport capability for radiation hydrodynamics in the Lagrangian frame. Verification problems testing the new capabilities have been simulated and analyzed, i.e. Marshak wave, Mach 45, Su-Olsen, picket fence, and crooked pipe, both in one and two dimensions. We also present results on two stretch goal problems: a simplified COAX high energy density physics experiment and a supernova shock simulation. Finally, we identify current limitations and future work needed to bring a full capability to the user community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Explorer 33 and 35 plasma observations of magnetosheath flow.

Mapping of the magnetosheath flow in the region -60 earth radii less than or equal to X sub SE less than or equal to -20 earth radii by using plasma data from the MIT plasma experiments on the earth-orbiting Explorer 33 and the lunar-orbiting Explorer 35. Temporal and spatial variations in the flow are separated by using 309 hours of simultaneous data from the two spacecraft, one spacecraft being in the solar wind and the other spacecraft in the magnetosheath. The overall flow pattern agrees well with the predictions of a simplified hydrodynamic theory. The measured ratio of magnetosheath to solar-wind densities is shown to have some dependence on the observed solar-wind sonic Mach number. This density ratio is consistently less than 1 in the region adjacent to the magnetotail boundary and thus reflects the irreversible nature of the flow across the bow shock.

Howe, H. C., Jr.↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2021

Post-explosion Evolution of Core-collapse Supernovae

Here we investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 explosion models during several seconds. We present a broad study based on three progenitors (11.2, 15, and 27 M ⊙ ), different neutrino heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind (NDW) that survives for several seconds. This indicates that NDWs are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 M ⊙ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

The growth of thermal instabilities in cooling flows

Numerical hydrodynamic simulations are used to determine the evolution of the condensation modes for thermal instabilities in 10 to the 6th to 10 to the 8th K gas. These simulations properly take into account the effects of nonequilibrium cooling below 10 to the 6th K. It is found that only perturbations greater than about 3 kpc develop shocks in the hot environment of cooling flows. These large perturbations evolve through three distinct phases: (1) the formation of a cold core, (2) a supersonic accretion phase onto the cold core which lasts for 10 to the 5th to 10 to the 6th yr with peak shock velocities about 120 km/s, and (3) a subsonic accretion phase. The X-ray emissivity and the optical-ultraviolet line emission associated with thermal instabilities are calculated. These calculations demonstrate that the optical emission from cooling flows probably does not arise from shocked condensations. The calculated optical emission lines are of a higher ionization level and a lower luminosity than the observations. Ultraviolet lines produced in the shocks may be detectable. The core must cool significantly after the shock has dissipated before star formation can commence.

David, Laurence P.↗
DOE OSTI2022

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Go with the Radflow: Thermal radiation experiments in the high energy-density regime [Slides]

Los Alamos National Laboratory is a nuclear weapons laboratory supporting our nation's defense. In support of the mission is a high energy-density physics program in which we design and execute experiments to study radiation-hydrodynamics phenomena and improve the predictive capability of our large-scale multi-physics software codes on our supercomputers. The Radflow Project is maturing a unique spectroscopic measurement that yields a spatially dependent supersonic radiation wavefront profile. The spectroscopic measurement is much more constraining than the typically used radiography, which measures the density changes due to the shock as the radiation wave goes subsonic. Recently, the spectroscopic measurement has been applied to advanced targets with a detailed heterogeneity representing a single realization of a stochastic medium. The Radflow Project also conducts opacity experiments to explain the current conundrum of the experimental iron opacity differing from theory, which is important in modeling the sun. More than that, though, we need to be able to answer these science questions and be able to validate our theory and codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗