Top Layer

The lubrication forces for $N$ (all) particles The lubrication forces for $N$ (all) particles

$$\bm{F} = \bm{L} \cdot \bm{U} ,$$ (6.1)

is calculated by

$${\tt calc\_lub\_\langle d\rangle\langle ver\rangle\ (sys, u, f)} ,$$ (6.2)

which is called by the top-level function

$${\tt solve\_\langle type\rangle\_lub\_\langle d\rangle\langle ver\rangle\ ()} ,$$ (6.3)

for mobility, resistance, and mixed problems (type) in 2- and 3-dimensions (d) with F, FT, and FTS versions (ver).

The function ${\tt calc\_lub\_\langle d\rangle\langle ver\rangle\ ()}$ has particle-loops. Under the periodic boundary condition, particles in the image cells are also taken into account.

In the loop for particle-particle pairs, either

$${\tt calc\_lub\_\langle ver\rangle\_2b()}$$ (6.4)

for monodisperse system or

$${\tt calc\_lub\_\langle ver\rangle\_2b\_poly()}$$ (6.5)

for polydisperse system is called.