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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1991

MPD work at MIT

MPD work at MIT is presented in the form of the view-graphs. The following subject areas are covered: the MIT program, its goals, achievements, and roadblocks; quasi one-dimensional modeling; two-dimensional modeling - transport effects and Hall effect; microscopic instabilities in MPD flows and modified two stream instability; electrothermal stability theory; separation of onset and anode depletion; exit plane spectroscopic measurements; phenomena of onset as performance limiter; explanations of onset; geometry effects on onset; onset at full ionization and its consequences; relationship to anode depletion; summary on self-field MPD; applied field MPD - the logical growth path; the case for AF; the challenges of AF MPD; and recommendations.

Martinez-Sanchez, Manuel↗
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 OSTI2016

Numerical investigation of initial condition effects on Rayleigh-Taylor instability with acceleration reversals

The influence of initial conditions on miscible incompressible baroclinically driven Rayleigh-Taylor instability undergoing nonuniform acceleration is explored computationally using an implicit large eddy simulation (ILES) technique. We consider here the particular case of evolution during multiple reversals of acceleration direction, where the flow is alternately statically stable or unstable. In the unstable phase, the flow is driven by the baroclinic release of potential energy, whereas in the stable phase, work is done against the density stratification with the energy exchange taking place by wavelike mechanisms. These dynamics are fundamentally different; here, we track the evolution of volume-averaged turbulent statistics that are most sensitive to changes in the distribution of spectral power and bandwidth of the initial conditions as the flow alternates between dynamical regimes due to acceleration reversal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Buoyancy instability of homologous implosions

With this study, I consider the hydrodynamic stability of imploding ideal gases as an idealized model for inertial confinement fusion capsules, sonoluminescent bubbles and the gravitational collapse of astrophysical gases. For oblate modes (short-wavelength incompressive modes elongated in the direction of the mean flow), a second-order ordinary differential equation is derived that can be used to assess the stability of any time-dependent flow with planar, cylindrical or spherical symmetry. Upon further restricting the analysis to homologous flows, it is shown that a monatomic gas is governed by the Schwarzschild criterion for buoyant stability. Under buoyantly unstable conditions, both entropy and vorticity fluctuations experience power-law growth in time, with a growth rate that depends upon mean flow gradients and, in the absence of dissipative effects, is independent of mode number. If the flow accelerates throughout the implosion, oblate modes amplify by a factor (2C) |N 0 |t i , where C is the convergence ratio of the implosion, N 0 is the initial buoyancy frequency and t i is the implosion time scale. If, instead, the implosion consists of a coasting phase followed by stagnation, oblate modes amplify by a factor exp(π|N 0 |t s ), where N 0 is the buoyancy frequency at stagnation and t s is the stagnation time scale. Even under stable conditions, vorticity fluctuations grow due to the conservation of angular momentum as the gas is compressed. For non-monatomic gases, this additional growth due to compression results in weak oscillatory growth under conditions that would otherwise be buoyantly stable; this over-stability is consistent with the conservation of wave action in the fluid frame. The above analytical results are verified by evolving the complete set of linear equations as an initial value problem, and it is demonstrated that oblate modes are the fastest-growing modes and that high mode numbers are required to reach this limit (Legendre mode ℓ ≳ 100 for spherical flows). Finally, comparisons are made with a Lagrangian hydrodynamics code, and it is found that a numerical resolution of ~30 zones per wavelength is required to capture these solutions accurately. This translates to an angular resolution of ~(12/ℓ)°, or ≲ 0.1° to resolve the fastest-growing modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Inviscid linear stability analysis of two vertical columns of different densities in a gravitational acceleration field

Here, we study the inviscid linear stability of a vertical interface separating two fluids of different densities and subject to a gravitational acceleration field parallel to the interface. In this arrangement, the two free streams are constantly accelerated, which means that the linear stability analysis is not amenable to Fourier or Laplace solution in time. Instead, we derive the equations analytically by the initial-value problem method and express the solution in terms of the well-known parabolic cylinder function. The results, which can be classified as an accelerating Kelvin–Helmholtz configuration, show that even in the presence of surface tension, the interface is unconditionally unstable at all wavemodes. This is a consequence of the ever increasing momentum of the free streams, as gravity accelerates them indefinitely. The instability can be shown to grow as the exponential of a quadratic function of time.

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

Effects of residual kinetic energy on yield degradation and ion temperature asymmetries in inertial confinement fusion implosions

The study of Rayleigh–Taylor instability in the deceleration phase of inertial confinement fusion implosions is carried out using the three-dimensional (3-D) radiation-hydrodynamic Eulerian parallel code DEC3D. In this paper, we show that the yield-over-clean is a strong function of the residual kinetic energy (RKE) for low modes. Our analytical models indicate that the behavior of larger hot-spot volumes observed in low modes and the consequential pressure degradation can be explained in terms of increasing the RKE. These results are derived using a simple adiabatic implosion model of the deceleration phase as well as through an extensive set of 3-D single-mode simulations using the code DEC3D. The effect of the bulk velocity broadening on ion temperature asymmetries is analyzed for different mode numbers ℓ = 1 -12. The jet observed in low mode ℓ = 1 is shown to cause the largest ion temperature variation in the mode spectrum. Finally, the vortices of high modes within the cold bubbles are shown to cause lower ion temperature variations than low modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Development of tearing instability in a current sheet forming by sheared incompressible flow

Sweet–Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal magnetohydrodynamics, a coordinate transformation reveals that the flow compresses and stretches perturbations. Here, a multiple scale formulation describes how linear tearing mode theory can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Thermosolutal convection during directional solidification

During solidification of a binary alloy at constant velocity vertically upward, thermosolutal convection can occur if the solute rejected at the crystal-melt interface decreases the density of the melt. We assume that the crystal-melt interface remains planar and that the flow field is periodic in the horizontal direction. The time-dependent nonlinear differential equations for fluid flow, concentration, and temperature are solved numerically in two spatial dimensions for small Prandtl numbers and moderately large Schmidt numbers. For slow solidification velocities, the thermal field has an important stabilizing influence: near the onset of instability the flow is confined to the vicinity of the crystal-melt interface. Further, for slow velocities, as the concentration increases, the horizontal wavelength of the flow decreases rapidly - a phenomenon also indicated by linear stability analysis. The lateral inhomogeneity in solute concentration due to convection is obtained from the calculations. For a narrow range of solutal Rayleigh numbers and wavelengths, the flow is periodic in time.

Mcfadden, G. B.↗
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 OSTI2019

Constitutive Behavior and Hot Workability of a Hot Isostatic Pressed Ti-22Al-25Nb Alloy during Hot Compression

A Ti-22Al-25Nb alloy was fabricated from prealloyed powders by hot isostatic pressing for 2 h at a temperature of 1050 °C and pressure of 100 MPa. The hot deformation behavior of the Ti-22Al-25Nb alloy was characterized by isothermal compression testing at deformation temperatures between 900 and 1060 °C and strain rates between 0.001 and 1 s{sup −1}. Based on the true stress–strain curves, a constitutive equation was constructed to describe the flow stress as a function of the strain rate and deformation temperature. Three-dimensional (3D) processing maps were developed based on dynamic material model theory using the stress flow data to identify the instability and optimization regions of the hot processing parameters. The results show that the flow stress decreases with increasing deformation temperature and decreasing strain rate, and the softening mechanisms are different under different deformation conditions. The apparent activation energies in the (α{sub 2} + β/B2 + O) and (α{sub 2} + B2) phase regions are calculated to be 865.177 and 590.661 kJ/mol, respectively, which suggests that the hot isostatic pressed (HIPed) Ti-22Al-25Nb alloy requires a high hot deformation activation energy. Combined with microstructural observations, the optimal processing domains are determined to be a temperature range of 1030-1060 °C and strain rate range of 0.01-0.1 s{sup −1} in the α{sub 2} + B2 phase region. Moreover, the results indicate that adiabatic shear bands and severe inhomogeneous deformation cause flow instability at lower temperatures (900-1005 °C) and higher strain rates (> 0.1 s{sup −1}).

36 MATERIALS SCIENCE↗
NASA NTRS1990

A review of instability and noise propagation in supersonic flows

Originally analytical and numerical models were to be developed for noise production in supersonic jets, wakes and free shear layers. While the effort was concentrated initially on these aspects, other topics were also pursued, most were of interest to the Jet Noise Group of the Aeroacoustics Branch. An overview is given of subjects reviewed and the investigations that were carried out. A significant effort was devoted to numerically predicting the flow field of a turbulent supersonic wall jet. This information is necessary for computing the pressure in the far field. The wall jet was selected because it represents a generic flow that can be associated with plug nozzle in supersonic engines. It combines the characteristic of a boundary layer with that of a free shear flow. The spatially evolving flow obtained using Dash's code would form the input for the stability analysis program. This analysis would determine the large scale instability wave within the flow. The far field pressure can be computed from the shape of the evolving large scale structure by asymptotic methods. Flow characteristics obtained from a program that analyses the turbulent downstream supersonic flow in a nozzle are described and compared with experimental results.

Daudpota, Q. Isa↗
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↗
DOE OSTI2024

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

RINGED ACCRETION DISKS: INSTABILITIES

We analyze the possibility that several instability points may be formed, due to the Paczyński mechanism of violation of mechanical equilibrium, in the orbiting matter around a supermassive Kerr black hole. We consider a recently proposed model of a ringed accretion disk, made up by several tori (rings) that can be corotating or counter-rotating relative to the Kerr attractor due to the history of the accretion process. Each torus is governed by the general relativistic hydrodynamic Boyer condition of equilibrium configurations of rotating perfect fluids. We prove that the number of the instability points is generally limited and depends on the dimensionless spin of the rotating attractor.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

A Note on the Wave Action Density of a Viscous Instability Mode on a Laminar Free-shear Flow

Using the assumptions of an incompressible and viscous flow at large Reynolds number, we derive the evolution equation for the wave action density of an instability wave traveling on top of a laminar free-shear flow. The instability is considered to be viscous; the purpose of the present work is to include the cumulative effect of the (locally) small viscous correction to the wave, over length and time scales on which the underlying base flow appears inhomogeneous owing to its viscous diffusion. As such, we generalize our previous work for inviscid waves. This generalization appears as an additional (but usually non-negligible) term in the equation for the wave action. The basic structure of the equation remains unaltered.

Balsa, Thomas F.↗