On the Long Time Simulation of Reaction-Diffusion Equations with Delay

For a consistent numerical method to be practically useful, it is widely accepted that it must preserve the asymptotic stability of the original continuous problem. However, in this study, we show that it may lead to unreliable numerical solutions in long time simulation even if a classical numerical method gives a larger stability region than that of the original continuous problem. Some numerical experiments on the reaction-diffusion equations with delay are presented to confirm our findings. Finally, some open problems on the subject are proposed.