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.

[1]  Doron Levy,et al.  A Third-Order Semidiscrete Central Scheme for Conservation Laws and Convection-Diffusion Equations , 2000, SIAM J. Sci. Comput..

[2]  Doron Levy,et al.  A Third-Order Semi-Discrete Central Scheme for Conservation Laws and Convection-Diffusion Equations , 2000, math/0002133.

[3]  Thomas Hagstrom,et al.  An efficient spectral method for ordinary differential equations with rational function coefficients , 1996, Math. Comput..

[4]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[5]  P. Moin,et al.  Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows , 1997 .

[6]  Evangelos A. Coutsias,et al.  Fundamental interactions of vortical structures with boundary layers in two-dimensional flows , 1991 .

[7]  J. Rasmussen,et al.  Formation and temporal evolution of the Lamb-dipole , 1997 .

[8]  A. Arakawa.  Computational design for long-term numerical integration of the equations of fluid motion: two-dimen , 1997 .

[9]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[10]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .

[11]  S. Orszag,et al.  High-order splitting methods for the incompressible Navier-Stokes equations , 1991 .

[12]  Self organization in 2D circular shear layers , 1994 .

[13]  Keith Bergeron,et al.  Dynamical properties of forced shear layers in an annular geometry , 2000, Journal of Fluid Mechanics.

[14]  G. J. F. van Heijst,et al.  On Chaplygin's investigations of two-dimensional vortex structures in an inviscid fluid , 1994, Journal of Fluid Mechanics.

[15]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[16]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[17]  H. Panofsky,et al.  Atmospheric Turbulence: Models and Methods for Engineering Applications , 1984 .

[18]  B. Scott.  Three-Dimensional Computation of Drift Alfven Turbulence , 1997 .

[19]  T. Bohr,et al.  Vortex merging and spectral cascade in two-dimensional flows , 1996 .

[20]  P. Gresho.  Incompressible Fluid Dynamics: Some Fundamental Formulation Issues , 1991 .

[21]  R. LeVeque.  Numerical methods for conservation laws , 1990 .