scholarly journals Online Coupled Generalized Multiscale Finite Element Method for the Poroelasticity Problem in Fractured and Heterogeneous Media

Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 298
Author(s):  
Aleksei Tyrylgin ◽  
Maria Vasilyeva ◽  
Dmitry Ammosov ◽  
Eric T. Chung ◽  
Yalchin Efendiev
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Heterogeneous Media ◽  
Coupled System ◽  
Coarse Grid ◽  
Basis Functions ◽  
Fine Grid ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Element Method

In this paper, we consider the poroelasticity problem in fractured and heterogeneous media. The mathematical model contains a coupled system of equations for fluid pressures and displacements in heterogeneous media. Due to scale disparity, many approaches have been developed for solving detailed fine-grid problems on a coarse grid. However, some approaches can lack good accuracy on a coarse grid and some corrections for coarse-grid solutions are needed. In this paper, we present a coarse-grid approximation based on the generalized multiscale finite element method (GMsFEM). We present the construction of the offline and online multiscale basis functions. The offline multiscale basis functions are precomputed for the given heterogeneity and fracture network geometry, where for the construction, we solve a local spectral problem and use the dominant eigenvectors (appropriately defined) to construct multiscale basis functions. To construct the online basis functions, we use current information about the local residual and solve coupled poroelasticity problems in local domains. The online basis functions are used to enrich the offline multiscale space and rapidly reduce the error using residual information. Only with appropriate offline coarse-grid spaces can one guarantee a fast convergence of online methods. We present numerical results for poroelasticity problems in fractured and heterogeneous media. We investigate the influence of the number of offline and online basis functions on the relative errors between the multiscale solution and the reference (fine-scale) solution.

scholarly journals Mixed Generalized Multiscale Finite Element Method for a Simplified Magnetohydrodynamics Problem in Perforated Domains

Computation ◽  
2020 ◽  
Vol 8 (2) ◽  
pp. 58
Author(s):  
Valentin Alekseev ◽  
Qili Tang ◽  
Maria Vasilyeva ◽  
Eric T. Chung ◽  
Yalchin Efendiev
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Model Problem ◽  
Mixed Finite Element ◽  
Basis Functions ◽  
Fine Grid ◽  
Perforated Domains ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Element Method

In this paper, we consider a coupled system of equations that describes simplified magnetohydrodynamics (MHD) problem in perforated domains. We construct a fine grid that resolves the perforations on the grid level in order to use a traditional approximation. For the solution on the fine grid, we construct approximation using the mixed finite element method. To reduce the size of the fine grid system, we will develop a Mixed Generalized Multiscale Finite Element Method (Mixed GMsFEM). The method differs from existing approaches and requires some modifications to represent the flow and magnetic fields. Numerical results are presented for a two-dimensional model problem in perforated domains. This model problem is a special case for the general 3D problem. We study the influence of the number of multiscale basis functions on the accuracy of the method and show that the proposed method provides a good accuracy with few basis functions.

scholarly journals Generalized Multiscale Finite Element Method for Elastic Wave Propagation in the Frequency Domain

Computation ◽  
2020 ◽  
Vol 8 (3) ◽  
pp. 63 ◽  
Author(s):  
Uygulana Gavrilieva ◽  
Maria Vasilyeva ◽  
Eric T. Chung
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Wave Propagation ◽  
Elastic Wave ◽  
Coarse Grid ◽  
Basis Functions ◽  
Interface Conditions ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Element Method

In this work, we consider elastic wave propagation in fractured media. The mathematical model is described by the Helmholtz problem related to wave propagation with specific interface conditions (Linear Slip Model, LSM) on the fracture in the frequency domain. For the numerical solution, we construct a fine grid that resolves all fracture interfaces on the grid level and construct approximation using a finite element method. We use a discontinuous Galerkin method for the approximation by space that helps to weakly impose interface conditions on fractures. Such approximation leads to a large system of equations and is computationally expensive. In this work, we construct a coarse grid approximation for an effective solution using the Generalized Multiscale Finite Element Method (GMsFEM). We construct and compare two types of the multiscale methods—Continuous Galerkin Generalized Multiscale Finite Element Method (CG-GMsFEM) and Discontinuous Galerkin Generalized Multiscale Finite Element Method (DG-GMsFEM). Multiscale basis functions are constructed by solving local spectral problems in each local domains to extract dominant modes of the local solution. In CG-GMsFEM, we construct continuous multiscale basis functions that are defined in the local domains associated with the coarse grid node and contain four coarse grid cells for the structured quadratic coarse grid. The multiscale basis functions in DG-GMsFEM are discontinuous and defined in each coarse grid cell. The results of the numerical solution for the two-dimensional Helmholtz equation are presented for CG-GMsFEM and DG-GMsFEM for different numbers of multiscale basis functions.

scholarly journals Mixed Generalized Multiscale Finite Element Method for Darcy-Forchheimer Model

Mathematics ◽  
2019 ◽  
Vol 7 (12) ◽  
pp. 1212
Author(s):  
Denis Spiridonov ◽  
Jian Huang ◽  
Maria Vasilyeva ◽  
Yunqing Huang ◽  
Eric T. Chung
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Heterogeneous Media ◽  
Mixed Finite Element ◽  
Picard Iteration ◽  
Fine Grid ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Forchheimer Model ◽  
Element Method

In this paper, the solution of the Darcy-Forchheimer model in high contrast heterogeneous media is studied. This problem is solved by a mixed finite element method (MFEM) on a fine grid (the reference solution), where the pressure is approximated by piecewise constant elements; meanwhile, the velocity is discretized by the lowest order Raviart-Thomas elements. The solution on a coarse grid is performed by using the mixed generalized multiscale finite element method (mixed GMsFEM). The nonlinear equation can be solved by the well known Picard iteration. Several numerical experiments are presented in a two-dimensional heterogeneous domain to show the good applicability of the proposed multiscale method.

scholarly journals An Online Generalized Multiscale Finite Element Method for Unsaturated Filtration Problem in Fractured Media

Mathematics ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1382
Author(s):  
Denis Spiridonov ◽  
Maria Vasilyeva ◽  
Aleksei Tyrylgin ◽  
Eric T. Chung
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Model Reduction ◽  
Basis Functions ◽  
Computational Domain ◽  
Filtration Problem ◽  
Spectral Problems ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Element Method

In this paper, we present a multiscale model reduction technique for unsaturated filtration problem in fractured porous media using an Online Generalized Multiscale finite element method. The flow problem in unsaturated soils is described by the Richards equation. To approximate fractures we use the Discrete Fracture Model (DFM). Complex geometric features of the computational domain requires the construction of a fine grid that explicitly resolves the heterogeneities such as fractures. This approach leads to systems with a large number of unknowns, which require large computational costs. In order to develop a more efficient numerical scheme, we propose a model reduction procedure based on the Generalized Multiscale Finite element method (GMsFEM). The GMsFEM allows solving such problems on a very coarse grid using basis functions that can capture heterogeneities. In the GMsFEM, there are offline and online stages. In the offline stage, we construct snapshot spaces and solve local spectral problems to obtain multiscale basis functions. These spectral problems are defined in the snapshot space in each local domain. To improve the accuracy of the method, we add online basis functions in the online stage. The construction of the online basis functions is based on the local residuals. The use of online bases will allow us to get a significant improvement in the accuracy of the method. We present results with different number of offline and online multisacle basis functions. We compare all results with reference solution. Our results show that the proposed method is able to achieve high accuracy with a small computational cost.

AN ENERGY-CONSERVING DISCONTINUOUS MULTISCALE FINITE ELEMENT METHOD FOR THE WAVE EQUATION IN HETEROGENEOUS MEDIA

2011 ◽  
Vol 03 (01n02) ◽  
pp. 251-268 ◽  
Author(s):  
ERIC T. CHUNG ◽  
YALCHIN EFENDIEV ◽  
RICHARD L. GIBSON
Keyword(s):  
Finite Element ◽  
Wave Equation ◽  
Heterogeneous Media ◽  
Coarse Grid ◽  
Multiscale Methods ◽  
Basis Functions ◽  
Coupling Mechanism ◽  
Time Step ◽  
Fine Grid ◽  
Multiscale Finite Element

Seismic data are routinely used to infer in situ properties of earth materials on many scales, ranging from global studies to investigations of surficial geological formations. While inversion and imaging algorithms utilizing these data have improved steadily, there are remaining challenges that make detailed measurements of the properties of some geologic materials very difficult. For example, the determination of the concentration and orientation of fracture systems is prohibitively expensive to simulate on the fine grid and, thus, some type of coarse-grid simulations are needed. In this paper, we describe a new multiscale finite element algorithm for simulating seismic wave propagation in heterogeneous media. This method solves the wave equation on a coarse grid using multiscale basis functions and a global coupling mechanism to relate information between fine and coarse grids. Using a mixed formulation of the wave equation and staggered discontinuous basis functions, the proposed multiscale methods have the following properties. • The total wave energy is conserved. • Mass matrix is diagonal on a coarse grid and explicit energy-preserving time discretization does not require solving a linear system at each time step. • Multiscale basis functions can accurately capture the subgrid variations of the solution and the time stepping is performed on a coarse grid. We discuss various subgrid capturing mechanisms and present some preliminary numerical results.

Sign in / Sign up

Export Citation Format

Share Document