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

A. Fortin and D. Yakoubi

Applied Mathematics and Computation, 2019, vol. 358, issue C, 330-347

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.

Keywords: Discontinuous Galerkin; High order method; Stiff ODEs; Adaptive time stepping (search for similar items in EconPapers)
Date: 2019
References: View references in EconPapers View complete reference list from CitEc
Citations: View citations in EconPapers (2)

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0096300319302875
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:apmaco:v:358:y:2019:i:c:p:330-347

DOI: 10.1016/j.amc.2019.04.011

Access Statistics for this article

Applied Mathematics and Computation is currently edited by Theodore Simos

More articles in Applied Mathematics and Computation from Elsevier
Bibliographic data for series maintained by Catherine Liu ().