کد مقاله | کد نشریه | سال انتشار | مقاله انگلیسی | نسخه تمام متن |
---|---|---|---|---|
501751 | 863641 | 2012 | 4 صفحه PDF | دانلود رایگان |

In the case of adiabatic electrons the gyrokinetic field equation for the electrostatic potential includes an averaging operator acting on flux surfaces. For realistic three-dimensional configurations, as e.g. in stellarator devices, the discretisation of this integro-differential equation leads to very large nearly dense matrices (full matrix approach) which typically cannot be stored in computer memory explicitly. A low memory consuming partly matrix-free approach, based on a preconditioned iterative matrix solver, has been developed where the Helmholtz part of the field equation is used in a matrix formulation while the averaging term is treated matrix-free. For matrices which could still be stored in memory explicitly, it is demonstrated that this approach is also much faster than the full matrix approach.
► Parallel solution of the gyrokinetic field equation with flux surface averaging.
► Standard matrix discretisation leads to full matrix.
► A new fast low memory solver was developed.
► This solver is based on a partly matrix-free iterative approach.
► It shows good scaling and is much faster than the standard approach.
Journal: Computer Physics Communications - Volume 183, Issue 2, February 2012, Pages 305–308