Imperial College London

Dr Eric E Keaveny

Faculty of Natural SciencesDepartment of Mathematics

Reader in Applied Mathematics
 
 
 
//

Contact

 

+44 (0)20 7594 2780e.keaveny

 
 
//

Location

 

741Huxley BuildingSouth Kensington Campus

//

Summary

 

Publications

Citation

BibTex format

@article{Bao:2018:10.1016/j.jcp.2018.08.021,
author = {Bao, Y and Rachh, M and Keaveny, E and Greengard, L and Donev, A},
doi = {10.1016/j.jcp.2018.08.021},
journal = {Journal of Computational Physics},
pages = {1094--1119},
title = {A fluctuating boundary integral method for Brownian suspensions},
url = {http://dx.doi.org/10.1016/j.jcp.2018.08.021},
volume = {374},
year = {2018}
}

RIS format (EndNote, RefMan)

TY  - JOUR
AB - We present a fluctuating boundary integral method (FBIM) for overdamped Brownian Dynamics (BD) of two-dimensional periodic suspensions of rigid particles of complex shape immersed in a Stokes fluid. We develop a novel approach for generating Brownian displacements that arise in response to the thermal fluctuations in the fluid. Our approach relies on a first-kind boundary integral formulation of a mobility problem in which a random surface velocity is prescribed on the particle surface, with zero mean and covariance proportional to the Green's function for Stokes flow (Stokeslet). This approach yields an algorithm that scales linearly in the number of particles for both deterministic and stochastic dynamics, handles particles of complex shape, achieves high order of accuracy, and can be generalized to three dimensions and other boundary conditions. We show that Brownian displacements generated by our method obey the discrete fluctuation–dissipation balance relation (DFDB). Based on a recently-developed Positively Split Ewald method Fiore et al. (2017) [24], near-field contributions to the Brownian displacements are efficiently approximated by iterative methods in real space, while far-field contributions are rapidly generated by fast Fourier-space methods based on fluctuating hydrodynamics. FBIM provides the key ingredient for time integration of the overdamped Langevin equations for Brownian suspensions of rigid particles. We demonstrate that FBIM obeys DFDB by performing equilibrium BD simulations of suspensions of starfish-shaped bodies using a random finite difference temporal integrator.
AU - Bao,Y
AU - Rachh,M
AU - Keaveny,E
AU - Greengard,L
AU - Donev,A
DO - 10.1016/j.jcp.2018.08.021
EP - 1119
PY - 2018///
SN - 0021-9991
SP - 1094
TI - A fluctuating boundary integral method for Brownian suspensions
T2 - Journal of Computational Physics
UR - http://dx.doi.org/10.1016/j.jcp.2018.08.021
UR - http://arxiv.org/abs/1709.01480v2
UR - https://www.sciencedirect.com/science/article/pii/S0021999118305448?via%3Dihub
UR - http://hdl.handle.net/10044/1/64112
VL - 374
ER -