Stabilized Galerkin for transient advection of differential forms

Holger Heumann (1, 2) , Ralf Hiptmair (3) , Cecilia Pagliantini (3)

1 CASTOR - Control, Analysis and Simulations for TOkamak Research
2 LJAD - Laboratoire Jean Alexandre Dieudonné
3 SAM - Seminar for Applied Mathematics [Zürich]

We deal with the discretization of generalized transient advection problems for differential forms on bounded spatial domains. We pursue an Eulerian method of lines approach with explicit timestepping. Concerning spatial discretization we extend the jump stabilized Galerkin discretization proposed in [ H. HEUMANN and R.HIPTMAIR, Stabilized Galerkin methods for magnetic advection, Math. Modelling Numer. Analysis, 47 (2013), pp.1713--1732] to forms of any degree and, in particular, advection velocities that may have discontinuities resolved by the mesh. A rigorous a priori convergence theory is established for Lipschitz continuous velocities, conforming meshes and standard finite element spaces of discrete differential forms. However, numerical experiments furnish evidence of the good performance of the new method also in the presence of jumps of the advection velocity.

Domaines: Mathématiques [math]

Holger Heumann, Ralf Hiptmair, Cecilia Pagliantini. Stabilized Galerkin for transient advection of differential forms. Discrete and Continuous Dynamical Systems - Series S, 2016, 9 (1), pp.185 - 214. ⟨10.3934/dcdss.2016.9.185⟩. ⟨hal-01248140⟩