Partager cette page :

Structure preserving schemes for Maxwell and PIC: conservation of charge and long-time stability

le 21 janvier 2015

14h00

ENS Rennes Bâtiment Sauvy, Salle 5 (rdc)
Plan d'accès

Séminaire de Martin Campos Pinto (UPMC) au groupe de travail "Applications des mathématiques"

Groupe de travail

/medias/photo/gt-applimaths_1418650100525-png

Lien vers la page Web de l'orateur

Résumé :
In this talk we propose a unified analysis for conforming and non-conforming finite element methods that provides a partial answer to the problem of preserving discrete divergence constraints when computing numerical solutions to the time-dependent Maxwell system. In particular, we identify a set of conditions that a Maxwell solver must satisfy to preserve a compatible discrete Gauss law and in the case of a pure Maxwell problem we give a characterization of a compatible source approximation operator that takes the form of a generalized commuting diagram leading to long-time error and stability estimates.

We next apply these findings by specifying charge conserving schemes for several classes of Galerkin methods such as (i) the usual curl-conforming finite elements, (ii) the centered discontinuous Galerkin (DG) scheme and (iii) a new conforming/non-conforming Galerkin method that shares several advantages with both conforming and DG schemes.

These results are also extended in the case where the Maxwell solver is coupled with a PIC scheme, and numerical results in 2d demonstrate some of the advantages of the proposed methods.
Thématique(s)
Recherche - Valorisation
Contact
Thibaut Deheuvels et Nicolas Crouseilles

Mise à jour le 7 mars 2016