Abstract
This article presents an analysis of the implementation of the adaptive mesh refinement (AMR) technique to one-dimensional three-phase flows in heterogeneous fractured media. Variations of matrix permeability and matrix–fracture transfer shape factor cause matrix saturation discontinuities. In order to avoid directly performing refining or coarsening operations on these matrix quantities, we suggest using two separated grid systems for the fracture and matrix equations, respectively. The fracture and matrix equations are decoupled by a special decomposition approach during the Newton-Raphson procedure for solving the resulting nonlinear equations. The numerical examples show that the proposed AMR technique is fast with good accuracy.
This work is supported by the National Hi-Tech Research and Development Program of China (863 Program No. 2006AA06Z228), the National Natural Science Foundation of China (No. 50674081), and the Program for New Century Excellent Talents in the University of China.