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:
-
- New York Univ. (NYU), NY (United States)
- 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}
}
Web of Science
Works referenced in this record:
Computation of MHD equilibrium of tokamak plasma
journal, March 1991
- Takeda, Tatsuoki; Tokuda, Shinji
- Journal of Computational Physics, Vol. 93, Issue 1
Axialsymmetrische magnetohydrodynamische Gleichgewichtskonfigurationen
journal, October 1957
- Lüst, R.; Schlüter, A.
- Zeitschrift für Naturforschung A, Vol. 12, Issue 10
The ITER divertor concept
journal, April 1995
- Janeschitz, G.; Borrass, K.; Federici, G.
- Journal of Nuclear Materials, Vol. 220-222
Engineering and design aspects related to the development of the ITER divertor
journal, March 1995
- Dietz, K.
- Fusion Engineering and Design, Vol. 27, Issue 1-2
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
Solving the Grad–Shafranov equation with spectral elements
journal, May 2014
- Howell, E. C.; Sovinec, C. R.
- Computer Physics Communications, Vol. 185, Issue 5
A mimetic spectral element solver for the Grad–Shafranov equation
journal, July 2016
- Palha, A.; Koren, B.; Felici, F.
- Journal of Computational Physics, Vol. 316
Quasi-static free-boundary equilibrium of toroidal plasma with CEDRES++: Computational methods and applications
journal, January 2015
- Heumann, H.; Blum, J.; Boulbe, C.
- Journal of Plasma Physics, Vol. 81, Issue 3
A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries
journal, April 2017
- Heumann, Holger; Rapetti, Francesca
- Journal of Computational Physics, Vol. 334
Application of Parallelized Multigrid Method to Solution of MHD Equilibrium With MPI
journal, March 2004
- Hanawa, T.; Ikuno, S.; Kamitani, A.
- IEEE Transactions on Magnetics, Vol. 40, Issue 2
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
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
ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria
journal, May 2015
- Lee, Jungpyo; Cerfon, Antoine
- Computer Physics Communications, Vol. 190
Meshless method for solving fixed boundary problem of plasma equilibrium
journal, July 2015
- Imazawa, Ryota; Kawano, Yasunori; Itami, Kiyoshi
- Journal of Computational Physics, Vol. 292
Numerical solution of ideal MHD equilibrium via radial basis functions collocation and moving least squares approximation methods
journal, June 2016
- Ghasemi, Maryam; Amrollahi, Reza
- Engineering Analysis with Boundary Elements, Vol. 67
Computation of fixed boundary tokamak equilibria using a method based on approximate particular solutions
journal, September 2015
- Nath, D.; Kalra, M. S.; Munshi, P.
- Computers & Mathematics with Applications, Vol. 70, Issue 5
Solution of Grad–Shafranov equation by the method of fundamental solutions
journal, February 2014
- Nath, D.; Kalra, M. S.
- Journal of Plasma Physics, Vol. 80, Issue 3
Accurate derivative evaluation for any Grad–Shafranov solver
journal, January 2016
- Ricketson, L. F.; Cerfon, A. J.; Rachh, M.
- Journal of Computational Physics, Vol. 305
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
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
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
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
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
Flow patterns around heart valves: A numerical method
journal, October 1972
- Peskin, Charles S.
- Journal of Computational Physics, Vol. 10, Issue 2
A priori error analysis for HDG methods using extensions from subdomains to achieve boundary conformity
journal, July 2013
- Cockburn, Bernardo; Qiu, Weifeng; Solano, Manuel
- Mathematics of Computation, Vol. 83, Issue 286
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
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
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
Reduction of stiffness and mass matrices
journal, February 1965
- Guyan, Robert J.
- AIAA Journal, Vol. 3, Issue 2
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
Conditions for superconvergence of HDG methods for second-order elliptic problems
journal, September 2012
- Cockburn, Bernardo; Qiu, Weifeng; Shi, Ke
- Mathematics of Computation, Vol. 81, Issue 279
A Novel Two-Grid Method for Semilinear Elliptic Equations
journal, January 1994
- Xu, Jinchao
- SIAM Journal on Scientific Computing, Vol. 15, Issue 1
Two-Grid Discretization Techniques for Linear and Nonlinear PDE s
journal, October 1996
- Xu, Jinchao
- SIAM Journal on Numerical Analysis, Vol. 33, Issue 5
Iterative Procedures for Nonlinear Integral Equations
journal, October 1965
- Anderson, Donald G.
- Journal of the ACM, Vol. 12, Issue 4
Convergence Analysis for Anderson Acceleration
journal, January 2015
- Toth, Alex; Kelley, C. T.
- SIAM Journal on Numerical Analysis, Vol. 53, Issue 2
Anderson Acceleration for Fixed-Point Iterations
journal, January 2011
- Walker, Homer F.; Ni, Peng
- SIAM Journal on Numerical Analysis, Vol. 49, Issue 4
“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
Analytical solutions to the Grad–Shafranov equation for tokamak equilibrium with dissimilar source functions
journal, September 1999
- Mc Carthy, P. J.
- Physics of Plasmas, Vol. 6, Issue 9
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
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
Axisymmetric MHD equilibrium solver with bicubic Hermite elements
journal, March 1992
- Lütjens, H.; Bondeson, A.; Roy, A.
- Computer Physics Communications, Vol. 69, Issue 2-3
Hydromagnetic equilibria and force-free fields
journal, September 1958
- Grad, Harold; Rubin, Hanan
- Journal of Nuclear Energy (1954), Vol. 7, Issue 3-4
ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria
journal, May 2015
- Lee, Jungpyo; Cerfon, Antoine
- Computer Physics Communications, Vol. 190
A mimetic spectral element solver for the Grad–Shafranov equation
journal, July 2016
- Palha, A.; Koren, B.; Felici, F.
- Journal of Computational Physics, Vol. 316
Works referencing / citing this record:
Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
text, January 2019
- Giacomini, Matteo; Sevilla, Ruben
- arXiv
A priori and a posteriori error analysis of an unfitted HDG method for semi-linear elliptic problems
preprint, January 2019
- Sánchez, Nestor; Sánchez-Vizuet, Tonatiuh; Solano, Manuel E.
- arXiv
A parallel cut-cell algorithm for the free-boundary Grad-Shafranov problem
preprint, January 2020
- Liu, Shuang; Tang, Qi; Tang, Xian-Zhu
- arXiv