Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid instabilities”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Convective instability when the temperature gradient and rotation vector are oblique to gravity. II - Real fluids with effects of diffusion

The linear stability analysis of Hathaway et al. (1979) is repeated for Boussinesq fluids with viscous and thermal diffusion. The fluid is confined between plane parallel boundaries and the rotation vector is oblique to gravity; this tilted rotation vector introduces a preference for roll-like disturbances whose axes are oriented north-south. The presence of a latitudinal temperature gradient produces a thermal wind shear which favors axisymmetric convective rolls if the gradient exceeds some critical value. The axisymmetric rolls are similar in many respects to the cloud bands on Jupiter provided they extend to a depth of about 15,000 km.

Hathaway, D. H.↗
NASA NTRS1999

Electric Field Induced Interfacial Instabilities

The study of the interface in a charge-free, critical and near-critical binary fluid in the presence of an externally applied electric field is presented. At sufficiently large fields, the interface between the two phases of the binary fluid should become unstable and exhibit an undulation with a predefined wavelength on the order of the capillary length. As the critical point is approached, this wavelength is reduced, potentially approaching length-scales such as the correlation length or critical nucleation radius. At this point the critical properties of the system may be affected. In this paper, the flat interface of a marginally polar binary fluid mixture is stressed by a perpendicular alternating electric field and the resulting instability is characterized by the critical electric field E(sub c) and the pattern observed. The character of the surface dynamics at the onset of instability is found to be strongly dependent on the frequency f of the field applied. The plot of E(sub c) vs. f for a fixed temperature shows a sigmoidal shape, whose low and high frequency limits are well described by a power-law relationship, E(sub c) = epsilon(exp zeta) with zeta = 0.35 and zeta = 0.08, respectively. The low-limit exponent compares well with the value zeta = 4 for a system of conducting and non-conducting fluids. On the other hand, the high-limit exponent coincides with what was first predicted by Onuki. The instability manifests itself as the conducting phase penetrates the non-conducting phase. As the frequency increases, the shape of the pattern changes from an array of bifurcating strings to an array of column-like (or rod-like) protrusions, each of which spans the space between the plane interface and one of the electrodes. For an extremely high frequency, the disturbance quickly grows into a parabolic cone pointing toward the upper plate. As a result, the interface itself changes its shape from that of a plane to that of a high sloping pyramid.

Kusner, Robert E.↗
NASA NTRS1989

Nonlinear generation of whistler waves by an ion beam

An electromagnetic hybrid code is used to simulate a new mechanism for whistler wave generation by an ion beam. First, a field-aligned ion beam becomes unstable to the electromagnetic ion/ion right-hand resonant instability which generates large amplitude MHD-like waves. These waves then trap the ion beam and increase its effective temperature anisotropy. As a result, the growth rates of the electron/whistler instability are significantly enhanced, and whistlers start to grow above the noise level. At the same time, because of the reduced parallel drift speed of the ion beam, the frequencies of the whistlers are also downshifted. Full simulations were performed to isolate and separately investigate the electron/ion whistler instability. The results are in agreement with the assumption of fluid electrons in the hybrid simulations and with the linear theory of the instability.

Akimoto, K.↗
NASA NTRS1988

Effects of eletron heating on the current driven electrostatic ion cyclotron instability and plasma transport processes along auroral field lines

Fluid simulations of the plasma along auroral field lines in the return current region have been performed. It is shown that the onset of electrostatic ion cyclotron (EIC) related anomalous resistivity and the consequent heating of electrons leads to a transverse ion temperature that is much higher than that produced by the current driven EIC instability (CDICI) alone. Two processes are presented for the enhancement of ion heating by anomalous resistivity. The anomalous resistivity associated with the turbulence is limited by electron heating, so that CDICI saturates at transverse temperature that is substantially higher than in the absence of resistivity. It is suggested that this process demonstrates a positive feedback loop in the interaction between CDICI, anomalous resistivity, and parallel large-scale dynamics in the topside ionosphere.

Ganguli, Supriya B.↗
NASA NTRS1985

Magnetic buoyancy instabilities incorporating rotation

Plane layer stability analysis is extended to incorporate the effects of uniform rotation. Detailed studies are made of the interchange, or 'axisymmetric modes' and of the undular or wavelike motions, considering both high and low frequency modes. The force due to rotation on the fluid motions are determined. The types of instability relevant to astrophysical problems, are also discussed with attention given to the occurrence of two distinct modes of instability in a bottom heavy field gradient. Some peculiar stability boundaries, which are due to the absence of a dominating force, are described.

Hughes, D. W.↗
NASA NTRS1990

Interfacial dynamics of two liquids under an oscillating gravitational field

The evolution of two miscible liquids meeting at an initially sharp interface inside a cavity under microgravity g-jitter conditions is studied numerically. The response of the interface and kinematics of the flowfield to various g-jitter accelerations and aspect ratio variations is shown. The interface region acts like a vortex source sheet, and it can be unstable to Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The vortices produced along the interface can serve as a stirring mechanism to promote local mixing. Below the critical Stokes-Reynolds number, the destabilization of the interface results in deformation into wavy structures. In some parameter regions, these structures oscillate in time; in others, they are quasi-steady. Above the critical Stokes-Reynolds number, 'chaotic' instability results, and the interface breaks into concentration pockets. The morphology of the initial breakup is similar to observed Rayleigh-Taylor instability. Subsequent mixing of the two fluids after the breakup of the interface is then very rapid.

Duval, Walter M. B.↗
NASA NTRS1991

Nonaxisymmetric instabilities in rotating shear flows - Internal gravity modes in stratified media and analogies with plane flows

The role of the internal gravity modes in mediating the growth of nonaxisymmetric instabilities is investigated by studying the instability of stratified incompressible differentially rotating fluid cylinders to global nonaxisymmetric modes. The results indicate that, in addition to a modified version of the well-known principal branch mediated by surface modes of the system (analogous to f-modes of stars), there exist unstable branches of the dispersion relation mediated by internal gravity modes of the system (similar to the g-modes of stars). These branches arise due to the interaction between the g-modes. It is shown that the maximum growth rate on one of the new branches can sometimes equal or exceed that on the principal branch, thus modifying the principal branch.

Ghosh, P.↗
NASA NTRS1992

Anomalous ion mixing within a Kelvin-Helmholtz vortex in a collisionless plasma

Anomalously fast ion mixing is observed in a hybrid code simulation (particle ions and a massless electron fluid) of the Kelvin-Helmholtz instability in a collisionless plasma. While the traditional view predicts that the ion mixing occurs on a time scale longer than the roll-up time T(v) of the Kelvin-Helmholtz vortex by a factor of (velocity shear scale length)/(ion Larmor radius), the simulation results show that the mixing occurs within a time scale comparable to T(v). It is shown that this fast mixing is due to the scattering of ions by unsteady electromagnetic fields.

Terasawa, T.↗
NASA NTRS2004

Torque Balances on the Taylor Cylinders in the Geomagnetic Data Assimilation

In this presentation we report on our continuing effort in geomagnetic data assimilation, aiming at understanding and predicting geomagnetic secular variation on decadal time scales. In particular, we focus on the effect of the torque balances on the cylindrical surfaces in the core co-axial with the Earth's rotation axis (the Taylor cylinders) on the time evolution of assimilated solutions. We use our MoSST core dynamics,model and observed geomagnetic field at the Earth's surface derived via Comprehensive Field Model (CFM) for the geomagnetic data assimilation. In our earlier studies, a model solution is selected randomly from our numerical database. It is then assimilated with the observations such that the poloidal field possesses the same field tomography on the core-mantel boundary (CMB) continued downward from surface observations. This tomography change is assumed to be effective through out the outer core. While this approach allows rapid convergence between model solutions and the observations, it also generates sevee numerical instabilities: the delicate balance between weak fluid inertia and the magnetic torques on the Taylor cylinders are completely altered. Consequently, the assimilated solution diverges quickly (in approximately 10% of the magnetic free-decay time in the core). To improve the assimilation, we propose a partial penetration of the assimilation from the CMB: The full-scale modification at the CMB decreases linearly and vanish at an interior radius r(sub a). We shall examine from our assimilation tests possible relationships between the convergence rate of the model solutions to observations and the cut-off radius r(sub a). A better assimilation shall serve our nudging tests in near future.

Kuang, Weijia↗
NASA NTRS2000

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗
NASA NTRS2002

Phase Transition in Dusty Plasmas: A Microphysical Description

Dust grains immersed in plasma discharges acquire a large negative charge and settle into a dust cloud at the edge of the sheath. In this region, the plasma ions stream toward the electrode at a velocity u approx. cs=(T(sub e)/m(sub i))(exp 1/2). Experimentally at sufficiently high gas pressure P, the random kinetic energy of the grains is damped by gas friction, and the grains are strongly coupled and self-organize into a crystalline configuration. For lower pressures despite the dissipation of grain kinetic energy to gas friction, the dust grains reach a steady-state kinetic temperature T(sub d) which is much larger than the temperature of any other component in the plasma. T(sub d) is so large that the dust acts like a fluid. We have used the dynamically shielded dust (DSD) model to simulate these physical processes. We find that the known experimental features are nicely reproduced in the simulations, and that additional features are revealed. In the figure we plot the variation of T(sub d) as P is continuously varied in a DSD code run. A marked difference is evident between the critical pressure P(sub m) for the melting transition as P is decreased, and the critical pressure P(sub c) for the condensation transition as P is increased. For P(sub m) is less than P is less than P(sub c), mixed phase states are seen. This hysteresis occurs because the instability which triggers melting is different from the instability that heats the dust in the fluid phase and inhibits freezing. At low pressure, the dust is subject to a two-stream instability with the ions. This instability is responsible for the high temperature of the dust at low pressure. The basic physics underlying the melting transition has been elucidated in a series of papers. We are developing a first-principles analytic approach to the melting transition, which embodies the same physics that is present in the DSD code.

Joyce, Glenn↗
NASA NTRS2019

Operational Stability Limits in Rotating Detonation Engine Numerical Simulations

An instability is described which arises in computational fluid dynamic (CFD) simulations of semi-idealized rotating detonation engines (RDE) configured with a throat at the exit. Its existence is verified by examining output from two independently developed CFD codes simulating the same configuration and producing solutions that agree well. The instability is shown to be thermo-acoustic in that a spatial integral of the product of pressure and heat release fluctuations develops a regular oscillation which grows in time. The instability can become severe enough to cause detonation failure. Its onset is shown to be closely linked to the size of the exit throat and the size of the inlet restriction; both parameters that strongly influence RDE performance. It is shown that the instability places a cap on ideal RDEperformance, but that an optimized exhaust throat and inlet restriction combination still yields substantial pressure gain. Other parametric sensitivities are also examined in terms ofinstability growth. These include axial length, inlet manifold pressure, and air-fuel ratio.

Paxson, Daniel E.↗
NASA NTRS1996

Electric Field Induced Interfacial Instabilities

The study of the interface in a charge-free, nonpolar, critical and near-critical binary fluid in the presence of an externally applied electric field is presented. At sufficiently large fields, the interface between the two phases of the binary fluid should become unstable and exhibit an undulation with a predefined wavelength on the order of the capillary length. As the critical point is approached, this wavelength is reduced, potentially approaching length-scales such as the correlation length or critical nucleation radius. At this point the critical properties of the system may be affected. In zero gravity, the interface is unstable at all long wavelengths in the presence of a field applied across it. It is conjectured that this will cause the binary fluid to break up into domains small enough to be outside the instability condition. The resulting pattern formation, and the effects on the critical properties as the domains approach the correlation length are of acute interest. With direct observation, laser light scattering, and interferometry, the phenomena can be probed to gain further understanding of interfacial instabilities and the pattern formation which results, and dimensional crossover in critical systems as the critical fluctuations in a particular direction are suppressed by external forces.

Kusner, Robert E.↗
NASA NTRS1985

Whirl/whip demonstration

Fluid flow in bearings and seals, set in motion by shaft rotation, generates dynamic forces which may result in a well recognized instability known as whirl and whip. These are lateral, forward precessional, self excited, subsynchronous vibrations in which the amplitude may vary from very small to nearly the limit of the bearing or seal clearances. Oil whirl in lubricated bearings, in particular, typically occurs at somewhat less than half rotative speed. As the rotative speed increases, the frequency relationship remains constant until the whirl frequency approaches the first balance resonance. Now the whirl is smoothly replaced by whip at a nearly constant frequency asymptotically approaching first balance resonance, independent of increasing rotative speed. Changes in bearing/seal radial loading can permit, prevent, or eliminate this instability. The oil whirl/whip rig demonstrates the effects of fluid dynamic forces generated by the rotating shaft. At low rotative speeds, this produces changes of the journal static equilibrium position within the bearing. The demonstrator shows the relationship between any load direction and the average journal equilibrium position. At higher rotative speeds, the instability threshold is observed as a function of unidirectional radial load, unbalance, and rotor configuration.

Grissom, R.↗
NASA NTRS1990

Buoyancy-driven instabilities and the nonlinear breakup of a sheared magnetic layer

Motivated by problems concerning the storage and subsequent escape of the solar magnetic field, a study has been made of how a magnetic layer embedded in a convectively stable atmosphere evolves due to axisymmetric instabilities driven by magnetic buoyancy. The initial equilibrium consists of a toroidal field sheared by a weaker poloidal component. The linear stability problem is investigated for both ideal and resistive MHD, and the nonlinear evolution is followed by numerical integration of the equations of motion. In all cases, the instability is greatly affected by the distribution and strength of the poloidal field. In particular, both the horizontal and vertical scales of the motions are controlled by the location of the surface on which the poloidal field vanishes: the resonant surface. In the nonlinear regime, a resonant surface close to the interface between the magnetized and field-free fluid leads to the localization of the instability, so that only a fraction of the magnetic region is disrupted by the motions. By contrast, a deeply seated resonant surface leads to the complete disruption of the layer and to the formation of large, helical magnetic fragments whose identity is preserved for the entire simulation.

Cattaneo, Fausto↗
NASA NTRS1993

Kinetic physics of the mirror instability

Kinetic mechanisms for the growth and saturation of the mirror instability are described using one-dimensional hybrid simulations. Two parameter regimes are considered. In the first regime, a relatively small ion anisotropy excites a slowly growing instability that produces small-amplitude waves; most ions respond to the waves as an adiabatic fluid. In the second regime a large anisotropy excites a rapidly growing instability that generates large-amplitude waves; the response of many ions in this case is nonadiabatic. The difference in ion response is due to the relative importance of two ion populations, resonant and nonresonant. Resonant ions, those ions with low velocities parallel to the background magnetic field, contribute to the growth of the instability as a result of their gyrointeractions with the noncoplanar component of the waves electric field and respond to the mirror waves nonadiabatically. Nonreasonant ions, those with large parallel velocities, respond as an adiabatic fluid. In both regimes, ion anisotropy is reduced by means of the magnetic mirror force; in the second regime, the anisotropy reduction reduces the free energy available for wave growth and leads to the saturation of the mirror instability.

Mckean, M. E.↗
NASA NTRS1985

Remarks on rotor stability (A contribution to discussion during symposium on instability in rotating machinery)

Influence of nonlinear factors in fluid lubricated bearings on rotor vibrations is discussed as are computation of threshold of stability for systems with multiple degrees of freedom. The reason for this investigation was to learn more about the influence of tilting pad bearings on vertical rotors. The eigenvalues, the stability threshold, the response on the Heaviside-function and the unbalance response were calculated for a rigid rotor, a Jeffcott rotor and for the finite element model of a real rotor. The comparison of the linear and nonlinear calculations shows the limit of the linear calculations. The computation is supposed to show, that the model of the Jeffcott rotor is not sufficient in most practical cases to determine the stability threshold. This was known for some years, when computed results were compared with the behavior of real machines. That's why a program was developed, which is roughly described.

Kraemer, E.↗