Engineering Papers⌕ Search

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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

On the stability of Galerkin methods for initial-boundary value problems for hyperbolic systems

The stability of approximating the solution of mixed initial-boundary value problems for hyperbolic systems by semidiscrete Galerkin methods is studied. It is shown that a particular straightforward Galerkin method yields an unstable approximation, and that this numerical instability is caused by an improper treatment of the boundary. Stable schemes are then presented, one of which differs from the unstable scheme only insofar as the treatment of the boundary is concerned. These stable schemes make use of a particular matrix which symmetrizes the differential system. It is therefore shown that the use of this matrix is crucial to the stability of the computations as well as for obtaining a priori bounds on the energy of the continuous system. This symmetrizing matrix is also related to the diagonalizing matrix for the system of hyperbolic equations and to the Liapunov matrix for the system of ordinary differential equations resulting from the application of Galerkin's method.

Gunzburger, M. D.↗
DOE OSTI2024

Not all that is β0 is β-function: the DGLAP resummation and the running coupling in NLO JIMWLK

Abstract We reanalyze the origin of the large transverse logarithms associated with the QCD one loopβfunction coefficient in the NLO JIMWLK Hamiltonian. We show that some of these terms are not associated with the running of the QCD coupling constant but rather with the DGLAP evolution. The DGLAP-like resummation of these logarithms is mandatory within the JIMWLK Hamiltonian, as long as the color correlation length in the projectile is larger than that in the target. This regime in fact covers the whole range of rapidities at which JIMWLK evolution is supposed to be applicable. We derive the RG equation that resums these logarithms to all orders inα s in the JIMWLK Hamiltonian. This is a nonlinear equation for the eikonal scattering matrixS(x). We solve this equation, and perform the DGLAP resummation in two simple cases: the dilute limit, where both the projectile and the target are far from saturation, and the saturated regime, where the target correlation length also determines its saturation momentum.

Physics↗
DOE OSTI2023

Determining the oxidation behavior of matrix graphite

This work presents the oxidation behavior of matrix graphite in air. Matrix graphite, graphite powder/flakes bonded by a small amount of non-graphitic carbon, surrounds coated fuel particles in order to form cylindrical fuel compacts (in prismatic core designs) or spheres (in pebble-bed reactor designs). This work focuses on oxidation tests conducted on two matrix graphite materials, one provided by Kairos Power and the other A3 matrix graphite. Some of the tests followed American Society for Testing and Materials (ASTM) oxidation testing standards using a vertical furnace system and others were performed in a thermogravimetric analyzer (TGA). It was determined that, at temperatures of 450 °C–700 °C, the oxidation rate of the Kairos matrix graphite follows the Arrhenius equation. In comparison with A3 matrix graphite, the Kairos matrix graphite shows better oxidation resistance at high temperatures (≥550 °C), but also a higher oxidation rate at low temperatures. Both the A3 matrix graphite and the Kairos matrix graphite materials may experience preferential oxidation of the partially graphitized binder. An oxygen penetration gradient was also observed when using the three characterization methods (i.e., optical microscope, x-ray tomography [XCT], and density profile by the lathe) enlisted in this research. In conclusion, the oxygen penetration depth increases with decreasing isothermal oxidation temperature, while the center of the oxidized samples (10% weight loss) remains almost untouched even at 500 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1974

Computer programs for calculation of matrix stability and frequency response from a state-space system description

FORTRAN computer subroutines stemming from requirements to process state variable system equations for systems of high order are presented. They find the characteristic equation of a matrix using the method of Danilevsky, the number of roots with positive real parts using the Routh-Horwitz alternate formulation, convert a state variable system description to a Laplace transfer function using the method of Bollinger, and evaluate that transfer function and obtain its frequency response. A sample problem is presented to demonstrate use of the subroutines.

Seidel, R. C.↗
DOE OSTI2020

Parallel-in-Time Simulation of Lindblad's Equation

Constructing fast quantum logic gates is critical to building a scalable quantum computer. We consider a qudit, a quantum version of a bit that can take an arbitrary number of states, coupled with a cavity. In this project, we wish to force the qudit to reach the 0-state, for any possible initial state. The coupled system changes in time according to Lindblad’s equation, an ordinary differential equation on the density matrix of the quantum system. Lindblad’s equation contains some parameters that we can control, so-called control functions. We seek control functions which force the qudit to the 0-state within 2 microseconds, which is much faster than what is currently done in practice. The search method is gradient descent, a numerical optimization method that uses gradient information to iteratively improve the control parameters. My contribution to this project is an attempt to speed up the computation of the gradient. It currently takes about 40 seconds to compute the gradient which involves solving a set of ODEs sequentially. Current supercomputers have thousands of cores, but sequential computations can only make use of 1 core at a time. We wish to divide up the work better, so that we can use many more cores at once. To this end, we have implemented the Multigrid Reduction in Time (MGRIT) algorithm. We perform a systematic parameter search on how to best apply this algorithm. Results indicate a 25 percent speed up for solving Lindblad’s equation and determining how close the final state is the 0-state.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

Edwards, Jack R.↗
NASA NTRS1982

A general algorithm for solving the algebraic Riccati equation

The generalized eigenvalue problem provides a suitable framework for reliable solutions of many system theoretic, control, and estimation problems. A general algorithm for solving the matrix algebraic Riccati equation (ARE) which utilizes a pencil structure is described here. This algorithm avoids unnecessary inversion of cost or transition matrices, making it a numerically sound way to solve for the gains and/or ARE with singular quadratic costs, for cases satisfying detectability and stabilizability conditions. Examples are solution with discrete dead-beat control, noiseless measurements in Kalman filters and time-delays in discrete-time systems, which cause difficulties in the Hamiltonian standard eigenvalue problem formulation. The ARE algorithm implementatiton and numerical examples are shown.

Walker, R. A.↗
NASA NTRS2002

A Numerical Scheme for Ordinary Differential Equations Having Time Varying and Nonlinear Coefficients Based on the State Transition Matrix

A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.

Bartels, Robert E.↗
NASA NTRS1993

Proposed framework for thermomechanical life modeling of metal matrix composites

The framework of a mechanics of materials model is proposed for thermomechanical fatigue (TMF) life prediction of unidirectional, continuous-fiber metal matrix composites (MMC's). Axially loaded MMC test samples are analyzed as structural components whose fatigue lives are governed by local stress-strain conditions resulting from combined interactions of the matrix, interfacial layer, and fiber constituents. The metallic matrix is identified as the vehicle for tracking fatigue crack initiation and propagation. The proposed framework has three major elements. First, TMF flow and failure characteristics of in situ matrix material are approximated from tests of unreinforced matrix material, and matrix TMF life prediction equations are numerically calibrated. The macrocrack initiation fatigue life of the matrix material is divided into microcrack initiation and microcrack propagation phases. Second, the influencing factors created by the presence of fibers and interfaces are analyzed, characterized, and documented in equation form. Some of the influences act on the microcrack initiation portion of the matrix fatigue life, others on the microcrack propagation life, while some affect both. Influencing factors include coefficient of thermal expansion mismatch strains, residual (mean) stresses, multiaxial stress states, off-axis fibers, internal stress concentrations, multiple initiation sites, nonuniform fiber spacing, fiber debonding, interfacial layers and cracking, fractured fibers, fiber deflections of crack fronts, fiber bridging of matrix cracks, and internal oxidation along internal interfaces. Equations exist for some, but not all, of the currently identified influencing factors. The third element is the inclusion of overriding influences such as maximum tensile strain limits of brittle fibers that could cause local fractures and ensuing catastrophic failure of surrounding matrix material. Some experimental data exist for assessing the plausibility of the proposed framework.

Halford, Gary R.↗
NASA NTRS1977

A study of the effect of radical load distributions on calibrated strain gage load equations

For several decades, calibrated strain gages have been used to measure loads on airplanes. The accuracy of the equations used to relate the strain gage measurements to the applied loads has been based primarily on the results of the load calibration. An approach is presented for studying the effect of widely varying load distributions on strain gage load equations. The computational procedure provides a link between the load calibration and the load to be measured in flight. A matrix approach to equation selection is presented, which is based on equation standard error, load distribution, and influence coefficient plots of the strain gage equations, and is applied to a complex, delta-wing structure.

Jenkins, J. M.↗
NASA NTRS1987

Effects of the oceans on polar motion: Extended investigations

Matrix formulation of the tide equations (pole tide in nonglobal oceans); matrix formulation of the associated boundary conditions (constraints on the tide velocity at coastlines); and FORTRAN encoding of the tide equations excluding boundary conditions were completed. The need for supercomputer facilities was evident. Large versions of the programs were successfully run on the CYBER, submitting the jobs from SUNY through the BITNET network. The code was also restructured to include boundary constraints.

Dickman, Steven R.↗
NASA NTRS1989

A family of permutations for concurrent factorization of block tridiagonal matrices

The inherent strong seriality of closely coupled systems is circumvented by defining a family of permutations for reordering equation sets whose matrix of coefficients is Hermitian block tridiagonal. The authors show how these permutations can be used to achieve relatively high concurrency in the Cholesky factorization of banded systems at the expense of introducing limited extra computations due to fill-in terms in the factors. Directed graphs are developed for the concurrent factorization of the transformed matrix of coefficients by the Cholesky algorithm. Expressions for speedup and efficiency are derived in terms of parameters of the permutation, set of equations, and machine architecture.

Utku, Senol↗
NASA NTRS1982

Integrating matrix formulations for vibrations of rotating beams including the effects of concentrated masses

By expressing partial differential equations of motion in matrix notation, utilizing the integrating matrix as a spatial operator, and applying the boundary conditions, the resulting ordinary differential equations can be cast into standard eigenvalue form upon assumption of the usual time dependence. As originally developed, the technique was limited to beams having continuous mass and stiffness properties along their lengths. Integrating matrix methods are extended to treat the differential equations governing the flap, lag, or axial vibrations of rotating beams having concentrated masses. Inclusion of concentrated masses is shown to lead to the same kind of standard eigenvalue problem as before, but with slightly modified matrices.

Lakin, W. D.↗
NASA NTRS1994

An implicit numerical scheme for the simulation of internal viscous flows on unstructured grids

The Navier-Stokes equations are solved numerically for two-dimensional steady viscous laminar flows. The grids are generated based on the method of Delaunay triangulation. A finite-volume approach is used to discretize the conservation law form of the compressible flow equations written in terms of primitive variables. A preconditioning matrix is added to the equations so that low Mach number flows can be solved economically. The equations are time marched using either an implicit Gauss-Seidel iterative procedure or a solver based on a conjugate gradient like method. A four color scheme is employed to vectorize the block Gauss-Seidel relaxation procedure. This increases the memory requirements minimally and decreases the computer time spent solving the resulting system of equations substantially. A factor of 7.6 speed up in the matrix solver is typical for the viscous equations. Numerical results are obtained for inviscid flow over a bump in a channel at subsonic and transonic conditions for validation with structured solvers. Viscous results are computed for developing flow in a channel, a symmetric sudden expansion, periodic tandem cylinders in a cross-flow, and a four-port valve. Comparisons are made with available results obtained by other investigators.

Jorgenson, Philip C. E.↗
NASA NTRS1989

Solving periodic block tridiagonal systems using the Sherman-Morrison-Woodbury formula

Many algorithms for solving the Navier-Stokes equations require the solution of periodic block tridiagonal systems of equations. By applying a splitting to the matrix representing this system of equations, it may first be reduced to a block tridiagonal matrix plus an outer product of two block vectors. The Sherman-Morrison-Woodbury formula is then applied. The algorithm thus reduces a periodic banded system to a non-periodic banded system with additional right-hand sides and is of higher efficiency than standard Thomas algorithm/LU decompositions.

Yarrow, Maurice↗
NASA NTRS1990

Mechanical behaviors of ceramic matrix composites with matrix cracking and fiber debonding

The purpose of this paper is to summarize the current research of the authors on the mechanical behaviors of ceramic matrix composites, including (1) the stress distributions in a composite with matrix cracking and interfacial debonding, (2) the critical strain for matrix cracking, and (3) the effects of fiber/matrix debonding and thermal residual stresses on the critical strain. The stress fields in both bonded and debonded regions are evaluated by taking into account thermal effects. An energy balance approach is followed to determine the critical strain for matrix cracking. From the general equation of the critical strain for matrix cracking, close form solutions have been obtained for two limiting cases: perfect bonding and complete debonding. Numerical solutions are given for the cases of partial fiber debonding and nonzero debonding energy. It is found that thermal residual stresses and the debonding energy have significant effects on the critical strain.

Kuo, Wen-Shyong↗
NASA NTRS1992

A three-dimensional nonlinear Timoshenko beam based on the core-congruential formulation

A three-dimensional, geometrically nonlinear two-node Timoshenkoo beam element based on the total Larangrian description is derived. The element behavior is assumed to be linear elastic, but no restrictions are placed on magnitude of finite rotations. The resulting element has twelve degrees of freedom: six translational components and six rotational-vector components. The formulation uses the Green-Lagrange strains and second Piola-Kirchhoff stresses as energy-conjugate variables and accounts for the bending-stretching and bending-torsional coupling effects without special provisions. The core-congruential formulation (CCF) is used to derived the discrete equations in a staged manner. Core equations involving the internal force vector and tangent stiffness matrix are developed at the particle level. A sequence of matrix transformations carries these equations to beam cross-sections and finally to the element nodal degrees of freedom. The choice of finite rotation measure is made in the next-to-last transformation stage, and the choice of over-the-element interpolation in the last one. The tangent stiffness matrix is found to retain symmetry if the rotational vector is chosen to measure finite rotations. An extensive set of numerical examples is presented to test and validate the present element.

Crivelli, Luis A.↗
NASA NTRS1984

Radiative transfer in spherical atmospheres

A method for defining spherical model atmospheres in radiative/convective and hydrostatic equilibrium is presented. A finite difference form is found for the transfer equation and a matrix operator is developed as the discrete space analog (in curvilinear coordinates) of a formal integral in plane geometry. Pressure is treated as a function of temperature. Flux conservation is maintained within the energy equation, although the correct luminosity transport must be assigned for any given level of the atmosphere. A perturbed integral operator is used in a complete linearization of the transfer and constraint equations. Finally, techniques for generating stable solutions in economical computer time are discussed.

Kalkofen, W.↗