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.
Electrohydrodynamic Rayleigh-Taylor bulk instability
Stability criterion for initially static, stratified fluid subject to electrical stress
Quantifying Instability Sources in Liquid Rocket Engines
Computational fluid dynamics methodology to predict the effects of combusting flows on acoustic pressure oscillations in liquid rocket engines (LREs) is under development. 'Me intent of the investigation is to develop the causal physics of combustion driven acoustic resonances in LREs. The crux of the analysis is the accurate simulation of pressure/density/sound speed in a combustor which when used by the FDNS-RFV CFD code will produce realistic flow phenomena. An analysis of a gas generator considered for the Fastrac engine will be used as a test validation case.
Surface Instability of Liquid Propellant under Vertical Oscillatory Forcing
Fluid motion in a fuel tank produced during thrust oscillations can circulate sub-cooled hydrogen near the liquid-vapor interface resulting in increased condensation and ullage pressure collapse. The first objective of this study is to validate the capabilities of a Computational Fluid Dynamics (CFD) tool, CFD-ACE+, in modeling the fundamental interface transition physics occurring at the propellant surface. The second objective is to use the tool to assess the effects of thrust oscillations on surface dynamics. Our technical approach is to first verify the CFD code against known theoretical solutions, and then validate against existing experiments for small scale tanks and a range of transition regimes. A 2D axisymmetric, multi-phase model of gases, liquids, and solids is used to verify that CFD-ACE+ is capable of modeling fluid-structure interaction and system resonance in a typical thrust oscillation environment. Then, the 3D mode is studied with an assumed oscillatory body force to simulate the thrust oscillating effect. The study showed that CFD modeling can capture all of the transition physics from solid body motion to standing surface wave and to droplet ejection from liquid-gas interface. Unlike the analytical solutions established during the 1960 s, CFD modeling is not limited to the small amplitude regime. It can extend solutions to the nonlinear regime to determine the amplitude of surface waves after the onset of instability. The present simulation also demonstrated consistent trends from numerical experiments through variation of physical properties from low viscous fluid to high viscous fluids, and through variation of geometry and input forcing functions. A comparison of surface wave patterns under various forcing frequencies and amplitudes showed good agreement with experimental observations. It is concluded that thrust oscillations can cause droplet formation at the interface, which results in increased surface area and enhanced heat transfer between the liquid and gas phases as the ejected droplets travel well into the warmer gas region.
Dissipation of plasma oscillations excited in a current-carrying plasma
Heating, dissipation, and containment of plasma in strong magnetic fluid, under conditions where bunched instabilities excited by escaping electrons develop
Electrohydrodynamic Rayleigh-Taylor bulk instability.
Electrohydrodynamic Rayleigh-Taylor bulk instability in initially static stratified fluid under electric stress
Nonthermal turbulent heating in the solar envelope.
It is shown that MHD pulses, in the form of fast magnetosonic waves or solitons, can produce a strong electron-ion coupling capable of maintaining electron-proton temperature equilibrium in the solar envelope. The mechanism producing the nonthermal heating is the fluid-like modified two-stream instability, which, since it is essentially independent of the electron-proton temperature ratio and the value of beta, becomes a prime candidate for the anomalous collisions required by the fluid models inside a distance less than 30 solar radii, in order to explain the dominant features of the solar-wind flow.
Electrostatic Kelvin-Helmholtz instability produced by a localized electric field perpendicular to an external magnetic field
The evolution of a plasma with a localized electric field perpendicular to an external magnetic field is shown to be dominated by the Kelvin-Helmholtz instability. For small ion gyroradius, the instability is similar to the fluid mode. When the ion gyroradius is an appreciable fraction of the spatial extent of the electric field, the plasma is not in equilibrium, and the initial drift profile relaxes. Subsequent evolution still leads to vortex flows.
Self-similar evolution of the nonlinear magnetic buoyancy instability
A new type of self-similar solution of ideal magnetohydrodynamics (MHD) in the nonlinear stage of the undular model (k parallel to B) of the magnetic buoyancy instability (the ballooning instability in fusion plasma physics or the Parker instability in astrophysics) is found through MHD simulation and theory. The linear theory developed agrees well with the simulation in the early (linear) stage. The nonlinear stages of the instability in the simulation show the self-similar evolution. One of the solutions obtained from the nonlinear analysis has the characteristics of nonlinear instability in Lagrangian coordinates; the fluid velocity and the Alfven speed on each magnetic loop increase exponentially with time, because the loop is evacuated by the field-aligned motion of matter resulting from gravitational acceleration. In the later stage of the nonlinear evolution, the solution property changes from exponential to power-law time dependence. The latter corresponds to a force-free expansion solution. The later saturation of the velocity increment is also discussed.
Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence
When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.
Influence of Shock Wave on the Flutter Behavior of Fan Blades Investigated
Modern fan designs have blades with forward sweep; a lean, thin cross section; and a wide chord to improve performance and reduce noise. These geometric features coupled with the presence of a shock wave can lead to flutter instability. Flutter is a self-excited dynamic instability arising because of fluid-structure interaction, which causes the energy from the surrounding fluid to be extracted by the vibrating structure. An in-flight occurrence of flutter could be catastrophic and is a significant design issue for rotor blades in gas turbines. Understanding the flutter behavior and the influence of flow features on flutter will lead to a better and safer design. An aeroelastic analysis code, TURBO, has been developed and validated for flutter calculations at the NASA Glenn Research Center. The code has been used to understand the occurrence of flutter in a forward-swept fan design. The forward-swept fan, which consists of 22 inserted blades, encountered flutter during wind tunnel tests at part speed conditions.
Investigation of Combustion Control in a Dump Combustor Using the Feedback Free Fluidic Oscillator
A feedback free fluidic oscillator was designed and integrated into a single element rocket combustor with the goal of suppressing longitudinal combustion instabilities. The fluidic oscillator uses internal fluid dynamics to create an unsteady outlet jet at a specific frequency. An array of nine fluidic oscillators was tested to mimic modulated secondary oxidizer injection into the combustor dump plane. The combustor has a coaxial injector that uses gaseous methane and decomposed hydrogen peroxide with an overall O/F ratio of 11.7. A sonic choke plate on an actuator arm allows for continuous adjustment of the oxidizer post acoustics enabling the study of a variety of instability magnitudes. The fluidic oscillator unsteady outlet jet performance is compared against equivalent steady jet injection and a baseline design with no secondary oxidizer injection. At the most unstable operating conditions, the unsteady outlet jet saw a 67% reduction in the instability pressure oscillation magnitude when compared to the steady jet and baseline data. Additionally, computational fluid dynamics analysis of the combustor gives insight into the flow field interaction of the fluidic oscillators. The results indicate that open loop high frequency propellant modulation for combustion control can be achieved through fluidic devices that require no moving parts or electrical power to operate.
Convective instability when the temperature gradient and rotation vector are oblique to gravity. I - Fluids without diffusion
A linear stability analysis of fluid layers under uniform rotation (generally oblique to gravity) which possess both vertical and horizontal temperature gradients is made by considering ideal fluids without diffusion within a Boussinesq approximation. This simplified configuration is used to assess the preferred convective modes as a function of latitude on a planet like Jupiter. The tilted rotation vector introduces a preference for roll-like disturbances with north-south orientations, while the horizontal temperature gradient produces a thermal wind shear which favors convective rolls oriented parallel to the flow in an east-west direction. It is found that the horizontal temperature gradient needed to produce a preference for the axisymmetric or east-west rolls increases with an increasing rotation rate and a decreasing latitude. The parameter values for Jupiter are estimated with the use of a simple radiative convective model, indicating a preference for axisymmetric rolls at nearly all latitudes if the convection zone depth is greater than about 200 km below the one atmosphere pressure level and convective roll characteristics which contribute to an equatorial acceleration.
Numerical simulation of the non-Newtonian mixing layer
This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.
International Symposium on Microgravity Science and Application (ISMSA), Beijing, China, May 10-13, 1993
Papers at this conference on Microgravity Science and Technology concerned fluid science and biotechnology. Some representative titles include the following: Hydrodynamic Instabilities in Thermocapillary Flow; Marangoni Convection in Immiscible Double Liquid Layers; Note on the Instability of Capillary Jet with Thermocapillarity; Residual Gravity Jitter Effects on Fluid Processes; Responsive Motion of Bubbles to Periodic g-jitter; and Protein Crystallization in Space.
Kinetic simulations of the Kelvin-Helmholtz instability at the magnetopause
Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.
A general theory for dynamic instability of tube arrays in crossflow
A general theory of fluidelastic instability for a tube array in crossflow is presented. Various techniques to obtain the motion-dependent fluid-force coefficients are discussed and the general instability characteristics are summarized. The theory is also used to evaluate the results of other mathematical models for crossflow-induced instability.
Electrohydrodynamic charge relaxation and interfacial perpendicular-field instability.
Small amplitude motions of plane interface between fluids stressed by initially perpendicular electric field, modeling fluids as ohmic conductors