Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TAYLOR 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2020

Unstable Invasion of Sedimenting Granular Suspensions

Here we investigate the development of mobility inversion and fingering when a granular suspension is injected radially between horizontal parallel plates of a cell filled with a miscible fluid. While the suspension spreads uniformly when the suspension and the displaced fluid densities are exactly matched, even a small density difference is found to result in a dense granular front which develops fingers with angular spacing that increase with granular volume fraction and decrease with injection rate. We show that the timescale over which the instability develops is given by the volume fraction dependent settling timescale of the grains in the cell. We then show that the mobility inversion and the nonequilibrium Korteweg surface tension due to granular volume fraction gradients determine the number of fingers at the onset of the instability in these miscible suspensions.

02 PETROLEUM↗
DOE OSTI2024

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Kinetic Kelvin-Helmholtz instability at a finite sized object

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the self-consistent flow around a two-dimensional obstacle with zero intrinsic magnetic field. Plasma outlfow from the obstacle is used to establish a boundary layer between the incoming solar wind and the outgoing plasma. Because the self-consistent flow solution, a velocity shear is naturally set up at this interface, and since the magnetic field for these simulations is transverse to this flow, the Kelvin-Helmholtz (K-H) instability can be excited at low-velocity shear. Simulations demonstrate the existence of the instability even near the subsolar location, which normally is thought to be stable to this instability. The apparent reason for this result is the overall time dependence at the boundary layer, which gives rise to a Rayleigh-Taylor like instability which provides seed perturbations for the K-H instability. These results are directly applicable to Venus, comets, artificial plasma releases, and laser target experiments. This result has potentially important ramifications for the interpretation of observational results as well as for an estimation of the cross-field transport. The results suggest that the K-H instability may play a role in dayside processes and the Venus ionopause, and may exist within the context of more general situations, for example, the Earth's magnetopause.

Thomas, V. A.↗
NASA NTRS1982

The investigation of flow instabilities on a rotating disk with curvature in the radial direction

The major objective is to explore any visible differences of the flow field with wall curvature of the test body, including possible interaction between Taylor-Gortler instabilities present along concave walls and the inflexional instabilities investigated here. An experimental study was conducted with emphasis placed on making visual observations and recording photographically the flow instabilities present under three different rotating bodies: a flat disk, a concave paraboloid, and a convex paraboloid. The data collected for the three test bodies lead to the conclusion that the wall curvature of the concave and convex paraboloids did not alter the observed flow field significantly from that observed on the flat disk.

Intemann, P. A.↗
DOE OSTI2022

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗
DOE OSTI2025

Anisotropic Turbulent Flows Observed in Above-the-loop-top Regions during Solar Flares

Abstract Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in three-dimensional (3D) magnetohydrodynamic (MHD) simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band extreme-ultraviolet images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh–Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the antisunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous 3D flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Xie, Xiaoyan (ORCID:0009000705827807)↗
DOE OSTI2025

Rotation of electrothermal-instability-driven overheating structure due to helically oriented surface magnetic field on a high-current-density aluminum rod

Experiments on the 1-MA, 100-ns-rise-time Mykonos Facility demonstrate rotation of electrothermal instability (ETI)-driven overheating structure on 1.00-mm-diameter, 10-nm-surface roughness, 99.999%-pure aluminum rods, which are pulsed with helically polarized surface magnetic field. Rods are machined to include pairs of 10-micron-scale quasi-hemispherical voids or “engineered defects (ED)” which provide the dominant current density perturbation from which ETI grows most rapidly. Experiments include an axial magnetic field component through the addition of a helically wound return-current electrode or “helical return can (HRC).” For a given HRC design, azimuthal field (B ɵ ) and axial field (B z ) components rise at a prescribed and fixed ratio, driving an increasing magnetic field of constant polarization at the rod's surface; most experiments generated surface magnetic field at a 15-degree field polarization angle (from horizontal) defined as ɸ B = arctan(B z /B ɵ ). ETI-driven emission patterns from individual ED are observed to rotate along ɸ B , while emission patterns from dielectric-coated ED pairs are shown to elongate and preferentially merge along ɸ B , in qualitative agreement with 3D-magnetohydrodynamic simulations. These data strongly support that for a randomized distribution of current density perturbations on a high-current density conductor, nearby perturbations will favorably merge about ɸ B , with the degree of merging increasing with current. Such observations offer fundamental new understanding of the seeding mechanisms of the helical magneto-Rayleigh Taylor (MRT) instabilities observed from axially magnetized magnetically driven imploding liners.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Double cylinder implosion experiments at the National Ignition Facility

Here, cylindrical implosion experiments are used to directly measure instability growth in a convergent geometry, providing a wealth of data for model validation. Double cylinders are a natural extension of the platform and enable measurements at a classically unstable interface, the outer surface of the inner cylinder, which experiences no ablative stabilization from the laser drive. However, the utility of this platform relies upon maintaining adequate axial uniformity of the inner cylinder during the implosion. Although previous smaller-scale double cylinder experiments exhibited acceptable levels of axial uniformity, radiation-hydrodynamics simulations of larger-scale double cylinders predict more axial non-uniformity induced by the impedance mismatch as the shock wraps around the axial ends of the inner cylinder. A mechanism to reduce axial non-uniformity in these larger double cylinder implosions is presented, and preliminary experimental data confirms the efficacy of the selected mitigation approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

On the control of electron heating for optimal laser radiation pressure ion acceleration

We study the onset of electron heating in intense laser–solid interactions and its impact on the spectral quality of radiation pressure accelerated ions in both hole boring and light sail regimes. Two- and three-dimensional particle-in-cell (PIC) simulations are performed over a wide range of laser and target parameters and reveal how the pulse duration, profile, polarization and target surface stability control the electron heating, the dominant ion acceleration mechanisms and the ion spectra. We find that the onset of strong electron heating is associated with the growth of the Rayleigh–Taylor-like instability at the front surface and must be controlled to produce high-quality ion beams, even when circularly polarized lasers are employed. We define a threshold condition for the maximum duration of the laser pulse that allows mitigation of electron heating and radiation pressure acceleration of narrow energy spread ion beams. The model is validated by three-dimensional PIC simulations, and the few experimental studies that reported low energy spread radiation pressure accelerated ion beams appear to meet the derived criteria. The understanding provided by our work will be important in guiding future experimental developments, for example for the ultrashort laser pulses becoming available at state-of-the-art laser facilities, for which we predict that proton beams with $\sim$ 150–250 MeV, $\sim$ 30% energy spread, and a total laser-to-proton conversion efficiency of $\sim$ 20% can be produced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Controlling morphology and improving reproducibility of magnetized liner inertial fusion experiments

X-ray imaging indicates magnetized liner inertial fusion (MagLIF) stagnation columns have a complicated quasi-helical structure with significant variations in x-ray brightness along the column. In this work, we describe MagLIF experiments aimed at controlling these stagnation structures by varying the initial liner geometry and composition. First, by varying the initial aspect ratio of the liner, we demonstrate a change in the stagnation structures that is consistent with helical magneto Rayleigh–Taylor (MRT) instabilities feedthrough from the outer-to-inner surfaces of the liner. Second, to minimize the seed for such instabilities, we incorporate a dielectric coating on the outer surface of the beryllium liner, which has previously been shown to reduce the growth of the electrothermal instability, a likely seed for MRT growth. Using this coating, we achieve a stagnation column with significantly reduced helical structure and axial variation in x-ray brightness. We discuss how this coating changes the evolution of structures through stagnation along with the spatial uniformity of neutron production. Finally, we show that these more uniform stagnations also result in improved reproducibility in stagnation temperatures and primary DD neutron yield.

Ampleford, D. J. (ORCID:0000000237422421)↗
DOE OSTI2026

Staged Z-pinch experiment at the double-EAGLE pulsed power facility

The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.

Physics - Plasma physics↗
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

Azimuthal magnetic field distribution in gas-puff Z -pinch implosions with and without external magnetic stabilization

An experimental study of the magnetic field distribution in gas-puff Z pinches with and without a preembedded axial magnetic field (B z0 ) is presented. Spatially resolved, time-gated spectroscopic measurements were made at the Weizmann Institute of Science on a 300 kA, 1.6 μs rise time pulsed-power driver. The radial distribution of the azimuthal magnetic field, B θ , during the implosion, with and without a preembedded axial magnetic field of B z0 = 0.26 T, was measured using Zeeman polarization spectroscopy. The spectroscopic measurements of B θ were consistent with the corresponding values of B θ inferred from current measurements made with a B-dot probe. One-dimensional magnetohydrodynamic simulations, performed with the code TRAC-II, showed agreement with the experimentally measured implosion trajectory, and qualitatively reproduced the experimentally measured radial B θ profiles during the implosion when B z0 = 0.26 T was applied. Simulation results of the radial profile of B θ without a preembedded axial magnetic field did not qualitatively match experimental results due to magneto-Rayleigh-Taylor (MRT) instabilities. Furthermore, our analysis emphasizes the importance of MRT instability mitigation when studying the magnetic field and current distributions in Z pinches. Discrepancies of the simulation results with experiment are discussed.

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↗