Abstract:The swelling of Ca-montmorillonite at elevated temperatures is important for many applications including geological disposal of radioactive waste, subsurface carbon sequestration, and shale gas exploration. However, the experimentally observed swelling behaviors of Ca-montmorillonite contacting liquid water and the temperature effects on the swelling pressure are not well understood. In this work, molecular dynamics simulations are carried out to study the swelling of Wyoming Ca-montmorillonite with a d-spacing ( d ) range of 1.40–4.00 nm at 5 MPa and various temperatures (298–500 K). The ClayFF and SPC are adopted for modeling Ca-montmorillonite and water, respectively. The simulation box is measured to be 11.15, 3.66, and 28.00 nm in the x- , y- , and z -direction. Atomistic pistons are used to control the bulk pressure of the water environment, and the implicit walls are implemented for preventing the ions from leaking from the pore into the water environment. The clay atoms are fixed during the simulation and the swelling pressure is calculated through dividing the force by the area. The equilibrium time is at least 20 ns and the production time falls in a range of 50–88 ns. The swelling pressure results show that for small d , high temperature reduces the magnitude of the oscillating curve of swelling pressure and also reduces the range of d where hydration force dominates the swelling pressure. This temperature effect is due to the weakened hydration force as evidenced from the weakened water density distributions inside the pore. For large d , high temperature reduces the swelling pressure, which is consistent with the experimental result, and increases the range of d where double layer force dominates the swelling pressure. The reduction of the swelling pressure can be explained by the enhanced ion correlation that reduces the double layer force according to the strong coupling theory, given that the calculated coupling parameters at higher temperatures are smaller. The swelling pressures are negative at elevated temperatures and large d , which prevents the clay from further swelling. However, the classical Poisson-Boltzmann (PB) equation predicts the positive double layer force since the ion correlation effect is not considered in the PB equation. Furthermore, the calculated swelling free energy curve shows that at 298 K and 5 MPa, it is difficult for Ca-montmorillonite to swell beyond a d -spacing of around 1.9 nm, which is in good agreement with the experimental result. The energy barrier for Ca-montmorillonite to swell to large d is larger than that for Na-montmorillonite, which means that it is more difficult for Ca-montmorillonite to swell to large d . This behavior is consistent with experimental observation and can be explained by the larger ion correlation effect in the Ca-montmorillonite system. These findings enhance the understanding of swelling of Ca-montmorillonite at elevated temperatures and could help to engineer better barrier materials for nuclear waste storage.

Molecular Dynamics Simulation of Hydrated Na-montmorillonite with Inorganic Salts Addition at High Temperature and High Pressure

Effects of Inorganic Salts Addition on Na-montmorillonite Clay at High Temperature and High Pressure: Insights from Molecular Dynamics Simulation

A Molecular Dynamics Study on the Structural and Mechanical Properties of Pyrophyllite and M-Montmorillonites (M = Na, K, Ca, and Ba)

Influence of inorganic and organic salts on the hydration mechanism of montmorillonite based on molecular simulation

Thermoelastic Parameter Αkt of Sodium Chloride at High Pressure and High Temperature

SWELLING BEHAVIOR OF MONTMORILLONITE, RHEOLOGICAL BEHAVIOR OF MONTMORILLONITE AQUEOUS SUSPENSION AND BLEND OF SUSPENSION AND EMULSION

A Molecular Dynamics Study on the Structural and Mechanical Properties of Hydrated Kaolinite System under Tension

Molecular Dynamics Study of Crystalline Swelling of Montmorillonite as Affected by Interlayer Cation Hydration

Early stage of swelling process of dehydrated montmorillonite through molecular dynamics simulation

Effects of Temperature on Structural Properties of Hydrated Montmorillonite: Experimental Study and Molecular Dynamics Simulation

Real-time dynamics in spin-1/2 chains with adaptive time-dependent density matrix renormalization group.

Evolution of Chemical Bonding and Crystalline Swelling-Shrinkage of Montmorillonite upon Temperature Changes Probed by in Situ Fourier Transform Infrared Spectroscopy and X-ray Diffraction

Molecular dynamics simulation of swelling properties of Ca-montmorillonite at high temperatures

Effect of Layer Charge Density on Hydration Properties of Montmorillonite: Molecular Dynamics Simulation and Experimental Study

Changes in the Interlayer Structure and Thermodynamics of Hydrated Montmorillonite under Basin Conditions: Molecular Simulation Approaches

Internal Forces Within the Layered Structure of Na-Montmorillonite Hydrates: Molecular Dynamics Simulation

Hydration and Mobility of Interlayer Ions of (Na x , Ca y )-Montmorillonite: A Molecular Dynamics Study

Nanoscale chemomechanical variations of montmorillonite induced by the specificity of counterions—An in situ XRD and AFM study

Elevated temperature effects (T > 100 °C) on the interfacial water and microstructure swelling of Na-montmorillonite

Atomic-level insights into the mechanism of saline-regulated montmorillonite (001)-salt droplet interface wetting: A molecular dynamics study