An official website of the United States government
The objective of this paper is to develop efficient numerical algorithms for the linear advection-diffusion equation in fractured porous media. A reduced fracture model is considered where the fractures are treated as interfaces between subdomains and the interactions between the fractures and the surrounding porous medium are taken into account. The model is discretized by a backward Euler upwind-mixed hybrid finite element method in which the flux variable represents both the advective and diffusive fluxes. The existence, uniqueness, as well as optimal error estimates in both space and time for the fully discrete coupled problem are established. Moreover, to facilitate different time steps in the fracture-interface and the subdomains, global-in-time, nonoverlapping domain decomposition is utilized to derive two implicit iterative solvers for the discrete problem. The first method is based on the time-dependent Steklov–Poincaré operator, while the second one employs the optimized Schwarz waveform relaxation (OSWR) approach with Ventcel-Robin transmission conditions. A discrete space-time interface system is formulated for each method and is solved iteratively with possibly variable time step sizes. The convergence of the OSWR-based method with conforming time grids is also proved. Finally, numerical results in two dimensions are presented to verify the optimal order of convergence of the monolithic solver and to illustrate the performance of the two decoupled schemes with local time-stepping on problems of high Péclet numbers.
more »
« less
Fast and Accuracy-Preserving Domain Decomposition Methods for Reduced Fracture Models with Nonconforming Time Grids
This paper is concerned with the numerical solution of compressible fluid flow in a fractured porous medium. The fracture represents a fast pathway (i.e., with high permeability) and is modeled as a hypersurface embedded in the porous medium. We aim to develop fast-convergent and accurate global-in-time domain decomposition (DD) methods for such a reduced fracture model, in which smaller time step sizes in the fracture can be coupled with larger time step sizes in the subdomains. Using the pressure continuity equation and the tangential PDEs in the fracture-interface as transmission conditions, three different DD formulations are derived; each method leads to a space-time interface problem which is solved iteratively and globally in time. Efficient preconditioners are designed to accelerate the convergence of the iterative methods while preserving the accuracy in time with nonconforming grids. Numerical results for two-dimensional problems with non-immersed and partially immersed fractures are presented to show the improved performance of the proposed methods.
more »
« less
- Award ID(s):
- 1912626
- PAR ID:
- 10473527
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Journal of Scientific Computing
- Volume:
- 96
- Issue:
- 1
- ISSN:
- 0885-7474
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This paper is concerned with the numerical solution of the flow problem in a fractured porous medium where the fracture is treated as a lower dimensional object embedded in the rock matrix. We consider a space-time mixed variational formulation of such a reduced fracture model with mixed finite element approximations in space and discontinuous Galerkin discretization in time. Different spatial and temporal grids are used in the subdomains and in the fracture to adapt to the heterogeneity of the problem. Analysis of the numerical scheme, including well-posedness of the discrete problem, stability and a priori error estimates, is presented. Using substructuring techniques, the coupled subdomain and fracture system is reduced to a space-time interface problem which is solved iteratively by GMRES. Each GMRES iteration involves solution of time-dependent problems in the subdomains using the method of lines with local spatial and temporal discretizations. The convergence of GMRES is proved by using the field-of-values analysis and the properties of the discrete space-time interface operator. Numerical experiments are carried out to illustrate the performance of the proposed iterative algorithm and the accuracy of the numerical solution.more » « less
-
In this work, we present a numerical method that provides accurate real-time detection for the widthsof the fractures in a fractured porous medium based on observational data on porous medium fluidmass and velocity. To achieve this task, an inverse problem is formulated by first constructing aforward formulation based on the reduced fracture model of the diffusion equation. A parameter estimationproblem is then performed online by utilizing a direct filter method. Numerical experimentsare carried out to demonstrate the accuracy of our method in approximating the target parameters.more » « less
-
In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn– Hilliard–Navier–Stokes equations in the free flow region and Cahn–Hilliard–Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions.more » « less
-
We consider the interaction between a free flowing fluid and a porous medium flow, where the free flowing fluid is described using the time dependent Stokes equations, and the porous medium flow is described using Darcy’s law in the primal formulation. To solve this problem numerically, we use a diffuse interface approach, where the weak form of the coupled problem is written on an extended domain which contains both Stokes and Darcy regions. This is achieved using a phase-field function which equals one in the Stokes region and zero in the Darcy region, and smoothly transitions between these two values on a diffuse region of width (ϵ) around the Stokes-Darcy interface. We prove convergence of the diffuse interface formulation to the standard, sharp interface formulation, and derive rates of convergence. This is performed by deriving a priori error estimates for discretizations of the diffuse interface method, and by analyzing the modeling error of the diffuse interface approach at the continuous level. The convergence rates are also shown computationally in a numerical example.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s10915-023-02251-0
