# zbMATH — the first resource for mathematics

Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme. (English) Zbl 1436.65202
Summary: We present a fast, high-order accurate and adaptive boundary integral scheme for solving the Stokes equations in complex – possibly nonsmooth – geometries in two dimensions. We apply the panel-based quadratures of Helsing and coworkers to evaluate to high accuracy the weakly-singular, hyper-singular, and super-singular integrals arising in the Nyström discretization, and also the near-singular integrals needed for flow and traction evaluation close to boundaries. The resulting linear system is solved iteratively via calls to a Stokes fast multipole method. We include an automatic algorithm to “panelize” a given geometry, and choose a panel order, which will efficiently approximate the density (and hence solution) to a user-prescribed tolerance. We show that this adaptive panel refinement procedure works well in practice even in the case of complex geometries with large number of corners, or close-to-touching smooth curves. In one example, for instance, a model 2D vascular network with 378 corners required less than 200K discretization points to obtain a 9-digit solution accuracy.
##### MSC:
 65N50 Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs 76D07 Stokes and related (Oseen, etc.) flows 76M15 Boundary element methods applied to problems in fluid mechanics
pvfmm
Full Text:
##### References:
 [1] Youngren, G. K.; Acrivos, A., Stokes flow past a particle of arbitrary shape: a numerical method of solution, J. Fluid Mech., 69, 377-403 (May 1975) [2] Youngren, G. K.; Acrivos, A., On the shape of a gas bubble in a viscous extensional flow, J. Fluid Mech., 76, 433-442 (August 1976) [3] Sorgentone, Chiara; Tornberg, Anna-Karin, A highly accurate boundary integral equation method for surfactant-laden drops in 3D, J. Comput. Phys., 360, 167-191 (2018) · Zbl 1391.76456 [4] Rahimian, Abtin; Veerapaneni, Shravan K.; Zorin, Denis; Biros, George, Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions, J. Comput. Phys., 298, 766-786 (2015) · Zbl 1349.76947 [5] Bryngelson, Spencer H.; Freund, Jonathan B., Global stability of flowing red blood cell trains, Phys. Rev. Fluids, 3, 7, Article 073101 pp. (2018) [6] Greengard, Leslie; Kropinski, Mary Catherine; Mayo, Anita, Integral equation methods for Stokes flow and isotropic elasticity in the plane, J. Comput. Phys., 125, 2, 403-414 (1996) · Zbl 0847.76066 [7] Ying, L.; Biros, G.; Zorin, D., A kernel-independent adaptive fast multipole method in two and three dimensions, J. Comput. Phys., 196, 2, 591-626 (2004) · Zbl 1053.65095 [8] Wang, Haitao; Lei, Ting; Li, Jin; Huang, Jingfang; Yao, Zhenhan, A parallel fast multipole accelerated integral equation scheme for 3D Stokes equations, Int. J. Numer. Methods Eng., 70, 7, 812-839 (2007) · Zbl 1194.76221 [9] Tornberg, Anna-Karin; Greengard, Leslie, A fast multipole method for the three-dimensional Stokes equations, J. Comput. Phys., 227, 3, 1613-1619 (2008) · Zbl 1290.76116 [10] Gimbutas, Zydrunas; Greengard, Leslie, Computational software: simple FMM libraries for electrostatics, slow viscous flow, and frequency-domain wave propagation, Commun. Comput. Phys., 18, 2, 516-528 (2015) · Zbl 1388.65003 [11] Malhotra, Dhairya; Biros, George, PVFMM: a parallel kernel independent FMM for particle and volume potentials, Commun. Comput. Phys., 18, 3, 808-830 (2015) · Zbl 1388.65169 [12] Saintillan, David; Darve, Eric; Shaqfeh, Eric S. G., A smooth particle-mesh Ewald algorithm for Stokes suspension simulations: the sedimentation of fibers, Phys. Fluids, 17, 3, Article 033301 pp. (2005) · Zbl 1187.76458 [13] Lindbo, Dag; Tornberg, Anna-Karin, Spectrally accurate fast summation for periodic Stokes potentials, J. Comput. Phys., 229, 23, 8994-9010 (2010) · Zbl 1282.76151 [14] Kumar, Amit; Graham, Michael D., Accelerated boundary integral method for multiphase flow in non-periodic geometries, J. Comput. Phys., 231, 20, 6682-6713 (2012) · Zbl 1284.76277 [15] af Klinteberg, Ludvig; Shamshirgar, Davoud Saffar; Tornberg, Anna-Karin, Fast Ewald summation for free-space Stokes potentials, Res. Math. Sci., 4, 1, 1 (2017) · Zbl 1360.65295 [16] Sangani, Ashok S.; Mo, Guobiao, An $$\mathcal{O}(\mathcal{N})$$ algorithm for Stokes and Laplace interactions of particles, Phys. Fluids, 8, 8, 1990-2010 (1996) · Zbl 1027.76627 [17] Zinchenko, A. Z.; Davis, R. H., An efficient algorithm for hydrodynamical interaction of many deformable drops, J. Comput. Phys., 157, 2, 539-587 (2000) · Zbl 0961.76055 [18] Wang, Xin; Kanapka, Joe; Ye, Wenjing; Aluru, Narayan R.; White, Jacob, Algorithms in FastStokes and its application to micromachined device simulation, IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst., 25, 2, 248-257 (2006) [19] Wang, Lei; Tlupova, Svetlana; Krasny, Robert, A treecode algorithm for 3D Stokeslets and stresslets, Adv. Appl. Math. Mech., 11, 737-756 (2019) [20] Rahimian, Abtin; Lashuk, Ilya; Veerapaneni, Shravan; Chandramowlishwaran, Aparna; Malhotra, Dhairya; Moon, Logan; Sampath, Rahul; Shringarpure, Aashay; Vetter, Jeffrey; Vuduc, Richard; Zorin, Denis; Biros, George, Petascale direct numerical simulation of blood flow on 200k cores and heterogeneous architectures, (Proceedings of the 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis, SC ‘10 (2010)), 1-11 [21] Nazockdast, Ehssan; Rahimian, Abtin; Zorin, Denis; Shelley, Michael, A fast platform for simulating semi-flexible fiber suspensions applied to cell mechanics, J. Comput. Phys., 329, 173-209 (2017) · Zbl 1406.92032 [22] Yan, Wen; Corona, Eduardo; Malhotra, Dhairya; Veerapaneni, Shravan; Shelley, Michael, A scalable computational platform for particulate Stokes suspensions (2019), preprint [23] Hulme, S. Elizabeth; Shevkoplyas, Sergey S.; Apfeld, Javier; Fontana, Walter; Whitesides, George M., A microfabricated array of clamps for immobilizing and imaging c. elegans, Lab Chip, 7, 11, 1515-1523 (2007) [24] Ishikawa, Takuji; Pedley, T. J., The rheology of a semi-dilute suspension of swimming model micro-organisms, J. Fluid Mech., 588, 399-435 (2007) · Zbl 1141.76464 [25] Barnett, Alex; Wu, Bowei; Veerapaneni, Shravan, Spectrally accurate quadratures for evaluation of layer potentials close to the boundary for the 2D Stokes and Laplace equations, SIAM J. Sci. Comput., 37, 4, B519-B542 (2015) · Zbl 1433.65323 [26] Ojala, Rikard; Tornberg, Anna-Karin, An accurate integral equation method for simulating multi-phase Stokes flow, J. Comput. Phys., 298, 145-160 (2015) · Zbl 1349.76635 [27] Helsing, J.; Ojala, R., On the evaluation of layer potentials close to their sources, J. Comput. Phys., 227, 2899-2921 (2008) · Zbl 1135.65404 [28] Helsing, J., Integral equation methods for elliptic problems with boundary conditions of mixed type, J. Comput. Phys., 228, 8892-8907 (2009) · Zbl 1177.65176 [29] Ojala, R., A robust and accurate solver of Laplace’s equation with general boundary conditions on general domains in the plane, J. Comput. Math., 30, 4, 433-448 (2012) · Zbl 1274.65310 [30] Kropinski, M. C.A., An efficient numerical method for studying interfacial motion in two-dimensional creeping flows, J. Comput. Phys., 171, 2, 479-508 (2001) · Zbl 1047.76572 [31] Helsing, J.; Jiang, S., On integral equation methods for the first Dirichlet problem of the biharmonic and modified biharmonic equations in nonsmooth domains, SIAM J. Sci. Comput., 40, 4, A2609-A2630 (2018) · Zbl 1398.65353 [32] Klöckner, Andreas; Barnett, Alexander; Greengard, Leslie; O’Neil, Michael, Quadrature by expansion: a new method for the evaluation of layer potentials, J. Comput. Phys., 252, 332-349 (2013) · Zbl 1349.65094 [33] Barnett, Alex H., Evaluation of layer potentials close to the boundary for Laplace and Helmholtz problems on analytic planar domains, SIAM J. Sci. Comput., 36, 2, A427-A451 (2014) [34] Helsing, J.; Holst, A., Variants of an explicit kernel-split panel-based Nyström discretization scheme for Helmholtz boundary value problems, Adv. Comput. Math., 41, 3, 691-708 (2015) · Zbl 1319.65117 [35] Carvalho, Camille; Khatri, Shilpa; Kim, Arnold D., Asymptotic analysis for close evaluation of layer potentials, J. Comput. Phys., 355, 327-341 (2018) · Zbl 1380.65396 [36] Rahimian, Abtin; Barnett, Alex; Zorin, Denis, Ubiquitous evaluation of layer potentials using quadrature by kernel-independent expansion, BIT Numer. Math., 58, 2, 423-456 (2018) · Zbl 1395.65151 [37] af Klinteberg, Ludvig; Tornberg, Anna-Karin, Adaptive quadrature by expansion for layer potential evaluation in two dimensions, SIAM J. Sci. Comput., 40, 3, A1225-A1249 (2018) · Zbl 1446.65204 [38] Pérez-Arancibia, Carlos; Faria, Luiz M.; Turc, Catalin, Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D, J. Comput. Phys., 376, 411-434 (2019) · Zbl 1416.65491 [39] Veerapaneni, Shravan K.; Rahimian, Abtin; Biros, George; Zorin, Denis, A fast algorithm for simulating vesicle flows in three dimensions, J. Comput. Phys., 230, 14, 5610-5634 (2011) · Zbl 1419.76475 [40] Kita, Eisuke; Kamiya, Norio, Error estimation and adaptive mesh refinement in boundary element method, an overview, Eng. Anal. Bound. Elem., 25, 7, 479-495 (2001) · Zbl 1008.65089 [41] Heuer, N.; Stephan, E. P., The hp-version of the boundary element method on polygons, J. Integral Equ. Appl., 8, 2, 173-212 (1996) · Zbl 0888.65121 [42] Bantle, Markus; Funken, Stefan, Efficient and accurate implementation of hp-BEM for the Laplace operator in 2D, Appl. Numer. Math., 95, 51-61 (2015), Fourth Chilean Workshop on Numerical Analysis of Partial Differential Equations (WONAPDE 2013) · Zbl 1320.65185 [43] Helsing, Johan, Solving integral equations on piecewise smooth boundaries using the RCIP method: a tutorial, Abstr. Appl. Anal., 2013 (2013), Article ID 938167 · Zbl 1328.65271 [44] Manas, Rachh; Serkh, Kirill, On the solution of Stokes equation on regions with corners, Commun. Pure Appl. Math. (2017), submitted for publication [45] Pozrikidis, Constantine, Boundary Integral and Singularity Methods for Linearized Viscous Flow (1992), Cambridge University Press · Zbl 0772.76005 [46] Hsiao, G. C.; Wendland, W. L., Boundary Integral Equations, Applied Mathematical Sciences, vol. 164 (2008), Springer · Zbl 1157.65066 [47] Hebeker, F.-K., Efficient boundary element methods for three-dimensional exterior viscous flows, Numer. Methods Partial Differ. Equ., 2, 273-297 (1986) · Zbl 0645.76035 [48] Kress, Rainer, Linear Integral Equations, Applied Mathematical Sciences, vol. 82 (1999), Springer · Zbl 0920.45001 [49] Rachh, M., : parallel Fortran code for the biharmonic FMM in 2D (2012) [50] Greenbaum, A.; Greengard, L.; Mayo, A., On the numerical solution of the biharmonic equation in the plane, Physica D, 60, 1-4, 216-225 (1992) · Zbl 0824.65117 [51] Trefethen, Lloyd N., Approximation Theory and Approximation Practice, vol. 128 (2013), SIAM · Zbl 1264.41001 [52] af Klinteberg, L.; Tornberg, A-K., Adaptive quadrature by expansion for layer potential evaluation in two dimensions, SIAM J. Sci. Comput., 40, 3, A1225-A1249 (2018) · Zbl 1446.65204 [53] af Klinteberg, L.; Barnett, A. H., Accurate quadrature of nearly singular line integrals in two and three dimensions by singularity swapping, BIT Numer. Math. (2019), submitted for publication [54] Pan, V. Y., How bad are Vandermonde matrices?, SIAM J. Matrix Anal. Appl., 37, 2, 676-694 (2016) · Zbl 1382.15008 [55] Tagliasacchi, Andrea, code (2017) [56] Gopal, Abinand; Trefethen, Lloyd N., Solving Laplace problems with corner singularities via rational functions, SIAM J. Numer. Anal., 57, 5, 2074-2094 (2019) · Zbl 1431.65223 [57] Alvarez, Yolanda; Cederlund, Maria L.; Cottell, David C.; Bill, Brent R.; Ekker, Stephen C.; Torres-Vazquez, Jesus; Weinstein, Brant M.; Hyde, David R.; Vihtelic, Thomas S.; Kennedy, Breandan N., Genetic determinants of hyaloid and retinal vasculature in zebrafish, BMC Dev. Biol., 7, 1, 114 (2007) [58] Barnett, Alex, : a MATLAB/octave unix tool to record true memory and CPU usage vs time (2018) [59] Barnett, Alex H.; Marple, Gary R.; Veerapaneni, Shravan; Zhao, Lin, A unified integral equation scheme for doubly periodic Laplace and Stokes boundary value problems in two dimensions, Commun. Pure Appl. Math., 71, 11, 2334-2380 (2018) · Zbl 1442.65451 [60] Sangani, Ashok S.; Mo, Guobiao, Inclusion of lubrication forces in dynamic simulations, Phys. Fluids, 6, 5, 1653-1662 (1994) · Zbl 0829.76024 [61] Kabacaoğlu, Gökberk; Biros, George, Optimal design of deterministic lateral displacement device for viscosity-contrast-based cell sorting, Phys. Rev. Fluids, 3, 12, Article 124201 pp. (2018) [62] Bonnet, Marc; Liu, Ruowen; Veerapaneni, Shravan, Shape optimization of Stokesian peristaltic pumps using boundary integral methods (2019), preprint [63] Helsing, Johan; Jiang, Shidong, On integral equation methods for the first Dirichlet problem of the biharmonic and modified biharmonic equations in nonsmooth domains, SIAM J. Sci. Comput., 40, 4, A2609-A2630 (2018) · Zbl 1398.65353 [64] Marple, Gary; Barnett, Alexander H.; Gillman, Adrianna; Veerapaneni, Shravan K., A fast algorithm for simulating multiphase flows through periodic geometries of arbitrary shape, SIAM J. Sci. Comput., 38, 5, B740-B772 (2016) · Zbl 1381.76229 [65] Zhang, Yabin; Gillman, Adrianna, A fast direct solver for boundary value problems on locally perturbed geometries, J. Comput. Phys., 356, 356-371 (2018) · Zbl 1380.65451 [66] Liu, Y., Fast Multipole Boundary Element Method: Theory and Applications in Engineering (2009), Cambridge University Press
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.