Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced basis method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Accessing the gluon momentum fraction of nucleons through the gradient flow

We calculate the gluon momentum fraction of the nucleon using lattice QCD, with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using 𝑁 𝑓 =2 +1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the $\overline{MS}$-scheme at 2 GeV. Our final result is ⟨𝑥⟩ 𝑔 ⁢(𝜇 =2 GeV) =0.482⁢(35), where we quote only statistical uncertainties.

Lattice QCD↗
DOE OSTI2023

Canonical and noncanonical Hamiltonian operator inference

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Collaboration on SBIR/TTO Phase II Project for Optical QC Device and Collaboration on SBIR/TTO Phase II-b Project for Advanced Quality Inspection Device Development: Cooperative Research and Development CRADA Number CRD-16-00652 (Final Report)

The overall goal of this program is to develop and commercialize a turnkey quality control solution for the entire PEM fuel cell manufacturing process including membrane, gas diffusion layers, catalyst, and assembled systems. This quality control solution is unique to each customer’s specific needs but includes a suite of in-line quality control systems for roll-to-roll manufacturing that can target thin, transparent membrane as well as opaque membrane, catalyst, and GDLs. The Phase II developed the CPNUVV system for thin, transparent membrane that operates using polarized filters to enhance defect resolution and determine thickness. However, PEM material manufacturers want a complete, turnkey solution for all components of the PEM fuel cell. In the Phase IIB, Mainstream will partner with NREL, and transition reflectance technology NREL developed (US Patent 9,234,843) to operate in real-time on a web-line to develop a total solution to PEM quality control. The device will identify and mark defects as well as monitor materials thickness in real-time to improve line efficiency and to reduce waste. The research performed at NREL under this CRADA increases the basis of knowledge about optical inspection methods, how they perform with regards to fuel cell component materials, and how they can be implemented in in-line, real-time configurations to provide quality inspection for roll-to-roll (R2R) manufacturing. Methods explored were found to be sensitive to catalyst loading in electrodes on both membrane and gas diffusion media substrates and to membrane thickness, across a broad range of thickness. The latter capability is novel and now patented and provides a previously unstudied and undemonstrated capability for R2R manufacturing of membranes.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
DOE OSTI2021

Reduction of the molecular hamiltonian matrix using quantum community detection

Abstract Quantum chemistry is interested in calculating ground and excited states of molecular systems by solving the electronic Schrödinger equation. The exact numerical solution of this equation, frequently represented as an eigenvalue problem, remains unfeasible for most molecules and requires approximate methods. In this paper we introduce the use of Quantum Community Detection performed using the D-Wave quantum annealer to reduce the molecular Hamiltonian matrix in Slater determinant basis without chemical knowledge. Given a molecule represented by a matrix of Slater determinants, the connectivity between Slater determinants (as off-diagonal elements) is viewed as a graph adjacency matrix for determining multiple communities based on modularity maximization. A gauge metric based on perturbation theory is used to determine the lowest energy cluster. This cluster or sub-matrix of Slater determinants is used to calculate approximate ground state and excited state energies within chemical accuracy. The details of this method are described along with demonstrating its performance across multiple molecules of interest and bond dissociation cases. These examples provide proof-of-principle results for approximate solution of the electronic structure problem using quantum computing. This approach is general and shows potential to reduce the computational complexity of post-Hartree–Fock methods as future advances in quantum hardware become available.

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

Reconstructing thermal quantum quench dynamics from pure states

Simulating the nonequilibrium dynamics of thermal states is a fundamental problem across scales from high-energy to condensed-matter physics. Quantum computers may provide a way to solve this problem efficiently. Preparing a thermal state on a quantum computer is challenging, but there exist methods to circumvent this by computing a weighted sum of time-dependent matrix elements in a convenient basis. Further, while the number of basis states can be large, in this paper we show that it can be reduced by simulating only the largest density matrix elements by weight, capturing the density matrix to a specified precision. Leveraging Hamiltonian symmetries enables further reductions. This approach paves the way to more accurate thermal-state dynamics simulations on near-term quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. Here, we select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.

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

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
DOE OSTI2021

Model Reduction by Generalized Falk Method for Efficient Field-Circuit Simulations

The Generalized Falk Method (GFM) for coordinate transformation, together with two model-reduction strategies based on this method, are presented for efficient coupled field-circuit simulations. Each model-reduction strategy is based on a decision to retain specific linearly-independent vectors, called trial vectors, to construct a vector basis for coordinate transformation. The reduced-order models are guaranteed to be stable and passive since the GFM is a congruence transformation of originally symmetric positive definite systems. We also show that, unlike the Pad´e-via-Lanczos (PVL) method, the GFM does not generate unstable positive poles while reducing the order of circuit problems. Further, the proposed GFM is also faster when compared to methods of the type Lanczos (or Krylov) that are already widely used in circuit simulations for electrothermal and electromagnetic problems. The concept of response participation factors is introduced for the selection of the trial vectors in the proposed model-reduction methods. Further, we present methods to develop simple equivalent circuit networks for the field component of the overall field-circuit system. The implementation of these equivalent circuit networks in circuit simulators is discussed. With the proposed model-reduction strategies, significant improvement on the efficiency of the generalized Falk method is illustrated for coupled field-circuit problems.

42 ENGINEERING↗
DOE OSTI2023

A Systematic, Polynomial-Cost Approach to Exact Correlation Energies (Final Technical Report)

The full configuration interaction (FCI) wave function provides the exact solution to the Schrödinger equation in a given basis set. While FCI is intractable due to exponential computational costs, the many-body expansion (or the method of increments) can reduce scaling to a low-order polynomial with system size. This project entailed advances in the incremental FCI (iFCI) approach, designed to allow iFCI to reach larger system sizes and maintain its intrinsic high accuracy. This document describes advances in solvers for iFCI, strategies to treat multiple charge and spin states, and virtual state management methods to reduce memory requirements. Overall, this project allows iFCI to correlate (for the first time) 142 valence electrons in 444 orbitals in a realistic model of a transition metal complex.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

High-Fidelity, Low-Dissipation/Symmetry-Preserving Numerical Scheme for Solving the Euler Equations with Unstructured, Metric-Based Mesh Adaptation

This work presents an overview of a high-fidelity compressible Euler solver that utilizes the continuous Galerkin (CG) method with added artificial numerical diffusion for stabilization to solve a variety of unsteady and steady benchmark inviscid flow problems. This work shows that discretizing the Euler equations with this CG approach and first order basis functions produces a cost-effective stencil as well as simple well-posed boundary conditions. We show through convergence testing with manufactured solutions that the reduced stencil of CG, combined with the low amount of artificial diffusion required when using the stabilization method outlined in this work, leads to stable and highly accurate results for a variety of unsteady and steady applications. When combined with the adaptive mesh refinement approach used for many of the cases in this work, our results show that the flow solver achieves even more accurate results. A variety of inviscid flow cases are presented in this work, including transient 2D cases with complex shock structures and several steady 3D airfoils sections with a constant span.

Doetsch, Kevin [ORNL] (ORCID:0000000267051705)↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2024

An Evaluation of The Dynamic Physical Security Risk Assessment Methodology for Fleet-Wide Applications

The requirements for U.S. nuclear power plants to maintain a large onsite physical security force contribute to their high operational costs. The cost of maintaining the current physical security posture is approximately 10% of the overall operation and maintenance budget for commercial nuclear power plants. The goal of the Light Water Reactor Sustainability (LWRS) program’s physical security pathway is to develop tools, methods, and technologies and provide the technical basis for an optimized physical security posture. The conservatisms built into current security postures may be analyzed and minimized to reduce security costs while still ensuring adequate security and operational safety. The research performed at Idaho National Laboratory within LWRS program’s physical security pathway has successfully developed a dynamic force-on-force modeling framework using various computer simulation tools and integrating them with the dynamic assessment Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool. This integrated process for physical security analysis is named Modeling and Analysis for Safety Security using Dynamic EMRALD Framework (MASS-DEF). This document provides an update on the progress in applying the MASS-DEF process to an operating commercial nuclear power plant as well as additional industry feedback regarding use of the tool for other physical security risk-informed topics. This report is only a summary of the progress and does not contain specific modeling results as those contain sensitive security information. Previous reports described how a user could integrate their plant-specific force-on-force models with the dynamic simulation tool EMRALD, model operator actions, and integrate with probabilistic risk assessment tools, such as CAFTA (Computer Aided Fault Tree Analysis System) or SAPHIRE (Systems Analysis Programs for Hands-on Integrated Reliability Evaluations), and with thermal-hydraulic tools, such as RELAP-5 or MAAP. Previous reports applied various combinations of available simulations codes with EMRALD using generic plant models to demonstrate how to perform the analysis. This report is an update the progress of applying the dynamic computational framework to an actual nuclear facility using their security scenarios and timelines. This report also provides an update to the procedural guidance for the MASS-DEF process and an overview of the generic models available for use by utilities. This report does not contain any plant’s sensitive information and/or safeguards information. This study’s purpose was to verify that the results achieved using generic models are similar to actual plant results and refine our guidance on the use of the framework. This assessment enables further analysis, such as what-if scenarios and staff-reduction evaluation, thereby optimizing physical security at plants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Double-Wiebe Function for Reactivity Controlled Compression Ignition Combustion Using Reformate Diesel

Abstract Reactivity controlled compression ignition (RCCI) combustion has previously been proposed as a method to achieve high fuel conversion efficiency and reduce engine emissions. A single-fuel RCCI combustion strategy can have decreased fuel system complexity by using a reformate fuel for port fuel injection and the parent fuel (diesel) for direct injection. This paper presents a one-dimensional computational model of a compression ignition engine with single-fuel RCCI. A Wiebe function is used to predict the combustion process by representing the mass fraction burned (MFB) on a crank angle resolved basis. One single-Wiebe function (SWF) and two double-Wiebe functions(DWFs) were fitted to experimentally derive MFB data using the least-square method. The fitted results were compared with MFBs calculated from experimental data to verify the accuracy. The SWF did not fully capture the MFB curve with high fidelity while the detailed DWF captured the MFB curve within a root mean square error of 1.4%. The reduced double-Wiebe function (RDWF) also resulted in a predicted combustion profile with similar accuracy. Hence, the RDWF was used in a GT-power thermodynamic study to understand the effects of the low-temperature heat release (LTHR) fraction and combustion phasing on combustion characteristics. At optimum phasing of 5–10 crank angle degree after the top dead center, increasing the LTHR fraction from 20% to 60% resulted in the fuel conversion efficiency increasing from 39.5% to 41.1%, thus suggesting that the reformate fuel-based RCCI strategy is viable to unlock improved combustion performance.

Energy & Fuels↗
DOE OSTI2024

Capturing the Relative‐Humidity‐Sensitive Gas–Particle Partitioning of Organic Aerosols in a 2D Volatility Basis Set

Abstract Aerosol water affects the physicochemical properties and mass concentration of organic aerosols (OA), but it is typically omitted by air quality, weather, and climate models. We compare two classes of simplified models to estimate the OA water uptake and gas–particle partitioning of organic compounds. One class uses a single‐hygroscopicity‐parameter (κ) approach while the other is based on the reduced‐complexity Binary Activity Thermodynamics (BAT) model. We show that a BAT‐based two‐dimensional volatility basis set (VBS) model always predicts a higher OA mass concentration at elevated relative humidity (RH), for example, ∼16% at 80% RH, than any variation of theκ‐based method considered—even when BAT‐VBS predicts a lower water uptake. The main reason being that the BAT‐VBS model captures variations in effective saturation mass concentration of organics (C*) with RH, a feature that other VBS methods lack. The BAT‐VBS framework offers an efficient, RH‐sensitive treatment for reduced‐complexity OA modeling.

Geology↗
DOE OSTI2021

Geometry optimization speedup through a geodesic approach to internal coordinates

We present a new geodesic-based method for geometry optimization in a basis set of redundant internal coordinates. Overall, our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to converge to a minimum. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to numerically solve an ordinary differential equation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗