Generalized discontinuous multiscale methods for flows in highly heterogeneous porous media.

摘要

This dissertation is devoted to the development, study and testing of numerical methods for elliptic and parabolic equations with heterogeneous coefficients. The motivation for this study is to meet the need for fast and robust methods for numerical upscaling and simulation of single and multi-phase fluid flow in highly heterogeneous porous media. We consider the multiscale model reduction technique in the framework of the discontinuous Galerkin (DG) and the hybridizable discontinuous Galerkin (HDG) finite element methods. First, we design multiscale finite element methods for second order elliptic equations by applying the symmetric interior penalty discontinuous Galekin finite element method. We propose two different types of finite element spaces on the coarse mesh within DG framework. The first type of spaces is based on a local spectral problem that uses an interior weighted L 2-norm and a boundary weighted L2-norm for computing the mass matrix. The second choice is based on generation of a snapshot space and subsequent selection of a subspace of a reduced dimension. Second, we develop multiscale model reduction methods within the HDG framework. We provide construction of several multiscale finite element spaces (related to the coarse-mesh edges) that guarantee a reasonable approximation on a reduced dimensional space of the numerical traces. In these approaches, we use local snapshot spaces and local spectral decomposition following the concept of Generalized Multiscale Finite Element Methods. We also provide a general framework for systematic construction of multiscale spaces. By using local snapshots we were able to add local features to the solution space and to avoid high dimensional representation of trace spaces. Further, we extend multiscale finite element methods within HDG method to nonlinear and/or time-dependent problems. These extensions demonstrate the potential of the proposed constructions for some advanced and more practical applications. For most of the proposed methods, we investigate their stability and derive error estimates for the approximate solutions. Furthermore we study the performance of all proposed methods on a representative number of numerical examples. In the numerical tests, we use various permeability data of highly heterogeneous porous media and contrasts ranging from 103 to 106. Since the exact solution is in general unknown, we first generate solutions on a very fine mesh and use them as reference solutions in our tests. The numerical results confirm the theoretical study of the accuracy of the proposed methods and their robustness with respect to the media contrast. Our numerical experiments also show that the proposed methods could be implemented in a practical and efficient way.