Maria Vasilyeva

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.

Multiscale model reduction for shale gas transport in fractured media

In this paper, we develop a multiscale model reduction technique that describes shale gas transport in fractured media. Due to the pore-scale heterogeneities and processes, we use upscaled models to describe the matrix. We follow our previous work (Akkutlu et al. Transp. Porous Media 107(1), 235–260, 2015), where we derived an upscaled model in the form of generalized nonlinear diffusion model to describe the effects of kerogen. To model the interaction between the matrix and the fractures, we use Generalized Multiscale Finite Element Method (Efendiev et al. J. Comput. Phys. 251, 116–135, 2013, 2015). In this approach, the matrix and the fracture interaction is modeled via local multiscale basis functions. In Efendiev et al. (2015), we developed the GMsFEM and applied for linear flows with horizontal or vertical fracture orientations aligned with a Cartesian fine grid. The approach in Efendiev et al. (2015) does not allow handling arbitrary fracture distributions. In this paper, we (1) consider arbitrary fracture distributions on an unstructured grid; (2) develop GMsFEM for nonlinear flows; and (3) develop online basis function strategies to adaptively improve the convergence. The number of multiscale basis functions in each coarse region represents the degrees of freedom needed to achieve a certain error threshold. Our approach is adaptive in a sense that the multiscale basis functions can be added in the regions of interest. Numerical results for two-dimensional problem are presented to demonstrate the efficiency of proposed approach.

Splitting schemes for poroelasticity and thermoelasticity problems

Abstract In this work, we consider the coupled systems of linear unsteady partial differential equations, which arise in the modelling of poroelasticity processes. Stability estimates of weighted difference schemes for the coupled system of equations are presented. Approximation in space is based on the finite element method. We construct splitting schemes and give some numerical comparisons for typical poroelasticity problems. The results of numerical simulation of a 3D problem are presented. Special attention is given to using high performance computing systems.

Mixed GMsFEM for second order elliptic problem in perforated domains

We consider a class of second order elliptic problems in perforated domains with homogeneous Neumann boundary condition. It is well-known that numerically solving these problems require a very fine computational mesh, and some model reduction techniques are therefore necessary. We will develop a new model reduction technique based on the generalized multiscale finite element method (GMsFEM). The GMsFEM has been applied successfully to second order elliptic problems in perforated domains with Dirichlet boundary conditions Chung et?al. (2015). However, due to the use of multiscale partition of unity functions, the same method cannot be applied to the case with Neumann boundary conditions. The aim of this paper is to develop a new mixed GMsFEM, based on a piecewise constant approximation for pressure and a multiscale approximation for velocity, giving a mass conservative method. The method can handle the Neumann boundary condition naturally. The multiscale basis functions for velocity are constructed by some carefully chosen local snapshot spaces and local spectral decompositions. The spectral convergence of the method is analyzed. Moreover, by using some local error indicators, the basis functions can be added locally and adaptively. We also consider an online procedure for the construction of new basis functions in the online stage in order to capture the distant effects. We will present some numerical examples to show the performance of the method.

Online adaptive local multiscale model reduction for heterogeneous problems in perforated domains

In this paper, we develop and analyze an adaptive multiscale approach for heterogeneous problems in perforated domains. We consider commonly used model problems including the Laplace equation, the elasticity equation, and the Stokes system in perforated regions. In many applications, these problems have a multiscale nature arising because of the perforations, their geometries, the sizes of the perforations, and configurations. Typical modeling approaches extract average properties in each coarse region, that encapsulate many perforations, and formulate a coarse-grid problem. In some applications, the coarse-grid problem can have a different form from the fine-scale problem, e.g. the coarse-grid system corresponding to a Stokes system in perforated domains leads to Darcy equations on a coarse grid. In this paper, we present a general offline/online procedure, which can adequately and adaptively represent the local degrees of freedom and derive appropriate coarse-grid equations. Our approaches start with the offline procedure, which constructs multiscale basis functions in each coarse region and formulates coarse-grid equations. We presented the offline simulations without the analysis and adaptive procedures, which are needed for accurate and efficient simulations. The main contributions of this paper are (1) the rigorous analysis of the offline approach, (2) the development of the online procedures and their analysis, and (3) the development of adaptive strategies. We present an online procedure, which allows adaptively incorporating global information and is important for a fast convergence when combined with the adaptivity. We present online adaptive enrichment algorithms for the three model problems mentioned above. Our methodology allows adding and guides constructing new online multiscale basis functions adaptively in appropriate regions. We present the convergence analysis of the online adaptive enrichment algorithm for the Stokes system. In particular, we show that the online procedure has a rapid convergence with a rate related to the number of offline basis functions, and one can obtain fast convergence by a sufficient number of offline basis functions, which are computed in the offline stage. The convergence theory can also be applied to the Laplace equation and the elasticity equation. To illustrate the performance of our method, we present numerical results with both small and large perforations. We see that only a few (1 or 2) online iterations can significantly improve the offline solution.

A Generalized Multiscale Finite Element Method for poroelasticity problems I

In this paper, we consider the numerical solution of poroelasticity problems that are of Biot type and develop a general algorithm for solving coupled systems. We discuss the challenges associated with mechanics and flow problems in heterogeneous media. The two primary issues being the multiscale nature of the media and the solutions of the fluid and mechanics variables traditionally developed with separate grids and methods. For the numerical solution we develop and implement a Generalized Multiscale Finite Element Method (GMsFEM) that solves problem on a coarse grid by constructing local multiscale basis functions. The procedure begins with construction of multiscale bases for both displacement and pressure in each coarse block. Using a snapshot space and local spectral problems, we construct a basis of reduced dimension. Finally, after multiplying by a multiscale partitions of unity, the multiscale basis is constructed in the offline phase and the coarse grid problem then can be solved for arbitrary forcing and boundary conditions. We implement this algorithm on two heterogeneous media and compute error between the multiscale solution with the fine-scale solutions. Randomized oversampling and forcing strategies are also tested.

A Generalized Multiscale Finite Element Method for poroelasticity problems II: nonlinear coupling

In this paper, we consider the numerical solution of some nonlinear poroelasticity problems that are of Biot type and develop a general algorithm for solving nonlinear coupled systems. We discuss the difficulties associated with flow and mechanics in heterogenous media with nonlinear coupling. The central issue being how to handle the nonlinearities and the multiscale scale nature of the media. To compute an efficient numerical solution we develop and implement a Generalized Multiscale Finite Element Method (GMsFEM) that solves nonlinear problems on a coarse grid by constructing local multiscale basis functions and treating part of the nonlinearity locally as a parametric value. After linearization with a Picard Iteration, the procedure begins with construction of multiscale bases for both displacement and pressure in each coarse block by treating the staggered nonlinearity as a parametric value. Using a snapshot space and local spectral problems, we construct an offline basis of reduced dimension. From here an online, parametric dependent, space is constructed. Finally, after multiplying by a multiscale partitions of unity, the multiscale basis is constructed and the coarse grid problem then can be solved for arbitrary forcing and boundary conditions. We implement this algorithm on a geometry with a linear and nonlinear pressure dependent permeability field and compute error between the multiscale solution with the fine-scale solutions.

On numerical homogenization of shale gas transport

In this paper, we study the numerical homogenization for the gas transport in organic rich shale with heterogeneous kerogen distribution. We consider organic rich shale as the domain with two subdomains: inorganic matrix and kerogen. The processes in both regions are described by nonlinear parabolic equations, which take into account the filtration, diffusion, and adsorption. We follow the work of Yucel Akkutlu et?al. (2015), where the authors develop homogenization techniques for the gas transport in organic rich shale. We use the framework of Yucel Akkutlu et?al. (2015) and develop numerical homogenization for two dimensional examples. We discuss local subgrid calculations and the approaches to compute the effective properties numerically. In our methods, the local problems use Dirichlet boundary conditions. We compare the properties of the fine-grid reference solution against those of the coarse grid model. Our approaches show a good agreement between macroscopic quantities.

Numerical Simulation of Thermoelasticity Problems on High Performance Computing Systems

In this work we consider the coupled linear system of equations for temperature and displacements which describes the thermoelastic behaviour of the body. For numerical solution we approximate our system using finite element method. As model problem for simulation we consider the thermomechanical state of the ceramic substrates with metallization, which are used for the manufacturing of light-emitting diode modules. The results of numerical simulation of the 3D problem in the complex geometric area are presented.

Multiscale model reduction for transport and flow problems in perforated domains

Abstract Convection-dominated transport phenomenon is important for many applications. In these applications, the transport velocity is often a solution of heterogeneous flow problems, which results to a coupled flow and transport phenomena. In this paper, we consider a coupled flow (Stokes problem) and transport (unsteady convection–diffusion problem) in perforated domains. Perforated domains (see Fig. 1) represent void space outside hard inclusions as in porous media, filters, and so on. We construct a coarse-scale solver based on Generalized Multiscale Finite Element Method (GMsFEM) for a coupled flow and transport. The main idea of the GMsFEM is to develop a systematic approach for computing multiscale basis functions. We use a mixed formulation and appropriate multiscale basis functions for both flow and transport to guarantee a mass conservation. For the transport problem, we use Petrov–Galerkin mixed formulation, which provides a stability. As a first approach, we use the multiscale flow solution in constructing the basis functions for the transport equation. In the second approach, we construct multiscale basis functions for coupled flow and transport without solving global flow problem. The novelty of this approach is to construct a coupled multiscale basis function. Numerical results are presented for computations using offline basis. We also present an algorithm for adaptively adding online multiscale basis functions, which are computed using the residual information. Numerical examples using online GMsFEM show the speed up of convergence.