 Solution

Let $t(n)$ denote traction: the force per unit area exerted across an oriented internal surface with unit normal $n$. Define ${\sigma }_{ij}$ as component $i$ of the traction on a coordinate plane whose normal is $e_j$. Balance of linear momentum on a shrinking tetrahedron makes the surface-force terms, of order area, dominate its inertia and body forces, of order volume. Resolving its oblique face into the three coordinate projected areas gives Cauchy's traction formulaOpposite faces have opposite tractions, as follows independently from a shrinking pillbox. This fixes both the index and normal conventions.

For a material volume $V(t)$, the Reynolds transport theorem says $d{\int }_{V(t)}fdV/dt={\int }_{V(t)}[f_t+\nabla \cdot (fu)]dV$. Apply it componentwise to $f=\rho u_i$ and use mass conservation. With no body force in this calculation, balance of linear momentum and the divergence theorem giveArbitrariness of the material volume yields $\boxed{\rho D{u}_i/Dt={\partial }_j{\sigma }_{ij}}$, with $D/Dt={\partial }_t+u_j{\partial }_j$.

For constant dynamic viscosity and incompressible velocity, ${\partial }_j{u}_j=0$, soThus $\boxed{\rho (u_t+u\cdot \nabla u)=-\nabla p+\mu {\nabla }^{2}u}$.

In fully developed steady pipe flow $u=w(x,y)e_z$, the convective acceleration vanishes. The transverse equations give $p_x=p_y=0$ and the axial equation gives $p_z=\mu (w_{xx}+w_{yy})$. One side depends only on $z$ and the other only on $(x,y)$, so both equal a constant. HenceIf $n$ is outward from the fluid cross-section, the wall-on-fluid axial shear per length is $\mu \oint {\partial }_{n}wds=\mu {\int }_A\Delta wdA=A{p}_z$. The fluid-on-pipe force is its negative:For forward pressure-driven flow $p_z<0$, this force points along the flow. Its value depends on cross-sectional area and pressure drop, even though the detailed velocity profile depends on the shape.