×

Factorized schemes of second-order accuracy for numerically solving unsteady problems. (English) Zbl 1359.65083

Summary: 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.

MSC:

65J08 Numerical solutions to abstract evolution equations
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs

Software:

mftoolbox; RODAS
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Ascher U. M., Numerical Methods for Evolutionary Differential Equations, Society for Industrial Mathematics, Philadelphia, 2008.; Ascher, U. M., Numerical Methods for Evolutionary Differential Equations (2008) · Zbl 1157.65048
[2] Baker G. A. and Graves-Morris P. R., Padé Approximants, Cambridge University Press, Cambridge, 1996.; Baker, G. A.; Graves-Morris, P. R., Padé Approximants (1996)
[3] Butcher J. C., Numerical Methods for Ordinary Differential Equations, John Wiley & Sons, Hoboken, 2008.; Butcher, J. C., Numerical Methods for Ordinary Differential Equations (2008) · Zbl 1167.65041
[4] Dekker K. and Verwer J. G., Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations, North-Holland, Amsterdam, 1984.; Dekker, K.; Verwer, J. G., Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations (1984) · Zbl 0571.65057
[5] Gear C. W., Numerical Initial Value Problems in Ordinary Differential Equations, Prentice Hall, Englewood Cliffs, 1971.; Gear, C. W., Numerical Initial Value Problems in Ordinary Differential Equations (1971) · Zbl 1145.65316
[6] Gustafsson B., High Order Difference Methods for Time Dependent PDE, Springer, Berlin, 2008.; Gustafsson, B., High Order Difference Methods for Time Dependent PDE (2008) · Zbl 1146.65064
[7] Hairer E. and Wanner G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Springer, Berlin, 2010.; Hairer, E.; Wanner, G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems (2010) · Zbl 1192.65097
[8] Higham N. J., Functions of Matrices: Theory and Computation, SIAM, Philadelphia, 2008.; Higham, N. J., Functions of Matrices: Theory and Computation (2008) · Zbl 1167.15001
[9] Hundsdorfer W. H. and Verwer J. G., Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations, Springer, Berlin, 2003.; Hundsdorfer, W. H.; Verwer, J. G., Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations (2003) · Zbl 1030.65100
[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.; LeVeque, R. J., Finite Difference Methods for Ordinary and Partial Differential Equations: Steady-State and Time-Dependent Problems (2007) · Zbl 1127.65080
[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.; Moler, C.; Van Loan, C., Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM Rev., 45, 1, 3-49 (2003) · Zbl 1030.65029
[12] Rakitskii Y. V., Ustinov S. M. and Chernorutskii I. G., Numerical Methods for Solving Stiff Systems (in Russian), Nauka, Moscow, 1979.; Rakitskii, Y. V.; Ustinov, S. M.; Chernorutskii, I. G., Numerical Methods for Solving Stiff Systems (1979)
[13] Samarskii A. A., The Theory of Difference Schemes, Marcel Dekker, New York, 2001.; Samarskii, A. A., The Theory of Difference Schemes (2001) · Zbl 0971.65076
[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.; Samarskii, A. A.; Gavrilyuk, I. P.; Makarov, V. L., Stability and regularization of three-level difference schemes with unbounded operator coefficients in Banach spaces, SIAM J. Numer. Anal., 38, 2, 708-723 (2001) · Zbl 1002.65101
[15] Samarskii A. A. and Gulin A. V., Stability of Difference Schemes (in Russian), Nauka, Moscow, 1973.; Samarskii, A. A.; Gulin, A. V., Stability of Difference Schemes (1973) · Zbl 0303.65079
[16] Samarskii A. A., Matus P. P. and Vabishchevich P. N., Difference Schemes with Operator Factors, Kluwer Academic, Dordrecht, 2002.; Samarskii, A. A.; Matus, P. P.; Vabishchevich, P. N., Difference Schemes with Operator Factors (2002) · Zbl 1018.65103
[17] Samarskii A. A. and Vabishchevich P. N., Computational Heat Transfer, Wiley, Chichester, 1995.; Samarskii, A. A.; Vabishchevich, P. N., Computational Heat Transfer (1995)
[18] Vabishchevich P. N., Factorized SM-stable two-level schemes, Comput. Math. Math. Phys. 50 (2010), no. 11, 1818-1824.; Vabishchevich, P. N., Factorized SM-stable two-level schemes, Comput. Math. Math. Phys., 50, 11, 1818-1824 (2010) · Zbl 1224.65240
[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.; Vabishchevich, P. N., Two-level finite difference scheme of improved accuracy order for time-dependent problems of mathematical physics, Comput. Math. Math. Phys., 50, 1, 112-123 (2010) · Zbl 1224.65184
[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.; Vabishchevich, P. N., SM stability for time-dependent problems, Numerical Methods and Applications (Borovets 2010), 29-40 (2011) · Zbl 1317.65180
[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.; Vabishchevich, P. N., Two level schemes of higher approximation order for time dependent problems with skew-symmetric operators, Comput. Math. Math. Phys., 51, 6, 1050-1060 (2011) · Zbl 1249.65175
[22] Vabishchevich P. N., SM-stability of operator-difference schemes, Comput. Math. Math. Phys. 52 (2012), no. 6, 887-894.; Vabishchevich, P. N., SM-stability of operator-difference schemes, Comput. Math. Math. Phys., 52, 6, 887-894 (2012) · Zbl 1274.65260
[23] Vabishchevich P. N., Explicit schemes for parabolic and hyperbolic equations, Appl. Math. Comput. 250 (2015), 424-431.; Vabishchevich, P. N., Explicit schemes for parabolic and hyperbolic equations, Appl. Math. Comput., 250, 424-431 (2015) · Zbl 1328.65185
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.