Stability Analysis and Optimization of Semi-Explicit Predictor–Corrector Methods
Abstract
:1. Introduction
2. Materials and Methods
2.1. Semi-Explicit and Semi-Implicit Adams–Bashforth Integration Formula
2.2. Stability Analysis
2.3. Algorithms for Generating the Minimal Finite-Difference Scheme
- Differential equations for which the corresponding right-hand side functions have the least number of occurrences of the different state variables should be calculated first;
- If two or more state equations contain the same number of different variables in the right-hand side functions, then it is worth choosing the one included in the other differential equations.
- A column is added to the feedback matrix end that contains the sums of ones in each row. This number defines the feedbacks number in each equation;
- The feedback matrix is sorted by the last column in ascending order;
- The row with the minimum number in the sum column is considered:
- (1)
- If this is the smallest number of all calculated amounts, i.e., the line with the minimum feedbacks number is unique, the index of the corresponding variable is added to the output array;
- (2)
- If there are several rows with a minimum feedbacks number, then new sums over all rows are considered for each variable’s feedback matrix. Columns corresponding to each of these variables are removed from the matrix. The variable index is chosen to add into the output array, for which at least one of the new calculated sums is the minimum among the sums calculated for all variables. If there are several such variables, then any of them is chosen, for example, the first in order;
- The column and the row corresponding to the variable whose index was added to the output array are removed from the feedback matrix;
- If there are no variables left in the feedback matrix, the algorithm terminates. Otherwise, execution resumes from step 1.
- From the array specifying the order of calculations in the corrector, a variable is chosen at the i-th index;
- For the chosen variable, using the feedback matrix, the indices of those variables for which there are feedbacks in the right-hand side function are determined;
- In the auxiliary array, in a loop, we go through the values by the indices defined in step 2;
- (1)
- If the value from the auxiliary array is 1, then we do nothing;
- (2)
- If the value from the auxiliary array is 0, then we change it to 1, and in the output array, we write the variable to which the considered index corresponds;
- If the value of the auxiliary array at the index determined in step 1 is 0, then change it to 1;
- Check the auxiliary array for content 0;
- (1)
- If the auxiliary array contains only ones, then the algorithm terminates;
- (2)
- If the auxiliary array contains zeros, then increase i by 1 and return to step 1.
3. Experimental Results
3.1. Three-Body Problem
3.2. Simulation of Coupled Oscillators Networks
4. Conclusions and Discussion
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Anderson, S.W.; Sedatole, K. Designing quality into products: The use of accounting data in new product development. Account. Horiz. 1998, 12, 213. [Google Scholar]
- Faragó, I.; Georgiev, K.; Havasi, Á.; Zlatev, Z. Efficient numerical methods for scientific applications: Introduction. Comput. Math. Appl. 2013, 65, 297–300. [Google Scholar] [CrossRef]
- Faragó, I.; Georgiev, K.; Havasi, Á.; Zlatev, Z. Efficient algorithms for large scale scientific computations: Introduction. Comput. Math. Appl. 2014, 67, 2085–2087. [Google Scholar] [CrossRef] [Green Version]
- Butusov, D.; Karimov, A.; Tutueva, A.; Kaplun, D.; Nepomuceno, E.G. The effects of Padé numerical integration in simulation of conservative chaotic systems. Entropy 2019, 21, 362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Jaffe, A. Ordering the universe: The role of mathematics. SIAM Rev. 1984, 26, 473–500. [Google Scholar] [CrossRef]
- Rahrovani, S.; Abrahamsson, T.; Modin, K. An efficient exponential integrator for large nonlinear stiff systems part 1: Theoretical investigation. Conf. Proc. Soc. Exp. Mech. Ser. 2014, 2, 259–268. [Google Scholar]
- Jackson, K.R. The numerical solution of large systems of stiff IVPs for ODEs. Appl. Numer. Math. 1996, 20, 5–20. [Google Scholar] [CrossRef]
- Petcu, D. Software issues in solving initial value problems for ordinary differential equations. Creat. Math. 2004, 13, 97–110. [Google Scholar]
- Butusov, D.N.; Ostrovskii, V.Y.; Tutueva, A.V. Simulation of dynamical systems based on parallel numerical integration methods. In Proceedings of the 2015 IEEE NW Russia Young Researchers in Electrical and Electronic Engineering Conference (EIConRusNW), St. Petersburg, Russia, 2–4 February 2015; pp. 56–59. [Google Scholar]
- Wolfram, S. An elementary introduction to the Wolfram langauge. In Wolfram Media, Incorporated; Wolfram Media, Inc.: Champaign, IL, USA, 2017. [Google Scholar]
- Butusov, D. Adaptive Stepsize Control for Extrapolation Semi-Implicit Multistep ODE Solvers. Mathematics 2021, 9, 950. [Google Scholar] [CrossRef]
- Faleichik, B.V. Minimal residual multistep methods for large stiff non-autonomous linear problems. J. Comput. Appl. Math. 2019, 389, 112498. [Google Scholar] [CrossRef] [Green Version]
- Lopez, S. A predictor-corrector time integration algorithm for dynamic analysis of nonlinear systems. Nonlinear. Dyn. 2020, 101, 1365–1381. [Google Scholar] [CrossRef]
- Aboubakr, A.; Shabana, A.A. Efficient and robust implementation of the TLISMNI method. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, American Society of Mechanical Engineers, Boston, MA, USA, 2–5 August 2015. [Google Scholar]
- Oghonyon, J.G.; Okunuga, S.A.; Omoregbe, N.A.; Agboola, O.O. A Computational Approach in Estimating the Amount of Pond Pollution and Determining the Long Time Behavioural Representation of Pond Pollution Model. Glob. J. Pure Appl. Math. 2015, 11, 2773–2785. [Google Scholar]
- Hamid, A.M.R.; Ahmad, R.R.; Din, U.K.S.; Ismail, A. Modified predictor-corrector multistep method using arithmetic mean for solving ordinary differential equations. In Proceedings of the AIP Conference Proceedings, American Institute of Physics, Putrajaya, Malaysia, 18–20 December 2012; Volume 1522, pp. 703–707. [Google Scholar]
- Biasa, T.; Majid, Z.A.; Suleiman, M. Predictor-corrector block iteration method for solving ordinary differential equations. Sains Malays. 2011, 40, 659–664. [Google Scholar]
- Tutueva, A.; Karimov, T.; Butusov, D. Semi-Implicit and Semi-Explicit Adams-Bashforth-Moulton Methods. Mathematics 2020, 8, 780. [Google Scholar] [CrossRef]
- Cellier, F.E.; Kofman, E. Continuous System Simulation; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
- Kannan, R.; Wang, Z.J. A study of viscous flux formulations for a p-multigrid spectral volume Navier stokes solver. J. Sci. Comput. 2009, 41, 165–199. [Google Scholar] [CrossRef]
- Kannan, R. A high order spectral volume method for elastohydrodynamic lubrication problems: Formulation and application of an implicit p-multigrid algorithm for line contact problems. Comput. Fluids 2011, 48, 44–53. [Google Scholar] [CrossRef]
- Butusov, D.N.; Ostrovskii, V.Y.; Karimov, A.I.; Andreev, V.S. Semi-explicit composition methods in memcapacitor circuit simulation. Int. J. Embed. Real-Time Commun. Syst. 2019, 10, 37–52. [Google Scholar] [CrossRef]
- Yang, L.; Yang, Q.; Chen, G. Hidden attractors, singularly degenerate heteroclinic orbits, multistability and physical realization of a new 6D hyperchaotic system. Comm. Nonlinear. Sci. Numer. Simul. 2020, 90, 105362. [Google Scholar] [CrossRef]
- Kapfer, E.M.; Stapor, P.; Hasenauer, J. Challenges in the calibration of large-scale ordinary differential equation models. IFAC-Pap. 2019, 52, 58–64. [Google Scholar] [CrossRef]
- Gowers, T.; Barrow-Green, J.; Leader, I. The Three-Body Problem. In The Princeton Companion to Mathematics; Princeton University Press: Princeton, NJ, USA, 2008; pp. 726–728. [Google Scholar]
- Marchal, C. The Three-Body Problem; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
- Hairer, E.; Nørsett, S.P.; Wanner, G. Solving Ordinary Differential Equations I: Nonstiff Problems; Springer: Berlin, Germany, 1993; Volume 8. [Google Scholar]
- Yang, Q.; Yang, L.; Ou, B. Hidden hyperchaotic attractors in a new 5D system based on chaotic system with two stable node-foci. Int. J. Bifurc. Chaos 2019, 29, 1950092. [Google Scholar] [CrossRef]
- Anishchenko, V.S.; Strelkova, G.I. Chimera structures in the ensembles of nonlocally coupled chaotic oscillators. Radiophys. Quantum Electron. 2019, 61, 659–671. [Google Scholar] [CrossRef]
- Talebi, B.; Abdi, A. Nordsieck representation of high order predictor-corrector Obreshkov methods and their implementation. CMDE 2019, 7, 16–27. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Tutueva, A.; Butusov, D. Stability Analysis and Optimization of Semi-Explicit Predictor–Corrector Methods. Mathematics 2021, 9, 2463. https://doi.org/10.3390/math9192463
Tutueva A, Butusov D. Stability Analysis and Optimization of Semi-Explicit Predictor–Corrector Methods. Mathematics. 2021; 9(19):2463. https://doi.org/10.3390/math9192463
Chicago/Turabian StyleTutueva, Aleksandra, and Denis Butusov. 2021. "Stability Analysis and Optimization of Semi-Explicit Predictor–Corrector Methods" Mathematics 9, no. 19: 2463. https://doi.org/10.3390/math9192463
APA StyleTutueva, A., & Butusov, D. (2021). Stability Analysis and Optimization of Semi-Explicit Predictor–Corrector Methods. Mathematics, 9(19), 2463. https://doi.org/10.3390/math9192463