Reactive transport processes in a porous medium will often both cause changes to the pore structure, via precipitation and dissolution of biomass or minerals, and be affected by these changes, via changes to the material's porosity and permeability. An understanding of the pore structure morphology and the changes to flow parameters during these processes is critical when modeling reactive transport. Commonly applied porosity-permeability relations in simulation models on the REV scale use a power-law relation, often with slight modifications, to describe such features; they are often used for modeling the effects of mineral precipitation and/or dissolution on permeability. To predict the reduction in permeability due to biomass growth, many different and often rather complex relations have been developed and published by a variety of authors. Some authors use exponential or simplified Kozeny-Carman relations. However, many of these relations do not lead to fundamentally different predictions of permeability alteration when compared to a simple power-law relation with a suitable exponent. Exceptions to this general trend are only few of the porosity-permeability relations developed for biomass clogging; these consider a residual permeability even when the pore space is completely filled with biomass. Other exceptions are relations that consider a critical porosity at which the porous medium becomes impermeable; this is often used when modeling the effect of mineral precipitation. This review first defines the scale on which porosity-permeability relations are typically used and aims at explaining why these relations are not unique. It shows the variety of existing approaches and concludes with their essential features.

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

We present a fully coupled soil-atmosphere model that includes radiation in the energy balance of the coupling conditions between the two domains. The model is able to describe evaporation processes under the influence of turbulence, surface roughness, and soil heterogeneities and focuses specifically on the influence of radiation on the mass and energy transport across the soil-atmosphere interface. It is shown that evaporation rates are clearly dominated by the diurnal cycle of solar irradiance. During Stage-I evaporation maximum temperatures are regulated due to evaporative cooling, but after a transition into Stage-II evaporation, temperatures rise tremendously. We compare two different soil types, a coarser, sandy soil and a finer, silty soil, and analyze evaporation rates, surface temperatures, and net radiation for three different wind conditions. The influence of surface undulations on radiation and evaporation is analyzed and shows that radiation can lead to different local drying patterns in the hills and the valleys of the porous medium, depending on the height of the undulations and on the direction of the Sun. At last a comparison of lysimeter measurement data to the numerical examples shows a good match for measured and calculated radiation values but evaporation rates are still overestimated in the model. Possible reasons for the discrepancy between measurement and model data are analyzed and are found to be uncertainties about the parameters close to the interface, which are decisive for determining evaporation rates. Several models are available that are able to describe evaporation with different levels of complexity, especially regarding simplifications made for the coupling conditions between the soil and the atmosphere.