Error estimation for numerical approximations of ODEs via composition techniques. Part I: One-step methods

Ahmad Deeb,Denys Dutykh
2024-09-02
Abstract:In this study, we introduce a refined method for ascertaining error estimations in numerical simulations of dynamical systems via an innovative application of composition techniques. Our approach involves a dual application of a basic one-step numerical method of order p in this part, and for the class of Backward Difference Formulas schemes in the second part [Deeb A., Dutykh D. and AL Zohbi M. Error estimation for numerical approximations of ODEs via composition techniques. Part II: BDF methods, Submitted, 2024]. This dual application uses complex coefficients, resulting outputs in the complex plane. The methods innovation lies in the demonstration that the real parts of these outputs correspond to approximations of the solutions with an enhanced order of p + 1, while the imaginary parts serve as error estimations of the same order, a novel proof presented herein using Taylor expansion and perturbation technique. The linear stability of the resulted scheme is enhanced compared to the basic one. The performance of the composition in computing the approximation is also compared. Results show that the proposed technique provide higher accuracy with less computational time. This dual composition technique has been rigorously applied to a variety of dynamical problems, showcasing its efficacy in adapting the time step,particularly in situations where numerical schemes do not have theoretical error estimation. Consequently, the technique holds potential for advancing adaptive time-stepping strategies in numerical simulations, an area where accurate local error estimation is crucial yet often challenging to obtain.
General Mathematics
What problem does this paper attempt to address?