IDEAS home Printed from https://ideas.repec.org/a/eee/apmaco/v358y2019icp330-347.html
   My bibliography  Save this article

An adaptive discontinuous Galerkin method for very stiff systems of ordinary differential equations

Author

Listed:
  • Fortin, A.
  • Yakoubi, D.

Abstract

We present a discontinuous Galerkin (DG) method for the numerical solution of stiff systems of ordinary differential equations (ODEs). We use a standard DG variational formulation with polynomials of degree k in each time interval. We show that the method is A-stable for every k. We then introduce a hierarchical Legendre finite element basis and we show that a whole family of approximations can be obtained simply by truncating the last p degrees of freedom from the computed solution. We show that these approximations converge to order k+1−p in L2-norm and to order k+1/2−p in supremum norm. We then show how this can be used to control the error and the time step length. We present numerical examples of solutions on very stiff problems and on stiff problems with very long time integration where the time step length can vary on many orders of magnitude.

Suggested Citation

  • Fortin, A. & Yakoubi, D., 2019. "An adaptive discontinuous Galerkin method for very stiff systems of ordinary differential equations," Applied Mathematics and Computation, Elsevier, vol. 358(C), pages 330-347.
  • Handle: RePEc:eee:apmaco:v:358:y:2019:i:c:p:330-347
    DOI: 10.1016/j.amc.2019.04.011
    as

    Download full text from publisher

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

    File URL: https://libkey.io/10.1016/j.amc.2019.04.011?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. Fortin, A. & Yakoubi, D., 2023. "A very high order discontinuous Galerkin method for the numerical solution of stiff DDEs," Applied Mathematics and Computation, Elsevier, vol. 443(C).
    2. Meng, Jian & Mei, Liquan, 2020. "Discontinuous Galerkin methods of the non-selfadjoint Steklov eigenvalue problem in inverse scattering," Applied Mathematics and Computation, Elsevier, vol. 381(C).

    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:apmaco:v:358:y:2019:i:c:p:330-347. 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: https://www.journals.elsevier.com/applied-mathematics-and-computation .

    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.