DuMux 3 – an open-source simulator for solving flow and transport problems in porous media with a focus on model coupling

Abstract

We present version 3 of the open-source simulator for flow and transport processes in porous media DuMu x . DuMu x 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 DuMu x , 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.

Experimental and Simulation Study on Validating a Numerical Model for CO2 Density-Driven Dissolution in Water

“…On the other hand, the buffering by the color indicator should act protractingly. A third point: as we have shown, the implementation of the discretization scheme (see Section 2.2.2) [20,21] can cause a delay of the onset time for coarse grids, since it requires some time to have the first discrete finite volume at the top of the water body at a high enough CO 2 concentration to trigger instability. While it is not easy to separate these effects quantitatively, we can still state that the agreement between experiment and simulations is very reasonable and the tendency towards longer onset times for lower CO 2 concentrations is consistently reproduced in all cases.…”

“…The Newton scheme chooses the time-step size according to its convergence with a user-controlled maximum time-step size. For further details on discretization, numerical solution methods, and their implementation, we refer to [20] or the handbook of DuMu x [21].…”

Call for Participation: Collaborative Benchmarking of Functional-Structural Root Architecture Models. The Case of Root Water Uptake

“…The flux is added as a source term in the corresponding segment in the root. The model is implemented in the open-source porous media simulator DuMu x (Flemisch et al, 2011;Koch et al, 2018a;Koch et al, 2019). The coupled system is solved with a fully coupled manner, using Newton's method, and monolithic linear solver (block-preconditioned stabilized bi-conjugate gradient solver) in each Newton iteration.…”

Modeling of methane migration from gas wellbores into shallow groundwater at basin scale

“…A fully implicit Euler time discretization scheme using the Newton method is applied to solve the linearized system of non-linear partial differential equations (Helmig 1997). The first time-step size is set to 0.01 s and the maximum time-step is 1e6 s. The time-step size is controlled by the number of iterations required by the Newton method to achieve convergence for the last time integration (Flemisch et al 2011;Koch et al 2020).…”