A fully coupled thermo-poro-mechanical finite element analysis to predict the thermal pressurization and thermally induced pore fluid flow in soil media

Tamizdoust, Mohammadreza Mir; Ghasemi‐Fare, Omid

doi:10.1016/j.compgeo.2019.103250

Cited by 59 publications

(10 citation statements)

References 38 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…TAMIZDOUST AND GHASEMI-FARE All the temperature-dependent and saturation-dependent parameters used in Equations 1 to 13 are presented in Table 1. The rest of parameters are calculated from the formulations provided by the International Association for the Properties of Water and Steam (IAPWS 2007) (Tamizdoust & Ghasemi-Fare, 2020;Wagner & Kretzschmar, 2008). In Table 1, the thermal conductivity of the medium depends on the volumetric liquid content, temperature, and clay fraction on the medium according to Bittelli et al (2015), and it has been incorporated in different numerical models to predict the nonisothermal evaporation process (Bittelli et al, 2008;Kroener et al, 2014;Trautz et al, 2015).…”

Section: 1029/2020wr027381mentioning

confidence: 99%

Utilization of Nonequilibrium Phase Change Approach to Analyze the Nonisothermal Multiphase Flow in Shallow Subsurface Soils

Tamizdoust

Ghasemi‐Fare

2020

Water Resources Research

Self Cite

View full text Add to dashboard Cite

The prediction of coupled nonisothermal multiphase flow in porous media has been the subject of many theoretical and experimental studies in the past half a century. In particular, the evaporation phenomenon from the shallow subsurface has been extensively studied based on the notion of equilibrium phase change between liquid water and water vapor (i.e., instantaneous phase change). One of the frequent assumptions in equilibrium phase change approach is that liquid water is hydraulically connected throughout the vadose zone. Furthermore, classical soil-water retention curves (e.g., van Genuchten model), which have been extensively used in the literature to model evaporation process, are only valid for high and intermediate saturation degrees. Although these limitations have been addressed and improved in separate studies, they have not yet been rigorously incorporated in the numerical modeling of nonisothermal multiphase flow in shallow subsurface of in-field soils. Therefore, the aim of this study is to investigate the coupled heat, liquid, and vapor flow in soil media through the Hertz-Knudsen-Schrage (HKS) phase change model and by incorporating a water retention model which captures the soil-water characteristics from full to oven-dried saturation degrees. A numerical model is developed and validated against the in-field experimental data. Reasonable agreements between the calculated and measured values of water contents at all depths, as well as the temperature, and cumulative evaporation are observed. Results also confirm that the contribution of the film flow in overall mass flow in the medium is required for accurate modeling and cannot be ignored.

show abstract

Section: 1029/2020wr027381mentioning

confidence: 99%

Utilization of Nonequilibrium Phase Change Approach to Analyze the Nonisothermal Multiphase Flow in Shallow Subsurface Soils

Tamizdoust

Ghasemi‐Fare

2020

Water Resources Research

Self Cite

View full text Add to dashboard Cite

show abstract

“…When the pore pressure and temperature gradients exist in the pore medium, the pore flow has heat exchange with the solid skeleton, which leads to the redistribution of temperature, seepage, and stress fields. This phenomenon is known as the "Dufour flow" and is widely used in shallow geothermal systems [3], enhanced geothermal systems (EGS) [4], and water conservancy projects [5].…”

Section: Introductionmentioning

confidence: 99%

Effect of Impounding on Bedrock Temperature of High Arch Dam Site: A Case Study of the Xiluodu Project, China

Zhang

Ren

Shen

et al. 2023

Water

View full text Add to dashboard Cite

The 300 m-level high arch dams built in China have different degrees of valley-narrowing deformation in the initial impoundment. This unexpected phenomenon may increase the long-term safety risk of high arch dams. There are several explanations for the valley-narrowing deformation phenomenon, and one of the explanations—that it is caused by a drop in bedrock temperature—is still under debate. In this paper, coupled thermo-hydromechanical (THM) simulations of the temporal and spatial evolutions of the bedrock temperature and seepage fields of the Xiluodu arch dam site were performed. The results showed that the mechanism of the bedrock temperature drop after impoundment can be divided into two types: short-term type and long-term. For the short-term type, the maximum temperature drop of 9 °C is caused by the heat exchange between the bedrock surface layer and the reservoir water during the initial impoundment. For the long-term type, the maximum temperature drop of 15 °C in front of the grouted curtain is caused by the long-term infiltration of cold water from upstream to downstream bypassing the grout curtain. The monitoring data showed that the valley-narrowing deformation mainly occurred in the initial impoundment, which corresponds to the short-term type. Notably, the valley-narrowing deformation may be overestimated if the long-term type of bedrock temperature drop is considered. Taken together, these findings demonstrate that the bedrock temperature change due to impoundment is dominated by the short-term type and can be applied to study the effect of temperature on valley-narrowing deformation after impounding.

show abstract

“…FEniCS is a popular and open‐source computing framework that enables automated solution of PDEs with great flexibility and efficiency by a collection of FE‐based libraries 9 . COMSOL Multiphysics is a commercial cross‐platform finite element software and efficient solver for coupled PDEs, which has been successfully applied in THM(C) topics such as heat transfer in saturated soil, 10 methane hydrate 11 and CO 2 injection 12 . Apart from models that are solely based on FEM, the Los Alamos National Laboratory (LANL) has developed the Finite Element Heat and Mass Transfer (FEHM) code where the finite volume method (FVM) is utilized for flow and mass balance while FEM for stress equilibrium 13 .…”

Section: Introductionmentioning

confidence: 99%

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

Zhou

Zhang

2023

Num Anal Meth Geomechanics

View full text Add to dashboard Cite

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanicalchemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. 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.

show abstract

A fully coupled thermo-poro-mechanical finite element analysis to predict the thermal pressurization and thermally induced pore fluid flow in soil media

Cited by 59 publications

References 38 publications

Utilization of Nonequilibrium Phase Change Approach to Analyze the Nonisothermal Multiphase Flow in Shallow Subsurface Soils

Utilization of Nonequilibrium Phase Change Approach to Analyze the Nonisothermal Multiphase Flow in Shallow Subsurface Soils

Effect of Impounding on Bedrock Temperature of High Arch Dam Site: A Case Study of the Xiluodu Project, China

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

Contact Info

Product

Resources

About