A robust ‘time-marching’ solver for one-dimensional nucleating steam flows