A Dual-Porosity-Stokes Model and Finite Element Method for Coupling Dual-Porosity Flow and Free Flow
- Award ID(s):
- 1418624
- PAR ID:
- 10249881
- Date Published:
- Journal Name:
- SIAM Journal on Scientific Computing
- Volume:
- 38
- Issue:
- 5
- ISSN:
- 1064-8275
- Page Range / eLocation ID:
- B710 to B739
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Summary 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.more » « less
-
Abstract The emergence of the “mush paradigm” has raised several questions for conventional models of magma storage and extraction: how are melts extracted to form eruptible liquid-rich domains? What mechanism controls melt transport in mush-rich systems? Recently, reactive flow has been proposed as a major contributing factor in the formation of high porosity, melt-rich regions. Yet, owing to the absence of accurate geochemical simulations, the influence of reactive flow on the porosity of natural mush systems remains under-constrained. Here, we use a thermodynamically constrained model of melt-mush reaction to simulate the chemical, mineralogical, and physical consequences of reactive flow in a multi-component mush system. Our results demonstrate that reactive flow within troctolitic to gabbroic mushes can drive large changes in mush porosity. For example, primitive magma recharge causes an increase in the system porosity and could trigger melt channelization or mush destabilization, aiding rapid melt transfer through low-porosity mush reservoirs.more » « less