Dear all,
I am trying to understand the lubrication correction
in Stokesian dynamics implemented in current developing version of ESPResSo.
According to David Schwoerer's thesis, the function "sd_compute_resistance_matrix_sparse()"
computes the lubrication correction described in equation (3.24) as R^{lc} = R^{lub} - R^{2b,ff}.
In addition, there is one comment in the code
referring R^{lub} to 'N.-Q. Nguyen and A. J. C. Ladd, PHYSICAL REVIEW E 66, 046708 (2002) equation (34)'.
But I still do not understand this function quite well, especially the terms containing the variable "ls = log(s) = log(|r|/a - 2)",
which I cannot find neither in section 3.1.2 in the thesis nor in Ladd's paper.
Would anyone like to give me more references about how ESPResSo calculates this correction?
Many thanks in advance
Lei