Core Concepts
The authors propose a novel approach to construct approximations for the Poisson-Nernst-Planck equations focusing on positivity preservation and mass conservation. They employ a projection method to satisfy physical constraints, resulting in a second-order Crank-Nicolson scheme that is both linear and efficient.
Abstract
The paper introduces a novel method for approximating the Poisson-Nernst-Planck equations with a focus on preserving positivity and conserving mass. By employing a projection step based on L2 minimization, the authors develop an efficient second-order Crank-Nicolson finite difference scheme. The study includes rigorous error estimates in L2 norm, demonstrating the effectiveness of the proposed method through numerical examples.
Stats
Based on initialization: ⟨p0h, 1⟩ = ⟨n0h, 1⟩ = max{⟨Ihp0, 1⟩, ⟨Ihn0, 1⟩} = M0.
Error bounds established under Assumption (A).
Estimates derived for local truncation errors Rkp, Rkn, Rkϕ.