Advection Calculations Using Spline Schemes

The calculation of advective transport in one and two dimensions is done by the method of characteristics using spline interpolating polynomials at the forward time step. The method has attractive accuracy, is unconditionally stable, applies equally well to both uniform and variable velocities and mesh widths, and can eliminate numerical oscillations without severely diffusing concentration gradients.