An open-source pipeline for solving continuous reaction-diffusion models in image-based geometries of porous media

Justina Stark,Ivo F. Sbalzarini
DOI: https://doi.org/10.1016/j.jocs.2023.102118
2023-10-16
Abstract: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.
Mathematical Software,Computational Physics
What problem does this paper attempt to address?
### Problems the paper attempts to solve The paper aims to address the challenges of simulating heterogeneous reaction - diffusion processes in porous media. Specifically, the paper proposes an open - source pipeline for efficiently solving continuous reaction - diffusion models in image - based porous media geometries. These models usually have complex and irregular geometries, which makes it very difficult to construct accurate numerical models. The main problems include: 1. **Modeling of complex geometries**: - Geometries in porous media are usually highly irregular and cannot be described by global parameterization. Therefore, algorithmic surface representation methods, such as surface triangulation or level - set methods, need to be used. 2. **Memory requirements for multi - scale simulations**: - Since the scale of the pore space in porous media differs by several orders of magnitude from the scale of the entire system, multi - scale simulations usually require a large amount of computer memory, limiting their feasibility and efficiency. 3. **Efficient numerical methods**: - Traditional numerical methods (such as the finite element method and the finite volume method) have high computational costs when dealing with complex geometries. Therefore, a new method that combines sparse block grids and the finite difference method needs to be developed to improve memory and computational efficiency. 4. **Validation in practical applications**: - This method needs to be validated in practical applications, such as the diffusion of fertilizers in soil, heat transfer in porous ceramic catalysts, etc., to prove its effectiveness and practicality. ### Solutions To address the above challenges, the paper proposes the following solutions: 1. **Level - set method**: - Use the level - set method to represent the complex geometries of porous media. This method can implicitly describe arbitrary geometries without the need for explicit surface meshes. 2. **Sparse block grid**: - Combine the sparse block grid data structure and allocate discrete points only in the regions of interest, thereby significantly reducing memory requirements and computational costs. 3. **GPU acceleration**: - Utilize GPU parallel computing technology to accelerate the simulation process through the OpenFPM framework and improve computational efficiency. 4. **Complete simulation process**: - Provide a complete workflow from 3D sample image segmentation to computational representation generation, then to solving PDEs and result visualization. ### Application examples The paper demonstrates the effectiveness of this method through two practical application examples: 1. **Heterogeneous diffusion of fertilizers in soil**: - Simulate the diffusion process of fertilizers in soil, taking into account the complex geometries of pores in the soil. 2. **Heat transfer in porous ceramic catalysts**: - Simulate the heat transfer process in porous ceramic catalysts, which is of great significance for the application of renewable energy technologies. ### Summary The paper proposes an efficient and accurate open - source pipeline for simulating heterogeneous reaction - diffusion processes in image - based porous media geometries. By combining the level - set method, sparse block grid, and GPU acceleration, this method can handle complex geometries and significantly reduce memory and computational requirements. The paper also validates the effectiveness and practicality of this method through practical application examples.