RUNGE-KUTTA METHODS FOR PARTIAL-DIFFERENTIAL EQUATIONS AND FRACTIONAL ORDERS OF CONVERGENCE

We apply Runge-Kutta methods to linear partial differential equations of the form u(t)(x, t) = L(x, partial derivative) u(x, t) + f(x, t). Under appropriate assumptions on the eigenvalues of the operator L and the (generalized) Fourier coefficients of f , we give a sharp lower bound for the order of convergence of these methods. We further show that this order is, in general, fractional and that it depends on the L(r)-norm used to estimate the global error. The analysis also applies to systems arising from spatial discretization of partial differential equations by finite differences or finite element techniques. Numerical examples illustrate the results.