Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process coupling”

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 469 records · Page 26

Quantifying physical parameterization uncertainties associated with land-atmosphere interactions in the WRF model over Amazon

The Weather Research and Forecasting (WRF) model can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations but subjected to uncertainties associated with model physical parameterizations. In this study, we propose a framework to quantify and reduce model physical parameterization uncertainties associated with surface fluxes and L-A coupling. An ensemble of WRF simulations with different physical schemes is used to simulate surface fluxes and land-atmosphere coupling strength over the Amazon region. The physical parameterizations investigated include cloud microphysics (MP), land surface processes (LSM), planetary boundary layer (PBL), surface layer (SL), and cumulus (CU). We perform 120 ensemble simulations using the WRF model and different combinations of six MPs, three LSMs, six PBLs and SLs and three CUs. The measurements from the GoAMAZON field campaign and satellite data are used to evaluate model performance. A Multi-way analysis of variance (ANOVA) approach is applied to quantify the relative importance of different physics processes on L-A coupling. The Tukey’s test is used to sort schemes that have no significant differences into one group. The suite of physics that result in the best simulations of the corresponding variables are selected based on the Taylor skill score. Results show that the relative importance of processes and their interaction vary with the variables of interest. For example, CU was the most important process in modulating soil moisture, 2 m-humidity, latent heat, and net radiation. LSM showed dominant effects on 2 m-temperature and also has the largest impact on sensible heat and the lifting condensation level. The best physical parameterization ensembles show much narrower ranges of the variables of interest than the priori ensemble. Results of this study show the roles of different physical processes in modulating L-A interactions, quantify model uncertainties from physical processes, and provide insights for improving the model physics parameterizations.

Wang, Chen↗

Automated processing of environmental transmission electron microscopy images for quantification of thin film dewetting and carbon nanotube nucleation dynamics

We report scalable production of carbon nanotubes (CNTs) requires catalysts and reaction conditions that provide high nucleation efficiency. In situ characterization methods such as environmental transmission electron microscopy (ETEM) can reveal fundamental mechanisms of synthesis, but to date have primarily provided qualitative observations on small sample sizes. Here, quantitative analysis is performed using high-resolution, high-rate video capture of ETEM experimentation coupled with automated image processing, involving computer vision algorithms and convolutional neural networks. By this approach, we detect distinct nanoparticle formation from an alumina-supported iron thin film and subsequent CNT nucleation from the nanoparticles. The statistical summary of particles in each video shows that, compared to a H 2 -only atmosphere, pretreatment of the catalyst with carbon added to the H 2 atmosphere results in a smaller average particle diameter, a 2-fold increase in particle density (to 5300 particles/μm 2 ), a 3-fold increase in CNT nucleation efficiency (to 92%), and more than a 5-fold increase in CNT density (to 4800/μm 2 ). Addition of carbon during exposure to H 2 is also more effective than NH 3 at dewetting the catalyst film and increasing the CNT nucleation efficiency, in spite of NH 3 being a stronger reducing agent for iron. Insights from this study are applicable to improving CNT yield and productivity in both batch-style and continuous processes.

25 ENERGY STORAGE↗

Graph-Based Representations and Applications to Process Simulation

Rapid and robust convergence of a process flowsheet is critical to enable large-scale simulations that address core scientific questions related to process design, optimization, and sustainability. However, due to the highly coupled and nonlinear nature of chemical processes, efficiently solving a flowsheet remains a challenge. In this work, we show that graph representations of the underlying physical phenomena in unit operations may help identify potential avenues to systematically reformulate the network of equations and enable more robust topology-based convergence of flowsheets. To this end, we developed graph abstractions of the governing equations of vapor-liquid and liquid-liquid equilibrium separation equipment. These graph abstractions consist of a mesh of interconnected variable nodes and equation nodes that are systematically generated through PhenomeNode, a new open-source library in Python developed in this study. We show that partitioning the graph into separate mass, energy, and equilibrium subgraphs can help decouple nonlinearities and guide decomposition algorithms. By employing the graph abstraction on an industrial separation process for separating glacial acetic acid from water, we implemented a new block decomposition scheme in BioSTEAM and demonstrated that this can accelerate convergence over a traditional sequential modular approach.

Distillation↗

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗

Reactive chemical transport simulations of geologic carbon sequestration: Methods and applications

Chemical reaction simulations are considerably used to quantitatively assess the long-term geologic carbon sequestration (GCS), such as CO 2 sequestration capacity estimations, leakage pathway analyses, enhanced oil recovery (EOR) efficiency studies, and risk assessments of sealing formations (caprocks), wellbores, and overlying underground water resources. All these require a deep understanding of the CO 2 - associated chemical reactions. To ensure long-term, safe CO 2 sequestration in the intended formations, modeling is the only way to plausibly assess the CO 2 flow, reaction, and transport over thousands of years. This review summarizes the multiple methodologies for describing homogeneous and heterogeneous chemical reaction patterns and multiscale application examples, the recent progress and current status of chemical reaction simulations for GCS, and the impact of such simulations on geological CO2 sequestration performance. Technical gaps and future challenges are also discussed for further study. The trends and challenges of such studies include: (1) the combination of coupled chemical, mechanical, and transport processes with calibrated experiments and associated uncertainty/risk assessments; (2) enhancement of the ability to simulate detailed geophysical and geochemical equations to mimic in situ conditions; and (3) characterization of multiscale subsurface systems with detailed conceptual models and assignment of suitable boundary conditions for field-scale sequestration fields. One major gap remaining is the current lack of accurate (and scale-justified) kinetic and equilibrium chemical reaction parameters under reservoir conditions. Advanced models that couple chemical, mechanical, and transport processes with scale-justified parameters, from lab to field-scale experiments, are required for quantitative assessments of sequestration capacity and the long-term safety of GCS projects.

58 GEOSCIENCES↗

BNF Radar b1 Data Processing Report: Spring 2025

The U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) User Facility supports atmospheric science through an integrated network of fixed and mobile observatories. These facilities collect continuous and campaign-based observations of atmospheric properties, with the goal of improving the representation of clouds, aerosols, precipitation, and radiation in Earth system models. The Bankhead National Forest (BNF) site, established as an ARM Mobile Facility (AMF) on 1 October 2024, is situated in a forested region of northern Alabama. Its strategic location in a southeastern U.S. environment characterized by complex terrain, diverse land cover, and frequent convective storms provides a valuable opportunity to examine coupled land-atmosphere processes under natural variability.

54 ENVIRONMENTAL SCIENCES↗

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: LDRD Final Project Poster

Electric-field assisted sintering (EFAS) manufacturing techniques reduce energy usage requirements, fabricate near-net-shape parts, and produce bulk nanostructured materials through rapid heating and cooling rates achievable through Joule heating. Process parameters of pulsed DC voltage and constant pressure are applied through stainless steel rams and graphite tooling to the ceramic yttrium oxide powder to be sintered. Our high-fidelity modeling code application advances manufacturing for extreme environments by enabling process-informed design and improving microstructure consistency through multiphysics multiscale process-structure-property-performance correlation simulations.

36 MATERIALS SCIENCE↗

The U.S. DOE ARM User Facility Establishes a New Site for Studies of Land–Aerosol–Cloud Interactions in the Southeastern United States

The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility has established a new site in the Bankhead National Forest (BNF) in northern Alabama that will gather data on how clouds, the land surface, and aerosols interact at a hierarchy of scales important to understanding and simulating the Earth system. Starting its operations in October 2024, the BNF site provides a multi-year opportunity for scientists to unravel complex land-atmosphere interactions. A suite of ground-based sensors, elevated tower-based instrumentation, and aerial facilities will enable scientists to investigate those interactions from within the canopy to the clouds. The Southeastern United States was recommended by the DOE ARM and its collaborators in the broader community as an important region to address their common scientific questions given the region’s abundant surface-forced convective clouds and mesoscale convective systems that pose ongoing challenges in Earth system models. The region is also home to significant terrain complexity and land-use heterogeneity that will unleash new understanding of anthropogenic and biogenic aerosol processes, boundary-layer aerosol-cloud interactions, and the interactions between the terrestrial ecosystem and coupled aerosol-cloud-radiation processes.

Aerosols↗

Flexible FlueCO2

Carbon dioxide (CO2) emission reductions remain a significant challenge on the path to clean energy. There are increasing legislative, social, and environmental factors motivating CO2 emissions reduction from power plants with carbon capture and storage (CCS). CCS in natural gas combined cycle (NGCC) power plants is critical to achieve a net-zero carbon electricity grid. Enhanced 45Q tax credits provide new incentives, but currently available technologies are unable to profitably operate in grids with deep variable renewable penetration which require flexible NGCC operation. Luna Labs has developed the FlueCO2 membrane to enable a profitable NGCC-CCS process. The FlueCO2 membrane couples steam transport across the membrane to CO2 transport in the opposite direction, enabling high capture efficiencies and low energy costs even at low CO2 concentrations. The dual-phase membrane can operate in the range of typical flue gas temperatures and pressures and does not require temperature or pressure cycling. Luna Labs’ FlueCO2 technology enables flexible and profitable operation of NGCC plants with lower capital investment and impact on electricity prices. In this Phase 1 project, Luna Labs utilized experimental testing, modeling, process simulation, and standardized costing methodologies to evaluate the techno-economic value of a 650 MW greenfield NGCC plant with FlueCO2 (NGCC-FlueCO2). Key design requirements for operation were established and plant performance under load-leveling conditions was validated through computational fluid dynamics and process modeling. Luna Labs developed a dynamic modeling tool which modeled plant operational modes across a variety of tax structures and electricity pricing scenarios to project the overall Net Present Value (NPV) of the NGCC-FlueCO2. FlueCO2 minimizes the impact of CCS integration on plant operation by integrating directly into the NGCC heat recovery steam generator (HRSG). By tapping into the plant’s low-pressure (LP) steam, operators can divert LP steam to the greenfield NGCC and/or CCS process in response to dynamic markets. Since FlueCO2 will not significantly affect HRSG (or NGCC) operation, CCS only turns off during peak power demand (>$250/MWh). Under baseload conditions, FlueCO2 lowers the capital (37%), energy (36%) and carbon capture (<$40/tonne) costs and can increase the overall plant lifetime NPV by approximately ~$1B in comparison with NGCC solvent-based capture reference cases (NETL Case 31B). Luna Labs has shared its costing tools with several interested partners and customers, which follows a generalizable approach to costing analysis.

Kelly, Jesse↗

Ionization-driven competitive (recovery) process in pre-damaged KTaO 3 : A brief review

The nuclear (S n ) and electronic (S e ) energy dissipation processes have been considered to be independent and largely uncorrelated, influencing our understanding of ion–solid interaction and damage processes in the last decades. Recently, however, it has become more generally accepted that S n and S e are coupled as they interact both in time and space. To decouple these processes, separating these effects in experiments using sequential dual-beam irradiations have become accepted as the logical path to advance the understanding of complex interactions between S e and pre-existing defects that may be created from displacement events. This experimental approach has been recently applied to studies of KTaO 3 to reveal new insights into this critical research topic. Here, we offer a forward-looking and comprehensive perspective on the fundamental coupling between Se and pre-existing defects in KTaO 3 . The origins behind the competitive two-stage phase transition process leading to damage healing are revealed and discussed. Furthermore, the evidence resulting from synergistic effects is also included for comparison. Additionally, our findings are rationalized using both Se and the ion velocity as key parameters. We highlight how the inelastic thermal spike (i-TS) calculations provide insights into the nature of this coupled process and further confirm that the ion velocity effect governs annealing kinetics. This work emphasizes that through the introduction of a small amount of local disorder in materials, MeV ion irradiation (i.e., not extreme S e ) may also be one additional option in subsequent material modification and functionalization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Elucidating the Link Between Alkali Metal Ions and Reaction-Transport Mechanisms in Cathode Electrodes for Alkali-ion Batteries

Our long-term goal is to improve the reliability of electrode materials and their ability to transport and store various metal ions for electrochemical energy storage applications. The main objective of this work was to investigate the intrinsic relationship between the role of alkali metal ions and electrochemically driven mechanical stability and kinetic properties of battery materials. The overall question was “What is the role of alkali metal ions on the electrochemical and mechanical behavior of cathode electrodes? Our guiding hypothesis was that intercalation of larger alkali metal ions (Na and K) inevitably alters the coupled transport-reaction processes during battery operation in organic electrolytes, leading to more intensive chemo-mechanical instabilities in cathode electrodes, resulting in rapid capacity fade. To validate the hypothesis, we experimentally characterized the reaction-transport processes and governing forces driving the instability of electrode materials in different alkali metal-ion environments. The project had three main tasks. The first one was to investigate intercalation-induced strains and associated stress generation, and their impact on structural deformations in composite cathode electrodes. The second task focused on identifying potential-dependent dynamic changes in the electrode-electrolyte interface in alkali metal ion batteries. The last task was focused on determining how larger alkali metal ions with slower diffusivity affect the transport-mechanics coupling at faster scan rates, compared to smaller ions with faster diffusivity in electrodes. We shortly provided the outcome of each task in the accomplishment section. This project produced 10 peer-reviewed publications (9 research papers and one review manuscript) and supported two Ph.D. students, who graduated from Oklahoma State University.

25 ENERGY STORAGE↗

Coordination and control – limits in standard representations of multi-reservoir operations in hydrological modeling

Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.

54 ENVIRONMENTAL SCIENCES↗

High-speed phase contrast imaging of spray breakup of jet fuels under combusting conditions

X-ray phase contrast imaging is conducted at 90,517 Hz on a fuel spray to characterize the breakup of the bulk fuel into ligaments and individual droplets in a single-sector, swirl-stabilized combustor under reacting conditions. Phase contrast imaging is used to qualitatively assess spray breakup for two fuels: Jet-A, which represents a normal viscosity fuel, and a blend of JP-5 (64%) and farnesane (36% by volume), which represents a high-viscosity fuel, at three fuel-flow rates and two air inlet temperatures. Time-series images of the fuels show that atomization occurs much more rapidly at higher fuel pressure with thinner ligaments being rapidly shredded. The high-viscosity fuel is qualitatively observed to have markedly longer, thicker ligaments and larger droplet diameters after the primary breakup than the normal viscosity fuel. Using an image processing and droplet quantification strategy, the mean droplet velocity magnitudes were found to decrease with decreasing fuel pressure and fuel-flow rate. The increase in preheat temperature from 323 to 370 K results in an increase in the mean droplet velocity magnitude and a decrease in numerical average and Sauter mean diameter, indicating improved atomization. The high-viscosity fuel was found to have larger average droplet diameters at each condition than the normal viscosity fuel. The increased preheat is found to enhance the combustion efficiency, which is hypothesized to couple with the improved atomization due to the higher combustion temperatures. This indicates the importance of examining spray atomization and breakup under realistic combusting conditions for relevant fuels as the heat release directly couples with the atomization process.

42 ENGINEERING↗

On the Spectral Evolution of Hot White Dwarf Stars. II. Time-dependent Simulations of Element Transport in Evolving White Dwarfs with STELUM

White dwarf stars are subject to various element transport mechanisms that can cause their surface composition to change radically as they cool, a phenomenon known as spectral evolution. In this paper, we undertake a comprehensive theoretical investigation of the spectral evolution of white dwarfs. First, we introduce STELUM, a new implementation of the stellar evolutionary code developed at the Université de Montréal. We provide a thorough description of the physical content and numerical techniques of the code, covering the treatment of both stellar evolution and chemical transport. Then, we present two state-of-the-art numerical simulations of element transport in evolving white dwarfs. Atomic diffusion, convective mixing, and mass loss are considered simultaneously as time-dependent diffusive processes and are fully coupled to the cooling. We first model the PG 1159-DO-DB-DQ evolutionary channel: a helium-, carbon-, and oxygen-rich PG 1159 star transforms into a pure-helium DB white dwarf due to gravitational settling and then into a helium-dominated, carbon-polluted DQ white dwarf through convective dredge-up. We also compute for the first time the full DO-DA-DC evolutionary channel: a helium-rich DO white dwarf harboring residual hydrogen becomes a pure-hydrogen DA star through the float-up process and then a helium-dominated, hydrogen-bearing DC star due to convective mixing. We demonstrate that our results are in excellent agreement with available empirical constraints. In particular, our DO-DA-DC simulation perfectly reproduces the lower branch of the bifurcation observed in the Gaia color–magnitude diagram, which can therefore be interpreted as a signature of spectral evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultrafast population and structural dynamics of a Ni-bipyridine photoredox catalyst reveal a significant deactivation pathway

The ultrafast excited state pathways and dynamics of NiII-bipyridine complexes influence the yield of photochemical processes involved in their catalytic cross-coupling reactions. Here we present ultrafast Ni K x-ray emission spectroscopy (XES) and x-ray solution scattering (XSS) of a NiII-bipyridine aryl halide complex, [Ni(t-Bubpy)(o-tol)Br], to quantify the excited state population dynamics and structural changes of the pre-catalyst. Due to the local spin-sensitivity of XES, the population dynamics of metal-to-ligand charge transfer (MLCT) and metal-centered (MC) excited states is established. A rapid ground state recovery pathway is newly identified, representing a significant deactivation pathway during photocatalysis. Furthermore, the pseudotetrahedral structure of the long-lived MC excited state is unambiguously identified and refined by XSS. The results advance our understanding of the ultrafast relaxation mechanisms that impact the photocatalytic mechanism and yield for NiII-bipyridine aryl halide cross-coupling catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effect of Many Modes on Self-Polarization and Photochemical Suppression in Cavities

The standard description of cavity-modified molecular reactions typically involves a single (resonant) mode, while in reality the quantum cavity supports a range of photon modes. In this work we demonstrate that as more photon modes are accounted for, physico-chemical phenomena can dramatically change, as illustrated by the cavity-induced suppression of the important and ubiquitous process of proton-coupled electron-transfer. Using a multi-trajectory Ehrenfest treatment for the photonmodes, we find that self-polarization effects become essential, and we introduce the concept of selfpolarization- modified Born-Oppenheimer surfaces as a new construct to analyze dynamics. As the number of cavity photon modes increases, the increasing deviation of these surfaces from the cavityfree Born-Oppenheimer surfaces, together with the interplay between photon emission and absorption inside the widening bands of these surfaces, leads to enhanced suppression. The present findings are general and will have implications for the description and control of cavity-driven physical processes of molecules, nanostructures and solids embedded in cavities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗