  We present a versatile open-source pipeline for simulating inhomogeneous
reaction-diffusion processes in highly resolved, image-based geometries of
porous media with reactive boundaries. Resolving realistic pore-scale
geometries in numerical models is challenging and computationally demanding, as
the scale differences between the sizes of the interstitia and the whole system
can lead to prohibitive memory requirements. The present pipeline combines a
level-set method with geometry-adapted sparse block grids on GPUs to
efficiently simulate reaction-diffusion processes in image-based geometries. We
showcase the method by applying it to fertilizer diffusion in soil, heat
transfer in porous ceramics, and determining effective diffusion coefficients
and tortuosity. The present approach enables solving reaction-diffusion partial
differential equations in real-world geometries applicable to porous media
across fields such as engineering, environmental science, and biology.
