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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Liquid neon heat transfer as applied to a 30 tesla cryomagnet

A 30-tesla magnet design is studied which calls for forced convection liquid neon heat transfer in small coolant channels. The design also requires suppressing boiling by subjecting the fluid to high pressures through use of magnet coils enclosed in a pressure vessel which is maintained at the critical pressure of liquid neon. This high pressure reduces the possibility of the system flow instabilities which may occur at low pressures. The forced convection heat transfer data presented were obtained by using a blowdown technique to force the fluid to flow vertically through a resistance heated, instrumented tube.

Papell, S. S.↗
NASA NTRS1992

Granular instability in fluidized beds: A small-over-large instability

For different particle diameters arranged in a gas fluidized bed, a new gravitational instability is predicted. The first order stability analysis yields conditions for growth of the most unstable wavenumber both as a function of particle diameter and dimensionless Bond number. Gravity acts to distribute the greater fluid mass associated with smaller particles, while surface tension stabilizes the entire bed against accelerations. The phenomenon, called the 'Small-over-Large' (SL) instability, evolves somewhat akin to a Rayleigh-Taylor mechanism and in addition to fundamental interest, carries implications for predicting solute redistribution and gravitational alignments in particle beds and chemical separation columns.

Noever, David A.↗
NASA NTRS2000

The Effect of a Rotating Magnetic Field on Flow Stability During Crystal Growth

The effect of a rotating magnetic field on the stability of flow in crystal growth configurations has been experimentally modeled using liquid gallium contained in a finite cylinder and heated from below. Several distinct flow regions were determined as a function of the Rayleigh and Hartmann numbers. At low values of the Rayleigh and Hartmann numbers, a region of stationary flow exists. As the rotating magnetic field is increased, the critical Rayleigh number bounding the stationary flow region can increase by a factor of 10. However, the rotating magnetic field itself induces an instability at a critical value of the Hartmann number independent of the Rayleigh number. In the stationary flow region, the rotating magnetic field can induce fluid motion with velocities several orders of magnitude larger than typical semiconductor crystal growth velocities. Thus, a rotating magnetic field can be used to achieve the benefits of forced convection without triggering deleterious instabilities.

Volz, Martin P.↗
NASA NTRS1985

Formation of Bright Terrain Bands on Ganymede by Unstable Lithosphere Extension

The surface of Ganymede, as seen in Voyager images, was greatly modified by tectonic activity. Geologic evidence is consistent with the formation of younger bright terrain by the emplacement of icy material into rift zones formed in the older dark terrain. To better understand the formation of both terrestrial rift zones and bright terrain bands on Ganymede, the mechanics of rift initiation by the growth of necking instabilities in an extending lithosphere are studied. In particular, how an initial, small perturbation in lithosphere thickness at the lithosphere/asthenosphere boundary is amplified by uniform horizontal extension is studied. The lithosphere is treated as a layer of power-law viscous fluid of thickness H overlying a weaker, viscous mantle substrate. In order to assess the relative importance of buoyant upwelling and mechanical instability in rift formation, both density and strength stratification were incorporated into the model.

Zuber, M. T.↗
NASA NTRS1992

Magnetic interchange instability of accretion disks

The nonlinear evolution of the magnetic interchange or buoyancy instability of a differentially rotating disk threaded by an ordered vertical magnetic field is investigated. A 2D ideal fluid in the equatorial plane of a central mass in the corotating frame of reference is considered as a model for the disk. If the rotation rate of the disk is Keplerian, the disk is found to be stable. If the vertical magnetic field is sufficiently strong, and the field strength decreases with distance from the central object, and thus the rotation of the disk deviates from Keplerian, if is found that an instability develops. The magnetic flux and disk matter expand outward in certain ranges of azimuth, while disk matter with less magnetic flux moves inward over the remaining range of azimuth, showing a characteristic development of an interchange instability.

Kaisig, M.↗
NASA NTRS1993

Flow instability in particle-bed nuclear reactors

A three-dimensional model of the stability of the particle-bed reactor is presented, in which the fluid has mobility in three dimensions. The model accurately represents the stability at low Re numbers as well as the effects of the cold and hot frits and of the heat conduction and radiation in the particle bed. The model can be easily extended to apply to the cylindrical geometry of particle-bed reactors. Exemplary calculations are carried out, showing that a particle bed without a cold frit would be subject to instability if operated at the high-temperature ratios used for nuclear rockets and at power densities below about 4 MW/l; since the desired power density for such a reactor is about 40 MW/l, the operation at design exit temperature but at reduced power could be hazardous. Calculations show however that it might be possible to remove the instability problem by appropriate combinations of cold and hot frits.

Kerrebrock, J. L.↗
DOE OSTI2024

Instabilities and Mixing in Inertial Confinement Fusion

By imploding fuel of hydrogen isotopes, inertial confinement fusion (ICF) aims to create conditions that mimic those in the Sun's core. This is fluid dynamics in an extreme regime, with the ultimate goal of making nuclear fusion a viable clean energy source. The fuel must be reliably and symmetrically compressed to temperatures exceeding 100 million degrees Celsius. After the best part of a century of research, the foremost fusion milestone was reached in 2021, when ICF became the first technology to achieve an igniting fusion fuel (thermonuclear instability), and then in 2022 scientific energy breakeven was attained. A key trade-off of the ICF platform is that greater fuel compression leads to higher burn efficiency, but at the expense of amplified Rayleigh–Taylor and Richtmyer–Meshkov instabilities and kinetic-energy-wasting asymmetries. In extreme cases, these three-dimensional instabilities can completely break up the implosion. Even in the highest-yielding 2022 scientific breakeven experiment, high-atomic-number (high-Z) contaminants were unintentionally injected into the fuel. Here we review the pivotal role that fluid dynamics plays in the construction of a stable implosion and the decades of improved understanding and isolated experiments that have contributed to fusion ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Instability during directional solidification - Gravitational effects

After an elementary introduction to solidification, recent research that attempts to understand the interaction of fluid flow with a crystal-melt interface under gravitational influence using linear stability analysis is reviewed. Recent numerical calculations of fluid flow and interface morphology during alloy solidification are briefly discussed.

Coriell, S. R.↗
NASA NTRS1989

On the Gortler vortex instability mechanism at hypersonic speeds

The linear instability of the hypersonic boundary layer on a curved wall is considered. As a starting point the viscosity of the fluid is taken to be a linear function of temperature and real-gas effects are ignored. It is shown that the flow is susceptible to Gortler vortices and that they are trapped in the logarithmically thin adjustment layer in which the temperature of the basic flow changes rapidly to its free stream value. The vortices decay exponentially in both directions away from this layer and are most unstable when their wavelength is comparable with the depth of the adjustment layer. The non-uniqueness of the neutral stability curve associated with incompressible Gortler 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. It is shown that in the hypersonic limit wall-cooling has a negligible effect on the stability of a fluid with a given value of the Chapman constant.

Hall, Philip↗
NASA NTRSDate not supplied

Surface Instability of Liquid Propellants in Microgravity During Pulsed Settling Operations

Pulsing reaction control system (RCS) thrusters, vent valves, or other propulsion devices can preserve propellant resources in space, but this operation also effectively introduces a vibration to the vehicle. When the vibration is perpendicular to the liquid propellant surface, Faraday waves may be generated at the liquid-vapor interface. These Faraday instabilities can perturb or break up the liquid surface of cryogenic tanks, leading to inefficiencies in thermal management or even ullage collapse. Drawing from theory and experiments, an engineering model defining the allowable design regions for pulsed settling in microgravity was assembled and verified with computational fluid dynamics (CFD) simulations. A traditional settling metric, the Bond number, was also overlaid in the model to indicate which duty cycles were insufficient to overcome surface tension and aggregate propellant. Mission planners and engineers can consult the tool to rapidly evaluate the stability of a liquid-vapor interface given the pulse frequency and the excitation acceleration. Expressions developed for Faraday waves induced by a sinusoidal forcing input at standard gravity were found to provide excellent predictive capabilities for pulsed, or rectangular, waveforms in the absence of a consistent gravitational acceleration. This study extends the usage of these equations to an alternative forcing function and microgravity environments for the purpose of estimating natural frequencies, surface mode shapes, surface wave amplitudes, and the onset of droplet ejection. CFD simulations with the Loci/STREAM-VoF (Volume of Fluid) solver were initially validated against experimental results in standard gravity. Discrete points on the design map were then investigated with CFD and confirmed that the engineering model reliably indicates surface stability and most Faraday wave characteristics without requiring higher-fidelity tools. The engineering model is highly extensible and can be adapted for various propellant fill fractions, fluid properties, and tank sizes.

Faraday Waves↗
NASA NTRS1987

The nonlinear interaction of Tollmien-Schlichting waves and Taylor-Goertler vortices in curved channel flows

It is known that a viscous fluid flow with curved streamlines can support both Tollmien-Schlichting and Taylor-Goertler instabilities. In a situation where both modes are possible on the basis of linear theory a nonlinear theory must be used to determine the effect of the interaction of the instabilities. The details of this interaction are of practical importance because of its possible catastrophic effects on mechanisms used for laminar flow control. This interaction is studied in the context of fully developed flows in curved channels. A part form technical differences associated with boundary layer growth the structures of the instabilities in this flow are very similar to those in the practically more important external boundary layer situation. The interaction is shown to have two distinct phases depending on the size of the disturbances. At very low amplitudes two oblique Tollmein-Schlichting waves interact with a Goertler vortex in such a manner that the amplitudes become infinite at a finite time. This type of interaction is described by ordinary differential amplitude equations with quadratic nonlinearities.

Hall, P.↗
NASA NTRS1987

Plasma transport in the auroral return current region

The classical and anomalous transport properties of a multifluid plasma consisting of H(+), O(+), and electron populations in the presence of auroral field-aligned return currents are investigated, using a multimoment fluid model with anomalous transport coefficients. The macroscopic effects of the electrostatic ion cyclotron (EIC) instability and of an EIC-related anomalous resistivity mechanism which heats the electrons are included in the present version of the model. The responses of the outflowing polar wind plasma to the application of current, with and without instabilities, are exhibited. The simulations show that the electron drift velocity corresponding to a return current of 0.65 micro-A/sq m is above the threshold for EIC waves. Downward electron heat flow competes with upward convection and adiabatic effects to determine the direction of the electron temperature anisotropy. Resistive electron heating lowers the critical drift velocity for marginal EIC stability and leads to enhanced ion heating.

Ganguli, Supriya B.↗
NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗
NASA NTRS1992

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗
NASA NTRS1992

Numerical simulation of jet noise

Jet noise and jet-induced structural loads have become key issues in the design of commercial and military aircraft. Computational Fluid Dynamics (CFD) can be of use in predicting the underlying jet shear-layer instabilities and, in conjunction with classical acoustic theory, jet noise. The computational issues involved in the resolution of high Reynolds number unsteady jet flows are addressed in this paper. Once these jet flows can be accurately resolved, it should be possible to use acoustic theory to extract, for example, the far-field jet noise. An assessment of future work and computational resources required for directly computing far-field jet noise is also presented.

Van Dalsem, W. R.↗
NASA NTRS1998

Application of Aeroelastic Solvers Based on Navier-Stokes Equations

A pre-release version of the Navier-Stokes solver (TURBO) was obtained from MSU. Along with Dr. Milind Bakhle of the University of Toledo, subroutines for aeroelastic analysis were developed and added to the TURBO code to develop versions 1 and 2 of the TURBO-AE code. For specified mode shape, frequency and inter-blade phase angle the code calculates the work done by the fluid on the rotor for a prescribed sinusoidal motion. Positive work on the rotor indicates instability of the rotor. The version 1 of the code calculates the work for in-phase blade motions only. In version 2 of the code, the capability for analyzing all possible inter-blade phase angles, was added. The version 2 of TURBO-AE code was validated and delivered to NASA and the industry partners of the AST project. The capabilities and the features of the code are summarized in Refs. [1] & [2]. To release the version 2 of TURBO-AE, a workshop was organized at NASA Lewis, by Dr. Srivastava and Dr. M. A. Bakhle, both of the University of Toledo, in October of 1996 for the industry partners of NASA Lewis. The workshop provided the potential users of TURBO-AE, all the relevant information required in preparing the input data, executing the code, interpreting the results and bench marking the code on their computer systems. After the code was delivered to the industry partners, user support was also provided. A new version of the Navier-Stokes solver (TURBO) was later released by MSU. This version had significant changes and upgrades over the previous version. This new version was merged with the TURBO-AE code. Also, new boundary conditions for 3-D unsteady non-reflecting boundaries, were developed by researchers from UTRC, Ref. [3]. Time was spent on understanding, familiarizing, executing and implementing the new boundary conditions into the TURBO-AE code. Work was started on the phase lagged (time-shifted) boundary condition version (version 4) of the code. This will allow the users to calculate non-zero interblade phase angles using, only one blade passage for analysis.

Keith, Theo G., Jr.↗
NASA NTRS1999

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗
NASA NTRS2009

Modeling Techniques for Evaluation the Effectiveness of Particle Damping in Turbomachinery

High power turbopumps are frequently used to supply propellants to the combustion chambers of rocket engines. Due to the high pressures and flow-rates required, turbopump components are subjected to harsh environments which include dynamic excitation due to random, sine, and acoustic vibration. Additionally, fluid-induced forces can couple with the dynamics of the structure resulting in flow induced instabilities (flutter). Structural response to these forms of excitation results in reduced fatigue life and increases the likelihood of an operational failure. Particle damping has been used successfully on vibration problems in the past by increasing the damping and therefore reducing the response to acceptable levels. Empirical methods have typically been employed to evaluate the performance of the particles in reducing the structural response. This report explores the use of finite element methods to estimate the effectiveness of particle damping in a typical non-rotating turbopump component. Axisymmetric harmonic models are used to estimate the increase in modal damping produced by the addition of particles in the cavity of an axisymmetric seal. Target modes of vibration are evaluated to quantify how the effective particle damping is altered by geometry changes in the seal design. A new method to predict the performance of particle dampers is developed and shown to provide more reasonable estimates of damping.

Ehrgott, R.↗