In crystalline and porous media, damage and healing refer to variations of mechanical and physical properties induced by pore or crack evolution. The gap between microscopic and macroscopic models makes it infeasible to uniquely characterize the pore- and crack- scale mechanisms that control deformation and flow regimes, predict percolation thresholds coupled to changes of rock stiffness, or relate crack rebonding time to stiffness and permeability healing time. The goal of our research is, therefore, to understand and predict chemo-mechanical damage and healing processes by coupling microstructural and poromechanical models. Fundamental scientific questions include: Why do pores and cracks heal? How long do mechanical and hydraulic recovery take? How much energy does healing require?
To establish closed forms that relate microstructure changes to macroscopic damage and healing evolution, we use both bottom-up and top-to-bottom approaches. In the former, moments of probability density functions of microstructure descriptors define fabric tensors that are used as dissipation variables at the macroscale. In the latter, homogenization schemes are built for polycrystals and for porous media with inclusions endowed with eigenstress. Fabric-enriched and matrix-inclusion models are inspired from and validated against laboratory experiments and image analyses. After calibration and validation at the material point, models are implemented in codes based on the Finite Element Method (FEM) and the eXtended Finite Element Method (XFEM). When the continuum assumption cannot hold, at a given damage threshold, a cohesive segment is inserted in replacement of a portion of the damage process zone. Rigorous calibrations procedures for the cohesive strength and the cohesive energy release rate ensure the balance of energy dissipated at the micro and macro scales. In order to simulate multi-scale fracture propagation driven by fluid injection, interpolation on fracture elements is enriched with jump functions for displacements, and with level-set-based distance functions for fluid pressure, which ensures that displacements are discontinuous across the fracture, but that the pressure field remains continuous.
A micromechanical model allowed predicting damage induced by a set of crack families subject to mixed mode propagation. A phenomenological approach was also proposed to relate micro-crack density, deformation and stiffness variations by deriving damage and healing dissipation potentials. Thermodynamic consistency is ensured, even when multiple damage and healing mechanisms are coupled. Upscaling was done by enriching the macroscopic formulation with fabric tensors defined as moments of probability of crack length, crack spacing, local porosity, pore. Additionally, self-consistent homogenization schemes were established to model viscous damage and pressure-solution driven healing in salt rock and implemented in the FEM to simulate long-term storage in caverns. We used a Mori-Tanaka homogenization scheme to capture the stress redistribution and subsequent damage induced by the expansion of biotite minerals as they convert to vermiculite over time. Biotite chemical weathering was accounted for by introducing inclusions' eigenstrains in the formulation. Simulations show that, at typical saprolite depth, micro-crack propagation can be controlled by weathering. Under constrained displacements, mechanical damage triggers earlier if the biotite volume fraction is high or if biotite minerals present a large orientation angle compared to the horizontal. A computational tool was developed to simulate the propagation of a discrete fracture within a continuum damage process zone. An XFEM was proposed to simulate multi-scale fracture propagation driven by fluid injection in transverse isotropic porous media. Intrinsic anisotropy was accounted for at the continuum scale, by using a non-local damage model in which four equivalent strains are defined to distinguish tension and compression, parallel and perpendicular to the bedding. Simulation results highlight the influence of intrinsic anisotropy on fluid driven fracture propagation.
Energy demand, environment protection and public health preservation raise pressing needs for safe and sustainable geological storage systems. The mining operations required to recover mineral resources, store energy underground and dispose of waste in deep geological layers involve coupled mechanical, physical and chemical rock microstructure changes. Our research findings can be used to recommend the conditions of moisture and temperature necessary to minimize damage and/or enhance healing in rocks and to design safe and sustainable geological storage systems. In addition to preventing subsidence, borehole instabilities and contaminant leakage, the proposed models will be applicable for predicting fault rupture during earthquakes, improving manufacturing methods and designing sustainable materials. Original contributions include: A theory to predict pore geometry evolution upon multi-physics damage and healing processes; Creative mathematical models to describe pore network topology with geometric variables that control damage and healing; Fundamental relationships between pore-scale healing time and macroscopic mechanical recovery time - a step forward to bridge poromechanics and damage mechanics; Innovative computational methods to predict mechanical instabilities and percolation thresholds upon damage and healing; Realistic multi-physics simulations of geological storage. The rigorous integration of topology, thermodynamics, poromechanics and continuum mechanics will transform the theory of damage and healing mechanics and provide a framework to interpret rock stress path history from topology descriptors.
- X. Shen, C. Arson, 2018. An isotropic self-consistent homogenization scheme for chemo-mechanical healing driven by pressure solution in halite, International Journal of Solids and Structures, DOI: 10.1016/j.ijsolstr.2018.11.010
- W. Jin, C. Arson, 2018. Anisotropic nonlocal damage model for materials with intrinsic transverse isotropy, International Journal of Solids and Structures, DOI: 10.1016/j.ijsolstr.2018.01.020
- W. Jin, C. Arson, 2017. Micromechanics based discrete damage model with multiple non-smooth yield surfaces: theoretical formulation, numerical implementation and engineering applications, International Journal of Damage Mechanics, DOI:10.1177/1056789517695872
- W. Jin, C. Arson, 2017. Discrete Wing Crack based Damage Model for Brittle Solids, International Journal of Solids and Structures, DOI:10.1016/j.ijsolstr.2016.12.025.
- A. Pouya, C. Zhu, C. Arson, 2017. Self-consistent micromechanical approach for damage accommodation in rock-like polycrystalline materials, International Journal of Damage Mechanics, DOI: 10.1177/1056789517747665
- W. Jin, H. Xu, C. Arson, S. Busetti, 2016. A Multi-Scale Computation Tool Coupling Mode II Fracture Propagation and Damage Zone Evolution, International Journal of Numerical and Analytical Methods in Geomechanics, DOI: 10.1002/nag.2553
- A. Pouya, C. Zhu, C. Arson, 2015. Micro-Macro Approach of Salt Viscous Fatigue under Cyclic Loading, Mechanics of Materials, DOI: 10.1016/j.mechmat.2015.10.009
- C. Zhu, A. Pouya, C. Arson, 2015. Micro-macro analysis and phenomenological modeling of salt viscous damage and application to salt caverns, Rock Mechanics and Rock Engineering, DOI: 10.1007/s00603-015-0832-
- C. Zhu, C. Arson, 2014. A model of damage and healing coupling halite thermomechanical behavior to microstructure evolution, Geotechnical and Geological Engineering, DOI: 10.1007/s10706-014-9797-9.
- C. Arson, 2014. Generalized Stress Variables in Continuum Damage Mechanics, Mechanics Research Communications, DOI: 10.1016/j.mechrescom.2014.06.006.