Coupled and Decoupled Stabilized Mixed Finite Element Methods for Nonstationary Dual-Porosity-Stokes Fluid Flow Model


In this paper, we propose and analyze two stabilized mixed finite element methods for the dual-porosity-Stokes model, which couples the free flow region and microfracture-matrix system through four interface conditions on an interface. The first stabilized mixed finite element method is a coupled method in the traditional format. Based on the idea of partitioned time stepping, the four interface conditions, and the mass exchange terms in the dual-porosity model, the second stabilized mixed finite element method is decoupled in two levels and allows a noniterative splitting of the coupled problem into three subproblems. Due to their superior conservation properties and convenience of the computation of flux, mixed finite element methods have been widely developed for different types of subsurface flow problems in porous media. For the mixed finite element methods developed in this article, no Lagrange multiplier is used, but an interface stabilization term with a penalty parameter is added in the temporal discretization. This stabilization term ensures the numerical stability of both the coupled and decoupled schemes. The stability and the convergence analysis are carried out for both the coupled and decoupled schemes. Three numerical experiments are provided to demonstrate the accuracy, efficiency, and applicability of the proposed methods.


Mathematics and Statistics

Research Center/Lab(s)

Center for High Performance Computing Research

Keywords and Phrases

Decoupled numerical methods; Dual-porosity-Stokes model; Horizontal wellbore; Mixed finite elements; Stabilization

International Standard Serial Number (ISSN)

0029-5981; 1097-0207

Document Type

Article - Journal

Document Version


File Type





© 2019 John Wiley & Sons Ltd, All rights reserved.

Publication Date

09 Nov 2019