Mesoscale Poroelasticity of Heterogeneous Media

The poromechanics of heterogeneous media is reformulated in a discrete framework using the lattice element method (LEM) that accounts for the presence of interfaces as well as local microtextural and elastic variations. The exchange of mechanical information between pore and solid(s) is captured by means of force field potentials for these domains, which eliminate the requirement of scale separability of continuum-based poromechanics approaches. In congruence with μVT and NPT ensembles of statistical mechanics, discrete expressions for Biot poroelastic coefficients are derived. Considering harmonic-type interaction potentials for each link, analytical expressions for both isotropic and transversely isotropic effective elasticity are presented. The theory is validated against continuum-based expressions of Biot poroelastic coefficients for porous media with isotropic and transversely isotropic elastic solid behavior

Capturing material toughness by molecular simulation: accounting for large yielding effects and limits

The inherent computational cost of molecular simulations limits their use to the study of nanometric systems with potentially strong size effects. In the case of fracture mechanics, size effects due to yielding at the crack tip can affect strongly the mechanical response of small systems. In this paper we consider two examples: a silica crystal for which yielding is limited to a few atoms at the crack tip, and a nanoporous polymer for which the process zone is about one order of magnitude larger. We perform molecular simulations of fracture of those materials and investigate in particular the system and crack size effects. The simulated systems are periodic with an initial crack. Quasi-static loading is achieved by increasing the system size in the direction orthogonal to the crack while maintaining a constant temperature. As expected, the behaviors of the two materials are significantly different. We show that the behavior of the silica crystal is reasonably well described by the classical framework of linear elastic fracture mechanics (LEFM). Therefore, one can easily upscale engineering fracture properties from molecular simulation results. In contrast, LEFM fails capturing the behavior of the polymer and we propose an alternative analysis based on cohesive crack zone models. We show that with a linear decreasing cohesive law, this alternative approach captures well the behavior of the polymer. Using this cohesive law, one can anticipate the mechanical behavior at larger scale and assess engineering fracture properties. Thus, despite the large yielding of the polymer at the scale of the molecular simulation, the cohesive zone analysis offers a proper upscaling methodology.MIT Energy InitiativeShell Oil CompanySchlumberger Limite

Fracture toughness of calcium-silicate-hydrate from molecular dynamics simulations

Concrete is the most widely manufactured material in the world. Its binding phase, calcium-silicate-hydrate (C-S-H), is responsible for its mechanical properties and has an atomic structure fairly similar to that of usual calcium silicate glasses, which makes it appealing to study this material with tools and theories traditionally used for non-crystalline solids. Here, following this idea, we use molecular dynamics simulations to evaluate the fracture toughness of C-S-H, inaccessible experimentally. This allows us to discuss the brittleness of the material at the atomic scale. We show that, at this scale, C-S-H breaks in a ductile way, which prevents one from using methods based on linear elastic fracture mechanics. Knowledge of the fracture properties of C-S-H at the atomic scale opens the way for an upscaling approach to the design of tougher cement paste, which would allow for the design of slender environment-friendly infrastructures, requiring less material

Fracture Mechanisms in Organic-Rich Shales: Role of Kerogen

In this work we study role of kerogen in the fracture properties of organic-rich shales and, in particular, in the ductility of shales. The presence of kerogen and clays in shale is known to increase the ductility. We propose here a multiscale approach to develop a fine understanding of shale ductility from the molecular scale. We develop and validate a methodology at the molecular scale that can capture the toughness and ductility of a material. We apply this methodology successfully to a silica polymorph and to a kerogen analog, and we confirm the significant ductility of kerogen. Interestingly the silica-kerogen interface exhibits a similar ductility, which is central for the properties of the heterogeneous shale. Finally, we consider a tentative upscaling considering the pull out phenomenon as a likely mechanism of fracture of the shale