SEARCH · Engineering Papers

Results for “linear algebra”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Analytic prediction of airplane equilibrium spin characteristics

The nonlinear equations of motion are solved algebraically for conditions for which an airplane is in an equilibrium spin. Constrained minimization techniques are employed in obtaining the solution. Linear characteristics of the airplane about the equilibrium points are also presented and their significance in identifying the stability characteristics of the equilibrium points is discussed. Computer time requirements are small making the method appear potentially applicable in airplane design. Results are obtained for several configurations and are compared with other analytic-numerical methods employed in spin prediction. Correlation with experimental results is discussed for one configuration for which a rather extensive data base was available. A need is indicated for higher Reynolds number data taken under conditions which more accurately simulate a spin.

Adams, W. M., Jr.↗

NASA NTRS1975

Some properties of extremal solutions of zero-sum differential games

Necessary conditions for infinite-time linear-quadratic two-player differential games with perfect state information are examined in the light of Willems' work on the algebraic Riccati equation. Necessary conditions for the nonzero-sum game do not always reduce precisely to those of the zero-sum game, when the weighting parameters of the cost index are taken as those of a zero-sum game. Formulation of zero-sum game problems in extended L2-spaces reveals the existence of nonstabilizing but infimal controls, and explicit frequency-domain conditions for existence of solutions and stabilizing solutions are established through the algebraic Riccati equation (ARE).

Johnson, T. L.↗

NASA NTRS1983

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity is studied in detail via the eigenstrain approach. A complete formulation and a treatment of both types of eigenstrains for equivalence between the inhomogeneity problem and the inclusion problem are given. This approach is shown to be mathematically identical to other approaches such as the direct volume integral formulation. Expanding the eigenstrains and applied strains in the polynomial form in the position vector and satisfying the equivalence conditions at every point, the governing simultaneous algebraic equations for the unknown coefficients in the eigenstrain expansion are derived. The elastodynamic field outside an ellipsoidal inhomogeneity in a linear elastic isotropic medium is given as an example. The angular and frequency dependence of the induced displacement field, as well as the differential and total cross sections are formally given in series expansion form for the case of uniformly distributed eigenstrains.

Fu, L. S.↗

NASA NTRS1984

Nonlinear evolution of magnetopause tearing modes

Since the magnetosheath plasma is highly turbulent, reconnection at the dayside magnetopause is likely to be temporally unsteady. The tearing mode can be viewed as a model for the unsteady development of a reconnecting magnetic topology. Magnetopause tearing occurs in the guide-field limit and has a wave packet spatial structure in the east-west direction. This paper solves for the nonlinear evolution of a single wavelength guide-field tearing mode including the effects of finite transit time on the Landau resonant electrons. Short wavelength modes evolve algebraically in time with perturbation amplitudes proportional to t-squared. Long wavelength modes are fully nonlinear, and the amplitude grows linearly in time.

Coroniti, F. V.↗

NASA NTRS1988

Recursive mass matrix factorization and inversion: An operator approach to open- and closed-chain multibody dynamics

This report advances a linear operator approach for analyzing the dynamics of systems of joint-connected rigid bodies.It is established that the mass matrix M for such a system can be factored as M=(I+H phi L)D(I+H phi L) sup T. This yields an immediate inversion M sup -1=(I-H psi L) sup T D sup -1 (I-H psi L), where H and phi are given by known link geometric parameters, and L, psi and D are obtained recursively by a spatial discrete-step Kalman filter and by the corresponding Riccati equation associated with this filter. The factors (I+H phi L) and (I-H psi L) are lower triangular matrices which are inverses of each other, and D is a diagonal matrix. This factorization and inversion of the mass matrix leads to recursive algortihms for forward dynamics based on spatially recursive filtering and smoothing. The primary motivation for advancing the operator approach is to provide a better means to formulate, analyze and understand spatial recursions in multibody dynamics. This is achieved because the linear operator notation allows manipulation of the equations of motion using a very high-level analytical framework (a spatial operator algebra) that is easy to understand and use. Detailed lower-level recursive algorithms can readily be obtained for inspection from the expressions involving spatial operators. The report consists of two main sections. In Part 1, the problem of serial chain manipulators is analyzed and solved. Extensions to a closed-chain system formed by multiple manipulators moving a common task object are contained in Part 2. To retain ease of exposition in the report, only these two types of multibody systems are considered. However, the same methods can be easily applied to arbitrary multibody systems formed by a collection of joint-connected regid bodies.

Rodriguez, G.↗

NASA NTRS1991

Effect of periodic accelerations on interface stability in a multilayered fluid configuration

The increasing number of research opportunities in a microgravity environment will benefit not only fundamental studies in fluid dynamics, but also technological applications such as those involving materials processing. In particular, fluid configurations that involve fluid-fluid interfaces would occur in a variety of experimental investigations. This work investigates the stability of a configuration involving fluid-fluid interfaces in the presence of a time-dependent (periodic) forcing. The fluid configuration is multilayered and infinite in extent. The analysis is linear and inviscid, and the acceleration vector is oriented perpendicular to each interface. A Floquent analysis is employed, and the resulting algebraic eigensystem is truncated. Nondimensional parameters appear in the algebraic system. A numerical study is performed to elucidate the regions of instability and the effect of parameter variation on the fluid configuration stability.

Lyell, M. J.↗

NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗

NASA NTRS1997

Complementary Reliability-Based Decodings of Binary Linear Block Codes

This correspondence presents a hybrid reliability-based decoding algorithm which combines the reprocessing method based on the most reliable basis and a generalized Chase-type algebraic decoder based on the least reliable positions. It is shown that reprocessing with a simple additional algebraic decoding effort achieves significant coding gain. For long codes, the order of reprocessing required to achieve asymptotic optimum error performance is reduced by approximately 1/3. This significantly reduces the computational complexity, especially for long codes. Also, a more efficient criterion for stopping the decoding process is derived based on the knowledge of the algebraic decoding solution.

Fossorier, Marc P. C.↗

DOE OSTI2020

Approximate inverse-based block preconditioners in poroelasticity

We focus on the fully implicit solution of the linear systems arising from a three-field mixed finite element approximation of Biot’s poroleasticity equations. The objective is to develop algebraic block preconditioners for the efficient solution of such systems by Krylov subspace methods. In this work, we investigate the use of approximate inverse-based techniques to decouple the native system of equations and obtain explicit sparse approximations of the Schur complements related to the physics-based partitioning of the unknowns by field type. Here, the proposed methods are tested in various numerical experiments including real-world applications dealing with petroleum and geotechnical engineering.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

On non-autonomous dynamical systems

In usual realistic classical dynamical systems, the Hamiltonian depends explicitly on time. In this work, a class of classical systems with time dependent nonlinear Hamiltonians is analyzed. This type of problems allows to find invariants by a family of Veronese maps. The motivation to develop this method results from the observation that the Poisson-Lie algebra of monomials in the coordinates and momenta is clearly defined in terms of its brackets and leads naturally to an infinite linear set of differential equations, under certain circumstances. To perform explicit analytic and numerical calculations, two examples are presented to estimate the trajectories, the first given by a nonlinear problem and the second by a quadratic Hamiltonian with three time dependent parameters. In the nonlinear problem, the Veronese approach using jets is shown to be equivalent to a direct procedure using elliptic functions identities, and linear invariants are constructed. For the second example, linear and quadratic invariants as well as stability conditions are given. Explicit solutions are also obtained for stepwise constant forces. For the quadratic Hamiltonian, an appropriated set of coordinates relates the geometric setting to that of the three dimensional manifold of central conic sections. It is shown further that the quantum mechanical problem of scattering in a superlattice leads to mathematically equivalent equations for the wave function, if the classical time is replaced by the space coordinate along a superlattice. The mathematical method used to compute the trajectories for stepwise constant parameters can be applied to both problems. It is the standard method in quantum scattering calculations, as known for locally periodic systems including a space dependent effective mass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Cauchy-type integral method for solving the linearized one-dimensional Vlasov-Poisson equation

Here, we present a method for solving the linearized Vlasov-Poisson equation, based on analyticity properties of the equilibrium and initial condition through Cauchy-type integrals, that produces algebraic expressions for the distribution and field, i.e., the solution is expressed without integrals. Standard extant approaches involve deformations of the Bromwich contour that give erroneous results for certain physically reasonable configurations or eigenfunction expansions that are misleading as to the temporal structure of the solution. Our method is more transparent, lacks these defects, and predicts previously unrecognized behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Xyce(™) Parallel Electronic Simulator

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗

DOE OSTI2022

Xyce(™) Parallel Electronic Simulator v.7.5

Source record↗

DOE OSTI2013

Xyce(™) Parallel Electronic Simulator v.7.6

Keiter, Eric↗

NASA NTRS1975

Controllability in nonlinear systems

An explicit expression for the reachable set is obtained for a class of nonlinear systems. This class is described by a chain condition on the Lie algebra of vector fields associated with each nonlinear system. These ideas are used to obtain a generalization of a controllability result for linear systems in the case where multiplicative controls are present.

Hirschorn, R. M.↗

NASA NTRS1996

Generalized Bezout's Theorem and its applications in coding theory

This paper presents a generalized Bezout theorem which can be used to determine a tighter lower bound of the number of distinct points of intersection of two or more curves for a large class of plane curves. A new approach to determine a lower bound on the minimum distance (and also the generalized Hamming weights) for algebraic-geometric codes defined from a class of plane curves is introduced, based on the generalized Bezout theorem. Examples of more efficient linear codes are constructed using the generalized Bezout theorem and the new approach. For d = 4, the linear codes constructed by the new construction are better than or equal to the known linear codes. For d greater than 5, these new codes are better than the known codes. The Klein code over GF(2(sup 3)) is also constructed.

Berg, Gene A.↗

NASA NTRSDate not supplied

Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D

The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.

High Lift↗

NASA NTRS2001

Glovebox Integrated Microgravity Isolation Technology (g-LIMIT): A Linearized State-Space Model

Vibration acceleration levels on large space platforms exceed the requirements of many space experiments. The Glovebox Integrated Microgravity Isolation Technology (g-LIMIT) is being built by the NASA Marshall Space Flight Center to attenuate these disturbances to acceptable levels. G-LIMIT uses Lorentz (voice-coil) magnetic actuators to levitate and isolate payloads at the individual experiment/sub-experiment (versus rack) level. Payload acceleration, relative position, and relative orientation measurements are fed to a state-space controller. The controller, in turn, determines the actuator Currents needed for effective experiment isolation. This paper presents the development of an algebraic, state-space model of g-LIMIT, in a form suitable for optimal controller design. The equations are first derived using Newton's Second Law directly, then simplified to a linear form for the purpose of controller design.

Hampton, R. David↗