Efficient exponential Rosenbrock methods till order four

Abstract

In a previous paper, a technique was described to avoid order reduction with exponential Rosenbrock methods when integrating initial boundary value problems with time-dependent boundary conditions. That requires to calculate some information on the boundary from the given data. In the present paper we prove that, under some assumptions on the coefficients of the method which are mainly always satisfied, no numerical differentiation is required to approximate that information in order to achieve order 4 for parabolic problems with Dirichlet boundary conditions. With Robin/Neumann ones, just numerical differentiation in time may be necessary for order 4, but none for order $\le 3$.

Furthermore, as with this technique it is not necessary to impose any stiff order conditions, in search of efficiency, we recommend some methods of classical orders 2, 3 and 4 and we give some comparisons with several methods in the literature, with the corresponding stiff order.