Hybrid parallelization of the XTOR-2F code for the simulation of two-fluid MHD instabilities in tokamaks

详细信息    查看全文

• 作者：Alain Marx ; ^{alain.marx@polytechnique.edu} ; Hinrich Lü ; tjens
• 关键词：Implicit ; Newton&ndash ; Krylov ; MPI/OpenMP ; SPIKE
• 刊名：Computer Physics Communications
• 出版年：2017
• 出版时间：March 2017
• 年：2017
• 卷：212
• 期：Complete
• 页码：90-99
• 全文大小：2339 K
• 卷排序：212

文摘

A hybrid MPI/OpenMP parallel version of the XTOR-2F code [Lütjens and Luciani, J. Comput. Phys. 229 (2010) 8130] solving the two-fluid MHD equations in full tokamak geometry by means of an iterative Newton–Krylov matrix-free method has been developed. The present work shows that the code has been parallelized significantly despite the numerical profile of the problem solved by XTOR-2F, i.e. a discretization with pseudo-spectral representations in all angular directions, the stiffness of the two-fluid stability problem in tokamaks, and the use of a direct LU decomposition to invert the physical pre-conditioner at every Krylov iteration of the solver. The execution time of the parallelized version is an order of magnitude smaller than the sequential one for low resolution cases, with an increasing speedup when the discretization mesh is refined. Moreover, it allows to perform simulations with higher resolutions, previously forbidden because of memory limitations.