SEARCH · Engineering Papers

Results for “modeling framework”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Dual-module humidity pump for efficient air dehumidification: Demonstration and performance limitations

Condensation dehumidification in conventional air conditioning technologies is energy-intensive, accounting for up to 50% of building cooling energy used in some climates. Selective vacuum membrane dehumidification (VMD) is one of the leading alternative dehumidification technologies due to its potential for significant energy savings, and the “dual-module humidity pump” is one of the most promising VMD concepts. Here, this work is the first to provide experimental proof-of-concept for the dual-module humidity pump system and provides the first thermodynamic modeling framework that accounts for realistic steady-state operating limitations, both of which are lacking in the current literature. Additionally, this work is the first to provide a system design solution that overcomes practical challenges associated with air accumulation in the vacuum channels. The experimental results in this work show that the current prototype can remove up to 45% of the water vapor in the air stream, and the vapor pressure difference in the vapor rejection module needs to be approximately 2–4 times greater than that of the dehumidification module in order to maintain balanced mass transfer. The thermodynamic model applied to typical air conditioning conditions shows that the ideal dehumidification (latent) COPs can reach up to 40, but practical COPs are limited to approximately 10. Furthermore, the model shows that the overall energy efficiency increases as the membrane air selectivity increases, though this improvement gradually starts to diminish when the membrane selectivity is increased above 10,000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Passive Energy-Saving Solutions for Clothes Dryers: A Modeling and Experimental Study for Improved Efficiency and Affordability

Dryers are integral appliances in modern households, yet their significant power consumption remains a critical challenge for reducing energy bills and upgrades. This work presents the outcomes of a comprehensive investigation aimed at lowering the operating energy costs of the dryer through modeling and experimental approaches. Specifically, the objectives of this research are twofold: (1) to reduce energy consumption without negatively impacting drying performance or time; and (2) to ensure affordability by developing retrofittable solutions characterized by low cost and a quick payback period. These strategies promise universal applicability to all dryer categories, encompassing both gas and electric models, by tackling core issues such as unnecessary heat loss and excess energy supply, both prevalent across dryer designs. Our methodology combines robust modeling frameworks and experimental validation to target energy efficiency improvements through three primary pathways: (1) effective heat loss management using ultralow-cost insulation materials tailored for dryer systems; (2) heat management across the drying cycle, facilitated by passive heat transfer mechanisms, and (3) optimization of heat supply to the load to ensure precise energy delivery, minimizing waste. These innovations are designed to seamlessly integrate with existing dryer configurations, providing a scalable and retrofittable solution that ensures affordability without requiring substantial redesigns or expensive components. Using these methods, the study demonstrates improvement opportunities in dryer energy consumption while maintaining the desired performance. The modeling component utilizes computational simulations to evaluate the thermal and energy performance of these innovations under varied operational conditions, providing foundational insights for experimental design. By addressing critical areas such as heat loss, energy recovery, and supply optimization, this work proposes impactful solutions for lowering the energy costs in domestic clothes dryers. The ultralow-cost, retrofittable nature of the proposed strategies ensures widespread adoption potential across diverse dryer categories, making significant strides toward affordable household energy practices.

Cheekatamarla, Praveen [ORNL] (ORCID:0000000248827↗

DOE OSTI2020

The Role of Interfacial Thermal Resistance in Li-Ion Battery Thermal Management: Preprint

Temperature critically affects the performance, life and safety of lithium-ion batteries. Therefore, it is essential to understand heat generation and dissipation within individual battery cells and battery packs to plan a proper thermal management strategy. One of the key challenges is that interfacial heat transfer of a battery unit is difficult to quantify. The steady-state absolute method and the transient laser-flash-diffusivity method were employed to measure heat conductivities of battery layer stacks and individual battery layer separately. Results show flash diffusivity method gives higher thermal conductivity at both cross-plane and in-plane directions. The difference is primarily caused by interfacial thermal resistance so that it can be estimated by steady-state and transient measurements. To investigate the effects of interfacial thermal transport beyond individual cell level, a multiphysics battery model is used. The model is built upon a multi-scale multi-domain modeling framework for battery packs that accounts for the interplay across multiple physical phenomena. Benefits of a battery module using thermal management materials are quantified through numerical experiments. During a thermal runaway event, it is found interfacial thermal resistance can mitigate thermal runaway in a battery module by significantly reducing heat transfer between cells.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Self-consistent investigation of density fueling needs on ITER and CFETR utilizing the new Pellet Ablation Module

Abstract Self-consistent modeling using the stability, transport, equilibrium, and pedestal (STEP) workflow in the OMFIT integrated modeling framework (predicting pedestal with EPED, core profiles with TGYRO, current profile with ONETWO, and EFIT for equilibrium) suggests ITER and future devices such as China Fusion Engineering Test Reactor (CFETR) Zhuang et al (2019 Nucl. Fusion 59 112010) will benefit from high-density operation (Greenwald limit fraction f g w ≈ 0.7−1.3). Regimes with an operational density near the Greenwald limit will likely need peaked density profiles so that the pedestal density remains below the Greenwald limit. Peaked density profiles can be achieved with the help of pellet injection. A flexible Pellet Ablation Module (PAM), which predicts the density source based on a comprehensive analytical pellet ablation model, has been developed for predicting pellet fueling for transport studies, and has been incorporated into the STEP workflow for predictive modeling. This workflow is applied to DIII-D and finds good agreement with experiments. On ITER the effect of pellet fueling is examined in an advanced inductive scenario, where a fusion gain of up to Q = 9 is predicted with strong central pellet fueling. On CFETR, with a mid-radius density source, an average of 1.5 × 10 22 electrons s −1 are required to achieve the density and temperature profiles necessary for the 1000 MW advanced scenario with a tritium burn-up fraction of ∼ 3 % .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗

DOE OSTI2025

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Using probability distribution function as a scaling approach to incorporate soil heterogeneity into biogeochemical models for greenhouse gas predictions (Final Technical Report)

The project investigated biogeochemical processes at terrestrial-aquatic interfaces (TAIs), focusing on soil microsite heterogeneity and its impact on greenhouse gas (GHG) fluxes. Using laboratory experiments, modeling, and data integration, researchers explored redox-driven microbial processes under fluctuating hydrological conditions. Key advancements included modifying the DAMM-GHG model to incorporateelectron acceptor availability and enhancing the AquaMEND model for improved microbial metabolism representation. Results highlighted microsite redox variability as a key driver of GHG fluxes, informing Earth system models. The project fostered interdisciplinary collaborations, student training, and the development of novel modeling frameworks to improve Earth'senergy budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Critically assessing sodium-ion technology roadmaps and scenarios for techno-economic competitiveness against lithium-ion batteries

Sodium-ion batteries have garnered notable attention as a potentially low-cost alternative to lithium-ion batteries, which have experienced supply shortages and price volatility for key minerals. Here we assess their techno-economic competitiveness against incumbent lithium-ion batteries using a modelling framework incorporating componential learning curves constrained by minerals prices and engineering design floors. We compare projected sodium-ion and lithium-ion price trends across over 6,000 scenarios while varying Na-ion technology development roadmaps, supply chain scenarios, market penetration and learning rates. Assuming that substantial progress can be made along technology roadmaps via targeted research and development, we identify several sodium-ion pathways that might reach cost-competitiveness with low-cost lithium-ion variants in the 2030s. In addition, we show that timelines are highly sensitive to movements in critical minerals supply chains—namely that of lithium, graphite and nickel. Our modelled outcomes suggest that being price advantageous against low-cost lithium-ion variants in the near term is challenging and increasing sodium-ion energy densities to decrease materials intensity is among the most impactful ways to improve competitiveness.

25 ENERGY STORAGE↗

DOE OSTI2025

A DECADE of dwarfs: first detection of weak lensing around spectroscopically confirmed low-mass galaxies

We present the first detection of weak gravitational lensing around spectroscopically confirmed dwarf galaxies, using the large overlap between DESI DR1 spectroscopic data and DECADE/DES weak lensing catalogs. A clean dwarf galaxy sample with well-defined redshift and stellar mass cuts enables excess surface mass density measurements in two stellar mass bins ($\log \rm{M}_*=[8.2, 9.2]~M_\odot$ and $\log \rm{M}_*=[9.2, 10.2]~M_\odot$), with signal-to-noise ratios of $5.6$ and $12.4$ respectively. This signal-to-noise drops to $4.5$ and $9.2$ respectively for measurements without applying individual inverse probability (IIP) weights, which mitigates fiber incompleteness from DESI's targeting. The measurements are robust against variations in stellar mass estimates, photometric shredding, and lensing calibration systematics. Using a simulation-based modeling framework with stellar mass function priors, we constrain the stellar mass-halo mass relation and find a satellite fraction of $\simeq 0.3$, which is higher than previous photometric studies but $1.5σ$ lower than $Λ$CDM predictions. We find that IIP weights have a significant impact on lensing measurements and can change the inferred $f_{\rm{sat}}$ by a factor of two, highlighting the need for accurate fiber incompleteness corrections for dwarf galaxy samples. Our results open a new observational window into the galaxy-halo connection at low masses, showing that future massively multiplexed spectroscopic observations and weak lensing data will enable stringent tests of galaxy formation models and $Λ$CDM predictions.

To, Chun-Hao [Chicago U., Astron. Astrophys. Ctr.;↗

DOE OSTI2025

Interacting mesons as degrees of freedom in a chiral model

Here, we study the equation of state of hot and dense hadronic matter using an extended chiral mean field (CMF) model framework where the addition is the inclusion of interactions of thermally excited mesons. This is implemented by calculating the in-medium masses of pseudoscalar and vector mesons, obtained through the explicit chiral symmetry-breaking and vector-interaction terms in the Lagrangian, respectively, prior to applying the mean-field approximation. As a result, the in-medium meson contributions generate a feedback term to the CMF’s equations of motion, which then modifies the equation of state. With this improvement, we quantify the effect on the equation of state of strongly interacting matter through comparisons with state-of-the-art lattice QCD results and other hadronic models like the hadron resonance gas model. We find that the results of the updated hadronic CMF model with an improved meson description (mCMF) provide a better agreement with lattice-QCD data for thermodynamic state variables across a wide range of temperatures and baryon chemical potentials.

phenomenology↗

DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗

DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗

DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗

DOE OSTI2026

Stochastic Thermo-Hydro Modeling and Neural Network Surrogate Development for Thermal Resource Assessment of the Galleries-to-Calories Geobattery

The Galleries-to-Calories Geobattery concept explores the use of abandoned coal mine workings for large-scale thermal energy transport and storage. The system involves injecting waste heat from a supercomputing facility into flooded mine galleries, where groundwater flow can store and transport thermal energy for potential recovery in downgradient district heating and cooling applications. To evaluate the feasibility and performance of the Geobattery under geological and operational uncertainty, we developed a suite of stochastic thermo-hydrological (TH) simulations using Monte Carlo sampling of key uncertain parameters (e.g., permeability, porosity, thermal conductivity, specific heat capacity) and operating conditions (e.g., injection rate, injection temperature). Results identified injection rate and temperature as the most influential parameters governing thermal front propagation, while the geometry of the room-and-pillar structure played a critical role in directing the extent and orientation of thermal advancement. Optimal combinations of material properties for maximizing heat recovery were also determined. To address the high computational cost of coupled-process stochastic modeling, we trained a neural network surrogate model on 24,000 physics-based realizations, achieving an R² > 0.99 and MAE < 0.1 for temperature predictions at monitoring locations. This surrogate enabled an additional 100,000 realizations for global sensitivity analysis and probabilistic thermal resource assessment. The integrated stochastic physics–surrogate modeling framework offers a computationally efficient tool for quantifying uncertainty, identifying key drivers, and informing early-stage design decisions for Geobattery systems.

15 - GEOTHERMAL ENERGY↗

DOE OSTI2023

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗