A class of analytical solutions for multidimensional multispecies diffusive transport coupled with precipitation‐dissolution reactions and porosity changes

One of the challenging problems in mathematical geosciences is the determination of analytical solutions of nonlinear partial differential equations describing transport processes in porous media. We are interested in diffusive transport coupled with precipitation‐dissolution reactions. Several numerical computer codes that simulate such systems have been developed. Analytical solutions, if they exist, represent an important tool for verification of numerical solutions. We present a methodology for deriving such analytical solutions that are exact and explicit in space and time variables. They describe transport of several aqueous species coupled to precipitation and dissolution of a single mineral in one, two, and three dimensions. As an application, we consider explicit analytical solutions for systems containing one or two solute species that describe the evolution of solutes and solid concentrations as well as porosity. We use one of the proposed analytical solutions to test numerical solutions obtained from two conceptually different reactive transport codes. Both numerical implementations could be verified with the help of the analytical solutions and show good agreement in terms of spatial and temporal evolution of concentrations and porosities.