Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow 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 109 records · Page 6

Theory of the Tertiary Instability and the Dimits Shift from Reduced Drift-Wave Models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $γ$ TI is derived within the modified Hasegawa-Wakatani model. We show that $γ$TI equals the primary-instability growth rate plus a term that depends on the local U''; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry-Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear |U'| is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Drift waves enstrophy, zonal flow, and nonlinear evolution of the modulational instability

The interaction of the drift wave (DW) turbulence and zonal flow (ZF) is investigated with the modified Hasegawa–Mima equation taking into account the backreaction of ZF velocity on DW turbulence. It is shown that the y-averaged enstrophy of DW turbulence and the velocity of ZF are intrinsically related. By utilizing this feature, a nonlinear stage of DW modulational instability is considered within the framework of the wave kinetic equation. It is shown that in this approximation, the nonlinear stage of the modulational instability results in the collapsing solutions, accompanied by the “wave breaking” phenomenon. Numerical simulations based on the Hasegawa–Mima equation show that for a weak DW turbulence, Φ ̃ = ( e φ ̃ / T e ) ( L n / ρ s ) ⪝ 1 , the collapsing-like features on both ZF and y-averaged enstrophy of DW turbulence decay in time and then re-emerge again at different locations. For the case of a strong DW turbulence, Φ ̃ > 1 , where nonlinear interactions of DW harmonics dominate, stable spatial structures of ZF and y-averaged enstrophy of DW turbulence emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Operando synchrotron X-ray analysis of melt pool dynamics in an Al-Sn immiscible alloy

The melt flow in an Al-50vol.% Sn immiscible alloy, produced by single-track laser melting of Al and Sn elemental powders, was studied in real time. High-speed synchrotron X-ray imaging was used to track Al and Sn liquids’ movements, and also to examine elemental distributions in the laser tracks, complimented by electron microscopy after solidification. Key aspects, including melt pool geometry, keyhole instability, and flow dynamics (flow pattern and velocity), were examined using digital image analysis. Relatively deeper melt pools formed at 400 W and 300 mm/s exhibited greater stability, with smooth surfaces, consistent outward flow, and minor vortices near the keyhole. In contrast, shallower pools produced at higher scanning speeds (>500 mm/s) demonstrated greater instability with increased surface waviness, and stronger velocity fluctuations, leading to numerous micro-vortices and increased Al-Sn heterogeneity. Velocity scale estimations, supported by experimental observations, examined the roles of vapour pressure, Marangoni effect, buoyancy, inertial, and surface tension forces in the flow. The results revealed that vapour pressure and mechanical waves dominated at high scanning speeds (shallow pools), while Marangoni forces were equally significant in deep pools at lower speeds (300 mm/s). Buoyancy was found to have minimal impact in both cases. Furthermore, the interaction between inertial and surface tension forces played a critical role in determining the degree of waviness of the pools’ surfaces. These findings offer valuable insights into melt pool dynamics during laser processing of immiscible alloys and other metallic systems using elemental powders, and provide guidance for developing high-fidelity computational fluid dynamics models.

immiscible alloys↗

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↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

Additive Manufacturing Interface Profiles for Jet Studies on the Vertical Shock Tube

Understanding the mixing between two fluids of disparate densities (i.e., variable density mixing) is important for the lab's mission. Experiments are required to validate our codes, such as RAGE, FLAG, and the BHR mix model. The role initial conditions (ICs) have on instability growth and the transition to turbulence is not completely understood and comparisons with experiments are needed. The Extreme Fluids team focuses on experiments at Los Alamos that utilize gases and low shock velocities (Mach <2 in air). These experiments have the advantage of high fidelity diagnostics capable of temporally and spatially resolved velocity and density measurements. One main difficulty is setting up a specific interface profile between the two gases. While attempts have been made to do so without the aid of a solid structure (i.e., a membrane), the results are limited in the range of possible shapes and reproducibility. Membranes provide the best way to provide a deterministic initial condition. However, fragments from the membrane can interfere with the flow and instability growth, as well as the diagnostics. In this project, we focused on exploring the materials and techniques available via modern additive manufacturing (i.e., 3D printing) to test if they could address some of the shortcomings of prior membrane designs. Our goal was to find a material that would have minimal impact on our experiment while also providing us with a means to produce complex initial conditions. This report details our work to explore the range of additive manufacturing capabilities and the jet-like structures we were able to generate using three-dimensional initial conditions.

36 MATERIALS SCIENCE↗