A New Higher-Order Godunov Method for General Hyperbolic Equations

A new numerical method of the higher-order Godunov-type is proposed. The spatial profile within a grid is interpolated with the cubic spline. The time integration is represented by the time evolution of the value and the first spatial derivative of the quantity. This scheme is successfully applied to a linear wave propagation with less diffusion retaining monotonicity and to the nonlinear hydrodynamic equations. The extension to higher dimensions and to more general hyperbolic equations including coupled hyperbolic-parabolic equations is straightforward.