DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation

Abstract

In axisymmetric fusion reactors, the equilibrium magnetic configuration can be expressed in terms of the solution to a semi-linear elliptic equation known as the Grad–Shafranov equation, the solution of which determines the poloidal component of the magnetic field. When the geometry of the confinement region is known, the problem becomes an interior Dirichlet boundary value problem. We propose a high order solver based on the Hybridizable Discontinuous Galerkin method. The resulting algorithm (1) provides high order of convergence for the flux function and its gradient, (2) incorporates a novel method for handling piecewise smooth geometries by extension from polygonal meshes, (3) can handle geometries with non-smooth boundaries and x-points, and (4) deals with the semi-linearity through an accelerated two-grid fixed-point iteration. The effectiveness of the algorithm is verified with computations for cases where analytic solutions are known on configurations similar to those of actual devices (ITER with single null and double null divertor, NSTX, ASDEX upgrade, and Field Reversed Configurations).

Authors:
ORCiD logo [1];  [2]
  1. New York Univ. (NYU), NY (United States)
  2. Univ. de Concepción (Chile)
Publication Date:
Research Org.:
New York Univ. (NYU), NY (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1610077
Alternate Identifier(s):
OSTI ID: 1635825
Grant/Contract Number:  
FG02-86ER53233
Resource Type:
Accepted Manuscript
Journal Name:
Computer Physics Communications
Additional Journal Information:
Journal Volume: 235; Journal ID: ISSN 0010-4655
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Hybridizable Discontinuous Galerkin (HDG); curved boundary; Grad–Shafranov; Anderson acceleration; plasma equilibrium; magnetohydrodynamics (MHD)

Citation Formats

Sánchez-Vizuet, Tonatiuh, and Solano, Manuel E. A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation. United States: N. p., 2018. Web. doi:10.1016/j.cpc.2018.09.013.
Sánchez-Vizuet, Tonatiuh, & Solano, Manuel E. A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation. United States. https://doi.org/10.1016/j.cpc.2018.09.013
Sánchez-Vizuet, Tonatiuh, and Solano, Manuel E. Fri . "A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation". United States. https://doi.org/10.1016/j.cpc.2018.09.013. https://www.osti.gov/servlets/purl/1610077.
@article{osti_1610077,
title = {A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation},
author = {Sánchez-Vizuet, Tonatiuh and Solano, Manuel E.},
abstractNote = {In axisymmetric fusion reactors, the equilibrium magnetic configuration can be expressed in terms of the solution to a semi-linear elliptic equation known as the Grad–Shafranov equation, the solution of which determines the poloidal component of the magnetic field. When the geometry of the confinement region is known, the problem becomes an interior Dirichlet boundary value problem. We propose a high order solver based on the Hybridizable Discontinuous Galerkin method. The resulting algorithm (1) provides high order of convergence for the flux function and its gradient, (2) incorporates a novel method for handling piecewise smooth geometries by extension from polygonal meshes, (3) can handle geometries with non-smooth boundaries and x-points, and (4) deals with the semi-linearity through an accelerated two-grid fixed-point iteration. The effectiveness of the algorithm is verified with computations for cases where analytic solutions are known on configurations similar to those of actual devices (ITER with single null and double null divertor, NSTX, ASDEX upgrade, and Field Reversed Configurations).},
doi = {10.1016/j.cpc.2018.09.013},
journal = {Computer Physics Communications},
number = ,
volume = 235,
place = {United States},
year = {Fri Sep 28 00:00:00 EDT 2018},
month = {Fri Sep 28 00:00:00 EDT 2018}
}

Journal Article:

Citation Metrics:
Cited by: 10 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Computation of MHD equilibrium of tokamak plasma
journal, March 1991


Axialsymmetrische magnetohydrodynamische Gleichgewichtskonfigurationen
journal, October 1957

  • Lüst, R.; Schlüter, A.
  • Zeitschrift für Naturforschung A, Vol. 12, Issue 10
  • DOI: 10.1515/zna-1957-1014

The ITER divertor concept
journal, April 1995


Engineering and design aspects related to the development of the ITER divertor
journal, March 1995


Isoparametric Bicubic Hermite Elements for Solution of the Grad-Shafranov Equation
journal, March 1991

  • Huysmans, G. T. A.; Goedbloed, J. P.; Kerner, W.
  • International Journal of Modern Physics C, Vol. 02, Issue 01
  • DOI: 10.1142/S0129183191000512

Solving the Grad–Shafranov equation with spectral elements
journal, May 2014


A mimetic spectral element solver for the Grad–Shafranov equation
journal, July 2016


Quasi-static free-boundary equilibrium of toroidal plasma with CEDRES++: Computational methods and applications
journal, January 2015


A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries
journal, April 2017


Application of Parallelized Multigrid Method to Solution of MHD Equilibrium With MPI
journal, March 2004


High-resolution magnetohydrodynamic equilibrium code for unity beta plasmas
journal, July 2006

  • Gourdain, P. -A.; Leboeuf, J. -N.; Neches, R. Y.
  • Journal of Computational Physics, Vol. 216, Issue 1
  • DOI: 10.1016/j.jcp.2005.12.005

A fast, high-order solver for the Grad–Shafranov equation
journal, June 2013

  • Pataki, Andras; Cerfon, Antoine J.; Freidberg, Jeffrey P.
  • Journal of Computational Physics, Vol. 243
  • DOI: 10.1016/j.jcp.2013.02.045

ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria
journal, May 2015


Meshless method for solving fixed boundary problem of plasma equilibrium
journal, July 2015


Numerical solution of ideal MHD equilibrium via radial basis functions collocation and moving least squares approximation methods
journal, June 2016


Computation of fixed boundary tokamak equilibria using a method based on approximate particular solutions
journal, September 2015


Solution of Grad–Shafranov equation by the method of fundamental solutions
journal, February 2014


Accurate derivative evaluation for any Grad–Shafranov solver
journal, January 2016


To CG or to HDG: A Comparative Study
journal, July 2011

  • Kirby, Robert M.; Sherwin, Spencer J.; Cockburn, Bernardo
  • Journal of Scientific Computing, Vol. 51, Issue 1
  • DOI: 10.1007/s10915-011-9501-7

Solving Dirichlet Boundary-value Problems on Curved Domains by Extensions from Subdomains
journal, January 2012

  • Cockburn, Bernardo; Solano, Manuel
  • SIAM Journal on Scientific Computing, Vol. 34, Issue 1
  • DOI: 10.1137/100805200

Unified Hybridization of Discontinuous Galerkin, Mixed, and Continuous Galerkin Methods for Second Order Elliptic Problems
journal, January 2009

  • Cockburn, Bernardo; Gopalakrishnan, Jayadeep; Lazarov, Raytcho
  • SIAM Journal on Numerical Analysis, Vol. 47, Issue 2
  • DOI: 10.1137/070706616

Optimal Isoparametric Finite Elements and Error Estimates for Domains Involving Curved Boundaries
journal, June 1986

  • Lenoir, M.
  • SIAM Journal on Numerical Analysis, Vol. 23, Issue 3
  • DOI: 10.1137/0723036

A finite element method for interface problems in domains with smooth boundaries and interfaces
journal, December 1996

  • Bramble, James H.; King, J. Thomas
  • Advances in Computational Mathematics, Vol. 6, Issue 1
  • DOI: 10.1007/BF02127700

Flow patterns around heart valves: A numerical method
journal, October 1972


A priori error analysis for HDG methods using extensions from subdomains to achieve boundary conformity
journal, July 2013


A projection-based error analysis of HDG methods
journal, March 2010

  • Cockburn, Bernardo; Gopalakrishnan, Jayadeep; Sayas, Francisco-Javier
  • Mathematics of Computation, Vol. 79, Issue 271
  • DOI: 10.1090/S0025-5718-10-02334-3

A Characterization of Hybridized Mixed Methods for Second Order Elliptic Problems
journal, January 2004

  • Cockburn, Bernardo; Gopalakrishnan, Jayadeep
  • SIAM Journal on Numerical Analysis, Vol. 42, Issue 1
  • DOI: 10.1137/S0036142902417893

Algorithm 949: MATLAB Tools for HDG in Three Dimensions
journal, June 2015

  • Fu, Zhixing; Gatica, Luis F.; Sayas, Francisco-javier
  • ACM Transactions on Mathematical Software, Vol. 41, Issue 3
  • DOI: 10.1145/2658992

Reduction of stiffness and mass matrices
journal, February 1965


To CG or to HDG: A Comparative Study in 3D
journal, July 2015

  • Yakovlev, Sergey; Moxey, David; Kirby, Robert M.
  • Journal of Scientific Computing, Vol. 67, Issue 1
  • DOI: 10.1007/s10915-015-0076-6

Conditions for superconvergence of HDG methods for second-order elliptic problems
journal, September 2012


A Novel Two-Grid Method for Semilinear Elliptic Equations
journal, January 1994

  • Xu, Jinchao
  • SIAM Journal on Scientific Computing, Vol. 15, Issue 1
  • DOI: 10.1137/0915016

Two-Grid Discretization Techniques for Linear and Nonlinear PDE s
journal, October 1996


Iterative Procedures for Nonlinear Integral Equations
journal, October 1965


Convergence Analysis for Anderson Acceleration
journal, January 2015

  • Toth, Alex; Kelley, C. T.
  • SIAM Journal on Numerical Analysis, Vol. 53, Issue 2
  • DOI: 10.1137/130919398

Anderson Acceleration for Fixed-Point Iterations
journal, January 2011

  • Walker, Homer F.; Ni, Peng
  • SIAM Journal on Numerical Analysis, Vol. 49, Issue 4
  • DOI: 10.1137/10078356X

“One size fits all” analytic solutions to the Grad–Shafranov equation
journal, March 2010

  • Cerfon, Antoine J.; Freidberg, Jeffrey P.
  • Physics of Plasmas, Vol. 17, Issue 3
  • DOI: 10.1063/1.3328818

Analytical solutions to the Grad–Shafranov equation for tokamak equilibrium with dissimilar source functions
journal, September 1999


A family of analytic equilibrium solutions for the Grad–Shafranov equation
journal, November 2007

  • Guazzotto, L.; Freidberg, J. P.
  • Physics of Plasmas, Vol. 14, Issue 11
  • DOI: 10.1063/1.2803759

Noncircular, finite aspect ratio, local equilibrium model
journal, April 1998

  • Miller, R. L.; Chu, M. S.; Greene, J. M.
  • Physics of Plasmas, Vol. 5, Issue 4
  • DOI: 10.1063/1.872666

Axisymmetric MHD equilibrium solver with bicubic Hermite elements
journal, March 1992


Hydromagnetic equilibria and force-free fields
journal, September 1958


ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria
journal, May 2015


A mimetic spectral element solver for the Grad–Shafranov equation
journal, July 2016


Works referencing / citing this record: