Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Design pattern”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Locking water molecules via ternary O–H⋯O intramolecular hydrogen bonds in perhydroxylated closo -dodecaborate

A multitude of applications related to perhydroxylated closo-dodecaborate B 12 (OH) 12 2− in the condensed phase are inseparable from the fundamental mechanisms underlying the high water orientation selectivity based on the base B 12 (OH) 12 2− . Herein, we directly compare the structural evolution of water clusters, ranging from monomer to hexamer, oriented by functional groups in the bases B 12 H 12 2− , B 12 H 11 OH 2− and B 12 (OH) 12 2− using multiple theoretical methods. A significant revelation is made regarding B 12 (OH) 12 2− : each additional water molecule is locked into the intramolecular hydrogen bond B–O–H ternary ring in an embedded form. This new pattern of water cluster growth suggests that B–(H–O)⋯H–O interactions prevail over the competition from water–hydrogen bonds (O⋯H–O), distinguishing it from the behavior observed in B 12 H 12 2− and B 12 H 11 OH 2− bases, in which competition arises from a mixed competing model involving dihydrogen bonds (B–H⋯H–O), conventional hydrogen bonds (B–(H–O)⋯H–O) and water hydrogen bonds (O⋯H–O). Through aqueous solvation and ab initio molecular dynamics analysis, we further demonstrate the largest water clusters in the first hydrated shell with exceptional thermodynamic stability around B 12 (OH) 12 2− . These findings provide a solid scientific foundation for the design of boron cluster chemistry incorporating hydroxyl-group-modified borate salts with potential implications for various applications.

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

FORCE-DISPATCHES Integration - Initial Demonstration

Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to improve economic viability under uncertain market and weather conditions. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Continuous efforts and investments from the IES programs have been made to expand and improve the versatility of the FORCE toolset in fiscal year 2022. Code-coupling and cross-tool communication have been important methods for improving this versatility. This report focuses on an additional workflow in the HERON tool for capacity and dispatch stochastic optimization through integration with the external tool Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES). DISPATCHES was primarily developed by the National Energy Technology Laboratory, in collaboration with other national laboratories, which included INL, universities, and industry partners. It is coupled to a library of algebraic models for specific plant components, and to a framework for stochastic optimization different from that provided in the current Risk Analysis Virtual Environment (RAVEN)-running-RAVEN algorithm in HERON. HERON currently conducts stochastic optimization via an outer-inner loop: it optimizes over variable capacity on the outer loop, and at each step within the capacity parameter space, conducts an inner optimization over scenarios (of market signals, demand, and/or weather patterns) and hourly dispatch throughout a user-specified number of years. On the other hand, DISPATCHES conducts stochastic optimization via an “all-at-once” strategy in which capacity variables are optimized at the same level as dispatch variables, as all scenarios are considered at once. The latter method works especially well for projects of limited size and project length, as the necessary computational power and memory increases with the number of variables and scenarios. The new capability to use the DISPATCHES workflow in HERON enhances standalone simulations by leveraging FORCE tools—namely, the economic metrics from the Tool for Economic Analysis (TEAL) and reduced-order model (ROM) sampling from RAVEN. The initial demonstration of the DISPATCHES workflow simulates an existing nuclear-case flowsheet within the DISPATCHES repository—this models a NPP with a secondary revenue stream for hydrogen production. Electrical output from the plant is converted to hydrogen via a proton-exchange membrane (PEM) electrolyzer, hydrogen tanks are used for storage, and an additional turbine is added for hydrogen combustion. Continued work regarding this FORCE-DISPATCHES integration will include automatic generation of DISPATCHES models from HERON inputs, offering analysts the option of using either the RAVEN-runsRAVEN or DISPATCHES workflow to solve technoeconomic optimization problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Hardware-in-the-loop Laboratory Performance Verification of Flexible Building Equipment in a Typical Commercial Building

This project aims to develop high-resolution equipment performance and occupant data that quantifies demand flexibility in typical commercial buildings. The dataset documented in this report includes comprehensive time-series measurements from hardware-in-the-loop (HIL) experiments conducted across multiple testbeds designed to simulate realistic operational environments for HVAC systems. Specifically, it captures minute-by-minute high-resolution data on the performance of various typical HVAC systems, including a variable-air-volume (VAV) air handling unit (AHU) system with chillers and an ice tank in the Intelligent Building Agents Laboratory (IBAL) at the National Institute of Standards and Technology (NIST), a two-stage air-source heat pump (ASHP) at the NIST, and a water-source heat pump (WSHP) at Texas A&M University (TAMU). The data were generated under a range of controlled conditions reflecting different grid scenarios and climatic influences, as well as various control strategies, building types, occupancy patterns, and occupant behaviors. This dataset provides DE-EE0009153 Final Report 5 detailed insights into the demand flexibility of these systems, including their response to grid signals, occupant behaviors, energy consumption patterns, and operational efficiency under different conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Climate impacts of particle physics

The pursuit of particle physics requires a stable and prosperous society. Today, our society is increasingly threatened by global climate change. Human-influenced climate change has already impacted weather patterns, and global warming will only increase unless deep reductions in emissions of CO$_2$ and other greenhouse gases are achieved. Current and future activities in particle physics need to be considered in this context, either on the moral ground that we have a responsibility to leave a habitable planet to future generations, or on the more practical ground that, because of their scale, particle physics projects and activities will be under scrutiny for their impact on the climate. In this white paper for the U.S. Particle Physics Community Planning Exercise ("Snowmass"), we examine several contexts in which the practice of particle physics has impacts on the climate. These include the construction of facilities, the design and operation of particle detectors, the use of large-scale computing, and the research activities of scientists. We offer recommendations on establishing climate-aware practices in particle physics, with the goal of reducing our impact on the climate. We invite members of the community to show their support for a sustainable particle physics field (https://indico.fnal.gov/event/53795/).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

FAIRification, Quality Assessment, and Missingness Pattern Discovery for Spatiotemporal Photovoltaic Data

The ongoing growth of the photovoltaic market has pushed the demand for power forecasting and performance evaluation for a huge population of PV power plants. Through access to a large number of time series data sets from different power plants, we have found common issues that impede the modeling process. Namely, the time series data are hard to transfer between groups due to differences in variable nomenclature, and the quality of the data sets can vary. We address the issue of variable nomenclature by FAIRifying spatiotemporal PV time series data. Through the creation of a solar power plant ontology, we propose standards for the naming and structure of metadata used to describe the data from these power plants. Using the structure from this ontology, we have developed both R and Python packages for the automation of the FAIRification process. We have also developed an R package that automates the analysis of the quality of a data set through the designation of letter grades. With access to large time series data sets across many power plants, we can utilize spatiotemporal coherence between the sites in order to improve the quality of our data. To solve the issue of data missingness, we propose the use of Spatiotemporal-GNN autoencoders to detect and impute missing values from a data set by utilizing data from power plants nearby.

14 SOLAR ENERGY↗
DOE OSTI2024

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2021

Maximizing oil production from water alternating gas (CO 2 ) injection into residual oil zones: The impact of oil saturation and heterogeneity

Residual oil zones (ROZs) are widespread reservoirs, characterized by oil at residual saturation, either underlying oil fields (brownfield) or lateral (greenfield) to such fields. These reservoirs have the potential to produce volumes of oil sufficiently significant to make appreciable impacts on the US’s oil reserves and associated incidental CO 2 sequestration. The objective of this study is to improve our understanding the impact of heterogeneous and low oil saturations, in brownfield ROZs, on the effectiveness of water alternating gas (WAG) injection strategies. ROZs occur in the Permian Basin and elsewhere, and operators are using CO 2 injection for enhanced oil recovery (EOR) in these zones. The consensus model for the formation of ROZs is that they were formed by the effect of faster regional aquifer flow, acting over millions of years. Both the magnitude of oil saturation and the spatial distribution of oil differ from water-flooded main pay zones (MPZs). To explore the most effective injection strategies, we conducted simulations of CO 2 injection into synthetic geologic reservoirs. These simulations focused on injection into reservoirs subject to either man-made waterflooding or long-term natural waterflooding. By exploring the impact of varying: oil saturation; well patterns; reservoir heterogeneity; and permeability anisotropy, we attempt to quantify the factors that most influence the effectiveness of WAG injection. WAG ratios (the ratio of injected water and CO 2 , in reservoir volumes) of interest are those that either minimize the net CO 2 utilization ratios or maximize oil production rates. In general, the most effective WAG ratios for ROZs, are consistently less than those observed undergoing CO 2 injection in the same geologic reservoir models after traditional (man-made) waterflooding. This work demonstrates that most favorable WAG ratios for oil production in ROZs are different from those in traditional MPZs because of oil saturation differences. Therefore, CO 2 injection into both zones or directly copying WAG injection designs from MPZs to ROZs might not maximize oil production.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Engineering Self-Assembled Domain Asymmetry in Solvent Vapor Annealed Block Copolymer–Homopolymer Blend Films

Thin-film block copolymer (BCP) self-assembly is a powerful approach to generate highly uniform nanopatterns across large areas, yet the symmetries of these nanopatterns are constrained by the relative volume occupied by each polymer block. Here, we present a conceptually new approach to circumvent this limitation by combining homopolymer (HP) blending with solvent vapor annealing (SVA), demonstrated here for ternary blends of a near-symmetric BCP with athermal, low molar mass HPs. Screening of BCP–HP interactions by a weakly selective solvent promotes entropically driven delocalization of HP throughout both blocks and enables assembly of metastable lamellae for volume fractions as high as 0.78. Subsequent brief thermal annealing induces HPs to withdraw to their enthalpically favored domains, sharpening domain interfaces on a time scale much shorter than pattern coarsening. This renders lamellar nanopatterns with tunable widths that can be transferred to other materials with high fidelity. Furthermore, the persistence of the metastable asymmetric lamellae upon thermal annealing is sensitive to film confinement, as they are preserved in submonolayer films but transition to horizontal cylinders in films more than one monolayer thick. SVA using a strongly selective solvent results in assembled morphologies aligned closely with expectations based on the total polymer blend composition, underscoring the key role of BCP–HP interactions in dictating domain asymmetry. Overall, this work details important principles for using SVA to mediate BCP–HP interactions in thin films, thereby presenting opportunities to engineer pathways for the assembly of well-ordered nanopatterns with designer feature asymmetry for lithographic or nanotexturing applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Investigating the origin of the far-field reflection interference fringe (RIF) of microdroplets

We show that the reflection interference fringe (RIF) is formed on a screen far away from the microdroplets placed on a prism-based substrate, which have low contact angles and thin droplet heights, caused by the dual convex–concave profile of the droplet, not a pure convex profile. The geometric formulation shows that the interference fringes are caused by the optical path difference when the reflected rays from the upper convex profile at the droplet–air interface interfere with reflection from the lower concave profile at oblique angles lower than the critical angle. Analytic solutions are obtained for the droplet height and the contact angle out of the fringe number and the fringe radius in RIF from the geometric formulation. Furthermore, the ray tracing simulation is conducted using the custom-designed code. The geometric formulation and the ray tracing show excellent agreement with the experimental observation in the relation between the droplet height and the fringe number and the relation between the contact angle and the fringe radius. This study is remarkable as the droplet's dual profile cannot be easily observed with the existing techniques. However, the RIF technique can effectively verify the existence of a dual profile of the microdroplets in a simple setup. In this work, the RIF technique is successfully developed as a new optical diagnostic technique to determine the microdroplet features, such as the dual profile, the height, the contact angle, the inflection point, and the precursor film thickness, by simply measuring the RIF patterns on the far-field screen.

42 ENGINEERING↗
DOE OSTI2026

Improving ICARUS track reconstruction algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino program at Fermilab. Its primary objective is to explore the possible existence of sterile neutrinos in the O(1 eV) mass range and to clarify the anomalies observed in the Liquid Scintillator Neutrino Detector and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber, capable of producing high-resolution 3D images and precise calorimetric measurements of ionizing particles. This technology allows for a detailed study of neutrino interactions across a broad energy range, from a few keV to several hundred GeV. The track reconstruction is achieved through a software framework that applies a series of pattern recognition algorithms, transforming raw detector signals into fully reconstructed event topologies. This process involves identifying interaction vertices, particle tracks, and electromagnetic showers within the TPC. However, in certain cases, these algorithms may mistakenly break a single particle track into several shorter segments, interpreting each as a distinct particle. Since track length is used to estimate the particle's energy, such fragmentation can result in an energy underestimation of several hundred MeV. Furthermore, when a track is split into multiple segments, the particle identification (which relies on analyzing the energy loss as a function of the residual range) may fail, potentially leading to the loss of the entire event. To mitigate this problem, we have developed a dedicated algorithm designed to identify and reconnect (“stitch”) the tracks that were erroneously divided into multiple segments.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗
DOE OSTI2025

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Effect of resonant magnetic perturbations including toroidal sidebands on magnetic footprints and fast ion losses in HL-2M

Externally applied resonant magnetic perturbations (RMPs), generated by magnetic coils located outside the plasma (referred to as RMP coils), provide an effective way to control the edge localized mode (ELM) in tokamak devices. Due to the discrete nature of the toroidal distribution of these window-frame coils, toroidal sidebands always exist together with the fundamental harmonics designed for ELM control. In this work, the MARS-F code (Liu et al 2000 Phys. Plasmas7 3681) is applied to investigate the detailed features of the RMP spectra considering both the dominant harmonic (n = 2) and the associated sideband (n = 6), and the impact of the combined fields on magnetic footprints as well as on the fast ion losses for a reference double-null scenario in the HL-2M device. It is found that the sum of the n = 2 and n = 6 RMP fields splits the footprint and widens the footprint area, as compared to the single-n (n = 2) harmonic case. The resistive plasma response breaks the up–down symmetry of the footprint pattern on the outer divertor plates, which is otherwise symmetric assuming vacuum RMP fields. Considering fast ion losses, a threshold value exists for the initially launched radial position of test particles, as well as for the RMP coil current, before the loss occurs. When the threshold criterion is satisfied, the combined n = 2 and n = 6 RMP fields enhance the fast ion loss rate by , as compared to that of the n = 2 component alone. These results illustrate the important role of the sideband of RMP fields on the magnetic footprints and fast ion losses in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Lanthanum-Promoted Electrocatalyst for the Oxygen Evolution Reaction: Unique Catalyst or Oxide Deconstruction?

A conventional performance metric for electrocatalysts that promote the oxygen evolution reaction (OER) is the current density at a given overpotential. However, the assumption that increased current density at lower overpotentials indicates superior catalyst design is precarious for OER catalysts in the working environment, as the crystalline lattice is prone to deconstruction and amorphization, thus greatly increasing the concentration of catalytic active sites. We show this to be the case for La 3+ incorporation into Co 3 O 4 . Powder X-ray diffraction (PXRD), Raman spectroscopy and extended X-ray absorption fine structure (EXAFS) reveal smaller domain sizes with decreased long-range order and increased amorphization for La-modified Co 3 O 4 . This lattice deconstruction is exacerbated under the conditions of OER as indicated by operando spectroscopies. The overpotential for OER decreases with increasing La 3+ concentration, with maximum activity achieved at 17% La incorporation. HRTEM images and electron diffraction patterns clearly show the formation of an amorphous overlayer during OER catalysis that is accelerated with La 3+ addition. O 1s XPS spectra after OER show the loss of lattice-oxide and an increase in peak intensities associated with hydroxylated or defective O-atom environments, consistent with Co(O) x (OH) y species in an amorphous overlayer. Furthermore, our results suggest that improved catalytic activity of oxides incorporated with La 3+ ions (and likely other metal ions) is due to an increase in the number of terminal octahedral Co(O) x (OH) y edge sites upon Co 3 O 4 lattice deconstruction, rather than enhanced intrinsic catalysis.

Catalysts↗
DOE OSTI2021

Atmospheric Carbon and Transport – America (ACT-America) Data Sets: Description, Management, and Delivery

The ACT-America project is a NASA Earth Venture Suborbital-2 mission designed to study the transport and fluxes of greenhouse gases. The open and freely available ACT-America data sets provide airborne in situ measurements of atmospheric carbon dioxide, methane, trace gases, aerosols, clouds, and meteorological properties, airborne remote sensing measurements of aerosol backscatter, atmospheric boundary layer height and columnar content of atmospheric carbon dioxide, tower-based measurements, and modeled atmospheric mole fractions and regional carbon fluxes of greenhouse gases over the Central and Eastern United States. In this study, we conducted 121 research flights during five campaigns in four seasons during 2016–2019 over three regions of the US (Mid-Atlantic, Midwest and South) using two NASA research aircraft (B-200 and C-130). We performed three flight patterns (fair weather, frontal crossings, and OCO-2 underflights) and collected more than 1,140 h of airborne measurements via level-leg flights in the atmospheric boundary layer, lower, and upper free troposphere and vertical profiles spanning these altitudes. We also merged various airborne in situ measurements onto a common standard sampling interval, which brings coherence to the data, creates geolocated data products, and makes it much easier for the users to perform holistic analysis of the ACT-America data products. Here, we report on detailed information of data sets collected, the workflow for data sets including storage and processing of the quality controlled and quality assured harmonized observations, and their archival and formatting for users. Finally, we provide some important information on the dissemination of data products including metadata and highlights of applications of ACT-America data sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Enhancing photonic systems using topology and non-Hermiticity

The broad goal of this project was to develop new analytical and numerical insights for how important photonic processes can be improved using recently discovered principles in topological and non-Hermitian physics. In particular, there are two recent discoveries that we aimed to harness to achieve this goal. First, it was discovered in condensed matter physics that crystalline symmetries can protect low-dimensional topologically protected states in lattices without the need for breaking time-reversal symmetry. These so-called ‘higher-order’ topological systems represent an important development for photonic systems, where it is very difficult to break time-reversal symmetry, and which had been previously thought necessary to realize topological phenomena. Second, the last decade has seen a significant amount of interest in phenomena which are unique to non-Hermitian systems, i.e., systems which do not conserve energy. For example, spatially patterned gain and loss can be used to realize exceptional points, which are degeneracies in a system’s spectrum where the system becomes defective, while the existence of radiative losses also enables a new route to confinement through bound states in the continuum. For such non-Hermitian phenomena, photonics again represents a critical platform, as photonic systems naturally lose light to their radiative environments, making them generally non-Hermitian, and it is also possible to incorporate additional gain or loss. Based on these broad principles, we pursued a range of projects to harness bound states in the continuum in a variety of different systems and architectures, develop real-space methods for classifying topological systems to yield better photonic design principles, and a novel Brillouin-based fiber laser for sensing strain.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Fast and Flexible Inference Framework for Continuum Reverberation Mapping Using Simulation-based Inference with Deep Learning

Continuum reverberation mapping (CRM) of active galactic nuclei (AGN) monitors multiwavelength variability signatures to constrain accretion disk structure and supermassive black hole (SMBH) properties. The upcoming Vera Rubin Observatory’s Legacy Survey of Space and Time will survey tens of millions of AGN over the next decade, with thousands of AGN monitored with almost daily cadence in the deep drilling fields. However, existing CRM methodologies often require long computation time and are not designed to handle such large amounts of data. In this paper, we present a fast and flexible inference framework for CRM using simulation-based inference (SBI) with deep learning to estimate SMBH properties from AGN light curves. We use a long short-term memory summary network to reduce the high dimensionality of the light curve data and then use a neural density estimator to estimate the posterior of SMBH parameters. Using simulated light curves, we find SBI can produce more accurate SMBH parameter estimation with 10 3 –10 5 times speed up in inference efficiency compared to traditional methods. The SBI framework is particularly suitable for wide-field CRM surveys as the light curves will have identical observing patterns, which can be incorporated into the SBI simulation. We explore the performance of our SBI model on light curves with irregular-sampled, realistic observing cadence and alternative variability characteristics to demonstrate the flexibility and limitation of the SBI framework.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Comparative Study of Layer Heating and Continuous Heating Methods on Prediction Accuracy of Residual Stresses in Selective Laser Melted Tube Samples

Thermal distortion and residual stresses are two important factors that affect the quality and reliability of steel parts manufactured by laser powder bed fusion (LPBF) processes. A cost-effective model for evaluation of those heat effects is needed to refine the manufacturing process and provides insights into the product design and heat treatment. In this study, the layer heating method and sophisticated track-layer scanning method were applied to simulate the thermo-mechanical response of IN625 tube parts built by LPBF. Based on the similarity of temperature field in each layer deposit, a swept mesh was constructed to perform the thermal analysis for top layer, with the rest of layers referring to the temperature by node number offsetting. A novel explicit finite element analysis code accelerated by graphics processing unit was used for the massive-element numerical analysis. The computational accuracy and efficiency of the layer heating and track-layer scanning methods were compared in detail. It is shown that layer heating method can efficiently capture the pattern of stress distribution with reasonable accuracy in stress magnitude. The grouped track-layer scanning method can predict the residual stress and strain more accurately at a higher cost (5 ~ 10×). The elastic strain distribution was compared with the measurement by X-ray diffraction, confirming the accuracy of residual stress prediction.

36 MATERIALS SCIENCE↗