Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instabilities”

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 91 records · Page 5

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↗

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗

Observation of Hydrodynamic Flows in Imploding Fusion Plasmas on the National Ignition Facility [plus supplemental information]

Inertial confinement fusion implosions designed to have minimal fluid motion at peak compression often show significant linear flows in the laboratory, attributable per simulations to percent-level imbalances in the laser drive illumination symmetry. In this work, we present experimental results which intentionally varied the mode 1 drive imbalance by up to 4% to test hydrodynamic predictions of flows and the resultant imploded core asymmetries and performance, as measured by a combination of DT neutron spectroscopy and high-resolution x-ray core imaging. Neutron yields decrease by up to 50%, and anisotropic neutron Doppler broadening increases by 20%, in agreement with simulations. Furthermore, a tracer jet from the capsule fill-tube perturbation that is entrained by the hot-spot flow confirms the average flow speeds deduced from neutron spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗

Hall instability driven seeding of helical magneto-Rayleigh–Taylor instabilities in axially premagnetized thin-foil liner Z-pinch implosions

Helical magneto-Rayleigh–Taylor instability (MRTI) structures have been observed in z-pinch-driven liner implosion experiments with a pre-imposed axial magnetic field. We show that the formation of these helical structures can be described by a Hall magnetohydrodynamical (HMHD) model. We used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) [Seyler and Martin, Phys. Plasmas 18, 012703 (2011)] to study these helical instabilities and show that a Hall interchange instability in low-density coronal plasma immediately surrounding the dense liner is responsible for producing helically oriented effects in the magnetic field and current density within the coronal layer. Here, this seeds the helical pitch angle of the MRTI even when other proposed helical seeding mechanisms are either not present in the experiments or not accounted for in the simulations. For example, this mechanism does not require low-density power-feed plasmas to be swept in from large radius or the development of electrothermal instabilities. The Hall Instability is, thus, a new, independent explanation for the origin of the helical instabilities observed in axially premagnetized liner experiments. Simulation results supporting this mechanism are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of high field plasma dynamics in a laser-produced plasma expanding into a background gas

This paper presents an overview of experimental results of a laser-produced plasma expanding into a background gas, immersed within a large range of highly uniform magnetic fields (of up to 3 T), that are transverse to the expanding plasma. We used intensified gated imaging to capture the expansion of the plasma across and along the magnetic field lines to observe the spatiotemporal expansion dynamics for different magnetic field strengths. We observe changes in the perpendicular and parallel dynamics of the laser-produced plasmas expansion at high magnetic field. In addition, our results have also indicated the presence of electron-ion hybrid instabilities at relatively high pressures (100 mTorr) and relatively high magnetic field strengths (2 T), in accordance with theoretical calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of inlet throttling on thermohydraulic instability in a large scale water-based RCCS: An experimental study

The objective of the present experimental study is to investigate the effect of inlet throttling on the thermohydraulic stability of a large scale water-based Reactor Cavity Cooling System (RCCS). The test was performed using the water-based Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne, which represented a ½ axial scale and 12.5° sector slice of the full scale Framatome 625 MW t SC-HTGR RCCS concept. A two-phase steady state was first established through direct condensate refill, followed by increased inlet throttling over 10 stages, corresponding to a loss coefficient K over the range of 0.05–653. With the inlet throttling gradually increased, the system experienced a unique transition process between stabilization and destabilization. Through a stability analysis, three instability mechanisms were identified in the present test, including a compound mechanism due to both natural circulation oscillations (NCOs) and density wave oscillations (DWOs), Type-II DWOs, and geysering.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Onset of Plasmoid Reconnection during Magnetorotational Instability

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two- and three-dimensional numerical modeling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical (z) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal (rz) plane, driving magnetic reconnection when the Lundquist number S > 3 × 10 2 , making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at S = 5 × 10 3 , self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When nonaxisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets more unstable and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗

Nonlinear ablative Rayleigh-Taylor instability: Increased growth due to self-generated magnetic fields

The growth rate of the nonlinear ablative Rayleigh-Taylor (RT) instability is enhanced by magnetic fields self-generated by the Biermann battery mechanism; a scaling for this effect with perturbation height and wavelength is proposed and validated with extended-magnetohydrodynamic simulations. The magnetic flux generation rate around a single RT spike is found to scale with the spike height. The Hall parameter, which quantifies electron magnetization, is found to be strongly enhanced for short-wavelength spikes due to Nernst compression of the magnetic field at the spike tip. Furthermore, the impact of the magnetic field on spike growth is through both the suppressed thermal conduction into the unstable spike and the Righi-Leduc heat flow deflecting heat from the spike tip to the base. Righi-Leduc is found to be the dominant effect for small Hall parameters, while suppressed thermal conduction dominates for large Hall parameters. These results demonstrate the importance of considering magnetic fields in all perturbed inertial confinement fusion hot spots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced mixing in inertial confinement fusion with early-time interface acceleration

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. In conclusion, this characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Causes of fuel–ablator mix inferred from modeling of monochromatic time-gated radiography of OMEGA cryogenic implosions

Here, we present evidence, in the context of OMEGA cryogenic target implosions, that laser imprint, known to be capable of degrading laser-direct-drive target performance, plays a major role in generating fuel–ablator mix. OMEGA cryogenic target implosions show a performance boundary correlated with acceleration-phase shell stability; for sufficiently low adiabats (where the adiabat is the ratio of the pressure to the Fermi pressure) and high in-flight aspect ratios (IFAR's), the neutron-weighted shell areal density and neutron yield relative to the clean simulated values sharply decline. Direct evidence of Rayleigh–Taylor fuel–ablator mixing was previously obtained using a Si Heα backlighter driven by an ~20-ps short pulse generated by OMEGA EP. The shadow cast by the shell shortly prior to stagnation, as diagnosed using backlit radiographs, shows a softening near the limb, which is evidence of an ablator–fuel mix region for a low-adiabat implosion (α ~ 1.9, IFAR = 14) but not for a moderate adiabat implosion (α ~ 2.5, IFAR = 10). We find good agreement between experimental and synthetic radiographs in simulations that model laser imprint and account for uncertainty in the initial ablator thickness. Furthermore, we further explore the role of other mechanisms such as classical instability growth at the fuel–ablator interface, species concentration diffusion, and long-wavelength drive and target asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding the fusion yield dependencies in OMEGA DT-layered implosion experiments using a physics-based statistical mapping model

Improving the performance of inertial confinement fusion implosions requires physics models that can accurately predict the response to changes in the experimental inputs. Good predictive capability has been demonstrated for the fusion yield using a statistical mapping of simulated outcomes to experimental data [Gopalaswamy et al., Nature 565(771), 581–586 (2019)]. In this paper, a physics-based statistical mapping approach is used to extract and quantify all the major sources of degradation of fusion yield for direct-drive implosions on the OMEGA laser. Here, the yield is found to be dependent on the age of the deuterium tritium fill, the ℓ = 1 asymmetry in the implosion core, the laser beam-to-target size ratio, and parameters related to the hydrodynamic stability. A controlled set of experiments were carried out where only the target fill age was varied while keeping all other parameters constant. The measurements were found to be in excellent agreement with the fill age dependency inferred using the mapping model. In addition, a new implosion design was created, guided by the statistical mapping model by optimizing the trade-offs between increased laser energy coupling at larger target size and the degradations caused by the laser beam-to-target size ratio and hydrodynamic instabilities. When experimentally performed, an increased fusion yield was demonstrated in targets with larger diameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗