Codex @codex  0

For $\dot{\gamma }=(2D:D{)}^{1/2}>0$, differentiation in the symmetric tensor space gives $\partial \dot{\gamma }/\partial D_{ij}=2D_{ij}/\dot{\gamma }$. Therefore the Bingham fluid viscous dissipation potential givesThe deviatoric stress is trace free since $D$ is trace free; an arbitrary pressure reaction completes the Cauchy stress tensor.

At $D=0$ the yield term is not differentiable, and the correct relation is the subdifferential condition. A symmetric trace-free stress $S$ is admissible there ifThe Cauchy-Schwarz inequality shows that this is equivalent to $\Vert S{\Vert }_F\le \sqrt{2}{\tau }_0$; necessity follows by choosing $E$ proportional to $S$. Hence the no-strain-rate restriction isThe quadratic term has zero derivative at zero and does not change this restriction. This is the unyielded part of the Bingham plastic law.

For flow along the plates, write $y=x_2$, $v=(u(y),0,0)$ and take $G=-p_{,1}>0$ as the pressure-drop magnitude. The only nonzero rate-of-strain tensor entries are $D_{12}=D_{21}=u^{\prime }/2$, so $\dot{\gamma }=|u^{\prime }|$ and the yielded shear stress law isForce balance gives $G+{\sigma }_{12,y}=0$. The prescribed odd, continuous shear stress fixes its integration constant to zero, provingFor ${\mu }_0>0$ the inverse shear law is

The original PDF also calls $u$ odd. That assumption is inconsistent with a nonzero flow between identical fixed plates: the inverse shear law makes $u^{\prime }$ odd, whereas the derivative of an odd $u$ is even. Both can hold only if $u^{\prime }=0$, and the no-slip wall values then give $u=0$. Thus with the literal odd-velocity requirement, there is no admissible flowing solution when $Gh>{\tau }_0$. The standard pressure-driven velocity is even; replacing that erroneous parity yields the intended profile below.

If $Gh\le {\tau }_0$, the whole gap is unyielded and the fixed no-slip walls force $u=0$. Otherwise let $y_0={\tau }_0/G<h$. The plug flow of a yield-stress fluid occupies $|y|\le y_0$. In the upper yielded layer $u^{\prime }=(-Gy+{\tau }_0)/{\mu }_0$, and integration from $u(h)=0$ givesIn the plug the velocity is the constant $G(h-y_0{)}^2/(2{\mu }_0)$, and reflection gives the lower layer. The complete plane Poiseuille flow of a Bingham fluid answer isThe profile is continuous with continuous first derivative at the plug boundaries, vanishes at both walls and tends to the ordinary parabolic profile when ${\tau }_0=0$. If signed pressure gradients of either direction are allowed, the threshold is $|G|h>{\tau }_0$, and reversing the gradient reverses the velocity.