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 modest extension of conventional WKB methods of solving second-order differential equations asymptotically is developed and applied to the problem of baroclinic instability. Accurate results are obtained for the Charney problem using simple expressions. More important, the present expressions can be used immediately for a wide variety of basic velocity profiles.
Accurate approximations are presented for the self-similar structures of nonradiating blast waves with adiabatic ions, isothermal electrons, and equation ion and electron temperatures at the shock. The cases considered evolve in cavities with power law ambient densities (including the uniform density case) and have negligible external pressure. The results provide the early time asymptote for systems with shock heating of electrons and strong thermal conduction. In addition, they provide analytical results against which two fluid numerical hydrodynamic codes can be checked.
The selection of the proper coordinate system in solving any fluid flow or heat transfer problem is a very important consideration. A new technique of moving mesh points in physical space is introduced so as to reduce the error in a computed asymptotic solution relative to that obtained using a fixed mesh. The technique has been used to solve the simple viscous Burgers' equation in one and two dimensions. Substantial error reductions are demonstrated. The treatment of boundary points and the effect of using different error criteria in generating grids are discussed.
The asymptotic response of the Charney baroclinic instability problem to a localized perturbation is determined using the formalism of Briggs (1964) and exploiting a recently obtained highly accurate WKB approximate dispersion relation (Lindzen and Rosenthal, 1981). Comparison is made with previous results for two-level and Eady models. Small scales and rapid growth characteristic of the initial stages of cyclogenesis are found and the linear dispersion relation, which can be obtained from observation of zonal wind and stability, emerges as a forecast tool for prediction of geographically local cyclogenesis.
Computer programs were developed to implement the computational scheme arising from Van Holten's asymptotic method for calculating airloads on a helicopter rotor blade in forward flight, and a similar technique which is based on a discretized version of the method. The basic outlines of the two programs are presented, followed by separate descriptions of the input requirements and output format. Two examples illustrating job entry with appropriate input data and corresponding output are included. Appendices contain a sample table of lift coefficient data for the NACA 0012 air foil and listings of the two programs.
Analytical approximations for the acoustic modes in a duct carrying a uniform core flow with a thin shear layer at the walls are developed using the Method of Matched Asymptotic Expansions. Both two-dimensional and cylindrical duct propagation are considered. Numerical results for eigenvalues calculated using the theory are presented for the two dimensional problem and compared with results from earlier analyses. It is found that the new approximations yield a significant increase in accuracy.
Asymptotic flows inside curved ducts of rectangular as well as polar cross sections have been analyzed using the Navier-Stokes equations in terms of the axial velocity and vorticity and the cross-flow stream function. Numerical solutions of the three second-order coupled elliptic partial differential equations governing this flow have been obtained using the multi-grid-SI scheme, to determine fine-grid results for configurations with strong secondary flow.
Methods that can be used in the numerical implementation of the uniform swallowtail approximation are described. An explicit expression for that approximation is presented to the lowest order, showing that there are three problems which must be overcome in practice before the approximation can be applied to any given problem. It is shown that a recently developed quadrature method can be used for the accurate numerical evaluation of the swallowtail canonical integral and its partial derivatives. Isometric plots of these are presented to illustrate some of their properties. The problem of obtaining the arguments of the swallowtail integral from an analytical function of its argument is considered, describing two methods of solving this problem. The asymptotic evaluation of the butterfly canonical integral is addressed.
During hot firing ground testing of the Space Shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.
The effort to develop classical methods to compute wall interference at transonic speeds is outlined. The two-dimensional theory and three-dimensional development are discussed. Also, some numerical application of the two-dimensional work are indicated. The basic advantages of the asymptotic theory are noted.
During hot firing ground testing of the Space shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.
The optimal control of structures which consist of composite materials or of perforated materials is discussed. Asymptotic formula, derived from the so-called homogenization theory, are presented which allow the replacement of very complicated problems by much simpler ones.
A system of parabolic nonlinear equations that describe the diffusion of a fully collisional plasma across a strong magnetic field is discussed. It is demonstrated that the solution to this system tends to a time asymptotic state which is of space-time separable form, theta(t)f(x). Furthermore, f(x) is independent of the initial conditions and theta(t) depends slightly on the initial conditions. The rate of decay of the temporal part is governed by a nonlinear eigenvalue problem. Since the equations are considered in a bounded domain we are able to analyze the effect of boundary conditions on the evolution of the system. Additional effects on radiation, heating, and particle injection can also be accounted for. Essential differences between the behavior of a fully-coupled system and a scalar equation are observed.
An investigation of high Reynolds number stationary instabilities in the boundary layer on a rotating disc is given. It is shown that in addition to the inviscid mode at high Reynolds numbers, there is a stationary short wavelength mode. This mode has its structure fixed by a balance between viscous and Coriolis forces and cannot be described by an inviscid theory. The asymptotic structure of the wavenumber and orientation of this mode is obtained. A similar analysis is given for the inviscid mode, the expansion procedure used is capable of taking nonparallel effects into account in a self consistent manner. The results are compared to numerical calculations and experimental observations.
Formal asymptotic expansions with differently scaled variables are matched to obtain a uniform approximation to the similarity solution of the shock-wedge diffraction problem.
The delegation of the labor of calculating higher-order terms in singular perturbation (SP) expansions to a computer by the use of MACSYMA is considered. The method of matched asymptotic expansions is studied in detail for two model SP problems: a model resembling the boundary layer equation with a small parameter multiplying the highest derivatives; and a turning-point problem. It is shown that MACSYMA has successfully performed the higher-order matching in both problems.
An asymptotic, large wave number approximation for the Helmholtz equation is derived. The theory is an extension of the geometric acoustic theory, and provides corrections to that theory in the form of multiplicative functions which satisfy parabolic equations. A simple example is used both to illustrate failure of the geometric theory for large propagation distances, and to show the improvement obtained by use of the new theory.
The paper considers the problem of how the momentum of ions created by electron impact ionization of neutrals moving at a speed v(0) perpendicular to the magnetic field through a background plasma is coupled to this plasma. It has been found that the plasma accelerates, and the relative velocity between neutrals and plasma decreases. If this decrease is rapid and large enough, the critical ionization velocity (CIV) phenomenon may turn off. Equations for the evolution of plasma density, electron and ion thermal energy, and plasma velocity have been derived. It was found that the CIV process reaches an asymptotic quasi-steady state, in which the ionization rate reaches a constant value which depends on the properties of the surrounding medium and the value of v(0).