View in EDS

Accurate calculation of the gradients of the equilibrium poloidal flux in tokamaks.

Saved in:
Bibliographic Details
Title: Accurate calculation of the gradients of the equilibrium poloidal flux in tokamaks.
Authors: Woo, M.1 (AUTHOR) mhwoo@kfe.re.kr, Jo, G.1 (AUTHOR) jogahyung@kfe.re.kr, Park, B.H.1 (AUTHOR) bhpark@kfe.re.kr, Aydemir, A.Y.1 (AUTHOR) yaydemir@gmail.com, Kim, J.-H1 (AUTHOR) yegakjh@kfe.re.kr
Source: Computer Physics Communications. Apr2026, Vol. 321, pN.PAG-N.PAG. 1p.
Subjects: Tokamaks, Derivatives (Mathematics), Partial differential equations, Frequency-domain analysis, Boundary element methods, Toroidal plasma, Fusion reactors
Abstract: This paper presents a novel method for calculating the first, second, and third derivatives of the equilibrium poloidal flux in different directions in tokamaks. The method is implemented in a new code called Equilibrium Derivative in Arbitrary Mesh (EDAM) which is designed for practical fusion applications. The spectral method is adopted along the boundary with evenly spaced angles, while unstructured triangular meshes are used inside the computational domain. A new boundary integral equation (BIE) is derived and solved numerically to obtain the first and higher-order derivatives at the boundary. Using GS equation, linear partial differential equations for the first and higher-order flux derivatives are then constructed and solved. Validation is performed using an analytical equilibrium constructed by Cicogna, which describes D-shaped plasmas with steep profiles near the boundary. The code demonstrates similar convergence rates for the first and higher-order derivatives, achieving second order accuracy. This new method has significant potential for practical fusion simulations, providing derivatives up to the third order with the required accuracy and precisely given values at any nodal points of the unstructured mesh. [ABSTRACT FROM AUTHOR]
Copyright of Computer Physics Communications is the property of Elsevier B.V. and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Engineering Source
Description
Abstract:This paper presents a novel method for calculating the first, second, and third derivatives of the equilibrium poloidal flux in different directions in tokamaks. The method is implemented in a new code called Equilibrium Derivative in Arbitrary Mesh (EDAM) which is designed for practical fusion applications. The spectral method is adopted along the boundary with evenly spaced angles, while unstructured triangular meshes are used inside the computational domain. A new boundary integral equation (BIE) is derived and solved numerically to obtain the first and higher-order derivatives at the boundary. Using GS equation, linear partial differential equations for the first and higher-order flux derivatives are then constructed and solved. Validation is performed using an analytical equilibrium constructed by Cicogna, which describes D-shaped plasmas with steep profiles near the boundary. The code demonstrates similar convergence rates for the first and higher-order derivatives, achieving second order accuracy. This new method has significant potential for practical fusion simulations, providing derivatives up to the third order with the required accuracy and precisely given values at any nodal points of the unstructured mesh. [ABSTRACT FROM AUTHOR]
ISSN:00104655
DOI:10.1016/j.cpc.2026.110022