@article {250, title = {Effect of Polydispersity of Clay Platelets on the Aggregation And Mechanical Properties of Clay at the Mesoscale}, journal = {Clays and Clay Minerals}, volume = {64}, year = {2016}, month = {Aug-2016}, pages = {425 - 437}, abstract = {

The results from mesoscale simulations of the formation and evolution of microstructure for assemblies of Na-smectite particles based on assumed size distributions of individual clay platelets are presented here. The analyses predicted particle arrangements and aggregation (i.e. platelets linked in face-face configurations) and are used to link geometric properties of the microstructure and mechanical properties of the particle assemblies. Interactions between individual ellipsoidal clay platelets are represented using the Gay-Berne potential based on atomistic simulations of the free energy between two Na-smectite clay-platelets in liquid water, following a novel coarse-graining method developed previously. The current study describes the geometric (aggregate thickness, orientation, and porosity) and elastic properties in the \‘jammed states\’ from the mesoscale simulations for selected ranges of clay particle sizes and confining pressures. The thickness of clay aggregates for monodisperse assemblies increases (with average stack thickness consisting of n = 3-8 platelets) with the diameter of the individual clay platelets and with the level of confining pressure. Aggregates break down at high confining pressures (50-300 atm) due to slippage between the platelets. Polydisperse simulations generate smaller aggregates (n = 2) and show much smaller effects of confining pressure. All assemblies show increased order with confining pressure, implying more anisotropic microstructure. The mesoscale simulations are also in good agreement with macroscopic compression behavior measured in conventional 1-D laboratory compression tests. The mesoscale assemblies exhibit cubic symmetry in elastic properties. The results for larger platelets (D = 1000 \Å) are in good agreement with Nano-indentation measurements on natural clays and shale samples.

}, issn = {0009-8604}, doi = {10.1346/CCMN.2016.0640407}, author = {Ebrahimi, Davoud and Andrew J. Whittle and Roland Jean-Marc Pellenq} } @article {113, title = {Mesoscale properties of clay aggregates from potential of mean force representation of interactions between nanoplatelets}, journal = {Journal of Chemical Physics }, volume = {140}, year = {2014}, month = {Apr-21-2014}, pages = {Article Number: 154309}, type = {Article}, abstract = {

Face-to-face and edge-to-edge free energy interactions of Wyoming Na-montmorillonite platelets were studied by calculating potential of mean force along their center to center reaction coordinate using explicit solvent (i.e., water) molecular dynamics and free energy perturbation methods. Using a series of configurations, the Gay-Berne potential was parametrized and used to examine the mesoscale aggregation and properties of platelets that are initially random oriented under isothermal-isobaric conditions. Aggregates of clay were defined by geometrical analysis of face-to-face proximity of platelets with size distribution described by a log-normal function. The isotropy of the microstructure was assessed by computing a scalar order parameter. The number of platelets per aggregate and anisotropy of the microstructure both increases with platelet plan area. The system becomes more ordered and aggregate size increases with increasing pressure until maximum ordered state at confining pressure of 50 atm. Further increase of pressure slides platelets relative to each other leading to smaller aggregate size. The results show aggregate size of (3-8) platelets for sodium-smectite in agreement with experiments (3-10). The geometrical arrangement of aggregates affects mechanical properties of the system. The elastic properties of the meso-scale aggregate assembly are reported and compared with nanoindentation experiments. It is found that the elastic properties at this scale are close to the cubic systems. The elastic stiffness and anisotropy of the assembly increases with the size of the platelets and the level of external pressure. (C) 2014 AIP Publishing LLC.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
Image of FIG. 4.
}, issn = {0021-9606}, doi = {10.1063/1.4870932}, author = {Ebrahimi, Davoud and Andrew J. Whittle and Roland Jean-Marc Pellenq} } @proceedings {251, title = {Mesoscale simulation of clay aggregate formation and mechanical properties}, journal = {3rd International Symposium on Geomechanics from Micro to Macro}, volume = {Geomechanics from Micro to Macro: Granular Matter}, year = {2014}, month = {Sep-2014}, pages = {539-544}, address = {SEP 01-03 2014 Univ Cambridge, Cambridge, ENGLAND}, abstract = {

This paper proposes a novel methodology for understanding the meso-scale aggregation of clay platelets in water. We use Molecular Dynamics simulations using the CLAYFF force fields to represent the interactions between two layers of Wyoming montmorillonite (Na-smectite) in bulk water. The analyses are used to establish the potential of mean force at different spacings between the layers for edge-to-edge and face-to-face interactions. This is accomplished by finding the change in free energy as a function of the separation distance between the platelets using thermodynamic perturbation theory with a simple overlap sampling method. These nanoscale results are then used to calibrate the Gay\–Berne (GB) potential that represents each platelet as a single-site ellipsoidal body. A coarse-graining upscaling approach then uses the GB potentials and molecular dynamics to represent the meso-scale aggregation of clay platelets (at submicron length scale). Results from meso-scale simulations obtain the equilibrium/jamming configurations for mono-disperse clay platelets. The results show aggregation for a range of clay platelets dimensions and pressures with mean stack size ranging from 3 to 8 platelets. The particle assemblies become more ordered and exhibit more pronounced elastic anisotropy at higher confining pressures. The results are in good agreement with previously measured nano-indentation moduli over a wide range of clay packing densities.

}, issn = {1434-5021}, doi = {10.1007/s10035-016-0655-8}, author = {Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Andrew J. Whittle}, editor = {Kenichi Soga and Krishna Kumar and Giovanna Biscontin} } @article {47, title = {Nanoscale Elastic Properties of Montmorillonite upon Water Adsorption}, journal = {Langmuir}, volume = {28}, year = {2012}, month = {Dec-11-2012}, pages = {16855{\textendash}16863}, abstract = {
Abstract Image

Smectites are an important group of clay minerals that experience swelling upon water adsorption. This paper uses molecular dynamics with the CLAYFF force field to simulate isothermal isobaric water adsorption of interlayer Wyoming Na-montmorillonite, a member of the smectite group. Nanoscale elastic properties of the clay\–interlayer water system are calculated from the potential energy of the model system. The transverse isotropic symmetry of the elastic constant matrix was assessed by calculating Euclidean and Riemannian distance metrics. Simulated elastic constants of the clay mineral are compared with available results from acoustic and nanoindentation measurements.

}, doi = {10.1021/la302997g}, author = {Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Andrew J. Whittle} } @proceedings {364, title = {Springer Series in Geomechanics and Geoengineering Multiscale and Multiphysics Processes in Geomechanics SIMULATION OF HYDRATION AND ELASTIC PROPERTIES OF MONTMORILLONITE USING MOLECULAR DYNAMICS}, journal = {International Workshop on Multiscale and Multiphysics Processes in Geomechanics}, volume = {Springer Series in Geomechanics and Geoengineering}, year = {2011}, month = {2011}, pages = {105 - 108}, publisher = {Springer Berlin Heidelberg}, address = {JUN 23-25 2010 Stanford Univ, Stanford, CA}, abstract = {

This article describes molecular-level simulations of the adsorption isotherm and elastic properties for a Wyoming Na-montmorillonite using the General Utility Lattice Program. The selected clay mineral includes isomorphous substitutions in both the silica tetrahedral and octahedral gibbsite sheets. The resulting negative charges are balanced by Na+ cations. The authors have investigated the importance of different interatomic potentials in understanding the adsorption of water molecules (at ambient temperature, 300K). A unique combination of NPT and Grand Canonical ensemble simulations using Core\–Shell potential shows intra-layer adsorption, while the same approach based on empirical CLAYFF potential produces only interlayer water adsorption. There is a correspondingly large difference in the theoretical elastic modulus normal to the clay sheets using the two different sets of force potentials. This article discusses the underlying reasons for these differences and reviews the experimental evidence supporting intra-layer water absorption in montmorillonite.

}, isbn = {978-3-642-19629-4}, issn = {1866-8755}, doi = {10.1007/978-3-642-19630-0_27}, author = {Roland Jean-Marc Pellenq and Andrew J. Whittle and Ebrahimi, Davoud}, editor = {Borja, Ronaldo I. and Wu, Wei} }