Towards an unconditionally stable numerical scheme for continuum dislocation transport