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.
An iterative method for solving nonsymmetric indefinite linear systems is proposed. The method involves the successive use of a modified version of the conjugate residual method. A numerical example is given to illustrate the method.
The purpose of this report is to establish linear, decoupled models of rigid body motion for the fixed wing configuration of the Rotor Systems Research Aircraft (RSRA). Longitudinal and lateral control surface fixed linear models were created from aircraft time histories using current system identification techniques. Models were obtained from computer simulation at 160 KCAS and 200 KCAS, and from flight data at 160 KCAS. Comparisons were performed to examine modeling accuracy, variation of dynamics with airspeed and correlation of simulation and flight data results. The results showed that the longitudinal and lateral linear models accurately predicted RSRA dynamics. The flight data results showed that no significant handling qualities problems were present in the RSRA fixed wing aircraft at the flight speed tested.
A brief description of the Model of the World Economy implemented at the Institute for Economic Analysis is presented, together with our experience in converting the software to vector code. For each time period, the model is reduced to a linear system of over 2000 variables. The matrix of coefficients has a bordered block diagonal structure, and we show how some of the matrix operations can be carried out on all diagonal blocks at once.
During the process of solving a mathematical model numerically, there is often a need to operate on a vector v by an operator which can be expressed as f(A) while A is NxN matrix (ex: exp(A), sin(A), A sup -1). Except for very simple matrices, it is impractical to construct the matrix f(A) explicitly. Usually an approximation to it is used. In the present research, an algorithm is developed which uses a polynomial approximation to f(A). It is reduced to a problem of approximating f(z) by a polynomial in z while z belongs to the domain D in the complex plane which includes all the eigenvalues of A. This problem of approximation is approached by interpolating the function f(z) in a certain set of points which is known to have some maximal properties. The approximation thus achieved is almost best. Implementing the algorithm to some practical problem is described. Since a solution to a linear system Ax = b is x= A sup -1 b, an iterative solution to it can be regarded as a polynomial approximation to f(A) = A sup -1. Implementing the algorithm in this case is also described.
Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.
Deriving algorithms for computations involving sparse matrices
Explore the source record for details and available documents.
Explore the source record for details and available documents.
FORTRAN coded out-of-core equation solvers that solve using direct methods symmetric banded systems of simultaneous algebraic equations. Banded, frontal and column (skyline) solvers were studied as well as solvers that can partition the working area and thus could fit into any available core. Comparison timings are presented for several typical two dimensional and three dimensional continuum type grids of elements with and without midside nodes. Extensive conclusions are also given.
Explore the source record for details and available documents.
Work aimed at modeling the charge voltage and current characteristics of nickel-cadmium cells subject to taper charge is presented. Work reported at previous NASA Battery Workshops has shown that the voltage of cells subject to constant current charge and discharge can be modeled very accurately with the equation: voltage = A + (B/(C-X)) + De to the -Ex where A, B, D, and E are fit parameters and x is amp-hr of charge removed during discharge or returned during charge. In a constant current regime, x is also equivalent to time on charge or discharge.
Explore the source record for details and available documents.
Geometrical methods for formation flying design based on the analytical solution to Hill's equations have been previously developed and used to specify desired relative motions in near circular orbits. By generating relationships between the vehicles that are intuitive, these approaches offer valuable insight into the relative motion and allow for the rapid design of satellite configurations to achieve mission specific requirements, such as vehicle separation at perigee or apogee, minimum separation, or a specific geometrical shape. Furthermore, the results obtained using geometrical approaches can be used to better constrain numerical optimization methods; allowing those methods to converge to optimal satellite configurations faster. This paper presents a set of geometrical relationships for formations in eccentric orbits, where Hill.s equations are not valid, and shows how these relationships can be used to investigate formation designs and how they evolve with time.
Linear functional differential equations trajectory optimization, proving maximal principle
Localization of linear integral equation eigenvalues with applications to linear ordinary differential equations
The role of electrostatic conic instability in the diffusion of spilled particles and the acceleration of the auroral plasma of the Earth are numerically investigated by means of quasi-linear equations. Linear increments to the conic instability under magnetospheric conditions are introduced. The quasi-linear equation is averaged along tubes of force taking into account the drift across the magnetic field lines. It is shown that the existence of a conic instability leads to small spills and isotropic distribution of particles, but is also accompanied by significant acceleration of electrons.
The equivalence of a considered integral equation form with an infinite system of linear equations is proved, and the localization of the eigenvalues of the infinite system is expressed. Error estimates are derived, and the problems of finding upper bounds and lower bounds for the eigenvalues are solved simultaneously.
Linearized and weakly nonlinear Boltzmann equation boundary value problems for gas between parallel plates, noting solution existence and uniqueness