Search

This paper presents a numerical procedure for computational modeling of moving interfaces between fluid and porous medium domains. To avoid the direct description of the interface boundary conditions at the interface between the fluid and porous medium domains, Darcy's law is used to simulate fluid flows in both the fluid and porous medium domains. For the purpose of effectively simulating the fluid flow using Darcy's law, the artificial permeability of the fluid domain is used to establish a permeability ratio between the artificial permeability in the fluid domain and the real permeability in the porous medium domain. Using the proposed permeability ratio, the ratio of fluid pressure gradient in the porous medium domain to that in the fluid domain can be appropriately simulated. To verify the proposed numerical procedure, analytical solutions have been derived for a benchmark problem, which can be simulated using the proposed numerical procedure. Comparison of the numerical solutions with the derived analytical solutions has demonstrated the correctness and accuracy of the proposed numerical procedure. Through applying the proposed numerical procedure to several examples associated with the fluid-porous medium interface propagation problems, the related numerical solutions have demonstrated that: (1) the proposed numerical procedure is capable of simulating the morphological instability of the fluid-porous medium interface in a coupled fluid flow-chemical dissolution system when the system is in a supercritical state; and (2) the permeability ratio can have a considerable effect on the evolved morphologies of the fluid-porous medium interface in the coupled fluid flow-chemical dissolution system. [ABSTRACT FROM AUTHOR]

In dealing with chemical-dissolution-front propagation problems in fluid-saturated porous media, the chemical dissolution front represented by the porosity of the medium may have a very steep slope (i.e., a very large porosity gradient) at the dissolution front, depending on the mineral dissolution ratio that is defined as the equilibrium concentration of the dissolved minerals in the pore-fluid to the solid molar density of the dissolvable minerals in the solid matrix. When the mineral dissolution ratio approaches zero, the theoretical value of the porosity gradient tends to infinity at the chemical dissolution front. Even for a very small value of the mineral dissolution ratio, which is very common in geochemical systems, the porosity gradient can be large enough to cause the solution hard to converge when the conventional finite element method is used to solve a chemical dissolution problem in a fluid-saturated porous medium where the pore-fluid is compressible. To improve the convergent speed of solution, a porosity-gradient replacement approach, in which the term involving porosity-gradient computation is replaced by a new term consisting of pore-fluid density-gradient and pressure-gradient computation, is first proposed and then incorporated into the finite element method in this study. Through comparing the numerical results obtained from the proposed approach with the theoretical solutions for a benchmark problem, it has been demonstrated that not only can the solution divergence be avoid, but also the accurate simulation results can be obtained when the proposed porosity-gradient replacement approach is used to solve chemical-dissolution-front propagation problems in fluid-saturated porous media including pore-fluid compressibility. [ABSTRACT FROM AUTHOR]