  We study a technique for verification of stress and pressure computations on
boundaries in flow simulations. We utilize existing experiments to provide
validation of the simulations. We show that this approach can reveal critical
flaws in simulation algorithms. Using the successful computational algorithms,
we examine Lamb's model for cylinder drag at low Reynolds numbers. We comment
on a discrepancy observed in an experimental paper, suggesting that the domain
size may be a contributing factor. Our simulations on suitably large domains
confirm Lamb's model. We highlight a paradox related to imposing Dirichlet
(Stokes) boundary conditions on polygonal approximations of the curved surface
using finite-element methods that are exactly divergence free. The
finite-element simulations provide very poor representations of drag when the
boundary conditions are imposed strongly. We demonstrate that relaxing the
boundary conditions using Nitsche's method restores high-order approximation.
