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: Molecular dynamics modeling of a partially saturated clay‐water system at finite temperature
Summary The mechanical and hydraulic properties of unsaturated clay under nonisothermal conditions have practical implications in geotechnical engineering applications such as geothermal energy harvest, landfill cover design, and nuclear waste disposal facilities. The water menisci among clay particles impact the mechanical and hydraulic properties of unsaturated clay. Molecular dynamics (MD) modeling has been proven to be an effective method in investigating clay structures and their hydromechanical behavior at the atomic scale. In this study, we examine the impact of temperature increase on the capillary force and capillary pressure of the partially saturated clay‐water system through high‐performance computing. The water meniscus formed between two parallel clay particles is studied via a full‐scale MD modeling at different elevated temperatures. The numerical results have shown that the temperature increase impacts the capillary force, capillary pressure, and contact angle at the atomic scale. The capillary force on the clay particle obtained from MD simulations is also compared with the results from the macroscopic theory. The full‐scale MD simulation of the partially saturated clay‐water system can not only provide a fundamental understanding of the impact of temperature on the interface physics of such system at the atomic scale, but also has practical implication in formulating physics‐based multiscale models for unsaturated soils by providing interface physical properties of such materials directly through high‐performance computing.  more » « less
Award ID(s):
1659932
PAR ID:
10453131
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical and Analytical Methods in Geomechanics
Volume:
43
Issue:
13
ISSN:
0363-9061
Page Range / eLocation ID:
p. 2129-2146
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The engineering problems involving clay under non‐isothermal conditions (e.g., geothermal energy harvest, landfill cover system, and nuclear waste disposal) are multiscale and multiphysics by nature. The nanoscale hydrodynamics of clay at elevated temperature is essential in developing a physics‐based multiscale model for clay under non‐isothermal conditions. The nonequilibrium molecular dynamics (NEMD) is a useful tool to study the nanoscale hydrodyndamics of clay. This article presents an NEMD modeling of hydrodynamics of clay nanopores at elevated temperatures. Water flow confined in pyrophyllite and montmorillonite clay nanopores is investigated. The nonequilibrium state is maintained by uniformly exerting an external force on each water molecule. The NEMD simulations have provided a molecular‐scale perspective of temperature effect on clay‐water density, water flow velocity, shear viscosity, clay‐water slip length, hydraulic conductivity, and clay‐water friction coefficient. The numerical results have shown a strong temperature dependence of fluid flow velocity, shear viscosity, clay‐water slip length, and hydraulic conductivity at the nanoscale. We have validated the applicability of cubic law in determining hydraulic conductivity at the nanopore scale at elevated temperatures. It is found from our numerical results that slip clay‐water boundary condition is an essential factor in properly determining nanoscale fluid flow velocity. By numerical examples, we also study the impact of nanopore size and clay layer thickness on the hydrodynamics of the clay‐water system. 
    more » « less
  2. Bardanis, M. (Ed.)
    This paper characterizes nanoscale soil-water retention mechanism of unsaturated clay through molecular dynamics simulation. Series of molecular dynamics simulations of clay at low degrees of saturation were conducted. Soil water was represented by a point cloud through the centre-of-massmethod. Water-air interface area was measured numerically by the alpha shape method. Spatial variation of water number density is characterized and used to determine the adsorbed water layer. The soil-water retention mechanism at the nanoscale was analysed by distinguishing adsorptive pressure and capillary pressure at different mass water contents and considering apparent interface area (water-air interface area per unit water volume). 
    more » « less
  3. 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
  4. In order to fully understand the thermo-hydro-mechanical behavior of the geotechnical infrastructures, the effects of temperature variations on soil properties and soil behavior have to be studied. Hydraulic conductivity, strength, volume change, moisture content, and pore pressure generation and dissipation rates depend on temperature variations. Thermal loading might induce excess pore water pressure and volumetric changes. Temperature changes in the fine-grained soils will cause expansion in water and soil particles. Since the coefficient of expansion for soil particles is much smaller than that for water, a generation of pore water pressure is expected. This thermally induced pore water pressure and then its dissipation during the relaxation period results in a time dependent consolidation. Thermal consolidation in fine grained soil is more dominant and can be irreversible in normally consolidated clay. However, the volumetric changes of highly over consolidated soil caused by temperature increment is reversible by temperature reduction. In this research, a modified consolidation testing device is used to study the effect of temperature increments (e.g., increasing step by step temperature increments to 80ºC) on the consolidation of fine grained soils. In another words the effect of temperature increments during the test on the consolidation process is studied. Time of applying the heating, target temperature, and initial void ratio are parameters affecting the rate and the amount of consolidation in the samples. 
    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