A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations. Commun. Comput. Phys. 2018;19(1):53-93. doi:10.4208/cicp.180615.280815a