A posteriori-based, local multilevel mesh refinement for the Darcy porous media flow problem

In this work we develop an a posteriori-based adaptive local multilevel mesh refinement strategy for the Darcy porous media flow problem discretized with the finite volume method. Our approach is based on coupling the FIC "flux interface correction" method Angot et al. (1992) and equilibrated flux a posteriori error estimates Vohralík et al. (2018). The choice of the FIC method is motivated by the fact that it is well-adapted to problems discretized in conservative form since that the FIC provides corrections by balancing fluxes computed from both coarse and fin grids across the interface. The equilibrated flux a posteriori estimates provide a guaranteed upper bound on the total error in the fluxes and take advantage of the conservative scheme in such a way that the estimates can be easily coded and cheaply evaluated. We use a posteriori estimates to detect automatically the zone of interest where the error is dominant in order to defined the local sub-grids. Numerous numerical experiments on practical problems illustrate the performance of our methodology. A comparison with a classical h-adaptive method is also carried out in this work.