The hydrodynamic interactions are implemented as Ball-Melrose lubrication terms, as detailed in Ball and Melrose . The Ball-Melrose lubrication terms assume a linearisation of interparticle dissipative forces, as a result of low "Reynoldsnumber" flow of colloidal particles in a hydrodynamic medium. This allows for the force equation to be solved using a resistance matrix to determine the most divergent terms for an N-body system such as the one under consideration.