Multiresolution Operator Decomposition for Flow Simulation in Fractured Porous Media.

Qingfu Zhang,Houman Owhadi,Jun Yao,Florian Schafer,Zhaoqin Huang,Yang Li

DOI: https://doi.org/10.1016/j.jcp.2018.12.032

IF: 4.645

2019-01-01

Journal of Computational Physics

Abstract:Fractures should be simulated accurately given their significant effects on whole flow patterns in porous media. But such high-resolution simulations impose severe computational challenges to numerical methods in the applications. Therefore, the demand for accurate and efficient coarse-graining techniques is increasing. In this work, a near-linear complexity multiresolution operator decomposition method is proposed for solving and coarse graining flow problems in fractured porous media. We use the Discrete Fracture Model (DFM) to describe fractures, in which the fractures are explicitly represented as (n−1)-dimensional elements. Using operator adapted wavelets, the solution space is decomposed into subspaces where DFM subsolutions can be computed by solving sparse and well-conditioned linear systems. By keeping only the coarse-scale part of the solution space, we furthermore obtain a reduced order model. We provide numerical experiments that investigate the accuracy of the reduced order model for different resolutions and different choices of medium.