Brittle anisotropic fracture propagation in quartz sandstone: insights from phase-field simulations

We developed a generalized multiphase-field modeling framework for addressing the problem of brittle fracture propagation in quartz sandstones at microscopic length scale. Within this numerical approach, the grain boundaries and crack surfaces are modeled as diffuse interfaces. The two novel aspects of the model are the formulations of (I) anisotropic crack resistance in order to account for preferential cleavage planes within each randomly oriented quartz grain and (II) reduced interfacial crack resistance for incorporating lower fracture toughness along the grain boundaries that might result in intergranular crack propagation. The presented model is capable of simulating the competition between inter- and transgranular modes of fracturing based on the nature of grain boundaries, while exhibiting preferred fracturing directions within each grain. In the full parameter space, the model can serve as a powerful tool to investigate the complicated fracturing processes in heterogeneous polycrystalline rocks comprising of grains of distinct elastic properties, cleavage planes, and grain boundary attributes. We demonstrate the performance of the model through the representative numerical examples.