BDDC for mixed-hybrid formulation of flow in porous media with combined mesh dimensions
Loading...
Links to Files
Author/Creator
Author/Creator ORCID
Date
2015-05-20
Type of Work
Department
Program
Citation of Original Publication
Šístek, Jakub; Březina, Jan; Sousedík, Bedřich; BDDC for mixed-hybrid formulation of flow in porous media with combined mesh dimensions; Numerical Linear Algebra with Applications, 22, 6, pages 903-929, 20 May, 2015; https://doi.org/10.1002/nla.1991
Rights
This item is likely protected under Title 17 of the U.S. Copyright Law. Unless on a Creative Commons license, for uses protected by Copyright Law, contact the copyright holder or the author.
This is the pre-peer reviewed version of the following article: Šístek, Jakub; Březina, Jan; Sousedík, Bedřich; BDDC for mixed-hybrid formulation of flow in porous media with combined mesh dimensions; Numerical Linear Algebra with Applications, 22, 6, pages 903-929, 20 May, 2015; https://doi.org/10.1002/nla.1991, which has been published in final form at https://onlinelibrary.wiley.com/doi/10.1002/nla.1991. This article may be used for noncommercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions.
This is the pre-peer reviewed version of the following article: Šístek, Jakub; Březina, Jan; Sousedík, Bedřich; BDDC for mixed-hybrid formulation of flow in porous media with combined mesh dimensions; Numerical Linear Algebra with Applications, 22, 6, pages 903-929, 20 May, 2015; https://doi.org/10.1002/nla.1991, which has been published in final form at https://onlinelibrary.wiley.com/doi/10.1002/nla.1991. This article may be used for noncommercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions.
Subjects
Abstract
We extend the balancing domain decomposition by constraints (BDDC) method to flows in porous media discretised by mixed-hybrid finite elements with combined mesh dimensions. Such discretisations appear when major geological fractures are modelled by one-dimensional or two-dimensional elements inside three-dimensional domains. In this set-up, the global problem and the substructure problems have a symmetric saddle-point structure, containing a ‘penalty’ block due to the combination of meshes. We show that the problem can be reduced by means of iterative substructuring to an interface problem, which is symmetric and positive definite. The interface problem can thus be solved by conjugate gradients with the BDDC method as a preconditioner. A parallel implementation of this algorithm is incorporated into an existing software package for subsurface flow simulations. We study the performance of the iterative solver on several academic and real-world problems. Numerical experiments illustrate its efficiency and scalability. Copyright © 2015 John Wiley & Sons, Ltd.