An adaptive hybridizable discontinuous Galerkin method for Darcy-Forchheimer flow in fractured porous media

In this paper, we consider an adaptive hybridizable discontinuous Galerkin (HDG) method based on the discrete fracture model for approximation of a single-phase flow in fractured porous media. We are interested in the case that the flow rate in fracture is large enough to justify the use of Forchheimer's law for modeling flow within the fracture, while Darcy's law is applied to the surrounding matrix. The HDG method could be designed to simulate the flow in porous media with reduced fractures which consist of many straight lines or planes. More specifically, we use piecewise polynomials of degree k to approximate the velocity and pressure in fracture and surrounding porous media. The existence and uniqueness of discrete solutions are proved by the Brouwer fixed point theorem, and an efficient and reliable a posteriori error estimator is obtained with respect to an energy norm. Moreover, the HDG scheme, the existence and uniqueness of discrete solutions, and the a posteriori error estimates are also extended to the problem with non-planar, embedded, and intersecting fractures. Finally, several numerical examples are provided to validate the performance of the obtained a posteriori error estimator.