Generating, from scratch, the near-field asymptotic forms of scalar resistance functions for two unequal rigid spheres in low Reynolds number flow

Abstract

The motion of rigid spherical particles suspended in a low Reynolds number fluid can be related to the forces, torques, and stresslets acting upon them by 22 scalar resistance functions, commonly notated X11A, X12A, Y11A, etc. Near-field asymptotic forms of these resistance functions were derived in Jeffrey and Onishi [“Calculation of the resistance and mobility functions for two unequal rigid spheres in low-Reynolds-number flow,” J. Fluid Mech. 139, 261–290 (1984)] and Jeffrey [“The calculation of the low Reynolds number resistance functions for two unequal spheres,” Phys. Fluids A 4, 16–29 (1992)]; these forms are now used in several numerical methods for suspension mechanics. However, the first of these important papers contains a number of small errors, which make it difficult for the reader to correctly evaluate the functions for parameters not explicitly tabulated. This short article comprehensively corrects these errors and adds formulas that were originally omitted from both papers, so that the reader can verify and implement the equations independently. The corrected expressions, rationalized and using contemporary nondimensionalization, are shown to match mid-field values of these scalars, which are calculated through an alternative method. A Python script to generate and evaluate these functions is provided.