A reduced basis finite element heterogeneous multiscale method for Stokes flow in porous media

• Published in: Computer Methods in Applied Mechanics and Engineering (ISSN: 0045-7825), vol. 307, p. 1-31
• Lausanne: Elsevier Science Sa, 2016

A reduced basis Darcy-Stokes finite element heterogeneous multiscale method (RB-DS-FE-HMM) is proposed for the Stokes problem in porous media. The multiscale method is based on the Darcy-Stokes finite element heterogeneous multiscale method (DS-FE-HMM) introduced in Abdulle and Budac (2015) that couples a Darcy equation solved on a macroscopic mesh, with missing permeability data extracted from the solutions of Stokes micro problems at each macroscopic quadrature point. To overcome the increasingly growing cost of repeatedly solving the Stokes micro problems as the macroscopic mesh is refined, we parametrize the microscopic solid geometry and approximate the infinite-dimensional manifold of parameter dependent solutions of Stokes problems by a low-dimensional space. This low-dimensional (reduced basis) space is obtained in an offline stage by a greedy algorithm and used in an online stage to compute the effective Darcy permeability at a cost independent of the microscopic mesh. The discretization of the parametrized Stokes problems relies on a Petrov-Galerkin formulation that allows for a stable and fast online evaluation of the required permeabilities. A priori and a posteriori estimates of the RB-DS-FE-HMM are derived and a residual-based adaptive algorithm is proposed. Two- and three-dimensional numerical experiments confirm the accuracy of the RB-DS-FE-HMM and illustrate the speedup compared to the DS-FE-HMM. (C) 2016 Elsevier B.V. All rights reserved.

Keywords: Stokes flow ; Darcy equation ; Numerical homogenization ; Reduced basis ; Adaptivity

Reference

• EPFL-ARTICLE-210480
• doi:10.1016/j.cma.2016.03.016
• View record in Web of Science

Record created on 2015-08-11, modified on 2017-05-12