Results for “thermal-hydraulic-mechanical coupling”

DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

DOE OSTI2026

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗

DOE OSTI2020

Prediction of the thermal-hydraulic-mechanical response of a geological repository at large scale and sensitivity analyses

To better understand the Thermal-Hydraulic-Mechanical (THM) response of the rock caused by the high-level radioactive waste (HLW) released heat in the Callovo-Oxfordian formation (COx) in the near-field and far-field areas, a series of coupled THM modelling is performed. This study presents a case study based on the French geological repository for HLW (Cigéo project) to better assess modelling of deep geological repositories (DGR) within the DECOVALEX-2019 framework. In this study, the proposed coupled THM model, implemented in the platform of COMSOL, is validated by comparing the modelled THM results for a point heat source in an infinite rock mass with an analytical solution. It is also validated by comparing its thermal component for the Base Case with accurate results calculated using another numerical method. To shorten the calculation time, the DGR is modelled using six cuboid blocks representing 162 HLW cells and only six central HLW cells of interest modelled in detail. This simplification is validated. This model can easily be used to model not only the near-field THM response of the rock but also the far-field THM response. The influence of the boundary conditions applied on the gallery wall, HLW cell walls, and external surfaces of the model and the influence of different vertical dimensions of the geometry are studied. The sensitivity analyses of the THM parameters of COx on the THM response at different locations are performed. Finally, the influence of using a 2-dimensional (2D) model to represent a 3-dimensional (3D) repository is also investigated.

58 GEOSCIENCES↗

DOE OSTI2025

An international study on THM modelling of the full-scale heater experiment at Mont Terri laboratory

We present results from an international model comparison study of the Full-Scale Emplacement (FE) experiment in Opalinus Clay at the Mont Terri Laboratory, Switzerland. Based on a provided parameter set the teams decided which parameters they adopted for their models, whether they considered the excavation and the ventilation phase in addition to the heating phase and if they included technical features like the shotcrete or the EDZ. The teams were able to reproduce the measured parameters temperature, relative humidity and pore pressure. The modelled results for temperature agree very closely between the teams especially in the sensors in Opalinus Clay. All teams were able to reproduce the redistribution of water in the bentonite backfill due to heating. The evolution of the relative humidity showed similar trends with differences in the intensity of the dry out effect. To model the pore pressure evolution is more complex because it comprises the full interaction of the coupled THM processes. The spread between the pore pressure modelled by the teams was larger, with some teams overestimating the pressure increase due to heating and some teams overestimating the extent of drainage. The agreement of modelled results with measurements improves with larger distance to the heater. We conclude that the EDZ and the shotcrete potentially influence the behaviour of the rock causing higher differences closer to the heater. Further research is needed to better implement those influences into the models. Based on the calibrated models, the future evolution of temperature, relative humidity and pore pressure was predicted over the next 10 years following a change of the heat power applied in 2023 and 2024. Again, the predicted temperatures agree very closely between the teams. Most teams do not expect an increase in relative humidity during the next 10 years after the initial dry-out.

58 GEOSCIENCES↗

DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2024

KOMPASS-II: Compaction of Crushed salt for Safe Containment – Phase 2

Long-term stable sealing elements are a basic component in the safety concept for a possible repository for heat-emitting radioactive waste in rock salt. The sealing elements will be part of the closure concept for drifts and shafts. They will be made from a welldefinied crushed salt in employ a specific manufacturing process. The use of crushed salt as geotechnical barrier as required by the German Site Selection Act from 2017 /STA 17/ represents a paradigm change in the safety function of crushed salt, since this material was formerly only considered as stabilizing backfill for the host rock. The demonstration of the long-term stability and impermeability of crushed salt is crucial for its use as a geotechnical barrier. The KOMPASS-II project, is a follow-up of the KOMPASS-I project and continues the work with focus on improving the understanding of the thermal-hydraulic-mechanical (THM) coupled processes in crushed salt compaction with the objective to enhance the scientific competence for using crushed salt for the long-term isolation of high-level nuclear waste within rock salt repositories. The project strives for an adequate characterization of the compaction process and the essential influencing parameters, as well as a robust and reliable long-term prognosis using validated constitutive models. For this purpose, experimental studies on long-term compaction tests are combined with microstructural investigations and numerical modeling. The long-term compaction tests in this project focused on the effect of mean stress, deviatoric stress and temperature on the compaction behavior of crushed salt. A laboratory benchmark was performed identifying a variability in compaction behavior. Microstructural investigations were executed with the objective to characterize the influence of pre-compaction procedure, humidity content and grain size/grain size distribution on the overall compaction process of crushed salt with respect to the deformation mechanisms. The created database was used for benchmark calculations aiming for improvement and optimization of a large number of constitutive models available for crushed salt. The models were calibrated, and the improvement process was made visible applying the virtual demonstrator.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗