Unstructured finite element solvers in gyrokinetic turbulence simulations of burning plasmas

M. F. Adams, Y. Nishimura, Emden Henson Van Emden Henson

研究成果: Article同行評審


Many plasma physics simulation applications utilize an unstructured finite element method (FEM) discretization of an elliptic operator, of the form -∇2u + αu = f; when α is equal to zero a "pure" Laplacian or Poisson equation results and when α is greater than zero a Helmholtz equation is produced. Discretized equations of this form, often 2D, occur in many tokamak fusion plasma, or burning plasma, applications - from MHD to Gyrokinetic codes. This report investigates the performance characteristics of basic classes of linear solvers (ie, direct, one-level iterative, and multilevel iterative methods) on 2D unstructured FEM problems of the form -∇2u + αu = f, with both α = 0 and α ≠ 0. The purpose of this work is to provide computational physicists guidelines as to appropriate linear solvers for their problems via detailed performance analysis, in terms of both the scalability and the constants in the solution times. We show, as expected, almost perfect scalability of multilevel methods and quantify the solution costs on a common computational platform - the IBM SP Power3.

頁(從 - 到)410-424
期刊Journal of Physics: Conference Series
出版狀態Published - 2005 一月 1

All Science Journal Classification (ASJC) codes

  • 物理與天文學 (全部)


深入研究「Unstructured finite element solvers in gyrokinetic turbulence simulations of burning plasmas」主題。共同形成了獨特的指紋。