Statistical Effective Diffusivity Estimation in Porous Media Using an Integrated On-site Imaging Workflow for Synchrotron Users

Transport in porous media plays an essential role for many physical, engineering, biological and environmental processes. Novel synchrotron imaging techniques and image-based models have enabled more robust quantification of geometric structures that influence transport through the pore space. However, image-based modelling is computationally expensive, and end users often require, while conducting imaging campaign, fast and agile bulk-scale effective parameter estimates that account for the pore-scale details. In this manuscript we enhance a pre-existing image-based model solver known as OpenImpala to estimate bulk-scale effective transport parameters. In particular, the boundary conditions and equations in OpenImpala were modified in order to estimate the effective diffusivity in an imaged system/geometry via a formal multi-scale homogenisation expansion. Estimates of effective pore space diffusivity were generated for a range of elementary volume sizes to estimate when the effective diffusivity values begin to converge to a single value. Results from OpenImpala were validated against a commercial finite element method package COMSOL Multiphysics (abbreviated as COMSOL). Results showed that the effective diffusivity values determined with OpenImpala were similar to those estimated by COMSOL. Tests on larger domains comparing a full image-based model to a homogenised (geometrically uniform) domain that used the effective diffusivity parameters showed differences below 2 % error, thus verifying the accuracy of the effective diffusivity estimates. Finally, we compared OpenImpala's parallel computing speeds to COMSOL. OpenImpala consistently ran simulations within fractions of minutes, which was two orders of magnitude faster than COMSOL providing identical supercomputing specifications. In conclusion, we demonstrated OpenImpala's utility as part of an on-site tomography processing pipeline allowing for fast and agile assessment of porous media processes and to guide imaging campaigns while they are happening at synchrotron beamlines.

Automated summary

This paper presents an image-based model solver called OpenImpala to estimate bulk-scale effective transport parameters in porous media. By modifying the boundary conditions and equations in OpenImpala, a formal multi-scale homogenisation expansion method is used to estimate the effective diffusivity in an imaged system/geometry. The results show that OpenImpala has faster computing speeds compared to COMSOL and accurately estimates the effective diffusivity values.