Rigorous vertical averaging for pseudo-3D modeling of precipitation obstructed fluid flow
Abstract
Computing effective hydraulic properties of porous materials with a dynamically changing pore morphology presents a complex and critical challenge in both theoretical and applied sciences. This paper presents a tool to effectively compute the intrinsic permeability k(x,t) for domains with both temporally and spatially varying pore structures, by utilizing a sample of snapshots of microfluidic experiments with calcium carbonate precipitation. To compute large domains at reasonable cost, a pseudo-3D Stokes solver with an additional viscous drag term is employed. In order to simulate the intricate domains with highly heterogeneous pore space, we derive vertically averaged governing equations including a vertically averaged drag term from scratch and without limiting assumptions. For exemplary precipitates (semi-spheres) of varying sizes within a rectangular channel with fixed solid boundaries, results are compared across multiple modeling approaches, each employing different drag force formulations. In addition, the performance of the derived models is examined for segments of an EICP experiment in terms of accuracy (comparison with 3D solvers) and computational efficiency. Finally, an entire experimental domain is used to demonstrate the capabilities offered by the solver thanks to its efficiency as well as accuracy and the permeability tensor k is computed and compared with experimental results. In all these applications, the vertically averaged drag term proves to be the optimal solution in terms of speed and accuracy, thus allowing for an efficient analysis of the non-uniformity of local fluid flow fluctuations resulting from alterations in the pore space.
Related
- Trade-offs between predictive accuracy and distributional fidelity in surrogate modeling of extreme precipitation
- Investigating hydrological modeling uncertainties in the Mediterranean region by combining precipitation and soil moisture products
- Added value of integrated satellite precipitation products for event-based flood modeling in semi-arid regions
- Increasing daily precipitation extremes despite declining annual totals in southern Europe: a modeling study on the effects of Mediterranean Sea warming
- Improved GPP modeling reveals rising precipitation influence in typical grasslands of Northern China's Hulun Lake Basin.
BibTeX
@misc{krach2026rigorous,
title = {{Rigorous vertical averaging for pseudo-3D modeling of precipitation obstructed fluid flow}},
author = {David Krach and Felix Weinhardt and Holger Class and Holger Steeb and Martin Schneider},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.39391},
eprint = {2609.39391},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.39391},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Rigorous vertical averaging for pseudo-3D modeling of precipitation obstructed fluid flow
AU - David Krach
AU - Felix Weinhardt
AU - Holger Class
AU - Holger Steeb
AU - Martin Schneider
PY - 2026
DA - 2026/09/30
PB - arXiv
DO - 10.48550/arxiv.2609.39391
UR - https://arxiv.org/abs/2609.39391
AB - Computing effective hydraulic properties of porous materials with a dynamically changing pore morphology presents a complex and critical challenge in both theoretical and applied sciences. This paper presents a tool to effectively compute the intrinsic permeability k(x,t) for domains with both temporally and spatially varying pore structures, by utilizing a sample of snapshots of microfluidic experiments with calcium carbonate precipitation. To compute large domains at reasonable cost, a pseudo-3D Stokes solver with an additional viscous drag term is employed. In order to simulate the intricate domains with highly heterogeneous pore space, we derive vertically averaged governing equations including a vertically averaged drag term from scratch and without limiting assumptions. For exemplary precipitates (semi-spheres) of varying sizes within a rectangular channel with fixed solid boundaries, results are compared across multiple modeling approaches, each employing different drag force formulations. In addition, the performance of the derived models is examined for segments of an EICP experiment in terms of accuracy (comparison with 3D solvers) and computational efficiency. Finally, an entire experimental domain is used to demonstrate the capabilities offered by the solver thanks to its efficiency as well as accuracy and the permeability tensor k is computed and compared with experimental results. In all these applications, the vertically averaged drag term proves to be the optimal solution in terms of speed and accuracy, thus allowing for an efficient analysis of the non-uniformity of local fluid flow fluctuations resulting from alterations in the pore space.
ER -