Reactive Dissolution in Porous Media

This dataset contains 32 numerical simulations produced by GeoChemFoam. The 2D simulation results are stored at 100 equally spaced time steps, resulting in time-evolving dissolution state maps of size 260 x 260. Each of these simulations conveys a particular rock sample with its own pore structure, as well as a distinct fluid trajectory over the course of the dissolution process. All simulations employ the same flow and reactive transport conditions, where both Péclet and kinetic numbers are equal to 1. Each HDF5 contains a Python dictionary of NumPy arrays, and comprises five different properties: C: the concentration of the acidic solution used in the dissolution; eps: the volume fraction of the pore space occupied by pore in each voxel; Ux: he magnitude and direction of flow of the acidic solution in the horizontal axis of a 2-D Cartesian plane; Uy: same as Ux, but for the vertical axis; time: the total time elapsed (in seconds) of the dissolution process at a particular time step.