is a state-of-the-art library for computing particle and volume potentials. For volume potentials, we use high order piecewise polynomial discretization on adaptive octrees and precomputed operators for computing near-interactions. Far-field interactions are computed using the KIFMM. It includes several new performance optimizations, support for GPU accelerators and distributed memory parallelism. It is one of the fastest libraries for computing solutions of elliptic PDEs on the unit cube with free-space and periodic boundary conditions.

Input vorticity field (left) and output velocity field (right) after convolution with Biot-Savert kernel.

Volume Integral Equation Solver

We use the PVFMM library along with iterative linear solvers to solve volume integral equation (VIE) formulations of variable coefficient elliptic PDEs.

Concentrated Vesicle Flows

Implemented algorithms for near-singular integration, collision handling and adaptive time-stepping to allow accurate, long time-scale simulations of concentrated vesicle suspensions in periodic flows.