  In this paper, we establish the estimates for the gradient and the
second-order partial derivatives for the Stokes flow in the presence of two
closely located strictly convex inclusions in dimension three. Moreover, the
blow-up rate of the gradient is showed to be optimal by a pointwise upper bound
and a lower bound in the narrowest region. We also show the optimal blow-up
rate of Cauchy stress tensor. In dimensions greater than three, the upper
bounds of the gradient are established. These results answer the questions
raised in [25].
