Fast solvers for heterogeneous saddle point problems arising from phasefield-models