Phase-field simulation of superconductor vortex clustering in the vicinity of ferromagnetic domain bifurcations