= Solution
Let $a=\lambda\dot\gamma$. For simple shear, the steady conformation equation is
$$
F\mathbf C-\mathbf I
=\lambda[(\nabla\mathbf u)\mathbf C
+\mathbf C(\nabla\mathbf u)^T].
$$
Its nonzero components are
$$
C_{yy}=C_{zz}=\frac1F,
\qquad
C_{xy}=\frac a{F^2},
\qquad
C_{xx}=\frac1F+\frac{2a^2}{F^3}.
$$
Since $\boldsymbol\tau=\eta_s\dot{\boldsymbol\gamma}+(\eta_p/\lambda)(F\mathbf C-\mathbf I)$,
$$
\boxed{\tau_{xy}=\left(\eta_s+\frac{\eta_p}{F}\right)\dot\gamma},
$$
$$
\boxed{\tau_{xx}=\frac{2\eta_p\lambda\dot\gamma^2}{F^2}},
\qquad
\tau_{yy}=\tau_{zz}=0.
$$
Thus the <FENE-P model> has
$$
\boxed{N_1=\frac{2\eta_p\lambda\dot\gamma^2}{F^2}>0,
\qquad N_2=0}.
$$
Back to article page
