Accuracy of Spectral and Finite Difference Schemes in 2D Advection Problems

In this paper we investigate the accuracy of two numerical procedures commonly used to solve 2D advection problems: spectral and finite difference (FD) schemes. These schemes are widely used, simulating, e.g., neutral and plasma flows. FD schemes have long been considered fast, relatively easy to implement, and applicable to complex geometries, but are somewhat inferior in accuracy compared to spectral schemes. Using two study cases at high Reynolds number, the merging of two equally signed Gaussian vortices in a periodic box and dipole interaction with a no-slip wall, we will demonstrate that the accuracy of FD schemes can be significantly improved if one is careful in choosing an appropriate FD scheme that reflects conservation properties of the nonlinear terms and in setting up the grid in accordance with the problem.