skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach

Abstract

We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. Key to the efficiency of the method is a technique for fast computation of the product of the blob-blob mobility matrix and a vector. For unbounded suspensions, we rely on existing analytical expressions for the Rotne-Prager-Yamakawa tensor combined with a fast multipole method (FMM) to obtain linear scaling in the number of particles. For suspensions sedimented against a single no-slip boundary, we use a direct summation on a graphical processing unit (GPU), which gives quadratic asymptotic scaling with the number of particles. For fully confined domains, such as periodicmore » suspensions or suspensions confined in slit and square channels, we extend a recently developed rigid-body immersed boundary method by B. Kallemov, A. P. S. Bhalla, B. E. Griffith, and A. Donev (Commun. Appl. Math. Comput. Sci. 11 (2016), no. 1, 79-141) to suspensions of freely moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid-body equations converges in a bounded number of iterations regardless of the system size. In our approach, each iteration only requires a few cycles of a geometric multigrid solver for the Poisson equation, and an application of the block-diagonal preconditioner, leading to linear scaling with the number of particles. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.« less

Authors:
 [1];  [2];  [1];  [3];  [3];  [1]
  1. New York Univ., New York, NY (United States)
  2. New York Univ., New York, NY (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
  3. Univ. of North Carolina, Chapel Hill, NC (United States)
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
National Science Foundation (NSF); USDOE
OSTI Identifier:
1393005
Grant/Contract Number:  
AC02-05CH11231
Resource Type:
Accepted Manuscript
Journal Name:
Communications in Applied Mathematics and Computational Science
Additional Journal Information:
Journal Volume: 11; Journal Issue: 2; Journal ID: ISSN 1559-3940
Publisher:
Mathematical Sciences Publishers
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; Stokes flow; colloidal suspensions; Stokesian dynamics; immersed boundary method

Citation Formats

Usabiaga, Florencio Balboa, Kallemov, Bakytzhan, Delmotte, Blaise, Bhalla, Amneet Pal Singh, Griffith, Boyce E., and Donev, Aleksandar. Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach. United States: N. p., 2016. Web. doi:10.2140/camcos.2016.11.217.
Usabiaga, Florencio Balboa, Kallemov, Bakytzhan, Delmotte, Blaise, Bhalla, Amneet Pal Singh, Griffith, Boyce E., & Donev, Aleksandar. Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach. United States. doi:10.2140/camcos.2016.11.217.
Usabiaga, Florencio Balboa, Kallemov, Bakytzhan, Delmotte, Blaise, Bhalla, Amneet Pal Singh, Griffith, Boyce E., and Donev, Aleksandar. Tue . "Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach". United States. doi:10.2140/camcos.2016.11.217. https://www.osti.gov/servlets/purl/1393005.
@article{osti_1393005,
title = {Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach},
author = {Usabiaga, Florencio Balboa and Kallemov, Bakytzhan and Delmotte, Blaise and Bhalla, Amneet Pal Singh and Griffith, Boyce E. and Donev, Aleksandar},
abstractNote = {We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. Key to the efficiency of the method is a technique for fast computation of the product of the blob-blob mobility matrix and a vector. For unbounded suspensions, we rely on existing analytical expressions for the Rotne-Prager-Yamakawa tensor combined with a fast multipole method (FMM) to obtain linear scaling in the number of particles. For suspensions sedimented against a single no-slip boundary, we use a direct summation on a graphical processing unit (GPU), which gives quadratic asymptotic scaling with the number of particles. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently developed rigid-body immersed boundary method by B. Kallemov, A. P. S. Bhalla, B. E. Griffith, and A. Donev (Commun. Appl. Math. Comput. Sci. 11 (2016), no. 1, 79-141) to suspensions of freely moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid-body equations converges in a bounded number of iterations regardless of the system size. In our approach, each iteration only requires a few cycles of a geometric multigrid solver for the Poisson equation, and an application of the block-diagonal preconditioner, leading to linear scaling with the number of particles. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.},
doi = {10.2140/camcos.2016.11.217},
journal = {Communications in Applied Mathematics and Computational Science},
number = 2,
volume = 11,
place = {United States},
year = {2016},
month = {1}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 15 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Simple and efficient representations for the fundamental solutions of Stokes flow in a half-space
journal, July 2015

  • Gimbutas, Z.; Greengard, L.; Veerapaneni, S.
  • Journal of Fluid Mechanics, Vol. 776
  • DOI: 10.1017/jfm.2015.302

Force-coupling method for particulate two-phase flow: Stokes flow
journal, January 2003


Fluctuating force-coupling method for simulations of colloidal suspensions
journal, July 2014


Analysis of the Brinkman equation as a model for flow in porous media
journal, January 1987

  • Durlofsky, L.; Brady, J. F.
  • Physics of Fluids, Vol. 30, Issue 11
  • DOI: 10.1063/1.866465

Calculation of Hydrodynamic Properties of Globular Proteins from Their Atomic-Level Structure
journal, February 2000


Inducing Propulsion of Colloidal Dimers by Breaking the Symmetry in Electrohydrodynamic Flow
journal, November 2015


Ewald sum of the Rotne–Prager tensor
journal, August 1986

  • Beenakker, C. W. J.
  • The Journal of Chemical Physics, Vol. 85, Issue 3
  • DOI: 10.1063/1.451199

Brownian Motion of Arbitrarily Shaped Particles in Two Dimensions
journal, November 2014

  • Chakrabarty, Ayan; Konya, Andrew; Wang, Feng
  • Langmuir, Vol. 30, Issue 46
  • DOI: 10.1021/la5037053

Generalization of the Rotne–Prager–Yamakawa mobility and shear disturbance tensors
journal, August 2013

  • Wajnryb, Eligiusz; Mizerski, Krzysztof A.; Zuk, Pawel J.
  • Journal of Fluid Mechanics, Vol. 731
  • DOI: 10.1017/jfm.2013.402

Slow viscous motion of a sphere parallel to a plane wall—I Motion through a quiescent fluid
journal, April 1967


A fast multipole method for the Rotne–Prager–Yamakawa tensor and its applications
journal, February 2013

  • Liang, Zhi; Gimbutas, Zydrunas; Greengard, Leslie
  • Journal of Computational Physics, Vol. 234
  • DOI: 10.1016/j.jcp.2012.09.021

Sedimentation of hard-sphere suspensions at low Reynolds number
journal, January 1999


Staggered Schemes for Fluctuating Hydrodynamics
journal, January 2012

  • Balboa, Florencio; Bell, John B.; Delgado-Buscalioni, Rafael
  • Multiscale Modeling & Simulation, Vol. 10, Issue 4
  • DOI: 10.1137/120864520

Application of the Langevin equation to fluid suspensions
journal, December 1975


The immersed boundary method
journal, January 2002


Efficient Variable-Coefficient Finite-Volume Stokes Solvers
journal, November 2014

  • Cai, Mingchao; Nonaka, Andy; Bell, John B.
  • Communications in Computational Physics, Vol. 16, Issue 5
  • DOI: 10.4208/cicp.070114.170614a

Hydrodynamic transport coefficients of random dispersions of hard spheres
journal, September 1990

  • Ladd, Anthony J. C.
  • The Journal of Chemical Physics, Vol. 93, Issue 5
  • DOI: 10.1063/1.458830

Calculation of the resistance and mobility functions for two unequal rigid spheres in low-Reynolds-number flow
journal, February 1984


Simulating Brownian suspensions with fluctuating hydrodynamics
journal, December 2015

  • Delmotte, Blaise; Keaveny, Eric E.
  • The Journal of Chemical Physics, Vol. 143, Issue 24
  • DOI: 10.1063/1.4938173

Brownian dynamics without Green's functions
journal, April 2014

  • Delong, Steven; Usabiaga, Florencio Balboa; Delgado-Buscalioni, Rafael
  • The Journal of Chemical Physics, Vol. 140, Issue 13
  • DOI: 10.1063/1.4869866

The method of regularized Stokeslets in three dimensions:  Analysis, validation, and application to helical swimming
journal, March 2005

  • Cortez, Ricardo; Fauci, Lisa; Medovikov, Alexei
  • Physics of Fluids, Vol. 17, Issue 3
  • DOI: 10.1063/1.1830486

Simulation of hydrodynamically interacting particles near a no-slip boundary
journal, November 2007

  • Swan, James W.; Brady, John F.
  • Physics of Fluids, Vol. 19, Issue 11
  • DOI: 10.1063/1.2803837

Brownian Dynamics Simulation of Rigid Particles of Arbitrary Shape in External Fields
journal, December 2002


Living Crystals of Light-Activated Colloidal Surfers
journal, January 2013


Fast Ewald summation for Stokesian particle suspensions: FAST EWALD SUMMATION FOR STOKESIAN PARTICLE SUSPENSIONS
journal, September 2014

  • Af Klinteberg, Ludvig; Tornberg, Anna-Karin
  • International Journal for Numerical Methods in Fluids, Vol. 76, Issue 10
  • DOI: 10.1002/fld.3953

Particle motion between parallel walls: Hydrodynamics and simulation
journal, October 2010

  • Swan, James W.; Brady, John F.
  • Physics of Fluids, Vol. 22, Issue 10
  • DOI: 10.1063/1.3487748

A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles
journal, December 1949

  • Brinkman, H. C.
  • Flow, Turbulence and Combustion, Vol. 1, Issue 1
  • DOI: 10.1007/BF02120313

Brownian dynamics of confined rigid bodies
journal, October 2015

  • Delong, Steven; Balboa Usabiaga, Florencio; Donev, Aleksandar
  • The Journal of Chemical Physics, Vol. 143, Issue 14
  • DOI: 10.1063/1.4932062

Variational Treatment of Hydrodynamic Interaction in Polymers
journal, June 1969

  • Rotne, Jens; Prager, Stephen
  • The Journal of Chemical Physics, Vol. 50, Issue 11
  • DOI: 10.1063/1.1670977

A highly accurate boundary treatment for confined Stokes flow
journal, August 2012


Hydrodynamic interactions in a suspension of spherical particles
journal, April 1988

  • Ladd, Anthony J. C.
  • The Journal of Chemical Physics, Vol. 88, Issue 8
  • DOI: 10.1063/1.454658

Block Preconditioners Based on Approximate Commutators
journal, January 2006

  • Elman, Howard; Howle, Victoria E.; Shadid, John
  • SIAM Journal on Scientific Computing, Vol. 27, Issue 5
  • DOI: 10.1137/040608817

An adaptive, formally second order accurate version of the immersed boundary method
journal, April 2007

  • Griffith, Boyce E.; Hornung, Richard D.; McQueen, David M.
  • Journal of Computational Physics, Vol. 223, Issue 1
  • DOI: 10.1016/j.jcp.2006.08.019

Fast Computation of Many-Particle Hydrodynamic and Electrostatic Interactions in a Confined Geometry
journal, April 2007


Flow-induced nonequilibrium self-assembly in suspensions of stiff, apolar, active filaments
journal, January 2016

  • Pandey, Ankita; Sunil Kumar, P. B.; Adhikari, R.
  • Soft Matter, Vol. 12, Issue 44
  • DOI: 10.1039/C6SM02104B

Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions
journal, February 2016


Slender-body theory for slow viscous flow
journal, June 1976


Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. I. Algorithmic details
journal, January 2005

  • Donev, Aleksandar; Torquato, Salvatore; Stillinger, Frank H.
  • Journal of Computational Physics, Vol. 202, Issue 2
  • DOI: 10.1016/j.jcp.2004.08.014

Brownian motion in confinement
journal, August 2003


Friction and mobility for colloidal spheres in Stokes flow near a boundary: The multipole method and applications
journal, February 2000

  • Cichocki, B.; Jones, R. B.; Kutteh, Ramzi
  • The Journal of Chemical Physics, Vol. 112, Issue 5
  • DOI: 10.1063/1.480894

NlogN method for hydrodynamic interactions of confined polymer systems: Brownian dynamics
journal, October 2006

  • Hernández-Ortiz, Juan P.; de Pablo, Juan J.; Graham, Michael D.
  • The Journal of Chemical Physics, Vol. 125, Issue 16
  • DOI: 10.1063/1.2358344

An immersed boundary method for rigid bodies
journal, January 2016

  • Kallemov, Bakytzhan; Bhalla, Amneet; Griffith, Boyce
  • Communications in Applied Mathematics and Computational Science, Vol. 11, Issue 1
  • DOI: 10.2140/camcos.2016.11.79

Dynamic density functional theory with hydrodynamic interactions and fluctuations
journal, June 2014

  • Donev, Aleksandar; Vanden-Eijnden, Eric
  • The Journal of Chemical Physics, Vol. 140, Issue 23
  • DOI: 10.1063/1.4883520

An immersed boundary method for Brownian dynamics simulation of polymers in complex geometries: Application to DNA flowing through a nanoslit with embedded nanopits
journal, January 2012

  • Zhang, Yu; de Pablo, Juan J.; Graham, Michael D.
  • The Journal of Chemical Physics, Vol. 136, Issue 1
  • DOI: 10.1063/1.3672103

The slow motion of a sphere through a viscous fluid towards a plane surface
journal, December 1961


The Rotne-Prager-Yamakawa approximation for periodic systems in a shear flow
journal, May 2014

  • Mizerski, Krzysztof A.; Wajnryb, Eligiusz; Zuk, Pawel J.
  • The Journal of Chemical Physics, Vol. 140, Issue 18
  • DOI: 10.1063/1.4871113

Stochastic simulations of DNA in flow: Dynamics and the effects of hydrodynamic interactions
journal, May 2002

  • Jendrejack, Richard M.; de Pablo, Juan J.; Graham, Michael D.
  • The Journal of Chemical Physics, Vol. 116, Issue 17
  • DOI: 10.1063/1.1466831

Brownian particles at different times scales: a new derivation of the Smoluchowski equation
journal, October 1992


Spectrally accurate fast summation for periodic Stokes potentials
journal, November 2010


Stochastic Eulerian Lagrangian methods for fluid–structure interactions with thermal fluctuations
journal, April 2011


Inertial coupling method for particles in an incompressible fluctuating fluid
journal, February 2014

  • Balboa Usabiaga, Florencio; Delgado-Buscalioni, Rafael; Griffith, Boyce E.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 269
  • DOI: 10.1016/j.cma.2013.10.029

Accelerated Stokesian Dynamics simulations
journal, November 2001


PyCUDA and PyOpenCL: A scripting-based approach to GPU run-time code generation
journal, March 2012


Collective Hydrodynamics of Swimming Microorganisms: Living Fluids
journal, January 2011


Modeling hydrodynamic self-propulsion with Stokesian Dynamics. Or teaching Stokesian Dynamics to swim
journal, July 2011

  • Swan, James W.; Brady, John F.; Moore, Rachel S.
  • Physics of Fluids, Vol. 23, Issue 7
  • DOI: 10.1063/1.3594790

Accelerated Stokesian dynamics: Brownian motion
journal, June 2003

  • Banchio, Adolfo J.; Brady, John F.
  • The Journal of Chemical Physics, Vol. 118, Issue 22
  • DOI: 10.1063/1.1571819

Propulsion of Microorganisms by Surface Distortions
journal, November 1996


Effect of confinement on DNA dynamics in microfluidic devices
journal, July 2003

  • Jendrejack, Richard M.; Schwartz, David C.; Graham, Michael D.
  • The Journal of Chemical Physics, Vol. 119, Issue 2
  • DOI: 10.1063/1.1575200

Tetratic order in the phase behavior of a hard-rectangle system
journal, February 2006


Large-scale simulation of steady and time-dependent active suspensions with the force-coupling method
journal, December 2015

  • Delmotte, Blaise; Keaveny, Eric E.; Plouraboué, Franck
  • Journal of Computational Physics, Vol. 302
  • DOI: 10.1016/j.jcp.2015.09.020

Dynamical simulations of sedimenting spheres
journal, February 1993

  • Ladd, Anthony J. C.
  • Physics of Fluids A: Fluid Dynamics, Vol. 5, Issue 2
  • DOI: 10.1063/1.858695

Dynamic simulation of hydrodynamically interacting suspensions
journal, October 1988

  • Brady, John F.; Phillips, Ronald J.; Lester, Julia C.
  • Journal of Fluid Mechanics, Vol. 195, Issue -1
  • DOI: 10.1017/S0022112088002411

The hydrodynamics of confined dispersions
journal, October 2011


A Gaussian-like immersed-boundary kernel with three continuous derivatives and improved translational invariance
journal, July 2016


Prediction of Hydrodynamic and Other Solution Properties of Rigid Proteins from Atomic- and Residue-Level Models
journal, August 2011


Modeling and simulation of active suspensions containing large numbers of interacting micro-swimmers
journal, June 2013