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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Reynolds Number Effects on the Richtmyer-Meshkov Instability

This presentation compares the results of two very different experimental studies of Richtmyer-Meshkov instability: shock tube experiments in which an air/SF6 interface is accelerated by a weak shock wave; and incompressible experiments in which a box containing two different density miscible liquids is impulsively accelerated by bouncing it off of a fixed coil spring. Both experiments are initiated with sinusoidal initial perturbations. The interface perturbation initially remains sinusoidal as it grows in amplitude, but eventually the interfacial vorticity concentrates into points, forming a row of line vortices of alternating sign. The Reynolds number based on vortex circulation ranges from 1,000 to 45,000 in these experiments. It is found that viscous effects have a large, quantifiable effect on the evolution of the individual vortices. The effects of viscosity on the overall perturbation amplitude, however, are small and will be compared to theory.

Niederhaus, Charles↗
NASA NTRS2022

Finite Difference Methods for Turbulence Simulations

The optimal finite difference discretization used in simulations of turbulent flows is influenced by both, the type of the scale resolving simulation (DNS or LES), as well as the flow-physics (hydrodynamic instabilities, shocks, acoustics, etc.) one expects to resolve. Insight into dispersion and dissipation error requirements for some common scale-resolving simulation scenarios help to highlight the issues faced in selecting a scheme.

Finite Difference Methods↗
DOE OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Two-dimensional simulations of the ion/ion acoustic instability and electrostatic shocks

A newly developed 2D electrostatic code with particle ions and Boltzmann electrons is used to investigate the details of the ion/ion acoustic instability and the structure of electrostatic shocks. The simulation results show that, for the parameters relevant to the plasma sheet boundary layer, the saturation mechanism of the ion/ion acoustic instability is ion trapping. It is also shown that the 2D structure of electrostatic shocks is considerably different from that suggested by previous 1D simulations. The main reason for this difference is the presence of shock reflected ions, which through the ion/ion acoustic instability lead to the generation of large amplitude waves in the upstream region propagating obliquely to the shock normal. These waves play an important role in the shock dissipation process.

Karimabadi, H.↗
DOE OSTI2024

Richtmyer–Meshkov instability when a shock is reflected for fluids with arbitrary equation of state

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

Napieralski, Mario (ORCID:0009000692344901)↗
NASA NTRS1997

The Wardle Instability in Interstellar Shocks: Gas Temperture and Line Emission - 2

We have modeled the gas temperature structure in unstable C-type shocks and obtained predictions for the resultant CO and H2 rotational line emissions, using numerical simulations of the Wardle instability. Our model for the thermal balance of the gas includes ion-neutral frictional heating; compressional heating; radiative cooling due to rotational and ro-vibrational transitions of the molecules CO, H2O, and H2; and gas-grain collisional cooling. We obtained results for the gas temperature distribution in-and H2 and CO line emission from-shocks of neutral Alfvenic Mach number 10 and velocity 20 or 40 km/ s in which the Wardle instability has saturated. Both two- and three-dimensional simulations were carried out for shocks in which the preshock magnetic field is perpendicular to the shock propagation direction, and a two-dimensional simulation was carried out for the case in which the magnetic field is obliquely oriented with respect to the shock propagation direction. Although the Wardle instability profoundly affects the density structure behind C-type shocks, most of the shock-excited molecular line emission is generated upstream of the region where the strongest effects of the instability are felt. Thus the Wardle instability has a relatively small effect on the overall gas temperature distribution in-and the emission-line spectrum from-C-type shocks, at least for the cases that we have considered. In none of the cases that we have considered thus far did any of the predicted emission-line luminosities change by more than a factor of 2.5, and in most cases the effects of instability were significantly smaller than that. Slightly larger changes in the line luminosities seem likely for three-dimensional simulations of oblique shocks, although such simulations have yet to be carried out and lie beyond the scope of this study. Given the typical uncertainties that are always present when model predictions are compared with real astronomical data, we conclude that Wardle instability does not imprint any clear observational signature on the shock-excited CO and H2 line strengths. This result justifies the use of one-dimensional steady shock models in the interpretation of observations of shock-excited line emission in regions of star formation. Our three-dimensional simulations of perpendicular shocks revealed the presence of warm filamentary structures that are aligned along the magnetic field, a result that is of possible relevance to models of water maser emission from C-type shocks.

Neufeld, David A.↗
NASA NTRS1993

The stability of cosmic-ray-dominated shocks - A secondary instability

Previous studies in one spatial dimension showed that the precursors of cosmic-ray-dominated shocks are unstable against traveling acoustic disturbances. Here we report that a secondary, Rayleigh-Taylor type instability can exist in cosmic-ray-dominated media influenced by the acoustic instability. Using the local WKB approximation, the growth rate of the secondary instability is shown to be comparable to that of the one-dimensional acoustic instability itself in the cases we have considered. The nonlinear development has been followed numerically with a two-dimensional PPM hydrodynamics code that also incorporates the two-fluid cosmic-ray energy equation. We show that the secondary instability may cause the precursor and postshock flows to become highly turbulent.

Ryu, Dongsu↗
NASA NTRS1994

Nonlinear instability of accelerating shock waves with application to supernovae

We consider the stability of a planar accelerating shock front in an exponential atmosphere, a situation that can be described by a self-similar solution. A previous linear instability analysis showed that there are three regimes, depending on the wavelength of the perturbation along the direction of the shock front: at long wavelengths the shock is unstable, at intermediate wavelengths it is overstable, and at short wavelenghts it is stable, where the characteristic length is set by the initial density scale height. We have carried out numerical simulations of such an accelerating shock front and have confirmed the results of the semianalytic linear analysis in the linear regime. In the nonlinear regime, the evolution again depends on the three wavelength domains. At long wavelengths, the instability continuously grows but the flow remains smooth, and at short wavelengths, the flow is stable. The intermediate wavelength, overstable regime shows more complex evolution. The growing linear oscillations saturate soon after entering the nonlinear regime and the oscillations continue with the same approximate period. The shock front develops moving points of intersection, which generate weak shock fronts and density and pressure structure in the immediate postshock flow. The density fluctuations, with a contrast of a factor of 2-3, become frozen into the downstream flow. The growth of the instability is slow, so substantial initial perturbations are needed; thesemay be present in Type II supernovae with red supergiant progenitor stars which have outer convective regions. The clumping in the outer supernova atmosphere may affect spectral line formation and may play a role in the formation of fast knots.

Luo, Ding↗
NASA NTRS1983

Temperatures of shock-induced shear instabilities and their relationship to fusion curves

New emission spectra for MgO and CaAl2Si2O8 (glass) are observed from 430 to 820 nm. Taken with previous data, it is suggested that transparent solids display three regimes of light emission upon shock compression to successively higher pressures: (1) characteristic radiation such as observed in MgO and previously in other minerals, (2) heterogeneous hot spot (greybody) radiation observed in CaAl2Si2O8 and previously in all transparent solids undergoing shock-induced phase transformations, and (3) blackbody emission observed in the high pressure phase regime in NaCl, SiO2, CaO, CaAl2Si2O8, and Mg2SiO4. The onset of the second regime may delineate the onset of shock-induced polymorphism whereas the onset of the third regime delineates the Hugoniot pressure required to achieve local thermal equilibrium in the shocked solid. It is also proposed that the hot spot temperatures and corresponding shock pressures determined in the second regime delineate points on the fusion curves of the high pressure phase.

Schmitt, D. R.↗
NASA NTRS1992

Kinetic structure of slow shocks - Effects of the electromagnetic ion/ion cyclotron instability

The structure of slow magnetosonic shocks in the low beta regime is analyzed with attention given to ion heating and the effects of waves upstream of the electromagnetic ion/ion cyclotron (EMIIC) instability. Shock formation is assessed by means of three methods - a relaxation method and two based on dynamic flow interactions - to determine the effects of initialization and boundary conditions on the formation. Good solutions are found with the piston method and the similar flow-flow method in which the plasma is injected from two boundaries to form two slow shocks. Plasma parameters and shock normal angle are found to be the key variables dictating the structure of the magnetosonic shocks. Four unique classes of resultant shock structures are described in which classical, steady, or nonsteady behavior is found. The analysis also yields insight into the relationship between EMIIC instability and ion dissipation.

Omidi, N.↗
DOE OSTI2022

Unveiling the origins of work-hardening enhancement and mechanical instability in laser shock peened titanium

Laser shock-peening (LSP) produces complex gradients in microstructures and residual stresses, which enables fabrication of metallic engineering components with superior mechanical properties. In this work, we use non-destructive and spatially resolved high-energy X-ray diffraction (HE-XRD) to investigate the influences of LSP on microstructure, surface topography, and residual stress for high purity titanium plates. LSP is found to produce large compressive in-plane residual stresses near the peened surface which monotonically decay to zero about 2.5 mm below the surface. These properties were also tracked during in-situ tensile loading, allowing stress partitioning and work-hardening rates to be measured as a function of depth. The surface region is found to have the highest work-hardening rate and remain mechanically stable up to sample failure. In addition, a crystal rigid rotation along a transversal axis of about 27 & nbsp; was observed at large applied strain near the surface region, which is attributed to the formation of denser dispersed shear bands. Furthermore, a sudden drop in tensile stress was observed at a depth of similar to 160 mu m from the LSP surface, indicating the formation of a localized mechanical instability zone (LMIZ). The origin of this LMIZ could be attributed to the increase in shear stress induced by a load transfer from LSP surface due to the formation of denser dispersed shear bands. These observations indicate a localized transition in deformation behavior has occurred from dispersed shear banding to LMIZ and homogeneous deformation.

36 MATERIALS SCIENCE↗
NASA NTRS1994

A Fully Implicit Time Accurate Method for Hypersonic Combustion: Application to Shock-induced Combustion Instability

A new fully implicit, time accurate algorithm suitable for chemically reacting, viscous flows in the transonic-to-hypersonic regime is described. The method is based on a class of Total Variation Diminishing (TVD) schemes and uses successive Gauss-Siedel relaxation sweeps. The inversion of large matrices is avoided by partitioning the system into reacting and nonreacting parts, but still maintaining a fully coupled interaction. As a result, the matrices that have to be inverted are of the same size as those obtained with the commonly used point implicit methods. In this paper we illustrate the applicability of the new algorithm to hypervelocity unsteady combustion applications. We present a series of numerical simulations of the periodic combustion instabilities observed in ballistic-range experiments of blunt projectiles flying at subdetonative speeds through hydrogen-air mixtures. The computed frequencies of oscillation are in excellent agreement with experimental data.

Yungster, Shaye↗
DOE OSTI2022

Design optimization for Richtmyer–Meshkov instability suppression at shock-compressed material interfaces

We report the Richtmyer–Meshkov instability (RMI) is a phenomenon that occurs at the interface of two substances of different densities due to an impulsive acceleration, such as a shock wave passing through this interface. Under these conditions, the instability can be seen as interface perturbations begin to grow into narrow jets or spikes of one substance that propagate into the other. In some cases, this interface may involve an elastic–plastic material, which can play a significant role in the development and behavior of the RMI. The ability to effectively control RMI jetting and spike growth is one major limiting factor in technological challenges, such as inertial confinement fusion, that involve using high-pressure shock waves to implode a fuel target. The propagation of RMI growth can lead to increased asymmetry in this implosion process and significantly reduce the obtained energy yield. We use hydrodynamics simulations of impactor shock-compression experiments and methods based in design optimization to suppress RMI spike growth by altering the geometry and other properties of a shock-compressed elastic–plastic material target that shares an interface with atmospheric air. These hydrodynamics simulations use an arbitrary Lagrangian–Eulerian method with a high-order finite element approach. Our results demonstrate that RMI suppression can be achieved by intentionally creating a separate upstream interface instability to counteract the growth of long narrow RMI spikes at an interface with initial perturbations.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Cyclotron drift instability in the bow shock.

It is argued that the cyclotron drift instability (coupled ion waves and electron Bernstein modes) can explain many structural features of the weak electrostatic turbulence observed in bow shock magnetic-field gradients and thus that it provides a more attractive speculation than the ion acoustic or Buneman instabilities previously suggested as weak turbulence sources.

Wu, C. S.↗
DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗
NASA NTRS1985

Electromagnetic ion beam instabilities - Hot beams at interplanetary shocks

This paper considers Vlasov instabilities driven by a very hot ion beam streaming parallel to a magnetic field. The linear theory of electromagnetic instabilities driven by such a beam in a homogeneous, collisionless, nonrelativistic plasma is used. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions are presented. If the thermal speed of the beam is greater than its drift speed, the two dominant modes are the right-hand and left-hand resonant ion beam instabilities. The parametric dependencies of the threshold drift speeds and growth rates for both modes are presented. In particular, the thresholds are shown to be sensitive functions of the beam temperature anisotropies. Arguments are presented that these two modes are driven by the suprathermal ion component at interplanetary shocks, and that the growth of these modes is sufficient to account for the amplitudes of MHD-like waves observed at such shocks.

Gary, S. P.↗
NASA NTRS1998

3D Simulations of the Richtmyer-Meshkov Instability with Re-Shock

We present results of inviscid simulations, in three dimensions, of Richtmyer-Meshkov instability for high incident shock Mach number. The growth rate of a single harmonic perturbation is quantified and compared with the results of a 2D calculation. Upon re-shock, the perturbation amplitude undergoes a phase reversal while the mean velocity of the interface is zero. Before re-shock the normalized growth rate of a 2D and 3D interface are nearly the same, but the growth rate after re-shock is significantly larger for the 3D than the 2D case. We also examine the evolution of multiple harmonic perturbations. Computational and parallelization issues of the simulation code will also be briefly discussed. The computations were done on the T3E at Pittsburgh Supercomputing Center.

Meiron, Daniel I.↗