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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Thresholds for the onset of fluid and magnetofluid turbulence

Linear stability calculations conducted in plasma physics are based on the theory of hydrodynamic stability of neutral fluids. The present investigation is concerned with the validity of procedures based on linear stability analysis in six much-studied situations. All have the property that the fluid goes from laminar to turbulent at critical values of some dimensionless number, such as the Reynolds number or the Rayleigh number. The question is, at what threshold values of the dimensionless number do the unstable motions set in, and can these thresholds be predicted by a linear analysis of the stability of the laminar state as the threshold is approached from the stable side. In three cases linear stability analysis clearly seems to fail. These cases include Plane Poiseuille Flow, Plane Couette Flow, and Cylindrical Pipe Flow (Hagen-Poiseuille Flow). Situations in which predictions provided by linear stability analysis are correct are related to Rotating Couette Flow, Thermally-Driven Convection, and Instability of Laminar Boundary Layers.

Montgomery, D.↗

NASA NTRS1985

A self-consistent linear-mode model of stellar convection

A normal-mode expansion of the linearized fluid equations in terms of small subset of spherical harmonics can provide a foundation for a physically motivated, self-consistent description of a solar-type convection zone. In the absence of dissipation, a second-order differential equation governs the radial dependence of the modes, so that interpretation of the effects on convection quantities of the normal-form 'potential well' is straightforward. The philosophy is quite different from the more recent work of Narasimha and Antia (1982): all envelopes presented here differ substantially from MLT envelopes, and therefore, from theirs, which are constructed to be consistent with MLT. The amplitude of all modes is set by a Kelvin-Helmholtz-('shear'-) instability argument unrelated to solar observations, with the result that the convection description may be considered to arise from 'first-hueristic-principles'. The thermodynamics modelled vaguely resemble the sun's, and more vigorously convective envelopes show some phenomena qualitatively like solar observations (e.g., atmospheric velocity spectra).

Macauslan, J.↗

NASA NTRS1989

A normal-mode approach to Jovian atmospheric dynamic

A nonlinear, quasi-geostrophic, baroclinic model of Jovian atmospheric dynamics is proposed, in which vertical variations of velocity are represented by a truncated sum over a complete set of orthogonal functions obtained by a separation of variables of the linearized quasi-geostrophic potential vorticity equation. A set of equations for the time variation of the mode amplitudes in the nonlinear case is then derived. It is shown that, for a planet with a neutrally stable, fluid interior instead of a solid lower boundary, the barotropic mode represents motions in the interior, and is not affected by the baroclinic modes. One consequence of this is that a normal-mode model with one baroclinic mode is dynamically equivalent to a one-layer model with solid lower topography. It is also shown that, for motions in Jupiter's cloudy lower troposphere, the stratosphere behaves nearly as a rigid lid, so that the normal-mode is applicable to Jupiter. The accuracy of the normal-mode model for Jupiter is tested using the following simple problems: (1) forced, vertically propagating Rossby waves, using two and three baroclinic modes, and (2) baroclinic instability, using two baroclinic modes. It is found that the normal-mode model provides qualitatively correct results, even with only a very limited number of vertical degrees of freedom.

Achterberg, Richard K.↗

NASA NTRS1996

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗

NASA NTRS1982

Three-dimensional stability of vortex arrays

The stability to three-dimensional disturbances of three classical steady vortex configurations in an incompressible inviscid fluid is studied in the limit of small vortex cross-sectional area and long axial disturbance wavelength. The configurations examined are the single infinite vortex row, the Karman vortex street of staggered vortices and the symmetric vortex street. It is shown that the single row is most unstable to a two-dimensional disturbance, while the Karman vortex street is most unstable to a three-dimensional disturbance over a significant range of street spacing ratios. The symmetric vortex street is found to be most unstable to three-dimensional or two-dimensional symmetric disturbances depending on the spacing ratio of the street. Short remarks are made concerning the relevance of the calculations to the observed instabilities in free shear layer, wake and boundary-layer type flows.

Robinson, A. C.↗

NASA NTRS1992

Perspectives on hypersonic viscous and nonequilibrium flow research

An attempt is made to reflect on current focuses in certain areas of hypersonic flow research by examining recent works and their issues. Aspects of viscous interaction, flow instability, and nonequilibrium aerothermodynamics pertaining to theoretical interest are focused upon. The field is a diverse one, and many exciting works may have either escaped the writer's notice or been abandoned for the sake of space. Students of hypersonic viscous flow must face the transition problems towards the two opposite ends of the Reynolds or Knudsen number range, which represents two regimes where unresolved fluid/gas dynamic problems abound. Central to the hypersonic flow studies is high-temperature physical gas dynamics; here, a number of issues on modelling the intermolecular potentials and inelastic collisions remain the obstacles to quantitative predictions. Research in combustion and scramjet propulsion will certainly be benefitted by advances in turbulent mixing and new computational fluid dynamics (CFD) strategies on multi-scaled complex reactions. Even for the sake of theoretical development, the lack of pertinent experimental data in the right energy and density ranges is believed to be among the major obstacles to progress in aerothermodynamic research for hypersonic flight. To enable laboratory simulation of nonequilibrium effects anticipated for transatmospheric flight, facilities capable of generating high enthalpy flow at density levels higher than in existing laboratories are needed (Hornung 1988). A new free-piston shock tunnel capable of realizing a test-section stagnation temperature of 10(exp 5) at Reynolds number 50 x 10(exp 6)/cm is being completed and preliminary tests has begun (H. Hornung et al. 1992). Another laboratory study worthy of note as well as theoretical support is the nonequilibrium flow experiment of iodine vapor which has low activation energies for vibrational excitation and dissociation, and can be studied in a laboratory with modest resources (Pham-Van-Diep et al. 1992).

Cheng, H. K.↗

NASA NTRS1984

An analysis of the stability of the compressible Ekman boundary layer

The linear stability problem for the compressible Ekman boundary layer common to rotating fluids is formulated and the stability properties determined numerically. Three classes of unstable waves are identified (called class A, B, and C), their properties are described. The class C waves have only recently been reported in the literature and are present only in compressible Ekman boundary layers. Most of the calculations presented here are for uranium hexafluoride gas; however, critical Reynolds numbers are also computed for air and ammonia gas. Compressibility is generally found to decrease the critical Reynolds number for each class of wave. A comparison of results for the three different gases shows the stability to be largely unaffected by changes in the gas properties. Maximum growth rate calculations for each wave show the class A and B waves to be the dominant instabilities.

Spall, J. R.↗

NASA NTRS1987

Straight cylindrical seal for high-performance turbomachines

A straight cylindrical seal configuration representing the seal for a high-performance turbopump (e.g., the space shuttle main engine fuel pump) was tested under static (nonrotating) conditions. The test data included critical mass flux and pressure profiles over a wide range of inlet temperatures and pressures for fluid nitrogen and fluid hydrogen with the seal in concentric and fully eccentric positions. The critical mass fluxes (or leakage rates) for the concentric and fully eccentric configurations were nearly the same when based on stagnation conditions upstream of the seal. The fully eccentric configuration pressure profiles of the gas and liquid were different. Further, the pressure differences between the maximum and the minimum clearance positions were highly dependent on the geometric conditions, the temperature, and the absolute pressure at both the inlet and the exit. The pressure differences were greatest in the inlet region. The results, although complex, tend to follow the corresponding-states principles for critical flows. Gaseous injection near the seal exit plane significantly altered the pressure profiles and could be used to control turbomachine instabilities.

Hendricks, Robert C.↗

NASA NTRS1991

Simulation of high-speed (orbital) releases of electron attachment materials in the ionosphere

The dynamics of ionospheric plasma irregularities produced by the release of electron attachment materials at orbital velocities across the geomagnetic field is studied. A two-dimensional electrostatic fluid model which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species transport is developed. Numerical simulations are performed to study the behavior at early and at late times after the release. At early times, of the order of or less than the attachment material neutral collision time, the negative ion cloud produced by the release may structure owing to the shear in the E x B velocity within the cloud. At high altitudes the cloud may bifurcate and form vortices on the back. At lower altitudes where ion-neutral collisional effects dominate, this structuring is suppressed. At late times, after a plasma depletion has formed due to neutralization chemistry, the cloud structures by the E x B interchange instability. Depending on the release altitude, the depletion structures by the collisional or inertial limit of this instability.

Scales, W. A.↗

NASA NTRS1993

Analytical description of the breakup of liquid jets in air

A viscous or inviscid cylindrical jet with surface tension in a vacuum tends to pinch due to the mechanism of capillary instability. Similarity solutions are constructed which describe this phenomenon as a critical time is encountered, for two physically distinct cases: inviscid jets governed by the Euler equations and highly viscous jets governed by the Stokes equations. In both cases the only assumption imposed is that at the time of pinching the jet shape has a radial length scale which is smaller than the axial length scale. For the inviscid case, we show that our solution corresponds exactly to one member of the one-parameter family of solutions obtained from slender jet theories and the shape of the jet is locally concave at breakup. For highly viscous jets our theory predicts local shapes which are monotonic increasing or decreasing indicating the formation of a mother drop connected to the jet by a thin fluid tube. This qualitative behavior is in complete agreement with both direct numerical simulations and experimental observations.

Papageorgiou, Demetrios T.↗

NASA NTRS2013

Characterization of Aerodynamic Interactions with the Mars Science Laboratory Reaction Control System Using Computation and Experiment

On August 5 , 2012, the Mars Science Laboratory (MSL) entry capsule successfully entered Mars' atmosphere and landed the Curiosity rover in Gale Crater. The capsule used a reaction control system (RCS) consisting of four pairs of hydrazine thrusters to fly a guided entry. The RCS provided bank control to fly along a flight path commanded by an onboard computer and also damped unwanted rates due to atmospheric disturbances and any dynamic instabilities of the capsule. A preliminary assessment of the MSL's flight data from entry showed that the capsule flew much as predicted. This paper will describe how the MSL aerodynamics team used engineering analyses, computational codes and wind tunnel testing in concert to develop the RCS system and certify it for flight. Over the course of MSL's development, the RCS configuration underwent a number of design iterations to accommodate mechanical constraints, aeroheating concerns and excessive aero/RCS interactions. A brief overview of the MSL RCS configuration design evolution is provided. Then, a brief description is presented of how the computational predictions of RCS jet interactions were validated. The primary work to certify that the RCS interactions were acceptable for flight was centered on validating computational predictions at hypersonic speeds. A comparison of computational fluid dynamics (CFD) predictions to wind tunnel force and moment data gathered in the NASA Langley 31-Inch Mach 10 Tunnel was the lynch pin to validating the CFD codes used to predict aero/RCS interactions. Using the CFD predictions and experimental data, an interaction model was developed for Monte Carlo analyses using 6-degree-of-freedom trajectory simulation. The interaction model used in the flight simulation is presented.

Schoenenberger, Mark↗

NASA NTRS2022

On the Stability and Phase Behavior of Titan's Subsurface Liquid Columns

On Titan, liquid hydrocarbon may stay in the subsurface porous permeable crust known as the alkanofer, analogous to water in Earth’s aquifer. In addition to pressure gradient, the subsurface liquid in alkanofers is subject to vertical compositional grading due to the gravity and temperature gradient. The common wisdom is that the liquid would normally stay underground in a stability established by the pressure that increases with depth as observed in aquifers on Earth. However, Titan’s liquids consist of nitrogen and hydrocarbons, mainly methane and ethane, the behavior of which is very sensitive to temperature and pressure. Consequently, the liquid density does not always increase with depth, thus may introduce a reverse density profile that leads to vertical convective instability of the liquid column. If reverse density profiles are present, capillary pressures arising from liquid trapped within small pores in the crust can help with the column stability. The liquid held in the capillaries can seal the space below it thus helping with the stability, unless the overpressure built from underneath becomes larger than the capillary pressure, which causes leakage to allow the liquids to seep upward from the deep. This situation is analogous to hydrocarbon seeps on Earth, where oil and natural gas escape the reservoir and flow slowly through network of cracks to the surface. An algorithm based on an extended Gibbs equation commonly used in petroleum reservoir engineering is employed in this work to produce pressure, density, and compositional profiles for evaluating the stability and phase behavior of Titan’s subsurface fluids.

Sugata Tan↗

NASA NTRS1993

Models of thermal/chemical boundary layer convection: Potential application to Venus

The upper boundary layer of Venus is comprised of at least two distinct chemical components, mantle and crust. Fluid dynamical models of convection within Venus' mantle were primarily of the thermal boundary layer type. Models assessing the ability of convective mantle flows to deform the crust were undertaken, but models exploring the effects of a variable thickness crust on mantle convection were largely lacking. A Venusian crust of variable thickness could couple back into, and alter, the mantle flow patterns that helped create it, leading to deformation mechanisms not predicted by purely thermal boundary layer convection models. This possibility is explored through a finite element model of thermal/chemical boundary layer convection. Model results suggest that a crust of variable thickness can serve as a mantle flow driver by perturbing lateral temperature gradients in the upper mantle. Resulting mantle flow is driven by the combination of free convective and nonuniform crustal distribution. This combination can lead to a flow instability manifest in the occurrence of episodic mantle lithosphere subduction initiated at the periphery of a crustal plateau. The ability of a light, near surface, chemical layer to potentially alter mantle flow patterns suggest that mantle convection and the creation and/or deformation of such a chemical layer may be highly nonseparable problems on time scales of 10(exp 8) years.

Lenardic, A.↗

NASA NTRS1994

Subduction and volatile recycling in Earth's mantle

The subduction of water and other volatiles into the mantle from oceanic sediments and altered oceanic crust is the major source of volatile recycling in the mantle. Until now, the geotherms that have been used to estimate the amount of volatiles that are recycled at subduction zones have been produced using the hypothesis that the slab is rigid and undergoes no internal deformation. On the other hand, most fluid dynamical mantle flow calculations assume that the slab has no greater strength than the surrounding mantle. Both of these views are inconsistent with laboratory work on the deformation of mantle minerals at high pressures. We consider the effects of the strength of the slab using two-dimensional calculations of a slab-like thermal downwelling with an endothermic phase change. Because the rheology and composition of subducting slabs are uncertain, we consider a range of Clapeyron slopes which bound current laboratory estimates of the spinel to perovskite plus magnesiowustite phase transition and simple temperature-dependent rheologies based on an Arrhenius law diffusion mechanism. In uniform viscosity convection models, subducted material piles up above the phase change until the pile becomes gravitationally unstable and sinks into the lower mantle (the avalanche). Strong slabs moderate the 'catastrophic' effects of the instabilities seen in many constant-viscosity convection calculations; however, even in the strongest slabs we consider, there is some retardation of the slab descent due to the presence of the phase change.

King, S. D.↗

NASA NTRS2004

Model of a Hollow Cathode Insert Plasma

A 2-D axisymmetric fluid model of the plasma in the insert region of a hollow cathode is presented. The level of sophistication included in the model is motivated in part by the need to determine quantitatively plasma fluxes to the emitter surface. The ultimate goal is to assess whether plasma effects can degrade the life of impregnated inserts beyond those documented throughout the 30-50 year history of vacuum cathode technologies. Results from simulations of a 1.2-cm diameter cathode operating at a discharge current of 25 A, and a gas flow rate of 5 sccm, suggest that approximately 10 A of electron current, and 3.5 A of ion current return to the emitter surface. The total emitted electron current computed by the model is about 35 A. Comparisons with plasma measurements suggest that anomalous heating of the plasma due to two-stream instabilities is possible near the orifice region. Solution to the heavy species energy equation, with classical transport and no viscous effects, predicts heavy species temperatures as high as 2640 K.

hollow cathodes↗

NASA NTRS1998

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗

NASA NTRS2005

Parametric Inlet Tested in Glenn's 10- by 10-Foot Supersonic Wind Tunnel

The Parametric Inlet is an innovative concept for the inlet of a gas-turbine propulsion system for supersonic aircraft. The concept approaches the performance of past inlet concepts, but with less mechanical complexity, lower weight, and greater aerodynamic stability and safety. Potential applications include supersonic cruise aircraft and missiles. The Parametric Inlet uses tailored surfaces to turn the incoming supersonic flow inward toward an axis of symmetry. The terminal shock spans the opening of the subsonic diffuser leading to the engine. The external cowl area is smaller, which reduces cowl drag. The use of only external supersonic compression avoids inlet unstart--an unsafe shock instability present in previous inlet designs that use internal supersonic compression. This eliminates the need for complex mechanical systems to control unstart, which reduces weight. The conceptual design was conceived by TechLand Research, Inc. (North Olmsted, OH), which received funding through NASA s Small-Business Innovation Research program. The Boeing Company (Seattle, WA) also participated in the conceptual design. The NASA Glenn Research Center became involved starting with the preliminary design of a model for testing in Glenn s 10- by 10-Foot Supersonic Wind Tunnel (10 10 SWT). The inlet was sized for a speed of Mach 2.35 while matching requirements of an existing cold pipe used in previous inlet tests. The parametric aspects of the model included interchangeable components for different cowl lip, throat slot, and sidewall leading-edge shapes and different vortex generator configurations. Glenn researchers used computational fluid dynamics (CFD) tools for three-dimensional, turbulent flow analysis to further refine the aerodynamic design.

Slater, John W.↗

NASA NTRSDate not supplied

Classical Flutter Analysis of X-57 Aircraft Mod II, III, and IV Configurations

To support airworthiness assessment of the X-57 electric propulsion demonstrator aircraft, the National Aeronautics and Space Administration (NASA) completed classical flutter analysis of the Mod II, Mod III, and Mod IV aircraft configurations. Aeroelastic analysis was primarily performed using the ZAERO analysis code developed by Zona Technology. Supplemental computational fluid dynamics (CFD) flutter analysis for the Mod III and Mod IV configurations was performed using NASA’s FUN3D software to verify the aeroelastic behavior at deployed flaps conditions. The NASTRAN finite element models (FEMs) used for aeroelastic analysis were developed by NASA. The final Mod II FEM utilized shell elements for the wing and fuselage, and beam elements for the tail and wing control surfaces. For Mod III and Mod IV aircraft models, shell element wing FEMs developed throughout the Mod III/IV wing design cycle were utilized. The fuselage and empennage FEM components were shared between the Mod II, III, and IV aircraft models. ZAERO analyses predicted flutter boundaries were >60% beyond the X-57 flight-test envelope never exceed airspeed (V NE ) for all configurations. Aeroelastic Reynolds-averaged Navier–Stokes (RANS) CFD simulations using FUN3D of the Mod III/IV configurations predicted flap mode instability did not occur until at least speeds of 60% higher than the Mod III/IV maximum flap deployment airspeed (V FE ).

Keerti K. Bhamidipati↗