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.
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.
A derivation of an explicit solution to the two point boundary-value problem of exoatmospheric guidance and trajectory optimization is presented. Fixed initial conditions and continuous burn, multistage thrusting are assumed. Any number of end conditions from one to six (throttling is required in the case of six) can be satisfied in an explicit and practically optimal manner. The explicit equations converge for off nominal conditions such as engine failure, abort, target switch, etc. The self starting, predictor/corrector solution involves no Newton-Rhapson iterations, numerical integration, or first guess values, and converges rapidly if physically possible. A form of this algorithm has been chosen for onboard guidance, as well as real time and preflight ground targeting and trajectory shaping for the NASA Space Shuttle Program.
Results of detailed numerical evaluation of turbulent and reacting three-dimensional compressible flows are presented. These solutions were obtained using a finite element algorithm to solve the system of explicitly non-linear partial differential equations governing three-dimensional, confined unidirectional turbulent flow of a compressible, reacting fluid. Assessments are made of solution accuracy and convergence, including detailed comparisons to experimental data.
A mapping transformation is derived for the inverse solution of nonlinear and linear integral equations of the types encountered in remote sounding studies. The method is applied to the solution of specific problems for the determination of the thermal and composition structure of planetary atmospheres from a knowledge of their upwelling radiance.
Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.
An application-independent geometric model within a data base framework should support the use of Boolean operators which allow the user to construct a complex model by appropriately combining a series of simple models. The use of these operators leads to the concept of implicitly and explicitly defined surfaces. With an explicitly defined model, the surface area may be computed by simply summing the surface areas of the bounding surfaces. For an implicitly defined model, the surface area computation must deal with active and inactive regions. Because the surface intersection problem involves four unknowns and its solution is a space curve, the parametric coordinates of each surface must be determined as a function of the arc length. Various subproblems involved in the general intersection problem are discussed, and the mathematical basis for their solution is presented along with a program written in FORTRAN IV for implementation on the IBM 370 TSO system.
Si3N4, AlN, and Al2O3 were used as basic constituents in a study of the pressureless sintering of beta prime-Si3N4 solid solution as a function of temperature. Y2O3-SiO2 additions were used to promote liquid-phase sintering. The sintered specimens were characterized with respect to density, microstructure, strength, oxidation, and thermal shock resistance. Density greater than 98 percent of theoretical was achieved by pressureless sintering at 1750 C. The microstructure consisted essentially of fine-grained beta prime-Si3N4 solid solution as the major phase. Modulus of rupture strengths up to 483 MPa were achieved at moderate temperature (1000 C), but decreased to 228 MPa at 1380 C. This substantial strength loss was attributed to a glassy grain boundary phase formed during cooling from the sintering temperature. The best oxidation resistance was exhibited by a composition containing 3 mol % Y2O3-SiO2 additives. Water quench thermal shock resistance was equivalent to that of reaction sintered silicon nitride but lower than hot-pressed silicon nitride.
The paper considers a thin elastic bidirectional symmetric composite laminate subjected to a uniform axial strain. All interlaminar stresses are determined by a zeroth-order perturbation analysis of the governing equations. The method of solution considers the laminate to be composed of an interior region (removed from the free edge) and a boundary layer region. The interior region is analyzed by dropping terms related to powers of the thickness-to-width ratio in the dimensionless elasticity equations. The boundary layer region is studied by introducing a stretching transformation to the governing equations. A uniform solution is then formed by satisfying Prandtl's matching principle of perturbation theory. Results for a four-layer bidirectional graphite-epoxy laminate are compared with existing finite-difference results. It is shown that the perturbation analysis provides more insight into the stress distribution near the free edge.
Exact analytical solutions are presented for the standard time-independent spherically symmetric convection-diffusion-adiabatic deceleration equation governing the transport of cosmic rays in the interplanetary medium for the case in which particles are produced with spherical symmetry at the sun. It is assumed that the solar-wind speed is constant and radial, and that the spatial diffusion coefficient has a power-law dependence on momentum. The Green's function describing the modulation of a monoenergetic production of particles is presented. The solutions provide a useful basis for the study of time-integrated properties of energetic solar-flare particle spectra.