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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
DOE OSTI2022

Understanding the role of flow dynamics in thermoacoustic combustion instability

Thermoacoustic combustion instability is one of the most challenging operational issues in several high-performance, low-emissions combustion technologies, including gas turbines, aircraft engines, rockets, and industrial boilers. Driven by the coupling between combustor acoustics and flame heat release rate fluctuations, thermoacoustic combustion instability can lead to reduced operability, increased emissions, and, in the most extreme cases, catastrophic failure of combustor components. The feedback loop between acoustics and combustion is often facilitated by fluid mechanic oscillations, referred to as “velocity coupling,” whereby acoustic oscillations drive flow fluctuations, which in turn create fluctuations in the flame. The character of these fluid mechanic oscillations is highly dependent on the structure of the flow field and the receptivity of the flow to external excitation. Combustor flow fields use features like fluid recirculation and shear to enhance flame holding and reduce emissions, but these are also the same features that can make the flow receptive to acoustic excitation or even drive self-excited oscillations. Here, in this paper, we discuss the basics of thermoacoustic instability with a focus on the role of hydrodynamic oscillations in typical combustor flows. To facilitate this discussion, we explore the hydrodynamic instability characteristics of several key combustor unit flows (wakes, swirling jets, etc.) and show how the hydrodynamic stability of a flow is an important consideration in determining a combustor’s propensity for thermoacoustic oscillations. Several examples of coupling between hydrodynamics and thermoacoustics are discussed to illustrate this important link. The paper concludes by discussing the potential for designing flow fields that are thermoacoustic instability resistant, either through a reduction in the receptivity of the flow or through nonlinear coupling mechanisms by which self-excited flow instabilities can suppress velocity-coupled combustion oscillations.

42 ENGINEERING↗
DOE OSTI2024

Cooling and instabilities in colliding radiative flows with toroidal magnetic fields

ABSTRACT We report on the results of a simulation-based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, is relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative losses is included in the calculation. We systematically vary plasma beta (βm) in the flows, the strength of the cooling (Λ0), and the exponent (α) of temperature dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a ‘spine’ is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depend on initial flow parameters and provide a simple analytical model that captures the basic features of the flow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Transition from wall modes to multimodality in liquid gallium magnetoconvection

Coupled laboratory-numerical experiments of Rayleigh-Bénard convection in liquid gallium subject to a vertical magnetic field are presented. The experiments are carried out in two cylindrical containers with diameter-to-height aspect ratio Γ = 1.0 and 2.0 at varying thermal forcing (Rayleigh numbers 10 5 ≲ Ra ≲ 10 8 ) and magnetic field strength (Chandrasekhar numbers 0 ≲ Ch ≲ 3 × 10 5 ). Laboratory measurements and numerical simulations confirm that magnetoconvection in our finite cylindrical tanks onsets via nondrifting wall-attached modes, in good agreement with asymptotic predictions for a semi-infinite domain. With increasing supercriticality, the experimental and numerical thermal measurements and the numerical velocity data reveal transitions between wall mode states with different azimuthal mode numbers and between wall-dominated convection to wall and interior multimodality. These transitions are also reflected in the heat transfer data, which combined with previous studies connect onset to supercritical turbulent behaviors in liquid metal magnetoconvection over a large parameter space. The gross heat transfer behaviors between magnetoconvection and rotating convection in liquid metals are compared and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The effect of shear flow on the resistive tearing instability

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-Ψ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-Ψ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

astrophysical plasmas↗
DOE OSTI2024

Saturation of Fishbone Instability by Self-Generated Zonal Flows in Tokamak Plasmas

Gyrokinetic simulations of the fishbone instability in DIII-D tokamak plasmas find that self-generated zonal flows can dominate the nonlinear saturation by preventing coherent structures from persisting or drifting in the energetic particle phase space when the mode frequency down-chirps. Results from the simulation with zonal flows agree quantitatively, for the first time, with experimental measurements of the fishbone saturation amplitude and energetic particle transport. Moreover, the fishbone-induced zonal flows are likely responsible for the formation of an internal transport barrier that was observed after fishbone bursts in this DIII-D experiment. Finally, gyrokinetic simulations of a related ITER baseline scenario show that the fishbone induces insignificant energetic particle redistribution and may enable high performance scenarios in ITER burning plasma experiments. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of zonal flows on ion temperature gradient instability in the scrape-off layer of a field-reversed configuration

Gyrokinetic simulations of long wavelength ion temperature gradient (ITG) turbulence in the scrape-off layer (SOL) of a field-reversed configuration (FRC) find that zonal flows are nonlinearly generated and are the dominant mechanism for the nonlinear saturation of the ITG instability. After the ITG saturation, zonal flows remain undamped and gradually suppress the turbulent transport to a very low level. In the simulations with collisions, collisional damping gradually reduces zonal flow amplitude to a lower level, which allows finite ITG turbulence intensity and ion heat transport in the SOL. The steady state turbulence intensity and ion heat transport are found to be proportional to the collision frequency. Furthermore, this favorable scaling suggests that minimizing collisions (e.g. increasing temperature, reducing impurity content, etc) and preserving toroidal symmetry could improve plasma confinement in the FRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hydrodynamic Mechanism for Clumping along the Equatorial Rings of SN1987A and Other Stars

An explanation for the origin and number of clumps along the equatorial ring of Supernova 1987A has eluded decades of research. Our linear analysis and hydrodynamic simulations of the expanding ring prior to the supernova reveal that it is subject to the Crow instability between vortex cores. Here, the dominant wave number is remarkably consistent with the number of clumps, suggesting that the Crow instability stimulates clump formation. Although the present analysis focuses on linear fluid flow, future nonlinear analysis and the incorporation of additional stellar physics may further elucidate the remnant structure and the evolution of the progenitor and other stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2023

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

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 OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2021

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↗