Modelling of a mineral dissolution front propagation is of interest in a wide range of scientific and engineering fields. The dissolution of minerals often involves complex physico-chemical processes at the solid–liquid interface (at nano-scale), which at the micro-to-meso-scale can be simplified to the problem of continuously moving boundaries. In this work, we studied the diffusion-controlled congruent dissolution of minerals from a meso-scale phase transition perspective. The dynamic evolution of the solid–liquid interface, during the dissolution process, is numerically simulated by employing the Finite Element Method (FEM) and using the phase–field (PF) approach, the latter implemented in the open-source Multiphysics Object Oriented Simulation Environment (MOOSE). The parameterization of the PF numerical approach is discussed in detail and validated against the experimental results for a congruent dissolution case of NaCl (taken from literature) as well as on analytical models for simple geometries. In addition, the effect of the shape of a dissolving mineral particle was analysed, thus demonstrating that the PF approach is suitable for simulating the mesoscopic morphological evolution of arbitrary geometries. Finally, the comparison of the PF method with experimental results demonstrated the importance of the dissolution rate mechanisms, which can be controlled by the interface reaction rate or by the diffusive transport mechanism.
Anisotropic solid–liquid interface kinetics in silicon: an atomistically informed phase-field model
