Stable coupling of the Yee scheme with a linear current model

Filipe da Silva (1) , Martin Campos Pinto (2) , Bruno Després (2) , Stephane Heuraux (3)

1 Instituto de Plasmas e Fusa ̃o Nuclear, Instituto Superior Te ́cnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
2 LJLL - Laboratoire Jacques-Louis Lions
3 IJL - Institut Jean Lamour

This work analyzes the stability of the Yee scheme for non stationary Maxwell's equations coupled with a linear current model. We show that the usual procedure may yield unstable scheme for physical situations that correspond to strongly magnetized plasmas in X-mode (TE) polarization. We propose to use first order clustered discretization of the vectorial product that gives back a stable coupling. We validate the schemes on some test cases representative of direct numerical simulations of X-mode in a magnetic fusion plasma.

Mots clés: Maxwell's equation; fusion plasma; stable scheme; FDTD; linear current model

Domaines: Analyse numérique [math.NA]; Physique mathématique [math-ph]; Physique mathématique [math-ph]

Filipe da Silva, Martin Campos Pinto, Bruno Després, Stephane Heuraux. Stable coupling of the Yee scheme with a linear current model. Journal of Computational Physics, 2015, 295, pp.24-45. ⟨hal-01016883⟩