Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reshock”

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.

50 records · Page 3

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Temporal evolution of scalar modes in Richtmyer–Meshkov instability of inclined interface using high-speed PIV and PLIF measurements at 60 kHz

Abstract The current work presents simultaneous, high-speed measurements at 60,000 fields per second of velocity and mole fraction using particle image velocimetry (PIV) and planar laser induced acetone-fluorescence in a Richtmyer–Meshkov instability of an inclined interface (Atwood number, At = 0.22). Specifically, around 2 ms of temporal evolution of the vortex structures and their associated scalar modes immediately following the interface-reshock interaction is presented. Two initial interface conditions are discussed—(a) a sharp, inclined ‘single mode’ interface and (b) a ‘multi-mode’ interface where small perturbations are imposed on the single mode case. A 2D wavelet decomposition of the scalar flow field shows a highly intermittent distribution of small-scale variance throughout the interface even at late times. These are correlated strongly with the vortex structures and local turbulence intensity, where each small-scale scalar mode is sandwiched between two co-rotating vortex structures. This indicates that the interstitial regions between the vortices are significant hotspots of entrainment, which is then dispersed by the induced, counter-flow velocity fields. The multimode case demonstrates similar organization at large scales, while the scalar field is much more homogeneous at smaller scales. These observations highlight the importance of capturing the early time vortex evolution to accurately estimate any late time intermittency, especially where deposition of intense vorticity on sharp interfaces is present.

Engineering↗

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Argon equation of state data to 1 TPa: Shock compression experiments and simulations

Argon is the most abundant noble gas on Earth and its noble, atomic fluid nature makes it an excellent candidate for comparison of experiment and theory at extreme conditions. For this study, we performed a combined computational and experimental study on shock compressed cryogenic liquid argon. Using Sandia's Z machine, we shock compressed liquid argon to 600 GPa and reshock states up to 950 GPa. Laser shock experiments at the Omega Laser facility extend the principal Hugoniot to 1000 GPa and provided temperature data along the principal Hugoniot. The plate impact experiments and laser shock experiments used well-characterized impedance matching standards and demonstrate consistent results between the two platforms over a common range. Density functional theory based molecular dynamics simulations provided additional data on the Hugoniot to 600 GPa. The combined experimental data and simulation results provide constraints on the development of new equation of state models at extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Developing quartz and molybdenum as impedance-matching standards in the 100-Mbar regime

We present quartz and molybdenum single-shock Hugoniot data to 62 Mbar and 120 Mbar, respectively. Quartz and molybdenum are attractive standard materials for Hugoniot equation-of-state impedance-matching studies. Quartz is advantageous because a reflecting shock front is supported above ~1 Mbar whose velocity can be tracked in situ to very high accuracy. Molybdenum's relatively high impedance makes it a more appropriate standard for other high-impedance material studies to avoid large departures from the standard principal Hugoniot upon release or reshock into the sample. Extensive measurements already exist for both materials up to about 10 Mbar and 10 eV, but new high-energy-density facilities and techniques have extended the experimentally accessible pressure range to >100 Mbar. The data presented here indicate which modern equation-of-state models can be reliably extended into newly-accessible experimental regimes.

36 MATERIALS SCIENCE↗

Two-length-scale turbulence model for self-similar buoyancy-, shock-, and shear-driven mixing

The three-equation k-L-a turbulence model [B. Morgan and M. Wickett, Three-equation model for the self-similar growth of Rayleigh-Taylor and Richtmyer-Meshkov instabilities," Phys. Rev. E 91 (2015)] is extended by the addition of a second length scale equation. It is shown that the separation of turbulence transport and turbulence destruction length scales is necessary for simultaneous prediction of the growth parameter and turbulence intensity of a Kelvin-Helmholtz shear layer when model coeficients are constrained by similarity analysis. Constraints on model coeficients are derived that satisfy an ansatz of self-similarity in the low-Atwood-number limit and allow the determination of model coeficients necessary to recover expected experimental behavior. The model is then applied in one-dimensional simulations of Rayleigh-Taylor, reshocked Richtmyer-Meshkov, Kelvin{Helmholtz, and combined Rayleigh-Taylor/Kelvin-Helmholtz instability mixing layers to demonstrate that the expected growth rates are recovered numerically. Finally, it is shown that model behavior in the case of combined instability is to predict a mixing width that is a linear combination of Rayleigh-Taylor and Kelvin-Helmholtz mixing processes.

42 ENGINEERING↗

NIF 3-pager FY24: Mshock [Slides]

Complete initial condition scan examining surface roughness alteration of larger mode growth (begin FY23 Q2). Previous reshock experiments have shown that we can measure the change in growth rate of a larger carrier mode as a diagnostic for the small-scale non-linear instability growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mshock HED Council Request FY24 [Slides]

ICF and HEDP experiments experience complicated loadings including successive shocks, but modeling and analytic work has mostly engaged the simpler shock or shock-reshock case. This is because the co propagating case is difficult to achieve with conventional (non-HED) drivers. Successive shocks will be challenging to model in BHR or modal model interpretations: 1) Ex: Consider the case where long wavelength modes re-invert coherently, but short wavelength nonlinear modes are spun up turbulently. 2) Necessitates more advanced diagnostics and analysis than simple “mixing layers.” High resolution spectral information (concomitant with Ω-EP campaigns) and higher-order moment analyses will be necessary to stress and validate the reduced BHR-type models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Cygnus loop - A detailed comparison of X-ray and optical emission

A comparison at a resolution of 17 arcsec is presented between optical emission from forbidden O III and forbidden S II and the thermal X-ray emission for a field on the southeast edge of the Cygnus Loop SNR. The relationship between optical and X-ray emission in the field is described, generalizing to the Loop as a whole when possible. Several possible explanations for the presence of bright X-ray emission in the vicinity of optical emission are presented and evaluated based on the data. These mechanisms involve evaporation, gradual variations and gradients in the density of the preshock intercloud medium, and additional compression of material which has already been heated to X-ray temperatures by the adiabatic blast wave. The additional compression could result from rapid deceleration of the blast wave itself or from reshocking by reflected or bow shocks around dense clouds. Implications of the observations for models of SNR evolution and the interstellar medium are discussed.

Cox, D. P.↗

The formation of interstellar jets by the convergence of supersonic conical flows

The interaction of a stellar wind with a nonuniform environment leads, under the assumption of steady state, to the formation of an ovoidal cavity with acute ending tips. The stellar wind recollected by the walls of the cavity ends up being deposited at the tips. Here, it is shown that this focusing effect leads to the formation of a narrow cylindrical stream or 'jet' of reshocked stellar wind matter moving directly away from the star with a large velocity. A 'typical' T Tauri star may produce jets with densities of 1000-10,000/cu cm, velocities about 100 km/s, widths about 0.00001-0.001 pc over a length of 0.001-0.01 pc. The opening angles are 7 deg or less.

Canto, J.↗

Evaluating turbulence models at high energy densities

The RESHOCK campaign at Lawrence Livermore National Laboratory has been working on increasing our understanding of the evolution of turbulent, unstable plasma interfaces which are applicable in various HED applications including ICF implosions. A common approach to model such interfaces involves the use of Renolds-Averaged-Navier-Stokes (RANS) models whose parameters have been constrained by theory and experiments. These parameter sets are not unique and the applicability of one set of tuned model parameters to systems with differing shock strengths or material densities has not been well tested. Our new NIF experiments, along with related studies of turbulent shear flows, provide data and model validation for plasma interfaces in the HED regime. Our experiments measure mixing-layer width at an unstable interface in the plasma regime and are specifically designed to challenge mix models of the Reynolds-Averaged-Navier-Stokes (RANS) type where turbulence is assumed to be fully developed, i.e. the flow exhibits a broad spectrum of length scales without memory of the initial condition.4 We utilize precise control of the initial interface conditions (densities on both sides of the interface, an initial rippled perturbation pattern consisting of equally weighted wavelengths between 10 and 20 microns, and initial material compositions) along with a repeatable drive history

36 MATERIALS SCIENCE↗

Impact of ablator thickness and laser drive duration on a platform for supersonic, shockwave-driven hydrodynamic instability experiments

Here, we discuss changes to a target design that improved the quality and consistency of data obtained through a novel experimental platform that enables the study of hydrodynamic instabilities in a compressible regime. The experiment uses a laser to drive steady, supersonic shockwave over well-characterized initial perturbations. Early experiments were adversely affected by inadequate experimental timescales and, potentially, an unintended secondary shockwave. These issues were addressed by extending the 4 x 10 13 W/cm 2 laser pulse from 19 ns to 28 ns, and increasing the ablator thickness from 185 µm to 500 µm. We present data demonstrating the performance of the platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗