Search

We develop a linearly-scaling variant of the Force Coupling Method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly-periodic geometry with either a single bottom wall or two walls (slit channel) in the aperiodic direction. Our spectrally-accurate Stokes solver uses the Fast Fourier Transform (FFT) in the periodic $xy$ plane and Chebyshev polynomials in the aperiodic $z$ direction normal to the wall(s). We decompose the problem into two problems. The first is a doubly-periodic subproblem in the presence of particles (source terms) with free-space boundary conditions in the $z$ direction, which we solve by borrowing ideas from a recent method for rapid evaluation of electrostatic interactions in doubly-periodic geometries [O. Maxian, R. P. Pel\'aez, L. Greengard and A. Donev, J. Chem. Phys. 154, 204107 (2021)]. The second is a correction subproblem to impose the boundary conditions on the wall(s). Instead of the traditional Gaussian kernel, we use the exponential of a semicircle kernel to model the source terms (body force) due to the presence of particles, and provide optimum values for the kernel parameters that ensure a given hydrodynamic radius with at least two digits of accuracy and rotational and translational invariance. The computation time of our solver, which is implemented in graphical processing units, scales linearly with the number of particles, and allows computations with about a million particles in less than a second for a sedimented layer of colloidal microrollers. We find that in a slit channel, a driven dense suspension of microrollers maintains the same two-layer structure as above a single wall, but moves at a substantially lower collective speed due to increased confinement.

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the Fluctuating Immersed Boundary (FIB) approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both cases the Immersed Boundary (IB) method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long range electrostatic interactions are computed by solving the Poisson equation, with short range corrections included using a novel immersed-boundary variant of the classical Particle-Particle Particle-Mesh (P3M) technique. Also included is a short range repulsive force based on the Weeks-Chandler-Andersen (WCA) potential. The new methodology is validated by comparison to Debye-H{\"u}ckel theory for ion-ion pair correlation functions, and Debye-H{\"u}ckel-Onsager theory for conductivity, including the Wein effect for strong electric fields. In each case good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid., Comment: 30 pages, 12 figures, 2 tables