A Positivity-Preserving Second-Order BDF Scheme for the Cahn-Hilliard Equation with Variable Interfacial Parameters. (2020). Communications in Computational Physics, 28(3), 967-998. https://doi.org/10.4208/cicp.OA-2019-0037