A Family of Block-Centered Schemes for Contaminant Transport Equations with Adsorption via Integral Method with Variational Limit

Abstract

This paper develops a class of high-order conservative schemes for contaminant transport with equilibrium adsorption, based on the Integral Method with Variational Limit on block-centered grids. By incorporating four parameters, the scheme can reproduce classical fourth-order compact schemes and further extend to sixth- and eighth-order accurate formulations, all within a unified framework. Under periodic boundary conditions, we analyze the stability, convergence, and mass conservation of the parameterized numerical scheme. Numerical experiments are then conducted to examine the impact of parameter variations on errors, explore the relationship between parameters and the fourth-, sixth-, and eighth-order schemes, and verify that the schemes' high-order accuracy aligns with theoretical predictions. To enhance the applicability of the proposed method, we further develop two fourth-order compact boundary treatments that ensure uniform accuracy between boundary and interior regions. Numerical results confirm the effectiveness of the proposed schemes across various adsorption models.