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 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Quasi-exactly-solvable confining solutions for spin-1 and spin-0 bosons in (1+1)-dimensions with a scalar linear potential

We point out a misleading treatment in the recent literature regarding confining solutions for a scalar potential in the context of the Duffin–Kemmer–Petiau theory. We further present the proper bound-state solutions in terms of the generalized Laguerre polynomials and show that the eigenvalues and eigenfunctions depend on the solutions of algebraic equations involving the potential parameter and the quantum number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A matrix-free approach for finite-strain hyperelastic problems using geometric multigrid

This paper investigates matrix-free algorithms for problems in quasi-static finite-strain hyperelasticity. Iterative solvers with matrix-free operator evaluation have emerged as an attractive alternative to sparse matrices in the fluid dynamics and wave propagation communities because they significantly reduce the memory traffic, the limiting factor in classical finite element solvers. Specifically, we study different matrix-free realizations of the finite element tangent operator and determine whether generalized methods of incorporating complex constitutive behavior might be feasible. Furthermore, in order to improve the convergence behavior of iterative solvers, we also propose a method by which to construct level tangent operators and employ them to define a geometric multigrid preconditioner. Additionally, the performance of the matrix-free operator and the geometric multigrid preconditioner is compared to the matrix-based implementation with an algebraic multigrid (AMG) preconditioner on a single node for a representative numerical example of a heterogeneous hyperelastic material in two and three dimensions. We find that matrix-free methods for finite-strain solid mechanics are very promising, outperforming linear matrix-based schemes by two to five times, and that it is possible to develop numerically efficient implementations that are independent of the hyperelastic constitutive law.

42 ENGINEERING↗

DOE OSTI2016

Scalable algorithms for three-field mixed finite element coupled poromechanics

We introduce a class of block preconditioners for accelerating the iterative solution of coupled poromechanics equations based on a three-field formulation. The use of a displacement/velocity/pressure mixed finite-element method combined with a first order backward difference formula for the approximation of time derivatives produces a sequence of linear systems with a 3 x 3 unsymmetric and indefinite block matrix. The preconditioners are obtained by approximating the two-level Schur complement with the aid of physically-based arguments that can be also generalized in a purely algebraic approach. A theoretical and experimental analysis is presented that provides evidence of the robustness, efficiency and scalability of the proposed algorithm. In conclusion, the performance is also assessed for a real-world challenging consolidation experiment of a shallow formation.

42 ENGINEERING↗

NASA NTRS1971

Similarity analysis of differential equations by Lie group.

Methods for transforming partial differential equations into forms more suitable for analysis and solution are investigated. The idea of Lie's infinitesimal contact transformation group is introduced to develop a systematic method which involves mostly algebraic manipulations. A thorough presentation of the application of this general method to the problem of similarity analysis in a broader sense - namely, the similarity between partial and ordinary differential equations, boundary value and initial value problems, and nonlinear and linear equations - is given with new and very general methods evolved for deriving the possible groups of transformations.

Na, T. Y.↗

NASA NTRS1992

Reduced basis technique for calculating sensitivity coefficients of nonlinear structural response

An efficient reduced basis technique is presented for calculating the sensitivity of nonlinear structural response to variations in the design variables. The structure is discretized by using two-field mixed finite element models. The vector of structural response and its sensitivity coefficients (derivatives with respect to design variables) are each expressed as a linear combination of a small number of basis (or global approximation) vectors. The Bubnov-Galerkin technique is then used to approximate each of the finite element equations governing the response and the sensitivity coefficients by a small number of algebraic equations in the amplitudes of these vectors. The path derivatives (derivatives of the response vector with respect to path parameters, e.g., load parameters) are used as basis vectors for approximating the response. A combination of the path derivatives and their derivatives with respect to the design variables is used for approximating the sensitivity coefficients. The potential of the proposed technique is discussed and its effectiveness is demonstrated by means of numerical examples of laminated composite plates subjected to mechanical and thermal loads.

Noor, Ahmed K.↗

NASA NTRS1996

Active Sonic Boom Control

A theory and simulation code are developed to study nonsteady sources as means to control sonic booms of supersonic aircraft. A key result is that the source of sonic boom pressure is not confined to the length of the aircraft but occupies an extensive segment of flight path. An aircraft in nonsteady flight functions as a synthetic aperture antenna, generating complex acoustic waves with no simple relation to instantaneous volume or lift distributions. The theory applies linear acoustics to slender nonsteady sources but requires no far field approximation. The solution for pressure contains a term not seen in Whitham's theory for sonic booms of distant supersonic aircraft. The term describes a pressure field that decays algebraically behind the Mach cone and, in the case of steady flight, integrates to a ground load equal to the weight of the aircraft. The algebraic term is separate from those that describe the sonic boom. Two nonsteady source phenomena are evaluated: periodic velocity changes (surge), and periodic longitudinal lift redistribution (slosh). Surge can attenuate a sonic boom and convert it into prolonged weak reverberation, but accelerations needed to produce the phenomenon seem too large for practical use. Slosh may be practical and can alter sonic booms but does not, on the average, result in boom attenuation. The conclusion is that active sonic boom abatement is possible in theory but maybe not practical

Crow, Steven C.↗

DOE OSTI2024

Dualities among massive, partially massless and shift symmetric fields on (A)dS

We catalog all the electromagnetic-like dualities that exist between free dynamical bosonic fields of arbitrary symmetry type and mass on (anti-) de Sitter space in all dimensions, including dualities among the partially massless and shift symmetric fields. This generalizes to all these field types the well known fact that a massless p-form is dual to a massless (D − p − 2)-form in D spacetime dimensions. In the process, we describe the structure of the Weyl modules (the spaces of local operators linear in the fields and their derivative relations) for all the massive, partially massless and shift symmetric fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2011

MueLu Multigrid Preconditioning Package

MueLu is intended for the research and development of multigrid algorithms used in the solution of sparse linear systems arising from systems of partial differential equations. The software provides multigrid source code, test programs, and short example programs to demonstrate the various interfaces for creating, accessing, and applying the solvers. MueLu currently provides an implementation of smoothed aggregation algebraic multigrid method and interfaces to many commonly used smoothers. However, the software is intended to be extensible, and new methods can be incorporated easily. MueLu also allows for advanced usage, such as combining multiple methods and segregated solves. The library supports point and block access to matrix data. All algorithms and methods in MueLu have been or will be published in the open scientific literature.

Tuminaro, Raymond↗

NASA NTRS1976

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

NASA NTRS1985

An efficient algorithm for estimating noise covariances in distributed systems

An efficient computational algorithm for estimating the noise covariance matrices of large linear discrete stochatic-dynamic systems is presented. Such systems arise typically by discretizing distributed-parameter systems, and their size renders computational efficiency a major consideration. The proposed adaptive filtering algorithm is based on the ideas of Belanger, and is algebraically equivalent to his algorithm. The earlier algorithm, however, has computational complexity proportional to p to the 6th, where p is the number of observations of the system state, while the new algorithm has complexity proportional to only p-cubed. Further, the formulation of noise covariance estimation as a secondary filter, analogous to state estimation as a primary filter, suggests several generalizations of the earlier algorithm. The performance of the proposed algorithm is demonstrated for a distributed system arising in numerical weather prediction.

Dee, D. P.↗

NASA NTRS1987

Stochastic control of infinite dimensional systems in Hilbert space: A factorization perspective

A factorization perspective on problems of optimal causal estimation and optimal causal control of linear stochastic systems defined on an infinite-dimensional Hilbert space is presented. A separation principle is derived for the case in which the system input/output map is generated by an abstract evolution operator. The factorization formalism allows for an essentially algebraic approach to these problems.

Milman, Mark M.↗

DOE OSTI2012

MueLu Multigrid Preconditioning Package

Source record↗

DOE OSTI2012

Non-linear realizations of conformal symmetry and effective field theory for the pseudo-conformal universe

The pseudo-conformal scenario is an alternative to inflation in which the early universe is described by an approximate conformal field theory on flat, Minkowski space. Some fields acquire a time-dependent expectation value, which breaks the flat space so(4,2) conformal algebra to its so(4,1) de Sitter subalgebra. As a result, weight-0 fields acquire a scale invariant spectrum of perturbations. The scenario is very general, and its essential features are determined by the symmetry breaking pattern, irrespective of the details of the underlying microphysics. In this paper, we apply the well-known coset technique to derive the most general effective lagrangian describing the Goldstone field and matter fields, consistent with the assumed symmetries. The resulting action captures the low energy dynamics of any pseudo-conformal realization, including the U(1)-invariant quartic model and the Galilean Genesis scenario. We also derive this lagrangian using an alternative method of curvature invariants, consisting of writing down geometric scalars in terms of the conformal mode. Using this general effective action, we compute the two-point function for the Goldstone and a fiducial weight-0 field, as well as some sample three-point functions involving these fields.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗

DOE OSTI2018

Interference in electron–molecule elastic scattering

General formulas describing the multiple scattering of electron by polyatomic molecules have been derived within the framework of the model of non-overlapping atomic potentials. These formulas are applied to different carbon molecules, both for fixed-in-space and randomly oriented molecules. Here, the molecular continuum wave function is represented as a plane wave plus a linear combination of the Green’s functions for free motion and the derivatives of these functions. Far from the target the electron spherical waves interfere as in the case of the Young’s double slits experiment. This interference manifests itself as diffraction oscillations in the differential and total cross sections for elastic electron scattering. The amplitude of electron scattering by a molecule is defined by the phase shifts for each of the atoms forming the target and its geometry. The numerical calculation of the scattering amplitude in closed form (rather than in the form of S-matrix expansion) is reduced to solving a system of algebraic equations. The number of atoms in a molecule and the atomic phase shifts for the orbital angular momenta included in the calculation define the number of equations.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1972

Heat transfer by condensation of low pressure metal vapors.

The film condensation of low pressure metal vapors on isothermal vertical flat plates or tubes is considered. The liquid film is treated as a thin layer in which the acceleration and pressure forces are negligible and across which the temperature distribution is linear. The average behavior of the vapor is found from the linearized one-dimensional vapor flow equations. In order to calculate the rate of condensation, a consistent distribution function for the vapor particles at the liquid-vapor interface is necessary and is determined. The result of the analysis is a set of algebraic equations from which one can predict the condensation rate of low pressure metal vapors. A large but continuous temperature decrease in the vapor is predicted and calculated.

Huang, Y. S.↗

NASA NTRS1980

A tensor approach to modeling of nonhomogeneous nonlinear systems

Model following control methodology plays a key role in numerous application areas. Cases in point include flight control systems and gas turbine engine control systems. Typical uses of such a design strategy involve the determination of nonlinear models which generate requested control and response trajectories for various commands. Linear multivariable techniques provide trim about these motions; and protection logic is added to secure the hardware from excursions beyond the specification range. This paper reports upon experience in developing a general class of such nonlinear models based upon the idea of the algebraic tensor product.

Yurkovich, S.↗

NASA NTRS1992

Analysis and experimental evaluation of a Stewart platform-based force/torque sensor

The kinematic analysis and experimentation of a force/torque sensor whose design is based on the mechanism of the Stewart Platform are discussed. Besides being used for measurement of forces/torques, the sensor also serves as a compliant platform which provides passive compliance during a robotic assembly task. It consists of two platforms, the upper compliant platform (UCP) and the lower compliant platform (LCP), coupled together through six spring-loaded pistons whose length variations are measured by six linear voltage differential transformers (LVDT) mounted along the pistons. Solutions to the forward and inverse kinematics of the force sensor are derived. Based on the known spring constant and the piston length changes, forces/torques applied to the LCP gripper are computed using vector algebra. Results of experiments conducted to evaluate the sensing capability of the force sensor are reported and discussed.

Nguyen, Charles C.↗