ABSTRACT In this article, we propose and analyze an efficient projection‐based structure‐preserving scheme for Patlak–Keller–Segel equations on surfaces. The mass lumping surface finite element method is adopted for the spatial discretization. The time integration is based on the second‐order Crank–Nicolson/Adams–Bashforth scheme via an projection post‐processing strategy and scalar auxiliary variable approach, which is proven to be bound/positivity preserving, mass conservative, and energy stable with the modified energy. Rigorous convergence analysis is presented, achieving an optimal spatial convergence rate. Various numerical examples are performed to verify these theoretical results and demonstrate the efficiency and robustness of the proposed scheme.
Li et al. (Sun,) studied this question.