A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations. (2018). Communications in Computational Physics, 19(1), 53-93. https://doi.org/10.4208/cicp.180615.280815a