Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

Tonatiuh Sánchez-Vizuet, Manuel E. Solano, Antoine J. Cerfon

Research output: Contribution to journalArticlepeer-review

Abstract

We propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

Original languageEnglish (US)
Article number107239
JournalComputer Physics Communications
Volume255
DOIs
StatePublished - Oct 2020

Keywords

  • Adaptive Hybridizable Discontinuous Galerkin (HDG)
  • Curved boundaries
  • Local mesh refinement
  • Plasma Equilibrium
  • Residual error estimator
  • Un-fitted mesh

ASJC Scopus subject areas

  • Hardware and Architecture
  • Physics and Astronomy(all)

Fingerprint Dive into the research topics of 'Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains'. Together they form a unique fingerprint.

Cite this