Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Decentralized Low-Rank State Estimation for Power Distribution Systems

This article considers the low-observability state estimation problem in power distribution networks and develops a decentralized state estimation algorithm leveraging the matrix completion methodology. Matrix completion has been shown to be an effective technique in state estimation that exploits the low dimensionality of the power system measurements to recover missing information. This technique can utilize an approximate (linear) load flow model, or it can be used with no physical models in a network where no information about the topology or line admittance is available. The direct application of matrix completion algorithms requires solving a semi-definite programming (SDP) problem, which becomes computationally challenging for large networks. We therefore develop a decentralized algorithm that capitalizes on the popular proximal alternating direction method of multipliers (proximal ADMM). The method allows us to distribute the computation among different areas of the network, leading to a scalable algorithm. By doing all computations at individual control areas and only communicating with neighboring areas, the algorithm eliminates the need for data to be sent to a central processing unit and thus increases efficiency and contributes to the goal of autonomous control of distribution networks. We illustrate the advantages of the proposed algorithm numerically using standard IEEE test cases.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

Deep Minima in the Triply Differential Cross Section for Ionization of Atomic Hydrogen by Electron and Positron Impact

We investigate ionization of atomic hydrogen by electron- and positron-impact. We apply the Coulomb–Born (CB1) approximation, various modified CB1 approximations, the three body distorted wave (3DW) approximation, and the time-dependent close-coupling (TDCC) method to electron-impact ionization of hydrogen. For electron-impact ionization of hydrogen for an incident energy of approximately 76.45 eV, we obtain a deep minimum in the CB1 triply differential cross section (TDCS). However, the TDCC for 74.45 eV and the 3DW for 74.46 eV gave a dip in the TDCS. For positron-hydrogen ionization (breakup) we apply the CB1 approximation and a modified CB1 approximation. We obtain a deep minimum in the TDCS and a zero in the CB1 transition matrix element for an incident energy of 100 eV with a gun angle of 56.13 °. Corresponding to a zero in the CB1 transition matrix element, there is a vortex in the velocity field associated with this element. For both electron- and positron-impact ionization of hydrogen the velocity field rotates in the same direction, which is anticlockwise. All calculations are performed for a doubly symmetric geometry; the electron-impact ionization is in-plane and the positron-impact ionization is out-of-plane.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗
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 OSTI2021

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Transparent window film with embedded nano-shades for thermoregulation

Every year, more than 20% of the energy consumption in the United States and more than 10% of the global consumption is used towards HVAC (Heating, Ventilating, Air-Conditioning) systems in buildings. Although continued efforts in advancing renewable energy, efficient appliances, and smart building systems are desired, the most critical factors that cause extensive indoor energy consumption are thermal leakage and thermal waste generation. We demonstrate a smart window film that allows natural lighting with a clear view while blocking oblique incident sunlight to make the building interior cooler without using electricity or generating heat. Building upon the concept of, often overlooked, traditional window shades and privacy window films, we emphasize and elevate the application and impact of using magnetically arranged nanoscale material for energy conservation purpose in buildings. By embedding vertically aligned Ni flakes (VANF) in a polymer matrix film, nanoscale mirror array serve as shades to block the high-angled incident sunlight, blocking approximately 73% of total solar irradiance and bringing down the indoor temperature by 1.92 °C during the daytime in a model house while providing nearly 90% transparency. Global application was simulated by applying the respective solar angle and intensity at selected latitudes, indicating their potential. This passively energy-conserving smart film is a low cost, maintenance free, and simple product for users and easy-to-manufacture product for the industry.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-Principles Phonon Quasiparticle Theory Applied to a Strongly Anharmonic Halide Perovskite

Understanding and predicting lattice dynamics in strongly anharmonic crystals is one of the long-standing challenges in condensed matter physics. Here we propose a first-principles method that gives accurate quasiparticle (QP) peaks of the phonon spectrum with strong anharmonic broadening. On top of the conventional first-order self-consistent phonon (SC1) dynamical matrix, the proposed method incorporates frequency renormalization effects by the bubble self-energy within the QP approximation. We apply the developed methodology to the strongly anharmonic α-CsPbBr 3 that displays phonon instability within the harmonic approximation in the whole Brillouin zone. While the SC1 theory significantly underestimates the cubic-to-tetragonal phase transition temperature (T c ) by more than 50%, we show that our approach yields T c = 404–423 K, in excellent agreement with the experimental value of 403 K. Furthermore, we also demonstrate that an accurate determination of QP peaks is paramount for quantitative prediction and elucidation of phonon linewidth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗
DOE OSTI2020

Performance of High Stopping Power Bismuth-Loaded Plastic Scintillators for Radiation Portal Monitors

Plastic scintillators are widely used in radiation portal monitors because of their low cost and availability in large sizes. However, due to their low density and low effective atomic number (Z), they offer low intrinsic efficiency and little spectroscopic information. The addition of high-Z constituents to these plastics can greatly increase both their total stopping power and the amount of photoelectric absorption, leading to full-energy deposition and thus useful gamma spectra. Herein, we present the performance of the latest formulation of Bi-loaded plastic scintillators showing their useful spectroscopic information up to relatively high energy (~1 MeV) due to their high stopping power compared to the current commercially available plastics. These Bi-loaded plastics use 20 weight percent (wt%) Bi-pivalate (8 wt% elemental Bi) dissolved in polyvinyltoluene (PVT) matrix and conventional fast fluors (~10 ns decay time). These Bi-loaded plastics achieve up to approximately 6000 photons/MeV and have been produced in sizes up to 17 in 3 . The performance of these Bi-loaded plastics is also demonstrated in the existing portal monitor hardware (Rapiscan Model TSA Trainer 770) showing the possibility to provide improved sensitivity as a drop-in replacement with continued scale-up.

42 ENGINEERING↗
DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Cyber-Resilient Frequency Control of Power Grids with Energy Storage Systems

The integration of synchronous generators and energy storage systems operated through communication networks introduces new challenges and vulnerabilities to the electric grid, where cyber attacks can corrupt sensor measurements or control inputs and interrupt functions such as frequency regulation. This paper proposes a defense methodology for the design of resilient operating constraints imposed on each generation and storage unit in order to prevent any attack sequence from driving the system's frequency to unsafe conditions. The resilient operating constraints are found by using ellipsoidal approximations of the reachable set of the power system, leading to a convex optimization problem with linear matrix inequalities. Numerical results in a single-area power system with synchronous generation and energy storage demonstrate how the resilient constraints provide security guarantees against any type of attack affecting frequency measurements or controller setpoints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗
DOE OSTI2023

Even spheres as joint spectra of matrix models

The Clifford spectrum is a form of joint spectrum for noncommuting matrices. This theory has been applied in photonics, condensed matter and string theory. In applications, the Clifford spectrum can be efficiently approximated using numerical methods, but this only is possible in low dimensional example. In this paper we examine the higher-dimensional spheres that can arise from theoretical examples. We also describe a constructive method to generate five real symmetric almost commuting matrices that have a K-theoretical obstruction to being close to commuting matrices. For this, we look to matrix models of topological electric circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first-principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of approximately 50 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. In conclusion, with greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

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

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

LLNL measurements of thermally irradiated HEU sample and saltwater-matrixed HEU sample

Two scintillation vials, kindly provided by the PNNL team, were received in the LLNL radiochemistry building on June 10, 2024; one containing approximately 2 mL dissolved thermally-irradiated HEU, and the other containing dried salt from irradiated seawater. A quick, semiquantitative screening count of the salt showed that the activation product activities were very low, with the highest activity being ~1 Bq 24 Na. The dissolved HEU solution was transferred to a Prindle vial (LLNL standard counting geometry) and weighed. Approximately 2 mL 3 M HNO 3 was added so that the solution completely covered the bottom of the vial and Gamanal software would be able to accurately generate an efficiency curve. This sample, called “PNNLFP24” is the irradiated HEU solution received, gravimetrically diluted by a factor of 1.7356 ± 0.0003, and was sent for quick gamma counting prior to further modification. Meanwhile the salt was quantitatively transferred to a 250 mL polyethylene bottle and dissolved in 160 mL 3 M HNO 3 . After counting, 50 µL PNNLFP24 was aliquoted for Resonance Ionization Mass Spectrometry (RIMS) analysis, and then 10.0 mL of the salt solution (6.03% of the total salt, gravimetrically) was added to the remaining PNNLFP24 solution to make “PNNLFPSW24”. After weighing, two 100 µL aliquots were removed for RIMS and DES, and 3 mL was aliquoted for microfluidic chemistry. The remaining “PNNLFPSW24” solution was weighed and proceeded for singles and coincidence counting on the MCBOS system.

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

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗