Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid 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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Ionization effects in three-dimensional solar granulation simulations

These numerical studies show that ionization influences both the transport and dynamical properties of compressible convection near the surface of the Sun. About two-thirds of the enthalpy transported by convective motions in the region of partial hydrogen ionization is carried as latent heat. The role of fast downflow plumes in total convective transport is substantially elevated by this contribution. Instability of the thermal boundary layer is strongly enhanced by temperature sensitive variations in the radiative properties of the fluid, and this provides a mechanism for plume initiation and cell fragmentation in the surface layers. As the plumes descend, temperature fluctuations and associated buoyancy forces are maintained because of the increased specific heat of the partially ionized material. This can result is supersonic vertical flows. At greater depths, ionization effects diminish, and the plumes are decelerated by significant entrainment of surrounding fluid.

Rast, Mark P.↗
NASA NTRS2001

Fundamental Processes of Atomization in Fluid-Fluid Flows

This report outlines the major results of the grant "Fundamental Processes of Atomization in Fluid-Fluid Flows." These include: 1) the demonstration that atomization in liquid/liquid shear flow is driven by a viscous shear instability that triggers the formation of a long thin sheet; 2) discovery of a new mode of interfacial instability for oscillatory two-layer systems whereby a mode that originates within the less viscous liquid phase causes interfacial deformation as the oscillation proceeds; 3) the demonstration that rivulet formation from gravity front occurs because the local front shape specified by gravity and surface tension changes from a nose to a wedge geometry, thus triggering a large increase in viscous resistance; and 4) extension of the studies on nonlinear wave evolution on falling films and in stratified flow, particularly the evolution towards large-amplitude solitary waves that tend to generate drops.

McCready, M. J.↗
NASA NTRS2001

Absolute And Convective Instability and Splitting of a Liquid Jet at Microgravity

The objective is to establish a definitive role of the capillary, viscous, and inertial forces at a liquid-gas interface in the absence of gravity by using the fluid dynamics problem of the stability of a liquid jet as a vehicle. The objective is achieved by reexamining known theories and new theories that can be verified completely only in microgravity. The experiments performed in the microgravity facility at NASA Glenn Research Center enable the verification of the theory with experimental data. Of particular interest are (1) to capture for the first time the image of absolute instability, (2) to elucidate the fundamental difference in the physical mechanism of the drop and spray formation from a liquid jet, and (3) to find the origin of the newly discovered phenomenon of jet splitting on earth and in space.

Lin, S. P.↗
NASA NTRS1999

The Effect of a Rotating Magnetic Field on Rayleigh-Benard Convection

A destabilizing vertical temperature gradient and a rotating magnetic field have been applied to a cylindrical column of liquid gallium. The convective flows which arise as a function of these driving parameters are identified. For small magnetic field strengths, a region of stationary flow is observed. This region is bounded by critical values of the Rayleigh and Hartmann numbers. As the rotating magnetic field is increased, the critical Rayleigh number can increase by more than a factor of 10. The rotating magnetic field itself induces an instability at a critical value of the Hartmann number independent of the Rayleigh number. The nature of the bifurcations (whether subcritical or supercritical) and the convective flows occurring at the critical Rayleigh numbers are dependent upon the Hartmann number. In the stationary flow region, the rotating magnetic field can induce fluid motion with velocities on the order of cm/s. Thus, a rotating magnetic field can be used to achieve the benefits of forced convection during crystal growth processes without triggering deleterious instabilities.

Volz, Martin P.↗
DOE OSTI2025

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗
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↗
NASA NTRS1978

An instability of finite amplitude circularly polarized Alfven waves

A demonstration is presented that a finite-amplitude circularly-polarized Alfven wave is generally unstable in a MHD fluid. The wave decays by a four-wave coupling process in which the daughter waves are forward propagating random density and magnetic fluctuations and a backward-propagating magnetic wave. For parameters typical of the solar corona and the solar wind (thermal to magnetic energy density ratios between 0.1 and 1, and values between 0.1 and 0.9 for the ratio of magnetic energy density of the initial Alfven wave to that of the background magnetic field), large decay rates are found.

Goldstein, M. L.↗
NASA NTRS1989

Influence of instabilities on plasma flow around a comet

A multifluid hydrodynamic approximation allowing for the relative motion along the magnetic field of the newly created ions and the original fluid is used to treat the ion-pickup process. Due to the processes characterized by these means, the ion tail of a comet may not be antisolar; the derivation from radial is anticipated to be largest for oxygen due to its ionization at the greatest distances. Other ions, created nearer the comet where flow speed is lower, should have smaller transverse velocities.

Kellogg, Paul J.↗
NASA NTRS1983

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

For the case of a baroclinic flow whose Richardson number, Ri, is of order unity, a three-dimensional linear stability analysis is conducted on the basis of a model for a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The Hadley cell basic state and stability analysis are both based on the Navier-Stokes and energy equations, and perturbations possessing zonal, meridional, and vertical structures are considered. An attempt is made to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl and Ekman numbers, as well as to finite growth rates and a wider range of the zonal wavenumber. In general, it is found that the symmetric modes of maximum growth are not purely symmetric, but have a weak zonal structure.

Antar, B. N.↗
NASA NTRS1989

Calculations of inlet distortion induced compressor flow field instability

Calculations of the onset of flow instability are carried out for low-speed multistage axial compressors operating with asymmetric inlet flow. The modeling of the fluid dynamic interaction between the spoiled and unspoiled sectors of the compressor is the most important feature of the calculation procedure. The calculations show that annulus average slope of the compressor pressure rise characteristic equal to zero is a useful approximate stability criterion for situations where the dynamics of the compressor flow field do not couple strongly to the compression system or the structure of the imposed distortion is not similar to that of the eigenmodes of the flow in the compressor annulus. This criterion is employed to investigate the relationship between the present model and the 'parallel compressor' model. Calculations are also presented for cases when compressor and compressor system are closely coupled, and situations in which the compressor is subjected to a rotating distortion. These first-of-a-kind computations, and the accompanying description of the physical mechanisms, show that the stability of the flow in the compressor can be adversely affected if the temporal or spatial structure of the distortion is such that resonant type responses can be evoked either from the compressor or from compressor/compression system interactions.

Chue, R.↗
NASA NTRS2022

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

Hypersonic flow↗
DOE OSTI2021

Thermonuclear neutron emission from a sheared-flow stabilized Z-pinch

We report the fusion Z-pinch experiment (FuZE) is a sheared-flow stabilized Z-pinch designed to study the effects of flow stabilization on deuterium plasmas with densities and temperatures high enough to drive nuclear fusion reactions. Results from FuZE show high pinch currents and neutron emission durations thousands of times longer than instability growth times. While these results are consistent with thermonuclear neutron emission, energetically resolved neutron measurements are a stronger constraint on the origin of the fusion production. This stems from the strong anisotropy in energy created in beam-target fusion, compared to the relatively isotropic emission in thermonuclear fusion. In dense Z-pinch plasmas, a potential and undesirable cause of beam-target fusion reactions is the presence of fast-growing, “sausage” instabilities. This work introduces a new method for characterizing beam instabilities by recording individual neutron interactions in plastic scintillator detectors positioned at two different angles around the device chamber. Histograms of the pulse-integral spectra from the two locations are compared using detailed Monte Carlo simulations. These models infer the deuteron beam energy based on differences in the measured neutron spectra at the two angles, thereby discriminating beam-target from thermonuclear production. An analysis of neutron emission profiles from FuZE precludes the presence of deuteron beams with energies greater than 4.65 keV with a statistical uncertainty of 4.15 keV and a systematic uncertainty of 0.53 keV. This analysis demonstrates that axial, beam-target fusion reactions are not the dominant source of neutron emission from FuZE. These data are promising for scaling FuZE up to fusion reactor conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Asymmetric fluctuations and self-folding of active interfaces

We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid–liquid interfaces.

active fluid↗
NASA NTRS1985

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

A three-dimensional, linear stability analysis of a baroclinic flow for Richardson number, Ri, of order unity is presented. The model considered is a thin horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the complete set of governing, nonlinear equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in a closed form. The stability analysis is also based on the complete set of equations; and perturbation possessing zonal, meridional, and vertical structures were considered. Numerical methods were developed for the stability problem which results in a stiff, eighth-order, ordinary differential eigenvalue problem. The previous work on three-dimensional baroclinic instability for small Ri was extended to a more realistic model involving the Prandtl number, sigma, and the Ekman number, E, and to finite growth rates and a wider range of the zonal wavenumber.

Antar, B. N.↗
NASA NTRS1980

Baroclinic instability of flows along sloping boundaries

In the two-layer quasi-geostrophic model with boundaries sloping perpendicular to the basic flow, the ratios of the slopes of the bottom and the top to that of the interface between the fluid layers in the basic state are important parameters in the expression of the growth rate of unstable waves. When Eady's (1949) model is extended to include sloping bottom and top boundaries, the growth rates of unstable waves depend on the ratios of the slopes of the bottom and the top to that of the isentropes of the basic state. For the Eady model with sloping bottom, an important parameter characterizing the instability is the ratio between the vertical and horizontal heat transports by the wave divided by the slope of the isentropes of the basic state. An interpretation of these ratios and their relations clarifies the stabilization of the system for large slopes, the variation of the wavelength of the most unstable wave with the bottom slope, and the destabilization of some short waves for negative bottom slopes. It is found that the most unstable wave of the system has zero vertical energy flux convergence at the sloping bottom.

Mechoso, C. R.↗
DOE OSTI2020

Pair plasma instability in homogeneous magnetic guide fields

Pair plasmas, collections of both matter and antimatter particles of equal mass, represent a paradigm for the study of basic plasma science, and many open questions exist regarding these unique systems. They are found in many astrophysical settings, such as gamma-ray bursts, and have recently also been produced in carefully designed laboratory experiments. A central research topic in plasma physics is instability; however, unlike their more common ion–electron siblings, pair plasmas are generally thought to be stable to cross field pressure gradients in homogeneous magnetic fields. It is shown here by means of kinetic full-f simulations that, when a pressure gradient is first established, the Gradient-driven Drift Coupling mode is destabilized and becomes turbulent. Force balance is eventually achieved by a combination of flattened pressure profiles due to turbulent transport and establishment of a magnetic field gradient, saturating the growth. During the unstable phase, key physics can be captured by a δf gyrokinetic description, where it is shown analytically and numerically that parallel particle motion results in a coupling of all electromagnetic field components. A fluid model derived therefrom accurately predicts linear eigenmodes and is used to resolve global profile effects. For laser-based electron–positron plasma experiments, prompt instability is predicted with growth times much shorter than plasma lifetimes. Similarly, growth rates are calculated for the planned APEX experiment as well as gamma-ray burst scenarios, suggesting that the instability may contribute to the early evolution of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

On the Goertler vortex instability mechanism at hypersonic speeds

The linear instability of the hypersonic boundary layer on a curved wall is considered. As a starting point real-gas effects are ignored and the viscosity of the fluid is taken to be related to the temperature either by Chapman's Law or by Sutherland's Law. It is shown that the flow is susceptible to Goertler vortices. If Chapman's Law is used the vortices are trapped in the logarithmically thin adjustment layer in which the temperature of the basic flow changes rapidly to its free stream value and the nonuniqueness of the neutral stability curve associated with incompressible Goertler vortices is shown to disappear at high Mach numbers if the appropriate 'fast' streamwise dependence of the instability is built into the disturbance flow structure. If, on the other hand, Sutherland's Law is used, the vortices are found to spread into an O(1) region and the concept of a unique neutral stability curve is not tenable because of the nonparallel effects. For both laws the leading order terms in the expansions of the Goertler number are independent of the wave number and are due to the curvature of the basic state.

Hall, P.↗