Study of spectral stability of generalized Runge-Kutta methods in the initial-boundary problems for the wave equation
DOI:
https://doi.org/10.7242/1999-6691/2015.8.2.10Keywords:
generalized Runge-Kutta methods, Newmark method, hyperbolic equation of the second order, spectral stability, initial-boundary problem, stability functionAbstract
A general algorithm is developed to study the spectral stability of generalized multi-stage Runge-Kutta methods (RKMs) of different orders of accuracy for time integration of the wave equation. The stability function is obtained to estimate the spectral stability of these methods. Spectral stability of various explicit and implicit generalized RKMs is investigated. The behavior of the introduced stability function in a certain generalized RKM is found to be the same as in the previously studied case for the transport equation. It is shown that all the explicit generalized RKMs are spectral unstable, and the implicit generalized RKMs are spectral stable. Moreover, implicit methods based on the formulas of Rado, Lobatto IIIC, Merseta and Burridge possess the false attenuation property (asymptotic stability), and the methods of Gauss-Legendre, Lobatto IIIA, Lobatto IIIB of all orders of accuracy do not possess this property. Using the proposed stability functions, the spectral stability of Newmark’s family of methods is investigated. Computations demonstrate that one of the Newmark methods is a special case of the generalized RKM, namely, the one-step Gauss-Legendre (midpoint) method. Other Newmark methods are spectral unstable or they have the false attenuation property. A comparison is carried out between the approximate solutions obtained in terms of different generalized RKMs and Newmark methods and the exact solution of the problem of free vibrations of a string that is in equilibrium before its movement evoked by a concentrated force, immediately removed at the initial instant of time. The best finding in respect of a “simplicity of realization -achieved accuracy” ratio is the numerical result obtained from the three-stage fourth order diagonally implicit Burridge method, because the complexity of its realization is approximately the same as the Newmark methods, and its accuracy is higher by two orders of magnitude. It has been found that the algorithm developed to study the spectral stability of generalized RKMs and all theoretical results can be transferred with no changes to the parabolic equations, which contain the second time derivatives of unknown functions and describe the dynamic behavior of flexible beams or plates.
Downloads
References
Samarskij A.A., Popov U.P. Raznostnye shemy gazovoj dinamiki. - M.: Editorial URSS, 2009. - 424 s.
2. Igumnov L.A., Ratausko A.U. Sagovyj metod obrasenia preobrazovania Laplasa na uzlah shemy Runge-Kutty // Problemy procnosti i plasticnosti. - 2013. - No 75-3. - S. 178-184.
3. Tihonov A.N., Samarskij A.A. Uravnenia matematiceskoj fiziki. - M.: Nauka, 1977. - 735 s.
4. Novikov E.A. Avnye metody dla zestkih sistem. - Novosibirsk: Nauka. Sibirskoe predpriatie RAN, 1997. - 195 s.
5. Firsov D.K. Ustojcivost’ avnyh shem resenia uravnenij Maksvella metodom kontrol’nyh ob"emov vysokogo poradka tocnosti // Vycislitel’nye metody i programmirovanie. - 2014. - T. 15, No 2. - S. 286-303.
6. Vedernikova E.U., Kornev A.A. Struktura ustojcivogo mnogoobrazia polnost’u neavnyh shem // Vycislitel’nye metody i programmirovanie. - 2013. - T. 14, No 1. - S. 44-49.
7. Bazenov V.G., Igoniceva E.V. Nelinejnye processy udarnogo vypucivania uprugih elementov konstrukcij v vide ortotropnyh obolocek vrasenia. - Niznij Novgorod: Izd-vo Nizegorodskogo universiteta, 1991. - 132 s.
8. Bazenov V.G., Pavlenkova E.V., Artem’eva A.A. Cislennoe resenie obobsennyh osesimmetricnyh zadac dinamiki uprugoplasticeskih obolocek vrasenia pri bol’sih deformaciah // Vycisl. meh. splos. sred. - 2012. - T. 5, No 4. - S. 427-434. DOI
9. De Vogelaere R. A method for the numerical integration of differential equations of second order without explicit first derivatives // J. Res. Nat. Bur. Stand. - 1955. - Vol. 54, no. 3. - P. 119-125. DOI
10. Berezin I.S., Zidkov N.P. Metody vycislenij. - M.: Fizmatgiz, 1959. - T. 2. - 620 s.
11. Levin V.A., Nadkrinicnyj L.V. Cislennoe issledovanie generacii voln na poverhnosti pri pogruzenii tverdogo tela v zidkost’ // Vycisl. meh. splos. sred. - 2011. - T. 4, No 1. - S. 65-73. DOI
12. Zienkiewicz O.C., Taylor R.L. The finite element method. - Oxford: Butterworth-Heinemann, 2000. - 707 p.
13. Kuz’min M.A., Lebedev D.L., Popov B.G. Procnost’, zestkost’, ustojcivost’ elementov konstrukcij. Teoria i praktikum. Rascety na procnost’ elementov mnogoslojnyh kompozitnyh konstrukcij: Uceb. posobie. - M.: Izd-vo MGTU im. N.E. Baumana, 2012. - 344 s.
14. Vu R.V.G., Uitmer E.A. Ustojcivost’ metodov De Vozela cislennogo integrirovania po vremeni // Raketnaa tehnika i kosmonavtika. - 1973. -No 10. - S. 97-100. DOI
15. Dekker K., Verver A. Ustojcivost’ metodov Runge-Kutty dla zestkih nelinejnyh differencial’nyh uravnenij. - M.: Mir, 1988. - 334 s.
16. Banjai L., Messner M., Schanz M. Runge-Kutta convolution quadrature for the boundary element method // Comput. Method. Appl. M. - 2012. - Vol. 245-246. - P. 90-101. DOI
17. Lipanov A.M., Karskanov S.A. Primenenie shem vysokogo poradka approksimacii pri modelirovanii processov tormozenia sverhzvukovyh tecenij v pramougol’nyh kanalah // Vycisl. meh. splos. sred. - 2013. - T. 6, No 3. - S. 292-299. DOI
18. Roman’kov A.S., Romenskij E.I. Metod Runge-Kutty/WENO dla rasceta uravnenij voln maloj amplitudy v nasysennoj uprugoj poristoj srede // SibZVM. - 2014. - T. 17, No 3. - S. 259-271.
19. Okuongae R.I., Ihile M.N.O. L(a)-ustojcivye neavnye metody Runge-Kutty peremennogo poradka so vtoroj proizvodnoj // SibZVM. - 2014. - T. 17, No 4. - S. 373-387.
20. Nemirovskij U.V., Ankovskij A.P. Obobsenie metodov Runge-Kutty i ih primenenie k integrirovaniu nacal’no-kraevyh zadac matematiceskoj fiziki // SibZVM. - 2005. - T. 8, No 1. - S. 57-76.
21. Ankovskij A.P. Issledovanie spektral’noj ustojcivosti obobsennyh metodov Runge-Kutty primenitel’no k cislennomu integrirovaniu nacal’noj zadaci dla uravnenia perenosa // Vycisl. meh. splos. sred. - 2014. - T. 7, No 3. - S. 279-294. DOI
22. Rabotnov U.N. Mehanika deformiruemogo tverdogo tela. - M.: Nauka, 1979. - 744 s.
23. Samarskij A.A. Teoria raznostnyh shem. - M.: Nauka, 1989. - 616 s.
24. Donnell L.G. Balki, plastiny i obolocki. - M.: Nauka, 1982. - 568 s.
25. Vlasov V.Z., Leont’ev N.N. Balki, plity i obolocki na uprugom osnovanii. - M.: Fizmatgiz, 1960. - 492 s.
Downloads
Published
Issue
Section
License
Copyright (c) 2015 Computational Continuum Mechanics
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.