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

Laser-induced fluorescence measurements of ion fluctuations in electron and ion presheaths

Electrostatic flow-driven instabilities have been predicted in the presheath region of low-temperature plasmas for both ion and electron rich sheaths. The ion-acoustic fluctuations that arise from these instabilities can influence transport through an enhanced collision rate. In this work, ion fluctuation spectra were measured using laser-induced fluorescence (LIF) in the presheath region for both sheath polarities. Furthermore, the non-invasive nature of LIF compared to that of probes is particularly valuable for these measurements because the sheath and presheath around a probe may generate their own flow-driven instabilities, which obscure the measurement. Measurements were made using a recently implemented field programmable gate array-based system able to measure ion fluctuation spectra up to 1 MHz using a two-point correlation function technique. Narrow bandwidth ion fluctuations were observed near 0:46 f pi (460 kHz) for both sheath polarities, where f pi is the ion plasma frequency. The observed fluctuations were significantly stronger in the electron sheath case and were observed over a larger volume of plasma. Fluctuations were measured in the ion sheath case at locations far enough from the sheath that linear theory predicts stability, suggesting that ion-acoustic instabilities generated in the relatively small unstable region near the boundary reflect from the sheath and perturb a much larger volume of the plasma. This is expected to affect probes with both ion and electron rich sheaths, meaning any dc biased probe may effectively act as an ion-acoustic wave antenna. These measurements are consistent with the recent theory and particle-in-cell results.

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

On the Cutting of Metals: A Mechanics Viewpoint

The mechanics of large-strain deformation in cutting of metals is discussed, primarily from viewpoint of recent developments in in situ analysis of plastic flow and microstructure characterization. It is shown that a broad range of deformation parameters can be accessed in chip formation—strains of 1–10, strain rates of 10–10 5 /s, and temperatures up to 0.7T m —and controlled. This range is far wider than achievable by any other single-stage, severe plastic deformation (SPD) process. The resulting extreme deformation conditions produce a rich variety of microstructures in the chip. Four principal types of chip formation—continuous, shear-localized, segmented, and mushroom-type—as elucidated first by Nakayama (1974, “The Formation of ‘Saw-Toothed Chip’ in Metal Cutting,” Proceedings of International Conference on Production Engineering, Tokyo, pp. 572–577) are utilized to emphasize the diverse plastic flow phenomena, especially unsteady deformation modes that prevail in cutting. These chip types are intimately connected with the underlying flow, each arising from a distinct mode and triggered by an instability phenomenon. The role of plastic flow instabilities such as shear banding, buckling, and fracture in mediating unsteady flow modes is expounded, along with consequences of the flow modes and chip types for the cutting. Sinuous flow is shown to be the reason why gummy (highly strain-hardening) metals, although relatively soft, are so difficult to cut. Synthesizing the various observations, a hypothesis is put forth that it is the stability of flow modes that determines the mechanics of cutting. Furthermore, this leads to a flow-stability phase diagram that could provide a framework for predicting chip types and process attributes.

42 ENGINEERING↗
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 OSTI2020

Different responses of the Rayleigh–Taylor type and resistive drift wave instabilities to the velocity shear

The effects of velocity shear on the unstable modes driven by the effective gravity (Rayleigh–Taylor and interchange) and resistive drift wave instabilities for inhomogeneous equilibrium fluid/plasma density are analyzed for the localized eigenmode problems. It is shown that the effect of the velocity shear drastically depends on the type of instability. Whereas the velocity shear can significantly suppress both Rayleigh–Taylor and interchange instabilities, it has only a weak impact on the growth rate of the resistive drift wave. This is directly related to the physical nature of these instabilities. For the Rayleigh–Taylor and interchange instabilities, the shear flow tilts the eddies of the stream functions, while for the resistive drift wave instability the shear flow simply shifts the eddies in the radial direction with no tilting. However, for a large velocity shear, the eigenmode solutions for resistive drift waves cease to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
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 OSTI2020

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗