Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model order reduction”

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 199 records · Page 11

Practical Considerations for Understanding Surface Reaction Mechanisms Involved in Heterogeneous Catalysis

Acquiring useful knowledge about the active site(s) of a catalyst, nature of reactant–catalyst interactions, nature of reactive intermediates, rate-determining step, reaction rate orders that affect various process parameters, and reaction mechanism as a whole is exceedingly challenging. This is especially true in the case of heterogeneous catalysts due to the complexity of the nature of surface active sites and their nonstatic behavior. Here, we present our perspective on differentiating between various surface reaction mechanisms in light of pioneering studies by leaders in the field, with the aim of clarifying some of the confusion associated with these complex mechanisms, especially the Eley–Rideal mechanism. Using bibliometric analysis, we identify and discuss the following four reactions that most commonly invoke the Eley– Rideal mechanism: H 2 activation, CO oxidation, esterification of alcohols by acids, and selective catalytic reduction (SCR) of NO x with NH 3 . Our analysis of studies utilizing well-suited experimental and computational methodologies for differentiating surface reaction mechanisms suggests that the above-mentioned four reactions do not occur via the Eley–Rideal mechanism. Instead, each reaction occurs via the Langmuir–Hinshelwood mechanism with nonidealities present. Lastly, we highlight practical considerations regarding select experimental (characterization methods and differential kinetics) and computational modeling that we believe can provide useful insights to accurately discern between the various possible reaction mechanisms in heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Solar Fuels Nexus: Molecules and Materials for Light-Driven Catalysis

The American Chemical Society (ACS) selects two groups of graduate students each year to plan and host a one-day symposium at each national meeting (both fall and spring).This year our Graduate Student Symposium Planning Committee (GSSPC), composed of seven students from four universities, proposal entitled “A Solar Fuels Nexus: Molecules and Materials for Light-Driven Catalysis” was selected for the “Crossroads in Chemistry” ACS Meeting that will take place March 23-26, 2023 in Indianapolis, IN. All members of our GSSPC are affiliated with the DOE Fuels from Sunlight Energy Innovation Hub, with two from the Liquid Sunlight Alliance (LiSA) and five from the Center for Hybrid Approaches in Solar Energy to Liquid Fuels (CHASE). Here we request funds to support this symposium. This symposium will highlight research progress and perspectives in the solar fuels generation field and seeks to advance the four priority research objectives (PROs) established by the Department of Energy’s Office of Basic Energy Sciences (DOE-BES) Roundtable Report that are also central to many research goals within LiSA and CHASE. The symposium will consist of research presentations from 10 invited senior researcher speakers on topics such as molecular catalyst design, computational modeling of electron transfer systems, microenvironmental effects on CO2 reduction and H2O oxidation catalysis, and intelligent design of semiconductor interfaces with ample time for discussions. These research topics fit very well with the Solar Photochemistry supported research areas of “light-driven electron and energy transfer in condensed phase and interfacial molecular systems,” “electrocatalysis and photocatalysis of solar fuels reactions,” and “semiconductor photoelectrochemistry.” More broadly, this symposium seeks to advance the DOE-BES’s mission to: “support fundamental research to understand, predict, and ultimately control matter and energy at the level of electrons, atoms, and molecules” by providing a diverse atmosphere where such research will be disseminated, discussed, and debated. There will be a strong focus on Diversity, Equity, and Inclusivity (DEI) in our symposium. Of our 10 speakers, 7 will be from underrepresented demographics in STEM, including 5 who identify as women. Furthermore, we have representatives from academia accompanied by one national lab scientist and one officer from the Office of Fossil Energy and Carbon Management at the DOE. All speakers will be holding a short DEI moment ahead of their talks. In order to support the career development of attending early career scientists, we will also be hosting a luncheon specifically for graduate students and postdocs to provide them opportunities to network with the distinguished speakers and other attendees. DOE funds for this symposium will be used to support the attendance and participation of 15 graduate students from US institutions by defraying travel and registration costs. These funds will promote engagement and conversation between early career scientists in the solar fuels field, while disseminating solar fuels research funded by and relevant to the DOE.

30 DIRECT ENERGY CONVERSION↗

Retrieving Top-k Hyperedge Triplets: Models and Applications

Complex systems frequently exhibit multi-way, rather than pairwise, interactions. These group interactions can- not be faithfully modeled as collections of pairwise interactions using graphs and instead require hypergraphs. However, methods that analyze hypergraphs directly, rather than via lossy graph reductions, remain limited. Hypergraph motifs hold promise in this regard, as motif patterns serve as building blocks for larger group interactions which are inexpressible by graphs. Recent work has focused on categorizing and counting hypergraph motifs based on the existence of nodes in hyperedge intersection regions. Here, we argue that the relative sizes of hyperedge inter- sections within motifs contain varied and valuable information. We propose a suite of efficient algorithms for finding top-k triplets of hyperedges based on optimizing the sizes of these intersection patterns. This formulation uncovers interesting local patterns of interaction, finding hyperedge triplets that either (1) are the least similar with each other, (2) have the highest pairwise but not groupwise correlation, or (3) are the most similar with each other. We formalize this as a combinatorial optimization problem and design efficient algorithms based on filtering hyperedges. Our comprehensive experimental evaluation shows that the resulting hyperedge triplets yield insightful information on real-world hypergraphs. Our approach is also orders of magnitude faster than a naive baseline implementation.

hypergraphs, motifs, Combinatorial Algorithms↗

Electrochemical Oxidation in Garnet-Type Solid Electrolyte by Formation of Point Defects

All-solid-state batteries hold greater promise for improving safety and energy density over conventional battery technology employing organic liquid electrolytes. One of the required features of a Li + conducting solid electrolyte is electrochemical stability, attained thermodynamically or kinetically, within the targeted operating voltage and temperature ranges. Therefore, understanding of the oxidative or reductive degradation mechanism is important to allow the design of stable solid electrolyte materials. This work contributes to building an understanding of the oxidative degradation mechanism in lithium solid electrolytes at cell operating conditions. Here, we have focused on resolving the oxidative decomposition mechanism of Al-doped lithium garnet Li 6.28 Al 0.24 La 3 Zr 2 O 12 (LLZO) as a state-of-the-art inorganic ceramic electrolyte. By combining experimental and computational analyses, we show that oxidation of LLZO occurs by simultaneous loss of oxygen and lithium from the structure, resulting in substoichiometric LLZO, at a moderate temperature (80 °C) and a high electrode potential (4.3 V vs Li/Li + ). Based on X-ray absorption and diffraction analyses, we find that the zirconium coordination shells in LLZO contract while the crystal structure experiences positive chemical strain upon electrochemical oxidation. The results from ex situ structural characterization of both the local structure and crystal symmetry are supported by a substoichiometric LLZO with lithium and oxygen vacancies, modeled by density functional theory (DFT) calculations. These chemical and structural changes in LLZO suppress effective lithium-ion conductivity by an order of magnitude. Formation of lithium and oxygen vacancies in LLZO upon electrochemical oxidation is different from prior thermodynamic predictions of phase decomposition of LLZO. The difference here is that the experiments were conducted at near-room temperature, which can hinder the kinetics of phase separation, and thus, the resultant LLZO solid electrolyte is still single-phase but substoichiometric in Li and O. In conclusion, these findings contribute an important degradation mechanism of the electrolyte, relevant for practical operational conditions of solid-state batteries.

36 MATERIALS SCIENCE↗

Spectroscopy and crystallography define carotenoid oxygenases as a new subclass of mononuclear non-heme Fe II enzymes

Carotenoid cleavage dioxygenases (CCDs) are non-heme Fe II enzymes that catalyze the oxidative cleavage of alkene bonds in carotenoids, stilbenoids, and related compounds. How these enzymes control the reaction of dioxygen (O 2 ) with their alkene substrates is unclear. Here, we apply spectroscopy in conjunction with X-ray crystallography to define the iron coordination geometry of a model CCD, CAO1 (Neurospora crassa carotenoid oxygenase 1), in its resting state and following substrate binding and coordination sphere substitutions. Resting CAO1 exhibits a five-coordinate (5C), square pyramidal Fe II center that undergoes steric distortion toward a trigonal bipyramidal geometry in the presence of piceatannol. Titrations with the O 2 -analog, nitric oxide, show a >100-fold increase in iron–nitric oxide affinity upon substrate binding, defining a crucial role for the substrate in activating the Fe II site for O 2 reactivity. The importance of the 5C Fe II structure for reactivity was probed through mutagenesis of the second-sphere Thr151 residue of CAO1, which occludes ligand binding at the sixth coordination position. A T151G substitution resulted in the conversion of the iron center to a six-coordinate state and a 135-fold reduction in apparent catalytic efficiency toward piceatannol compared with the wildtype enzyme. Substrate complexation resulted in partial six-coordinate to 5C conversion, indicating solvent dissociation from the iron center. Additional substitutions at this site demonstrated a general functional importance of the occluding residue within the CCD superfamily. Taken together, these data suggest an ordered mechanism of CCD catalysis occurring via substrate-promoted solvent replacement by O 2 . CCDs thus represent a new class of mononuclear non-heme Fe II enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon Tetrachloride Degradation Results for 200-ZP-1 Operable Unit

Carbon tetrachloride (CT) contamination in the 200-ZP-1 Operable Unit (OU) at the Hanford Site originated from large-volume discharges to the subsurface during plutonium production operations between 1955 and 1973. Contamination migrated through more than 70 meters of unsaturated sediment to reach the underlying unconfined aquifer, where it persists as a large and complex groundwater plume. The 200-ZP-1 OU Record of Decision (ROD) requires that groundwater CT concentrations be reduced to 3.4 µg/L within 125 years. Current groundwater modeling projections estimate that the existing pump-and-treat, even when combined with monitored natural attenuation (specifically hydrolysis), will not achieve this target within the designated timeframe. A fundamental contributor to this shortfall is the extremely slow rate of CT hydrolysis under Hanford aquifer conditions, which has been estimated to have a half-life of 630 years. If faster-acting biotic and abiotic degradation processes are operating within the aquifer, their contribution to CT mass reduction could have a meaningful impact. However, site-specific measurements of these processes and their rates have not previously been performed. This report documents the results of a two-phase laboratory investigation designed to characterize and quantify the capacity of site-specific 200-ZP-1 OU sediments and groundwater to support natural attenuation of CT through biotic and abiotic pathways. In this context, degradation capacity is defined as the intrinsic potential of the subsurface matrix to transform CT under optimized, controlled conditions. System capacity is evaluated in two ways: (1) as rate-limited capacity, which establishes the maximum kinetic velocity of CT transformation and is measured using half-lives and first order rate constants; and (2) as mass limited capacity, which defines the total contaminant mass the batch experimental system can degrade before reactants are exhausted, representing the maximum amount of contaminant the microbial community and reactive mineral phases can transform under the experimental conditions.

abiotic degradation↗

Symmetry instability induced by topological phase transitions

The symmetry-topology interplay dictates how to define order parameters and classify material ordered phases. However, current understanding of this interplay has been predominately approached from a one-sided perspective, with topological states being classified within the constraints imposed by specific fixed symmetries. Here, in this study, we complete this full circle by demonstrating spontaneous symmetry breaking that results from a periodic alteration of topological phases induced by light in a centrosymmetric Dirac material ZrTe 5 . The distinguishing feature is the observation of robust correlation and striking anomalies in the fluence and temperature dependence of key transport parameters. First, both shift current 𝐽 𝑠 and displacement current 𝐽 𝑑 , arising from interband transition and infrared phonon driving, respectively, along with charge carrier pumping, exhibit similar behaviors. Second, they all peak at similar low pump fluence, followed by a subsequent reduction as the fluence further increases. This behavior cannot be explained by conventional energetically allowed, direct excitations. Third, all the three observables exhibit anomalies when they approach the topological phase transition temperature. These results highlight the unique low-energy pumping behaviors in ZrTe 5 , characterized by reversible fluence dependence and a “hinge-like” interaction that connects various electronic and lattice observables, including phonons, charge carriers, and currents. Our findings, supported by model analysis, provide key insights into the fragility of crystalline (inversion) and time-reversal symmetries during the dynamics of topological phase transitions. This fragility drives spontaneous symmetry breaking, evidenced by the synchronized emergence of off-resonant infrared phonons and broken-symmetry photocurrents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)

Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.

36 MATERIALS SCIENCE↗

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗

Quantifying Capital Cost Reduction Pathways for Advanced Nuclear Reactors

The framework developed in this study is provided both as an excel sheet (https://inl.gov/content/uploads/2023/11/Nuclear-Reactor-Cost-Reduction-Pathway-Spreadsheet-Tool.xlsx) and a Python (Jupyter) Notebook (link: https://github.com/accert-dev/ACCERT/tree/main/Cost%20Reduction). Capital cost considerations are one of the primary inhibitors to the large-scale deployment of nuclear power plants. While it is widely accepted that first units will likely be expensive and relatively uncompetitive, it is reasonable to expect that subsequent units, built in relative quick succession, will be cheaper as they benefit from the so-called “learning effects”. However, the large degree of uncertainty associated with this parameter renders it challenging for first movers to invest in the first few expensive units. To resolve this impasse, the U.S. Department of Energy’s Advanced Nuclear Liftoff study advocated for the formation of large, committed order books of plants of the same technology to spread the costs across several units and kickstart the nuclear supply chain. The study also advocated best practices for avoiding overruns and keeping reactors on budget. This report builds on these key recommendations by attempting to quantify specific pathways toward cost reduction for nuclear energy. A capital cost estimation framework was built to untangle the effect of learning into a subset of key cost drivers, referred to as “levers”. Collectively, the choice of these levers is intended to reflect the decision-making of high-level stakeholders like plant owners and the government. In addition to the size of the firm orderbook, these levers included (a) cost drivers that are most often attributed to cost overruns such as architect/engineering (A/E) proficiency, construction proficiency, procurement service proficiency, design completion prior to the start of construction, and design maturity, and (b) cost reduction strategies such as modular construction, cross-site standardization, safety classification of the reactor building, and of the balance of plant. Two advanced reactor designs were leveraged as use cases and bottom-up cost estimates made with assumptions consistent with a well-executed first-of-a-kind project (WE-FOAK, i.e., almost no overruns) were used as baselines for the models. Cost correlations were surveyed from the literature to determine the impact of important variables on projected timelines and costs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Packages of Distributed Energy Technologies Demonstrating Demand Flexibility at Community Scale

The combination of increased electric load growth across all sectors, deferred electrical infrastructure investment, and other factors resulting in variable electric power supply, has created technical challenges to maintaining a resilient and reliable grid. Many federal, regional, and local efforts are in play to modernize the electric grid, including advancing building technologies and distributed energy resources (DERs) that are utilizing smarter controls to become responsive to both occupant and grid needs. This report reviews ten pilot projects demonstrating how groups of buildings combined with behind-the-meter (BTM) DERs such as electric vehicle (EV) charging, battery storage, flexible HVAC and domestic hot water systems, and photovoltaic systems can reliably and cost effectively provide grid services. Each of the ten pilot projects aim to deliver both energy efficiency and demand flexibility (DF) while supporting load growth. The ten demonstration teams are piloting flexible DER packages across diverse communities of residential and commercial buildings to address a variety of regional grid needs. The outcomes of these pilot projects will be used to inform future scaling through utility program development. This paper characterizes the ten teams, showcasing the decision-making process used by each group to develop their packages (Section 2), the grid services they plan to deliver (Section 3), the types of DER packages selected for deployment within building sectors (Section 4) and trends between building sector, DER types, and grid services In order to achieve community scale benefits, the pilot projects must utilize aggregated control mechanisms for coordinating buildings and DERs together. Several types of coordinated control architectures have evolved amongst the teams, influenced by use type, existing market conditions, and integration type. Three coordinated controls architectures have been characterized, highlighting their use cases, benefits, challenges, and tradeoffs in their design. These insights can aid utilities, control vendors, and developers in scaling community-level energy systems (Paul, 2024). Ultimately, the technology packages selected by the ten teams will be coordinated to provide power system services, also known as grid services. Insights from these demonstrations will be useful for grid operators, regulators, aggregators and other stakeholders as they look to deploy demand flexible resources as grid services in the future. The grid services that each team is targeting for demonstration are described in Section 3 and Section 4. Methods for evaluating the grid services have been described in the paper Metrics for Evaluating Grid Service Provision from Communities of Grid-interactive and Efficient Buildings and other DER (MacDonald, 2023). To identify technology packages for demonstration, Section 2 shows that project teams used a range of analysis approaches, including building energy modeling, AMI data analysis, cost-benefit frameworks, and utility pilot data. Some teams emphasized technical modeling to quantify grid impacts and demand reduction potential, while others prioritized economic evaluations, stakeholder input, or exploratory pilots to inform deployment decisions. This diversity reflects the need to tailor selection methods to project goals, available data, and organizational context. Section 5 discusses trends between the DER technologies deployed and the grid service provisions from each team. Residential buildings (multifamily and single family) lean towards technologies that enhance energy efficiency (e.g. weatherization upgrades, smart thermostats) and onsite power generation integration (e.g. solar PV). Commercial building demonstrations prioritize technologies that ensure operational reliability (e.g. battery storage) and centralized energy management systems and optimization solutions. Teams that are deploying controllable storage-based technologies are more likely to provide grid services that require a near real-time response. Teams incorporating load shifting technologies like smart thermostats with HEMs are likely to include energy markets participation and customer bill management offerings. Campus demonstrations are adopting diverse sets of DERs to emphasize renewable generation, paired with centralized control. This section also describes technologies that were considered during project planning but ultimately excluded from final deployment. These demonstrations reveal that effective DER package design should be tailored to building type, customer segment, and construction vintage. Multifamily buildings benefit from centralized HVAC upgrades and supervisory controls, while single-family homes are well-suited for individualized technologies like solar, storage, and smart home energy monitors. Commercial and campus settings prioritize EMIS integration and load optimization. New construction enables cost-effective integration of DER-ready infrastructure, whereas retrofits require deployments aligned with owner and tenant value streams. For utility program planners, early coordination with developers and building owners, paired with segmented and modular program offerings, can improve adoption, scalability, and grid impact.

24 POWER TRANSMISSION AND DISTRIBUTION↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

ACRRF High-Bay Dose Calculations using MCNP (Part A)

Analytical tools and models have been developed as a starting point for directly assessing dose in the Annular Core Research Reactor Facility (ACRRF) due to reactor operation. Key results include peak dose along the Central Cavity (CC) centerline (beamline) at the cavity level, dose throughout the High-Bay (HB), and dose on the facility roof for partially-shielded reactor operation where the 4” insert is removed from the CC Shield Plug (SP). Model results in the beamline are benchmarked against measured doses from passive dosimetry evaluations. Personnel total (neutron and gamma) dose in the ACRRF HB is calculated using Monte Carlo N-Particle (MCNP). Various CC and SP configurations are analyzed, including unshielded (no SP) and partially shielded (SP installed but 4” insert removed). Novel application of Variance Reduction (VR) techniques, namely the Surface Source Write (SSW) and Surface Source Read (SSR) capabilities in MCNP, enable impressive resolution (in a Monte Carlo modeling sense) of dose throughout much the facility. The VR techniques reduce stochastic error for challenging tallies, with more advanced techniques explored in the companion to this report (Part B) [1]. Supplementary studies (including a verification analysis) and pedagogic evaluations in Part B involve neutron spectra, angular distributions, and the dose impact of facility characteristics. With the SP 4” insert removed and the Lead-Boron (44”) Bucket (LB–44) in the reactor cavity, Total Effective Dose (TED) within the CC beamline is ≈140 rem per 300 MJ of reactor yield (or 3900 rem per hour at 100% Steady-State (SS) power). With no SP (unshielded) and a Free-Field (FF) cavity, TED within the beamline is ≈610 rem per 300 MJ (or 17000 rem per hour at 100% SS power). Due to the predicted collimation of radiation by the reactor pool (and partial SP, if present), beamline dose is much greater than the scattered radiation field surrounding the cavity and reactor tank. Comparisons are made to beamline dosimetry measurements to validate the model. Model predictions agree reasonably well (⪅10%) with measured quantities of neutron fluence, gamma fluence, and spectral metrics. Away from the beamline, comparisons made to previous dose measurements in the HB agreement within an order of magnitude.

61 RADIATION PROTECTION AND DOSIMETRY↗

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗

Managing Subsurface Pressure Buildup and Interference in Commercial-Scale CO 2 Storage Project with Proximal Injection Wells

Large-scale decarbonization using carbon capture and storage (CCS) is likely to involve many commercial-scale CO 2 storage projects located in close proximity to each other. This close proximity raises concerns over pressure interference among the storage projects. Pressure interference between injection and storage efforts can reduce the practicable CO 2 storage resource and force wells to inject CO 2 at a lower rate to avoid the fracture pressure thresholds per United States Environmental Protection Agency (EPA) Class VI well regulations to preserve injection and confining zone integrity and potentially mitigate against inducing seismic activity. These analyses employ numerical full-physics reservoir modeling to evaluate how pressure buildup fronts and CO 2 plumes evolve under commercial-scale injection volumes of CO 2 in which multiple storage sites located in close proximity occur in tandem. The simulation models mimic injection at pseudo basin-scale and assume homogeneous saline formation(s) as storage targets with a pair of upper and lower homogeneous seal layer/s. These analyses specifically investigate the efficacy of two basin-wide reservoir pressure management strategies in addressing the technical challenges associated with pressure buildup and CO 2 plume commingling. The strategies explored include: 1) enlarging the area of injection well spacing (WS) and 2) storing CO 2 in a stacked sequence (SSS) of saline formations compared to a single formation. The storage and confining zones properties assumed were held common across the scenarios, unless specified otherwise. Analyses results show that after injecting 4 million tons per year for 30 years using 4 separate wells (each injecting 1 million metric tons per year), the radius of CO 2 plume extends to a mere 3 km or less from injection wells. Meanwhile, the radius of pressure buildup ranges on the order of tens to a few hundreds of kilometers, depending on the magnitude of pressure buildup threshold that one would use to define the front. CO 2 plume commingling from different injection wells appears to occur 50 years post-injection, especially under scenarios with narrowly spaced (i.e., < 5 km apart) injection well locations. Findings from sensitivity cases on the well spacing suggest that storage formations modeled would require different well spacing to avoid fracture pressure thresholds. For instance, modeled storage formations with high fracture gradients (i.e., 0.8 psi/ft) would need less than 5–km well spacing, whereas those with lower fracture gradients (i.e., 0.7 psi/ft) would need approximately 20–km well spacing, based on assumed modeling parameters. Under stacked injection, the pressure challenges (described above) still exist but are more alleviated due to distributing the same injection volume across more available reservoir volume. These analyses demonstrate that stacked-sequence storage can effectively address the challenges, while still providing the same target CO 2 storage volumes and allowing a large number of storage projects to be deployed in the same basin by better utilizing the available storage resource across different reservoir depths. Among cases modeled, the resulting pressure buildup front is most suppressed when each storage project distributes injection volumes over several wells, each of which injects a portion of the total CO 2 across the stacked sequence. This strategy results in the smallest CO 2 aerial footprint amongst scenarios evaluated but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 42 percent compared to the commercial-scale single-formation storage), the result of which is crucial to maintain caprock integrity. The findings presented by this research draw attention to the importance of greater coordination among storage operators and regulatory stakeholders to foster the upscaling and deployment of CCS. These analyses provide insights into required decision-making when considering multi-project deployment in a shared basin. Because these analyses evaluate a very specific geologic situation, they bear further investigations across other geologic situations.

42 ENGINEERING↗

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗

Comparison of Sub-Ppm Instrument Response Suggests Higher Detection Limits Could Be Used to Quantify Methane Emissions from Oil and Gas Infrastructure

Quantifying and controlling fugitive methane emissions from oil and gas facilities remains essential for addressing climate goals, but the costs associated with monitoring millions of production sites remain prohibitively expensive. Current thinking, supported by measurement and simple dispersion modelling, assumes single-digit parts-per-million instrumentation is required. To investigate instrument response, the inlets of three trace-methane (sub-ppm) analyzers were collocated on a facility designed to release gas of known composition at known flow rates between 0.4 and 5.2 kg CH 4 h –1 from simulated oil and gas infrastructure. Methane mixing ratios were measured by each instrument at 1 Hertz resolution over nine hours. While mixing ratios reported by a cavity ring-down spectrometer (CRDS)-based instrument were on average 10.0 ppm (range 1.8 to 83 ppm), a mid-infrared laser absorption spectroscopy (MIRA)-based instrument reported short-lived mixing ratios far larger than expected (range 1.8 to 779 ppm) with a similar nine-hour average to the CRDS (10.1 ppm). We suggest the peaks detected by the MIRA are likely caused by a micrometeorological phenomenon, where vortex shedding has resulted in heterogeneous methane plumes which only the MIRA can observe. Further analysis suggests an instrument like the MIRA (an optical-cavity-based instrument with cavity size ≤10 cm3 measuring at ≥2 Hz with air flow rates in the order of ≤0.3 slpm at distances of ≤20 m from the source) but with a higher detection limit (25 ppm) could detect enough of the high-concentration events to generate representative 20 min-average methane mixing ratios. Even though development of a lower-cost, high-precision, high-accuracy instrument with a 25 ppm detection threshold remains a significant problem, this has implications for the use of instrumentation with higher detection thresholds, resulting in the reduction in cost to measure methane emissions and providing a mechanism for the widespread deployment of effective leak detection and repair programs for all oil and gas infrastructure.

03 NATURAL GAS↗

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗