Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock instability”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Assessing the performance of MagLIF with 3D MHD simulations

The Magnetized Liner Inertial Fusion (MagLIF) experimental platform at Sandia National Laboratories has realized a number of recent improvements in electrical current and laser preheat coupling but fusion yields are below expectations, based on modeling with clean (that is, simulations with no interfacial instabilities and no mix) two-dimensional (2D) magneto-hydrodynamics (MHD) simulations, by a factor of 5 or more. However, magnetized liner implosions are known to produce helical magneto-Rayleigh–Taylor (MRT) modes and complex stagnation structures that cannot be directly modeled in 2D. This paper presents the results of 3D HYDRA MHD simulations, including helical MRT, showing that degradation from these instabilities can readily reproduce experimental yields, but it is difficult to simultaneously match the yield along with inferred fuel temperature, pressure, and burn history (assessed with x-rays), which are larger in the experiments considered here. Additional analysis methods and future experiments are proposed to help address the discrepancies. The 3D simulations also show, without substantive improvements to stability, changing the applied Bz or preheat is unlikely to increase performance of MagLIF to the same degree as clean 2D simulations. Finally, the first 3D HYDRA simulations with the Hall term show the ab initio production of helical MRT and produce a clear change in stagnation morphology compared to pre-seeded simulations without Hall. In spite of the differences, the simulations still produce comparable fusion performance but with longer, lower power, neutron yield history.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Flow visualization in high-energy-density laboratory plasmas

Flow visualization is often essential for evaluating high-energy-density (HED) laboratory plasmas, examining their dynamics, and assessing design strategies. However, imaging at fine phenomenological flow scales is a major challenge. Here, a zone-plate–based radiography system tailored for applications in HED physics is demonstrated. The burst-mode radiography system combines a dual-flash laser-generated x-ray source, a high-magnification zone plate optic, and a multiple frame hybrid complementary metal-oxide-semiconductor (hCMOS) sensor for image detection. Operating at a magnification of 42 x⁠, a pair of radiographs is measured in rapid succession, each time gated by a 4.75 keV flash of x rays lasting 100 ps. The system achieves a spatial resolution below 3 μm. The potential of this approach for applications in HED physics is exemplified by observing small-scale spike morphologies at a pusher-foam interface impacted by a strong, laser-driven blast wave.

Flow instabilities↗
DOE OSTI2025

Optimizing the Smoothness and Thickness Uniformity of Thin-Film Parylene-N Vapor-Deposited Coatings for Inertial Confinement Fusion Experiments

Polymer coatings with submicrometer smoothness and constant thickness are a required component in a variety of inertial confinement fusion experiments. Smoothness is important for minimizing Rayleigh-Taylor-driven hydrodynamic instabilities, and uniform thickness is important for uniform shock propagation and shell convergence, both of which are critical phenomena that affect the experiment. The preferred polymer coating method is to vapor deposit the parylene-N polymer because it provides nominally smooth conformal coatings. As the coating thickness exceeds 5 µm, however, dome-shaped nodular growth defects develop and the thickness will vary by up to 17% over a distance of 3 cm. This study presents a deterministic method for achieving uniform film thicknesses with ±2% variability over 3 cm and a predictive method to control the thickness to within 5% of the desired value. A coating smoothness of ∼50 nm rms, measured over 40 000 µm 2 , was achieved by adding additional surfaces near the substrates. This additional area improved the thickness uniformity, an effect that is attributed to the low sticking coefficient of the parylene monomer.

chemical vapor deposition (CVD)↗
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↗
NASA NTRS1973

Research in plasma physics

Three aspects of barium ion cloud dynamics are discussed. First, the effect of the ratio of ion cloud conductivity to background ionospheric conductivity on the motion of barium ion clouds is investigated and compared with observations of barium ion clouds. This study led to the suggestion that the conjugate ionosphere participates in the dynamics of barium ion clouds. Second, analytic work on the deformation of ion clouds is presented. Third, a linearized stability theory was extended to include the effect of the finite extent of an ion cloud, as well as the effect of the ratio of ion cloud to ionospheric conductivities. The stability properties of a plasma with contra-streaming ion beams parallel to a magnetic field are investigated. The results are interpreted in terms of parameters appropriate for collisionless shock waves. It is found that this particular instability can be operative only if the up-stream Alfven Mach number exceeds 5.5.

Source record↗
NASA NTRS1990

Real-time measurement of materials properties at high temperatures by laser produced plasmas

Determination of elemental composition and thermophysical properties of materials at high temperatures, as visualized in the context of containerless materials processing in a microgravity environment, presents a variety of unusual requirements owing to the thermal hazards and interferences from electromagnetic control fields. In addition, such information is intended for process control applications and thus the measurements must be real time in nature. A new technique is described which was developed for real time, in-situ determination of the elemental composition of molten metallic alloys such as specialty steel. The technique is based on time-resolved spectroscopy of a laser produced plasma (LPP) plume resulting from the interaction of a giant laser pulse with a material target. The sensitivity and precision were demonstrated to be comparable to, or better than, the conventional methods of analysis which are applicable only to post-mortem specimens sampled from a molten metal pool. The LPP technique can be applied widely to other materials composition analysis applications. The LPP technique is extremely information rich and therefore provides opportunities for extracting other physical properties in addition to the materials composition. The case in point is that it is possible to determine thermophysical properties of the target materials at high temperatures by monitoring generation and transport of acoustic pulses as well as a number of other fluid-dynamic processes triggered by the LPP event. By manipulation of the scaling properties of the laser-matter interaction, many different kinds of flow events, ranging from shock waves to surface waves to flow induced instabilities, can be generated in a controllable manner. Time-resolved detection of these events can lead to such thermophysical quantities as volume and shear viscosities, thermal conductivity, specific heat, mass density, and others.

Kim, Yong W.↗
NASA NTRS1990

Physics of the outer heliosphere; Proceedings of the 1st COSPAR Colloquium, Warsaw, Poland, Sept. 19-22, 1989

This volume includes chapters on spectroscopic data on the local interstellar medium and the related XUV radiation background; solar UV backscatter on neutral Galactic gases; entry and dynamics of Galactic and anomalous cosmic rays in the heliosphere; distant solar wind plasma, magnetic field, and solar energetic particles; the interaction of solar wind with the interstellar medium; and deep-space missions. Papers are on the characteristics and the velocity structure of the local interstellar medium, the scattering of solar UV on local neutral gases, solar effects on underground muons at 570 hg/sq cm, and cosmic rays and magnetosonic instabilities of solar wind flow near the heliospheric shock wave. Consideration is also given to the shock heating of the solar wind plasma, solar wind vortex flow in the outer heliosphere, the signature of a viscous interaction at the heliopause, the expected beams of energetic neutral atoms in the outer heliosphere, and the Pioneers 10 and 11 deep-space missions.

Grzedzielski, S.↗
NASA NTRS1992

A new class of galactic discrete gamma ray sources: Chaotic winds of massive stars

We propose a new class of galactic discrete gamma-ray sources, the chaotic, high mass-loss-rate winds from luminous early-type stars. Early-type stellar winds are highly unstable due to intrinsic line-driven instabilities, and so are permeated by numerous strong shocks. These shocks can accelerate a small fraction of thermal electrons and ions to relativistic energies via the first-order Fermi mechanism. A power-law-like photon spectrum extending from keV to above 10 MeV energies is produced by inverse Compton scattering of the extremely abundant stellar UV photons by the relativistic electrons. In addition, a typical pi(sup 0)-decay gamma-ray spectrum is generated by proton-ion interactions in the densest part of the winds.

Chen, Wan↗
NASA NTRS1993

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗
NASA NTRS1994

The peculiar type II supernova 1993J in M81: Transition to the nebular phase

We present optical spectra of the bright, peculiar Type II supernova 1993J in M81 spanning the first 14 months of its existence, revealing its transition to the nebular phase. Unlike the case in normal Type II supernovae, during the first 2-10 months the H-alpha emission line gradually becomes less prominent relative to other features such as (O I) lambda lambda 6300, 6364 and (Ca II) lambda lambda 7291, 7324, as we had predicted based on early-time (tau less than or approximately equal to 2 months) spectra. The nebular spectrum resembles those of the Type Ib/Ic supernovae 1985F and 1987M, although weak H-alpha emission is easily visible even at late times in SN 1993J. At tau = 8 months a close similarity is found with the spectrum of SN 1987K, the only other Type II supernova known to have undergone such a metamorphosis. The emission lines are considerably broader than those of normal Type II supernovae at comparable phases, consistent with the progenitor having lost a majority of its hydrogen envelope prior to exploding. Consequently, there is now little doubt that Type Ib, and probably Type Ic, supernovae result from core collapse in stripped, massive stars; models of the chemical evolution of galaxies in which these subtypes are ascribed to exploding white dwarfs must be appropriately modified. Although all of the emission lines in spectra of SN 1993J fade roughly exponentially for a considerable time, the fading of H-alpha begins to slow down at tau approximately = 8 months, and in the interval tau = 10-14 months its flux is constant, or even slightly rising in the wings of the line. This behavior, together with the box-like shape and great breadth (full width at half maximum (FWHM) approximately = 17 000 km/s) of the line profile, suggests that the H-alpha emission is being produced by the high-velocity outer layer of hydrogen ejecta interacting with circumstellar gas released by the progenitor prior to its explosion. A similar phenomenon has previously been seen at later phases in several Type II supernovae, most notably SN 1980K. Bumps (FWHM approximately = 1000 km/s, amplitude approximately = 20%) in the H-alpha profile are probably indicative of Rayleigh-Taylor instabilities in the cool gas behind the reverse shock. A very narrow component (unresolved, FWHM less than or approximately equal to 200 km/s) of H-alpha at the symmetric velocity of SN 1993J may instead be produced by a superposed H II region, or perhaps by recombination in a large circumstellar shell or ring that was ionized during the first few hours after outburst. In the near future the spectrum of SN 1993J should become increasingly dominated by broad H-alpha emission.

Filippenko, Alexei V.↗
NASA NTRS1998

Development of an Aeroelastic Analysis Including a Viscous Flow Model

Under this grant, the flutter versions of the three-dimensional Navier-Stokes aeroelastic code (TURBO-AE) have been developed. The TURBO-AE aeroelastic code allows flutter calculations for a fan, compressor, or turbine blade row. This code models a vibrating three-dimensional bladed disk geometry and the associated unsteady flow (including shocks, and viscous effects) to calculate the aeroelastic instability using a work-per-cycle or eigenvalue approach. Steady and unsteady calculations have been performed with the 3D Navier-Stokes (TURBO) code to gain familiarity with the code. Calculations were done with the NASA Rotor-67 configuration which has been used as a standard test case by turbomachinery researchers. Additional calculations have been performed for the NASA/GE EEE fan and a proprietary configuration.

Keith, Theo G., Jr.↗
NASA NTRS2001

Development of an Aeroelastic Analysis Including a Viscous Flow Model

Under this grant, Version 4 of the three-dimensional Navier-Stokes aeroelastic code (TURBO-AE) has been developed and verified. The TURBO-AE Version 4 aeroelastic code allows flutter calculations for a fan, compressor, or turbine blade row. This code models a vibrating three-dimensional bladed disk configuration and the associated unsteady flow (including shocks, and viscous effects) to calculate the aeroelastic instability using a work-per-cycle approach. Phase-lagged (time-shift) periodic boundary conditions are used to model the phase lag between adjacent vibrating blades. The direct-store approach is used for this purpose to reduce the computational domain to a single interblade passage. A disk storage option, implemented using direct access files, is available to reduce the large memory requirements of the direct-store approach. Other researchers have implemented 3D inlet/exit boundary conditions based on eigen-analysis. Appendix A: Aeroelastic calculations based on three-dimensional euler analysis. Appendix B: Unsteady aerodynamic modeling of blade vibration using the turbo-V3.1 code.

Keith, Theo G., Jr.↗
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
NASA NTRS2022

Revolutionizing Our Understanding of Particle Energization in Space Plasmas Using on-Board Wave-Particle Correlator Instrumentation

A leap forward in our understanding of particle energization in plasmas throughout the heliosphere is essential to answer longstanding questions in heliophysics, including the heating of the solar corona, acceleration of the solar wind, and energization of particles that lead to observable phenomena, such as the Earth’s aurora. The low densities and high temperatures of typical heliospheric environments lead to weakly collisional plasma conditions. Under these conditions, the energization of particles occurs primarily through collisionless interactions between the electromagnetic fields and the individual plasma particles with energies characteristic of a particular interaction. To understand how the plasma heating and particle acceleration impacts the macroscopic evolution of the heliosphere, impacting phenomena such as extreme space weather, it is critical to understand these collisionless wave-particle interactions on the characteristic ion and electron kinetic timescales. Such understanding requires high-cadence measurements of both the electromagnetic fields and the three-dimensional particle velocity distributions. Although existing instrument technology enables these measurements, a major challenge to maximize the scientific return from these measurements is the limited amount of data that can be transmitted to the ground due to telemetry constraints. A valuable, but underutilized, approach to overcome this limitation is to compute on-board correlations of the maximum-cadence field and particle measurements to improve the sampling time by several orders of magnitude. Here we review the fundamentals of the innovative field-particle correlation technique, present a formulation of the technique that can be implemented as an on-board wave-particle correlator, and estimate results that can be achieved with existing instrumental capabilities for particle velocity distribution measurements.

plasma heating↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
NASA NTRS2003

Two-Fluid Interface Instability Being Studied

The interface between two fluids of different density can experience instability when gravity acts normal to the surface. The relatively well known Rayleigh-Taylor (RT) instability results when the gravity is constant with a heavy fluid over a light fluid. An impulsive acceleration applied to the fluids results in the Richtmyer-Meshkov (RM) instability. The RM instability occurs regardless of the relative orientation of the heavy and light fluids. In many systems, the passing of a shock wave through the interface provides the impulsive acceleration. Both the RT and RM instabilities result in mixing at the interface. These instabilities arise in a diverse array of circumstances, including supernovas, oceans, supersonic combustion, and inertial confinement fusion (ICF). The area with the greatest current interest in RT and RM instabilities is ICF, which is an attempt to produce fusion energy for nuclear reactors from BB-sized pellets of deuterium and tritium. In the ICF experiments conducted so far, RM and RT instabilities have prevented the generation of net-positive energy. The $4 billion National Ignition Facility at Lawrence Livermore National Laboratory is being constructed to study these instabilities and to attempt to achieve net-positive yield in an ICF experiment.

Niederhaus, Charles E.↗