pith. sign in

arxiv: 2605.30805 · v1 · pith:634MXJS3new · submitted 2026-05-29 · ⚛️ physics.flu-dyn · cs.NA· math.NA

A scalable Ewald-free BIE framework for periodic Stokes flow via hierarchical proxy sums

Tianyue Li , Dhairya Malhotra , Shravan Veerapaneni This is my paper

Pith reviewed 2026-06-28 21:26 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn cs.NAmath.NA
keywords Stokes flowboundary integral equationsperiodic geometriesproxy sourceshierarchical summationEwald-freeStokesletfast multipole method
0
0 comments X

The pith

The periodization precomputation for periodic Stokes flow depends only on box geometry and reuses across kernels for O(N) accuracy.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a boundary integral equation method for particulate Stokes flow in periodic geometries using only the free-space Green's function. It introduces a hierarchical proxy sum based on kernel-independent fast multipole method sources to handle far-field periodic images without Ewald summation or solving an extended system. The method ensures convergence of the sum via a net-force-zero condition on the particle surfaces. The key feature is that the periodization setup depends solely on the periodic box geometry, making it independent of the flow kernel and the specific surfaces inside the box. This allows the same precomputation to be applied to Stokeslet, stresslet, and rotlet kernels, achieving high-order accuracy with linear computational cost when paired with adaptive surface discretizations.

Core claim

Proxy sources placed on equivalent surfaces of the kernel-independent FMM form the auxiliary basis, and contributions from far image boxes are captured by a hierarchical proxy sum made absolutely convergent by a net-force-zero compatibility condition. The resulting periodization precomputation depends only on the periodic-box geometry, independent of the kernel and of the surfaces inside the box, and is reused verbatim across the Stokeslet, stresslet, and rotlet.

What carries the argument

hierarchical proxy sum: a sum over hierarchical proxy sources on equivalent surfaces that captures periodic image contributions and is rendered convergent by the net-force-zero condition

Load-bearing premise

The net-force-zero compatibility condition is sufficient to make the hierarchical proxy sum absolutely convergent for periodic pipes, wall-bounded doubly-periodic, and triply-periodic geometries.

What would settle it

Running the hierarchical proxy sum for a configuration with nonzero net force and observing whether the sum diverges or fails to converge to the expected periodic solution would test the assumption.

Figures

Figures reproduced from arXiv: 2605.30805 by Dhairya Malhotra, Shravan Veerapaneni, Tianyue Li.

Figure 1
Figure 1. Figure 1: Three representative periodic geometries considered in this work. (a) Singly-periodic: a converging [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Left: A cubic domain B tiled periodically in 2D. The contributions from image boxes beyond the nearest [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Left: Construction of a multipole expansion for a source box. The potential is first evaluated at [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left: M2L translation at level 0. The operator MM2L, 0 evaluates the potential u d 0 at the incoming check points x d 0 around the root box B from the level-0 image copies of the multipole proxy sources (shown in red), as given by eq. (21). Right: M2L translation at level 1, followed by L2L translation to level 0. The operator MM2L, 1 evaluates the potential u d 1 at the level-1 incoming check points x d 1… view at source ↗
Figure 5
Figure 5. Figure 5: Hierarchical computation of the periodic far field for 1D periodization. Once the multipole expansions [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Self-convergence in Nf for singly-, doubly-, and triply-periodic pressure-driven flow past a spherical obstacle in a cylindrical channel (left), wall-confined toroidal loops (middle), and spherical suspensions (right), for several panel counts Np. Np = 18, Nf = 96 serves as the reference solution in each case. Spectral convergence is observed before plateaus caused by saturation of error from parameter Np.… view at source ↗
Figure 7
Figure 7. Figure 7: Convergence of the periodized solver for the 25-sphere singly-periodic system. (a) Joint dependence of [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Precomputation time for the Stokes single-layer (S [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Disturbance-velocity streamlines in triply-periodic polydisperse sphere suspensions of 25 (left), 400 [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Weak scaling efficiency (a) and strong scaling efficiency (b) for setup time (blue circles) and per [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
read the original abstract

Particulate Stokes flow in confined, periodic geometries underlies a broad class of problems in biophysics, microfluidics, and the rheology of complex fluids. Boundary integral equation (BIE) methods are a natural tool for such problems, but existing periodization schemes rely either on periodic Green's functions, which are restrictive for complex confining geometries, or on free-space schemes that solve auxiliary proxy strengths alongside the surface densities in an extended linear system whose cost scales unfavorably in three dimensions. We present a BIE framework for three-dimensional particulate Stokes flow in periodic pipes with circular cross-sections, wall-bounded doubly-periodic, and triply-periodic geometries that uses only the free-space Green's function and avoids both Ewald summation and the extended linear system. Proxy sources placed on equivalent surfaces of the kernel-independent FMM (KIFMM) form the auxiliary basis, and contributions from far image boxes are captured by a hierarchical proxy sum made absolutely convergent by a net-force-zero compatibility condition. The resulting periodization precomputation depends only on the periodic-box geometry, independent of the kernel and of the surfaces inside the box, and is reused verbatim across the Stokeslet, stresslet, and rotlet. Combined with high-order adaptive surface discretizations, the method achieves high-order accuracy at $\mathcal{O}(N)$ cost with a single layer of image boxes in the near field. Numerical examples on dense polydisperse suspensions with thousands of particles and on flow through complex periodic channels, together with strong and weak scaling studies, demonstrate efficient performance on systems with millions of degrees of freedom on distributed-memory architectures.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a boundary integral equation framework for three-dimensional particulate Stokes flow in periodic pipes, wall-bounded doubly-periodic, and triply-periodic geometries. It employs only the free-space Green's function together with hierarchical proxy sums constructed from KIFMM-style equivalent surfaces; a net-force-zero compatibility condition is asserted to render the proxy sums absolutely convergent, enabling a kernel-independent periodization precomputation that is reused across the Stokeslet, stresslet, and rotlet and combined with high-order adaptive surface discretizations to achieve O(N) cost using only a single layer of image boxes in the near field. Numerical examples on dense polydisperse suspensions and complex channels plus scaling studies are offered in support.

Significance. If the central claims hold, the work supplies a practical route to high-order, scalable simulations of confined periodic Stokes flow that avoids both Ewald summation and the cost of extended linear systems. The kernel- and surface-independent precomputation, together with the reported strong/weak scaling to millions of degrees of freedom on distributed architectures, would constitute a genuine technical contribution to the field.

major comments (2)
  1. [Abstract] Abstract: the assertion that the net-force-zero compatibility condition renders the hierarchical proxy sum "absolutely convergent" (and thereby permits truncation after a single image layer) is load-bearing for the claimed high-order accuracy and kernel independence. No explicit proof, decay-rate analysis, or uniform error bound is supplied for the stresslet and rotlet (whose far-field decay differs from the Stokeslet), nor is absolute versus conditional convergence verified numerically for the three distinct geometries.
  2. [Numerical examples] Numerical examples section: the reported error studies and scaling plots do not isolate the truncation error incurred by the single-layer image-box approximation when the proxy sum is applied to the stresslet or rotlet; without such targeted tests the claim that geometry-dependent errors are absorbed by the precomputation cannot be assessed.
minor comments (2)
  1. The definition of the hierarchical proxy sum and the precise placement of equivalent surfaces would benefit from an additional schematic or pseudocode block for readers outside the KIFMM literature.
  2. A short table summarizing the precomputation cost (in terms of box geometry only) versus the per-solve cost would clarify the claimed separation of concerns.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the work's potential contribution and for the detailed, constructive comments. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the assertion that the net-force-zero compatibility condition renders the hierarchical proxy sum "absolutely convergent" (and thereby permits truncation after a single image layer) is load-bearing for the claimed high-order accuracy and kernel independence. No explicit proof, decay-rate analysis, or uniform error bound is supplied for the stresslet and rotlet (whose far-field decay differs from the Stokeslet), nor is absolute versus conditional convergence verified numerically for the three distinct geometries.

    Authors: We agree that the manuscript would benefit from an explicit decay-rate analysis and numerical verification to support the absolute-convergence claim for all kernels. The net-force-zero condition cancels the leading far-field term of the Stokeslet; the stresslet and rotlet possess faster natural decay, and the same compatibility condition extends the argument. We will add a dedicated subsection containing the decay-rate analysis, uniform error bounds, and numerical checks of absolute versus conditional convergence for the Stokeslet, stresslet, and rotlet across the three geometries. revision: yes

  2. Referee: [Numerical examples] Numerical examples section: the reported error studies and scaling plots do not isolate the truncation error incurred by the single-layer image-box approximation when the proxy sum is applied to the stresslet or rotlet; without such targeted tests the claim that geometry-dependent errors are absorbed by the precomputation cannot be assessed.

    Authors: We concur that isolating the single-layer truncation error for the stresslet and rotlet would strengthen the assessment of the precomputation. We will augment the numerical-examples section with targeted experiments that apply the proxy sum separately to each kernel, quantify the truncation error under the single-image-layer approximation, and confirm that geometry-dependent contributions remain absorbed by the precomputed corrections while high-order accuracy is retained. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation relies on external KIFMM and net-force-zero condition with independent numerical validation

full rationale

The paper's central claims rest on the mathematical assertion that net-force-zero renders the hierarchical proxy sum (built from KIFMM equivalent surfaces) absolutely convergent for the listed geometries, combined with reuse of a geometry-only precomputation across kernels. No step reduces by definition or construction to its own inputs; the compatibility condition is an external assumption whose sufficiency is tested via reported scaling studies rather than assumed. No self-citations are load-bearing, no parameters are fitted then renamed as predictions, and the periodization independence is presented as a consequence of the proxy construction rather than a tautology. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

Framework depends on kernel-independent FMM proxy placement and the net-force-zero condition for convergence; no free parameters are described in the abstract.

axioms (1)
  • domain assumption Net-force-zero compatibility condition ensures absolute convergence of the hierarchical proxy sum from far image boxes.
    Invoked explicitly to avoid Ewald summation while using only free-space kernels.
invented entities (1)
  • Hierarchical proxy sum no independent evidence
    purpose: Capture far-field periodic image contributions via KIFMM proxy sources without extending the linear system.
    Introduced as the core periodization mechanism.

pith-pipeline@v0.9.1-grok · 5836 in / 1352 out tokens · 32244 ms · 2026-06-28T21:26:56.427120+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Fast summation on rectangular cuboids with arbitrary periodicity in the DMK framework

    math.NA 2026-06 unverdicted novelty 6.0

    DMK extended to rectangular cuboids with arbitrary periodicity via localized octree evaluations on cubical tilings and Fourier-space root-level summation with truncated kernels for reduced periodicity.

Reference graph

Works this paper leans on

43 extracted references · 29 canonical work pages · cited by 1 Pith paper

  1. [1]

    Fast Ewald summation for Stokesian particle suspensions.International Journal for Numerical Methods in Fluids, 76(10):669–698, 2014

    Ludvig af Klinteberg and Anna-Karin Tornberg. Fast Ewald summation for Stokesian particle suspensions.International Journal for Numerical Methods in Fluids, 76(10):669–698, 2014. doi: 10.1002/fld.3953. 18

    work page doi:10.1002/fld.3953 2014

  2. [2]

    Fast Ewald sum- mation for free-space Stokes potentials.Research in the Mathematical Sciences, 4(1):1, 2017

    Ludvig af Klinteberg, Davoud Saffar Shamshirgar, and Anna-Karin Tornberg. Fast Ewald sum- mation for free-space Stokes potentials.Research in the Mathematical Sciences, 4(1):1, 2017. doi: 10.1186/s40687-016-0092-7

    work page doi:10.1186/s40687-016-0092-7 2017

  3. [3]

    Fast sum- mation of Stokes potentials using a new kernel-splitting in the DMK framework.arXiv preprint arXiv:2509.21471, 2025

    Ludvig af Klinteberg, Leslie Greengard, Shidong Jiang, and Anna-Karin Tornberg. Fast sum- mation of Stokes potentials using a new kernel-splitting in the DMK framework.arXiv preprint arXiv:2509.21471, 2025

    arXiv 2025

  4. [4]

    Active turbulence.An- nual Review of Condensed Matter Physics, 13(1):143–170, 2022

    Ricard Alert, Jaume Casademunt, and Jean-François Joanny. Active turbulence.An- nual Review of Condensed Matter Physics, 13(1):143–170, 2022. doi: 10.1146/ annurev-conmatphys-082321-035957

    2022

  5. [5]

    Fast Ewald summation for Stokes flow with arbitrary periodicity.Journal of Computational Physics, 493:112473, 2023

    Joar Bagge and Anna-Karin Tornberg. Fast Ewald summation for Stokes flow with arbitrary periodicity.Journal of Computational Physics, 493:112473, 2023. doi: 10.1016/j.jcp.2023.112473

    work page doi:10.1016/j.jcp.2023.112473 2023

  6. [6]

    A new integral representation for quasi-periodic scattering problems in two dimensions.BIT Numerical Mathematics, 51(1):67–90, 2011

    Alex Barnett and Leslie Greengard. A new integral representation for quasi-periodic scattering problems in two dimensions.BIT Numerical Mathematics, 51(1):67–90, 2011. doi: 10.1007/ s10543-010-0297-x

    2011

  7. [7]

    Barnett, Gary Marple, Shravan Veerapaneni, and Lin Zhao

    Alex H. Barnett, Gary Marple, Shravan Veerapaneni, and Lin Zhao. A unified integral equa- tion scheme for doubly periodic Laplace and Stokes boundary value problems in two dimensions. Communications on Pure and Applied Mathematics, 71(11):2334–2380, 2018

    2018

  8. [8]

    exponentialofsemicircle

    Alexander H. Barnett, Jeremy Magland, and Ludvig af Klinteberg. A parallel nonuniform fast Fouriertransformlibrarybasedonan“exponentialofsemicircle” kernel.SIAM Journal on Scientific Computing, 41(5):C479–C504, 2019. doi: 10.1137/18M120885X

    work page doi:10.1137/18m120885x 2019

  9. [9]

    Leonard Berman and Leslie Greengard

    C. Leonard Berman and Leslie Greengard. A renormalization method for the evaluation of lattice sums.Journal of Mathematical Physics, 35(11):6036–6048, 1994. doi: 10.1063/1.530726

    work page doi:10.1063/1.530726 1994

  10. [10]

    Min Hyung Cho and Alex H. Barnett. Robust fast direct integral equation solver for quasi-periodic scattering problems with a large number of layers.Optics Express, 23(2):1775–1799, 2015. doi: 10.1364/OE.23.001775

    work page doi:10.1364/oe.23.001775 2015

  11. [11]

    Darling and Dino Di Carlo

    Eric M. Darling and Dino Di Carlo. High-throughput assessment of cellular mechani- cal properties.Annual Review of Biomedical Engineering, 17:35–62, 2015. doi: 10.1146/ annurev-bioeng-071114-040545

    2015

  12. [12]

    X.-J. Fan, N. Phan-Thien, and R. Zheng. Completed double layer boundary element method for periodic suspensions.Zeitschrift für angewandte Mathematik und Physik, 49(2):167–193, 1998. doi: 10.1007/s000330050214

    work page doi:10.1007/s000330050214 1998

  13. [13]

    Foss and John F

    David R. Foss and John F. Brady. Brownian dynamics simulation of hard-sphere colloidal disper- sions.Journal of Rheology, 44(3):629–651, 2000. doi: 10.1122/1.551104

    work page doi:10.1122/1.551104 2000

  14. [14]

    Jonathan B. Freund. Numerical simulation of flowing blood cells.Annual Review of Fluid Mechan- ics, 46:67–95, 2014. doi: 10.1146/annurev-fluid-010313-141349

    work page doi:10.1146/annurev-fluid-010313-141349 2014

  15. [15]

    A fast direct solver for quasiperiodic scattering problems

    Adrianna Gillman and Alex Barnett. A fast direct solver for quasiperiodic scattering problems. Journal of Computational Physics, 248:309–322, 2013

    2013

  16. [16]

    Greengard and V

    L. Greengard and V. Rokhlin. A fast algorithm for particle simulations.Journal of Computational Physics, 73(2):325–348, 1987. doi: 10.1016/0021-9991(87)90140-9

    work page doi:10.1016/0021-9991(87)90140-9 1987

  17. [17]

    Integral equation methods for Stokes flow in doubly-periodic domains.Journal of Engineering Mathematics, 48(2):157–170, 2004

    Leslie Greengard and Mary Catherine Kropinski. Integral equation methods for Stokes flow in doubly-periodic domains.Journal of Engineering Mathematics, 48(2):157–170, 2004. doi: 10.1023/ B:ENGI.0000011923.59797.92. 19

    arXiv 2004

  18. [18]

    On the periodic fundamental solutions of the Stokes equations and their application to viscous flow past a cubic array of spheres.Journal of Fluid Mechanics, 5(3):317– 328, 1959

    Hidenori Hasimoto. On the periodic fundamental solutions of the Stokes equations and their application to viscous flow past a cubic array of spheres.Journal of Fluid Mechanics, 5(3):317– 328, 1959. doi: 10.1017/S0022112059000222

    work page doi:10.1017/s0022112059000222 1959

  19. [19]

    A dual-space multilevel kernel-splitting framework for discrete and continuous convolution.Communications on Pure and Applied Mathematics, 78(5):1086–1143,

    Shidong Jiang and Leslie Greengard. A dual-space multilevel kernel-splitting framework for discrete and continuous convolution.Communications on Pure and Applied Mathematics, 78(5):1086–1143,

  20. [20]

    doi: 10.1002/cpa.22240

    work page doi:10.1002/cpa.22240

  21. [21]

    McFaul, Simon P

    Chao Jin, Sarah M. McFaul, Simon P. Duffy, Xiaoyan Deng, Peyman Tavassoli, Peter C. Black, and Hongshen Ma. Technologies for label-free separation of circulating tumor cells: from historical foundations to recent developments.Lab on a Chip, 14(1):32–44, 2014. doi: 10.1039/C3LC50625H

    work page doi:10.1039/c3lc50625h 2014

  22. [22]

    Spectrally accurate fast summation for periodic Stokes potentials.Journal of Computational Physics, 229(23):8994–9010, 2010

    Dag Lindbo and Anna-Karin Tornberg. Spectrally accurate fast summation for periodic Stokes potentials.Journal of Computational Physics, 229(23):8994–9010, 2010. doi: 10.1016/j.jcp.2010. 08.026

    work page doi:10.1016/j.jcp.2010 2010

  23. [23]

    Fast and spectrally accurate Ewald summation for 2- periodic electrostatic systems.The Journal of Chemical Physics, 136(16):164111, 2012

    Dag Lindbo and Anna-Karin Tornberg. Fast and spectrally accurate Ewald summation for 2- periodic electrostatic systems.The Journal of Chemical Physics, 136(16):164111, 2012. doi: 10. 1063/1.4704177

    2012

  24. [24]

    Yuxiang Liu and Alex H. Barnett. Efficient numerical solution of acoustic scattering from doubly- periodic arrays of axisymmetric objects.Journal of Computational Physics, 324:226–245, 2016. doi: 10.1016/j.jcp.2016.08.011

    work page doi:10.1016/j.jcp.2016.08.011 2016

  25. [25]

    Scientific computing template library

    Dhairya Malhotra. Scientific computing template library. URLhttps://sctl.readthedocs.io/ en/latest/index.html

  26. [26]

    Efficient convergent boundary integral methods for slender bodies.Journal of Computational Physics, 503:112855, 2024

    Dhairya Malhotra and Alex Barnett. Efficient convergent boundary integral methods for slender bodies.Journal of Computational Physics, 503:112855, 2024. doi: 10.1016/j.jcp.2024.112855

    work page doi:10.1016/j.jcp.2024.112855 2024

  27. [27]

    PVFMM: A parallel kernel independent FMM for particle and volume potentials.Communications in Computational Physics, 18(3):808–830, 2015

    Dhairya Malhotra and George Biros. PVFMM: A parallel kernel independent FMM for particle and volume potentials.Communications in Computational Physics, 18(3):808–830, 2015. doi: 10.4208/cicp.020215.150515sw

    work page doi:10.4208/cicp.020215.150515sw 2015

  28. [28]

    Algorithm 967: A distributed-memory fast multipole method for volume potentials.ACM Transactions on Mathematical Software, 43(2):1–27, 2016

    Dhairya Malhotra and George Biros. Algorithm 967: A distributed-memory fast multipole method for volume potentials.ACM Transactions on Mathematical Software, 43(2):1–27, 2016. doi: 10. 1145/2898349

    2016

  29. [29]

    Marple, Alex Barnett, Adrianna Gillman, and Shravan Veerapaneni

    Gary R. Marple, Alex Barnett, Adrianna Gillman, and Shravan Veerapaneni. A fast algorithm for simulating multiphase flows through periodic geometries of arbitrary shape.SIAM Journal on Scientific Computing, 38(5):B740–B772, 2016. doi: 10.1137/15M1043066

    work page doi:10.1137/15m1043066 2016

  30. [30]

    Simulation methods for particulate flows and concentrated suspensions.Annual Review of Fluid Mechanics, 49:171–193, 2017

    Martin Maxey. Simulation methods for particulate flows and concentrated suspensions.Annual Review of Fluid Mechanics, 49:171–193, 2017. doi: 10.1146/annurev-fluid-122414-034408

    work page doi:10.1146/annurev-fluid-122414-034408 2017

  31. [31]

    Cytoplasmic flows as signatures for the mechanics of mitotic positioning.Molecular Biology of the Cell, 28(23): 3261–3270, 2017

    Ehssan Nazockdast, Abtin Rahimian, Daniel Needleman, and Michael Shelley. Cytoplasmic flows as signatures for the mechanics of mitotic positioning.Molecular Biology of the Cell, 28(23): 3261–3270, 2017. doi: 10.1091/mbc.E16-02-0108

    work page doi:10.1091/mbc.e16-02-0108 2017

  32. [32]

    The physics of dense suspen- sions.Annual Review of Condensed Matter Physics, 13(1):97–117, 2022

    Christopher Ness, Ryohei Seto, and Romain Mari. The physics of dense suspen- sions.Annual Review of Condensed Matter Physics, 13(1):97–117, 2022. doi: 10.1146/ annurev-conmatphys-031620-105938

    2022

  33. [33]

    A fast method for imposing periodic boundary conditions on arbitrarily-shaped lattices in two dimensions.Journal of Compu- tational Physics, 474:111792, 2023

    Ruqi Pei, Travis Askham, Leslie Greengard, and Shidong Jiang. A fast method for imposing periodic boundary conditions on arbitrarily-shaped lattices in two dimensions.Journal of Compu- tational Physics, 474:111792, 2023. doi: 10.1016/j.jcp.2022.111792

    work page doi:10.1016/j.jcp.2022.111792 2023

  34. [34]

    Pozrikidis.Boundary Integral and Singularity Methods for Linearized Viscous Flow

    C. Pozrikidis.Boundary Integral and Singularity Methods for Linearized Viscous Flow. Cambridge University Press, Cambridge, UK, 1992. 20

    1992

  35. [35]

    Pozrikidis

    C. Pozrikidis. Computation of Stokes flow in two-dimensional periodic systems.Journal of Fluid Mechanics, 324:107–131, 1996

    1996

  36. [36]

    David Saintillan, Eric Darve, and Eric S. G. Shaqfeh. A smooth particle-mesh Ewald algorithm for Stokes suspension simulations: The sedimentation of fibers.Physics of Fluids, 17(3):033301, 2005. doi: 10.1063/1.1862262

    work page doi:10.1063/1.1862262 2005

  37. [37]

    Timothy W. Secomb. Blood flow in the microcirculation.Annual Review of Fluid Mechanics, 49: 443–461, 2017. doi: 10.1146/annurev-fluid-010816-060302

    work page doi:10.1146/annurev-fluid-010816-060302 2017

  38. [38]

    Mu Wang and John F. Brady. Spectral Ewald acceleration of Stokesian dynamics for polydisperse suspensions.Journal of Computational Physics, 306:443–477, 2016. doi: 10.1016/j.jcp.2015.11.042

    work page doi:10.1016/j.jcp.2015.11.042 2016

  39. [39]

    Flexibly imposing periodicity in kernel independent FMM: A multipole-to-local operator approach.Journal of Computational Physics, 355:214–232, 2018

    Wen Yan and Michael Shelley. Flexibly imposing periodicity in kernel independent FMM: A multipole-to-local operator approach.Journal of Computational Physics, 355:214–232, 2018. doi: 10.1016/j.jcp.2017.11.012

    work page doi:10.1016/j.jcp.2017.11.012 2018

  40. [40]

    A scalable computational platform for particulate Stokes suspensions.Journal of Computational Physics, 416:109524, 2020

    Wen Yan, Eduardo Corona, Dhairya Malhotra, Shravan Veerapaneni, and Michael Shelley. A scalable computational platform for particulate Stokes suspensions.Journal of Computational Physics, 416:109524, 2020. doi: 10.1016/j.jcp.2020.109524

    work page doi:10.1016/j.jcp.2020.109524 2020

  41. [41]

    A kernel-independent adaptive fast multipole algo- rithm in two and three dimensions.Journal of Computational Physics, 196(2):591–626, 2004

    Lexing Ying, George Biros, and Denis Zorin. A kernel-independent adaptive fast multipole algo- rithm in two and three dimensions.Journal of Computational Physics, 196(2):591–626, 2004. doi: 10.1016/j.jcp.2003.11.021

    work page doi:10.1016/j.jcp.2003.11.021 2004

  42. [42]

    G. K. Youngren and A. Acrivos. Stokes flow past a particle of arbitrary shape: a numerical method of solution.Journal of Fluid Mechanics, 69(2):377–403, 1975. doi: 10.1017/S0022112075001486

    work page doi:10.1017/s0022112075001486 1975

  43. [43]

    A. A. Zick and G. M. Homsy. Stokes flow through periodic arrays of spheres.Journal of Fluid Mechanics, 115:13–26, 1982. doi: 10.1017/S0022112082000627. 21

    work page doi:10.1017/s0022112082000627 1982