Engineering Papers⌕ Search

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.

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

At least 667 records · Page 37

UNCERTAINTIES IN ATOMIC DATA AND THEIR PROPAGATION THROUGH SPECTRAL MODELS. I

We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data. We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of O III and Fe II and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe II].

79 ASTRONOMY AND ASTROPHYSICS↗

Compressed Continuous Computation v. 12/20/2016

A library for performing numerical computation with low-rank functions. The (C3) library enables performing continuous linear and multilinear algebra with multidimensional functions. Common tasks include taking "matrix" decompositions of vector- or matrix-valued functions, approximating multidimensional functions in low-rank format, adding or multiplying functions together, integrating multidimensional functions.

Gorodetsky, Alex↗

Development of an unsteady aerodynamic analysis for finite-deflection subsonic cascades

An unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, was developed to predict aerodynamic forces and moments associated with free vibration or flutter phenomena in the fan, compressor, or turbine stages of modern jet engines. Based on the assumption of small amplitude blade motions, the unsteady flow is governed by linear equations with variable coefficients which depend on the underlying steady low. These equations were approximated using difference expressions determined from an implicit least squares development and applicable on arbitrary grids. The resulting linear system of algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., noniterative) solution. The solution procedure was extended to treat blades with rounded or blunt edges at incidence relative to the inlet flow.

Verdon, J. M.↗

On substructuring algorithms and solution techniques for the numerical approximation of partial differential equations

Substructuring methods are in common use in mechanics problems where typically the associated linear systems of algebraic equations are positive definite. Here these methods are extended to problems which lead to nonpositive definite, nonsymmetric matrices. The extension is based on an algorithm which carries out the block Gauss elimination procedure without the need for interchanges even when a pivot matrix is singular. Examples are provided wherein the method is used in connection with finite element solutions of the stationary Stokes equations and the Helmholtz equation, and dual methods for second-order elliptic equations.

Gunzburger, M. D.↗

Iterative solution of large, sparse linear systems on a static data flow architecture - Performance studies

The applicability of static data flow architectures to the iterative solution of sparse linear systems of equations is investigated. An analytic performance model of a static data flow computation is developed. This model includes both spatial parallelism, concurrent execution in multiple PE's, and pipelining, the streaming of data from array memories through the PE's. The performance model is used to analyze a row partitioned iterative algorithm for solving sparse linear systems of algebraic equations. Based on this analysis, design parameters for the static data flow architecture as a function of matrix sparsity and dimension are proposed.

Reed, D. A.↗

A "Kanes's Dynamics" Model for the Active Rack Isolation System

Many microgravity space-science experiments require vibratory acceleration levels unachievable without active isolation. The Boeing Corporation's Active Rack Isolation System (ARIS) employs a novel combination of magnetic actuation and mechanical linkages, to address these isolation requirements on the International Space Station (ISS). ARIS provides isolation at the rack (international Standard Payload Rack, or ISPR) level. Effective model-based vibration isolation requires (1) an appropriate isolation device, (2) an adequate dynamic (i.e., mathematical) model of that isolator, and (3) a suitable, corresponding controller. ARIS provides the ISS response to the first requirement. This paper presents one response to the second, in a state-space framework intended to facilitate an optimal-controls approach to the third. The authors use "Kane's Dynamics" to develop an state-space, analytical (algebraic) set of linearized equations of motion for ARIS.

Hampton, R. David↗

A "Kane's Dynamics" Model for the Active Rack Isolation System

Many microgravity space science experiments require vibratory acceleration levels unachievable without active isolation. The Boeing Corporation's Active Rack Isolation System (ARIS) employs a novel combination of magnetic actuation and mechanical linkages to address these isolation requirements on the International Space Station (ISS). ARIS provides isolation at the rack (International Standard Payload Rack (ISPR)) level. Effective model-based vibration isolation requires: (1) an appropriate isolation device, (2) an adequate dynamic (i.e., mathematical) model of that isolator, and (3) a suitable, corresponding controller. ARIS provides the ISS response to the first requirement. This paper presents one response to the second, in a state space framework intended to facilitate an optimal-controls approach to the third. The authors use "Kane's Dynamics" to develop a state-space, analytical (algebraic) set of linearized equations of motion for ARIS.

Hampton, R. D.↗

On the conformal symmetry of exceptional scalar theories

The DBI and special galileon theories exhibit a conformal symmetry at unphysical values of the spacetime dimension. We find the Lagrangian form of this symmetry. The special conformal transformations are non-linearly realized on the fields, even though conformal symmetry is unbroken. Commuting the conformal transformations with the extended shift symmetries, we find new symmetries, which when taken together with the conformal and shift symmetries close into a larger algebra. For DBI this larger algebra is the conformal algebra of the higher dimensional bulk in the brane embedding view of DBI. For the special galileon it is a real form of the special linear algebra. We also find the Weyl transformations corresponding to the conformal symmetries, as well as the necessary improvement terms to make the theories Weyl invariant, to second order in the coupling in the DBI case and to lowest order in the coupling in the special galileon case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Recursive Linearization Of Dynamical Models Of Manipulators

Report shows how spatial-operator algebra developed by authors used to derive linearized mathematical models of dynamics of robotic manipulators. Aspects of spatial-operator algebra also described in "Robot Control Based On Spatial-Operator Algebra" (NPO-17918), "Unified Formulation Of Dynamics Of Robot Arms" (NPO-18040), and "Spatial-Operator Algebra For Robotic Manipulators" (NPO-17770).

Jain, Abhinandan↗

A linear programming manual

Computer solutions of linear programming problems are outlined. Information covers vector spaces, convex sets, and matrix algebra elements for solving simultaneous linear equations. Dual problems, reduced cost analysis, ranges, and error analysis are illustrated.

Tuey, R. C.↗

AMG2023

The AMG2023 benchmark solves two diffusion problems with a linear solver preconditioned with algebraic multigrid. The code only contains a driver, a Makefile, and a documentation file. It requires an installation of the open source software library hypre that needs to be downloaded elsewhere and is not included here. Its purpose is to benchmark linear solver performance on high performance computers.

Li, Ruipeng↗

SCF and CI calculations of the dipole moment function of ozone

The constant and linear terms in a Taylor series expansion of the dipole moment function of the ground state of ozone are calculated with Cartesian Gaussian basis sets ranging in quality from minimal to double zeta plus polarization. Results are presented at both the self-consistent field and configuration-interaction levels. Although the algebraic signs of the linear dipole moment derivatives are all established to be positive, the absolute magnitudes of these quantities, as well as the infrared intensities calculated from them, vary considerably with the level of theory.

Curtiss, L. A.↗

Frequency analysis via the method of moment functionals

Several variants are presented of a linear-in-parameters least squares formulation for determining the transfer function of a stable linear system at specified frequencies given a finite set of Fourier series coefficients calculated from transient nonstationary input-output data. The basis of the technique is Shinbrot's classical method of moment functionals using complex Fourier based modulating functions to convert a differential equation model on a finite time interval into an algebraic equation which depends linearly on frequency-related parameters.

Pearson, A. E.↗

Frequency analysis via the method of moment functionals

Several variants are presented of a linear-in-parameters least squares formulation for determining the transfer function of a stable linear system at specified frequencies given a finite set of Fourier series coefficients calculated from transient nonstationary input-output data. The basis of the technique is Shinbrot's classical method of moment functionals using complex Fourier based modulating functions to convert a differential equation model on a finite time interval into an algebraic equation which depends linearly on frequency-related parameters.

Pearson, A. E.↗

Pfaffian systems and feedback linearization/obstructions

The concept of Cartan-Vessiot filtered Lie algebra/module is defined in terms of previous feedback linearizations constructed through canonical forms of Pfaffian systems associated with nonlinear systems. Nonlinearization/Cartan-Vessiot algebras are regarded as deformations of linear ones. It is shown how certain problems of interest in nonlinear system theory translate over to algebraic problems involving this sort of algebraic structure.

Hermann, R.↗

Quantum stochastic calculus associated with quadratic quantum noises

We first study a class of fundamental quantum stochastic processes induced by the generators of a six dimensional non-solvable Lie †-algebra consisting of all linear combinations of the generalized Gross Laplacian and its adjoint, annihilation operator, creation operator, conservation, and time, and then we study the quantum stochastic integrals associated with the class of fundamental quantum stochastic processes, and the quantum Itô formula is revisited. The existence and uniqueness of solution of a quantum stochastic differential equation is proved. The unitarity conditions of solutions of quantum stochastic differential equations associated with the fundamental processes are examined. The quantum stochastic calculus extends the Hudson-Parthasarathy quantum stochastic calculus.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗