A pseudospectral algorithm for three-dimensional magnetohydrodynamic simulation (Q761327)

An algorithm for the solution of the three-dimensional resistive magnetohydrodynamic equations in toroidal geometry is presented. The algorithm employs the pseudospectral method for approximation in the two periodic coordinates, and finite differences in the radial direction. Efficient fast Fourier transforms are used to communicate between configuration and Fourier space. Leapfrog time advancement is used for advective terms. Diffusion terms are treated implicitly to avoid severe time step restrictions. Sample cases are presented, and a comparison of the method with standard finite difference techniques is presented and discussed.