A new thermo-mechanical coupled DEM model with non-spherical grains for thermally induced damage of rocks

Thermally induced damage often occurs in rocks in geophysical systems. Discrete element method (DEM) is a useful tool to model this thermo-mechanical coupled process owing to its explicit representation of fracture initiation and propagation. However, the previous DEM models for this are mostly based on spherical discrete elements, which are not able to capture all consequences (e.g. high ratio of compressive to tensile strength) of real rocks (e.g. granite) composed of complex-geometry grains. In order to overcome this intrinsic limitation, we present a new model allowing to mimick thermally induced damage of brittle rock with non-spherical grains. After validations, the new model is used to study thermal gradient cracking with a special emphasis on the effects from rock heterogeneity. The obtained fracture initiation and propagation are consistent with experimental observations, which demonstrates the ability of current model to reproduce the thermally induced damage of rocks. Meanwhile, the results show that rock heterogeneity influences thermal gradient cracking significantly, and more micro cracks uniformly scattering around the borehole are induced in the heterogeneous sample, which is not good for applications such as nuclear waste disposal. The present model provides a promising approach at micro-scale to explore mechanisms of thermally induced damage of rocks in geological engineering.