Computational Simulation of Chemical Dissolution-Front Instability Problems Associated with Radially Divergent Flow in Fluid-Saturated Porous Media

This paper proposes a computational simulation procedure for simulating chemical dissolution-front instability problems, in which radially divergent flow is involved in fluid-saturated porous media. In the proposed computational simulation procedure, a combination of the finite element and finite difference methods is used to simulate a chemical dissolution-front instability problem involving radially divergent flow, while a new algorithm is used to apply a small perturbation to the problem. Particular attention is paid on simulating low-order modes of an unstable circular chemical dissolution-front propagating in a fluid-saturated porous medium, in which dissolvable materials only occupy a small part, so that the final porosity is remarkably smaller than unity when dissolvable materials are completely dissolved in the chemical dissolution system. To verify the proposed computational simulation procedure, analytical solutions for a benchmark chemical dissolution-front instability problem involving radially divergent flow are derived in a purely mathematical manner. The related computational simulation results have demonstrated that: (1) the proposed computational simulation procedure is correct and useful for simulating chemical dissolution-front instability problems, which are associated with both stable and unstable chemical dissolution systems involving radially divergent flow in fluid-saturated porous media; (2) the simulated shapes of the second-order, third-order and fifth-order modes associated with an unstable chemical dissolution-fronts are respectively an ellipse, a star of three angles and a star of five angles in the unstable chemical dissolution system involving radially divergent flow in the fluid-saturated porous medium; (3) although the heterogeneity of a porous medium can affect the propagation speed of a chemical dissolution-front, it does not affect the low-order mode shape in the unstable chemical dissolution system involving radially divergent flow in the fluid-saturated heterogeneous porous medium.