A numerical method for curvature driven boundary migration