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.
A theoretical investigation of the influence of the aerodynamic torque on a Skylab-type satellite is presented. Only the case of free flow of uniform velocity and completely diffuse reflection is considered. It is found by linear stability analysis that the satellite has no stable equilibrium.
A nonlinear proton distribution function that is an exact stationary solution of the nonlinear Vlasov equation and Maxwell's equations and which supports a single nonlinear transverse Alfven (ion cyclotron) wave that is circularly polarized and nondispersive is proposed for most of the observations during high-speed solar wind streams. This nonlinear distribution removes the strong Alfven wave instability, inconsistent with the persistence of the observed proton distribution functions in high-speed streams, found by the linear stability analysis. Model temperature anisotropies and drift velocities of the two spatially inhomogeneous bi-Maxwellian components are consistent with typical proton velocity distributions measured in high-speed streams at 1 AU. Two derived relations for each of the wave number and the phase velocity of the wave are obeyed within experimental uncertainties by two typical proton measurements. Our model also predicts that the alpha particle bulk flow velocity exceeds the proton particle bulk flow velocity, as is observed.
A linear stability analysis of a baroclinic zonal current contained between two parallel rigid boundaries is presented. Curvature is included by performing the analysis on a beta b-plane and viscosity by allowing for the effects of Ekman layers on the rigid boundaries. A two-layer model is used. This calculation was carried out to assist in the design of a spherical model of the general circulation of the earth's atmosphere for Spacelab. In the low-gravity environment on an orbiting vehicle, a dominant radial dielectric body force, analogous to planetary gravity, can be achieved over a volume of liquid held between two concentric spheres. The results show the Eady short wavelength cutoff, and long wavelength cutoffs due to Ekman damping and curvature.
This paper is concerned with analyses of nonlinear waves in solid propellant rockets. Most attention is given to an approximate technique which inexpensively provides results which appear to be quite accurate at least up to amplitudes of ten percent. The connection with linear stability analysis is shown. Primarily to study nonlinear stability, or triggering, the method is extended to third order in the amplitude of wave motion; no explicit results have been obtained. Application of the approximate method to the behavior of pulses is described.
Some problems of attitude stability of spinning satellites are treated in a rigorous manner. With certain restrictions, linearized stability analysis correctly predicts the attitude stability of spinning satellites, even in the critical cases of the Liapunov-Poincare stability theory.
The linear stability analysis for the stratified flow between two rotating circular cylinders is formulated. Two approaches for the stability analysis are presented. The first approach results in an algebraic eigenvalue problem, while the second results in an initial value problem for the perturbation function. The advantages and disadvantages of both approaches are discussed and a preferable numerical solution technique is outlined.
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.
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.
The two-dimensional steady and time dependent properties of plane Poiseuille and plane Couette flows are analyzed using iterative techniques and full numerical simulation of the Navier-Stokes equations. It is shown that the finite-amplitude two-dimensional states investigated are strongly unstable to very small three-dimensional perturbations. It is also shown, through full numerical simulation, that this explosive secondary instability can explain the subcritical transitions that occur in real flows. Finally, it is shown that the three-dimensional instability can be analyzed by a linear stability analysis of a two-dimensional flow consisting of the basic parallel flow and a steady (or quasi-steady) finite-amplitude two-dimensional cellular motion.
A quantitative description of flame stabilization in stagnation-point flow is proposed. Asymptotic and stability analyses are made for a flame model where the density of the gas is assumed to be constant and the reaction zone is assumed to be narrow and concentrated over the flame front. It is shown that, if blowing is sufficiently strong, the corrugations disappear and a plane flame results. The phenomena cannot be fully described by means of classical linear stability analysis.
The evolution of the unstable solar atmosphere into the nonlinear phase, in response to various perturbations, is followed. The initial dynamic evolution of the atmosphere follows the predictions of linear stability analysis. In the nonlinear phase, rapid changes are confined to the transition region; these changes are manifested as a propagation of the transition region through the plasma, i.e., chromospheric evaporation or condensation. Global evolution therefore proceeds on the coronal conductive time scale. The rate of propagation of the transition region is determined by the imbalance between the energy supplied by thermal conduction from the corona and radiative cooling within the transition region itself. Flow velocities in the lower corona during evaporation or condensation are, in the cases studied, of order 3 km/s. The observed dynamic evolution is consistent with the existence of relatively long-lived coronal loops whose brightnesses vary on the evaporative time scale.
In a Benard-type convection problem, the stationary flows of an infinite layer of fluid lying between two rigid horizontal walls and heated uniformly from below are determined. As the temperature difference across the layer increases beyond a certain value, other convective motions appear. These motions areoften cellular in character in that their streamlines are confined to certain well-defined cells having, for example, the shape of rolls or hexagons. A selection principle that explains why hexagonal cells seem to be preferred for certain ranges of the parameters is formulated. An operator-theoretical formulation of one generalized Bernard problem is given. The infinite dimensional problem is reduced to one of solving a finite dimensional system of equations, namely, the selection equations. These equations are solved and a linearized stability analysis of the resultant stationary flows is presented.
Completed and ongoing research activities are discussed briefly, including a three-dimensional, linear stability analysis of the baroclinic Hadley cell and a numerical model of the baroclinic flow between two rotating concentric spheres. This model simulates axisymmetric flow in the Atmospheric General Circulation Experiment configuration. A computer code designed to solve the strongly nonlinear stability problem for the Eady basic state is mentioned.
An elecrostatic finite-sized particle simulation model is used to study the early time-scale phenomena asociated with an impulsively injected plasma, expanding radially, normal to a strong ambient magnetic field. Results of the simulation show the early formation of a radial polarization electric field due to ion-electron charge separation which causes electrons to drift in the azimuthal direction. Velocity shear within this motion gives rise to a Kelvin-Helmholtz (dioxotron) instability, creating a radial fluted pattern as the plasma expands. In the nonlinear stage of instability, azimuthal variations in the electric field are observed which cause electrons to drift across the magnetic field, thereby reducing the space charge created by energetic ion expansion. The results of a linear stability analysis reveal that for a decreasing amount of charge separation in the plasma, the number of unstable azimuthal modes at maximum growth rate increases. Also, as the thickness of the electron ring increases, a fewer number of unstable modes develops for a given amount of charge separation. It is concluded that electrons become unstable at the very early time scale, within an ion gyroperiod, and that unstable electrons lead to modification of the ion dynamics.
For the case of a baroclinic flow whose Richardson number, Ri, is of order unity, a three-dimensional linear stability analysis is conducted on the basis of a model for a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The Hadley cell basic state and stability analysis are both based on the Navier-Stokes and energy equations, and perturbations possessing zonal, meridional, and vertical structures are considered. An attempt is made to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl and Ekman numbers, as well as to finite growth rates and a wider range of the zonal wavenumber. In general, it is found that the symmetric modes of maximum growth are not purely symmetric, but have a weak zonal structure.
A three-dimensional, linear stability analysis of a baroclinic flow for Richardson number, Ri, of order unity is presented. The model considered is a thin horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the complete set of governing, nonlinear equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in a closed form. The stability analysis is also based on the complete set of equations; and perturbation possessing zonal, meridional, and vertical structures were considered. Numerical methods were developed for the stability problem which results in a stiff, eighth-order, ordinary differential eigenvalue problem. The previous work on three-dimensional baroclinic instability for small Ri was extended to a more realistic model involving the Prandtl number, sigma, and the Ekman number, E, and to finite growth rates and a wider range of the zonal wavenumber.
Recent experiments have shown that the presence of a vertical buoyancy-driven flow adjacent to an initially cylindrical crystal-melt interface may produce a time-dependent helical deformation of the interface, with a rotation period ranging from several minutes to many hours, depending upon the width of the melt. The temperature distribution is such that the interface is expected to be morphologically stable in the absence of fluid flow. A linear stability analysis reveals that the instability is due to a coupling between a basic hydrodynamic instability in the buoyant flow and the deformable boundary separating the two phases. The crystal-melt interface lowers the critical Grashof number of an analogous rigid-walled system by an order of magnitude for succinonitrile with a Prandtl number P = 22.8; furthermore, the hydrodynamic mode that is actually destabilized by the interface is not the least stable mode in the rigid-walled system for P = 22.8. The results show that the instability may be regarded either as a rather large alteration of a basic hydrodynamic instability by the crystal-melt interface, or as a significant modification of the morphological stability of the interface by the presence of the buoyant flow.
During solidification of a binary alloy at constant velocity vertically upward, thermosolutal convection can occur if the solute rejected at the crystal-melt interface decreases the density of the melt. We assume that the crystal-melt interface remains planar and that the flow field is periodic in the horizontal direction. The time-dependent nonlinear differential equations for fluid flow, concentration, and temperature are solved numerically in two spatial dimensions for small Prandtl numbers and moderately large Schmidt numbers. For slow solidification velocities, the thermal field has an important stabilizing influence: near the onset of instability the flow is confined to the vicinity of the crystal-melt interface. Further, for slow velocities, as the concentration increases, the horizontal wavelength of the flow decreases rapidly - a phenomenon also indicated by linear stability analysis. The lateral inhomogeneity in solute concentration due to convection is obtained from the calculations. For a narrow range of solutal Rayleigh numbers and wavelengths, the flow is periodic in time.