Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix Theory”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Gravitational form factors of glueballs in Yang-Mills theory

This work presents preliminary results of the first determination of the energy-momentum tensor form factors of the scalar glueball, referred to as gravitational form factors (GFFs). The calculation has been carried out in lattice Yang-Mills theory at a single lattice spacing. Using variationally optimized operators, the matrix elements are extracted from ratios of three-point functions to two-point functions. The glueball GFFs and their kinematic dependence are compared to those of other hadrons from previous calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

ORACLS: A system for linear-quadratic-Gaussian control law design

A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.

Armstrong, E. S.↗
NASA NTRS1976

Residual stresses in polymer matrix composite laminates

Residual stresses in composites are induced during fabrication and by environmental exposure. The theory formulated can describe the shrinkage commonly observed after a thermal expansion test. Comparison between the analysis and experimental data for laminates of various material systems indicates that the residual stress-free temperature can be lower than the curing temperature, depending on the curing process. Effects of residual stresses on ply failure including the acoustic emission characteristics are discussed.

Hahn, H. T.↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1997

Exact-Output Tracking Theory for Systems with Parameter Jumps

We consider the exact output tracking problem for systems with parameter jumps. Necessary and sufficient conditions are derived for the elimination of switching-introduced output transient. Previous works have studied this problem by developing a regulator that maintains exact tracking through parameter jumps (switches). Such techniques are, however, only applicable to minimum-phase systems. In contrast, our approach is applicable to non-minimum-phase systems and it obtains bounded but possibly non-causal solutions. If the reference trajectories are generated by an exosystem, then we develop an exact-tracking controller in a feed-back form. As in standard regulator theory, we obtain a linear map from the states of the exosystem to the desired system state which is defined via a matrix differential equation. The constant solution of this differential equation provides asymptotic tracking, and coincides with the feedback law used in standard regulator theory. The obtained results are applied to a simple flexible manipulator with jumps in the pay-load mass.

Devasia, Santosh↗
NASA NTRS1989

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗
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.↗
NASA NTRS2001

The Oxidation Kinetics of Continuous Carbon Fibers in a Cracked Ceramic Matrix Composite

Experimental observations and results suggest two primary regimes as a function of temperature, i.e., diffusion and reaction controlled kinetics. Thermogravimetric analysis of carbon fiber in flowing oxygen gave an activation energy of 64.1 kJ/mol in the temperature range of 500 to 600 C and an apparent activation energy of 7.6 kJ/mol for temperatures from 600 to 1400 C. When C/SiC composite material was unstressed, matrix effects at temperatures from 900 to 1400 C protected the internal fibers. When under stress, self-protection was not observed. Increasing the stress from 10 to 25 ksi caused a 67 to 82 percent reduction in times to failure at temperatures from 750 to 1500 C. Based on experimental results, observation, and theory, a finite difference model was developed, which simulates the diffusion of oxygen into a matrix crack that is bridged by carbon fibers. The model allows the influence of important variables on oxidation kinetics to be studied systematically, i.e., temperature, reaction rate constant, diffusion coefficient, environment, and sample geometry.

Halbig, Michael C.↗
NASA NTRS1996

Exact-Output Tracking Theory for Systems with Parameter Jumps

In this paper we consider the exact output tracking problem for systems with parameter jumps. Necessary and sufficient conditions are derived for the elimination of switching-introduced output transient. Previous works have studied this problem by developing a regulator that maintains exact tracking through parameter jumps (switches). Such techniques are, however, only applicable to minimum-phase systems. In contrast, our approach is applicable to nonminimum-phase systems and obtains bounded but possibly non-causal solutions. If the reference trajectories are generated by an exo-system, then we develop an exact-tracking controller in a feedback form. As in standard regulator theory, we obtain a linear map from the states of the exo-system to the desired system state which is defined via a matrix differential equation. The constant solution of this differential equation provides asymptotic tracking, and coincides with the feedback law used in standard regulator theory. The obtained results are applied to a simple flexible manipulator with jumps in the pay-load mass.

Devasia, Santosh↗
NASA NTRS2003

Analytical Modeling of the High Strain Rate Deformation of Polymer Matrix Composites

The results presented here are part of an ongoing research program to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric matrix materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical 5 plasticity theory definitions of effective stress and effective plastic strain are modified by applying variations of the Drucker-Prager yield criterion. To verify the revised formulation, the shear and tensile deformation of a representative toughened epoxy is analyzed across a wide range of strain rates (from quasi-static to high strain rates) and the results are compared to experimentally obtained values. For the analyzed polymers, both the tensile and shear stress-strain curves computed using the analytical model correlate well with values obtained through experimental tests. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. In the micromechanics, the unit cell is divided up into a number of independently analyzed slices, and laminate theory is then applied to obtain the effective deformation of the unit cell. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite (composed using the representative polymer analyzed for the correlation of the polymer constitutive equations) for several fiber orientation angles across a variety of strain rates. The computed values compare favorably to experimentally obtained results.

Goldberg, Robert K.↗
DOE OSTI2026

Quantum Ising model on (2+1)-dimensional anti–de Sitter space using tensor networks

We study the quantum Ising model on (2+1)-dimensional anti-de Sitter space using matrix product states (MPS) and matrix product operators (MPOs). We explore the bulk phase diagram of the theory on regular tessellations of hyperbolic space with coordination number seven and find disordered and ordered phases separated by a phase transition. We find that the boundary-boundary spin correlation function exhibits power law scaling deep in the disordered phase of the Ising model consistent with holography. At the critical point, we find the boundary entanglement entropy scales logarithmically with subsystem size but away from this, we see a linear scaling. In comparison, the full system exhibits a volume law scaling, which is expected in chaotic and/or highly connected systems. We also measure out of time ordered correlators (OTOCs) to explore the scrambling behavior of the theory.

Quantum spin models↗
DOE OSTI2024

Evolution of Efimov states

The Efimov phenomenon manifests itself as an emergent discrete scaling symmetry in the quantum three-body problem. In the unitarity limit, it leads to an infinite tower of three-body bound states with energies forming a geometric sequence. Here in this work, we study the evolution of these so-called Efimov states using relativistic scattering theory. We identify them as poles of the three-particle S matrix and trace their trajectories in the complex energy plane as they evolve from virtual states through bound states to resonances. We dial the scattering parameters toward the unitarity limit and observe the emergence of the universal scaling of energies and couplings - a behavior known from the non-relativistic case. Interestingly, we find that Efimov resonances follow unusual, cyclic trajectories accumulating at the three-body threshold and then disappear at some values of the two-body scattering length. We propose a partial resolution to this ?missing states? problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

T -matrix analysis of static Wilson line correlators from lattice QCD at finite temperature

Here, we utilize a previously constructed thermodynamic T -matrix approach to the quark-gluon plasma (QGP) to derive constaints on its input by using results on Wilson line correlators (WLCs) of a static quark-antiquark pair from 2 + 1-flavor lattice-QCD (lQCD) computations with realistic pion mass. The self-consistent T -matrix results, which include previous constraints from the lQCD equation of state in the light-parton sector, can describe the lQCD data for WLCs fairly well once refinements of its driving kernel are applied. In particular, the input potential requires less screening than what has been inferred from previous T-matrix analyses. Pertinent predictions for the spectral and transport properties of the QGP are discussed, including the spatial diffusion coefficient for heavy quarks; the latter turns out to have a rather weak temperature dependence, in approximate agreement with recent unquenched lQCD results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Neutrinoless Double Beta Decay of 136 Xe and Related Nuclear Structure Studies (Final Scientific Report)

The search for neutrinoless double-beta decay (0nbb) explores new physics by directly probing the unknown mass scale and possible Majorana nature of the neutrino. The nEXO experiment, one of the two leading proposed ton-scale projects in the US, has a projected sensitivity to the 136 Xe 0nbb half-life of 10 28 years. The interpretation of possible signals in this next generation of experiments would be complicated by significant variations in theoretical calculations of the nuclear matrix elements (NMEs) of the decay, and nuclear structure measurements testing those theories can help address that uncertainty. We have 136 Xe(p,n) 136 Cs reaction at TUNL to deduce the level-scheme of states in 136 Cs through which the lowest-lying 1+ state decays. We find that over 99% of such decays will proceed through at least one isomeric state with a lifetime of order 100ns, which would enable large xenon detectors to employ a delayed-coincidence technique to search for charge-exchange processes including solar neutrino interactions. We have carried out measurements at TUNL of the 134 Xe( 3 He,n) 136 Ba and 136 Xe( 3 He,n) 138 Ba reactions to low-lying 0 + states in the residuals to probe the BCS assumption for QRPA NME calculations for the initial and final nuclei in 136 Xe 0nbb. While the analyses of these reactions is not yet complete and is ongoing, our initial results indicate tension with the BCS assumption. We have assembled a thermosyphon cooling R&D system for nEXO which has demonstrated up to 500 W of cooling in vacuum. We have developed a simple thermosyphon cooling simulation for nEXO and planned laboratory tests with the R&D system to benchmark it.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Application of Contraction Mappings to the Control of Nonlinear Systems

The theoretical and applied aspects of successive approximation techniques are considered for the determination of controls for nonlinear dynamical systems. Particular emphasis is placed upon the methods of contraction mappings and modified contraction mappings. It is shown that application of the Pontryagin principle to the optimal nonlinear regulator problem results in necessary conditions for optimality in the form of a two point boundary value problem (TPBVP). The TPBVP is represented by an operator equation and functional analytic results on the iterative solution of operator equations are applied. The general convergence theorems are translated and applied to those operators arising from the optimal regulation of nonlinear systems. It is shown that simply structured matrices and similarity transformations may be used to facilitate the calculation of the matrix Green functions and the evaluation of the convergence criteria. A controllability theory based on the integral representation of TPBVP's, the implicit function theorem, and contraction mappings is developed for nonlinear dynamical systems. Contraction mappings are theoretically and practically applied to a nonlinear control problem with bounded input control and the Lipschitz norm is used to prove convergence for the nondifferentiable operator. A dynamic model representing community drug usage is developed and the contraction mappings method is used to study the optimal regulation of the nonlinear system.

Killingsworth, W. R., Jr.↗
NASA NTRS1973

Application of modern control theory to the design of optimum aircraft controllers

The synthesis procedure presented is based on the solution of the output regulator problem of linear optimal control theory for time-invariant systems. By this technique, solution of the matrix Riccati equation leads to a constant linear feedback control law for an output regulator which will maintain a plant in a particular equilibrium condition in the presence of impulse disturbances. Two simple algorithms are presented that can be used in an automatic synthesis procedure for the design of maneuverable output regulators requiring only selected state variables for feedback. The first algorithm is for the construction of optimal feedforward control laws that can be superimposed upon a Kalman output regulator and that will drive the output of a plant to a desired constant value on command. The second algorithm is for the construction of optimal Luenberger observers that can be used to obtain feedback control laws for the output regulator requiring measurement of only part of the state vector. This algorithm constructs observers which have minimum response time under the constraint that the magnitude of the gains in the observer filter be less than some arbitrary limit.

Power, L. J.↗