This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
MIXED FINITE ELEMENT-BASED FULLY CONSERVATIVE METHODS FOR SIMULATING WORMHOLE PROPAGATION *JISHENG KOU † , SHUYU SUN ‡ , AND YUANQING WU § Abstract. Wormhole propagation during reactive dissolution of carbonates plays a very important role in the product enhancement of oil and gas reservoir. Because of high velocity and nonuniform porosity, the Darcy-Forchheimer model is applicable for this problem instead of conventional Darcy framework. We develop a mixed finite element scheme for numerical simulation of this problem, in which mixed finite element methods are used not only for the Darcy-Forchheimer flow equations but also for the solute transport equation by introducing an auxiliary flux variable to guarantee full mass conservation. In theoretical analysis aspects, based on the cut-off operator of solute concentration, we construct an analytical function to control and handle the change of porosity with time; we treat the auxiliary flux variable as a function of velocity and establish its properties; we employ the coupled analysis approach to deal with the fully coupling relation of multivariables. From this, the stability analysis and a priori error estimates for velocity, pressure, concentration and porosity are established in different norms. Numerical results are also given to verify theoretical analysis and effectiveness of the proposed scheme.