IDEAS home Printed from https://ideas.repec.org/a/eee/matcom/v64y2004i5p507-520.html
   My bibliography  Save this article

Stable schemes for partial differential equations: the one-dimensional reaction–diffusion equation

Author

Listed:
  • Teixeira, João

Abstract

Reaction–diffusion equations are fundamental in modelling several natural phenomena. The numerical schemes used to solve these equations often suffer from numerical stability problems. In this paper, a new type of algorithm to solve the diffusion equation in a stable and explicit manner is extended to the reaction–diffusion equation. The new scheme imposes a fixed value for the stability coefficient below the stability limit, and uses this information in order to determine a new grid. The values of the variables at this grid are then obtained by interpolation from the original grid. The scheme is applied to the linear single kinetic reaction–diffusion equation and to the classical Fisher equation. Different possibilities of extending the new scheme to the reaction–diffusion equations are discussed. It is shown that, for the linear case, including both terms (reaction and diffusion) in the computation of the new grid gives more accurate results and is more correct than just including the diffusion term. To solve the non-linear Fisher equation, a fractional-step method, where the reaction and diffusion terms are solved separately, is chosen. The new scheme provides realistic results when compared with analytic solutions.

Suggested Citation

  • Teixeira, João, 2004. "Stable schemes for partial differential equations: the one-dimensional reaction–diffusion equation," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 64(5), pages 507-520.
  • Handle: RePEc:eee:matcom:v:64:y:2004:i:5:p:507-520
    DOI: 10.1016/j.matcom.2003.10.001
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0378475403001599
    Download Restriction: Full text for ScienceDirect subscribers only

    File URL: https://libkey.io/10.1016/j.matcom.2003.10.001?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Witek, Marcin L. & Teixeira, Joao & Flatau, Piotr J., 2008. "On stable and explicit numerical methods for the advection–diffusion equation," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 79(3), pages 561-570.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:matcom:v:64:y:2004:i:5:p:507-520. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    We have no bibliographic references for this item. You can help adding them by using this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: http://www.journals.elsevier.com/mathematics-and-computers-in-simulation/ .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.