Periodic Fast Multipole Method

by   Ruqi Pei, et al.

A new scheme is presented for imposing periodic boundary conditions on unit cells with arbitrary source distributions. We restrict our attention here to the Poisson, modified Helmholtz, Stokes and modified Stokes equations. The approach extends to the oscillatory equations of mathematical physics, including the Helmholtz and Maxwell equations, but we will address these in a companion paper, since the nature of the problem is somewhat different and includes the consideration of quasiperiodic boundary conditions and resonances. Unlike lattice sum-based methods, the scheme is insensitive to the unit cell's aspect ratio and is easily coupled to adaptive fast multipole methods (FMMs). Our analysis relies on classical "plane-wave" representations of the fundamental solution, and yields an explicit low-rank representation of the field due to all image sources beyond the first layer of neighboring unit cells. When the aspect ratio of the unit cell is large, our scheme can be coupled with the nonuniform fast Fourier transform (NUFFT) to accelerate the evaluation of the induced field. Its performance is illustrated with several numerial examples.


page 1

page 2

page 3

page 4


An Efficient Discontinuous Galerkin Scheme for Simulating Terahertz Photoconductive Devices with Periodic Nanostructures

Photoconductive devices (PCDs) enhanced with nanostructures have shown a...

Penalization method for the Navier-Stokes-Fourier system

We apply the method of penalization to the Dirichlet problem for the Nav...

Convergence study of IB methods for Stokes equations with non-periodic boundary conditions

Peskin's Immersed Boundary (IB) model and method are among the most popu...

Universal AMG Accelerated Embedded Boundary Method Without Small Cell Stiffness

We develop a universally applicable embedded boundary finite difference ...

Force Analysis for Interactions beyond the Closest Neighbor in a Periodic Structure

Periodic structures are a type of metamaterials in which their physical ...

An optimal complexity spectral method for Navier–Stokes simulations in the ball

We develop a spectral method for solving the incompressible generalized ...

Level-set based design of Wang tiles for modelling complex microstructures

Microstructural geometry plays a critical role in a response of heteroge...