TY - CHAP
T1 - A data parallel pseudo-spectral semi-implicit magnetohydrodynamics code
T2 - High-Performance Computing and Networking
Y1 - 1997
A1 - Keppens, R.
A1 - Poedts, S.
A1 - Meijer, P. M.
A1 - Goedbloed, J. P.
ED - Hertzberger, B.
ED - Sloot, P.
AB - The set of eight nonlinear partial differential equations of magnetohydrodynamics (MHD) is used for time dependent simulations of three-dimensional (3D) fluid flow in a magnetic field. A data parallel code is presented, which integrates the MHD equations in cylindrical geometry, combining a semi-implicit time integration with a pseudo-spectral treatment of the poloidal and longitudinal directions. The semi-implicit method is devised to lift the severe CFL-condition imposed by the fastest waves. In the radial direction, we use centered finite differences on a staggered mesh. Together with the semi-implicit method, this leads to tridiagonal systems to be solved for each 2D Fourier mode. The parallelism is required to fully resolve small-scale dynamics in MHD simulations at affordable CPU costs. It is obtained by performing the 2D FFTs in an embarrassingly parallel way, and by solving the tridiagonal systems using a pipelined elimination algorithm. We discuss the scalability of the full code for a CM-Fortran implementation.
JF - High-Performance Computing and Networking
T3 - Lecture Notes in Computer Science
VL - 1225
SN - 0302-97433-540-62898-3
U5 - aae040c05d2543fd123eddd4e1a02f4e
ER -