Efficient FMM accelerated vortex methods in three dimensions via the Lamb-Helmholtz decomposition

  • Authors:

  • Nail A. Gumerov;Ramani Duraiswami

  • Affiliations:

  • Computer Science and UMIACS, University of Maryland, College Park, United States and Fantalgo, LLC, Elkridge, MD, United States;Computer Science and UMIACS, University of Maryland, College Park, United States and Fantalgo, LLC, Elkridge, MD, United States

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

Quantified Score

Hi-index 31.45

Visualization

Abstract

Vortex methods are used to efficiently simulate incompressible flows using Lagrangian techniques. Use of the FMM (Fast Multipole Method) allows considerable speed up of both velocity evaluation and vorticity evolution terms in these methods. Both equations require field evaluation of constrained (divergence free) vector valued quantities (velocity, vorticity) and cross terms from these. These are usually evaluated by performing several FMM accelerated sums of scalar harmonic functions. We present a formulation of vortex methods based on the Lamb-Helmholtz decomposition of the velocity in terms of two scalar potentials. In its original form, this decomposition is not invariant with respect to translation, violating a key requirement for the FMM. One of the key contributions of this paper is a theory for translation for this representation. The translation theory is developed by introducing ''conversion'' operators, which enable the representation to be restored in an arbitrary reference frame. Using this form, efficient vortex element computations can be made, which need evaluation of just two scalar harmonic FMM sums for evaluating the velocity and vorticity evolution terms. Details of the decomposition, translation and conversion formulae, and sample numerical results are presented.