SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1990

A model for predicting damage induced fatigue life of laminated composite structural components

This paper presents a model for predicting the life of laminated composite structural components subjected to fatigue induced microstructural damage. The model uses the concept of continuum damage mechanics, wherein the effects of microcracks are incorporated into a damage dependent lamination theory instead of treating each crack as an internal boundary. Internal variables are formulated to account for the effects of both matrix cracks and internal delaminations. Evolution laws for determining the damage variables as functions of ply stresses are proposed, and comparisons of predicted damage evolution are made to experiment. In addition, predicted stiffness losses, as well as ply stresses are shown as functions of damage state for a variety of stacking sequences.

Allen, David H.↗

NASA NTRS1998

Attitude/attitude-rate estimation from GPS differential phase measurements using integrated-rate parameters

A sequential filtering algorithm is presented for attitude and attitude-rate estimation from Global Positioning System (GPS) differential carrier phase measurements. A third-order, minimal-parameter method for solving the attitude matrix kinematic equation is used to parameterize the filter's state, which renders the resulting estimator computationally efficient. Borrowing from tracking theory concepts, the angular acceleration is modeled as an exponentially autocorrelated stochastic process, thus avoiding the use of the uncertain spacecraft dynamic model. The new formulation facilitates the use of aiding vector observations in a unified filtering algorithm, which can enhance the method's robustness and accuracy. Numerical examples are used to demonstrate the performance of the method.

Oshman, Yaakov↗

NASA NTRS1997

Vector Observation-Aided/Attitude-Rate Estimation Using Global Positioning System Signals

Oshman, Yaakov↗

NASA NTRS2004

Photoionization of Benzene and Small Polycyclic Aromatic Hydrocarbons in Ultraviolet-Processed Astrophysical Ices: A Computational Study

We employed density functional theory (DFT) calculations to model the photoionization behavior of benzene and small polycyclic aromatic hydrocarbons when they are embedded in a matrix of water ice in order to investigate issues raised by recent experimental work by Gudipati and Allamandola. The ionization energies of benzene, naphthalene, anthracene, and pyrene were found to be lowered by 1.5-2.1 eV in water ice. Low-lying vertical electronic excitation energies were computed with time-dependent DFT for both neutral and ionized species and are found in both cases to be remarkably unaffected by the ice matrix. Chemical behavior in ultraviolet-photoprocessed ices is also discussed, with a focus on electron recombination and pathways leading to phenol and analogous products.

Woon, D. E.↗

DOE OSTI2026

Position Operators in Terms of Converging Finite-Dimensional Matrices and Their Intertwining with Geometry, Transport, and Gauge

The position operator 𝑟̂ appears as 𝑖∂ 𝑝 in wave mechanics, while its matrix form (e.g., under a Bloch basis) is well known diverging in diagonals, causing difficulties in basis transformation, observable yielding, etc. We aim to find a convergent r-matrix (CRM) to improve the existing divergent r-matrix (DRM), and investigate its influence at both the conceptual and the application levels. A key modification is increasing the familiar substitution of 𝑟̂ by 𝑖∂ 𝑝 to 𝑖∑ 𝑗 ∂ 𝑘 𝑗 , namely the N-th Weyl algebra. Resolving the divergence makes r-matrix rigorously defined, and we are able to show r-matrix is distinct from a spin matrix in terms of its defining principles, transformation behavior, and the observable it yields. Conceptually, the CRM fills the logical gap between the r-matrix and the Berry connection (this unremarked vagueness has caused the diagonal divergence). In application, we focus on transport, and discover that the Hermitian matrix is not identical with the associative Hermitian operator, i.e., 𝑟 𝑚,𝑛 = 𝑟$^∗_{𝑛,𝑚}$ ⇎ 𝑟̂ = 𝑟̂ † , which subtly affects the celebrated Berry curvature formula for adiabatic current. We also discuss how such a non-representation CRM can contribute to building a unified transport theory.

Weyl algebras↗

DOE OSTI2022

One-loop matching for quark dipole operators in a gradient-flow scheme

The quark chromoelectric dipole (qCEDM) operator is a CP-violating operator describing, at hadronic energies, beyond-the-standard-model contributions to the electric dipole moment of particles with nonzero spin. In this paper we define renormalized dipole operators in a regularization-independent scheme using the gradient flow, and we perform the matching at one loop in perturbation theory to renormalized operators of the same and lower dimension in the more familiar MS scheme. We also determine the matching coefficients for the quark chromo-magnetic dipole operator (qCMDM), which contributes for example to matrix elements relevant to CP-violating and CP-conserving kaon decays. The calculation provides a basis for future lattice QCD computations of hadronic matrix elements of the qCEDM and qCMDM operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Bootstrapping gauge theories

We consider asymptotically free gauge theories with gauge group S U ( N c ) and N f quarks with mass m q ≪ Λ QCD that undergo chiral symmetry breaking and confinement. We propose a bootstrap method to compute the S matrix of the pseudo-Goldstone bosons (pions) that dominate the low-energy physics. For the important case of N c = 3 , N f = 2 , a numerical implementation of the method gives the phase shifts of the S 0 , P 1 and S 2 waves in good agreement with experimental results. The method incorporates gauge theory information ( N c , N f , m q , Λ QCD ) by using the form-factor bootstrap recently proposed by Karateev, Kuhn and Penedones together with a finite energy version of the Shifman-Vainshtein-Zakharov (SVZ) sum rules. At low energy we impose constraints from chiral symmetry breaking. The only low-energy numerical inputs are the pion mass m π and the quark and gluon condensates. Published by the American Physical Society 2024

He, Yifei (ORCID:0000000213666157)↗

DOE OSTI2025

Xenon–metal pair formation in UO 2 investigated using DFT + U

A recent experimental study on a spent uranium dioxide (UO 2 ) fuel sample from Belgium Reactor 3 identified a unique pair structure formed by the noble metal phase (NMP) and fission gas [xenon (Xe)] precipitate. However, the fundamental mechanism behind this structure remains unclear. The present study aims to provide an understanding of the interaction between five different metal precipitates [molybdenum (Mo), ruthenium (Ru), palladium (Pd), technetium (Tc), and rhodium (Rh)] and the Xe fission gas atoms in UO 2 , by using density functional theory (DFT) in combination with the Hubbard U correction to compute the formation energies involved. All DFT + U calculations were performed with occupation matrix control to ensure antiferromagnetic ordering of UO 2 . The calculated formation and binding energies of the Xe and solid fission products in the NMP reveal that these metal precipitates form stable pair structures with Xe. Notably, the formation energy of Xe–metal pairs is lower than that of the isolated single defects in all instances, with Pd and Mo showing the most favorable binding energy, likely accounting for the observed pair structure formation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Neutral excitations of quantum Hall states: A density matrix renormalization group study

We use the dynamical structure factors of the quantum Hall states at ν = 1 / 3 and 1 / 2 in the lowest Landau level to study their excitation spectrum. Using the density matrix renormalization group in combination with the time-dependent variational principle on an infinite cylinder geometry, we extract the low-energy properties. At ν = 1 / 3 , a sharp magnetoroton mode and the two-roton continuum are present and the finite-size effects can be understood using the fractional charge of the quasiparticle. At ν = 1 / 2 , we find low-energy modes with linear dispersion and the static structure factor s ¯ ( q ) ~ ( q ℓ ) 3 in the limit q ℓ → 0 . The properties of these modes agree quantitatively with the predictions of the composite-fermion theory placed on the infinite cylinder.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Toward ab-initio nuclear theory calculations of 𝛿 𝐶

We propose a theory framework to study the isospin-symmetry breaking correction 𝛿 C in superallowed nuclear 𝛽 decays, crucial for the precise determination of |𝑉 𝑢⁢𝑑 |. Based on a general assumptions of the isovector dominance in isospin-symmetry breaking interactions, we construct a set of functions 𝐹 𝑇 𝑧 which involve nuclear matrix elements of isovector monopole operators and the nuclear Green's function. Via the functions 𝐹 𝑇 𝑧 , a connection of 𝛿 C to measurable electroweak nuclear radii is established, providing an experimental gauge of the theory accuracy of 𝛿 C . Here, we outline a strategy to perform ab initio calculations of 𝐹 𝑇 𝑧 based on the Lanczos algorithm, and discuss its similarity with other nuclear-structure-dependent inputs in nuclear 𝛽 decays.

Ab initio calculations↗

DOE OSTI2023

Three-body scattering and quantization conditions from S -matrix unitarity

Not Available

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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.↗

DOE OSTI2021

Approaching the self-dual point of the sinh-Gordon model

One of the most striking but mysterious properties of the sinh-Gordon model (ShG) is the b → 1/b self-duality of its S-matrix, of which there is no trace in its Lagrangian formulation. Here b is the coupling appearing in the model’s eponymous hyperbolic cosine present in its Lagrangian, cosh(bΦ). In this paper we develop truncated spectrum methods (TSMs) for studying the sinh-Gordon model at a finite volume as we vary the coupling constant. We obtain the expected results for b << 1 and intermediate values of b, but as the self-dual point b = 1 is approached, the basic application of the TSM to the ShG breaks down. We find that the TSM gives results with a strong cutoff Ec dependence, which disappears according only to a very slow power law in Ec. Standard renormalization group strategies — whether they be numerical or analytic — also fail to improve upon matters here. We thus explore three strategies to address the basic limitations of the TSM in the vicinity of b = 1. In the first, we focus on the small-volume spectrum. We attempt to understand how much of the physics of the ShG is encoded in the zero mode part of its Hamiltonian, in essence how ‘quantum mechanical’ vs ‘quantum field theoretic’ the problem is. In the second, we identify the divergencies present in perturbation theory and perform their resummation using a supra-Borel approximate. In the third approach, we use the exact form factors of the model to treat the ShG at one value of b as a perturbation of a ShG at a different coupling. In the light of this work, we argue that the strong coupling phase b > 1 of the Lagrangian formulation of model may be different from what is naïvely inferred from its S-matrix. In particular, we present an argument that the theory is massless for b > 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1980

Multicyclic control of a helicopter rotor considering the influence of vibration, loads, and control motion

Weighted multiple linear regression is used to establish a transfer function matrix relationship between higher harmonic control inputs and transducer vibration outputs for a controllable twist rotor. Data used in the regression were taken from the test of a KAMAN controllable twist rotor conducted in the Ames Research Center's 40- by 80-Foot Wind Tunnel in June 1977. Optimal controls to minimize fixed system vibrational levels are calculated using linear quadratic regulatory theory with a control deflection penalty included in the performance criteria. Control sensitivity to changes in control travel, forward speed, and lift and propulsive forces is examined. It is found that the linear transfer matrix is a strong function of forward speed and a weak function of lift and propulsive force. An open-loop strategy is proposed for systems with limited control travel.

Brown, T. J.↗

DOE OSTI2019

Improved reconstruction stability for chemical shift encoded hyperpolarized 13 C magnetic resonance spectroscopic imaging using k‐t spiral acquisitions

Purpose A multiecho, field of view (FOV)‐oversampled k‐t spiral acquisition and direct iterative decomposition of water and fat with echo asymmetry and least‐squares estimation reconstruction is demonstrated to improve the stability of hyperpolarized 13 C magnetic resonance spectroscopic imaging (MRSI) in the presence of signal ambiguities attributed to low‐SNR (signal‐to‐noise‐ratio) species, local uncertainties in metabolite peaks, and echo‐to‐echo signal inconsistencies. Theory k‐t spiral acquisitions redistribute readout points to be more densely spaced radially in k‐space by acquiring an FOV and matrix that are oversampled by η . These more densely spaced spiral turns constitute effective intraspiral echoes and can supplement conventional interspiral echoes to improve spectral separation and reduce spectral cross‐talk to better resolve 13 C‐labeled species for spectroscopic imaging. Methods Digital simulations and imaging phantom experiments were performed for a range of interspiral echo spacings and η using multiecho, k‐t spiral acquisitions. Image spectral cross‐talk artifacts were evaluated both qualitatively and quantitatively as the percent error in measured metabolite ratios. In vivo murine experiments evaluated the feasibility of multiecho, k‐t spiral [1‐ 13 C]pyruvate MRSI to reduce spectral cross‐talk for 3 scenarios of different expected reconstruction stability. Results Digital simulations and imaging phantom experiments both demonstrated reduced or comparable image spectral cross‐talk and percent errors in measured metabolite ratios with increasing η and better choices of echo spacings. In vivo images displayed markedly reduced spectral cross‐talk in lactate images acquired with η = 7 versus η = 1. Conclusion The precision of hyperpolarized 13 C metabolic imaging and quantification in the presence of low‐SNR species, local uncertainties in metabolite resonances, and echo‐to‐echo signal inconsistencies can be improved with the use of FOV‐oversampled k‐t spiral acquisitions.

Macdonald, Erin B.↗

DOE OSTI2023

Interactions of πK, ππK and KKπ systems at maximal isospin from lattice QCD

We study the interactions of systems of two and three nondegenerate mesons composed of pions and kaons at maximal isospin using lattice QCD, specifically $π^+K^+, π^+π^+K^+$ and $K^+K^+π^+$. Utilizing the stochastic LapH method, we determine the spectrum of these systems on two CLS $N_f$ = 2 + 1 ensembles with pion masses of 200 MeV and 340 MeV, and include many levels in different momentum frames. We constrain the K matrices describing two- and three-particle interactions by fitting the spectrum to the results predicted by the finite-volume formalism, including up to p waves. This requires also results for the $π^+π^+$ and $K^+K^+$ spectrum, which have been obtained previously on the same configurations. We explore different fitting strategies, comparing fits to energy shifts with fits to energies boosted to the rest frame, and also comparing simultaneous global fits to all relevant two- and three-particle channels to those where we first fit twoparticle channels and then add in the three-particle information. We provide the first determination of the three-particle K matrix in $π^+π^+K^+$ and $K^+K^+π^+$ systems, finding statistically significant nonzero results in most cases. We include s and $\textit{p}$ waves in the K matrix for $π^+K^+$ scattering, finding evidence for an attractive p-wave scattering length. We compare our results to Chiral Perturbation Theory, including an investigation of the impact of discretization errors, for which we provide the leading order predictions obtained using Wilson Chiral Perturbation Theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Lattice QCD calculation of the pion distribution amplitude with domain wall fermions at physical pion mass

We present a direct lattice QCD calculation of the x-dependence of the pion distribution amplitude (DA), which is performed using the quasi-DA in large momentum effective theory on a domain-wall fermion ensemble at physical quark masses and spacing a ≈ 0.084 fm. The bare quais-DA matrix elements are renormalized in the hybrid scheme and matched to $\overline{MS}$ with a subtraction of the leading renormalon in the Wilson-line mass. For the first time, we include threshold resummation in the perturbative matching onto the light-cone DA, which resums the large logarithms in the soft gluon limit at next-to-next-to-leading log. The resummed results show controlled scale-variation uncertainty within the range of momentum fraction x ϵ [0.25, 0.75] at the largest pion momentum P z ≈ 1.85 GeV. In addition, we apply the same analysis to quasi-DAs from a highly-improved-staggered-quark ensemble at physical pion mass and a = 0.076 fm. By comparison we find with 2σ confidence level that the DA obtained from chiral fermions is flatter and lower near x = 0.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Interlayer fusion bonding of semi-crystalline polymer composites in extrusion deposition additive manufacturing

This work focuses on the evolution of interlayer fracture toughness properties of fiber-reinforced, semi-crystalline polymers in the extrusion deposition additive manufacturing (EDAM) process. Further, this work bridges the gap between the additive process conditions (time-temperature history) and the effective layer-to-layer fracture properties developed within a printed component. This is the first step to predict delamination that can occur during printing, during cooling to room temperature after printing, and during service performance of an additively manufactured geometry. A phenomenological model is developed for fusion bonding of semi-crystalline polymer matrix composites by coupling the interdiffusion of polymer chains with the evolution of polymer crystallinity. While the interdiffusion is captured by reptation theory of polymer dynamics, the evolution of crystallinity is modeled by phenomenological crystallization kinetics and crystal melting dynamics. Further, a methodology is developed to determine the critical strain energy release rate, G IC of the interlayer interface and experiments are conducted utilizing the double cantilever beam fracture test geometry. In conclusion, predictions of GIC as a function of thermal history are compared with experiments.

36 MATERIALS SCIENCE↗