Central Schemes for Porous Media Flows

We are concerned with central differencing schemes for solving scalar hyperbolic conservation laws arising in the simulation of multiphase flows in heterogeneous porous media. We compare the Kurganov-Tadmor, 2000 semi-discrete central scheme with the Nessyahu-Tadmor, 1990 central scheme. The KT scheme uses more precise information about the local speeds of propagation together with integration over nonuniform control volumes, which contain the Riemann fans. These methods can accurately resolve sharp fronts in the fluid saturations without introducing spurious oscillations or excessive numerical diffusion. We first discuss the coupling of these methods with velocity fields approximated by mixed finite elements. Then, numerical simulations are presented for two-phase, two-dimensional flow problems in multi-scale heterogeneous petroleum reservoirs. We find the KT scheme to be considerably less diffusive, particularly in the presence of high permeability flow channels, which lead to strong restrictions on the time step selection; however, the KT scheme may produce incorrect boundary behavior

Reactive Flows in Deformable, Complex Media

Many processes of highest actuality in the real life are described through systems of equations posed in complex domains. Of particular interest is the situation when the domain is changing in time, undergoing deformations that depend on the unknown quantities of the model. Such kind of problems are encountered as mathematical models in the subsurface, material science, or biological systems.The emerging mathematical models account for various processes at different scales, and the key issue is to integrate the domain deformation in the multi-scale context. The focus in this workshop was on novel techniques and ideas in the mathematical modelling, analysis, the numerical discretization and the upscaling of problems as described above

Python framework for HP adaptive discontinuous Galerkin methods for two phase flow in porous media

In this paper we present a framework for solving two-phase flow problems in porous media. The discretization is based on a Discontinuous Galerkin method and includes local grid adaptivity and local choice of polynomial degree. The method is implemented using the new Python frontend Dune-FemPy to the open source framework Dune. The code used for the simulations is made available as Jupyter notebook and can be used through a Docker container. We present a number of time stepping approaches ranging from a classical IMPES method to a fully coupled implicit scheme. The implementation of the discretization is very flexible allowing to test different formulations of the two-phase flow model and adaptation strategies

Local-global splitting for spatiotemporal-adaptive multiscale methods

We present a novel spatiotemporal-adaptive Multiscale Finite Volume (MsFV) method, which is based on the natural idea that the global coarse-scale problem has longer characteristic time than the local fine-scale problems. As a consequence, the global problem can be solved with larger time steps than the local problems. In contrast to the pressure-transport splitting usually employed in the standard MsFV approach, we propose to start directly with a local-global splitting that allows to locally retain the original degree of coupling. This is crucial for highly non-linear systems or in the presence of physical instabilities. To obtain an accurate and efficient algorithm, we devise new adaptive criteria for global update that are based on changes of coarse-scale quantities rather than on fine-scale quantities, as it is routinely done before in the adaptive MsFV method. By means of a complexity analysis we show that the adaptive approach gives a noticeable speed-up with respect to the standard MsFV algorithm. In particular, it is efficient in case of large upscaling factors, which is important for multiphysics problems. Based on the observation that local time stepping acts as a smoother, we devise a self-correcting algorithm which incorporates the information from previous times to improve the quality of the multiscale approximation. We present results of multiphase flow simulations both for Darcy-scale and multiphysics (hybrid) problems, in which a local pore-scale description is combined with a global Darcy-like description. The novel spatiotemporal-adaptive multiscale method based on the local-global splitting is not limited to porous media flow problems, but it can be extended to any system described by a set of conservation equations

Python Framework for HP Adaptive Discontinuous Galerkin Method for Two Phase Flow in Porous Media

In this paper we present a framework for solving two phase flow problems in porous media. The discretization is based on a Discontinuous Galerkin method and includes local grid adaptivity and local choice of polynomial degree. The method is implemented using the new Python frontend Dune-FemPy to the open source framework Dune. The code used for the simulations is made available as Jupyter notebook and can be used through a Docker container. We present a number of time stepping approaches ranging from a classical IMPES method to fully coupled implicit scheme. The implementation of the discretization is very flexible allowing for test different formulations of the two phase flow model and adaptation strategies.Comment: Keywords: DG, hp-adaptivity, Two-phase flow, IMPES, Fully implicit, Dune, Python, Porous media. 28 pages, 9 figures, various code snippet

Reactive Flows in Deformable, Complex Media

Many processes of highest actuality in the real life are described through systems of equations posed in complex domains. Of particular interest is the situation when the domain is variable, undergoing deformations that depend on the unknown quantities of the model. Such kind of problems are encountered as mathematical models in the subsurface, or biological systems. Such models include various processes at different scales, and the key issue is to integrate the domain deformation in the multi-scale context. Having this as the background theme, this workshop focused on novel techniques and ideas in the analysis, the numerical discretization and the upscaling of such problems, as well as on applications of major societal relevance today