Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “micro scale models”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Chemistry and Dynamics of Supercritical Carbon Dioxide and Methane in the Slit Pores of Layered Silicates

In the mid 2010s, high-pressure diffraction and spectroscopic tools opened a window into the molecular-scale behavior of fluids at the conditions of many CO 2 sequestration and shale/tight gas reservoirs - conditions where CO 2 and CH 4 are present as variably wet supercritical fluids. Integrating high-pressure spectroscopy and diffraction with molecular modeling has revealed much about the ways that supercritical CO 2 and CH 4 behave in reservoir components, particularly in the slit-shaped micro- and meso-pores of layered silicates (phyllosilicates) abundant in caprocks and shales. This Account summarizes how supercritical CO 2 and CH 4 behave in the slit pores of swelling phyllosilicates as functions of the H 2 O activity, framework structural features, and charge-balancing cation properties at 90 bar and 323 K – conditions similar to a reservoir at ~1 km depth. Slit pores containing cations with large radii, low hydration energy, and large polarizability readily interact with CO 2 , allowing CO 2 and H 2 O to adsorb and co-exist in these interlayer pores over a wide range of fluid humidities. In contrast, cations with small radii, high hydration energy, and low polarizability weakly interact with CO 2 leading to reduced CO 2 uptake and a tendency to exclude CO 2 from interlayers when H 2 O is abundant. The reorientation dynamics of confined CO 2 depends on the interlayer pore height, which is strongly influenced by the cation properties, framework properties, and fluid humidity. The silicate structural framework also influences CO 2 uptake and behavior - for example, smectites with increasing F-for-OH substitution in the framework take up greater quantities of CO 2 . Reactions that trap CO 2 in carbonate phases have been observed in thin H 2 O films near smectite surfaces, including a dissolution-reprecipitation mechanism when the edge surface area is high and an ion exchange-precipitation mechanism when the interlayer cation can form a highly insoluble carbonate. In contrast, supercritical CH 4 does not readily associate with cations, does not react with smectites, and is only incorporated into interlayer slit mesopores when (i) the pore has a $z$-dimension large enough to accommodate CH 4 , (ii) the smectite has low charge and (iii) the H 2 O activity is low. In conclusion, the adsorption and displacement of CH 4 by CO 2 and vice-versa have been studied on the molecular scale in one shale, but opportunities remain to examine behavioral details in this more complicated, slit-pore inclusive system.

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

An FEM-Based Peridynamic Model for Failure Analysis of Unidirectional Fiber-Reinforced Laminates

To predict the mixed damage modes of unidirectional fiber-reinforced polymer (FRP) laminates under dynamic loading, an FEM-based peridynamic model is introduced in this paper. Based on its geometric structure and material composition, a long fiber lamina is considered a transversely isotropic medium as a result of homogenization at the meso-scale. The laminated structure is modeled by stacking surface mesh layers with arbitrary fiber angles along the thickness direction. The peridynamic bonds between Gauss points connect the separated elements. These bonds are classified as inner-layer bonds and inter-layer bonds. To represent the anisotropy of a laminate, the micro-elastic modulus of the inner-layer and inter-layer bonds is calculated from the anisotropic engineering material constants separately. To capture complex failure behaviors of laminate structures, an empirical damage model is proposed for the tension/ compression breakage of peridynamic bonds. This damage model can control the in-plane and delamination failure process. Finally, benchmark tests are conducted to validate the elastic response of laminates under dynamic loading. In terms of damage analysis, the proposed model can capture the complex damage modes and resistive force of laminate structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

If We Build Them, They Will Run: Automated HPC Apps Deployment and Profiling with eBPF in Cloud

The high performance computing (HPC) community is in a period of transition. The rise of AI/ML coupled with a changing landscape of resources deems portability a new metric of performance, and methods to move between on-premises and cloud environments and assess compatibility are paramount. Here we design and test a strategy for bridging the gap between traditional HPC and Kubernetes environments – first containerizing applications, providing automated orchestration to run studies, and packaging the setup with automated means to assess performance using low overhead eXtended Berkeley Packet Filter (eBPF) programs. We first assess different designs for eBPF collection, demonstrating a tradeoff between number of programs deployed on a node and overhead added. We develop 5 low overhead eBPF programs that combine with streaming ML models to assess CPU, futex, TCP, shared memory, and file access across four different builds of an HPC application for CPU and GPU. We use eBPF data to generate insights into the possible underlying etiology of scaling issues. We then assess compatibility of a well-known benchmark, HPCG, across matrices of micro-architectures and optimization levels (217 containers across 24 instance types and over 7500 runs). We provide to the community 30 applications to deploy in our automated setup and perform a scaling study from 4 to a maximum of 256 nodes for both CPU and GPU applications. Finally, we use our gained knowledge about performance to generate compatibility artifacts that are used by a newly developed Kubernetes controller to intelligently select instance type based on optimizing a figure of merit. Along with insights to scaling in this environment with a collection of applications and templates to work from, we provide an overall strategy for approaching HPC application deployment and image selection based on compatibility in cloud.

Computer science↗
DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward a Holistic Performance Evaluation of Large Language Models Across Diverse AI Accelerators

Artificial intelligence (AI) methods have become critical in scientific applications to help accelerate scientific discovery. Large language models (LLMs) are being considered a promising approach to address some challenging problems because of their superior generalization capabilities across domains. The effectiveness of the models and the accuracy of the applications are contingent upon their efficient execution on the underlying hardware infrastructure. Specialized Al accelerator hardware systems have recently become available for accelerating Al applications. However, the comparative performance of these AI accelerators on large language models has not been previously studied. In this paper, we systematically study LLMs on multiple AI accelerators and GPUs and evaluate their performance characteristics for these models. We evaluate these systems with (i) a micro-benchmark using a core transformer block, (ii) a GPT-2 model, and (iii) an 1,I,M-driven science use case, GenSLM. We present our findings and analyses of the models' performance to better understand the intrinsic capabilities of AI accelerators. Furthermore, our analysis takes into account key factors such as sequence lengths, scaling behavior, and sensitivity to gradient accumulation steps.

Emani, Murali↗
DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗
DOE OSTI2024

Self-Organized Stress Distributions in Polycrystalline Materials [Dissertation]

Understanding stress distributions in solid materials is complicated by the fact that most materials are polycrystalline in nature, with each crystal having an elastically anisotropic reaction to force. This study is to gain a better understanding on how external forces placed on a polycrystal are related to internal reactions within and between the grains. The hypothesis is that the stress distribution in porous to fully dense materials are self-organized based on strong contacts between and within the individual grains created by force chains. Force chains, commonly known in loaded granular materials, and could be the phenomenon that connect micro to macro deformation. Scale bridging measurements conducted through Raman spectroscopy, Atomic Force Microscopy, and Digital Image Correlation will be used to create stress maps, modulus maps, and elastic strain maps across a variety of geological and pharmaceutical polycrystals. When possible, the resulting maps will be compared to current homogenization schemes and a full field models. Finite element modeling will be used to assess whether the patterning seen in the experimental map is a reasonable approximation based on the orientation data of the samples used. A minimum of three publications is projected to be accomplished focusing each on a different method to experimentally test and analyze stress distributions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bridging the length scales on mechanical property evaluation (Final Report)

The development of small-scale mechanical testing in combination with microstructural investigation is of great interest to the nuclear materials community for both materials development and monitoring applications. Dramatically reducing the sample sizes to reduce radioactivity and obtaining mechanical properties of irradiated samples is truly intriguing. Moreover, such studies promise a range of benefits including cost reduction, fundamental insight in structure-property relationships, increased statistics on less sample material, and reinvestigation of prior irradiated and tested reactor samples while simultaneously enabling the generation of mechanical test data on ion beam irradiated materials with limited penetration depths. Small scale materials testing on sub-sized samples has been studied for several decades, though it has only been after the development of micro-testing based on Focused Ion Beam (FIB) sample manufacturing in that orders of magnitude smaller samples could really be investigated in a quantitative manner. In recent years, small scale mechanical testing techniques at a number of length scales has been developed for both unirradiated and irradiated (ion and neutron) materials. Technological advances made in this field have enabled ex-situ and in-situ transmission electron microscopy (TEM) and scanning electron microscopy (SEM) examination, thus leading to more accurate measurements as well as additional mechanistic information. A recent review of the benefits of these techniques show that these techniques are at a stage to tackle multi-scale ranges of materials investigations and can be utilized to obtain fundamental science-based understanding of nuclear materials. Considering the tremendous advances made, one can see how small-scale mechanical testing techniques combined with modeling can enable true small scale to bulk scale mechanical property correlations. However, for the engineering community to adapt this approach fully, one needs to demonstrate that a) that these techniques can produce results with high fidelity and reproducibility, b) generate engineering stress-strain data that one can utilize to understand bulk behavior, and c) generate new insight into relevant phenomena fostering the true understanding of radiation damage and microstructure in materials for nuclear applications. It is the objective of this proposal to bridge the length scale between macro- and micro- scale mechanical testing of unirradiated and irradiated materials. This involves the development and demonstration of procedures for multi-scale mechanical testing that enable high fidelity reproducibility of data and the generation bulk property data from small-scale mechanical tests. Through this, the proposal aims to enhance the confidence in the obtained data at the smaller length scales and enhance the insight provided from these techniques for bulk scale applications on both unirradiated and irradiated materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Advancing Detector R&D for High-Pressure Gaseous Argon TPCs in Precision Neutrino Physics

High-pressure gaseous argon time projection chambers (HPgTPCs) represent an emerging detector paradigm for neutrino physics, combining increased target density with the intrinsic tracking and low thresholds of gaseous detectors. This approach enables detailed reconstruction of exclusive final states, improved particle identification, and sensitivity to low-energy and rare processes — capabilities that are increasingly central to precision oscillation measurements and searches for beyond-the-Standard-Model signatures. This abstract presents an overview of ongoing detector R&D toward high-pressure gaseous argon TPC operation, with emphasis on micro-pattern gas detector (MPGD) charge amplification in argon-based mixtures. We report experimental characterization of triple-GEM structures at pressures relevant for neutrino applications, including studies of multiplication factor scaling, stability, and operational voltage envelopes across gas admixtures. Measurements performed at the TOAD and GORG test stands at Fermilab help define viable amplification and electronics noise regimes in conditions where higher density imposes stricter constraints on signal formation. These results provide essential input to the optimization of high-pressure gaseous argon detectors for future neutrino experiments, including near-detector concepts such as ND-GAr in DUNE Phase II near detector upgrade. More broadly, this program shows how dedicated detector R&D can expand the precision frontier in neutrino physics by enabling complementary reconstruction capabilities beyond conventional detectors.

McConnell, Brenna [Indiana U.] (ORCID:000900041138↗
DOE OSTI2021

Micro on a macroscale: relating microbial-scale soil processes to global ecosystem function

ABSTRACT Soil microorganisms play a key role in driving major biogeochemical cycles and in global responses to climate change. However, understanding and predicting the behavior and function of these microorganisms remains a grand challenge for soil ecology due in part to the microscale complexity of soils. It is becoming increasingly clear that understanding the microbial perspective is vital to accurately predicting global processes. Here, we discuss the microbial perspective including the microbial habitat as it relates to measurement and modeling of ecosystem processes. We argue that clearly defining and quantifying the size, distribution and sphere of influence of microhabitats is crucial to managing microbial activity at the ecosystem scale. This can be achieved using controlled and hierarchical sampling designs. Model microbial systems can provide key data needed to integrate microhabitats into ecosystem models, while adapting soil sampling schemes and statistical methods can allow us to collect microbially-focused data. Quantifying soil processes, like biogeochemical cycles, from a microbial perspective will allow us to more accurately predict soil functions and address long-standing unknowns in soil ecology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Soil Carbon Saturation: What Do We Really Know?

Managing soils to increase organic carbon storage presents a potential opportunity to mitigate and adapt to global change challenges, while providing numerous co-benefits and ecosystem services. However, soils differ widely in their potential for carbon sequestration, and knowledge of biophysical limits to carbon accumulation may aid in informing priority regions. Consequently, there is great interest in assessing whether soils exhibit a maximum capacity for storing organic carbon, particularly within organo–mineral associations given the finite nature of reactive minerals in a soil. While the concept of soil carbon saturation has existed for over 25 years, recent studies have argued for and against its importance. Here, we summarize the conceptual understanding of soil carbon saturation at both micro- and macro-scales, define key terminology, and address common concerns and misconceptions. We review methods used to quantify soil carbon saturation, highlighting the theory and potential caveats of each approach. Critically, we explore the utility of the principles of soil carbon saturation for informing carbon accumulation, vulnerability to loss, and representations in process-based models. We highlight key knowledge gaps and propose next steps for furthering our mechanistic understanding of soil carbon saturation and its implications for soil management.

Environmental sciences↗
DOE OSTI2020

Nano-engineering the material structure of preferentially oriented nano-graphitic carbon for making high-performance electrochemical micro-sensors

Direct synthesis of thin-film carbon nanomaterials on oxide-coated silicon substrates provides a viable pathway for building a dense array of miniaturized (micron-scale) electrochemical sensors with high performance. However, material synthesis generally involves many parameters, making material engineering based on trial and error highly inefficient. Here, we report a two-pronged strategy for producing engineered thin-film carbon nanomaterials that have a nano-graphitic structure. First, we introduce a variant of the metal-induced graphitization technique that generates micron-scale islands of nano-graphitic carbon materials directly on oxide-coated silicon substrates. A novel feature of our material synthesis is that, through substrate engineering, the orientation of graphitic planes within the film aligns preferentially with the silicon substrate. This feature allows us to use the Raman spectroscopy for quantifying structural properties of the sensor surface, where the electrochemical processes occur. Second, we find phenomenological models for predicting the amplitudes of the redox current and the sensor capacitance from the material structure, quantified by Raman. Our results indicate that the key to achieving high-performance micro-sensors from nano-graphitic carbon is to increase both the density of point defects and the size of the graphitic crystallites. Our study offers a viable strategy for building planar electrochemical micro-sensors with high-performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Fracture Dynamics in Silicon Anode Solid-State Batteries

Solid-state batteries (SSBs) with silicon anodes could enable improved safety and energy density compared to lithium-ion batteries. However, degradation arising from the massive volumetric changes of silicon anodes during cycling is not well understood in solid-state systems. Here, we use operando X-ray computed microtomography to reveal micro- to macro-scale chemo-mechanical degradation processes of silicon anodes in SSBs. Mud-type channel cracks driven by biaxial tensile stress form across the electrode during delithiation. We also find detrimental cracks at the silicon/ solid electrolyte interface that form due to local reaction competition between neighboring domains of different sizes. Continuum phase-field damage modeling quantifies stress-driven channel cracking and shows that the lithiated silicon stress state is critical for determining the extent of interfacial fracture. This work reveals the mechanisms that govern SSBs compared to conventional lithium-ion batteries and provides guidelines for engineering chemomechanically resilient electrodes for high-energy batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Preliminary Thermal-hydraulics and Selected Safety Analysis for Holos-Quad Reactor Design

This report presents thermal hydraulics and safety analyses of the 10 MWe Holos-Quad micro-reactor design developed by HolosGen LLC. These analyses were executed under steady-state and station black-out (SBO) transient conditions using the System Analysis Module (SAM), with a focus on ensuring peak fuel temperatures do not exceed the safety thresholds of 1,250°C for steady state normal operating conditions and 1,600°C for transient conditions. For steady-state temperature predictions, a coupled 1D fluid-to-3D solid heat conduction methodology was used. The steady-state simulations included a model of the central core region where power peaking occurs to predict the maximum steady-state fuel temperature, and a coarse-mesh full-core model to provide the initial condition for the consequent transient SBO decay heat removal simulation. The maximum fuel temperature predicted from the central core region simulations under the normal operating condition was 1,222°C, which is below the design limit 1,250°C. To simulate a SBO transient, a full core model was necessary because of the non-axial symmetric core design and heat loss from the outer core structures to the environment. To reduce the required computational costs, a novel scaled subassembly approach was used to model the full core domain. With the full core steady-state temperature distribution provided as the initial condition, the passive decay heat removal of the core was simulated with the scaled subassembly full core model under SBO transient conditions. This analysis showed that the peak fuel temperatures of the Holos-Quad core remained below their steady-state values during the transient, therefore much lower than the safety limit 1,600°C. This is mainly due to the large thermal inertia of the graphite matrix, relatively low power density, and the large surface-to-volume-ratio of the core. The preliminary analyses presented in this report confirm the inherent safety of the Holos-Quad micro-reactor thermal-hydraulics design, as the peak fuel temperatures under both normal and off-normal operational conditions remain within design and safety limits.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Pore-Scale Simulation of Spectral Induced Polarization

A new model was developed to simulate spectral induced polarization (SIP) in porous media. The model is based on a pore-network approximation of the soil or rock, in which bonds between pores are represented by equivalent circuits. In the circuit model, electrolytic and surface conduction occur in parallel, with the electrolyte behaving as a single resistor, and the mineral surface behaving as a resistor and capacitor in series. Based on the pore geometry, fluid conductivity, and mineral surface conductivity, an impedance network is generated. SIP spectra are calculated by solving the sinusoidal steady-state electrical problem for a range of frequencies. The new model provides a tool to bridge the gap between (1) our fundamental understanding of polarization processes at the scale of a single pore or grains and (2) existing simulation capabilities for macroscale systems that do not account for pore-scale properties. The model is implemented in Matlab and can solve two-dimensional and three-dimensional (2D and 3D) networks. Examples for regular cubic lattices are presented, but the method and code are flexible with respect to pore geometries and allow for future leveraging of existing open-source capabilities for generation of pore networks and extraction of network geometries from micro imaging of core samples.

58 GEOSCIENCES↗
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↗