skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, May 16 until 2:00 AM ET on Saturday, May 17 due to maintenance. We apologize for the inconvenience.


Title: A Dual-Porosity-Stokes Model and Finite Element Method for Coupling Dual-Porosity Flow and Free Flow
Award ID(s):
1418624
PAR ID:
10249881
Author(s) / Creator(s):
; ; ; ; ;
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
  1. 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
  2. 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