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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1967

Matrix Analysis of Longitudinal and Torsional Vibrations in Nonuniform Multibranch Beams

Since longitudinal modes and frequencies provide basic data for dynamic analyses of arbitrary beam-like structures and since closed-form solutions for the modes are generally not feasible to obtain, an approximate method is developed for computing the natural frequencies and the corresponding mode shapes for a variable-section, unconstrained multibranch beam. A lumped mass analogy employing influence coefficients is used to represent the beam. The simultaneous equations of motion for the lumped mass system are derived in matrix form and algebraically manipulated to yield a classical eigenvalue equation solvable by standard procedures. The orthogonality relationship of the natural modes is derived and used to form the basis of an orthogonal sweeping process for determination of modes above the fundamental. Numerical examples including an application to a solid-fuel launch system are presented. Also, a detailed discussion is devoted to the theoretical verifications of the approximate modes and frequencies.

Wingate, Robert T.↗

NASA NTRS2001

Durability and Damage Development in Woven Ceramic Matrix Composites

Damage development in woven SiC/SiNC ceramic matrix composites (CMC's) under tensile and cyclic loading both at room and elevated temperatures have been investigated for the exhaust nozzle of high-efficient turbine engines. The ultimate strength, failure strain, proportional limit and modulus data at a temperature range of 23 to 1250 C are generated. The tensile strength of SiC/SiNC woven composites have been observed to increase with increased temperatures up to 1000 C. The stress/strain plot shows a pseudo-yield point at 25 percent of the failure strain (epsilon(sub r)) which indicates damage initiation in the form of matrix cracking. The evolution of damage beyond 0.25 epsilon(sub f), both at room and elevated temperature comprises multiple matrix cracking, interfacial debonding, and fiber pullout. Although the nature of the stress/strain plot shows damage-tolerant behavior under static loading both at room and elevated temperature, the life expectancy of SiC/SiNC composites degrades significantly under cyclic loading at elevated temperature. This is mostly due to the interactions of fatigue damage caused by the mechanically induced plastic strain and the damage developed by the creep strain. The in situ damage evolutions are monitored by acoustic event parameters, ultrasonic C-scan and stiffness degradation. Rate equations for modulus degradation and fatigue life prediction of ceramic matrix composites both at room and elevated temperatures are developed. These rate equations are observed to show reasonable agreement with experimental results.

Haque, A.↗

DOE OSTI2025

Erratum: Optical absorption and luminescence of 𝛼−LiV 2 ⁢O 5 from the Bethe-Salpeter equation [Phys. Rev. B 112, 035132 (2025)]

This is an erratrum to a previously published paper, reporting that the previously made choice of the method of calculating the optical matrix elements was inaccurate.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Parallel algorithms for computation of the manipulator inertia matrix

The development of an O(log2N) parallel algorithm for the manipulator inertia matrix is presented. It is based on the most efficient serial algorithm which uses the composite rigid body method. Recursive doubling is used to reformulate the linear recurrence equations which are required to compute the diagonal elements of the matrix. It results in O(log2N) levels of computation. Computation of the off-diagonal elements involves N linear recurrences of varying-size and a new method, which avoids redundant computation of position and orientation transforms for the manipulator, is developed. The O(log2N) algorithm is presented in both equation and graphic forms which clearly show the parallelism inherent in the algorithm.

Amin-Javaheri, Masoud↗

NASA NTRS1980

On time discretizations for spectral methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.

Gottlieb, D.↗

DOE OSTI2021

Spatiotemporal Crossover between Low- and High-Temperature Dynamical Regimes in the Quantum Heisenberg Magnet

The stranglehold of low temperatures on fascinating quantum phenomena in one-dimensional quantum magnets has been challenged recently by the discovery of anomalous spin transport at high temperatures. Whereas both regimes have been investigated separately, no study has attempted to reconcile them. For instance, the paradigmatic quantum Heisenberg spin-$1/2$ chain falls at low-temperature within the Tomonaga-Luttinger liquid framework, while its high-temperature dynamics is superdiffusive and relates to the Kardar-Parisi-Zhang universality class in $1+1$ dimensions. This work aims at reconciling the two regimes. Building on large-scale matrix product state simulations, we find that they are connected by a temperature-dependent spatiotemporal crossover. In this work, as the temperature $T$ is reduced, we show that the onset of superdiffusion takes place at longer length and time scales $\propto 1/T$. This prediction has direct consequences for experiments including nuclear magnetic resonance: it is consistent with earlier measurements on the nearly ideal Heisenberg $S=1/2$ chain compound Sr$_2$CuO$_3$ yet calls for new and dedicated experiments.

1-dimensional spin chains↗

DOE OSTI2021

Characterizing Topological Order with Matrix Product Operators

One of the most striking features of gapped quantum phases that exhibit topological order is the presence of long-range entanglement that cannot be detected by any local order parameter. The formalism of projected entangled-pair states is a natural framework for the parameterization of gapped ground state wavefunctions which allows one to characterize topological order in terms of the virtual symmetries of the local tensors that encode the wavefunction. In their most general form, these symmetries are represented by matrix product operators acting on the virtual level, which leads to a set of algebraic rules characterizing states with topological quantum order. This construction generalizes the concepts of $\mathsf G$- and twisted injectivity; the corresponding matrix product operators encode all topological features of the theory and provide a complete picture of the ground state manifold on the torus. We show how the string-net models of Levin and Wen fit within this formalism and in doing so provide a particularly intuitive interpretation of the pentagon equation for F-symbols as the pulling of matrix product operators through the string-net tensor network. Our approach paves the way to finding novel topological phases beyond string nets and elucidates the description of topological phases in terms of entanglement Hamiltonians and edge theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1986

Simplified composite micromechanics for predicting microstresses

A unified set of composite micromechanics equations is summarized and described. This unified set is for predicting the ply microstresses when the ply stresses are known. The set consists of equations of simple form for predicting three-dimensional stresses (six each) in the matrix, fiber, and interface. Several numerical examples are included to illustrate use and computational effectiveness of the equations in this unified set. Numerical results from these examples are discussed with respect to their significance on microcrack formation and, therefore, damage initiation in fiber composites.

Chamis, C. C.↗

NASA NTRS1986

Simplified composite micromechanics for predicting microstresses

Chamis, Christos C.↗

NASA NTRS1987

Simplified composite micromechanics for predicting microstresses

Chamis, Christos C.↗

NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗

NASA NTRS1990

Melt migration in a silicate liquid-olivine system - An experimental test of compaction theory

An experimentally derived melt migration profile is compared with those obtained from compaction theory. A couple is prepared in which a source-for-melt (a disk composed of a silicate glass and olivine) was placed in contact with a sink-for-melt (a disk of polycrystalline olivine) in order to induce melt migration under laboratory conditions. It is noted that melt infiltrates into the sink disk along triple junctions driven by capillary forces. In order to analyze the melt migration profile in terms of compaction theory, the equations developed by McKenzie (1984) to describe porous flow of a liquid in a deformable matrix were modified. The bouyancy force term is replaced with a capillary force term and governing equations are then solved numerically with the initial boundary conditions specified by the experimental design. Also, it is determined that compaction theory provides a good description of the experimental results provided that the permeability of these solid-liquid materials increases linearly with increasing liquid fraction.

Riley, G. N., Jr.↗

DOE OSTI2022

On the novel 3-D neutron transport kinetic tRAPID algorithm and its validation

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, based on the Multi-stage Response-function Transport (MRT) methodology, allows for real-time simulation of nuclear systems based on 3-D continuous-energy particle transport. RAPID's steady-state (criticality) neutron transport algorithm is based on the Fission Matrix (FM) method, and has been extensively verified and validated against computational benchmarks and experiments. This paper introduces the novel 3-D time-dependent transport algorithm that has been implemented into the code, tRAPID, and its validation using the JSI TRIGA Mark-II reactor. tRAPID accurately and efficiently calculates neutron kinetics parameters (such as β{sub eff}, l{sub eff} , Λ, α{sub Rossi}) and 3-D time-dependent neutron fission source distribution and neutron importances for both prompt and delayed neutrons. tRAPID is used to simulate a rod insertion experiment performed at the JSI TRIGA Mark-II reactor, during which signals from four fission chambers at four different locations in the core were collected. The results demonstrate how tRAPID is capable of calculating detailed and accurate results with only a minimal use computational resources and time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2021

Multi-Frequency NMM3D Simulations of Wave Propagation in Vegetation for Remote Sensing of Soil Moisture

To investigate the feasibility of using multi-frequencies to fur- ther improve the soil moisture retrieval, the recently devel- oped hybrid method is used to perform full-wave scattering simulations of a wheat field at three different frequencies of L-, S- and C- bands. In the hybrid method, the multiple scattering within a single wheat plant is first captured using the T-matrix of vec- tor cylindrical waves calculated from the full-wave solutions of HFSS. In the second step, the scattering among different plants is considered using the Foldy-Lax equations of multiple scattering theory (MST) and the extracted T-matrix. The transmission of microwaves through wheat field at the L-, S- and C-bands are calculated using the hybrid method. Results show that: (1) the transmission obtained from the full-wave simulations is much larger than those computed from the radiative transfer equations (RTE) model, (2) the hybrid method transmission has a weaker frequency dependence than those of the RTE, and (3) the attenuation caused by the vegetation layer would saturate with an increase in frequency.

Yueh, Simon↗

NASA NTRS1998

Preconditioning the Helmholtz Equation for Rigid Ducts

An innovative hyperbolic preconditioning technique is developed for the numerical solution of the Helmholtz equation which governs acoustic propagation in ducts. Two pseudo-time parameters are used to produce an explicit iterative finite difference scheme. This scheme eliminates the large matrix storage requirements normally associated with numerical solutions to the Helmholtz equation. The solution procedure is very fast when compared to other transient and steady methods. Optimization and an error analysis of the preconditioning factors are present. For validation, the method is applied to sound propagation in a 2D semi-infinite hard wall duct.

Baumeister, Kenneth J.↗

NASA NTRS1989

An optimal momentum management controller for the Space Station

This paper presents a new sequential design procedure for determining an optimal CMG (control moment gyro) momentum management and attitude control system for the Space Station Freedom. First, the Space Station equations of motion are linearized and uncoupled, and the associated state space equations are defined. Next, a new sequential procedure is used for the development of a continuous quadratic regulator with eigenvalue placement in a specified region of the complex plane. The regional pole assignment method is utilized since it is best suited for tradeoffs between eigenvalue locations and robustness with respect to parameter variations, sensor failures, implementation accuracies and gain reductions. The matrix sign function is used for solving the Riccati equations which appear in the design procedure. Simulation results are given which show that the resultant design provides desired system performance.

Sunkel, J. W.↗

DOE OSTI2025

Bottomonium suppression in pNRQCD and open quantum system approach

By employing the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory within an open quantum system framework, we derive a Lindblad equation governing the evolution of the heavy-quarkonium reduced density matrix, accurate to next-to-leading order (NLO) in the ratio of the state's binding energy to the medium's temperature [1]. The derived NLO Lindblad equation provides a more reliable description of heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological applications, we numerically solve this equation using the quantum trajectories algorithm. By averaging over Monte Carlo-sampled quantum jumps, we obtain solutions without truncation in the angular momentum quantum number of the considered states. Our analysis highlights the importance of quantum jumps in the nonequilibrium evolution of bottomonium states within the quark-gluon plasma [2]. Additionally, we demonstrate that the quantum regeneration of singlet states from octet configurations is essential to explain experimental observations of bottomonium suppression. The heavy-quarkonium transport coefficients used in our study align with recent lattice QCD determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗