A fast numerical method for computing doublyperiodic regularized Stokes flow in 3D
Abstract
A new numerical method for computing threedimensional Stokes flow driven by a doublyperiodic array of regularized forces is presented. The method is based on deriving an analytical representation of a regularized Green's function in Fourier space. Then only an inverse fast Fourier transform (inverse FFT) has to be computed to determine the fluid velocity on a grid in the physical domain. The velocity at other points can be interpolated from this grid. Accuracy is verified by comparing numerical results to a solution that is independent of the method. Although the regularized forces lead to a smooth velocity field, the Green's function may contain rapid transitions that are not captured properly on a coarse grid. In that case, an Ewald splitting technique is used to compute the gridresolved part of the flow using an inverse FFT and a sum in physical space for the localized part of the velocity. The splitting parameter can be chosen as small as a few grid cells, which makes the sum in physical space converge extremely fast. We present numerical examples that demonstrate that fact. In some cases, when the grid size is sufficiently small compared to the regularization parameter, the Ewald splitting is not needed.
 Publication:

Journal of Computational Physics
 Pub Date:
 February 2014
 DOI:
 10.1016/j.jcp.2013.10.032
 Bibcode:
 2014JCoPh.258....1C