free_solution = lambda kB, T, eta, Rh: kB * T / (6 * _np.pi * eta * Rh)