SEARCH · Search NASA

Results for “Rayleigh-Taylor 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Mitigation of magneto-Rayleigh-Taylor instability growth in a triple-nozzle, neutron-producing gas-puff Z pinch

The gas-puff Z-pinch is a well-known source of x-rays and/or neutrons, but it is highly susceptible to the magneto-Rayleigh-Taylor instability (MRTI). Approaches to MRTI mitigation include density profile tailoring, in which nozzles are added or modified to alter the acceleration trajectory, and axial pre-magnetization, in which perturbations are smoothed out via magnetic field line tension. In this work, we present two-dimensional magnetohydrodynamic simulations of loads driven by an 850 kA,160 ns driver that suggest these mitigation strategies can be additive. The initial axial magnetic field, B z0 , to stabilize a 2.5-cm-radius Ne gas liner imploding onto an on-axis deuterium target can be reduced from 0.7 T to 0.3 T by adding a second liner with a radius of 1.25 cm. Because MRTI mitigation tends to increasingly lower yield with higher B z0 , the use of a lower field is advantageous. Here, we predict a reduction in yield penalty from > 100× with the single liner to < 10× with a double liner. A premagnetized, triple nozzle gas puff could therefore be an attractive source for intense neutrons or other fusion applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Development and assessment of models for turbulent Rayleigh-Taylor mixing using the macroscopic forcing method

Reynolds-Averaged Navier Stokes (RANS) simulations are a popular method for designing ICF experiments, and accurate mixing models are crucial for these simulations to give good predictions. To this end, the present work seeks to demonstrate the Macroscopic Forcing Method (MFM) as a tool for both improving existing RANS models as well as assessing RANS model forms. First, MFM analysis from Lavacot et al. (Phys. Rev. Fluids, 2025) is used to develop the k–L–F model, an extension of the k–L model of Dimonte and Tipton (Phys. Fluids, 2006) that incorporates nonlocality through addition of a turbulent species flux transport equation. MFM is then applied to the k–L–F model along with the k–L and BHR–4 models to assess their forms and compare the model-implied eddy diffusivity moments to those measured from high-fidelity simulations. Furthermore, the analysis reveals that models incorporating nonlocality (k–L–F and BHR–4) match the high-fidelity simulation data better than purely local models (k–L), both in terms of mean fields and eddy diffusivity moments. However, all of the considered RANS models struggle to match temporal moments at high Atwood numbers, highlighting the importance of temporal nonlocality in these regimes and the need for additional improvement even among models incorporating nonlocality.

general physics↗

DOE OSTI2026

Dynamic strength of iron under pressure-temperature conditions of Earth’s inner core

Iron (Fe) is a primary constituent of terrestrial planetary cores, yet its rheological properties under extreme conditions remain uncertain. Here we present direct measurements of Fe strength at 310-430 GPa pressures and 3700-5800 K temperatures, obtained using Rayleigh-Taylor (RT) instability experiments at the National Ignition Facility. Single-crystal α-Fe samples with [001] and [111] orientations are shock-ramp compressed past the α-ε transition along paths approaching Earth’s inner core conditions. We find that ε-Fe derived from [001] α-Fe is consistently stronger (11-20 GPa) than that from [111] α-Fe (8-18 GPa), contrary to the trend at ambient conditions. Large-scale molecular dynamics simulations reproduce this atypical strength anisotropy and attribute it to microstructural evolution during the phase transition and subsequent ε-phase plasticity. Ripple growth analysis further constrains viscosities of 100-170 Pa·s under the driven conditions. Our results provide experimental benchmarks for Fe rheology at inner-core conditions, with implications for seismic anisotropy and the geodynamo.

Condensed-matter physics↗

DOE OSTI2025

An experimental and computational study of thin-layer Rayleigh–Taylor instability development during deceleration with and without an externally applied magnetic field

The importance of mitigating the Rayleigh–Taylor instability (RTI) in inertial confinement fusion (ICF) is critical to successfully achieve high gain fusion yield. Consequently, understanding the seed mechanisms of RTI and the potential evolution of RTI in ICF relevant conditions is crucial. Single feature perturbations consistently demonstrate non-linear RTI evolution, for which an experimental platform on OMEGA-EP is developed. Manufacturing defects introduced into the target design require exploration of unanticipated changes to RTI development and an identification of targets that will still render quantifiable physics results. Consequently, it is presented that the inherent 3D nature of experimental targets necessitates 3D modeling for accurate design work and predictive modeling of experimental targets, especially when high resolution imaging diagnostics, like Fresnel Zone Plates, are utilized. A study of the morphology of the RTI evolution due to changing initial conditions and the presence of an externally applied magnetic field are also explored. Experimental data show thin-layer RTI morphology comparable to resistive magneto hydrodynamic 3D results. A discussion on the impacts of an externally applied magnetic field makes the case for continued efforts to observe a magnetic field's impact on RTI morphology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Design of a high-resolution Rayleigh-Taylor experiment with the Crystal Backlighter Imager on the National Ignition Facility

The Rayleigh-Taylor (RT) instability affects a vast range of High Energy Density (HED) length scales, spanning from supernova explosions (10 13 m) to inertial confinement fusion (10 -6 m). In inertial confinement fusion, the RT instability is known to induce mixing or turbulent transition, which in turn cools the hot spot and hinders ignition. The fine-scale features of the RT instability, which are difficult to image in HED physics, may help determine if the system is mixing or is transitioning to turbulence. Earlier diagnostics lacked the spatial and temporal resolution necessary to diagnose the dynamics that occur along the RT structure. A recently developed diagnostic, the Crystal Backlighter Imager (CBI), [1,2] can now produce an x-ray radiograph capable of resolving the fine-scale features expected in these RT unstable systems. This paper describes an experimental design that adapts a well-characterized National Ignition Facility (NIF) platform to accommodate the CBI diagnostic. Simulations and synthetic radiographs highlight the resolution capabilities of the CBI in comparison to previous diagnostics. Furthermore, the improved resolution of the system can provide new observations to study the RT instability's involvement in mixing and the transition to turbulence in the HED regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Local wave-number model for inhomogeneous two-fluid mixing

We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Sensitivity of the dynamic-shell target to laser drive nonuniformities

The dynamic-shell concept for inertial confinement fusion (ICF) uses an initially homogeneous target and a carefully shaped laser pulse to form a shell and implode it. The laser pulse consists of a series of pickets that drive shocks into the target. The first few shocks converge inwards and rebound from the center of the target, creating an expanding, low-density plasma. Subsequent shocks are launched into the expanding plasma and eventually coalesce to form a shell, which is then imploded with a traditional ICF laser pulse. This study describes radiation-hydrodynamic simulations that investigate the sensitivity of dynamic-shell targets to imperfections in the laser drive. A one-dimensional (1D) study looks at mistiming and power variations in the pickets and a two-dimensional (2D) study examines irradiation perturbations imposed by the laser-beam geometry. Simulations show that less than ~2% power imbalance or 200 ps timing variation in the pickets is sufficient to keep the yield above 90% of the maximum. Additionally, the 2D simulations show that 72 or more beams are required to keep irradiation nonuniformities low enough to obtain fusion yields close to that of 1D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Measuring scale-dependent shape anisotropy by coarse-graining: Application to inhomogeneous Rayleigh-Taylor turbulence

Here, we generalize the “filtering spectrum” to probe scales along different directions by spatial coarse-graining. This multidimensional filtering spectrum quantifies the spectral content of flows that are not necessarily homogeneous. From multidimensional spectral information, we propose a simple metric for shape anisotropy at various scales. The method is applied to simulations of two-dimensional (2D) and 3D Rayleigh-Taylor (RT) turbulence, which is inhomogeneous and anisotropic. We show that 3D RT has clear shape anisotropy at large scales with approximately 4:3 vertical to horizontal aspect ratio but tends toward isotropy at small scales as expected. In sharp contrast, we find that RT in 2D simulations, which are still the main modeling framework for many applications, is isotropic at large scales and its shape anisotropy increases at smaller scales where structures tend to be horizontally elongated. While this may be surprising, it is consistent with recent results that large-scale isotropy in 2D RT is due to the generation of a large-scale overturning circulation via an upscale cascade, while small-scale anisotropy is due to the stable stratification resultant from such overturning and the inefficient mixing in 2D.

42 ENGINEERING↗

DOE OSTI2025

Atwood effects on nonlocality of the scalar transport closure in Rayleigh-Taylor mixing

The importance of nonlocality is assessed in modeling mean scalar transport for turbulent Rayleigh-Taylor (RT) mixing at different Atwood numbers. Building on the two-dimensional incompressible work of Lavacot et al. [J. Fluid Mech. 985, A47 (2024)], the present work extends the macroscopic forcing method to variable density problems in three-dimensional space to measure moments of the generalized eddy diffusivity kernel in RT mixing for increasing Atwood numbers (𝐴 = 0.05, 0.3, 0.5, 0.8). It is found that as 𝐴 increases, (1) the eddy diffusivity moments become asymmetric and (2) the higher-order eddy diffusivity moments become larger relative to the leading-order diffusivity, indicating that nonlocality becomes more important at higher 𝐴. There is a particularly strong temporal nonlocality at higher 𝐴, suggesting stronger history effects. In conclusion, the implications of these findings for closure modeling for finite-Atwood RT are discussed.

general physics↗

DOE OSTI2024

Determining Exact RANS Operators with the Macroscopic Forcing Method (Final Report)

This report contains a compilation of key results and findings from the ACT project “Determining Exact RANS Operators with the Macroscopic Forcing Method.” The Macroscopic Forcing Method (MFM), a numerical tool for determining closure operators, is used to measure eddy diffusivity moments in Rayleigh-Taylor (RT) instability. It is first applied to low-Atwood 2D RT instability; that work is then extended to 3D RT at different finite Atwood numbers. It is found that nonlocality is important for modeling the mean scalar transport closure operator in RT mixing. Additionally, work is done to improve the statistical convergence of MFM for chaotic problems like RT mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Suppression of Richtmyer-Meshkov Instability via Special Pairs of Shocks and Phase Transitions

The classical Richtmyer-Meshkov instability (RMI) is a hydrodynamic instability characterizing the evolution of an interface following shock loading. In contrast to other hydrodynamic instabilities such as Rayleigh-Taylor, it is known for being unconditionally unstable: regardless of the direction of shock passage, any deviations from a flat interface will be amplified. In this article, we show that for negative Atwood numbers, there exist special sequences of shocks which result in a nearly perfectly suppressed instability growth. Here, we demonstrate this principle computationally and experimentally with stepped fliers and phase transition materials. A fascinating immediate corollary is that in specific instances, a phase-transitioning material may self-suppress RMI.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Onset of Rayleigh-Bénard convection at fuel-ablator interfaces in inertial confinement fusion simulations

Here, we demonstrate the onset of Rayleigh-Bénard convection (RBC) near the fuel-ablator interface of inertial confinement fusion (ICF) implosion simulations for the scenario in which Rayleigh-Taylor or Richtmyer-Meshkov instabilities are suppressed. This convective heating within the cold fuel can degrade the adiabat and therefore performance of ICF implosions even without mixing of fuel and ablator. Our work suggests additional considerations to realize design adiabats in otherwise ideal ICF implosions: first, via control of the temperature differences between ablator and fuel and, second, control of the fuel acceleration which we show has a complex interplay with the RBC dynamics.

Physics - Plasma physics↗

DOE OSTI2022

Deceleration-stage Rayleigh–Taylor growth in a background magnetic field studied in cylindrical and Cartesian geometries

Experiments have identified the Rayleigh–Taylor (RT) instability as one of the greatest obstacles to achieving inertial confinement fusion. Consequently, mitigation strategies to reduce RT growth and fuel–ablator mixing in the hotspot during the deceleration phase of the implosion are of great interest. In this work, the effect of seed magnetic fields on deceleration-phase RT growth are studied in planar and cylindrical geometries under conditions relevant to the National Ignition Facility (NIF) and Omega experiments. The magnetohydrodynamic (MHD) and resistive-MHD capabilities of the FLASH code are used to model imploding cylinders and planar blast-wave-driven targets. Realistic target and laser parameters are presented that suggest the occurrence of morphological differences in late-time RT evolution in the cylindrical NIF case and a measurable difference in spike height of single-mode growth in the planar NIF case. The results of this study indicate the need for target designs to utilize an RT-unstable foam–foam interface in order to achieve sufficient magnetic field amplification to alter RT evolution. Benchmarked FLASH simulations are used to study these magnetic field effects in both resistive and ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Post-Shot Report for OMEGA Double Cylinders (CylDRT 22B)

The direct-drive double cylinder experimental platform is a high-energy-density (HED) science platform designed to image an imploding cylindrical target. The target consists of a directly-driven outer cylinder and a shock-and-collision driven inner cylinder. The purpose of this platform is to study hydrodynamic instability growth on the inner cylinder, the outer surface of which is classically Rayleigh-Taylor unstable during the acceleration phase. We present results from recent experiments at the OMEGA laser facility. These experiments were designed as a proof-of-principle for the platform, using the same cylinder exterior dimensions and direct-drive beam configuration as previous single cylinder experiments. In these experiments, three sets of targets were fielded: no machined perturbations (smooth), a sinusoidal mode-10 perturbation on the outer surface of the inner cylinder, and a mode-20 perturbation on the outer surface of the inner cylinder. The primary diagnostic was a gated x-ray framing camera which imaged the backlit inner cylinder on-axis for sixteen frames over a time window of 1 ns for each shot. A second side-lit radiograph captured one image per shot, diagnosing axial uniformity. In this report we include an overview of the results from both the backlighter and the sidelighter diagnostics. We discuss at length the experimental analysis process. Finally, we present the results of the smooth target implosion trajectory and compare them to post-shot simulations. We see favorable agreement between simulation and experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Magnetic-Field Effects in Unstable High-Energy-Density Plasmas

The influence of magnetic fields in high-energy-density (HED) systems is a widely unexplored field of research that can play an important role in magnetized inertial confinement fusion (ICF) and in astrophysical systems. Hydrodynamic instabilities are presently believed to be the dominant degradation mechanism in ICF implosions. Therefore, understanding the influence of an imposed B-field on RT growth is important for future mitigation strategies as well as fundamental physics of ICF implosions that utilize B-fields for enhanced performance. The Crab Nebula, one of the most observed objects in our universe, contains elongated spikes of material created by hydrodynamic instabilities, however these spikes do not break-up and become turbulent as they evolve late in time. One hypothesis for this behavior is that B-fields around these spikes prevent vortex generation typical in these types of systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Amplification of laser imprint in the presence of strong, externally imposed, target-normal magnetic fields

Here, an experiment was performed on the OMEGA EP Laser System to investigate the effects of magnetic fields on laser imprint. A 30 µ⁢m thick CH target was driven by a single beam delivering 3.0 kJ of UV energy in a 5 ns square pulse without smoothing by spectral dispersion. The Rayleigh-Taylor amplification of this beam's imprint on the target surface was monitored using face-on gated x-ray radiography from a Gd backlighter. Magnetic fields of up to 45 T were applied normal to the target surface. Analysis of the resulting radiographs shows a 60 ± 13% increase in the spectrally resolved surface perturbation amplitudes, consistent at all times and for all unsaturated frequencies. This consistency indicates that the increase in perturbation amplitudes was due to a change in the initial amplitudes rather than Rayleigh-Taylor growth. This is supported by the trajectory of individual modes and the inferred time-averaged perturbation growth rates. Laser imprint is therefore inferred to have increased due to strong magnetic fields that remain in the conduction zone of the target, suppressing off-axis electron motion and limiting the effects of thermal smoothing.

direct drive↗