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