Discontinuous Galerkin methods for plasma physics in the scrape-off layer of tokamaks

JOURNAL OF COMPUTATIONAL PHYSICS

卷:274

Article

Michoski, C.¹  Meyerson, D.²  Isaac, T.³  Waelbroeck, F.² 

[1] Univ Texas Austin, Dept Aerosp Engn & Engn Mech, Inst Computat Engn & Sci, Austin, TX 78712 USA

[2] Univ Texas Austin, Dept Phys, Inst Fus Studies, Austin, TX 78712 USA

[3] Univ Texas Austin, Inst Computat Engn & Sci, Austin, TX 78712 USA

关键词: Nodal/modal;    Discontinuous Galerkin;    Unified;    Elliptic;    Hyperbolic;    Parabolic;    Plasma physics;    Fusion;    Tokamaks;    Verification;    Blobs;    Modons;    Scrape-off layer;    Interior penalty;    Turbulent saturation;    Upwinding;    Artificial diffusion;    Massively parallel;   

DOI  :  10.1016/j.jcp.2014.06.058

来源: Elsevier

【 摘 要 】

A new parallel discontinuous Galerkin solver, called ArcOn, is developed to describe the intermittent turbulent transport of filamentary blobs in the scrape-off layer (SOL) of fusion plasma. The model is comprised of an elliptic subsystem coupled to two convection-dominated reaction-diffusion-convection equations. Upwinding is used for a class of numerical fluxes developed to accommodate cross product driven convection, and the elliptic solver uses SIPG, NIPG, IIPG, Brezzi, and Bassi-Rebay fluxes to formulate the stiffness matrix. A novel entropy sensor is developed for this system, designed for a space-time varying artificial diffusion/viscosity regularization algorithm. Some numerical experiments are performed to show convergence order on manufactured solutions, regularization of blob/streamer dynamics in the SOL given unstable parameterizations, long-time stability of modon (or dipole drift vortex) solutions arising in simulations of drift-wave turbulence, and finally the formation of edge mode turbulence in the scrape-off layer under turbulent saturation conditions. Published by Elsevier Inc.

【 授权许可】

Free