Numerical solution of hamiltonian systems in reaction-diffusion by symplectic difference schemes

https://doi.org/10.1016/0021-9991(91)90280-XGet rights and content

Abstract

Discrete models in time and space of Fishers equation, ∂u∂t = 2u∂x2 + f(u), in reaction diffusion are numerous in mathematical biology (Weinberger, SIAM J. Math. Anal. 13, 353 (1982) and the references therein). For f(u) =u(1−u) and no dissipation, May (Nature 261, 459 (1976)), using the Euler discretization of the time derivative, found stable solutions (period 2 in time) provided the time step satisfies 2 < k ⩽ √6, the linearized stability for period 1 solutions being 0 < k ⩽ 2. When the dissipation term in discretised form is added to May's ordinary difference scheme, it is shown by 4., 6., and Sleeman (Proc. Roy. Soc. London Ser. A 425, 17 (1989)) that the stable period 2 in time solutions persist. Here it is shown (Sleeman, op. cit.), that when the dissipation term in continuous form is added to May's difference equation, solutions period 2 in time for each value of x satisfy a Hamiltonian system in space. The latter, being non integrable, is solved numerically by Symplectic difference schemes constructed to maintain the values of the Hamiltonain energy up to large values of the space variable (1., 2.). The shape of the solution, in calculations involving 200,000 space steps, is shown to depend crucially on the type and location of the fixed points of the Hamiltonian system in phase space at the position of the initial data at x=0 relative to these fixed points.

References (12)

  • K. Fend

    J. Comput. Math.

    (1986)
  • K. Fling et al.
  • P.C. Fife et al.

    Arch. Rat. Mech. Anal.

    (1977)
  • D.F. Griffiths et al.

    IMA J. Numer. Anal.

    (1988)
  • R.M. May

    Nature

    (1976)
  • A.R. Mitchell et al.

    Numerical Analysis

There are more references available in the full text version of this article.
View full text