A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers

  • Authors:

  • Wim M. van Rees;Anthony Leonard;D. I. Pullin;Petros Koumoutsakos

  • Affiliations:

  • Chair of Computational Science, ETH Zürich, CH-8092, Switzerland;Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena, CA 91125, USA;Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena, CA 91125, USA;Chair of Computational Science, ETH Zürich, CH-8092, Switzerland

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.47

Visualization

Abstract

We present a validation study for the hybrid particle-mesh vortex method against a pseudo-spectral method for the Taylor-Green vortex at Re"@C=1600 as well as in the collision of two antiparallel vortex tubes at Re"@C=10,000. In this study we present diagnostics such as energy spectra and enstrophy as computed by both methods as well as point-wise comparisons of the vorticity field. Using a fourth order accurate kernel for interpolation between the particles and the mesh, the results of the hybrid vortex method and of the pseudo-spectral method agree well in both flow cases. For the Taylor-Green vortex, the vorticity contours computed by both methods around the time of the energy dissipation peak overlap. The energy spectrum shows that only the smallest length scales in the flow are not captured by the vortex method. In the second flow case, where we compute the collision of two anti-parallel vortex tubes at Reynolds number 10,000, the vortex method results and the pseudo-spectral method results are in very good agreement up to and including the first reconnection of the tubes. The maximum error in the effective viscosity is about 2.5% for the vortex method and about 1% for the pseudo-spectral method. At later times the flows computed with the different methods show the same qualitative features, but the quantitative agreement on vortical structures is lost.