English

A parallel cut-cell algorithm for the free-boundary Grad-Shafranov problem

Numerical Analysis 2021-08-31 v2 Numerical Analysis Computational Physics

Abstract

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad-Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. The algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad-Shafranov solver.

Keywords

Cite

@article{arxiv.2012.06015,
  title  = {A parallel cut-cell algorithm for the free-boundary Grad-Shafranov problem},
  author = {Shuang Liu and Qi Tang and Xian-Zhu Tang},
  journal= {arXiv preprint arXiv:2012.06015},
  year   = {2021}
}

Comments

24 pages, 11 figures, SIAM Journal on Scientific Computing

R2 v1 2026-06-23T20:53:19.301Z