A Positivity-Preserving Second-Order BDF Scheme for the Cahn-Hilliard Equation with Variable Interfacial Parameters. Commun. Comput. Phys. 2020;28(3):967-998. doi:10.4208/cicp.OA-2019-0037