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 973 records · Page 54

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Skylab experiment M-092 - Results of the first manned mission

Blood pressure, heart rate, and percentage increase in leg volume were continuously recorded for the Inflight Lower Body Negative Pressure (LBNP) experiment conducted during the first manned Skylab mission. Data were collected over a 5-month preflight period as well as at approximately 3-day intervals throughout the mission. Individual variations in cardiovascular responses to LBNP during the preflight period continued to be demonstrated inflight, and measurements of the calf showed that a large volume of fluid was shifted out of the legs early in the flight. A much greater increase in leg volume occurred during inflight LBNP than in preflight tests. Resting heart rates tended to be low early in the flight and to increase slightly as it progressed, while resting blood pressure varied. The LBNP protocol was a greater stress inflight, and the tests had to be stopped early on three occasions due to impending syncopal reactions. Inflight LBNP responses seemed to predict the degree of postflight orthostatic intolerance. Postflight responses to LBNP during the first 48 hours were characterized by marked elevations in heart rate, instability in blood pressure, and considerably elevated systolic and diastolic pressures both at rest and during stress.

Johnson, R. L.↗

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. 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. 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. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. 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, slip boundary conditions, and the supersonic 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, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the 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 one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. 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. 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, S. R.↗

NASA NTRS1988

Inflow/outflow boundary conditions and global dynamics of spatial mixing layers

The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed.

Buell, J. C.↗

NASA NTRS1993

Finite element calculations of transonic flutter in cascades

Aeroelastic stability and response of turbomachinery cascades operating in the transonic regime are explored using a new computational approach. The two-dimensional unsteady Euler equations are used to model the transonic cascade flows, and the cascade blades are modeled as thin plates. To minimize the discretization incompatibility, both the fluid and the solid domains are discretized using finite element procedures. The fluid and the structural models are coupled at the element level by introducing a channel theory, which is analogous to a strip theory in the classical sense. The system of space-discretized equations is marched forward in time using an explicit four-stage Runge-Kutta scheme. Results of free vibration analyses of typical blades are presented and compared to published finite element and Rayleigh-Ritz results. Aeroelastic stability and response characteristics of several transonic cascades are investigated using the new quasi-3D model. These results demonstrate the present 3D finite element blade model is capable of capturing aeroelastic instabilities involving higher modes, including (plate-type) camber bending modes.

Hsiao, Chingteng↗

NASA NTRS1983

On the stability of the unsteady boundary layer on a cylinder oscillating transversely in a viscous fluid

The stability of the two-dimensional flow induced by the tranverse oscillation of a cylinder in a viscous fluid is investigated in both the linear and weakly nonlinear regime. The major assumption that is made to simplify the problem is that the oscillation frequency is large in which case an unsteady boundary layer is set up on the cylinder. Results are given for cylinders of elliptic cross section and it is found that for any given eccentricity the most dangerous configuration is when the cylinder oscillates parallel to its minor axis. Some discussion of nonlinear effects is also given and for the circular cylinder it is shown that the steady streaming boundary layer of the basic flow is significantly altered by the instability.

Hall, P.↗

NASA NTRS1984

On the stability of the unsteady boundary layer on a cylinder oscillating transversely in a viscous fluid

The stability of the two-dimensional flow induced by the transverse oscillation of a cylinder in a viscous fluid is investigated in both the linear and weakly nonlinear regime. The major assumption that is made to simplify the problem is that the oscillation frequency is large in which case an unsteady boundary layer is set up on the cylinder. Results are given for cylinders of elliptic cross section and it is found that for any given eccentricity the most dangerous configuration is when the cylinder oscillates parallel to its minor axis. Some discussion of nonlinear effects is also given and for the circular cylinder it is shown that the steady streaming boundary layer of the basic flow is significantly altered by the instability.

Hall, P.↗

NASA NTRS1987

Spacelab experiments on convection in a rotating spherical shell with radial gravity

Data from the geophysical fluid flow cell experiment on Spacelab 3 revealed the diverse forms of fluid motion that can occur in geometrically simple models of global convection with rotation and radial gravity. The following types of convection were observed: midlatitude waves interacting with the low-latitude columnar convection rolls or 'banana cells', spiral waves' near the poles when latitudinal heating gradients are present on the bounding surfaces, and 'triangular waves' coupling midlatitude and equatorial disturbances under similar differential heating. It is believed that concepts based on the Taylor-Proudman theorem have a central role in many flows realized in the rapidly rotating experiments when H is approximately 0; however, the resulting banana cells are subject to secondary instabilities that produce intricate time dependence and eventually turbulent flows in which only hints of the simpler patterns are recognizable.

Toomre, Juri↗

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↗