Guanglian Li

Generalized multiscale finite element methods for problems in perforated heterogeneous domains

Complex processes in perforated domains occur in many real-world applications. These problems are typically characterized by physical processes in domains with multiple scales. Moreover, these problems are intrinsically multiscale and their discretizations can yield very large linear or nonlinear systems. In this paper, we investigate multiscale approaches that attempt to solve such problems on a coarse grid by constructing multiscale basis functions in each coarse grid, where the coarse grid can contain many perforations. In particular, we are interested in cases when there is no scale separation and the perforations can have different sizes. In this regard, we mention some earlier pioneering works, where the authors develop multiscale finite element methods. In our paper, we follow Generalized Multiscale Finite Element Method (GMsFEM) and develop a multiscale procedure where we identify multiscale basis functions in each coarse block using snapshot space and local spectral problems. We show that with a few basis functions in each coarse block, one can approximate the solution, where each coarse block can contain many small inclusions. We apply our general concept to (1) Laplace equation in perforated domains; (2) elasticity equation in perforated domains; and (3) Stokes equations in perforated domains. Numerical results are presented for these problems using two types of heterogeneous perforated domains. The analysis of the proposed methods will be presented elsewhere.

Generalized Multiscale Finite Element Methods. Nonlinear Elliptic Equations

In this paper we use the Generalized Multiscale Finite Element Method (GMsFEM) framework, introduced in [20], in order to solve nonlinear elliptic equations with high-contrast coefficients. The proposed solution method involves linearizing the equation so that coarse-grid quantities of previous solution iterates can be regarded as auxiliary parameters within the problem formulation. With this convention, we systematically construct respective coarse solution spaces that lend themselves to either continuous Galerkin (CG) or discontinuous Galerkin (DG) global formulations. Here, we use Symmetric Interior Penalty Discontinuous Galerkin approach. Both methods yield a predictable error decline that depends on the respective coarse space dimension, and we illustrate the effectiveness of the CG and DG formulations by offering a variety of numerical examples.

A generalized multiscale finite element method for the Brinkman equation

In this paper we consider the numerical upscaling of the Brinkman equation in the presence of high-contrast permeability fields. We develop and analyze a robust and efficient Generalized Multiscale Finite Element Method (GMsFEM) for the Brinkman model in two dimensions. In the fine grid, we use mixed finite element method with the velocity and pressure being continuous piecewise quadratic and piecewise constant finite element spaces, respectively. Using the GMsFEM framework we construct suitable coarse-scale spaces for the velocity and pressure that yield a robust mixed GMsFEM. We develop a novel approach to construct a coarse approximation for the velocity snapshot space and a robust small offline space for the velocity space. The stability of the mixed GMsFEM and a priori error estimates are derived. A variety of two-dimensional numerical examples are presented to illustrate the effectiveness of the algorithm.

On the Decay Rate of the Singular Values of Bivariate Functions

In this work, we establish a new truncation error estimate of the singular value decomposition (SVD) for a class of Sobolev smooth bivariate functions

Homogenization of High-Contrast Brinkman Flows

\kappa {\,\in\,} L^2(\Omega,H^s(D))

Multiscale Modeling of High Contrast Brinkman Equations with Applications to Deformable Porous Media.

,

Low-Rank Approximation to Heterogeneous Elliptic Problems

s{\,\...

Upscaled HDG Methods for Brinkman Equations with High-Contrast Heterogeneous Coefficient

Modeling porous flow in complex media is a challenging problem. Not only is the problem inherently multiscale but, due to high contrast in permeability values, flow velocities may differ greatly throughout the medium. To avoid complicated interface conditions, the Brinkman model is often used for such flows [O. Iliev, R. Lazarov, and J. Willems, Multiscale Model. Simul., 9 (2011), pp. 1350--1372]. Instead of permeability variations and contrast being contained in the geometric media structure, this information is contained in a highly varying and high-contrast coefficient. In this work, we present two main contributions. First, we develop a novel homogenization procedure for the high-contrast Brinkman equations by constructing correctors and carefully estimating the residuals. Understanding the relationship between scales and contrast values is critical to obtaining useful estimates. Therefore, standard convergence-based homogenization techniques [G. A. Chechkin, A. L. Piatniski, and A. S. Shamev, Homogen...

On the degree of ill-posedness of multi-dimensional magnetic particle imaging

Simulating porous media flows has a wide range of applications. Often, these applications involve many scales and multi physical processes. A useful tool in the analysis of such problems in that of homogenization as an averaged description is derived circumventing the need for complicated simulation of the fine scale features. In this work, we recall recent developments of homogenization techniques in the application of flows in deformable porous media. In addition, homogenization of media with high-contrast. In particular, we recall the main ideas of the homogenization of slowly varying Stokes flow and summarize the results of [4]. We also present the ideas for extending these techniques to high-contrast deformable media [3]. These ideas are connected by the modeling of multiscale fluid-structure interaction problems.

A Convergent Adaptive Finite Element Method for Cathodic Protection

In this work, we investigate the low-rank approximation of elliptic problems in heterogeneous media by means of Kolmogrov