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.

No abstract

Microscopic velocity field measurements inside a regular porous medium adjacent to a low Reynolds number channel flow. Physics of Fluids, 31(4), [042001].

Coupled system of free flow above a permeable medium appear in a wide range of technical and environmental applications. Prominent examples including two-phase flow range from fuel-cell water manage-

Modeling coupled systems of free flow adjacent to a porous medium by means of fully resolved Navier–Stokes equations is limited by the immense computational cost and is thus only feasible for relatively small domains. Coupled, hybrid-dimensional models can be much more efficient by simplifying the porous domain, e.g., in terms of a pore-network model. In this work, we present a coupled pore-network/free-flow model taking into account pore-scale slip at the local interfaces between free flow and the pores. We consider two-dimensional and three-dimensional setups and show that our proposed slip condition can significantly increase the coupled model’s accuracy: compared to fully resolved equidimensional numerical reference solutions, the normalized errors for velocity are reduced by a factor of more than five, depending on the flow configuration. A pore-scale slip parameter $$\beta _{{{{\rm pore}}}}$$ β pore required by the slip condition was determined numerically in a preprocessing step. We found a linear scaling behavior of $$\beta _{{{{\rm pore}}}}$$ β pore with the size of the interface pore body for three-dimensional and two-dimensional domains. The slip condition can thus be applied without incurring any run-time cost. In the last section of this work, we used the coupled model to recalculate a microfluidic experiment where we additionally exploited the flat structure of the micromodel which permits the use of a quasi-3D free-flow model. The extended coupled model is accurate and efficient.

No abstract

Carbon dioxide density-driven dissolution in a water-filled laboratory flume of the dimensions 60 cm length, 40 cm height, 1 cm thickness, was visualized using a pH-sensitive color indicator. We focus on atmospheric pressure conditions, like in caves where CO 2 concentrations are typically higher. Varying concentrations of carbon dioxide were applied as boundary conditions at the top of the experimental setup, leading to the onset of convective fingering at differing times. The data were used to validate a numerical model implemented in the numerical simulator DuMu x . The model solves the Navier-Stokes equations for density-induced water flow with concentration-dependent fluid density and a transport equation, including advective and diffusive processes for the carbon dioxide dissolved in water. The model was run in 2D, 3D, and pseudo-3D on two different grids. Without any calibration or fitting of parameters, the results of the comparison between experiment and simulation show satisfactory agreement with respect to the onset time of convective fingering, and the number and the dynamics of the fingers. Grid refinement matters, in particular, in the uppermost part where fingers develop. The 2D simulations consistently overestimated the fingering dynamics. This successful validation of the model is the prerequisite for employing it in situations with background flow and for a future study of karstification mechanisms related to CO 2 -induced fingering in caves.

A discretization is proposed for models coupling free flow with anisotropic porous medium flow.Our approach employs a staggered grid finite volume method for the Navier-Stokes equations in the free flow subdomain and a MPFA finite volume method to solve Darcy flow in the porous medium. After appropriate spatial refinement in the free flow domain, the degrees of freedom are conveniently located to allow for a natural coupling of the two discretization schemes. In turn, we automatically obtain a more accurate description of the flow field surrounding the porous medium.Numerical experiments highlight the stability and applicability of the scheme in the presence of anisotropy and show good agreement with existing methods, verifying our approach.