skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Comparison of thermo-poroelastic and thermo-poroelastoplastic constitutive models to analyze THM process in clays
Thermal pore pressurization in soil media has been investigated for the past few decades. It has been shown that temperature variations may significantly affect thermal pore pressure in clay soils confined deep into the ground. Moreover, thermal loading may lead to stress change and thermal deformation. Thermo-poroelastic and advance thermo-poroelastoplastic constitutive models have been formulated and incorporated numerically to simulate the thermo-hydro-mechanical process. However, the accurate response of soil media during THM process has not been completely understood. Although numerical modelling reasonably predicts the experimental observations, they still could not be used to completely justify the field observations. In this study, the main features of the thermo-poroelastic model are incorporated in a thermo-poroelastoplastic constitutive model (ACMEG-T) to further investigate the effect of different thermal and hydraulic properties on thermo-hydro-mechanical (THM) response of the soil media.  more » « less
Award ID(s):
1804822
PAR ID:
10206115
Author(s) / Creator(s):
;
Editor(s):
McCartney, J.S.; Tomac, I.
Date Published:
Journal Name:
E3S Web of Conferences
Volume:
205
ISSN:
2267-1242
Page Range / eLocation ID:
04008
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Temperature variations in low permeable soil (e.g. clay) induce pore pressure, which is known as thermal pressurization. Previous research showed that thermal pressurization highly depends on thermal pressurization coefficient. This coefficient depends on the soil type and changes with temperature due to temperature dependency of thermal expansion coefficient of water. Thermal pressurization is often investigated through thermo-hydro-mechanical (THM) numerical modeling. THM process, with respect to thermal loading, has been examined in the literature to justify the field observations by incorporating advanced thermo-mechanical constitutive models. However, result of numerical simulations using advanced thermo-elastoplastic models still show some discrepancies with experimental and field observations. In this study, the assessment of thermal pressurization in Boom clay is scrutinized through employing a relatively simple while practical thermo-poroelastic finite element model with careful consideration of the temperature-dependent thermal, hydraulic, and mechanical properties of the medium and saturating fluid (i.e. water). The numerical model is carried out using COMSOL Multiphysics and the results of the numerical simulations are compared and validated with the ATLAS project, a large-scale experimental facility in Belgium. The results confirm that thermal and hydraulic coupling parameters are the key factors to change thermal pressurization. 
    more » « less
  2. The mechanical behavior of unsaturated porous media under non-isothermal conditions plays a vital role in geo-hazards and geo-energy engineering (e.g., landslides triggered by fire and geothermal energy harvest and foundations). Temperature increase can trigger localized failure and cracking in unsaturated porous media. This article investigates the shear banding and cracking in unsaturated porous media under non-isothermal conditions through a thermo–hydro–mechanical (THM) periporomechanics (PPM) paradigm. PPM is a nonlocal formulation of classical poromechanics using integral equations, which is robust in simulating continuous and discontinuous deformation in porous media. As a new contribution, we formulate a nonlocal THM constitutive model for unsaturated porous media in the PPM paradigm in this study. The THM meshfree paradigm is implemented through an explicit Lagrangian meshfree algorithm. The return mapping algorithm is used to implement the nonlocal THM constitutive model numerically. Numerical examples are presented to assess the capability of the proposed THM mesh-free paradigm for modeling shear banding and cracking in unsaturated porous media under non-isothermal conditions. The numerical results are examined to study the effect of temperature variations on the formation of shear banding and cracking in unsaturated porous media. 
    more » « less
  3. Abstract. Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need to reexamine the assumption of predetermined frictional breakdown energy common in dynamic rupture modeling and to better understand the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes. 
    more » « less
  4. Summary This study presents a thermo‐hydro‐mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA‐FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker‐Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi‐static loading conditions. Temperature dependency of the surface tension is implemented in the soil‐water retention curve. Nonuniform rational B‐splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA‐based framework offers a viable alternative for solving THM problems in unsaturated soils. 
    more » « less
  5. Climate change is expected to alter the statistics of extreme events including rainfall storms, floods, droughts, and heatwaves. Climate-adaptive geotechnical structures warrant a quantitative assessment of the impacts of emerging and projected extreme patterns on the short and long-term behaviors of earthen structures. Furthermore, long-term changes to soil carbon and moisture due to non-extreme climate events should also be considered. While several large-scale studies have been conducted to evaluate various aspects of climate change, there is a clear gap in the state of knowledge in terms of assessing the resilience of geotechnical structures to changes in climatic trends (e.g., warmer climate, protracted droughts, intensified extreme precipitations, and sea level rise). The majority of the aforementioned climatic trends pose multi-physics problems involving thermo-hydro-mechanical (THM) processes in partially saturated soils and earthen structures. This review paper discusses how soil-atmospheric interactions and extreme event patterns in a changing climate can alter soil properties and loading conditions, affecting the performance of partially saturated geotechnical structures. We speculate how changes in climatic trends may weaken partially saturated earthen structures through strength reduction, drying, soil desiccation cracking, shrinkage, microbial oxidation of soil organic matter, fluctuation in the ground water table, land and surface erosion, and highly dynamic pore pressure changes. Each of these weakening processes is primarily induced by variations in the soil moisture and temperature. Finally, we discuss potential modes of failure imposed on partially saturated earthen structures by climatic trends. 
    more » « less