THE GLOBAL DYNAMICS OF DISCRETE SEMILINEAR PARABOLIC EQUATIONS

A class of scalar semilinear parabolic equations possessing absorbing sets, a Lyapunov functional, and a global attractor are considered. The gradient structure of the problem implies that, provided all steady states are isolated, solutions approach a steady state as t -* oo. The dynamical properties of various finite difference and finite element schemes for the equations are analysed. The existence of absorbing sets, bounded independently of the mesh size, is proved for the numerical methods. Discrete Lyapunov functions are constructed to show that, under appropriate conditions on the mesh parameters, numerical orbits approach steady state solutions as discrete time increases. However, it is shown that insufficient spatial resolution can introduce deceptively smooth spurious steady solutions and cause the stability properties of the true steady solutions to be incorrectly represented. Furthermore, it is also shown that the explicit Euler scheme introduces spurious solutions with period 2 in the timestep. As a result, the absorbing set is destroyed and there is initial data leading to blow up of the scheme, however small the mesh parameters are taken. To obtain stabilization to a steady state for this scheme, it is necessary to restrict the timestep in terms of the initial data and the space step. Implicit schemes are constructed for which absorbing sets and Lyapunov functions exist under restrictions on the timestep that are independent of initial data and of the space step; both one-step and multistep (BDF) methods are studied.

[1]  A. R. Mitchell,et al.  Stable periodic bifurcations of an explicit discretization of a nonlinear partial differential equation in reaction diffusion , 1988 .

[2]  Andrew M. Stuart,et al.  The Dynamics of the Theta Method , 1991, SIAM J. Sci. Comput..

[3]  Jack K. Hale,et al.  Upper semicontinuity of attractors for approximations of semigroups and partial differential equations , 1988 .

[4]  Daniel B. Henry.  Geometric Theory of Semilinear Parabolic Equations , 1989 .

[5]  David Hoff,et al.  Stability and Convergence of Finite Difference Methods for Systems of Nonlinear Reaction-Diffusion Equations , 1978 .

[6]  G. R. Shubin,et al.  Multiple Solutions and Bifurcation of Finite Difference Approximations to Some Steady Problems of Fluid Dynamics , 1981 .

[7]  Andrew M. Stuart,et al.  Runge-Kutta methods for dissipative and gradient dynamical systems , 1994 .

[8]  Charles M. Elliott,et al.  The Cahn-Hilliard Model for the Kinetics of Phase Separation , 1989 .

[9]  Jie Shen,et al.  Convergence of approximate attractors for a fully discrete system for reaction-diffusion equations , 1989 .

[10]  A. R. Humphries.  Spurious solutions of numerical methods for initial value problems , 1993 .

[11]  Andrew M. Stuart.  Linear Instability Implies Spurious Periodic Solutions , 1989 .

[12]  R. Temam.  Infinite Dimensional Dynamical Systems in Mechanics and Physics Springer Verlag , 1993 .

[13]  Edriss S. Titi,et al.  Dissipativity of numerical schemes , 1991 .

[14]  J. Carr,et al.  Metastable patterns in solutions of ut = ϵ2uxx − f(u) , 1989 .

[15]  S. G. N. Khalsa.  Finite element approximation of a reaction-diffusion equation. I. Application of topological techniques to the analysis of asymptotic behavior of the semidiscrete approximations , 1986 .

[16]  J. Rappaz,et al.  On numerical approximation in bifurcation theory , 1990 .

[17]  J. Smoller.  Shock Waves and Reaction-Diffusion Equations , 1983 .

[18]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[19]  Hiroshi Matano,et al.  Asymptotic Behavior and Stability of Solutions of Semilinear Diffusion Equations , 1979 .

[20]  Jack K. Hale,et al.  Lower semicontinuity of attractors of gradient systems and applications , 1989 .

[21]  E. F. Infante,et al.  A Bifurcation Problem for a Nonlinear Partial Differential Equation of Parabolic Type. , 1974 .