On the use of local maximum entropy approximants for Cahn–Hilliard phase-field models in 2D domains and on surfaces