High accuracy cubic spline finite difference approximation for the solution of one-space dimensional non-linear wave equations

In this paper, we propose a new three-level implicit nine point compact cubic spline finite difference formulation of order two in time and four in space directions, based on cubic spline approximation in x-direction and finite difference approximation in t-direction for the numerical solution of one-space dimensional second order non-linear hyperbolic partial differential equations. We describe the mathematical formulation procedure in details and also discuss how our formulation is able to handle wave equation in polar coordinates. The proposed method when applied to a linear hyperbolic equation is also shown to be unconditionally stable. Numerical results are provided to justify the usefulness of the proposed method.