A numerical study of Marangoni flows in zone-refined silicon crystals