A method for the computation of open-boundary non-linear transient magnetic fields in axisymmetry