A Phase-Field Discrete Element Method to study chemo-mechanical coupling in granular materials

Alexandre Sac-Morane,Manolis Veveakis,Hadrien Rattez

DOI: https://doi.org/10.1016/j.cma.2024.116900

2023-10-24

Abstract:This paper presents an extension of the discrete element method using a phase-field formulation to incorporate grain shape and its evolution. The introduction of a phase variable enables an effective representation of grain geometry and facilitates the application of physical laws, such as chemo-mechanical couplings, for modeling shape changes. These physical laws are solved numerically using the finite element method coupled in a staggered scheme to the discrete element model. The efficacy of the proposed Phase-Field Discrete Element Model (PFDEM) is demonstrated through its ability to accurately capture the real grain shape in a material subjected to dissolution only and compute the stress evolution. It is then applied to model the phenomenon of pressure solution involving dissolution and precipitation in granular materials at the microscale and enables to reproduce the creep response observed experimentally. This framework contributes to the enhanced understanding and simulation of complex behaviors in granular materials and sedimentary rocks for many geological processes like diagenesis or earthquake nucleation.

Materials Science,Geophysics

What problem does this paper attempt to address?

The paper aims to address the complex interactions between chemical and mechanical processes in soil and rock, particularly in the chemo-mechanical coupling phenomena of granular materials. Specifically, the objectives of the paper include:

1. **Introducing the Phase Field Discrete Element Method (PFDEM)**: An extended Discrete Element Method (DEM) is developed, incorporating particle shapes and their evolution through phase field formulations. This method effectively represents particle geometries and facilitates the application of physical laws to simulate shape changes.

2. **Investigating the Impact of Particle Shape**: Conducting oedometer tests under acid injection conditions to highlight the influence of particle shape on experimental results.

3. **Simulating Pressure Solution Phenomena**: Simulating pressure solution phenomena through micro-scale dissolution and precipitation processes, and reproducing the creep response observed in experiments. This helps to better understand and simulate the complex behavior of granular materials and sedimentary rocks, as well as many geological processes such as diagenesis or earthquake nucleation.

Through these efforts, the paper aims to enhance the understanding and simulation capabilities of the complex behavior in granular materials and sedimentary rocks, especially in engineering applications such as underground storage and geothermal energy utilization, as well as geological processes like diagenesis and earthquake nucleation.