© 2014 John Wiley & Sons, Ltd. This paper primarily deals with the computational aspects of chemical dissolution-front instability problems in two-dimensional fluid-saturated porous media under non-isothermal conditions. After the dimensionless governing partial differential equations of the non-isothermal chemical dissolution-front instability problem are briefly described, the formulation of a computational procedure, which contains a combination of using the finite difference and finite element method, is derived for simulating the morphological evolution of chemical dissolution fronts in the non-isothermal chemical dissolution system within two-dimensional fluid-saturated porous media. To ensure the correctness and accuracy of the numerical solutions, the proposed computational procedure is verified through comparing the numerical solutions with the analytical solutions for a benchmark problem. As an application example, the verified computational procedure is then used to simulate the morphological evolution of chemical dissolution fronts in the supercritical non-isothermal chemical dissolution system. The related numerical results have demonstrated the following: (1) the proposed computational procedure can produce accurate numerical solutions for the planar chemical dissolution-front propagation problem in the non-isothermal chemical dissolution system consisting of a fluid-saturated porous medium; (2) the Zhao number has a significant effect not only on the dimensionless propagation speed of the chemical dissolution front but also on the distribution patterns of the dimensionless temperature, dimensionless pore-fluid pressure, and dimensionless chemical-species concentration in a non-isothermal chemical dissolution system; (3) once the finger penetrates the whole computational domain, the dimensionless pore-fluid pressure decreases drastically in the non-isothermal chemical dissolution system.
|Journal||International Journal for Numerical Methods in Engineering|
|Publication status||Published - 2015|