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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1970

Galactic astronomy; Proceedings of the Second Summer Institute for Astronomy and Astrophysics, State University of New York, Stony Brook, N.Y., June 19-July 17, 1968. Volumes 1 & 2.

Basic problems on the structure of our galaxy are presented. The essentially spiral structure of the galaxies is examined by a theory of the spiral shape, and by a study of the spiral structure with the 21-cm line of neutral hydrogen. Irregular spirals are discussed, and some problems that arise as a result of large deviations from the regular spiral structure are treated. The density wave theory of spiral structure is compared with the observed spiral structure of the galaxies. A section on the star migration technique, which is used to test the density wave theory of galactic spirals, is included. The interstellar medium is described in terms of such phenomena as supernova explosions and cosmic rays, both of which contribute to the energy balance of the medium. Other topics that are germaine to a description of the interstellar medium are the UV spectrum, kinematics of interstellar clouds, and ionization equilibrium. The implications of spiral galactic shock pattern on star formation along spiral arms are investigated, and by means of a blast wave model, it is shown that a hydrodynamic shock of galactic scale exists in galaxy M82.

Chiu, H.-Y.↗
DOE OSTI2024

A Shock Stabilization of the HLLC Riemann Solver for the Carbuncle Instability

The HLLC approximate Riemann solver improves upon the HLL Riemann solver by resolving contact discontinuities. This is a particularly desirable property for multi-material codes in which problems usually contain material interfaces. However, the HLLC solver is known to suffer from the carbuncle phenomenon, a numerical instability most apparent at grid-aligned shocks in multi-dimensional simulations. Many problems of interest, including high energy-density physics applications, require the accurate resolution of both material interfaces and hydrodynamic shocks. A variety of methods have been developed to cure this instability, with varying degrees of complexity. The objective of this work is to describe a simple approach to modify the HLLC Riemann solver and prevent the carbuncle instability. Here, the method is then demonstrated for assorted two-dimensional test problems known to exhibit the shock instability. The performance of the new solver is compared with that of the standard HLL and HLLC Riemann solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Solid-to-plasma transition of polystyrene induced by a nanosecond laser pulse within the context of inertial confinement fusion

Laser Direct Drive (LDD) Inertial Confinement Fusion (ICF) involves irradiating a spherical target of thermonuclear fuel coated with an ablator, usually made of polystyrene. Laser energy absorption near the target surface leads to matter ablation, hydrodynamic shocks, and ultimately capsule implosion. The conservation of spherical symmetry is crucial for implosion efficiency, yet spatial modulations in laser intensity can induce nonuniformities, causing the laser imprint phenomenon. Understanding laser imprint, especially considering the initial solid state, is essential for advancing LDD ICF. A first microscopic model of solid-to-plasma transition was built in 2019, accounting for laser absorption in the solid state with a band structure based ionization model. This model has been improved to include chemical fragmentation and a more accurate description of electron collision frequency in various matter states. The latest development involves assessing the model reliability by comparing theoretical predictions with experimental observations. Despite the success of this approach, questions remain, leading to further investigations and observations under different irradiation conditions. This work presents an experiment with a nanosecond pulse, taking into account hydrodynamic effects, and measures transmission dynamics over the entire laser beam area to observe two-dimensional effects. As a result, the objective is to adapt the theoretical model, couple it with a hydrodynamic code, and observe additional effects related to the initial solid state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Far-infrared emission lines of CO and OH in the Orion-KL molecular shock

Observations of far infrared rotational emission lines which arise in the shocked gas associated with Orion-Kl are presented, including detections of the CO J = 34 yields 33, J = 31 yields 30, J = 26 yields 25, and OH sup 2 PI sub (3/2) J sup P = 7/2(-) yields 5/2(+) emission lines, as well as improved measurements of the CO J = 22 yields 21 and OH sup 2 PI sub (3/2) J = 5/2 yields 3/2 lines. These lines are observed to have velocity widths of Del V approx. 20 to 30 km/sec, somewhat less than either the 2 micro H sub 2 lines or the high velocity plateau component of the millimeter wave CO lines seen in this object. An H sub 2 column density of approx. 3 x 10 to the 21st power, a total mass of approx. 1 solar mass and characteristic temperature and density T approx. 750 K and approx. 2 x 10 to the 6th power per cu cm can be derived from the CO intensities. The density is too low by at least an order of magnitude for the observed infrared H sub 2 and far infrared CO emission to be accounted for by a purely hydrodynamic shock, and support is lent to hydromagnetic shock models. From the present measurements, the relative abundance of CO is estimated to be CO H sub 2 = 1.2 x .0001, corresponding to 20% of the cosmic abundance of C existing in the form of CO. The average relative abundance of OH in the shocked gas is O/H sub 2 or = 5 x 10 to the -7th power. An upper limit to the intensity of the HD J = 1 yields 0 line is used to derive an upper limit of tau or = 3 for the D/H relative abundance in the Orion cloud core.

Watson, D. M.↗
NASA NTRS1985

Far-infrared emission lines of CO and OH in the Orion-KL molecular shock

Observations of far infrared rotational emission lines which arise in the shocked gas associted with Orion-Kl are presented, including detections of the CO J = 34 yields 33, J = 31 yields 30, J = 26 yields 25, and OH sup 2 PI sub (3/2) J sup P = 7/2(-) yields 5/2(+) emission lines, as well as improved measurements of the CO J = 22 yields 21 and OH sup 2 PI sub (3/2) J = 5/2 yields 3/2 lines. These lines are observed to have velocity widths of Del V approx. 20 to 30 km/sec, somewhat less than either the 2 micro H sub 2 lines or the high velocity plateau component of the millimeter wave CO lines seen in this object. An H sub 2 column density of aprox. 3 x 10 to the 21st power, a total mass of approx. 1 solar mass and characteristic temperature and density T approx. 750 K and approx. 2 x 10 to the 6th power per cu cm can be derived from the CO intensities. The density is too low by at least an order of magnitude for the observed infrared H sub 2 and far infrared CO emission to be accounted for by a purely hydrodynamic shock, and support is lent to hydromagnetic shock models. From the present measurements, the relative abundance of CO is estimated to be CO H sub 2 = 1.2 x .0001, corresponding to 20 percent of the cosmic abundance of C existing in the form of CO. The average relative abundance of OH in the shocked gas is O/H sub 2 or = 5 x 10 to the -7th power. An upper limit to the intensity of the HD J - 1 yields 0 line is used to derive an upper limit of tau or = 3 for the D/H relative abundance in the Orion cloud core.

Watson, D. M.↗
DOE OSTI2022

Validation of ablation model for polyethylene using pulsed x-ray and proton exposures

The surface erosion of polyethylene is interrogated using pulsed x rays at the Z Machine (Sandia National Laboratories) and with proton beams at the Gamble II generator (Naval Research Laboratory) to validate a coupled model for volumetric thermal ablation, photoionization, finite-rate decomposition, and molecular recombination of radicals. The intense radiation pulses (up to ∼1014W/m2 over tens of nanoseconds) are used to generate one-dimensional vapor flows with low ionization fractions and a simplified geometry compared to typical laser ablation, allowing for evaluation of the model under local thermal equilibrium conditions. Areal momentum carried by the ensuing uniaxial hydrodynamic shock is used to indicate the extent of ablation. The threshold fluence for ablation is found to be in close correspondence with the bulk melt transition, and reasonable agreement with the model is obtained for peak temperatures in polyethylene up to 5500 K and heating rates up to 1011K/s where thermal decomposition reactions are also active.

Moore, N. W. (ORCID:0000000272259674)↗
DOE OSTI2023

Proton acceleration in an overdense hydrogen plasma by intense CO 2 laser pulses with nonlinear propagation effects in the underdense pre-plasma

We report on proton acceleration from intense CO 2 laser-irradiated hydrogen plasmas at near-critical densities, with the density gradient steepened by Nd:YAG laser ablation-driven hydrodynamic shocks. While the experimental results, such as the quasi-monoenergetic proton spectra and their scaling with respect to the laser energy, are generally in agreement with the simulations, certain laser shots produced significantly higher proton energies than anticipated during the experiment. The increased proton energy may be linked to nonlinear propagation effects in the steepened plasma density ramp before the critical surface, including relativistic self-focusing and, for the case of temporally-structured laser pulses observed in the experiment, focusing of the trailing pulse through the plasma channel formed by the leading pulse 25 ps ahead. The occurrence of channel focusing in the underdense hydrogen plasma is supported by a subsequent pump-probe experiment with a dark-field imaging technique, where the formation of ion channels was observed after the passage of an intense CO 2 laser pulse

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spray: Sparse Reductions of Arrays in OPEN MP

We present SPRAY, an open-source header-only C++ library for sparse reductions of arrays. SPRAY is meant for applications in which a large array is collaboratively updated by multiple threads using an associative and commutative operation such as +=. Especially when each thread accesses only parts of the array, SPRAY can perform significantly better than OPENMP's built-in reduction clause or atomic updates, while also using less memory than the former. SPRAY provides both an easy-to-use interface that can serve as a drop-in replacement for OPENMP reductions and a selection of reducer objects that accumulate the final result in different thread-safe ways. We demonstrate SPRAY through multiple test cases including the LULESH shock hydrodynamics code and a transpose-matrix-vector multiplication for sparse matrices stored in CSR format. SPRAY reductions outperform built-in OPENMP reductions consistently, in some cases improving run time and memory overhead by 20X, and even beating domain-specific approaches such as Intel MKI, by over 2X in some cases. Furthermore, SPRAY reductions have a minimal impact on the code base, requiring only a few lines of source code changes. Once in place, SPRAY reduction schemes can be switched easily, allowing performance portability and tuning opportunities by separating performance-critical implementation details from application code.

OpenMP↗
DOE OSTI2024

Introduction to Special Issue on the Early History of Nuclear Fusion

This introductory paper to the special issue of Fusion Science and Technology commemorates early research on fusion conducted at Los Alamos (the singular entity denoted Los Alamos Laboratory/Los Alamos Scientific Laboratory/Los Alamos National Laboratory at different times is designated “Los Alamos” in this paper) in support of the eventual H-bomb program. We survey the historical origins of the thermonuclear program, what was known of fusion reactions at the outbreak of the war, and the remarkable breakthroughs involving particularly the prospect of deuterium-tritium (DT) reactions conducted during the war, and we summarize the papers in this volume. Much of the nuclear fusion technical history presented herein has not been previously reported. Papers describe aspects of fusion science during these days, on shock hydrodynamics and on electron-radiation coupling, and on nuclear physics including the discoveries of resonances in both the DT cross section and in the lithium tritium-breeding cross section. Three papers follow our colleague Mark Paris’s finding Arthur Ruhlig’s 1938 paper on the first observation of DT fusion: one on how it influenced subsequent Manhattan Project research, another on a modern calculation of that historic experiment, and a third that has repeated the experiment using modern experimental capabilities. Other papers discuss how the first H-bomb test, Ivy Mike, led to the discovery of the new elements einsteinium and fermium and how the DT fusion processes played a key role in our universe’s development after the Big Bang. We also present a paper that analyzes the pioneering Cambridge University 1934 experiment by Marcus Oliphant, Paul Harteck, and Ernest Rutherford where deuterium-deuterium fusion was first observed and that describes how Ernest Lawrence missed identifying fusion in 1933. Finally, we present a summary of early concepts for controlled fusion energy that grew out of wartime discussions at Los Alamos. The papers show how J. Robert Oppenheimer played a leading technical role in the early developments of the H-bomb, before his later opposition—our first paper in this issue addresses the U.S. Department of Energy’s 2022 vacation of the earlier 1954 decision to revoke his security clearance.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1974

Possible evidence for structured acceleration of cosmic rays on a galactic scale from recent gamma-ray observations

Recent data from SAS-2 on the galactic gamma-ray line flux as a function of longitude reveal a broad maximum in the region below 30 deg. These data, as unfolded here, imply that the low-energy (1-10 GeV) galactic cosmic-ray flux varies with the radial distance from the galactic center and is about an order of magnitude higher than the local value in a toroidal region for radial distances between 4 and 5 kpc. We further show that this enhancement can be plausibly accounted for by Fermi acceleration and compression caused by a hydrodynamic shock driven by the expanding gas in the 3-kpc arm and invoked in some versions of galactic structure theory.

Stecker, F. W.↗
NASA NTRS1975

Comparison of magnetospheres and radio emissions of Jupiter with earth

The magnetosphere and radio emission of Jupiter is compared with those of the earth. It was predicted that Jupiter would have a Van Allen belt at a radius such that its magnetic field strength would be about equal to that in earth's Van Allen belt and that Jupiter's moon Io travels in the Van Allen belt. Because of Io's low conductivity, plasma sweeping past hits Io, producing a turbulent plasma proboscis which forms hydrodynamic shocks. These shocks travel down the magnetic field lines to the Jovian magnetosphere where they stimulate electron cyclotron emission and free radical spin-flip emission. The free radicals likely to exist abundantly and the richness of the likely decametric frequencies resulting from the many g values of the free radicals are discussed.

Libby, L. M.↗
NASA NTRS1981

Effect of thermal conduction and radiation on the dynamics of a flaring coronal loop

A systematic study of the effects of thermal conduction and radiation (Cox-Tucker type) on the dynamic response of a coronal loop is presented. The study is based on a one-fluid hydrodynamic description of a magnetic flux loop subject to a sudden localized heating, simulating the flare energy release. Numerical results are obtained for the heating of a full loop at one end and symmetrical heating at both ends. The results show that: (1) the basic response is the propagation of a hydrodynamic shock; (2) thermal conduction (Coulomb type) has significant effects; and radiation has only minor influence unless the density exceeds 10 to the 11th per cu cm. It is noted that the uniform flaring of a loop is associated with heating of long duration. Short pulse-like heating results in the nonuniform brightening and propagation of luminous fronts similar to that sometimes observed in flare loops.

Wu, S. T.↗
NASA NTRS1984

The laboratory simulation of unmagnetized supernova remnants Absence of a blast wave

Supernova remnants are experimentally simulated by irradiating spherical targets with eight-beam carbon dioxide laser in a chamber containing finite amounts of neutral gas, the gas being ionized by radiation from the hot target. The expansion velocities of the target plasmas are approximately the same as the expansion velocities of supernova ejecta and the experiment is successfully scaled to the case of a supernova remnant in an unmagnetized, low-density, interstellar medium. No sweep-up of the ambient plasma is detected, indicating that no hydrodynamic shock wave is formed to couple the target ejecta to the ambient gas. The experiment implies that if supernova ejecta couple to the interstellar medium, magnetic-field effects may be crucial to the physical description.

Borovsky, J. E.↗
NASA NTRS1987

Collapse of 9 solar mass stars

General relativistic hydrodynamical calculations of the collapse of O + Ne + Mg cores of a 9 solar mass star are reported. Collapse is induced by rapid electron captures as the O + Ne + Mg is burned to nuclear statistical equilibrium. The high entropy in the core after burning leads to a large abundance of free protons which readily capture electrons. This leads to large neutrino losses and a correspondingly small infalling homologous core. The hydrodynamic shock thus forms at a small mass point. The shock stalls before reaching the edge of the O + Ne + Mg core and thereby fails to produce a successful supernova explosion by the direct mechanism. No enhancement in the shock energy due to nuclear burning is found.

Baron, E.↗
NASA NTRS1992

Relativistic two-moment neutrino transport

We implement a general relativistic neutrino transport scheme appropriate to problems in stellar collapse and neutron star formation. In this scheme the spectral energy and momentum balance equations are solved as a coupled set. We perform transport calculations on two frozen stellar background models neglecting the feedback of the neutrinos on the matter. We compare with schemes that ignore either relativity or momentum balance or both. The work performed by the neutrino stress on the matter is taken into account and has an important effect. In particular, if there is a hydrodynamical shock in the system, this work contribution boosts the energy deposition behind the shock. This may enhance the chances for a successful neutrino driven supernova explosion.

Cernohorsky, J.↗
NASA NTRS1996

Interplanetary and Interstellar Dust Particles: Reentry Heating and Capture in Aerogel

Studies of interplanetary and interstellar dust can provide significant information on the evolution of the solar system or stars, respectively. However, for reliable analysis it is crucial to know how the particles have been modified during reentry (in the case of interplanetary dust particles, or [DP's) and impact into the capture medium. In the case of stratospheric capture, particles will be heated by atmospheric drag. Subsequent capture of the particles will result in heating, ablation, accretion of the capture medium and possible fragmentation. Modeling of these processes is a useful way of assessing their effects on the interpretation of the compositional data for these particles. Previous work on reentry heating has shown that heat diffusion alone cannot adequately account for temperature gradients observed in IDP's. In fact, for any reasonable thermal parameters, calculations show the particles to be nearly isothermal. Here we extend those studies to include phase transitions. These preliminary results are promising and show a significant temperature gradient for a 100 micron (diameter) particle. The actual capture of the particles in silica aerogel is being modeled using a comprehensive shock hydrodynamic code (called CTH). Various options of this code were investigated to attempt to make the most appropriate choice of methods of impact, equations of state, and processes of energy transfer from capture material to particle. The initial calculations with the code used only 'reasonable' estimates for the physical parameters of silica aerogel. Through the literature searches and personal contacts with the knowledgeable scientists, the best possible mechanical and thermal data have been made available for these computations.

Szydlik, Paul P.↗
NASA NTRS1997

Gamma-rays from supernova remnants and the signatures of diffusive shock acceleration

A nonlinear shock acceleration model which generates non-thermal proton distributions and includes a self-consistent determination of shock hydrodynamics, is considered. Gamma ray spectra are obtained for supernova remnants, allowing for the cessation of acceleration to high energies due to the finite ages and the sizes of the remnants. Gamma ray spectral cutoffs can be observed in the TeV range for reasonable remnant parameters and deviations from power law behavior are found at all energies from 1 MeV to cutoff. Correlated observations by the International Gamma Ray Astrophysics Laboratory, Whipple and other instruments may provide stringent constraints to understanding supernova remnants.

Baring, Matthew G.↗