Design of Stable Adaptive Step Size Controller for ODE Numerical Solutions Based on Condition Number Optimization

Jia Dong

doi:10.54097/678fjv81

Authors

Jia Dong School of Tianjin Foreign Studies University, Tianjin, 300204, China

DOI:

https://doi.org/10.54097/678fjv81

Keywords:

Ordinary Differential Equations; Adaptive Step Size; Condition Number Optimization; Numerical Stability; Jacobian Matrix; Control Theory.

Abstract

A numerical method to solve ordinary differential equations (ODEs) is used in computational science for many kinds of complex dynamic systems in all fields of engineering and science. Conventional adaptive step-size controllers are mainly concerned with bounding the local truncation error to maintain solution accuracy. However, for stiff systems, implicit numerical methods need to solve non-linear algebraic equations at each step and usually employ Newton-Raphson iterations. Therefore, at this time, the condition number of the iteration matrix affects both the stability and speed of convergence in practice. A new adaptive step-size controller that considers the condition number optimisation in the control loop is proposed in this paper. Continuously monitor and bound the condition number of the Jacobian-derived iteration matrix to adjust the step size dynamically in the proposed dual-loop controller, thus avoiding divergence of Newton iteration and reducing the amplification of floating-point errors. Based on the above theoretical analysis and control system construction, it has been shown that adding condition number constraints can enhance both the stability and efficiency of numerical integration for stiff ODEs. The above results provide an all-encompassing system for the development of high-accuracy and algebraically stable numerical solvers.

Downloads

Download data is not yet available.

References

[1] Gustafsson, K. (1991). Control theoretic techniques for stepsize selection in explicit Runge Kutta methods. ACM Transactions on Mathematical Software, 17(4), 533 554. https://doi.org/10.1145/210232.210242

[2] Hairer, E., Nørsett, S. P., & Wanner, G. (1993). Solving ordinary differential equations I: Nonstiff problems (2nd ed.). Springer Verlag. https://doi.org/10.1007/978 3 540 78862 1

[3] Hairer, E., & Wanner, G. (1996). Solving ordinary differential equations II: Stiff and differential algebraic problems (2nd ed.). Springer Verlag. https://doi.org/10.1007/978 3 642 05221 7

[4] Söderlind, G. (2002). Automatic control and adaptive time stepping. Numerical Algorithms, 31(1 4), 281 310. https://doi.org/10.1023/A:1021160023092

[5] Higham, N. J. (2002). Accuracy and stability of numerical algorithms (2nd ed.). SIAM. https://doi.org/10.1137/1.9780898718027

[6] Söderlind, G., & Wang, L. (2006). Adaptive time stepping and computational stability. Journal of Computational and Applied Mathematics, 185(2), 225 243. https://doi.org/10.1016/j.cam.2005.03.008

[7] Hindmarsh, A. C. (1983). ODEPACK, a systematized collection of ODE solvers. In R. S. Stepleman et al. (Eds.), Scientific Computing (pp. 55 64). North Holland.

[8] Ascher, U. M., & Petzold, L. R. (1998). Computer methods for ordinary differential equations and differential algebraic equations. SIAM. https://doi.org/10.1137/1.9781611971392

[9] Kennedy, C. A., & Carpenter, M. H. (2003). Additive Runge Kutta schemes for convection diffusion reaction equations. Applied Numerical Mathematics, 44(1 2), 139 181. https://doi.org/10.1016/S0168 9274(02)00138 1

[10] Rackauckas, C., & Nie, Q. (2017). DifferentialEquations.jl: A performant and feature rich ecosystem for solving differential equations in Julia. Journal of Open Research Software, 5(1), Article 15. https://doi.org/10.5334/jors.151