Evaluating the Evans function: Order reduction in numerical methods

Malham, Simon; Niesen, Jitse

doi:10.1090/S0025-5718-07-02016-9

Evaluating the Evans function: Order reduction in numerical methods
HTML articles powered by AMS MathViewer

by Simon Malham and Jitse Niesen PDF

Math. Comp. 77 (2008), 159-179 Request permission

Abstract:

We consider the numerical evaluation of the Evans function, a Wronskian-like determinant that arises in the study of the stability of travelling waves. Constructing the Evans function involves matching the solutions of a linear ordinary differential equation depending on the spectral parameter. The problem becomes stiff as the spectral parameter grows. Consequently, the Gauss–Legendre method has previously been used for such problems; however more recently, methods based on the Magnus expansion have been proposed. Here we extensively examine the stiff regime for a general scalar Schrödinger operator. We show that although the fourth-order Magnus method suffers from order reduction, a fortunate cancellation when computing the Evans matching function means that fourth-order convergence in the end result is preserved. The Gauss–Legendre method does not suffer from order reduction, but it does not experience the cancellation either, and thus it has the same order of convergence in the end result. Finally we discuss the relative merits of both methods as spectral tools.

References

Andrei L. Afendikov and Thomas J. Bridges, Instability of the Hocking-Stewartson pulse and its implications for three-dimensional Poiseuille flow, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 457 (2001), no. 2006, 257–272. MR 1848089, DOI 10.1098/rspa.2000.0665
J. Alexander, R. Gardner, and C. Jones, A topological invariant arising in the stability analysis of travelling waves, J. Reine Angew. Math. 410 (1990), 167–212. MR 1068805
Nairo D. Aparicio, Simon J. A. Malham, and Marcel Oliver, Numerical evaluation of the Evans function by Magnus integration, BIT 45 (2005), no. 2, 219–258. MR 2176193, DOI 10.1007/s10543-005-0001-8
Thomas J. Bridges, Gianne Derks, and Georg Gottwald, Stability and instability of solitary waves of the fifth-order KdV equation: a numerical framework, Phys. D 172 (2002), no. 1-4, 190–216. MR 1946769, DOI 10.1016/S0167-2789(02)00655-3
L. Q. Brin, Numerical testing of the stability of viscous shock waves, Ph.D. thesis, Indiana University, 1998.
Leon Q. Brin, Numerical testing of the stability of viscous shock waves, Math. Comp. 70 (2001), no. 235, 1071–1088. MR 1710652, DOI 10.1090/S0025-5718-00-01237-0
Philip J. Davis and Philip Rabinowitz, Numerical integration, Blaisdell Publishing Co. [Ginn and Co.], Waltham, Mass.-Toronto, Ont.-London, 1967. MR 0211604
Ilan Degani and Jeremy Schiff, RCMS: right correction Magnus series approach for oscillator ODEs, J. Comput. Appl. Math. 193 (2006), no. 2, 413–436. MR 2229552, DOI 10.1016/j.cam.2005.07.001
Kenth Engø, On the construction of geometric integrators in the RKMK class, BIT 40 (2000), no. 1, 41–61. MR 1759033, DOI 10.1023/A:1022362117414
John W. Evans, Nerve axon equations. IV. The stable and the unstable impulse, Indiana Univ. Math. J. 24 (1974/75), no. 12, 1169–1190. MR 393891, DOI 10.1512/iumj.1975.24.24096
M. Facão and D. F. Parker, Stability of screening solitons in photorefractive media, Phys. Rev. E 68 (2003), no. 1, 016610.
R. A. Fisher, The wave of advance of advantageous genes, Ann. Eugenics 7 (1937), 353–369.
C. González, A. Ostermann, and M. Thalhammer, A second-order Magnus-type integrator for nonautonomous parabolic problems, J. Comput. Appl. Math. 189 (2006), no. 1-2, 142–156. MR 2202970, DOI 10.1016/j.cam.2005.04.036
Leon Greenberg and Marco Marletta, Numerical solution of non-self-adjoint Sturm-Liouville problems and related systems, SIAM J. Numer. Anal. 38 (2001), no. 6, 1800–1845. MR 1856233, DOI 10.1137/S0036142999358743
V. V. Gubernov, G. N. Mercer, H. S. Sidhu, and R. O. Weber, Evans function stability of non-adiabatic combustion waves, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 460 (2004), no. 2048, 2415–2435. MR 2073926, DOI 10.1098/rspa.2004.1285
E. Hairer and G. Wanner, Solving ordinary differential equations. II, 2nd ed., Springer Series in Computational Mathematics, vol. 14, Springer-Verlag, Berlin, 1996. Stiff and differential-algebraic problems. MR 1439506, DOI 10.1007/978-3-642-05221-7
Marlis Hochbruck and Christian Lubich, On Magnus integrators for time-dependent Schrödinger equations, SIAM J. Numer. Anal. 41 (2003), no. 3, 945–963. MR 2005189, DOI 10.1137/S0036142902403875
Arieh Iserles, A first course in the numerical analysis of differential equations, Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, 1996. MR 1384977
Arieh Iserles, On the global error of discretization methods for highly-oscillatory ordinary differential equations, BIT 42 (2002), no. 3, 561–599. MR 1931887, DOI 10.1023/A:1022049814688
Arieh Iserles, Hans Z. Munthe-Kaas, Syvert P. Nørsett, and Antonella Zanna, Lie-group methods, Acta numerica, 2000, Acta Numer., vol. 9, Cambridge Univ. Press, Cambridge, 2000, pp. 215–365. MR 1883629, DOI 10.1017/S0962492900002154
Lucas Jódar and Marco Marletta, Solving ODEs arising from non-selfadjoint Hamiltonian eigenproblems, Adv. Comput. Math. 13 (2000), no. 3, 231–256. MR 1772541, DOI 10.1023/A:1018954227133
Todd Kapitula and Björn Sandstede, Edge bifurcations for near integrable systems via Evans function techniques, SIAM J. Math. Anal. 33 (2002), no. 5, 1117–1143. MR 1897705, DOI 10.1137/S0036141000372301
A. Kolmogorov, I. Petrovsky, and N. Piskunov, Étude de l’équation de la diffusion avec croissance de la quantité de matiére et son application à un probléme biologique, Bull. Univ. d’État à Moscou, Série Internationale, A 1 (1937), 1–25.
Wilhelm Magnus, On the exponential solution of differential equations for a linear operator, Comm. Pure Appl. Math. 7 (1954), 649–673. MR 67873, DOI 10.1002/cpa.3160070404
P. C. Moan, Efficient approximation of Sturm-Liouville problems using Lie-group methods, Tech. Report 1998/NA11, DAMTP, University of Cambridge, UK, 1998.
P. C. Moan and J. Niesen, Convergence of the Magnus series, Submitted, 2006. arXiv:math/0609198
Cleve Moler and Charles Van Loan, Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM Rev. 45 (2003), no. 1, 3–49. MR 1981253, DOI 10.1137/S00361445024180
J. D. Murray, Mathematical biology, Biomathematics, vol. 19, Springer-Verlag, Berlin, 1989. MR 1007836, DOI 10.1007/978-3-662-08539-4
J. Niesen, On the global error committed when evaluating the Evans function numerically, Tech. Report HWM 06/43, Dept of Mathematics, Heriot–Watt University, 2006. arXiv:math/0611855
Robert L. Pego and Michael I. Weinstein, Eigenvalues, and instabilities of solitary waves, Philos. Trans. Roy. Soc. London Ser. A 340 (1992), no. 1656, 47–94. MR 1177566, DOI 10.1098/rsta.1992.0055
A. Prothero and A. Robinson, On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations, Math. Comp. 28 (1974), 145–162. MR 331793, DOI 10.1090/S0025-5718-1974-0331793-2
John D. Pryce, Numerical solution of Sturm-Liouville problems, Monographs on Numerical Analysis, The Clarendon Press, Oxford University Press, New York, 1993. Oxford Science Publications. MR 1283388
Björn Sandstede, Stability of travelling waves, Handbook of dynamical systems, Vol. 2, North-Holland, Amsterdam, 2002, pp. 983–1055. MR 1901069, DOI 10.1016/S1874-575X(02)80039-X
Jonathan Swinton and John Elgin, Stability of travelling pulse solutions to a laser equation, Phys. Lett. A 145 (1990), no. 8-9, 428–433. MR 1052867, DOI 10.1016/0375-9601(90)90307-A
David Terman, Stability of planar wave solutions to a combustion model, SIAM J. Math. Anal. 21 (1990), no. 5, 1139–1171. MR 1062397, DOI 10.1137/0521063