Engineering Papers⌕ Search

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.

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

At least 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2023

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1991

Review of LFTs, LMIs, and mu

The authors present a tutorial overview of linear fractional transformations (LFTs) and the role of the structured singular value, mu, and linear matrix inequalities (LMIs) in solving LFT problems. The authors first introduce the notation for LFTs and briefly discuss some of their properties. They then describe mu and its connections with LFTs. They focus on two standard notions of robust stability and performance, mu stability and performance and Q stability and performance, and their relationship is discussed. Comparisons with the L1 theory of robust performance with structured uncertainty are considered.

Doyle, John↗
NASA NTRS1993

Linear and nonlinear pattern selection in Rayleigh-Benard stability problems

A new algorithm is introduced to compute finite-amplitude states using primitive variables for Rayleigh-Benard convection on relatively coarse meshes. The algorithm is based on a finite-difference matrix-splitting approach that separates all physical and dimensional effects into one-dimensional subsets. The nonlinear pattern selection process for steady convection in an air-filled square cavity with insulated side walls is investigated for Rayleigh numbers up to 20,000. The internalization of disturbances that evolve into coherent patterns is investigated and transient solutions from linear perturbation theory are compared with and contrasted to the full numerical simulations.

Davis, Sanford S.↗
DOE OSTI2026

Affinity and Treatment of LDR Organics in Low Activity Waste by Cementitious Materials – 26109

The physiochemical interactions between organic contaminants, minerals, and polymer additives in grout/cementitious materials are being identified. This research supports a proposed Resource Conservation and Recovery Act (RCRA) treatment variance relevant to the Land Disposal Restriction (LDR) organics present in Hanford tank waste. It provides the technical basis to address the aspects documented by the Environmental Protection Agency in a draft memo necessary to establish that organics can be immobilized/stabilized in a cementitious matrix. The organic species may interact with cementitious minerals, including slag, fly ash, cement, and other components/dopants like activated carbon. A list of 132 reasonably expected LDR organics associated with Hanford tank waste (RPP-RPT-63493, Rev 1a) was screened by functional groups, octanol-water partitioning coefficients, and detection frequency in Hanford tank waste samples. A subset of 15 compounds spanning the identified properties was tested using sorption and leachate tests. These compounds were subjected to traditional batch sorption tests and leaching tests on/from Cast Stone cementitious material with and without activated carbon (a potential organic adsorption additive). The interactions between organic compounds were observed in traditional batch sorption and modified TCLP leachate tests for surrogate Cast Stone material with and without activated carbon addition. Specifically, phthalic acid (negatively charged) and 4-chloroaniline (polar) sorbed strongly to unmodified Cast Stone, while a group of three phenolic compounds (2,4,6-Trichlorophenol, Pentachlorophenol, o-Cresol) sorbed to and were retained by Cast Stone with a 1% activated carbon amendment. Mixtures of these compounds demonstrated visible retention of the organics within the matrix with a quantifiable affinity of the organics to the cement. This novel testing supports the notion that organics within Hanford tank waste would be retained onto the cement phase with the addition of simple materials easily added to the formulation as part of the treatment process. To further provide additional support to this claim, models employing Flory-Huggins theory in conjunction with macromolecular and organic molecule solubility parameters and species activity coefficients are being utilized to better understand these interactions. Within this context, Van Krevelen (VK) and Hansen organic and polymer solubility parameters were used to inform the Flory-Huggins interaction parameter, allowing for a quantifiable method of evaluating polymer-solvent compatibility and, by extension, interactions with proposed binding additives (e.g., activated carbon). These models and measurement should test the hypothesis for our systems understanding and applicability of Flory-Huggins theory, and whether the prediction of organics in the polymer phase versus the aqueous phase can be confidently made using this strategy. The estimation of the organic molecular activity coefficient in these high salt-saturated aqueous systems is another factor that will be assessed. This project aims to predict and understand complicated organic interactions in these heterogenous systems, guiding optimal treatment and stabilization methodologies.

Jolin, Will [Savannah River National Laboratory (S↗
DOE OSTI2022

One and Two Dimensional Quantum Lattice Algorithms for Maxwell Equations in Inhomogeneous Scalar Dielectric Media I: Theory

A quantum lattice algorithm (QLA) is developed for Maxwell equations in scalar dielectric media using the Riemann-Silberstein representation on a Cartesian grid. For x-dependent and y-dependent dielectric inhomogeneities, the corresponding QLA requires a minimum of 8 qubits/spatial lattice site. This is because the corresponding Pauli spin matrices have off-diagonal components which permit the local collisional entanglement of these qubits. However, z-dependent inhomogeneities require a QLA with a minimum of 16 qubits/lattice site since the Pauli spin matrix σz is diagonal. For 2 dimensional inhomogeneities, one can readily couple the 8-8 qubit schemes for x-y variations. z-x and y-z variations can be treated by either a 16-8 qubit scheme or a 16-16 qubit representation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48 Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. In this work, we use our scheme to compute the matrix element that governs the neutrinoless double beta decay of 48 Ca to 48 Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76 Ge, 130 Te, and 136 Xe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗
NASA NTRS1975

Boron/aluminum graphite/resin advanced fiber composite hybrids

Fabrication feasibility and potential of an adhesively bonded metal and resin matrix fiber-composite hybrid are determined as an advanced material for aerospace and other structural applications. The results show that using this hybrid concept makes possible a composite design which, when compared with nonhybrid composites, has greater transverse strength, transverse stiffness, and impact resistance with only a small penalty on density and longitudinal properties. The results also show that laminate theory is suitable for predicting the structural response of such hybrids. The sequence of fracture modes indicates that these types of hybrids can be readily designed to meet fail-safe requirements.

Chamis, C. C.↗
NASA NTRS1991

The generation of O(1S) from the dissociative recombination of O2(+)

The multichannel quantum defect theory (MQDT) method and large scale wave functions are applied to the calculation of the cross sections and rates for dissociative recombination of O2(+) along the 1Sigma-u(+) dissociative potential. Indirect dissociative recombination is accounted for by simultaneously including both the vibronic and electronic coupling to the intermediate Rydberg resonances. An enhanced MQDT approach involving a second-order K matrix is described. Cross sections and rates for the lowest three vibrational levels of the ion are reported. The shapes of the cross sections are discussed in terms of Fano's profile index. It is found that, for each of the three ion vibrational levels, the intermediate Rydberg resonances reduce the dissociative recombination rate below the direct recombination rate. Just above threshold, resonances with centers below threshold play an important role.

Guberman, Steven L.↗
NASA NTRS1992

The transition from columnar to equiaxed dendritic growth in proeutectic, low-volume fraction copper, Pb-Cu alloys

Lead, 17.1, 11.2, and 5 volume fraction copper (14, 9, and 4 wt pct Cu) alloys have been directionally solidified at constant growth velocities ranging from 1 to 100 microns/s. Serially increasing the growth velocity within this range results in a graded microstructural transition from fully columnar, albeit segregated, copper dendrites in a lead matrix to one consisting only of equiaxed grains. The imposed velocity necessary to effect fully equiaxed growth is found to drop rapidly as the volume fraction of copper is decreased. Factors which complicate the controlled, directional solidification of these alloys are discussed and the experimental results are interpreted in view of, and seen to be in qualitative agreement with, Hunt's theory on the transition from columnar to equiaxed growth of dendrites.

Kim, Shinwoo↗
NASA NTRS1994

The effects of physical aging at elevated temperatures on the viscoelastic creep on IM7/K3B

Physical aging at elevated temperature of the advanced composite IM7/K3B was investigated through the use of creep compliance tests. Testing consisted of short term isothermal, creep/recovery with the creep segments performed at constant load. The matrix dominated transverse tensile and in-plane shear behavior were measured at temperatures ranging from 200 to 230 C. Through the use of time based shifting procedures, the aging shift factors, shift rates and momentary master curve parameters were found at each temperature. These material parameters were used as input to a predictive methodology, which was based upon effective time theory and linear viscoelasticity combined with classical lamination theory. Long term creep compliance test data was compared to predictions to verify the method. The model was then used to predict the long term creep behavior for several general laminates.

Gates, Thomas S.↗
DOE OSTI2025

A flavor of SO(10) unification with a spinor Higgs

We investigate Higgs Parity unification — a realization of SO(10) grand unification based on the Higgs Parity mechanism in which the Standard Model (SM) Higgs resides in a spinor representation. The theory has an intermediate left-right symmetric stage where the SU(2)R symmetry breaking scale is fixed by the vanishing of the SM Higgs quartic coupling. The strong CP problem is solved by parity. Gauge coupling unification successfully predicts αs(MZ) to within 1%. The spinor Higgs naturally leads to a seesaw origin for SM flavor observables. We identify a novel mechanism where large mixing of third generation fermions with additional heavy vector-like fermions accounts for the anarchical nature of the PMNS matrix and the lack of hierarchy in the neutrino mass spectrum, relative to the up-quarks. A fit to quark and lepton masses and mixings, with a minimal parameter set, predicts 1) A testable relation between the top quark mass and αs(MZ) which is about (1 – 2)σ from current best fit values, 2) The order of magnitude of the baryon asymmetry of the universe, via leptogenesis from second-generation right-handed neutrino decays. 3) The proton decay and the neutron EDM are likely observable in next generation experiments, and 4) A normal ordered neutrino mass spectrum where 0νββ decay and the mass of the lightest neutrino are out of reach of next generation experiments.

Baryo-and Leptogenesis↗
DOE OSTI2024

Polarized 𝐽/𝜓 production in semi-inclusive DIS at large 𝑄 2 : Comparing quark fragmentation and photon-gluon fusion

We compare the relative importance of different mechanisms for polarized 𝐽/𝜓 production in semi-inclusive deep inelastic scattering processes at large 𝑄 2 . The transverse momentum dependent (TMD) factorization framework and nonrelativistic quantum chromodynamics are used to study the leading contributions from light quark fragmentation to polarized 𝐽/𝜓, and compared to direct production via photon-gluon fusion, which can proceed through color-singlet as well as color-octet mechanisms. We identify kinematic regimes where light quark fragmentation dominates, allowing for the extraction of the 3 𝑆$^{[8]}_{1}$ matrix element, as well as regimes where photon gluon fusion dominates, suggesting that the gluon TMD parton distribution function can be probed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Correspondence between theory and observations of polar motion

The focus of this study is to write in terms of reported values the standard theoretical equation describing the earth's wobble, namely, the linearized conservation of angular momentum equation known as the Liouville equation. A relation between the locations of the rotation pole and the celestial ephemeris pole is first obtained to rewrite this equation in terms of the location of the reported pole. This is achieved in the time domain by considering the properties of the time-dependent transformation matrix that relates components of a position vector in the rotating, body-fixed frame to its components in the celestial, space-fixed frame.

Gross, Richard S.↗
DOE OSTI2020

Theoretical study of the electronic, thermodynamic, and thermo-conductive properties of y-LiAlO2 with 6Li isotope substitutions for tritium production

Lithium aluminate has attracted researchers’ interests due to its wide applications. The structural, electronic, optical, lattice phonon thermodynamic, and thermo-conductive properties of ?-LiAlO2 with 6Li-isotope substitution are investigated by density functional theory and lattice phonon dynamics. The calculated results show that ?-LiAlO2 possesses a wide band-gap of 4.63 eV. The valence band is mainly formed by p orbitals of Li, O, and Al. The calculated lattice thermal conductivity of ?-LiAlO2 is found to be in good agreement over a wide range of temperatures with the available experimental data. The calculated dielectric matrix shows anisotropic behavior along z and x (or y) axes. Due to the wide band-gap, the optical conductivity of ?-LiAlO2 is mainly contributed from the electron hopping between valence-conduction bands. Substitution of 7Li with 6Li isotope in ?-LiAlO2 leads to observable differences in the lattice phonon frequencies at higher frequencies (>8THz) and to slight changes in corresponding thermal conductivity and the infrared and Raman spectra which are in good agreement with the measured data. The relevance of these properties for tritium production through 6Li absorbing neutron in nuclear reactors is discussed.

ab initio thermodynamics, ?-LiAlO2 and 6Li-isotope↗