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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Advanced stability analysis for laminar flow control

Five classes of problems are addressed: (1) the extension of the SALLY stability analysis code to the full eighth order compressible stability equations for three dimensional boundary layer; (2) a comparison of methods for prediction of transition using SALLY for incompressible flows; (3) a study of instability and transition in rotating disk flows in which the effects of Coriolis forces and streamline curvature are included; (4) a new linear three dimensional instability mechanism that predicts Reynolds numbers for transition to turbulence in planar shear flows in good agreement with experiment; and (5) a study of the stability of finite amplitude disturbances in axisymmetric pipe flow showing the stability of this flow to all nonlinear axisymmetric disturbances.

Orszag, S. A.↗
DOE OSTI2015

Renormalized theory of ion temperature gradient instability of the magnetic-field-aligned plasma shear flow with hot ions

The developed kinetic theory for the stability of a magnetic-field-aligned (parallel) shear flow with inhomogeneous ion temperature [Mikhailenko et al., Phys. Plasmas 21, 072117 (2014)] predicted that a kinetic instability arises from the coupled reinforcing action of the flow velocity shear and ion temperature gradient in the cases where comparable ion and electron temperatures exist. In the present paper, the nonlinear theory was developed for the instability caused by the combined effects of ion-temperature-gradient and shear-flow (ITG–SF). The level of the electrostatic turbulence is determined for the saturation state of the instability on the basis of the nonlinear dispersion equation, which accounts for a nonlinear scattering of ions by the developed turbulence in a sheared flow. The renormalized quasilinear equation for the ion distribution function, which accounts for the turbulent scattering of ions by ITG–SF driven turbulence, was derived and employed for the estimation of the turbulent ion viscosity, the anomalous ion thermal conductivity, and anomalous ion heating rate at the saturation state of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Final Technical Report for DOE DE-FG02-05ER54831 "Laboratory Studies of Dynamos."

Laboratory Studies of Dynamos: Executive Summary. The self-generation of magnetic fields by astrophysical bodies like planets, stars, accretion disks, galaxies, and even galaxy clusters arises due to a mechanism referred to as a homogeneous dynamo. It is quite simple to demonstrate the generation of a magnetic fi eld from a rotating copper disk coupled with a coil of wire, a device known as the homopolar dynamo. The device works like a magnetic fi eld ampli er with a feedback circuit: the differential rotation of a metal disk past an infinitesimally small seed magnetic field induces currents in the disk which, when coupled to a coil winding, can amplify the field until it becomes strong enough to slow the rotation of the disk. What is remarkable is that the same type of circuit may be achieved in a flowing conducting fluid such as a liquid metal in the case of planetary dynamos or a plasma in the case of astrophysical dynamos. The complexity of describing planetary and stellar dynamos despite their ubiquity and the plethora of observational data from the Earth and the Sun motivates the demonstration of a laboratory homogenous dynamo. To create a homogenous dynamo, one first needs a su fficiently large, fast flow of a highly conducting fluid that the velocity shear in the fluid can bend magnetic field lines. With a high Rm-flow, the magnetic fi eld can be ampli ed by the stretching action provided by di fferential rotation. The other critical ingredient is a flow geometry that provides feedback so that the ampli ed eld reinforces the initial in nitesimal seed field - a mechanism that recreates the feedback provided by the coil of wire in the homopolar dynamo. In the Madison Dynamo Experiment, this combination of magnetic ampli cation and feedback is feasible in the simple geometry of two counter-rotating helical vortices in a 1 meter-diameter spherical vessel lled with liquid sodium. For an optimal helical pitch of the flow the threshold for exciting a dynamo is predicted from laminar flow modeling to be at peak flow speeds of 5 m/s. Liquid metals tend to have viscosities similar to that of water yielding inviscid flows. Whereas the timescale for the dynamo instability is on the resistive dissipation time, the timescale for hydrodynamic instability of the shear layer is quite short meaning that the shear layer required to generate the magnetic eld is broken up by Kelvin-Helmholtz instabilities. The eddies generated by large-scale flow drive instabilities at progressively smaller scale giving rise to a cascade of turbulent eddies driven at the largest scale of the experiment. The major contribution of the Madison Dynamo Experiment has been quantifying the role this turbulence plays in the generation of magnetic elds. Overall, the Madison Dynamo Experiment has now operated for about 1 decade and carried out experiments related to magnetic fi eld generation by turbulent flows of liquid metal. The principle thrust of research and indeed the main scienti fic outcomes are related to how turbulent flows create and transport magnetic fi elds.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state: A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state

The development of numerical approaches to perform direct numerical simulations of compressible multiphase flows has been an active field of research for several years. Proper treatment of fluid interfaces is crucial as important physics occur in this infinitesimally small region. Furthermore, the compressibility of the fluid requires proper treatment of discontinuities. Artificial diffusivity is among a number of methods widely used for compressible flows. This study develops a general form of consistent artificial diffusion fluxes and extends the localized artificial diffusivity method for high-order central schemes to solve multiphase flows with an interface-capturing method. These fluxes ensure an oscillation-free interface for pressure, velocity, and temperature without employing a sharpening technique. Moreover, the high-order representation of all scales in the flow helps capture the wide range of instabilities inherent in these flows. The goal is to develop an approach capable of performing high-fidelity simulations supported by physics-driven validation. lastly, this is achieved by solving the five-equation model with the stiffened-gas equation of state using the proposed method for multicomponent and multiphase flows on a variety of 1D and 2D problems.

42 ENGINEERING↗
NASA NTRS1983

Instabilities of dynamic thermocapillary liquid layers. I Convective instabilities

A steady shear flow is set up in a planar liquid layer bounded below by a rigid plate, and above by an interface with a passive gas, by imposing a temperature gradient along the layer and driving the motion by thermocapillarity. The thermal-convective instabilities to which this dynamic state is susceptible include stationary longitudinal rolls involving the classical Marangoni instability, and unsteady hydrothermal waves, whose new instability mechanism derives its energy from the horizontal temperature gradients. The thermal stability characteristics for liquid layers with and without return flow profiles are presently given as functions of the Prandtl number for the liquid and the Biot number for the interface, and comparisons are made with available experimental observations.

Smith, M. K.↗
NASA NTRS1990

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗
NASA NTRS1997

Entry Dispersion Analysis for the Stardust Comet Sample Return Capsule

Stardust will be the first mission to return samples from beyond the Earth-Moon system. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study analyzes the entry, descent, and landing of the returning sample capsule. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by an accelerometer and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) shows that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust program limits. Additionally, the size of the resulting 3-sigma landing ellipse is 83.5 km in downrange by 29.2 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗
NASA NTRS1999

Entry Trajectory Issues for the Stardust Sample Return Capsule

The Stardust mission was successfully launched on February 7, 1999. It will be the first mission to return samples from a comet. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study describes the analysis of the entry, descent, and landing of the returning sample capsule utilizing the final, launch configuration capsule mass properties. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by a gravity-switch and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) predicts that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust mission limits. Additionally, the size of the resulting 3-sigma landing ellipse is 60.8 km in downrange by 19.9 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗
NASA NTRS1999

Entry Trajectory Issues for the Stardust Sample Return Capsule

The Stardust mission was successfully launched on February 7, 1999. It will be the first mission to return samples from a comet. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study describes the analysis of the entry, descent, and landing of the returning sample capsule utilizing the final, launch configuration capsule mass properties. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by a gravity-switch and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) predicts that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust mission limits. Additionally, the size of the resulting 3-sigma landing ellipse is 60.8 km in downrange by 19.9 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗
DOE OSTI2013

On secondary and tertiary instability in electromagnetic plasma microturbulence

Zonal flows, widely accepted to be the secondary instability process leading to the nonlinear saturation of ion temperature gradient modes, are shown to grow at higher rates relative to the linear mode amplitude as the plasma pressure β is increased—thus, confirming that zonal flows become increasingly important in the turbulent dynamics at higher β. At the next level of nonlinear excitation, radial corrugations of the distribution function (zonal flow, zonal density, and zonal temperature) are demonstrated to modify linear growth rates moderately through perturbed-field, self-consistent gradients: on smaller scales, growth rates are reduced below the linear rate. In particular, excitation of kinetic ballooning modes well below their usual threshold is not to be expected under normal conditions. These findings strengthen the theory of the non-zonal transition [M. J. Pueschel et al., Phys. Rev. Lett. 110, 155005 (2013)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗
NASA NTRS1984

Morphological stability in the presence of fluid flow in the melt

Recent experiments have shown that the presence of a vertical buoyancy-driven flow adjacent to an initially cylindrical crystal-melt interface may produce a time-dependent helical deformation of the interface, with a rotation period ranging from several minutes to many hours, depending upon the width of the melt. The temperature distribution is such that the interface is expected to be morphologically stable in the absence of fluid flow. A linear stability analysis reveals that the instability is due to a coupling between a basic hydrodynamic instability in the buoyant flow and the deformable boundary separating the two phases. The crystal-melt interface lowers the critical Grashof number of an analogous rigid-walled system by an order of magnitude for succinonitrile with a Prandtl number P = 22.8; furthermore, the hydrodynamic mode that is actually destabilized by the interface is not the least stable mode in the rigid-walled system for P = 22.8. The results show that the instability may be regarded either as a rather large alteration of a basic hydrodynamic instability by the crystal-melt interface, or as a significant modification of the morphological stability of the interface by the presence of the buoyant flow.

Mcfadden, G. B.↗
DOE OSTI2016

STOCHASTIC PARTICLE ACCELERATION IN TURBULENCE GENERATED BY MAGNETOROTATIONAL INSTABILITY

We investigate stochastic particle acceleration in accretion flows. It is believed that magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation and without back reaction to the field. Our results show that the CRs randomly gain or lose their energy through interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: the diffusion coefficient in the direction of the unperturbed flow is about 20 times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration works efficiently for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction undergo unphysical runaway acceleration by the Lorentz transformation, which needs to be taken into account with special care.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Turbulence modeling in aerodynamic shear flows - Status and problems

This paper briefly summarizes the status and problems of turbulence modeling for aerodynamical applications. For complex flows the 'approach of choice' is (increasingly) full second-order (Reynolds stress equation) closure. These closures have not yet developed to anywhere near their full potential, significant further research is required especially regarding length-scale equations, representation of pressure-strain correlations, and wall region treatments. Recent developments in computer capability, algorithms, numerical simulations, theory and quantitative flow visualization should assist in and hasten this research. Several problem areas such as shock interaction and discrete dynamic instabilities of turbulent flows may require mega-to-large eddy simulation or theoretical adjuncts.

Bushnell, D. M.↗
DOE OSTI2012

PLANETESIMAL FORMATION AT THE BOUNDARY BETWEEN STEADY SUPER/SUB-KEPLERIAN FLOW CREATED BY INHOMOGENEOUS GROWTH OF MAGNETOROTATIONAL INSTABILITY

We have studied formation of planetesimals at a radial pressure bump in a protoplanetary disk created by radially inhomogeneous magnetorotational instability (MRI), through three-dimensional resistive MHD simulations including dust particles. In our previous papers, we showed that the inhomogeneous MRI developing in non-uniform structure of magnetic field or magnetic resistivity can transform the local gas flow in the disk to a quasi-steady state with local rigid rotation that is no longer unstable against the MRI. Since the outer part of the rigid rotation is super-Keplerian flow, a quasi-static pressure bump is created and dust concentration is expected there. In this paper, we perform simulations of the same systems, adding dust particles that suffer gas drag and modulate gas flow via the back-reaction of the gas drag (dust drag). We use {approx}O(10{sup 7}) super-particles, each of which represents {approx}O(10{sup 6})-O(10{sup 7}) dust particles with sizes of centimeter to meter. We have found that the dust drag suppresses turbulent motion to decrease the velocity dispersion of the dust particles while it broadens the dust concentrated regions to limit peaky dust concentration, compared with the simulation without the dust drag. We found that the positive effect for the gravitational instability (GI), reduction in the velocity dispersion, dominates over the negative one, suppression in particle concentration. For meter-size particles with the friction time {tau}{sub f} {approx_equal} 1/{Omega}, where {Omega} is Keplerian frequency, the GI of the dust particles that may lead to planetesimal formation is expected. For such a situation, we further introduced the self-gravity of dust particles to the simulation to demonstrate that several gravitationally bound clumps are actually formed. Through analytical arguments, we found that planetesimal formation from meter-sized dust particles is possible at {approx}5 AU, if dust spatial density is a few times larger than that in the minimum mass solar nebula.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

On the thermal instability of galactic and cluster halos

The paper presents a detailed study of thermal instabilities in cooling flows associated with galaxies and clusters of galaxies. In the case of purely radiation-driven accretion onto a central object such as the cD galaxy M87, it is found that the gas is largely subject to overstability, rather than to monotonic instability. If thermal conductivity is taken into account, the flow is stabilized on scales of several kiloparsecs, even if the conductivity is appreciably reduced (e.g., about 1 percent) with respect to the Spitzer value. In no case are the globular perturbations (i.e., perturbations with comparable radial and azimuthal dimensions) found to be monotonically unstable. The paper presents numerical solutions of the local dispersion relation for the cooling flow in M87 and discusses the possible consequences of the results for a correct understanding of cooling flows.

Malagoli, A.↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗