desirable for the credibility of our approach) it is then natural to specify $\epsilon$ in the form $\sqrt{\lambda}R_e^{\beta}$ and to demand that $\beta \ge 1$. In particular, the choice $\beta=1$ arranges that $$\Gamma_{HA} =  O(R_e^2/\log R_e)\ ,\ \ \Gamma_{NL} = O(R_e/\log^2{R_e})$$ and we have established that $\Gamma = o(R_e)$. It also decays exponentially at infinity.

The drag on the cylinder
========================

The calculation of the drag can be done by integrating the pointwise force on the cylinder over its surface. The pointwise force has two components. One involves the local rate of strain in the fluid, and the other is the pressure force. The first requires a purely local calculation, but the second requires an integration of the pressure equation from infinity to the cylinder. A question is how this second part can be carried out [*without*]{} any [*global*]{} representation of the flow field. For the [*specific*]{} case of low Reynolds number calculations with a certain symmetry there are ways around the problem that we shall exploit presently. To make these matters clear, we shall summarize a first-principles calculation of the force using the results for cylindrical polar coordinates for the rate of strain as given by Appendix A of 