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

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

A Petrov-Galerkin reduced basis approximation of the Stokes equation in parameterized geometries

We present a Petrov-Galerkin reduced basis (RB) approximation for the parameterized Stokes equation. Our method, which relies on a reduced solution space and a parameter-dependent test space, is shown to be stable (in the sense of Babuska) and algebraically stable '(a bound on the condition number of the online system can be established). Compared to other stable RB methods that can also be shown to be algebraically stable, our approach is among those with the smallest online time cost and it has general applicability to linear non-coercive problems without assuming a saddle-point structure. (C) 2015 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved

An adaptive finite element heterogeneous multiscale method for Stokes flow in porous media

A finite element heterogeneous multiscale method is proposed for solving the Stokes problem in porous media. The method is based on the coupling of an effective Darcy equation on a macroscopic mesh with unknown permeabilities recovered from micro finite element calculations for Stokes problems on sampling domains centered at quadrature points in each macro element. The numerical method accounts for nonperiodic microscopic geometry that can be obtained from a smooth deformation of a reference pore sampling domain. The computational work is nevertheless independent of the small size of the pore structure. A priori error estimates reveal that the overall accuracy of the numerical scheme is limited by the regularity of the solutions of the Stokes microproblems. This regularity is low for a typical situation of nonconvex microscopic pore geometries. We therefore propose an adaptive scheme with micro- macro mesh refinement driven by residual-based indicators that quantify both the macro- and microerrors. A posteriori error analysis is derived for the new method. Two- and three-dimensional numerical experiments confirm the robustness and the accuracy of the adaptive method

A discontinuous Galerkin reduced basis numerical homogenization method for fluid flow in porous media

We present a new conservative multiscale method for Stokes flow in heterogeneous porous media. The method couples a discontinuous Galerkin finite element method (DG-FEM) at the macroscopic scale for the solution of an effective Darcy equation with a Stokes solver at the pore scale to recover effective permeabilities at macroscopic quadrature points. To avoid costly computation of numerous Stokes problems throughout the macroscopic computational domain, the pore geometry is parametrized and a model order reduction algorithm is used to select representative microscopic geometries. Accurate Stokes solutions and related permeabilities are obtained for these representative geometries in an offline stage. In an online stage, the DG-FEM is computed with permeabilities recovered at the desired macroscopic quadrature points from the precomputed Stokes solutions. The multiscale method is shown to be mass conservative at the macro scale and the computational cost for the online stage is similar to the cost of solving a single scale Darcy problem. Numerical experiments for two and three dimensional problems illustrate the efficiency and the performance of the proposed method

MATHICSE Technical Report : Multiscale model reduction methods for flow in heterogeneous porous media

In this paper we provide a general framework for model reduction methods applied to fluid flow in porous media. Using reduced basis and numerical homogenization techniques we show that the complexity of the numerical approximation of Stokes flow in heterogeneous media can be drastically reduced. The use of such acomputational framework is illustrated at several model problems such as two and three scale porous media

MATHICSE Technical Report : A reduced basis finite element heterogeneous multiscale method for Stokes flow in porous media

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 [A. Abdulle, O. Budáč, Multiscale Model. Simul. 13 (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

MATHICSE Technical Report : Multiscale methods and model order reduction for flow problems in three-scale porous media

A new multiscale method combined with model order reduction is proposed for flow problems in three-scale porous media. We derive an effective three-scale model that couples a macroscopic Darcy equation, a mesoscopic Stokes-Brinkman equation, and a microscopic Stokes equation. A corresponding three-scale numerical method is then derived using the finite element discretization with numerical quadrature, where the macroscopic and mesoscopic permeability is upscaled at quadrature points from mesoscopic and microscopic problems, respectively. The computational cost of solving numerous mesoscopic and microscopic flow problems is further reduced by applying a Petrov–Galerkin reduced basis method at the mesocopic and microscopic scales. As there is nonatural way to obtain an affine decomposition of the mesoscopic problems, which is instrumental for the efficiency of the model order reduction, we derive a mesoscopic solver that makes use of empirical interpolation techniques. A priori and a posteriori error estimates are derived for the new method that is also tested numerically to corroborate the theoretical convergence rates and illustrate its efficiency