Multiscale extended finite element method for deformable fractured porous media

被引:17
|
作者
Xu, Fanxiang [1 ]
Hajibeygi, Hadi [2 ]
Sluys, Lambertus J. [1 ]
机构
[1] Delft Univ Technol, Fac Civil Engn & Geosci, Mat Mech Management & Design, POB 5048, NL-2600 GA Delft, Netherlands
[2] Delft Univ Technol, Fac Civil Engn & Geosci, Dept Geosci & Engn, POB 5048, NL-2600 GA Delft, Netherlands
基金
美国国家科学基金会;
关键词
Fractured porous media; Extended finite element; Multiscale; Geomechanics; Scalable iterative solver; CRACK-GROWTH; VOLUME METHOD; FLOW; MODEL; SIMULATION; HOMOGENIZATION; SOLVER; XFEM; PROPAGATION; FORMULATION;
D O I
10.1016/j.jcp.2021.110287
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Deformable fractured porous media appear in many geoscience applications. While the extended finite element method (XFEM) has been successfully developed within the computational mechanics community for accurate modeling of deformation, its application in geoscientific applications is not straightforward. This is mainly due to the fact that subsurface formations are heterogeneous and span large length scales with many fractures at different scales. To resolve this limitation, in this work, we propose a novel multiscale formulation for XFEM, based on locally computed enriched basis functions. The local multiscale basis functions capture heterogeneity of th e porous rock properties, and discontinuities introduced by the fractures. In order to preserve accuracy of these basis functions, reduced-dimensional boundary conditions are set as localization condition. Using these multiscale bases, a multiscale coarse-scale system is then governed algebraically and solved. The coarse scale system entails no enrichment due to the fractures. Such formulation allows for significant computational cost reduction, at the same time, it preserves the accuracy of the discrete displacement vector space. The coarse-scale solution is finally interpolated back to the fine scale system, using the same multiscale basis functions. The proposed multiscale XFEM (MS-XFEM) is also integrated within a two stage algebraic iterative solver, through which error reduction to any desired level can be achieved. Several proof-of-concept numerical tests are presented to assess the performance of the developed method. It is shown that the MS-XFEM is accurate, when compared with the fine-scale reference XFEM solutions. At the same time, it is significantly more efficient than the XFEM on fine-scale resolution, as it significantly reduces the size of the linear systems. As such, it develops a promising scalable XFEM method for large-scale heavily fractured porous media. (c) 2021 The Author(s). Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). <comment>Superscript/Subscript Available</comment
引用
收藏
页数:19
相关论文
共 50 条