Abstract
Schemes with the second-order approximation in time are considered for numerically solving the Cauchy problem for an evolutionary equation of first order with a self-adjoint operator. The implicit two-level scheme based on the Padé polynomial approximation is unconditionally stable. It demonstrates good asymptotic properties in time and provides an adequate evolution in time for individual harmonics of the solution (has spectral mimetic (SM) stability). In fact, the only drawback of this scheme is the necessity to solve an equation with an operator polynomial of second degree at each time level. We consider modifications of these schemes, which are based on solving equations with operator polynomials of first degree. Such computational implementations occur, for example, if we apply the fully implicit two-level scheme (the backward Euler scheme). A three-level modification of the SM-stable scheme is proposed. Its unconditional stability is established in the corresponding norms. The emphasis is on the scheme, where the numerical algorithm involves two stages, namely, the backward Euler scheme of first order at the first (prediction) stage and the following correction of the approximate solution using a factorized operator. The SM-stability is established for the proposed scheme. To illustrate the theoretical results of the work, a model problem is solved numerically.
Funding source: Russian Foundation for Basic Research
Award Identifier / Grant number: 14-01-00785
Funding statement: This work was supported by the Russian Foundation for Basic Research (project 14-01-00785).
References
[1] Ascher U. M., Numerical Methods for Evolutionary Differential Equations, Society for Industrial Mathematics, Philadelphia, 2008. 10.1137/1.9780898718911Search in Google Scholar
[2] Baker G. A. and Graves-Morris P. R., Padé Approximants, Cambridge University Press, Cambridge, 1996. 10.1017/CBO9780511530074Search in Google Scholar
[3] Butcher J. C., Numerical Methods for Ordinary Differential Equations, John Wiley & Sons, Hoboken, 2008. 10.1002/9780470753767Search in Google Scholar
[4] Dekker K. and Verwer J. G., Stability of Runge–Kutta Methods for Stiff Nonlinear Differential Equations, North-Holland, Amsterdam, 1984. Search in Google Scholar
[5] Gear C. W., Numerical Initial Value Problems in Ordinary Differential Equations, Prentice Hall, Englewood Cliffs, 1971. Search in Google Scholar
[6] Gustafsson B., High Order Difference Methods for Time Dependent PDE, Springer, Berlin, 2008. Search in Google Scholar
[7] Hairer E. and Wanner G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Springer, Berlin, 2010. Search in Google Scholar
[8] Higham N. J., Functions of Matrices: Theory and Computation, SIAM, Philadelphia, 2008. 10.1137/1.9780898717778Search in Google Scholar
[9] Hundsdorfer W. H. and Verwer J. G., Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations, Springer, Berlin, 2003. 10.1007/978-3-662-09017-6Search in Google Scholar
[10] LeVeque R. J., Finite Difference Methods for Ordinary and Partial Differential Equations: Steady-State and Time-Dependent Problems, Society for Industrial Mathematics, Philadelphia, 2007. 10.1137/1.9780898717839Search in Google Scholar
[11] Moler C. and Van Loan C., Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM Rev. 45 (2003), no. 1, 3–49. 10.1137/S00361445024180Search in Google Scholar
[12] Rakitskii Y. V., Ustinov S. M. and Chernorutskii I. G., Numerical Methods for Solving Stiff Systems (in Russian), Nauka, Moscow, 1979. Search in Google Scholar
[13] Samarskii A. A., The Theory of Difference Schemes, Marcel Dekker, New York, 2001. 10.1201/9780203908518Search in Google Scholar
[14] Samarskii A. A., Gavrilyuk I. P. and Makarov V. L., Stability and regularization of three-level difference schemes with unbounded operator coefficients in Banach spaces, SIAM J. Numer. Anal. 38 (2001), no. 2, 708–723. 10.1137/S0036142999357221Search in Google Scholar
[15] Samarskii A. A. and Gulin A. V., Stability of Difference Schemes (in Russian), Nauka, Moscow, 1973. Search in Google Scholar
[16] Samarskii A. A., Matus P. P. and Vabishchevich P. N., Difference Schemes with Operator Factors, Kluwer Academic, Dordrecht, 2002. 10.1007/978-94-015-9874-3Search in Google Scholar
[17] Samarskii A. A. and Vabishchevich P. N., Computational Heat Transfer, Wiley, Chichester, 1995. Search in Google Scholar
[18] Vabishchevich P. N., Factorized SM-stable two-level schemes, Comput. Math. Math. Phys. 50 (2010), no. 11, 1818–1824. 10.1134/S0965542510110059Search in Google Scholar
[19] Vabishchevich P. N., Two-level finite difference scheme of improved accuracy order for time-dependent problems of mathematical physics, Comput. Math. Math. Phys. 50 (2010), no. 1, 112–123. 10.1134/S0965542510010100Search in Google Scholar
[20] Vabishchevich P. N., SM stability for time-dependent problems, Numerical Methods and Applications (Borovets 2010), Lecture Notes in Comput. Sci. 6046, Springer, Berlin (2011), 29–40. 10.1007/978-3-642-18466-6_3Search in Google Scholar
[21] Vabishchevich P. N., Two level schemes of higher approximation order for time dependent problems with skew-symmetric operators, Comput. Math. Math. Phys. 51 (2011), no. 6, 1050–1060. 10.1134/S0965542511060170Search in Google Scholar
[22] Vabishchevich P. N., SM-stability of operator-difference schemes, Comput. Math. Math. Phys. 52 (2012), no. 6, 887–894. 10.1134/S0965542512060140Search in Google Scholar
[23] Vabishchevich P. N., Explicit schemes for parabolic and hyperbolic equations, Appl. Math. Comput. 250 (2015), 424–431. 10.1016/j.amc.2014.10.124Search in Google Scholar
© 2017 by De Gruyter
