On generalized multiscale methods for flow in complex porous media and their applications

Abstract

In this document, the Generalized Multiscale Finite Element Method (GMsFEM) is studied, which deals with constructing multiscale spectral basis functions designed for high-contrast multiscale problems. The multiscale basis functions are built from the product of the eigenvectors, computed from a local spectral problem and a partition of unity over the study domain. The eigenvalues detect essential features of the solutions that are not captured by the initial multiscale basis functions. This document reviews the general convergence study where the error estimates are written in terms of the eigenvalues associated with the eigenvectors not used in the construction. Error analysis involves local and global norms that measure the convergence speed of the expansion of the solution in terms of local eigenvectors; this is achieved with a careful choice of the initial multiscale basis functions and the configuration of the eigenvalue problems. Two novel important numerical applications are presented: the first is the free-boundary dam problem posed on a heterogeneous high-contrast medium, where we introduce a fictitious time variable that motivates an adequate time discretization that can be understood as a fixed-point iteration. For the steady-state solution, we use the duality method to deal with the multivalued nonlinear terms involved; then, efficient approximations of pressure and saturation are calculated using the GMsFEM method. The second application is the solution of a parabolic equation. Here implementing time discretizations, such as finite differences or exponential integrators in the presence of a high contrast coefficient, it may not be practical in because each time iteration one needs the computation of matrix operators involving very large and extremely ill-conditioned sparse matrices. The GMsFEM is essential since it allows obtaining the solution of the problem more simply, allowing to combine the GMsFEM with the method of exponential integrators in time to get a good approximation of the final temporary solution. (Texto tomado de la fuente)

En este documento se estudia el M´etodo de Elementos Finitos Multiescala Generalizados (GMsFEM), el cual trata de la construcci´on de funciones base espectrales multiescala que est´an dise˜nadas para problemas de alto contraste. Las funciones base multiescala se construyen a partir del producto entre los vectores propios, construidos a partir de un problema espectral local y una partici´on de la unidad sobre el dominio de estudio. Los valores propios detectan caracter´ısticas importantes de las soluciones que no son capturadas por las funciones base multiescala iniciales. En este trabajo, se presenta un estudio de convergencia donde las estimaciones de error son generales, y est´an escritas en t´erminos de los valores propios asociados a los vectores propios no utilizados en la construcci´on. El an´alisis de errores implica normas locales y globales que miden la descomposici´on de la expansi´on de la soluci´on en t´erminos de vectores propios locales, esto se logra con una elecci´on cuidadosa de las funciones de base multiescala iniciales y la configuraci´on de los problemas de valores propios. Se presentan dos aplicaciones num´ericas importantes: la primera, es el problema de represa con frontera libre planteado sobre un medio heterog´eneo de alto contraste, donde introducimos una variable de tiempo ficticia que motiva una discretizaci´on de tiempo adecuada que puede entenderse como una iteraci´on de punto fijo a la soluci´on de estado estacionario, y usamos el m´etodo de dualidad para tratar con los t´erminos no lineales multivaluados involucrados; luego, se calculan aproximaciones eficientes de la presi´on y la saturaci´on usando el m´etodo GMsFEM. La segunda aplicaci´on es la soluci´on de una ecuaci´on parab´olica donde al implementar discretizaciones de tiempo como diferencias finitas o integradores exponenciales sobre un coeficiente de alto contraste, puede no ser pr´actico porque cada iteraci´on de tiempo necesita el c´alculo de operadores matriciales que involucran matrices dispersas, muy grandes y mal condicionadas; es por esto que el GMsFEM es importante ya que permite la obtenci´on de la soluci´on del problema de una forma m´as sencilla, permitiendo combinar GMsFEM con el m´etodo de integradores exponenciales en el tiempo para obtener una buena aproximaci´on de la soluci´on temporal final