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 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Orbit averaging coherent states: holographic three-point functions of AdS giant gravitons

We study correlation functions of two AdS giant gravitons in AdS 5 × S 5 and a BPS supergravity mode using holography. In the gauge theory these are described by BPS correlators of Schur polynomials of fully-symmetric representations and a single trace operator. We find full agreement between the semiclassical gravity and gauge theory computations at large N, for both diagonal and off-diagonal structure constants. Our analysis in $\mathcal{N}$ = 4 SYM provides a simpler derivation to the results in the literature, and it can be readily generalized to operators describing bound states of AdS giant gravitons as well as bubbling geometries.

1/N Expansion↗

DOE OSTI2025

Non-perturbative quarkonium dissociation rates in strongly coupled quark-gluon plasma

Heavy quarks and quarkonia are versatile probes of the transport properties of the hot QCD medium produced in ultra-relativistic heavy-ion collisions (URHICs). A robust description of heavy-flavor transport coefficients requires a microscopic approach that treats the open and hidden heavy-flavor sectors on the same footing. Here, we employ the quantum many-body T -matrix formalism to evaluate the dissociation rates of heavy quarkonia in the quark-gluon plasma (QGP). The basic ingredient is the heavy-light T -matrix, which utilizes a nonperturbative driving kernel constrained by lattice-QCD data. Its resummation in a ladder series provides a much enhanced interaction strength compared to a previously used perturbative coupling to the quasiparticle partons in the QGP. The in-medium quarkonium properties, particularly their temperature-dependent binding energies, are obtained from selfconsistent calculations with the same interaction kernel, including interference effects (also referred to as the imaginary part of the heavy-quark potential) as well as off-shell parton spectral functions. We systematically investigate the interplay of these effects and elaborate on the connections to the dipole approximation used in effective field theory.

effective field theories of QCD↗

DOE OSTI2026

Optimality of Gradient-MUSIC for Spectral Estimation

We introduce the Gradient-MUSIC algorithm for estimating the unknown frequencies and amplitudes of a nonharmonic signal from noisy time samples. While the classical MUSIC algorithm performs a computationally expensive search over a fine grid, Gradient-MUSIC is significantly more efficient and eliminates the need for discretization over a fine grid by using optimization techniques. It coarsely scans the 1D landscape to find initialization simultaneously for all frequencies followed by parallelizable local refinement via gradient descent. We also analyze its performance when the noise level is sufficiently small and the signal frequencies are separated by at least 8π/m, where π/m is the standard resolution of this problem. Even though the 1D landscape is nonconvex, we prove a global convergence result for Gradient-MUSIC: coarse scanning provably finds suitable initialization and gradient descent converges at a linear rate. In addition to convergence results, we also upper bound the error between the true signal frequencies and amplitudes with those found by Gradient-MUSIC. For example, if the noise has $\ell^\infty$ norm at most ϵ, then the frequencies and amplitudes are recovered up to error at most Cϵ/m and Cϵ respectively, which are minimax optimal in m and ϵ. Our theory can also handle stochastic noise with performance guarantees under nonstationary independent Gaussian noise. Our main approach is a comprehensive geometric analysis of the landscape, a perspective that has not been explored before.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2014

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part I: Linear Theory

Conical shell theory and piston theory aerodynamics are used to study the aeroelastic stability of the thermal protection system (TPS) on the NASA Hypersonic Inflatable Aerodynamic Decelerator (HIAD). Structural models of the TPS consist of single or multiple orthotropic conical shell systems resting on several circumferential linear elastic supports. The shells in each model may have pinned (simply-supported) or elastically-supported edges. The Lagrangian is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the equations of motion. The natural modes of vibration and aeroelastic stability boundaries are found by calculating the eigenvalues and eigenvectors of a large coefficient matrix. When the in-flight configuration of the TPS is approximated as a single shell without elastic supports, asymmetric flutter in many circumferential waves is observed. When the elastic supports are included, the shell flutters symmetrically in zero circumferential waves. Structural damping is found to be important in this case. Aeroelastic models that consider the individual TPS layers as separate shells tend to flutter asymmetrically at high dynamic pressures relative to the single shell models. Several parameter studies also examine the effects of tension, orthotropicity, and elastic support stiffness.

Goldman, Benjamin D.↗

DOE OSTI2024

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗

DOE OSTI2022

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1986

Kalman filtering, smoothing and recursive robot arm forward and inverse dynamics

The inverse and forward dynamics problems for multi-link serial manipulators are solved by using recursive techniques from linear filtering and smoothing theory. The pivotal step is to cast the system dynamics and kinematics as a two-point boundary-value problem. Solution of this problem leads to filtering and smoothing techniques identical to the equations of Kalman filtering and Bryson-Frazier fixed time-interval smoothing. The solutions prescribe an inward filtering recursion to compute a sequence of constraint moments and forces followed by an outward recursion to determine a corresponding sequence of angular and linear accelerations. In addition to providing techniques to compute joint accelerations from applied joint moments (and vice versa), the report provides an approach to evaluate recursively the composite multi-link system inertia matrix and its inverse. The report lays the foundation for the potential use of filtering and smoothing techniques in robot inverse and forward dynamics and in robot control design.

Rodriguez, G.↗

NASA NTRS1989

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗

DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

Search for New Physics via EFT in collisions containing a top quark and a boosted Z or Higgs boson

We present a generator-level sensitivity study of physics beyond the Standard Model within the framework of Effective Field Theory (EFT). The targeted signal processes are single top quark production in association with a Lorentz-boosted Higgs or Z boson (tHq and tZq) in proton-proton collisions at a center-of-mass energy of 13.6 TeV. Event selection requires one electron or muon, a high-transverse-momentum, large-radius (AK8) jet reconstructing hadronic H → b̄b or Z → b̄b decays, at least two distinct, non-overlapping small-radius (AK4) jets with one originating from a b quark, and a missing transverse energy greater than 30 GeV. Signal and background samples are generated via Monte Carlo simulation using MadGraph with MadJax to perform EFT matrix-element calculations. Deviations from the Standard Model are parametrized through dimension-six EFT operators that probe modified top quark interactions. This generator-level study identifies kinematic observables with enhanced sensitivity to EFT operators and motivates future reconstruction-level analyses of tHq and tZq channels.

Meshramkar, Daniel [Baylor U., Waco]↗

DOE OSTI2026

Search for New Physics via EFT in collisions containing a top quark and a boosted Z or Higgs boson

Meshramkar, Daniel [Baylor U., Waco]↗

NASA NTRS2018

Coupled Thermo-Mechanical Micromechanics Modeling of the Influence of Thermally Grown Oxide Layer in an Environmental Barrier Coating System

Environmental Barrier Coatings (EBCs) have emerged as a promising means of protecting silicon based ceramic matrix composite (CMC) components for high temperature applications (e.g., aircraft engines). EBCs are often used to protect an underlying material (substrate) such as silicon carbide from extreme thermal/chemical environments. In a typical CMC/EBC system, an EBC may or may not be adhered to an underlying substrate with a bond coat (e.g., silicon). Irrespective, systems that utilize EBCs are susceptible to a number of failure modes including oxidation/delamination, recession, chemical attack and dissolution, thermomechanical degradation, erosion, and foreign object damage. Current work at NASA Glenn Research Center is aimed at addressing these failure modes in EBC systems and developing robust analysis tools to aid in the design process. The Higher-Order Theory for Functionally Graded Materials (HOTFGM), a precursor to the High-Fidelity Generalized Method of Cells micromechanics approach, was developed to investigate the coupled thermo-mechanical behavior of functionally graded composites and will be used herein to assess the development and growth of a low-stiffness thermally grown oxide (TGO) layer in EBC/CMC systems without a silicon bond coat. To accomplish this a sensitivity study is conducted to examine the influence of uniformly and nonuniformly grown oxide layer on the associated driving forces leading to mechanical failure (spallation) of EBC layer when subjected to isothermal loading.

thermally grown oxide↗

DOE OSTI2022

Capturing the ground state of uranium dioxide from first principles: Crystal distortion, magnetic structure, and phonons

Uranium dioxide (UO 2 ) remains a formidable challenge for first-principles approaches due to the complex interplay among spin-orbit coupling, Mott physics, magnetic ordering, and crystal distortions. Here we use DFT+U to explore UO 2 at zero temperature, incorporating all the aforementioned phenomena. The technical challenge is to navigate the many metastable electronic states produced by DFT+U, which is accomplished using f-orbital occupation matrix control to search for the ground state. We restrict our search to the high-symmetry ferromagnetic phase, including spin-orbit coupling, which produces a previously unreported occupation matrix. This newfound occupation matrix is then used as an initialization to explore the broken symmetry phases. We find the oxygen cage distortion of the 3k antiferromagnetic state to be in excellent agreement with experiments, and both the spin-orbit coupling and the Hubbard U are critical ingredients. Further, we demonstrate that only select phonon modes have a strong dependence on the Hubbard U, whereas magnetic ordering has only a small influence overall. We perform measurements of the phonon dispersion curves using inelastic neutron scattering, and our calculations show good agreement when using reasonable values of U. The quantitative success of DFT+U warrants exploration of thermal transport and other observables within this level of theory.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Modelling of thick composites using a layerwise laminate theory

The layerwise laminate theory of Reddy (1987) is used to develop a layerwise, two-dimensional, displacement-based, finite element model of laminated composite plates that assumes a piecewise continuous distribution of the tranverse strains through the laminate thickness. The resulting layerwise finite element model is capable of computing interlaminar stresses and other localized effects with the same level of accuracy as a conventional 3D finite element model. Although the total number of degrees of freedom are comparable in both models, the layerwise model maintains a 2D-type data structure that provides several advantages over a conventional 3D finite element model, e.g. simplified input data, ease of mesh alteration, and faster element stiffness matrix formulation. Two sample problems are provided to illustrate the accuracy of the present model in computing interlaminar stresses for laminates in bending and extension.

Robbins, D. H., Jr.↗

DOE OSTI2021

Theoretical study of CH 3 Cl-N 2 line shapes in the ν 1 band. Line mixing effects in QR doublets and QQ k sub-branches

Line coupling and line mixing effects of CH 3 Cl lines in the ν 1 band perturbed by N 2 have been investigated. We have taken into account the non-diagonality of the exp(–S 2 ) operator within specified line spaces as well as the k-degeneracy of the transitions (due to the double degeneracy of the j,k levels with k ≠ 0). These transitions should be considered as doublets coupled by the line mixing process. A new problem appears in the calculation when the atom-atom potential model is introduced. In order to overcome this difficulty, a pragmatic approach is proposed. Comparisons of theoretically calculated matrix elements of W with measurements of QR(j,k) doublets as well as some QQ k sub-branches, which are strongly affected by line mixing, have been made. Finally, the results show that the formalism improved in this way leads to rather accurate predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A Novel NASICON‐Type Na 4 MnCr(PO 4 ) 3 Demonstrating the Energy Density Record of Phosphate Cathodes for Sodium‐Ion Batteries

Abstract Sodium‐ion batteries (SIBs) have attracted incremental attention as a promising candidate for grid‐scale energy‐storage applications. To meet practical requirements, searching for new cathode materials with high energy density is of great importance. Herein, a novel Na superionic conductor (NASICON)‐type Na 4 MnCr(PO 4 ) 3 is developed as a high‐energy cathode for SIBs. The Na 4 MnCr(PO 4 ) 3 nanoparticles homogeneously embedded in a carbon matrix can present an extraordinary reversible capacity of 160.5 mA h g −1 with three‐electron reaction at ≈3.53 V during the Na + extraction/insertion process, realizing an unprecedentedly high energy density of 566.5 Wh kg −1 in the phosphate cathodes for SIBs. It is intriguing to reveal the underlying mechanism of the unique Mn 2+ /Mn 3+ , Mn 3+ /Mn 4+ , and Cr 3+ /Cr 4+ redox couples via X‐ray absorption near‐edge structure spectroscopy. The whole electrochemical reaction undergoes highly reversible single‐phase and biphasic transitions with a moderate volume change of 7.7% through in situ X‐ray diffraction and ex situ high‐energy synchrotron X‐ray diffraction. Combining density functional theory (DFT) calculations with the galvanostatic intermittent titration technique, the superior performance is ascribed to the low ionic‐migration energy barrier and desirable Na‐ion diffusion kinetics. The present work can offer a new insight into the design of multielectron‐reaction cathode materials for SIBs.

Zhang, Jian↗

DOE OSTI2020

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗