Modeling of a magnetohydrodynamics free surface problem arising in Czochralski crystal growth