摘要
In this paper,ETD3-Padéand ETD4-PadéGalerkin finite element methods are proposed and analyzed for nonlinear delayed convection-diffusion-reaction equations with Dirichlet boundary conditions.An ETD-based RK is used for time integration of the corresponding equation.To overcome a well-known difficulty of numerical instability associated with the computation of the exponential operator,the Padéapproach is used for such an exponential operator approximation,which in turn leads to the corresponding ETD-Padéschemes.An unconditional L^(2) numerical stability is proved for the proposed numerical schemes,under a global Lipshitz continuity assumption.In addition,optimal rate error estimates are provided,which gives the convergence order of O(k^(3)+h^(r))(ETD3-Padé)or O(k^(4)+h^(r))(ETD4-Padé)in the L^(2)norm,respectively.Numerical experiments are presented to demonstrate the robustness of the proposed numerical schemes.
基金
NSFC 11971047(Q.Huang)and NSF DMS-2012669(C.Wang).