Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 355, Issue -, Pages 1-21Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2017.11.002
Keywords
Porous media; Darcy flow; Domain decomposition; Multiscale approximation; Robin boundary conditions; Optimized Schwarz; Two-Lagrange-multiplier method
Funding
- Sao Paulo Research Foundation FAPESP [2013/07375-0]
- Brazilian National Research and Technology Council (CNPq)
- Brazilian oil company Petrobras
- FAPESP [2016/08115-0, 2015/02649-0]
- NSF-DMS [1514808]
- Science Without Borders/CNPq-Brazil grant
- UT Dallas
- Direct For Mathematical & Physical Scien
- Division Of Mathematical Sciences [1514808] Funding Source: National Science Foundation
Ask authors/readers for more resources
A multiscale mixed method aiming at the accurate approximation of velocity and pressure fields in heterogeneous porous media is proposed. The procedure is based on a new domain decomposition method in which the local problems are subject to Robin boundary conditions. The domain decomposition procedure is defined in terms of two independent spaces on the skeleton of the decomposition, corresponding to interface pressures and fluxes, that can be chosen with great flexibility to accommodate local features of the underlying permeability fields. The well-posedness of the new domain decomposition procedure is established and its connection with the method of Douglas et al. (1993) [12], is identified, also allowing us to reinterpret the known procedure as an optimized Schwarz (or Two-Lagrange-Multiplier) method. The multiscale property of the new domain decomposition method is indicated, and its relation with the Multiscale Mortar Mixed Finite Element Method (MMMFEM) and the Multiscale Hybrid-Mixed (MHM) Finite Element Method is discussed. Numerical simulations are presented aiming at illustrating several features of the new method. Initially we illustrate the possibility of switching from MMMFEM to MHM by suitably varying the Robin condition parameter in the new multiscale method. Then we turn our attention to realistic flows in high-contrast, channelized porous formations. We show that for a range of values of the Robin condition parameter our method provides better approximations for pressure and velocity than those computed with either the MMMFEM and the MHM. This is an indication that our method has the potential to produce more accurate velocity fields in the presence of rough, realistic permeability fields of petroleum reservoirs. (c) 2017 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available