SEARCH · Engineering Papers
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.
The paper presents a computational analysis of liquid rocket combustion instability. Consideration is given to both a fully nonlinear unsteady calculation as well as a new CFD-based linearized stability analysis. An analytical solution for the linear stability problem in a constant area combustion chamber with uniform mean flow is developed to verify the numerical analyses.
The onset of instabilities in a fluid contained in a rotating hemispherical shell, driven by thermal gradients imposed upon the hemispherical boundaries and by a spherically symmetric radial body force, is numerically studied. Computations are presented for a range of Taylor and thermal Rossby numbers. The analysis indicates the presence of an instability dependent upon the spherically radial gravity alone when the warmest temperatures are at the pole and an additional centrifugal buoyant instability for weak imposed gravity and fast rotation when the temperature decreases poleward.
The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.
An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, $\sigma_\rho/\bar{\rho}\sim1$) and strong shocks (shock Mach numbers $\bar{M}_s\gtrsim10$). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.
High-speed visible imaging of sub-microsecond electric explosion of wires at the low specific energy deposition threshold reveals three distinct modes of wire failure as capacitor charge voltage and energy deposition are increased. For 100 micron diameter gold-plated tungsten wires of 2 cm length, deposited energies of 1.9 eV/atom produces a liquid column that undergoes hydrodynamic breakup into droplets with radii of order of wire diameter on timescales of 200 μs. Instability growth, column breakup, and droplet coalescence follow classic Rayleigh-Plateau predictions for instability of viscous fluid column. Above 3.2 eV/atom of deposited energy, wires are seen to abruptly transition to an expanding mixture of micron scale liquid-droplets and vapor within one frame (less than 3.33 μs), which has been termed ‘phase explosion’ in literature. Between these two limits, at 2.5 eV/atom of deposited energy, wire radius is unchanged for the first 10 μs before the onset of a rapid expansion and disintegration that resembles homogenous nucleation of mechanically unstable bubbles. Here, thermodynamic calculations are presented that separate cases by temperature obtained during heating: below boiling point, near boiling point, and exceeding boiling point.
The footward shift of blood volume and compensatory autonomic reflex responses initiated during the assumption of the upright posture in terrestrial gravity is a common feature in humans. When regular exposure to upright posture is removed, the headward redistribution of blood induces numerous physiological adaptations which compromise this normal response and the development of low blood pressure upon standing (orthostatic hypotension) can ensue. Such microgravity conditions can be produced by prolonged exposure to spaceflight or prolonged bed rest. Since the reduction of blood and plasma volume during exposure to microgravity has been associated with orthostatic instability following spaceflight and bed rest, it has been a reasonable assumption that hypovolemia may be a primary contributing factor to the development of orthostatic hypotension. This view was supported by the observation that the attempt to restore vascular volume in astronauts by drinking saline solutions just prior to re-entry had proven effective in reducing orthostatic instability after spaceflights of short duration. However, as the duration of spaceflight lengthens, orthostatic instability persists despite fluid loading procedures, suggesting that mechanisms other than hypovolemia may contribute to orthostatic hypotension following prolonged exposure to microgravity conditions. Recent evidence has been generated from spaceflight and groundbase experiments that supports the notion that changes in autonomic baroreflexes that control cardiac and vascular responses during orthostatic challenges may be affected by longer periods of microgravity exposure and can contribute to postflight orthostatic instability.
This paper addresses the critical need for the development of intelligent or assisting software tools for the scientist who is working in the initial problem formulation and mathematical model representation stage of research. In particular, examples of that representation in fluid dynamics and instability theory are discussed. The creation of a mathematical model that is ready for application of certain solution strategies requires extensive symbolic manipulation of the original mathematical model. These manipulations can be as simple as term reordering or as complicated as discovery of various symmetry groups embodied in the equations, whereby Backlund-type transformations create new determining equations and integrability conditions or create differential Grobner bases that are then solved in place of the original nonlinear PDEs. Several examples are presented of the kinds of problem formulations and transforms that can be frequently encountered in model representation for fluids problems. The capability of intelligently automating these types of transforms, available prior to actual mathematical solution, is advocated. Physical meaning and assumption-understanding can then be propagated through the mathematical transformations, allowing for explicit strategy development.
Electrostatic heating of solar wind ions from instability in two fluid solar wind model leading to preferential heating of protons over electrons
A general formulation of some effects of swirl on turbulent mixing is given. The basis for the analysis is that momentum transport is enhanced by turbulence resulting from rotational instability of the fluid field. An appropriate form for the turbulent eddy viscosity is obtained by mixing length type arguments. The result takes the form of a corrective factor that is a function of the swirl and acts to increase the eddy viscosity. The factor is based upon the initial mixing conditions implying that the rotational turbulence decays in a manner similar to that of free shear turbulence. Existing experimental data for free jet combustion are adequately matched by using the modifying factor to relate the effects of swirl on eddy viscosity. The model is extended and applied to the supersonic combustion of a ring jet of hydrogen injected into a constant area annular air stream. The computations demonstrate that swirling the flow could: (1) reduce the burning length by one half, (2) result in more uniform burning across the annulus width, and (3) open the possibility of optimization of the combustion characteristics by locating the fuel jet between the inner wall and center of the annulus width.
A fluid mechanical convective instability is examined under the conditions found in the solar corona. Linearized density perturbations are shown to grow as they are carried outward by the solar wind. The non-linear instability may be proposed as a cause of coronal structures in nonmagnetic stars, or in regions of the solar corona with weak magnetic fields.
Fluid waves and instabilities are considered critical to the evolution of protoplanetary nebulae, particularly for their roles in mass, angular momentum, and energy transport. A number have been identified, however, notably absent, is an influential wave commonly found in planetary atmospheres and oceans: the planetary Rossby wave (PRW). Since, in the Earth's atmosphere, the PRW is of primary importance in shaping large-scale meteorological phenomena, it is reasonable to consider whether it might have similar importance in the protoplanetary nebula. The thrust of the research project this summer (1998) was to determine whether a nebular analog to the PRW is viable, a so-called nebular Rossby wave (NRW), and if so, to explore possible ramifications of this wave to the evolution of the nebula. This work was carried out primarily by S. Davis, J. Cuzzi and me, with significant discussions with P. Cassen. We believe we have established a good case for the NRW and as a result believe we have opened up a new and possibly interesting line of research in regard to the nebular development, in particular with regard to zonal jet formation, a potent accretion mechanism, and possible ties to vortex formation. The standard model of the protoplanetary nebula consists of a large disk of gas with about 1% entrained dust gravitationally bound to a large central mass, m(sub c) i.e., the protostar. The planet-forming region of the disk extends to roughly 100 A.U. in radius. Disk thickness, H, is believed to be on the order of 10-100 times less than disk radius. Disk lifetime is on the order of a million years.
Richtmyer-Meshkov (RM) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and eventually become a turbulent flow. It is closely related to Rayleigh-Taylor instability, which is the instability of a planar interface undergoing constant acceleration, such as caused by the suspension of a heavy fluid over a lighter one in the earth's gravitational field. Like the well-known Kelvin-Helmholtz instability, RM instability is a fundamental hydrodynamic instability which exhibits many of the nonlinear complexities that transform simple initial conditions into a complex turbulent flow. Furthermore, the simplicity of RM instability (in that it requires very few defining parameters), and the fact that it can be generated in a closed container, makes it an excellent test bed to study nonlinear stability theory as well as turbulent transport in a heterogeneous system. However, the fact that RM instability involves fluids of unequal densities which experience negligible gravitational force, except during the impulsive acceleration, requires RM instability experiments to be carried out under conditions of microgravity. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The impulsive acceleration is generated by bouncing a rectangular tank containing two different density liquids off a retractable vertical spring. The initial perturbation is produced prior to release by oscillating the tank in the horizontal direction to produce a standing wave. The instability evolves in microgravity as the tank travels up and then down the vertical rails of a drop tower until hitting a shock absorber at the bottom. Planar Laser Induced Fluorescence (PLIF) is employed to visualize the flow. PLIF images are captured by a video camera that travels with the tank. Figure 1 is as sequence of images showing the development of the instability from the initial sinusoidal disturbance far into the nonlinear regime which is characterized by the appearance of mushroom structures resulting from the coalescence of baroclinic vorticity produced by the impulsive acceleration. At later times in this sequence the vortex cores are observed to become unstable showing the beginnings of the transition to turbulence in this flow. The amplitude of the growing disturbance after the impulsive acceleration is measured and found to agree well with theoretical predictions. The effects of Reynolds number (based on circulation) on the development of the vortices and the transition to turbulence are also determined.
The effect of negative aerodynamic drag in an ideal fluid subject to convective instability is considered. It is shown that a cylinder moving in such a fluid is propelled forward in its motion by the convective forces and that the characteristic acceleration time is comparable to the onset time of convective motions in the fluid. It is suggested that convective propulsion plays an important role in the dynamics of flux tubes extending through the surface of the sun. The suppression of the upward heat flow in a Boussinesq convective cell with free upper and lower boundaries by a downdraft is then analyzed. Application to the solar convection zone indicates that downdrafts of 1 to 2 km/s at depths of 1000 to 4000 km beneath the visible surface of the sun are sufficient to reduce the upward heat flux to a small fraction of the ambient value.
The use of couplings for high speed turbocompressors or pumps is essential to transmit power from the driver. Typical couplings are either of the lubricated gear or dry diaphragm type design. Gear couplings have been the standard design for many years and recent advances in power and speed requirements have pushed the standard design criteria to the limit. Recent test stand and field data on continuous lube gear type couplings have forced a closer examination of design tolerances and concepts to avoid operational instabilities. Two types of mechanical instabilities are reviewed in this paper: (1) entrapped fluid, and (2) gear mesh instability resulting in spacer throw-out onset. Test stand results of these types of instabilities and other directly related problems are presented together with criteria for proper coupling design to avoid these conditions. An additional test case discussed shows the importance of proper material selection and processing and what can happen to an otherwise good design.
The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is γ τ A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of γ τ A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.
Abstract The resistance to volumetric deformations displayed by polymer networks is largely due to secondary and tertiary interactions between neighboring polymer chains. These interactions are both entropic and enthalpic in nature but are fundamentally different from the entropic forces that resist shearing in these networks. In this paper, we introduce a new depiction of elastomers as a crosslinked Van der Waals fluid. Starting from first principles, we develop constitutive equations that are implemented in a continuum model as well as a discrete network model. Our models predict that the failure of polymer networks may be driven by an instability in the underlying polymer bulk ‘fluid’ or by the breaking of polymer chains, depending on the loading path taken. The results of this study indicate that material failure in elastomers exposed to a purely triaxial state, such as in a poker chip experiment, may be driven by an entirely different mode of instability than those deformed in pure shear, such as in a uniaxial tension experiment.
We consider the superfluid phase of a specific renormalizable relativistic quantum field theory. We prove that, within the regime of validity of perturbation theory and of the superfluid effective theory, there are consistent and regular vortex solutions where the superfluid’s velocity field as traditionally defined smoothly interpolates between zero and arbitrarily large superluminal values. We show that this solution is free of instabilities and of superluminal excitations. We show that, in contrast, a generic vortex solution for an ordinary fluid does develop an instability if the velocity field becomes superluminal. All this questions the characterization of a superfluid velocity field as the actual velocity of “something”.
A device for controlling fluid flow. The device includes an arc generator coupled to electrodes. The electrodes are placed adjacent a fluid flowpath such that upon being energized by the arc generator, an arc filament plasma adjacent the electrodes is formed. In turn, this plasma forms a localized high temperature, high pressure perturbation in the adjacent fluid flowpath. The perturbations can be arranged to produce vortices, such as streamwise vortices, in the flowing fluid to control mixing and noise in such flows. The electrodes can further be arranged within a conduit configured to contain the flowing fluid such that when energized in a particular frequency and sequence, can excite flow instabilities in the flowing fluid. The placement of the electrodes is such that they are unobtrusive relative to the fluid flowpath being controlled.