Martin Schneider

Abstract

In this paper, we address the expensive computational cost resulting from limited time-step sizes during numerical simulations of two-phase flow in porous media using dynamic pore-network models. To overcome this issue, we propose a numerical method for dynamic pore-network models using a fully implicit approach. The proposed method introduces a regularization strategy considering the historical fluid configuration at the pore throat, which smooths the discontinuities in local conductivity caused by invasion and snap-off events. The results demonstrate the superiority of the proposed method in terms of accuracy, efficiency and consistency in comparison with other numerical schemes. With similar computational cost, determined by time-step sizes and number of Newton iterations, the developed method in this work yields more accurate results compared to similar schemes presented in the literature. Additionally, our results highlight the enhanced robustness of the our scheme, as it exhibits reduced sensitivity to variations in time-step sizes.

Abstract

We present a numerical scheme for the solution of nonlinear mixed-dimensional PDEs describing coupled processes in embedded tubular network system in exchange with a bulk domain. Such problems arise in various biological and technical applications such as in the modeling of root-water uptake, heat exchangers, or geothermal wells. The nonlinearity appears in form of solution-dependent parameters such as pressure-dependent permeability or temperature-dependent thermal conductivity. We derive and analyze a numerical scheme based on distributing the bulk-network coupling source term by a smoothing kernel with local support. By the use of local analytical solutions, interface unknowns and fluxes at the bulk-network interface can be accurately reconstructed from coarsely resolved numerical solutions in the bulk domain. Numerical examples give confidence in the robustness of the method and show the results in comparison to previously published methods. The new method outperforms these existing methods in accuracy and efficiency. In a root water uptake scenario, we accurately estimate the transpiration rate using only a few thousand 3D mesh cells and a structured cube grid whereas other state-of-the-art numerical schemes require millions of cells and local grid refinement to reach comparable accuracy.

Abstract

Flow in fractured porous media is of high relevance in a variety of geotechnical applications, given the fact that they ubiquitously occur in nature and that they can have a substantial impact on the hydraulic properties of rock. As a response to this, an active field of research has developed over the past decades, focusing on the development of mathematical models and numerical methods to describe and simulate the flow processes in fractured rock. In this work, we present a comparison of different finite volume-based numerical schemes for the simulation of flow in fractured porous media by means of numerical experiments. Two novel vertex-centered approaches are presented and compared to well-established numerical schemes in terms of convergence behavior and their performance on benchmark studies taken from the literature. The new schemes show to produce results that are in good agreement with those of established methods while leading to less degrees of freedom on unstructured simplex grids.

Abstract

We present version 3 of the open-source simulator for flow and transport processes in porous media DuMux. DuMux is based on the modular C++ framework Dune (Distributed and Unified Numerics Environment) and is developed as a research code with a focus on modularity and reusability. We describe recent efforts in improving the transparency and efficiency of the development process and community-building, as well as efforts towards quality assurance and reproducible research. In addition to a major redesign of many simulation components in order to facilitate setting up complex simulations in DuMux, version 3 introduces a more consistent abstraction of finite volume schemes. Finally, the new framework for multi-domain simulations is described, and three numerical examples demonstrate its flexibility.

Abstract

We present a new method for modeling tissue perfusion on the capillary scale. The microvasculature is represented by a network of one-dimensional vessel segments embedded in the extra-vascular space. Vascular and extra-vascular space exchange fluid over the vessel walls. This exchange is modeled by distributed sources using smooth kernel functions for the extra-vascular domain. It is shown that the proposed method may significantly improve the approximation of the exchange flux, in comparison with existing methods for mixed-dimension embedded problems. Furthermore, the method exhibits better convergence rates of the relevant quantities due to the increased regularity of the extra-vascular pressure solution. Numerical experiments with a vascular network from the rat cortex show that the error in the approximation of the exchange flux for coarse grid resolutions may be decreased by a factor of 3. This may open the way for computing on larger network domains, where a fine grid resolution cannot be achieved in practical simulations due to constraints in computational resources, for example in the context of uncertainty quantification.

Abstract

In this work we present an abstract finite volume discretization framework for incompressible immiscible two-phase flow through porous media. A priori error estimates are derived that allow us to prove the existence of discrete solutions and to establish the proof of convergence for schemes belonging to this framework. In contrast to existing publications the proof is not restricted to a specific scheme and it assumes neither symmetry nor linearity of the flux approximations. Two nonlinear schemes, namely a nonlinear two-point flux approximation and a nonlinear multipoint flux approximation, are presented, and some properties of these schemes, e.g. saturation bounds, are proven. Furthermore, the numerical behavior of these schemes (e.g. accuracy, coercivity, efficiency or saturation bounds) is investigated for different test cases for which the coercivity is checked numerically.

Abstract

In this paper, we present a fast streamline-based numerical method for the two-phase flow equations in high-rate flooding scenarios for incompressible fluids in heterogeneous and anisotropic porous media. A fractional flow formulation is adopted and a discontinuous Galerkin method (DG) is employed to solve the pressure equation. Capillary effects can be neglected in high-rate flooding scenarios. This allows us to present an improved streamline approach in combination with the one-dimensional front tracking method to solve the transport equation. To handle the high computational costs of the DG approximation, domain decomposition is applied combined with an algebraic multigrid preconditioner to solve the linear system. Special care at the interior interfaces is required and the streamline tracer has to include a dynamic communication strategy. The method is validated in various two- and three-dimensional tests, where comparisons of the solutions in terms of approximation of flow front propagation with standard fully implicit finite-volume methods are provided.