High order numerical algorithms based on biquadratic spline collocation for two-dimensional parabolic partial differential equations

We report a new algorithm for solving linear parabolic partial differential equations in two space dimension. The algorithm employs optimal biquadratic spline collocation for space discretization and modified trapezoidal rule for time discretization. We need to solve a block tridiagonal linear system at each time step, and obtain an approximate solution with error (Formula presented.) at space-time grid points. We analyse the stability of the new algorithm, and present a stability enhanced variant. Moreover, we give an acceleration strategy based on spectral deferred correction, and the theoretical accuracy can be increased to (Formula presented.), where k is the number of correction loops. We also analyse the stability for the accelerated algorithms. Numerical experiments are attached to demonstrate the effectiveness of the new algorithms.