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 379 records · Page 21

Bearing Tester Data Compilation Analysis, and Reporting and Bearing Math Modeling

The magnitude and direction of fluid induced torques and forces on the 57 mm bearing cage is considered to be a contributing factor in possible cage instabilities that can produce intermittent high heating in the bearing. Analyses of the fluid forces and torques are presented. Heat generated by viscous fluid work was estimated for two flow diverter configurations and a coolant flow of 10 lbs/sec to support the thermal evaluation of the LOX Bearing Materials Tester. Results of the analysis of the LOX turbopump turbine end bearings are discussed. Coolant velocities for the no. 4 LOX turbopump turbine end bearings were estimated as a function of shaft speed and coolant flow rate. Contact angles and track width data were developed for the 57 mm bearing as functions of shaft speed, and axial and radial loads. The Advanced Dynamics of Rolling Elements (ADORE) computer program was installed on the MSFC UNIVAC 1100 and a test case successfully run. Both the text output and the plotting output were verified. The Bearing Seal and Materials Tester - Test Condition Data Base was developed. The parametric analysis of the operating characteristics of the LOX turbopump pump end bearing using the 45 mm bearing thermal model was begun.

Source record↗

Filtering of non-linear instabilities

For Courant numbers larger than one and cell Reynolds numbers larger than two, oscillations and in some cases instabilities are typically found with implicit numerical solutions of the fluid dynamics equations. This behavior has sometimes been associated with the loss of diagonal dominance of the coefficient matrix. It is shown here that these problems can in fact be related to the choice of the spatial differences, with the resulting instability related to aliasing or nonlinear interaction. Appropriate 'filtering' can reduce the intensity of these oscillations and in some cases possibly eliminate the instability. These filtering procedures are equivalent to a weighted average of conservation and non-conservation differencing. The entire spectrum of filtered equations retains a three-point character as well as second-order spatial accuracy. Burgers equation has been considered as a model. Several filters are examined in detail, and smooth solutions have been obtained for extremely large cell Reynolds numbers.

Khosla, P. K.↗

The zero-frequency ion ring instability

The electrostatic zero-frequency ion ring instability with wave vector perpendicular to a uniform magnetic field B is examined through linear and second-order theory as well as by computer simulation. In the simulation ions are taken as magnetized particles; the electrons are described as a massless fluid subject to E x B motion. Saturation of the instability is primarily due to broadening of the ion ring distribution. A second-order theory provides an approximate criterion for the saturation amplitude, as does a simple trapping argument. Thus, for the simulation presented here, both quasi-linear and trapping effects contribute to saturation.

Gary, S. P.↗

On the evolution of centrifugal instabilities within curved incompressible mixing layers

It is known that certain configurations which possess curvature are prone to a class of instabilities which their 'flat' counterparts will not support. The main thrust of the study of these centrifugal instabilities has concentrated on curved solid boundaries and their effect on the fluid motion. In this article attention is shifted towards a fluid-fluid interface observed within a curved incompressible mixing layer. Experimental evidence is available to support the conjecture that this situation may be subjected to centrifugal instabilities. The evolution of modes with wavelengths comparable with the layer's thickness is considered and the high Taylor/Gortler number regime is also discussed which characterizes the ultimate fate of the modes.

Otto, S. R.↗

Self-Contained Automated Methodology for Optimal Flow Control

This paper describes a self-contained, automated methodology for active flow control which couples the time-dependent Navier-Stokes system with an adjoint Navier-Stokes system and optimality conditions from which optimal states, i.e., unsteady flow fields and controls (e.g., actuators), may be determined. The problem of boundary layer instability suppression through wave cancellation is used as the initial validation case to test the methodology. Here, the objective of control is to match the stress vector along a portion of the boundary to a given vector; instability suppression is achieved by choosing the given vector to be that of a steady base flow. Control is effected through the injection or suction of fluid through a single orifice on the boundary. The results demonstrate that instability suppression can be achieved without any a priori knowledge of the disturbance, which is significant because other control techniques have required some knowledge of the flow unsteadiness such as frequencies, instability type, etc. The present methodology has been extended to three dimensions and may potentially be applied to separation control, re-laminarization, and turbulence control applications using one to many sensors and actuators.

Joslin, Ronald D.↗

Beam-plasma interactions in a positive ion-negative ion plasma

An electron-free plasma consisting of negative ions /SF6(-)/ and positive ions /Ar(+)/, and negligible neutral-ion collision frequencies has been created in the laboratory. This plasma has a mass ratio of approximately 3.5-similar to many computer particle-in-cell simulated systems. A fluid description of this positive and negative ion confinement (PANIC) plasma is given and compared to experimental measurements of a beam-plasma instability for both beam species and a wide range of beam energies. The fluid dispersion relation and most growing modes are predicted to be insensitive to many parameters of the PANIC beam-plasma system, and found to the consistent with the data.

Intrator, T.↗

Advanced Booster Liquid Engine Combustion Stability

Combustion instability is a phenomenon in liquid rocket engines caused by complex coupling between the time-varying combustion processes and the fluid dynamics in the combustor. Consequences of the large pressure oscillations associated with combustion instability often cause significant hardware damage and can be catastrophic. The current combustion stability assessment tools are limited by the level of empiricism in many inputs and embedded models. This limited predictive capability creates significant uncertainty in stability assessments. This large uncertainty then increases hardware development costs due to heavy reliance on expensive and time-consuming testing.

Tucker, Kevin↗

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗

Discretization limits of lattice‐Boltzmann methods for studying immiscible two‐phase flow in porous media

Summary Digital images of porous media often include features approaching the image resolution length scale. The behavior of numerical methods at low resolution is therefore important even for well‐resolved systems. We study the behavior of the Shan‐Chen (SC) and Rothman‐Keller (RK) multicomponent lattice‐Boltzmann models in situations where the fluid‐fluid interfacial radius of curvature and/or the feature size of the medium approaches the discrete unit size of the computational grid. Various simple, small‐scale test geometries are considered, and a drainage test is also performed in a Bentheimer sandstone sample. We find that both RK and SC models show very high ultimate limits: in ideal conditions the models can simulate static fluid configuration with acceptable accuracy in tubes as small as three lattice units across for RK model (six lattice units for SC model) and with an interfacial radius of curvature of two lattice units for RK and SC models. However, the stability of the models is affected when operating in these extreme discrete limits: in certain circumstances the models exhibit behaviors ranging from loss of accuracy to numerical instability. We discuss the circumstances where these behaviors occur and the ramifications for larger‐scale fluid displacement simulations in porous media, along with strategies to mitigate the most severe effects. Overall we find that the RK model, with modern enhancements, exhibits fewer instabilities and is more suitable for systems of low fluid‐fluid miscibility. The shortcomings of the SC model seem to arise predominantly from the high, strongly pressure‐dependent miscibility of the two fluid components.

Li, Zhe↗

Streaming Instability with Multiple Dust Species – I. Favourable Conditions for the Linear Growth

A recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth time-scale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles tomutually validate and study the unstablemodes of the instability in more detail.We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution [Ts, min, Ts, max] and the total solid-to-gas mass ratio ε and introduce a third dimension with the slope q of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the ε–Ts, max space, while this boundary is not appreciably sensitive to q or Ts, min.With a wide range of dust sizes present in the disc (e.g. Ts, min  10−3), the growth rate in the slow-growth regime decreases as more dust species are considered.With a narrow range of dust sizes (e.g. Ts, max/Ts, min = 5), on the other hand, the growth rate in most of the ε–Ts, max space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by ε  1 or Ts, max  1, which may represent the favourable conditions for planetesimal formation.

Zhaohuan Zhu↗

Electrohydrodynamically coupled Kelvin-Helmholtz instability

An inviscid model of the surface dynamics at the interface between two parallel-flowing insulating dielectric fluids, subject to a tangential electric field with a weak gradient, is developed. The two flows are oppositely directed, resulting in an electrohydrodynamically coupled Kelvin-Hemholtz instability. A simple stability criterion is used successfully to correlate experimental data obtained with a test configuration.

Perry, M. P.↗

Rotating disk transition due to isolated roughness with intense acoustic irradiation

The low transition Reynolds numbers associated with boundary-layers which exhibit crossflow instability have made them an important area of study in fluid mechanics. An excellent test bed for the study of this flow is a flat disk rotating in still fluid. An experimental investigation of stability and transition of a rotating disk boundary-layer has been conducted. The effects of isolated roughness elements and acoustic forcing, both alone and in combination, were considered. The flow showed no receptivity to intense monochromatic acoustic forcing except in the presence of sizable roughness. A case was investigated where the growth of the primary instability, seeded by an isolated roughness element, was positively altered in the presence of sound, leading to a downstream movement of transition. In agreement with the work of several investigators, the results of this study suggest that receptivity to long wavelength phenomena, such as traveling sound waves, is contingent upon the presence of strong localized pressure gradients (such as those introduced by roughness) to reduce the length scales of the disturbance to those which are consistent with unstable wavelengths in the flow field.

Waitz, Ian A.↗

Convection effects at solid-liquid interfaces: Influence of gravity

The stability of the flow between two vertical, infinite, rigid, coaxial cylinders held at different temperatures is analyzed by linear stability theory. For a Prandtl number of 22.8 and a radius ratio of 0.02, the flow is unstable to an axisymmetric perturbation at a critical Grashof number of 2150; the wave speed of the instability is comparable to the maximum velocity of the unperturbed flow. When the rigid outer cylinder-fluid interface is replaced by a crystal-melt interface which can change shape, two new modes of instability occur at lower Grashof numbers.

Glicksman, M. E.↗

Magnetic Control of Concentration Gradient in Microgravity

A report describes a technique for rapidly establishing a fluid-concentration gradient that can serve as an initial condition for an experiment on solutal instabilities associated with crystal growth in microgravity. The technique involves exploitation of the slight attractive or repulsive forces exerted on most fluids by a magnetic-field gradient. Although small, these forces can dominate in microgravity and therefore can be used to hold fluids in position in preparation for an experiment. The magnetic field is applied to a test cell, while a fluid mixture containing a concentration gradient is prepared by introducing an undiluted solution into a diluting solution in a mixing chamber. The test cell is then filled with the fluid mixture. Given the magnetic susceptibilities of the undiluted and diluting solutions, the magnetic-field gradient must be large enough that the magnetic force exceeds both (1) forces associated with the flow of the fluid mixture during filling of the test cell and (2) forces imposed by any residual gravitation and fluctuations thereof. Once the test cell has been filled with the fluid mixture, the magnetic field is switched off so that the experiment can proceed, starting from the proper initial conditions.

Leslie, Fred↗

Deleterious Thermal Effects due to Randomized Flow Paths in Pebble Bed, and Particle Bed Style Reactors

Reactor fuel rod surface area that is perpendicular to coolant flow direction (+S) i.e. perpendicular to the P creates areas of coolant stagnation leading to increased coolant temperatures resulting in localized changes in fluid properties. Changes in coolant fluid properties caused by minor increases in temperature lead to localized reductions in coolant mass flow rates leading to localized thermal instabilities. Reductions in coolant mass flow rates result in further increases in local temperatures exacerbating changes to coolant fluid properties leading to localized thermal runaway. Unchecked localized thermal runaway leads to localized fuel melting. Reactor designs with randomized flow paths are vulnerable to localized thermal instabilities, localized thermal runaway, and localized fuel melting.

Moran, Robert P.↗

Use of rotation to suppress thermosolutal convection in directionally solidified binary alloys

Effects of rotation on onset of convection during plane-front directional solidification of Pb-Sn and the pseudobinary system mercury cadmium telluride (Hg(1-x)Cd(x)Te), and on dendritic solidification of Pb-Sn have been studied by means of linear stability analysis. Incorporating Coriolis and centrifugal accelerations into the momentum equation of Coriell et al., we find that under realistic processing conditions, a large degree of stabilization can be achieved using modest rotation rates for both Pb-Sn and mercury cadmium telluride. At a growth velocity of 5 micron/sec and nominal liquid-side temperature gradient of 200 K/cm in Pb-Sn, rotation at 500 rpm results in a hundredfold increase in the critical Sn concentration. Large increases in the maximum allowable growth velocity at fixed melt composition are also attainable with modest rotation rates. The effect is amplified under conditions of reduced gravitational acceleration. For Hg(1-x)Cd(x)Te, we have also studied the nonrotating case. The key differences are due to the existence of a composition range for Hg(1-x)Cd(x)Te in which the melt density has a local maximum as a function of temperature. When the melt solidifies by cooling from below, the liquid density may initially increase with distance above the interface, before ultimately decreasing as the melt temperature increases above the value at which the local density maximum occurs. In contrast to the Pb-Sn case where density depends monotonically on temperature and composition, for Hg(1-x)Cd(x)Te there exists a critical value of the growth velocity above which plane-front solidification is unstable for all bulk CdTe mole fractions. Again, rotation leads to significant inhibition of onset. We identify the predicted stabilization with the Taylor-Proudman mechanism by which rotation inhibits thermal convection in a single-component fluid heated from below. In a binary liquid undergoing solidification, rotation inhibits the onset of buoyancy-driven convection, and has no effect on the short-wavelength morphological instability. At large growth velocities, the plane-front interface between liquid and solid becomes unstable with respect to a morphological instability and solidification occurs dendritically, with a mushy zone of dendrites and interdendritic fluid separating the solid from the melt. For the Pb-Sn system, rotation substantially suppresses the onset of convection in the mushy zone and in the overlying liquid, holding open the promise that rotation can suppress freckling and other macrosegregation defects.

Pearlstein, Arne J.↗

Variable-density buoyancy-driven turbulence with asymmetric initial density distribution

The effects of different initial density distributions on the evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at low (0.05) and high (0.75) Atwood numbers are studied by using high-resolution direct numerical simulations. HVDT aims to mimic the acceleration-driven Rayleigh–Taylor and shock-driven Richtmyer–Meshkov instabilities and reveals new physics that arise from variable-density effects on the turbulent mixing. Here, the initial amounts of pure light and pure heavy flows are altered primarily to mimic the variable-density turbulence at the different locations of the Rayleigh–Taylor and Richtmyer–Meshkov instabilities’ mixing layers where the amounts of the mixing fluids are not equal. It is found that for the low Atwood number cases, the asymmetric initial density distribution has limited effects on both global and local flow evolution for HVDT. However, at high Atwood number, both global flow evolution and the local flow structures are strongly affected by the initial composition ratio. The flow composed of more light fluid reaches higher turbulent levels and the local statistics reach their fully-developed behavior earlier in the time evolution. In contrast, during the late time decay, where most of the flow is well-mixed, all parameters become independent of the initial composition ratio for both low and high Atwood number cases.

42 ENGINEERING↗