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

Nonlinear Matrix Approximation with Radial Basis Function Components

We introduce and investigate matrix approximation by decomposition into a sum of radial basis function (RBF) components. An RBF component is a generalization of the outer product between a pair of vectors, where an RBF function replaces the scalar multiplication between individual vector elements. Even though the RBF functions are positive definite, the summation across components is not restricted to convex combinations and allows us to compute the decomposition for any real matrix that is not necessarily symmetric or positive definite. We formulate the problem of seeking such a decomposition as an optimization problem with a nonlinear and non-convex loss function. Several modern versions of the gradient descent method, including their scalable stochastic counterparts, are used to solve this problem. We provide extensive empirical evidence of the effectiveness of the RBF decomposition and that of the gradient-based fitting algorithm. While being conceptually motivated by singular value decomposition (SVD), our proposed nonlinear counterpart outperforms SVD by drastically reducing the memory required to approximate a data matrix with the same L2 error for a wide range of matrix types. For example, it leads to 2 to 6 times memory save for Gaussian noise, graph adjacency matrices, and kernel matrices. Moreover, this proximity-based decomposition can offer additional interpretability in applications that involve, e.g., capturing the inner low-dimensional structure of the data, retaining graph connectivity structure, and preserving the acutance of images.

Rebrova, Elizaveta↗
DOE OSTI2020

Risk-Informed ATF and FLEX Analysis for an Enhanced Resilient BWR Under Design-Basis and Beyond-Design-Basis Accidents

This report documents the activities performed by the Idaho National Laboratory (INL) during fiscal year (FY) 2020 for the U.S. Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program, Risk-Informed System Analysis (RISA) Pathway, Enhanced Resilient Plant (ERP) Systems research. The purpose of the RISA Pathway research and development is to support plant owner-operator decisions with the aim to improve the economics, reliability, and maintain the high levels of safety of current nuclear power plants over periods of extended plant operations. The concept of ERP refers to the combinations of accident-tolerant fuel (ATF), optimal use of diverse and flexible coping strategy (FLEX), enhancements to plant components and systems, and the incorporation of augmented or new passive cooling systems, as well as improved fuel cycle efficiency. The objective of the ERP research effort is to use the RISA methods and toolkit in industry applications, including methods development and early demonstration of technologies, in order to enhance existing reactors’ safety features (both active and passive) and to substantially reduce operating costs through risk-informed approaches to plant design modifications to the plant and their characterization. One main focus of the FY 2020 efforts documented in this report was to extend the analyses conducted in FYs 2018 and 2019 for a pressurized water reactor (PWR) to a boiling water reactor (BWR). The same analysis process, risk analysis approaches, and analysis tools as in the previous work for PWR were used for a generic BWR with near-term ATF cladding (i.e., Iron-Chromium-Aluminum [FeCrAl] cladding and Chromium [Cr]-coated cladding) designs under the postulated station blackout (SBO) and medium loss-of-coolant (MLOCA) accident scenarios. In addition, a FLEX model was developed and incorporated into a generic BWR probabilistic risk assessment (PRA) model using the INL-developed software tool, Systems Analysis Programs for Hands-on Integrated Reliability Evaluations (SAPHIRE), to assess the risk impact from FLEX. The other main focus of the FY 2020 efforts was to advance analysis methods, including developing dynamic approach for FLEX human reliability analysis (HRA) using the INL-developed software tool, Event Modeling Risk Assessment using Linked Diagrams (EMRALD), as well as developing a multicriterion benefit evaluation (MCBE) method for evaluating costs and benefits of safety enhancements in nuclear power plants (NPPs). As a case study, the MCBE method was applied to evaluate the costs and benefits brought by FLEX implementation.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

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

A catalogue of white dwarfs in Gaia EDR3

ABSTRACT We present a catalogue of white dwarf candidates selected from Gaia Early Data Release 3 (EDR3). We applied several selection criteria in absolute magnitude, colour, and Gaia quality flags to remove objects with unreliable measurements while preserving most stars compatible with the white dwarf locus in the Hertzsprung–Russell diagram. We then used a sample of over 30 000 spectroscopically confirmed white dwarfs and contaminants from the Sloan Digital Sky Survey (SDSS) to map the distribution of these objects in the Gaia absolute magnitude–colour space. Finally, we adopt the same method presented in our previous work on Gaia Data Release 2 (DR2) to calculate a probability of being a white dwarf (PWD) for ≃1.3 million sources that passed our quality selection. The PWD values can be used to select a sample of ${\simeq} 359\,000$ high-confidence white dwarf candidates. We calculated stellar parameters (effective temperature, surface gravity, and mass) for all these stars by fitting Gaia astrometry and photometry with synthetic pure-H, pure-He, and mixed H–He atmospheric models. We estimate an upper limit of 93 per cent for the overall completeness of our catalogue for white dwarfs with G ≤ 20 mag and effective temperature (Teff) > 7000 K, at high Galactic latitudes (|b| > 20°). Alongside the main catalogue we include a reduced proper motion extension containing ${\simeq} 10\,200$ white dwarf candidates with unreliable parallax measurements that could, however, be identified on the basis of their proper motion. We also performed a cross-match of our catalogues with SDSS Data Release 16 (DR16) spectroscopy and provide spectral classification based on visual inspection for all resulting matches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High- temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes ad- vantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High-temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Parallelization and Performance Portability in Hydrodynamics Codes

With the eve of Exascale computing, performance and portability are at the forefront of all scientific codes. Adding more cores and more energy to a system is no longer a sustainable way to achieve performance, and extra effort must now be made to improve performance in all areas of code and code development. Using hydrodynamic codes as a basis, this work explores numerous techniques to achieve performance in different ways. Adaptive mesh refinement (AMR) is a necessary technique to improve memory optimization in mesh-based simulations. However it is invasive and conventionally difficult to integrate into existing applications, so we present a new branch of AMR to create a smooth transition to these optimizations, which not only improves performance, but also greatly reduces developer effort. We introduce the concept of this improvement as Phantom-Cell AMR, and assess theoretically the improvements, as well as present an application of its use. Other work included involves and investigation into an efficient data structure that ensures optimal memory layout for cache performance, with a target of making codes performant and portable across all architectures. All of the work targets both performance and portability, not just on CPU hardware, but specifically across GPU architectures. Parallel performance is key to all of the methods presented, but the research makes a great effort to improve the portability of all applications to prepare for current high performance computing systems and those on the horizon.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Drug discovery efforts at George Mason University

With over 39,000 students, and research expenditures in excess of $200 million, George Mason University (GMU) is the largest R1 (Carnegie Classification of very high research activity) university in Virginia. Mason scientists have been involved in the discovery and development of novel diagnostics and therapeutics in areas as diverse as infectious diseases and cancer. Below are highlights of the efforts being led by Mason researchers in the drug discovery arena. To enable targeted cellular delivery, and non-biomedical applications, Veneziano and colleagues have developed a synthesis strategy that enables the design of self-assembling DNA nanoparticles (DNA origami) with prescribed shape and size in the 10 to 100 nm range. The nanoparticles can be loaded with molecules of interest such as drugs, proteins and peptides, and are a promising new addition to the drug delivery platforms currently in use. The investigators also recently used the DNA origami nanoparticles to fine tune the spatial presentation of immunogens to study the impact on B cell activation. These studies are an important step towards the rational design of vaccines for a variety of infectious agents. To elucidate the parameters for optimizing the delivery efficiency of lipid nanoparticles (LNPs), Buschmann, Paige and colleagues have devised methods for predicting and experimentally validating the pKa of LNPs based on the structure of the ionizable lipids used to formulate the LNPs. These studies may pave the way for the development of new LNP delivery vehicles that have reduced systemic distribution and improved endosomal release of their cargo post administration. To better understand protein-protein interactions and identify potential drug targets that disrupt such interactions, Luchini and colleagues have developed a methodology that identifies contact points between proteins using small molecule dyes. The dye molecules noncovalently bind to the accessible surfaces of a protein complex with very high affinity, but are excluded from contact regions. When the complex is denatured and digested with trypsin, the exposed regions covered by the dye do not get cleaved by the enzyme, whereas the contact points are digested. The resulting fragments can then be identified using mass spectrometry. The data generated can serve as the basis for designing small molecules and peptides that can disrupt the formation of protein complexes involved in disease processes. For example, using peptides based on the interleukin 1 receptor accessory protein (IL-1RAcP), Luchini, Liotta, Paige and colleagues disrupted the formation of IL-1/IL-R/IL-1RAcP complex and demonstrated that the inhibition of complex formation reduced the inflammatory response to IL-1B. Working on the discovery of novel antimicrobial agents, Bishop, van Hoek and colleagues have discovered a number of antimicrobial peptides from reptiles and other species. DRGN-1, is a synthetic peptide based on a histone H1-derived peptide that they had identified from Komodo Dragon plasma. DRGN-1 was shown to disrupt bacterial biofilms and promote wound healing in an animal model. The peptide, along with others, is being developed and tested in preclinical studies. Other research by van Hoek and colleagues focuses on in silico antimicrobial peptide discovery, screening of small molecules for antibacterial properties, as well as assessment of diffusible signal factors (DFS) as future therapeutics. The above examples provide insight into the cutting-edge studies undertaken by GMU scientists to develop novel methodologies and platform technologies important to drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Tribal Revegetation Project Final Project Report: 92-Acre Area, Area 5 Radioactive Waste Management Complex, Nevada National Security Site, Nevada

Nuwu (Southern Paiute), Newe (Western Shoshone), and Nuumu (Owens Valley Paiute) are linguistically related, Numic-speaking peoples who are part of the broader Uto-Aztecan language group. Numic peoples view the land as a holistic, living, sentient being with feelings and purpose. The land is personified with human characteristics and it needs to be experienced to be understood through “learning by doing.” Numic peoples do not support ground disturbing activities within their ancestral lands, including activities tied to the storage of low-level radioactive waste or classified materials on the NNSS, which they view as culturally inappropriate. These deep-rooted ancestral connections are the impetus for reinforcing Numic responsibility for healing disturbed areas by integrating respect and patience with consistent Tribal interaction. Tribal Ecological Knowledge (TEK) is the science of describing Tribal approaches for understanding natural resources. Numic TEK is embedded in traditional teachings learned incrementally over time though experience and it evolves through lessons learned and responses to environmental changes over millennia. Therefore, TEK can broaden and enhance Western scientific knowledge associated with revegetation, especially in highly disturbed areas. The project blended TEK with Western scientific ecological methods to create a vegetative cover within test plots on the 92-Acre Area located at the Radioactive Waste Management Complex (RWMC) located in Area 5 on the Nevada National Security Site (NNSS). The vegetated test plots were systematically created for the Department of Energy (DOE) in tandem with the existing Federal Facilities Agreement and Consent Order (FFACO) with the Nevada Department of Environmental Protection (NDEP). Three previous contractor-lead attempts at revegetation, one targeting full cover revegetation and two targeting test plot revegetation, did not achieve the anticipated results at this location. When presented to the 16 American Indian Tribal nations and affiliated groups with cultural and historical ties to the NNSS, the group appointed a Tribal Revegetation Committee (TRC) that included six expert Tribal knowledge holders to collaborate with an ethnoecologist/cultural anthropologist and two biologists. Project design, planning, seed and outplant selection, spiritual land preparation, and methodology were guided by the TRC and an ethnoecologist/cultural anthropologist from Portland State University (PSU) and biologists from Desert Research Institute (DRI). Using TEK, the TRC recommended a specific seed mixture that contained nine native plant species and three species of outplants. The revegetation effort included preparing and planting thirty 10 m × 10 m (32.8 ft × 32.8 ft) seeded plots, twelve of which also included outplants; and eight 10 m × 100 m (32.8 ft × 328 ft) plots that only received outplants, all atop a waste cell cover. The TRC and the project team creatively adapted TEK with Western scientific methods so that the planned revegetation efforts could occur within the safety and security parameters of the RWMC. Test plots were subjected to one of five soil treatments with varying combinations of straw or mulch applications, soil amendment, and/or outplant planting and one of two watering treatments (watered or unwatered). Planting was divided into two events: one in the fall season and another during the subsequent spring season based on TEK and a corresponding recommendation from the TRC. Monitoring and spiritual management was conducted by the TRC to evaluate plant progress on a monthly basis (in conjunction with the biologist and anthropologist) during each respective growing season for a period of three years after planting. This approach allowed Tribal members the opportunity to conduct traditional blessings and other culturally appropriate activities to restore balance to the land in accordance with Tribal protocols. Following TEK-guided methods, successful plant establishment from seed stock and outplants was observed in many plots. Overall, plots planted in the spring, as recommended by the TRC, showed higher rates of outplant survival and native seedling emergence than those planted in the fall. This finding is significant because it is contrary to the original guidance and previous approaches provided for planting in this region. The TRC believes the frequent co-occurrence of native seedlings near surviving outplants indicates an important symbiotic relationship understood by Tribal communities and overlooked by others. Watered outplants displayed much higher survivability than unwatered plants, even after watering was reduced after the plants were established. Soil amendments and mulch created higher densities of native plants from seed. Many native seedlings showed significant delays in germination, which is considered a normal adaptation to desert climates. Some native plants did not germinate until the third growing year, whereas others germinated during the first growing year, which demonstrates the complexity of the desert environment. Evidence of native insects, reptiles, mammals, and birds, as well as native plants that were not part of the planted species, were noted and considered culturally significant. Despite the presence of non-native plants, native outplants continued to thrive and the incidence of native plant germination from seed increased over time. These successful revegetation results where previous efforts were unsuccessful reinforce the importance of integrating regionally appropriate, TEK-guided methodology with Western science to achieve positive results and the necessity of integrating Tribal involvement in all stages of the revegetation effort. Expanded approaches coupled with Tribal knowledge and tools from Western science addressed a complex problem tied to revegetating atop a low-level radioactive waste cell. The level of Tribal participation serves as a progressive model for building collaborative relationships and addressing ecological challenges on the NNSS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells (Final Report)

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co2O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly4 increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO 2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly 4 increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO 2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Removing Basis Set Incompleteness Error in Finite-Temperature Electronic Structure Calculations: Two-Electron Systems

We investigate the basis-set-size dependence for quantities related to interacting electrons in the canonical ensemble. Calculations are performed using exact diagonalization (finite temperature full configuration interaction method) on two-electron model systems–the uniform electron gas (UEG) and the helium atom. Our data reproduce previous observations of a competition for how the internal energy converges between the ground-state correlation energy and the high-temperature kinetic energy. We explore how this can be related to component parts of the internal energy including kinetic, exchange, and correlation energies and show there is surprising nuance in how this can be broken down into mostly monotonically converging quantities. We also show that separation of the free energy into a free energy with/without correlation allows for monotonic convergence with basis set size due to the variational principle. We find that the free energy convergence matches the previously observed convergence properties of the internal energy. We discuss the free energy divergence that happens when converging a finite basis analytical hydrogen atom to the complete basis set limit and compare this to the energies of a helium atom in a large periodic box. Reducing the box size, we saw convergence trends for the helium atom that were similar to the UEG.

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

Machine-learning based model reduction for partial differential equations

We develop a novel synergistic approach between model reduction and machine learning. The specific goal of this project is to aid in the construction of reduced order models for basis functions that are custom-made to represent the solution of partial differential equations. Partial differential equations (PDEs) are one of the main mathematical tools for describing physical phenomena. However, due to either efficiency or necessity, for many real-world problems, we are interested in constructing reduced order models (ROMs) which focus only on the explicit computation of subsets of the active spatio-temporal scales in the problem, while treating the interaction with the rest of the scales approximately. The task of accurate representation of such interactions (usually called memory terms) constitutes a vast area of research known as model reduction. PI Stinis has significant expertise in the construction of ROMs for complex systems. In addition, in recent work with the project key participant Qadeer, they have utilized machine learning to acquire custom-made basis functions (CBFs) to expand the solutions of PDEs. In the proposed work, we will merge the two concepts by constructing ROMs for subsets of the CBFs needed to represent the solution of a PDE. Specifically, we will use the Mori-Zwanzig model reduction formalism to construct ROMs for subsets of CBFs for nonlinear PDEs of various complexity, as well as investigate the usage of CBFs in the spectral vanishing viscosity method for problems that can form shocks in finite time. The outcome of the research is aimed to be proof-of-concept about a novel synergistic approach between model reduction and machine learning, thus advancing the field of scientific machine learning. Such a capability will benefit the efficient modeling of physical systems appearing in various areas of interest to the DOE.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗