SEARCH · Engineering Papers

Results for “Matrix equations”

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 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗

DOE OSTI2020

Non-Hermitian approach for quantum plasmonics

Herein, we examine the limits of applicability of a simple non-Hermitian model for exciton/plasmon interactions in the presence of dissipation and dephasing. The model can be used as an alternative to the more complete Lindblad density matrix approach and is computationally and conceptually simpler. We find that optical spectra in the linear regime can be adequately described by this approach. The model can fail, however, under continuous optical driving in some circumstances. In the case of two quantum dots or excitons interacting with a plasmon, the model can also describe coherences and entanglement qualitatively when both dissipation and dephasing are present and quantitatively in the limit with no dephasing. The approach, within a single excitation manifold, is also applied to assess the role of disorder for 50 quantum dots interacting with a plasmon, where we find that, on average, large enough disorder can help stabilize the ensemble average of the open quantum system toward a dark quasi-steady-state much faster than without disorder. Furthermore, while such single excitation manifold calculations in this size limit can readily be done with either the non-Hermitian or Lindblad forms, as one goes to larger Hilbert space sizes, the computational and storage advantages of the non-Hermitian approach can become more useful. Published under license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Implementation of 3-D isoparametric finite elements on supercomputer for the formulation of recursive dynamical equations of multi-body systems

A computational algorithmic procedure is developed and implemented for the dynamic analysis of a multibody system with rigid/flexible interconnected bodies. The algorithm takes into consideration the large rotation/translation and small elastic deformations associated with the rigid-body degrees of freedom and the flexibility of the bodies in the system respectively. Versatile three-dimensional isoparametric brick elements are employed for the modeling of the geometric configurations of the bodies. The formulation of the recursive dynamical equations of motion is based on the recursive Kane's equations, strain energy concepts, and the techniques of component mode synthesis. In order to minimize CPU-intensive matrix multiplication operations and speed up the execution process, the concepts of indexed arrays is utilized in the formulation of the equations of motion. A spin-up maneuver of a space robot with three flexible links carrying a solar panel is used as an illustrative example.

Shareef, N. H.↗

DOE OSTI2020

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗

NASA NTRS1989

An integrating matrix formulation for buckling of rotating beams including the effects of concentrated masses

The integrating matrix technique of computational mechanics is extended to include the effects of concentrated masses. The stability of a flexible rotating beam with discrete masses is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. In this problem, the beam is subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam. To determine the eigenvalues from which stability is assessed, the differential equations of motion are solved numerically by combining the extended integrating matrix method with an eigenanalysis. Stability boundaries for a discrete mass representation of a uniform beam are shown to asymptotically approach the stability boundaries for the corresponding continuous mass beam as the number of concentrated masses is increased. An error in the literature is also noted for the discrete mass problem concerning the behavior of the critical rotational speed for inplane buckling as the radius of rotation of the clamped end of the beam is reduced.

Lakin, William D.↗

NASA NTRS1972

The calculation of the eigenvalues and eigenfunctions of Mathieu's equation

The eigenfunctions of Mathieu's equation are expanded in trigonometric series, and the resulting eigenvalue problem is cast in matrix form. This matrix is found to be a symmetric, triagonal matrix, and the eigenvalues are computed using the bisection method. The eigenfunction expansion coefficients are obtained by the standard recursion method. This computational technique for the eigenvalues and eigenfunctions of Mathieu's equation is both rapid and accurate.

Hodge, D. B.↗

NASA NTRS1973

A penny shaped crack in a filament-reinforced matrix. 2: The crack problem

The elastostatic interaction problem between a penny-shaped crack and a slender inclusion or filament in an elastic matrix was formulated. For a single filament as well as multiple identical filaments located symmetrically around the crack the problem is shown to reduce to a singular integral equation. The solution of the problem is obtained for various geometries and filament-to-matrix stiffness ratios, and the results relating to the angular variation of the stress intensity factor and the maximum filament stress are presented.

Pacella, A. H.↗

NASA NTRS1988

Laminate behavior for SiC fiber-reinforced reaction-bonded silicon nitride matrix composites

The room temperature mechanical properties of SiC fiber reinforced reaction-bonded silicon nitride matrix composite laminates (SiC/RBSN) have been measured. The laminates contained approx 30 volume fraction of aligned 142-micron diameter SiC fiber in a porous RBSN matrix. Three types of laminate studied were unidirectional: (1) (0) sub 8, (2) (10) sub 8, and (3) (45) sub 8, and (90) sub 8; cross plied laminates (0 sub 2/90 sub 2); and angle plied laminates: (+45 sub 2/-45 sub 2). Each laminate contained eight fiber plies. Results of the unidirectionally reinforced composites tested at various angles to the reinforcement direction indicate large anisotropy in in-plane properties. In addition, strength properties of these composites along the fiber direction were independent of specimen gage length and were unaffected by notches normal to the fiber direction. Splitting parallel to the fiber at the notch tip appears to be the dominant crack blunting mechanism responsible for notch insensitive behavior of these composites. In-plane properties of the composites can be improved by 2-D laminate construction. Mechanical property results for (0 sub 2/90 sub 2)sub s and (+45/-45 sub 2) sub s laminates showed that their matrix failure strains were similar to that for (0) sub 8 laminates, but their primary elastic moduli, matrix cracking strengths, and ultimate composite strengths were lower. The elastic properties of unidirectional, cross-ply, and angle-ply composites can be predicted from modified constitutive equations and laminate theory. Further improvements in laminate properties may be achieved by reducing the matrix porosity and by optimizing the bond strength between the SiC fiber and RBSN matrix.

Rhatt, R. T.↗

NASA NTRS1990

Laminate behavior for SiC fiber-reinforced reaction-bonded silicon nitride matrix composites

The room temperature mechanical properties of SiC fiber reinforced reaction-bonded silicon nitride matrix composite laminates (SiC/RBSN) have been measured. The laminates contained approx 30 volume fraction of aligned 142-micron diameter SiC fiber in a porous RBSN matrix. Three types of laminate studied were unidirectional: (1) (0) sub 8, (2) (10) sub 8, and (3) (45) sub 8, and (90) sub 8; cross plied laminates (0 sub 2/90 sub 2); and angle plied laminates: (+45 sub 2/-45 sub 2). Each laminate contained eight fiber plies. Results of the unidirectionally reinforced composites tested at various angles to the reinforcement direction indicate large anisotropy in in-plane properties. In addition, strength properties of these composites along the fiber direction were independent of specimen gage length and were unaffected by notches normal to the fiber direction. Splitting parallel to the fiber at the notch tip appears to be the dominant crack blunting mechanism responsible for notch insensitive behavior of these composites. In-plane properties of the composites can be improved by 2-D laminate construction. Mechanical property results for (0 sub 2/90 sub 2) sub s and (+45/-45 sub 2) sub s laminates showed that their matrix failure strains were similar to that for (0) sub 8 laminates, but their primary elastic moduli, matrix cracking strengths, and ultimate composite strengths were lower. The elastic properties of unidirectional, cross-ply, and angle-ply composites can be predicted from modified constitutive equations and laminate theory. Further improvements in laminate properties may be achieved by reducing the matrix porosity and by optimizing the bond strength between the SiC fiber and RBSN matrix.

Bhatt, Ramakrishna T.↗

NASA NTRS1991

Closed-form solutions for linear regulator-design of mechanical systems including optimal weighting matrix selection

This paper addresses the restriction of Linear Quadratic Regulator (LQR) solutions to the algebraic Riccati Equation to design spaces which can be implemented as passive structural members and/or dampers. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical systems. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist. Some examples of simple spring mass systems are shown to illustrate key points.

Hanks, Brantley R.↗

NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗

NASA NTRS1988

Optical matrix-vector processing for computational fluid dynamics

An optical processor to solve partial differential equations for computational fluid dynamics applications is considered. This application is new and original for optical processors. The algorithms that are used are optical realizations of the Newton-Raphson method for nonlinear equations and a new optical LU direct decomposition and Gauss-Seidel iterative solution to the resultant linear algebraic equations. These algorithms are used to solve Burger's equation (a specific form of the momentum equation in fluid dynamics). The nonlinear equations provide 1-D velocity data at each time step. Simulation results of optical processing with these algorithms on computational fluid dynamics data is included.

Perlee, Caroline J.↗

NASA NTRS1991

Simplified micromechanical equations for thermal residual stress analysis of coated fiber composites

The fabrication of metal matrix composites poses unique problems to the materials engineer. The large thermal expansion coefficient (CTE) mismatch between the fiber and matrix leads to high tensile residual stresses at the fiber/matrix (F/M) interface which could lead to premature matrix cracking during cooldown. Fiber coatings could be used to reduce thermal residual stresses. A simple closed form analysis, based on a three phase composite cylinder model, was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. Parametric studies showed that the tensile thermal residual stresses at the F/M interface were very sensitive to the CTE and thickness of the interphase layer. The modulus of the layer had only a moderate effect on tensile residual stresses. For a silicon carbide titanium aluminide composite, the tangential stresses were 20 to 30 pct. larger than the axial stresses, over a wide range of interphase layer properties, indicating a tendency to form radial matrix cracks during cooldown. Guidelines for the selection of appropriate material properties of the fiber coating were also derived in order to minimize thermal residual stresses in the matrix during fabrication.

Naik, R. A.↗

DOE OSTI2021

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Applying transfer matrix method to the estimation of the modal characteristics of the NASA Mini-Mass Truss

It is beneficial to use a distributed parameter model for large space structures because the approach minimizes the number of model parameters. Holzer's transfer matrix method provides a useful means to simplify and standardize the procedure for solving the system of partial differential equations. Any large space structures can be broken down into sub-structures with simple elastic and dynamical properties. For each single element, such as beam, tether, or rigid body, we can derive the corresponding transfer matrix. Combining these elements' matrices enables the solution of the global system equations. The characteristics equation can then be formed by satisfying the appropriate boundary conditions. Then natural frequencies and mode shapes can be determined by searching the roots of the characteristic equation at frequencies within the range of interest. This paper applies this methodology, and the maximum likelihood estimation method, to refine the modal characteristics of the NASA Mini-Mast Truss by successively matching the theoretical response to the test data of the truss. The method is being applied to more complex configurations.

Shen, Ji-Yao↗

NASA NTRS1993

Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

Walker, Kevin P.↗

NASA NTRS2008

Non-Boltzmann Modeling for Air Shock-Layer Radiation at Lunar-Return Conditions

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.

Johnston, Christopher O.↗

NASA NTRS1996

Performance of NASA Equation Solvers on Computational Mechanics Applications

This paper describes the performance of a new family of NASA-developed equation solvers used for large-scale (i.e. 551,705 equations) structural analysis. To minimize computer time and memory, the solvers are divided by application and matrix characteristics (sparse/dense, real/complex, symmetric/nonsymmetric, size: in-core/out of core) and exploit the hardware features of current and future computers. In this paper, the equation solvers, which are written in FORTRAN, and are therefore easily transportable, are shown to be faster than specialized computer library routines utilizing assembly code. Twenty NASA structural benchmark models with NASA solver timings reside on World Wide Web with a challenge to beat them.

Storaasli, Olaf O.↗