An Open-Source MATLAB Implementation for elastic analysis of heterogeneous materials using the Extended Multiscale Finite Element Method

2021 ◽  
Author(s):  
Hui Liu ◽  
Dongsheng Yang ◽  
Jingkai Wu ◽  
Yonggang Zheng ◽  
Hongwu Zhang
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Open Source ◽  
Heterogeneous Materials ◽  
Elastic Analysis ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Matlab Implementation ◽  
Element 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.

Sign in / Sign up

Export Citation Format

Share Document