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 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Primary zone dynamics in a gas turbine combustor

Fluid mechanical investigations simulating the flow in the primary zone of a gas turbine combustor are presented using three generic test rigs: (1) rotating pipe yielding a swirling jet of air; (2) primary zone model with a single swirler and various primary jet configurations, operated with air; and (3) two rectangular models of a (stretched-out) annular combustor with five swirlers in the backwall and with various primary jet configurations, one operated with air and the other with water. Concentration measurements are obtained using laser sheet imaging techniques and velocity measurements using a laser Doppler velocimeter. The results show recirculation zones, intense mixing, instabilities of the interacting jets and the presence of large random vortical motions. The flowfields are shown to exhibit bimodal behavior, have asymmetries despite symmetrical geometry and inlet conditions and display strong jet/swirler and swirler/swirler interactions.

Sullivan, J. P.↗

NASA NTRS2003

Suppression of Buoyancy in Gaseous Media at High Temperatures

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side-walls (Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no matter how large the temperature (density) difference becomes at a given pressure even under normal gravity! This unexpected behavior is primarily attributed to highly temperature-sensitive kinematic viscosity which counteracts the tendency toward instability and dampens convection by making the gas more viscous at higher temperatures. This compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature beyond which a system will tend to be more stable as the hot-surface temperature increases.

Gokoglu, Suleyman A.↗

NASA NTRS2016

Heating and Acceleration of Solar Wind Ions by Turbulent Wave Spectrum in Inhomogeneous Expanding Plasma

Near the Sun (< 10Rs) the acceleration, heating, and propagation of the solar wind are likely affected by the background inhomogeneities of the magnetized plasma. The heating and the acceleration of the solar wind ions by turbulent wave spectrum in inhomogeneous plasma is studied using a 2.5D hybrid model. The hybrid model describes the kinetics of the ions, while the electrons are modeled as massless neutralizing fluid in an expanding box approach. Turbulent magnetic fluctuations dominated by power-law frequency spectra, which are evident from in-situ as well as remote sensing measurements, are used in our models. The effects of background density inhomogeneity across the magnetic field on the resonant ion heating are studied. The effect of super- Alfvenic ion drift on the ion heating is investigated. It is found that the turbulent wave spectrum of initially parallel propagating waves cascades to oblique modes, and leads to enhanced resonant ion heating due to the inhomogeneity. The acceleration of the solar wind ions is achieved by the parametric instability of large amplitude waves in the spectrum, and is also affected by the inhomogeneity. The results of the study provide the ion temperature anisotropy and drift velocity temporal evolution due to relaxation of the instability. The non-Maxwellian velocity distribution functions (VDFs) of the ions are modeled in the inhomogeneous solar wind plasma in the acceleration region close to the Sun.

Ofman, Leon↗

NASA NTRS1998

Microgravity Experiments and Numerical Modeling of Rotating Buoyant Convection in a Spherical Shell with Latitudinal Thermal Gradients

The results of numerical model calculations are compared with space laboratory experiments for rotating, thermally driven flow in a hemispherical shell, and various flow regimes are described. A spherically symmetric body force, analogous to gravity, is imposed in the radial direction through the use of a dielectric fluid and an electrostatic potential difference across the gap. The spherical boundaries are maintained at constant temperature profiles (with the pole and the inner sphere being the warmer surfaces), and the equatorial wall is an insulator. Typical parameter combinations result in highly nonlinear, but laminar, flow. For weak enough buoyant forcing, the flow is axisymmetric. It consists of a single meridional cell, rising in warm latitudes (in this case, near the pole) and sinking in cool latitudes, with prograde flow in the equatorial region and near the inner hemisphere for other latitudes, and retrograde motion elsewhere. For fast rotation, the first transition due to the instability of this simple flow is also axisymmetric, consisting of rings of convection in the polar region, which propagate poleward. The first nonaxisymmetric convection occurs at Rayleigh numbers which increase with rotation rate. The form of the convection near the transition also depends upon the rate of rotation. Selected flow patterns near the transition as well as those beyond it are studied numerically. For those cases where there exist laboratory experiments with which to compare, the numerical and experimental results agree very well.

Leslie, Fred W.↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

NASA NTRS1974

A zero-gravity demonstration of the collision and coalescence of water droplets

The mechanics of the collision and coalescence of liquid droplets is one of the main research areas in the fields of nuclear physics, astrophysics, meteorology and fluid mechanics. The crew members on the Skylab 3 and 4 missions were requested to perform demonstrations of the collision and coalescence of water droplets under the low gravity environment at orbital altitude. In Skylab 4 two water droplets with equal volumes, 30 cu cm each, were used. A dark colored droplet (contaminated with grape drink) moving with a velocity of 3.14 cm/sec collided with a stationary pink colored droplet (contaminated with strawberry drink) and coalescence occurred. Theoretical models are proposed to study the various stages of the collision-coalescence processes. Special considerations are concentrated in the investigation of the bounce-coalescence and coalescence-instability processes. The surface tension of the coalesced droplets was calculated to be 52 dynes/cm in perfect agreement with laboratory measurements made after the flight using a reproduction of the liquids.

Hung, R. J.↗

NASA NTRS1996

Performance Analysis of Two Early NACA High Speed Propellers with Application to Civil Tiltrotor Configurations

The helicopter industry is vigorously pursuing development of civil tiltrotors. One key to efficient high speed performance of this rotorcraft is prop-rotor performance. Of equal, if not greater, importance is assurance that the flight envelope is free of aeroelastic instabilities well beyond currently envisioned cuise speeds. This later condition requires study at helical tip Match numbers well in excess of 1.0. Two 1940's 'supersonic' propeller experiments conducted by NACA have provided an immensely valuable data bank with which to study prop-rotor behavior at transonic and supersonic helical tip Mach numbers. Very accurate 'blades alone' data were obtained by using nearly an infinite hub. Tabulated data were recreated from the many thrust and power figures and are included in two Appendices to this report. This data set is exceptionally well suited to re-evaluating classical blade element theories as well as evolving computational fluid dynamic (CFD) analyses. A limited comparison of one propeller's experimental results to a modem rotorcraft CFD code is made. This code, referred to as TURNS, gives very encouraging results. Detailed analysis of the performance data from both propellers is provided in Appendix A. This appendix quantifies the minimum power required to produce usable prop-rotor thrust. The dependence of minimum profile power on Reynolds number is quantified. First order compressibility power losses are quantified as well and a first approximation to design air-foil thickness ratio to avoid compressibility losses is provided. Appendix A's results are applied to study high speed civil tiltrotor cruise performance. Predicted tiltrotor performance is compared to two turboprop commercial transports. The comparison shows that there is no fundamental aerodynamic reason why the rotorcraft industry could not develop civil tiltrotor aircraft which have competitive cruise performance with today's regional, turboprop airlines. Recommendations for future study that will insure efficient prop-rotor performance to well beyond 400 knots are given.

Harris, Franklin D.↗

NASA NTRS1999

The Considere Condition and Rapid Stretching of Linear and Branched Polymer Melts

We analyze the onset of "necking" and subsequent filament failure during the transient uniaxial elongation of viscoelastic fluid samples in extensional rheometers. In the limit of rapid elongation (such that no molecular relaxation occurs), the external work applied is all stored elastically and the Considere criterion originally developed in solid mechanics can be used to quantitatively predict the critical Hencky strain to failure. By comparing the predictions of the Doi-Edwards model for linear homopolymer melts with those of the "Pom-Pom" model for prototypical branched melts we show that the critical strain to failure in rapid elongation of a rubbery material is intimately linked to the molecular topology of the chain, especially the degree of chain branching. The onset of necking instability is monotonically shifted to larger Hencky strains as the number of branches is increased. Numerical computations at finite Deborah numbers also show that there is an optimal range of deformation rates over which homogeneous extensions can be maintained to large strain. We also consider other rapid homogeneous stretching deformations, such as biaxial and planar stretching, and show that the degree of stabilization afforded by inclusion of material with long-chain branching is a sensitive function of the imposed mode of deformation.

McKinley, Gareth H.↗

NASA NTRS1962

STABILITY OF FILM BOILING TWO-PHASE FLOW IN CRYOGENIC SYSTEMS

The stability of film boiling two phase flow is investigated referring in particular to saturated cryogenic liquids. This type of fluid motion is unstable, such that a range of wave numbers will always be present in which disturbances are amplified. Further the mode of operation of the energy supply system may favour temperature fluctuations of large amplitude if the circuit does not provide feed-back compensation. The theoretically predicted stabilization due to capillarity is reduced considerably in liquefied gases with low surface tension, such as hydrogen. This effect causes the transition to turbulent flow to occur at low Reynolds numbers, and increases the heat transport in the turbulent regime. A comparison with experimental data, (horizontal cylinder), shows that the "critical" film Reynolds number ūδ/𝜈, at which instabilities become effective, is lowered from approximately 200, (O2), to about 30, (H2). Similar results have been found with other configurations. Therefore application of laminar film boiling solutions is restricted to small objects, entrance regions, or stagnation points. Turbulent correlations given so far are valid only for a certain group of substances with similar behavior.

LIQUID FLOW↗

NASA NTRSDate not supplied

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

Detonation↗

NASA NTRSDate not supplied

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

detonation↗

NASA NTRS2003

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗

NASA NTRS2002

Damping of Quasi-stationary Waves Between Two Miscible Liquids

Two viscous miscible liquids with an initially sharp interface oriented vertically inside a cavity become unstable against oscillatory external forcing due to Kelvin-Helmholtz instability. The instability causes growth of quasi-stationary (q-s) waves at the interface between the two liquids. We examine computationally the dynamics of a four-mode q-s wave, for a fixed energy input, when one of the components of the external forcing is suddenly ceased. The external forcing consists of a steady and oscillatory component as realizable in a microgravity environment. Results show that when there is a jump discontinuity in the oscillatory excitation that produced the four-mode q-s wave, the interface does not return to its equilibrium position, the structure of the q-s wave remains imbedded between the two fluids over a long time scale. The damping characteristics of the q-s wave from the time history of the velocity field show overdamped and critically damped response; there is no underdamped oscillation as the flow field approaches steady state. Viscous effects serve as a dissipative mechanism to effectively damp the system. The stability of the four-mode q-s wave is dependent on both a geometric length scale as well as the level of background steady acceleration.

Duval, Walter M. B.↗

NASA NTRS1982

Vortex simulation of an inviscid shear layer

The accuracy of the vortex-blob method was tested by simulating a free-shear-layer instability, Kirchhoff's elliptical vortex, and a circular vortex. The main numerical parameters in the vortex-blob method are the density of the vortices, and the distribution of vorticity within each vortex core. The growth rate of a periodic unstable mode of the shear layer was calculated numerically and compared with the exact result. The error is only a few percent for about 10 rows of vortex blobs. The error is reduced by decreasing the spacing between vortices and, correspondingly, the core size. In the simulation of the motion of the elliptical vortex, the rotation of the boundary, without change of shape, and the circular particle paths of the vortical fluid were well simulated. For the circular vortex, optimum sets of parameters were obtained by comparing them with the exact velocity. The results are consistent with convergence theories of the vortex-blob method. In particular, second-order convergence is observed with a Gaussian core from velocity calculation.

Nakamura, Y.↗

NASA NTRS1992

Fluid flow in a jet and the Ca XIX line profiles observed during solar flares

The present study considers two types of fluid models as methods for reproducing the blueshifts and line broadening observed in soft X-ray lines during the rise of solar flares. The fluid models representing the laminar flow of material in a jet and through a pipe were used to derive the velocity at each location in the flow. The optimal values for the free parameters in the jet model were within the ranges expected for solar flare conditions. The jet was much smaller than typical flare volumes estimated from spatial images, but the resulting densities were consistent with the upper limit of densities derived using density-sensitive line ratios. This supports previous reports of small filling factors and is interpreted as evidence for a chromospheric origin for the flow. Stability analysis indicates that, although the laminar jet can reproduce the observations, instabilities may cause the flow to become turbulent, and that this laminar treatment is not strictly valid.

Bornmann, P. L.↗

NASA NTRSDate not supplied

Full Stall Simulations of a Redesigned Ventilation Fan for the ISS

The concept of a stall is studied rigorously in the aerospace industry. From a design standpoint, instabilities such as stall are undesirable– operation in the stall regime has a tremendous impact on aerodynamic performance as well as structural integrity. In extreme cases, operating in stall conditions can cause failure. Further, stall can cause a loss of lift on aircraft wings or a loss of thrust in aircraft engines. In any case, stall continues to be a topic of interest in the aerospace industry. The present work aims to analyze the stall characteristics of a ventilation fan that was recently designed for the International Space Station (ISS). Although the ventilation fan has a rotor-stator design, this paper considers a rotor-only configuration. The FUN3D Computational Fluid Dynamic (CFD) solver developed by NASA Langley Research Center was used to simulate the operational characteristics of the ventilation fan. FUN3D solves the Unsteady Reynolds-Averaged Naiver-Stokes (URANS) equations using implicit time marching and a dynamic overset grid. The FUN3D solver was originally written for exterior flow fields; however, this work represents an extension of the FUN3D solver to turbomachinery or interior flow fields. The FUN3D results for the rotor-only ventilation fan accurately captured the operational characteristics inherent to compressors– the results shared similar performance trends when compared to the experimental results for the rotor-stator case. A peak adiabatic efficiency of 96% occurred at a MFR of 105.2 CFM, the minimum aerodynamically stable point. A computationally stable stall occurred at a mass flow rate of 43.8 CFM where the adiabatic efficiency dropped to 69%.

CFD↗

NASA NTRS2020

Flow Simulations for Martian and Lunar Lander Plume-Surface Interaction Prediction

Landing vehicle propulsive systems introduce several risks associated with plume/surface interactions (PSI) during the final stages of descent. For example, view obscuration during descent due to the liberation of dust driven by the rocket plume’s interaction with the lunar surface was reported on several Apollo missions. Plume-driven debris impacts on vehicle and nearby assets as well as landing instabilities due to crater formation are also potential risks associated with plume/surface interactions. The importance of studying these plume/surface interactions has been elevated due to NASA’s mission to land humans on the Moon by 2024. Existing and under-development predictive simulation tools are currently being brought to application readiness under a project awarded by the NASA Space Technology Mission Directorate (STMD) Game Changing Development (GCD) program. This tool development is divided into four tasks: (1) simulating plume structure in low-pressure environments, (2) simulating crater development and ejecta, (3) regolith particle phase modeling, and (4) gas-particle interaction modeling. Recent efforts under task 1 include demonstrating existing computational fluid dynamics (CFD) tools for simulating plume structure in low pressure environments similar to those found on Mars. The results of these demonstrative simulations are compared to available experimental data. Emerging capability being developed under task 2 is also presented, demonstrating simulations of crater formation in Lunar conditions.

Thomas P Shurtz↗

NASA NTRS2015

Advancing the State-of-the-Practice for Liquid Rocket Engine Injector Design

Current shortcomings in both the overall injector design process and its underlying combustion stability assessment methodology are rooted in the use of empirically based or low fidelity representations of complex physical phenomena and geometry details that have first order effects on performance, thermal environments and combustion stability. The result is a design and analysis capability that is often inadequate to reliably arrive at a suitable injector design in an efficient manner. Specifically, combustion instability has been particularly difficult to predict and mitigate. Large hydrocarbon-fueled booster engines have been especially problematic in this regard. Where combustion instability has been a problem, costly and time-consuming redesign efforts have often been an unfortunate consequence. This paper presents an overview of a recently completed effort at NASA Marshall Space Flight Center to advance the state-of-the-practice for liquid rocket engine injector design. Multiple perturbations of a gas-centered swirl coaxial (GCSC) element that burned gaseous oxygen and RP-1 were designed, assessed for combustion stability, and tested. Three designs, one stable, one marginally unstable and one unstable, were used to demonstrate both an enhanced overall injector design process and an improved combustion stability assessment process. High-fidelity results from state-of-the-art computational fluid dynamics CFD simulations were used to substantially augment and improve the injector design methodology. The CFD results were used to inform and guide the overall injector design process. They were also used to upgrade selected empirical or low-dimensional quantities in the ROCket Combustor Interactive Design (ROCCID) stability assessment tool. Hot fire single element injector testing was used to verify both the overall injector designs and the stability assessments. Testing was conducted at the Air Force Research Laboratory and at Purdue University. Companion papers provide details of the overall injector design process, full- and sub-scale testing, ROCCID-based stability assessments and the CFD simulations.

Tucker, P. K.↗