Fully discrete semi-Lagrangian methods for advection of differential forms

Holger Heumann, Ralf Hiptmair, Kun Li, Jinchao Xu

Research output: Contribution to journalArticlepeer-review

3 Scopus citations

Abstract

We study the discretization of linear transient transport problems for differential forms on bounded domains. The focus is on unconditionally stable semi-Lagrangian methods that employ finite element approximation on fixed meshes combined with tracking of the flow map. We derive these methods as finite element Galerkin approach to discrete material derivatives and discuss further approximations leading to fully discrete schemes. We establish comprehensive a priori error estimates, in particular a new asymptotic estimate of order O(hr+1 τ -1/2) for the L2-error of semi-Lagrangian schemes with exact L2-projection. Here, h is the spatial meshwidth, τ denotes the timestep, and r is the (full) polynomial degree of the piecewise polynomial discrete differential forms used as trial functions. Yet, numerical experiments hint that the estimates may still be sub-optimal for spatial discretization with lowest order discrete differential forms.

Original languageEnglish (US)
Pages (from-to)981-1007
Number of pages27
JournalBIT Numerical Mathematics
Volume52
Issue number4
DOIs
StatePublished - Dec 2012

All Science Journal Classification (ASJC) codes

  • Software
  • Computer Networks and Communications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Fully discrete semi-Lagrangian methods for advection of differential forms'. Together they form a unique fingerprint.

Cite this