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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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.↗
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.↗
NASA NTRS2011

Effects of Inlet Distortion on Aeromechanical Stability of a Forward-Swept High-Speed Fan

Concerns regarding noise, propulsive efficiency, and fuel burn are inspiring aircraft designs wherein the propulsive turbomachines are partially (or fully) embedded within the airframe; such designs present serious concerns with regard to aerodynamic and aeromechanic performance of the compression system in response to inlet distortion. Separately, a forward-swept high-speed fan was developed to address noise concerns of modern podded turbofans; however this fan encounters aeroelastic instability (flutter) as it approaches stall. A three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is applied to analyze and corroborate fan performance with clean inlet flow. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions, is modified and then applied in a computational study to preliminarily assess the effects of inlet distortion on aeroelastic stability of the fan. Computational engineering application and implementation issues are discussed, followed by an investigation into the aeroelastic behavior of the fan with clean and distorted inlets.

Herrick, Gregory P.↗
NASA NTRS1999

Extensional Flow of a Polystyrene Boger Fluid Through a 4:1:4 Axisymmetric Contraction/Expansion

The creeping flow of a dilute (0.025 wt%) monodisperse polystyrene/polystyrene Boger fluid through a 4:1:4 axisymmetric contraction/expansion is experimentally observed for a wide range of Deborah numbers. Pressure drop measurements across the orifice plate show a large extra pressure drop that increases monotonically with Deborah number above the value observed for a similar Newtonian fluid at the same flow rate. This enhancement in the dimensionless pressure drop is not associated with the onset of a flow instability, yet it is not predicted by existing steady-state or transient numerical computations with simple dumbbell models. It is conjectured that this extra pressure drop is the result of an additional dissipative contribution to the polymeric stress arising from a stress-conformation hysteresis in the strong non-homogeneous extensional flow near the contraction plane. Such a hysteresis has been independently measured and computed in recent studies of homogeneous transient uniaxial stretching of PS/PS Boger fluids. Flow visualization and velocity field measurements using digital particle image velocimetry (DPIV) show large upstream growth of the corner vortex with increasing Deborah number. At large Deborah numbers, the onset of an elastic instability is observed, first locally as small amplitude fluctuations in the pressure measurements, and then globally as an azimuthal precessing of the upstream corner vortex accompanied by periodic oscillations in the pressure drop across the orifice.

Rothstein, Jonathan P.↗
NASA NTRS1994

The proton cyclotron instability and the anisotropy/B inverse correlation

Spacecraft observations in the strongly compressed subsolar magnetosheath show an inverse correlation between the proton temperature anistrophy (T(sub perpendicular to p)/T(sub parallel to p) greater than 1 where perpendicular to and parallel to denote directions perpendicular and parallel to the background magnetic field) and the parallel proton beta (beta(sub parallel to p). This manuscript uses one-dimensional hybrid simulations of the proton cyclotron anisotropy instability in homogeneous electron-proton plasmas to study this correlation which may represent a limited closure relation for fluid theories of anisotropic space plasmas. The emphasis is on driven simulation which increase the temperature anisotropy by periodically reducing the magnetic-field-aligned velocities of the protons. The late-time states from ensembles of both initial value and driven simulations yield very similar expressions for the proton anisotropy/beta(sub parallel to P) inverse correlation, and provide a basis for explaining differences between sheath observations from different spacecraft. The driven simulations also yield expressions for the maximum instability growth rate and the fluctuating field energy as functions of beta (sub parallel to p) and a parameter characterizing the anisotropy driver.

Gary, S. Peter↗
NASA NTRS1989

Numerical investigation of the control of the secondary instability process in boundary layers

A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.

Kral, L. D.↗