The last advection scheme is very accurate in space (fourth-order accuracy) but less accurate in time (first-order accuracy) if the velocity field for the current moment of time (i.e. A-configuration velocity field) is used for B, C and D points. An
alternative algorithm is to use the Runge–Kutta scheme in both space and time: all markers are first displaced to their B points and then material properties are reinterpolated to Eulerian nodes and a new B-configuration velocity field is computed by solving the momentum and continuity equations, this B-velocity field is then used for moving markers to their C points for computing a C-configuration velocity