A locally field-aligned discontinuous Galerkin method for anisotropic wave equations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:408

Article

Dingfelder, Benedict¹  Hindenlang, Florian J.¹ 

[1] Max Planck Inst Plasma Phys, Boltzmannstr 2, D-85748 Garching, Germany

关键词: Discontinuous Galerkin;    Field-alignment;    Non-conforming;    Eigenvalue problem;    Plasma physics;    Stellarator;   

DOI  :  10.1016/j.jcp.2020.109273

来源: Elsevier

【 摘 要 】

In magnetized plasmas of fusion devices the strong magnetic field leads to highly anisotropic physics where solution scales along field lines are much larger than perpendicular to it. Hence, regarding both accuracy and efficiency, a numerical method should allow to address parallel and perpendicular resolutions independently. In this work, we consider the eigenvalue problem of a two-dimensional anisotropic wave equation with variable coefficients which is a simplified model of linearized ideal magnetohydrodynamics. For this, we propose to use a mesh that is aligned with the magnetic field and choose to discretize the problem with a discontinuous Galerkin method which naturally allows for non-conforming interfaces. First, we analyze the eigenvalue spectrum of a constant coefficient anisotropic wave equation, and demonstrate that this approach improves the accuracy by up to seven orders of magnitude, if compared to a non-aligned method with the same number of degrees of freedom. In particular, the results improve for eigenfunctions with high mode numbers. We also apply the method to compute the eigenvalue spectrum of the associated anisotropic wave equation with variable coefficients of flux surfaces of a Stellarator configuration. We benchmark the results against a spectral code. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free