A Modified Reaction–Diffusion–Advection Framework for Subsurface Co₂ Sequestration: Derivation, Non-Dimensionalization, and DuFort–Frankel Discretization of a Porous-Media Transport Model
Owuor Lucas Otieno *
Department of Mathematics, Kibabii University, Bungoma, Kenya.
Elijah Sifuna Wafula
Department of Mathematics, Kibabii University, Bungoma, Kenya.
*Author to whom correspondence should be addressed.
Abstract
Geological storage of CO₂ is an important approach for long-term climate change mitigation, but reliable prediction of subsurface plume behaviour requires mathematical models that clearly link physical assumptions to governing equations and numerical implementation. This study presents a derivation-only framework for a modified reaction–diffusion–advection model describing one-dimensional CO₂ concentration dynamics in a porous geological formation. Starting from mass conservation, the formulation incorporates pressure-driven advection, effective diffusion, linearised geochemical removal, injection input and caprock leakage in a single governing partial differential equation. The model is then expressed in dimensionless form to identify the Péclet number, Damköhler number and leakage number as controlling parameters for the relative influence of advective migration, diffusive spreading, geochemical trapping and containment loss. A DuFort–Frankel finite-difference discretisation is subsequently derived to obtain an explicit three-level update formula for numerical implementation. The resulting formulation provides a transparent theoretical basis for interpreting the roles of transport, reaction and leakage in simplified geological carbon storage systems. As the study does not include numerical experiments, benchmark comparison or field validation, the framework should be regarded as a foundation for future computational testing rather than a validated predictive tool.
Keywords: Carbon dioxide sequestration, Reaction–diffusion–advection equation, DuFort–Frankel scheme, Damköhler number, Péclet number, porous media, geological storage, CCS, numerical methods