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

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.

Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 298
Author(s):  
Aleksei Tyrylgin ◽  
Maria Vasilyeva ◽  
Dmitry Ammosov ◽  
Eric T. Chung ◽  
Yalchin Efendiev

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.


Mathematics ◽  
2019 ◽  
Vol 7 (12) ◽  
pp. 1212
Author(s):  
Denis Spiridonov ◽  
Jian Huang ◽  
Maria Vasilyeva ◽  
Yunqing Huang ◽  
Eric T. Chung

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.


Mathematics ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1382
Author(s):  
Denis Spiridonov ◽  
Maria Vasilyeva ◽  
Aleksei Tyrylgin ◽  
Eric T. Chung

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.


Computation ◽  
2020 ◽  
Vol 8 (3) ◽  
pp. 63 ◽  
Author(s):  
Uygulana Gavrilieva ◽  
Maria Vasilyeva ◽  
Eric T. Chung

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.


Sign in / Sign up

Export Citation Format

Share Document