Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix decomposition”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of Hydrogen Permeation Barrier and Matrix Interaction Protection Coatings for TZM Enclosure of Hydride Moderators (Microreactor Program M3 Report)

It is important to develop a high-performance enclosure solution for hydride-based microreactor moderation components to effectively suppress hydride thermal decomposition and prevent hydrogen loss in harsh in-reactor environments. Currently, the DOE-NE Microreactor Program (MR program) is investigating the titanium-zirconium-molybdenum (TZM) canned yttrium hydride (YH 2-x ) solution, which provides acceptable hydrogen retainment performance up to ~700 °C. However, improvements in the enclosure are necessary for increasing long-term operation temperature beyond 700 °C and enhancing short-term survivability at even higher temperatures during power transients. One potential solution is the addition of a dedicated hydrogen permeation barrier, which can potentially improve the enclosure's performance. Furthermore, the hydride moderator components will be inserted into the microreactor matrix, which may consist of materials that are not chemically inert against TZM. Therefore, the TZM might need to be protected from the external environment through the application of another surface protection coating. To address these critical issues, two types of barrier coatings optimized for TZM-based moderator enclosure has been developed. A multilayer coating approach has been used to identify TZM-compatible barrier materials. Behavior of these coating designs against H 2 permeation and long-term high-temperature graphite matrix interactions has been examined in detail. Additionally, the developed coatings' behavior when exposed to both high-temperature thermal cycling and heavy-ion irradiation has also been investigated. Overall, the development of these barrier coatings is essential for achieving maturity for the metal hydride moderator technology. The results in this report will provide crucial insights into the behavior of the developed coatings and their potential for improving the long-term performance of TZM-based moderator enclosures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗
DOE OSTI2026

Insights into Native Single-Atom Electrocatalyst Site Structures

Single-atom electrocatalysts consisting of metal atoms embedded in a carbon matrix are promising next-generation catalysts for green hydrogen production and utilization, CO2 reduction, low-temperature CO oxidation, ammonia production, plastic decomposition, and electrochemical energy storage. The origins of activity and stability for the single-atom sites are still debatable, however, because of constrained insights into their local structure resulting from idealized models and experiments derived from a large number of individual sites. Insights into structural variations around single atomic sites are therefore critical for the continued development of these next-generation catalysts. While electron microscopy commonly provides atomic-scale information about these materials, the beam sensitivity of individual sites makes structural determination by conventional low-voltage (60 keV) techniques challenging. Here, we introduce ultralow-voltage electron ptychography, performed at 30 keV, that enables determination of the lattice structure around individual metal sites in a well-defined single-atom electrocatalyst system while essentially eliminating knock-on structural modifications. Pairing these atomic-scale, site-specific measurements with computational methods will broaden our understanding of the activity and stability of these materials, which will accelerate the development of the next generation of catalysts.

Zachman, Michael [ORNL] (ORCID:0000000319101357)↗
DOE OSTI2021

Uniform Subspace Correction Preconditioners for Discontinuous Galerkin Methods with hp-Refinement

In this paper, we develop subspace correction preconditioners for discontinuous Galerkin (DG) discretizations of elliptic problems with hp-refinement. These preconditioners are based on the decomposition of the DG finite element space into a conforming subspace, and a set of small nonconforming edge spaces. The conforming subspace is preconditioned using a matrix-free low-order refined technique, which in this work, we extend to the hp-refinement context using a variational restriction approach. The condition number of the resulting linear system is independent of the granularity of the mesh h, and the degree of the polynomial approximation p. The method is amenable to use with meshes of any degree of irregularity and arbitrary distribution of polynomial degrees. Furthermore, numerical examples are shown on several test cases involving adaptively and randomly refined meshes, using both the symmetric interior penalty method and the second method of Bassi and Rebay (BR2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Matrix ensembles with global symmetries and ’t Hooft anomalies from 2d gauge theory

The Hilbert space of a quantum system with internal global symmetry $\textit{G}$ decomposes into sectors labelled by irreducible representations of $\textit{G}$. If the system is chaotic, the energies in each sector should separately resemble ordinary random matrix theory. We show that such “sector-wise” random matrix ensembles arise as the boundary dual of two- dimensional gravity with a $\textit{G}$ gauge field in the bulk. Within each sector, the eigenvalue density is enhanced by a nontrivial factor of the dimension of the representation, and the ground state energy is determined by the quadratic Casimir. We study the consequences of ’t Hooft anomalies in the matrix ensembles, which are incorporated by adding specific topological terms to the gauge theory action. The effect is to introduce projective representations into the decomposition of the Hilbert space. Finally, we consider ensembles with $\textit{G}$ symmetry and time reversal symmetry, and analyze a simple case of a mixed anomaly between time reversal and an internal $\mathbb{Z}$ 2 symmetry.

2D Gravity↗
DOE OSTI2022

When Do Extended Physics-Informed Neural Networks (XPINNs) Improve Generalization?

Physics-informed neural networks (PINNs) have become a popular choice for solving high-dimensional partial differential equations (PDEs) due to their excellent approximation power and generalization ability. Recently, extended PINNs (XPINNs) based on domain decomposition methods have attracted considerable attention due to their effectiveness in modeling multiscale and multiphysics problems and their parallelization. However, theoretical understanding of their convergence and generalization properties remains unexplored. In this study, we take an initial step towards understanding how and when XPINNs outperform PINNs. Specifically, for general multilayer PINNs and XPINNs, we first provide a prior generalization bound via the complexity of the target functions in the PDE problem and a posterior generalization bound via the posterior matrix norms of the networks after optimization. Moreover, based on our bounds, we analyze the conditions under which XPINNs improve generalization. Concretely, our theory shows that the key building block of XPINN, namely, the domain decomposition, introduces a tradeoff for generalization. On the one hand, XPINNs decompose the complex PDE solution into several simple parts, which decreases the complexity needed to learn each part and boosts generalization. On the other hand, decomposition leads to less training data being available in each subdomain, and hence such a model is typically prone to overfitting and may become less generalizable. Empirically, we choose five PDEs to show when XPINNs perform better than, similar to, or worse than PINNs, hence demonstrating and justifying our new theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Forced Oscillation Source Location in Power Systems Using MVMD-assisted DEF in TF Plane

Fast and accurate forced oscillation source location (FOSL) is critical for mitigating forced oscillations in power systems. Here, to overcome the shortcomings of traditional dissipating energy flow (DEF) in decomposing the forced oscillation components from the measurement responses, a multivariate variational mode decomposition (MVMD)-assisted FOSL is proposed in this paper to locate the forced oscillation sources from measurements in time-frequency (TF) plane. The multi-channel measurement matrix of each generator is formed, then the multi-channel modes are decomposed from the formed multi-channel measurement matrix using the MVMD. Further, the intrinsic mode functions (IMFs) associated with the forced oscillation mode are distinguished from the decomposed modes according to the relative energy weights. With the distinguished IMFs, the MVMD-assisted DEF calculation method is developed to locate the forced oscillation source. Finally, the performance of the proposed method is evaluated using the simulation data of the WECC 179-bus test system and IEEE-NASPI Oscillation Source Location Contest as well as the actual PMU data of ISO New England. The results validated the effectiveness and accuracy of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗
DOE OSTI2024

From phase decomposition to evaporation: A multi-modal evaluation of thermally degraded model lightweight high-entropy alloy

Lightweight high-entropy alloys (LHEAs) have the potential to replace conventional lightweight materials due to their superior mechanical properties and thermal stability. However, the thermal degradation pattern of LHEAs from phase decomposition to evaporation is not clear. We develop a new Al-based dual phase (FCC + HCP) LHEA—AlTi 0.45 CuZn, and further investigate its thermal degradation behavior for potential high-temperature structural applications. Using multimodal advanced characterization techniques such as differential scanning calorimetry/thermogravimetric analysis, scanning/transmission electron microscopy, and synchrotron X-ray diffraction/pair distribution function (XRD/PDF), a sequence of thermal degradation events beyond the thermal phase stability limit—between 250 and 360 °C—is observed. These include phase decomposition at ~360 °C, Zn evaporation at ~750 °C, and LHEA melting at 880 °C which results in ~25% cumulative weight loss. The formation of Al-Ti phase off the AlTi 0.45 CuZn matrix is due to the largest negative mixing enthalpy for Al-Ti than other binary pairs. Similarly, Zn evaporation from AlTi 0.45 CuZn LHEA is due to its faster evaporation rate than other constituent elements. The high-resolution synchrotron XRD and PDF results support the aforementioned observations; in addition, they reveal local atomic arrangements, local strain, and sluggish grain growth in the LHEA. Among other LHEAs of close density range (5.55 ≤ ρ ≤ 5.85 g/cc), the investigated LHEA exhibits outstanding nano-indentation hardness values due to the coupled grain size effect and HCP phase strengthening of the FCC matrix. As the search for LHEAs for lightweight applications grows, this study shows the potential use of AlTi 0.45 CuZn LHEA for structural applications even at elevated temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Tunable quantum logic gate on photonic qubits with a ladder emitter

We describe how a ladder emitter can implement a tunable quantum logic gate on photonic qubits encoded in the frequency basis. The ground-to-first excited state of the ladder emitter interacts with the control photon, and the first-to-second excited state transition interacts with the target photon. By controlling the relative detuning between the target photon and the first-to-second excited state transition of the ladder emitter, we enable any controlled-phase operation from 0 to π. We derive analytical formulas for the performance of the gate through the S-matrix formalism as well as describe the mechanism intuitively. This gate is deterministic, does not utilize any active control, and needs only a single ladder emitter, enabling low-footprint and more efficient decomposition of quantum circuits, especially the quantum Fourier transform. We suggest multiple potential systems for physical realization of our proposal, such as lanthanide ions embedded in Purcell-enhanced cavities. We expect these results to motivate further interest in photonic quantum information processing with designer emitters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Nonperturbative random matrix model of N = 2 JT supergravity

It is shown how to nonperturbatively define a random matrix model that captures key physics of N = 2 Jackiw-Teitelboim supergravity, going well beyond the perturbative topological expansion defined recently by Turiaci and Witten. A decomposition into an infinite family of certain multicritical models is derived, leading to the definition of a nonlinear ordinary differential equation from which the physics may be computed. Bogomol’nyi-Prasad-Sommerfield (BPS) states are naturally described by the model. The nonperturbative completions of the spectral densities for non-BPS multiplets are readily extracted. Published by the American Physical Society 2024

Johnson, Clifford V. (ORCID:0000000189645830)↗
DOE OSTI2020

Molecular Structure and Electron Affinity of Metal-Solvent Complexes: Insights from Density Functional Theory Simulations

A molecular level understanding of the structure and energetics of the monovalent and divalent metal ion complexes is of great importance for development of next-generation batteries. Here, Density Functional Theory (DFT) simulations at the ωb97xD/6-31 + G(d,p) level of theory are performed to investigate the interaction of metal ions (Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ ) with 26 organic solvent molecules. The reduction energetics (electron affinity and reduction potential) and structural responses of the solvent molecules and the molecular complexes are discussed. The DFT calculations are carried out to investigate the structure, energetics and electron affinities of chelated complexes of water (H 2 O), tetrahydrofuran (THF) and di-methoxy ethane (DME) solvent molecules. Additionally, ab initio dynamic simulations (AIMD) at 298 K using atom centered density matrix propagation (ADMP) formalism are performed to understand the spontaneous structure formation upon electron attachment of the metal ion-solvent complexes. The ADMP simulations indicate the decomposition of Mg + -(DME) 3 complex via cleavage of C–O bond of one of the three DME molecules indicating irreversible decomposition of DME in the presence of the Mg + radical. We believe that the data collected as part of this investigation serves as a library of fundamental knowledge towards a deeper understanding of the electrode-electrolyte interfacial reactions.

25 ENERGY STORAGE↗
DOE OSTI2024

Faster Tensor Network Decoding for Topological Quantum Codes

We present a fast and Bayes-optimal-approximating tensor network decoder for planar quantum LDPC codes based on the tensor renormalization group algorithm, originally proposed by Levin, and Nave. By precomputing the renormalization group flow for the null syndrome, we need only recompute tensor contractions in the causal cone of the measured syndrome at the time of decoding. This allows us to achieve an overall runtime complexity of ($pnχ^6$) where p is the depolarizing noise rate, and χ is the cutoff value used to control singular value decomposition approximations used in the algorithm. We apply our decoder to the surface code in the code capacity noise model and compare its performance to the original matrix product state (MPS) tensor network decoder introduced by Bravyi, Suchara, and Vargo. The MPS decoder has a p-independent runtime complexity of $\mathcal{O}(nχ^3)$ resulting in significantly slower decoding times compared to our algorithm in the low-p regime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Limiting Retained Austenite Decomposition in Quenched and Tempered Steels: Influences of Rapid Tempering and Silicon

Tempering reactions are critical to microstructure and property control in martensitic steels. Here, retained austenite decomposition and cementite precipitation are monitored using Mössbauer spectroscopy in 4340 and 300-M steel under conventional and rapid tempering conditions. Tempering times are compared at a constant tempered hardness by increasing tempering temperatures associated with short time conditions to achieve equivalent matrix softening to that of longer tempering times. Time-temperature combinations that provide equivalent tempered hardness generated microstructures with similar dislocation densities and cementite precipitation fractions; these mechanisms are controlled by self-diffusion. However, systematic differences in retained austenite content were observed at a given degree of softening, where shorter tempering times exhibited higher levels of retained austenite compared to more conventional conditions. At low temperatures, the differences in retained austenite preservation between explored time-temperature conditions are attributed to corresponding differences in carbon diffusion distance (in austenite), the controlling diffusional process of retained austenite decomposition. At higher temperatures, retained austenite decomposition exhibits C-curve kinetic behavior in 4340. Thus, reduced thermodynamic driving force for cementite and ferrite formation at higher temperature is believed to play a role in restricting retained austenite decomposition within some short-time, high temperature tempering regimes. The addition of silicon pushes cementite precipitation and retained austenite decomposition to higher temperatures, although retained austenite decomposition is suppressed to a greater extent than cementite precipitation. Potential is illustrated for coupling rapid tempering with silicon alloying to produce appreciably tempered martensite (~490 HV) with relatively less retained austenite decomposition compared to conventional tempering conditions.

36 MATERIALS SCIENCE↗