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.
Iterative methods for the approximate solution of the nonlinear state estimation problem are investigated in which the solution is retained in the form of a finite series of Chebyshev polynomials. Algorithms are presented which allow the state to be estimated from observational data in either the batch or the sequential form. The advantages of these techniques are discussed.
The conical flow solution for axisymmetric supersonic flow past cones has been found to be virtually independent of the ratio of specific heats when normalized in a certain way. A simple rational approximation to this flow is derived. The important singularities and the limiting behavior of the solution are also discussed.
The exact solution of the Lifshitz equations governing the cosmological evolution of an initial fluctuation is presented. Lifshitz results valid for squares of the sound velocity equal to zero and 1/3 are extended in closed form to any equation of state where the pressure equals the total energy density times the square of the sound velocity. The solutions embody all the results found previously for special cases of the square of the sound velocity. It is found that the growth of any initial fluctuation is only an exponential function of time with an exponent of not more than 4/3 and is insufficient to produce galaxies unless the initial fluctuation is very large. A possible way to produce very large initial fluctuations by modifying the equation of state by including gravitational interactions is also examined. It is found that a phase transition can occur at baryonic density of 1 nucleon per cubic Planck length or equivalently, at a time of about 10 to the -43rd power sec. At those early times, the masses allowed by causality requirements are too small to be of interest in galaxy formation.
This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.
A numerical investigation of vortex breakdown has been undertaken in an attempt to understand its properties, and the mechanisms responsible for it. Solutions of the full steady axisymmetric Navier-Stokes equations for breakdown in an unconfined viscous vortex have been obtained for core Reynolds numbers up to 200, for a two-parameter family of assumed upstream velocity distributions. Diffusion and convection of vorticity away from the vortex core, and the strong coupling between the circumferential and axial velocity fields in highly-swirling flows, are shown to lead to stagnation and reversal of the axial flow near the axis. The various theories of vortex breakdown are considered in light of the present numerical solutions.
The ionization structure of hydrogen and helium in dusty nebulae is discussed. General equations for pure hydrogen nebulae have been modified to include helium. An approximate analytical solution is presented for Stroemgren radii of hydrogen and helium in the absence of dust. It is shown that these results, with simple modifications, are also applicable to dusty nebulae where the effective absorption cross section of dust grains varies slowly with frequency in the range from 1000 to 200 A. No analytical solutions are possible if this cross section varies rapidly with frequency. In this case, however, simple coupled differential equations are derived which can easily be solved numerically. Approximate analytical expressions are given for evaluating the variation of the fraction of ionizing radiation absorbed by dust and the ratio of the volume emission measures of He II to H II regions with the spectrum of the ionizing source, helium abundance, and absorption properties of dust. The effects of dust on the He III zone are discussed. The present results are restricted to spherically symmetric nebulae, but nonuniform gas and dust distributions and clumpiness can be taken into account by the general results.
A simple solution to the problem of the acoustical coupling between a rectangular structure, its air content, and an external noise source is presented. This solution is a mathematical expression for the normalized acoustic pressure inside the structure. Numerical results for the sound-pressure response for a specified set of parameters are also presented.
A numerical method is presented for solving the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to three-dimensional supersonic laminar jet flow issuing parallel with a supersonic free stream. A coordinate transformation is introduced which maps the boundaries at infinity into a finite computational domain in order to eliminate difficulties associated with the imposition of free-stream boundary conditions. Results are presented for an approximate circular jet, a square jet, varying aspect ratio rectangular jets, and interacting square jets. The solution behavior varies from axisymmetric to nearly two-dimensional in character. For cases where comparisons of the present results with those obtained from shear layer calculations could be made, agreement was good.
Methods are given for computing the covariance matrix for the solution vector of an equality-constrained least squares problem. The methods are matched to the solution algorithms given in the book, 'Solving Least Squares Problems.'
A mathematically rigorous derivation is given of first order corrections to multi-impulse approximations to the solutions to space flight optimization problems with bang-bang control. The rocket was subjected to an inverse square gravitational force and to a thrust force with constant magnitude. The mass decreased linearly with time. An optimal impulsive solution was obtained for a problem with given initial and final conditions. The method was then used to obtain first-order corrections to the initial values of the costate variables. Indications are given on how the theory may be extended to higher order corrections. The theory was applied to intercept and rendezvous problems.
An experiment on melting and directional crystallization of an antimony (Sb) doped germanium silicon (GeSi) solid solution was designed for the Apollo-Soyuz Test Project (ASTP) to study the possibility of using zero-g conditions for obtaining solid-solution monocrystals with uniformly distributed components. Crystallization in the zero-g environment did not occur under ideal stationary growth and segregation conditions. Crystallization under zero-g conditions revealed the heterogeneous nature of Si and Sb distribution in the cross sections of crystals. The presence of the radial thermal gradient in the multipurpose furnace could be one of the reasons for such Si and Sb distribution. The structure of space-grown crystals correlates with the nature of heterogeneities of Si and Sb distribution in crystals. The type of surface morphology and the contour observed in space-grown crystals were never observed in ground-based crystals and indicate the absence of wetting of the graphitized walls of the ampoule by the melt during melting and crystallization.
A method is presented for automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multi-connected two-dimensional region containing any number of arbitrarily shaped bodies. No restrictions are placed on the shape of the boundaries, which may even be time-dependent, and the approach is not restricted in principle to two dimensions. With this procedure the numerical solution of a partial differential system may be done on a fixed rectangular field with a square mesh with no interpolation required regardless of the shape of the physical boundaries, regardless of the spacing of the curvilinear coordinate lines in the physical field, and regardless of the movement of the coordinate system in the physical plane. A number of examples of coordinate systems and application thereof to the solution of partial differential equations are given. The FORTRAN computer program and instructions for use are included.
A solution is presented to the radiative transfer of the solar irradiation through a turbid atmosphere, based on the single-scattering approximation, i.e., an assumption that a photon that underwent scattering either leaves the top of the atmosphere or strikes the surface. The solution depends on a special idealization of the scattering phase function of the aerosols. The equations developed are subsequently applied to analyze quantitatively the enhancement of the surface irradiation and the enhancement of the scattered radiant emittance as seen from above the atmosphere, caused by the surface reflectance and atmospheric back scattering. An order of magnitude error analysis is presented.
An experiment performed during the Apollo-Soyuz flight investigated a novel process for growing single crystals from solutions. In zero gravity, reactant solutions were allowed to diffuse slowly toward each other through a region of pure solvent in which they underwent a chemical reaction to form crystals of an insoluble material. The process is analogous to the growth of crystals in gels in normal gravity. The experiment was successful and yielded crystals of about the expected size, number, and quality. The results indicate that with further refinement the process should yield large, high quality crystals of a variety of materials.
Some properties are derived for the solutions of the variational equations of a class of dynamical systems. It is shown that under rather general conditions, the matrix of the linearized Lagrangian equations of motion have an important property for which the word 'skew-symplectic' has been introduced. It is also shown that the fundamental matrix of solutions is 'symplectic', the word symplectic being used here in a more general sense than in the classical literature. Two consequences of the symplectic property are that the fundamental matrix is easily invertible and that the eigenvalues appear in reciprocal pairs. The effect of coordinate transformations is also analyzed; in particular, the change from Lagrangian to canonical systems.
In the Boussinesq approximation, single-mode equations describing thermal convection are constructed by expanding the fluctuating velocity and temperature fields in a complete set of functions (or planforms) of the horizontal coordinates and retaining just one term. Numerical solutions of the single-mode equations are investigated, chief consideration being given to hexagonal planforms. Extensive surveys of steady solutions are presented for various Rayleigh numbers, Prandtl numbers, and horizontal wavenumbers. The dependences on Rayleigh number and Prandtl number at very large Rayleigh number are in satisfactory agreement with the results of asymptotic expansions.
As an alternative to pure finite-difference techniques, the Method of Integral Relations (MIR) is formulated for three-dimensional boundary-layer flows with separation, and applied to two problems. The first concerns an incompressible laminar flow over a plate with an attached cylinder, where the potential solution of the flow round the circular cylinder is used as the boundary condition. In the second application the solution to the three-dimensional boundary-layer flow over blunt bodies is discussed and results for an ellipsoid of revolution at 30 deg incidence are presented. In both cases the results are found to be in satisfactory agreement with those obtained by finite-difference methods.
The problem of calculating the steady-state free-electron energy distribution in a hydrogen gas is considered in order to study departures of that distribution from a Maxwellian at sufficiently low degrees of ionization. A model kinetic equation is formulated and solved analytically for the one-particle electron distribution function in a steady-state partially ionized hydrogen gas, and it is shown that the formal solution can be accurately approximated by using the WKB method. The solutions obtained indicate that the high-energy tail of the distribution is susceptible to distortion by imbalanced inelastic collisions for ionization fractions not exceeding about 0.1 and that such departures from a Maxwellian can lead to significant changes in the collisional excitation and ionization rates of ground-state hydrogen atoms. Expressions for the electron-hydrogen collision rates are derived which explicitly display their dependence on the hydrogen departure coefficients. The results are applied in order to compare self-consistent predictions with those based on the a priori assumption of a Maxwellian distribution for models of the thermal ionization equilibrium of hydrogen in the optically thin limit, spectral-line formation by a gas consisting of two-level atoms, and radiative transfer in finite slabs by a gas of four-level hydrogen atoms.