Positivity Preserving and Mass Conservative Projection Method for Poisson-Nernst-Planck Equation

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.