Skip to content

Constant Torque Through Every Shell

Irodov 1.333DifficultyDifficulty 3 of 5UniquenessUniqueness 3 of 5

Problem Irodov 1.333

Viscous fluid of viscosity η\eta fills the annular gap between two long, coaxial cylinders of radii R1R_1 and R2R_2, with R1<R2R_1<R_2. The inner cylinder (R1R_1) is held fixed while the outer one (R2R_2) is driven at a constant angular velocity ω2\omega_2, and the resulting flow of the fluid between them is laminar. The friction force per unit area on any cylindrical fluid surface of radius rr is given as σ=ηr (∂ω/∂r)\sigma=\eta r\,(\partial\omega/\partial r), where ω(r)\omega(r) is the fluid’s own local angular velocity. Using this, find (a) how the fluid’s angular velocity ω\omega depends on rr for R1≤r≤R2R_1\le r\le R_2, and (b) the moment, per unit length of the cylinders, of the friction forces acting on the outer cylinder.

Top view down the shared axis. The inner cylinder (radius R_1, ticked rim) is held; the outer (radius R_2) is driven at ω_2. The dashed circle is the variable shell at radius r, carrying tangential velocity v(r).

Top view down the shared axis. The inner cylinder (radius R1R_1, ticked rim) is held; the outer (radius R2R_2) is driven at ω2\omega_2. The dashed circle is the variable shell at radius rr, carrying tangential velocity v(r)v(r).

Given
  • η\eta
  • R1R_1, R2R_2
  • ω2\omega_2
  • σ=ηr ∂ω/∂r\sigma = \eta r\,\partial\omega/\partial r
Find
(a) ω(r)\omega(r) for R1≤r≤R2R_1\le r\le R_2; (b) NN, the moment of friction on the outer cylinder per unit length

Solution

Step 1The torque crossing every shell is the same constant

The torque per unit length crossing a cylindrical shell of radius rr is the shear stress σ(r)\sigma(r) times the shell’s area per unit length (2πr2\pi r) times the moment arm (rr again). Setting this equal to a constant MM and substituting the given σ\sigma,

M=σ(r)⋅2πr⋅r=2πηr3dωdr.(1)M = \sigma(r)\cdot 2\pi r \cdot r = 2\pi\eta r^3\frac{d\omega}{dr}. \tag{1}
Step 2Integrate the constant-torque condition

Dividing through by 2πηr32\pi\eta r^3 and integrating once gives ω\omega up to two constants.

dωdr=M2πη 1r3ω(r)=−M4πη1r2+C.(2) \begin{aligned} \frac{d\omega}{dr} &= \frac{M}{2\pi\eta}\,\frac{1}{r^3} \\ \omega(r) &= -\frac{M}{4\pi\eta}\frac{1}{r^2} + C. \end{aligned} \tag{2}
Step 3Apply the two boundary conditions

The inner cylinder is held, ω(R1)=0\omega(R_1)=0; the fluid at the outer wall shares that cylinder’s own speed by no-slip, ω(R2)=ω2\omega(R_2)=\omega_2.

−M4πη1R12+C=0−M4πη1R22+C=ω2\begin{align*}-\frac{M}{4\pi\eta}\frac{1}{R_1^2} + C &= 0 \\ -\frac{M}{4\pi\eta}\frac{1}{R_2^2} + C &= \omega_2\end{align*}

Subtracting the first equation from the second eliminates CC and isolates MM.

M4πη(1R12−1R22)=ω2\frac{M}{4\pi\eta}\left(\frac{1}{R_1^2}-\frac{1}{R_2^2}\right) = \omega_2
M=4πηω2R12R22R22−R12M = 4\pi\eta\omega_2\dfrac{R_1^2R_2^2}{R_2^2-R_1^2}
(3)

The first boundary equation then gives the remaining constant,

C=M4πηR12=ω2R22R22−R12.C = \frac{M}{4\pi\eta R_1^2} = \frac{\omega_2 R_2^2}{R_2^2-R_1^2}.
Step 4Assemble ω(r)\omega(r)

Substituting MM from (3) and CC back into the general solution (2),

ω(r)=−14πηr2(4πηω2R12R22R22−R12)+ω2R22R22−R12=ω2R22R22−R12(1−R12r2).\begin{align*}\omega(r) &= -\frac{1}{4\pi\eta r^2}\left(4\pi\eta\omega_2\frac{R_1^2R_2^2}{R_2^2-R_1^2}\right) + \frac{\omega_2 R_2^2}{R_2^2-R_1^2} \\ &= \omega_2\frac{R_2^2}{R_2^2-R_1^2}\left(1-\frac{R_1^2}{r^2}\right).\end{align*}

Taking R12R_1^2 out of the bracket gives the profile in its symmetric form.

ω(r)=ω2R12R22R22−R12(1R12−1r2)\omega(r) = \omega_2\dfrac{R_1^2R_2^2}{R_2^2-R_1^2}\left(\dfrac{1}{R_1^2}-\dfrac{1}{r^2}\right)
(4)
Step 5The moment on the outer cylinder

Because the torque transmitted across every shell is the same constant MM of (1), the moment the fluid exerts on the outer cylinder, at r=R2r=R_2, is that very same MM: no further integration is needed.

N=4πηω2R12R22R22−R12N = 4\pi\eta\omega_2\dfrac{R_1^2R_2^2}{R_2^2-R_1^2}
(5)
Check

As R1→0R_1\to0 (the inner cylinder shrinks to the axis), the prefactor R12R22/(R22−R12)→R12R_1^2R_2^2/(R_2^2-R_1^2)\to R_1^2, so ω(r)→ω2(1−R12/r2)→ω2\omega(r)\to\omega_2\bigl(1-R_1^2/r^2\bigr)\to\omega_2 at every fixed rr: with no inner boundary to hold it back, the fluid simply co-rotates rigidly with the outer wall, and the shear ∂ω/∂r→0\partial\omega/\partial r\to0 everywhere, consistent with N→4πηω2R12→0N\to4\pi\eta\omega_2 R_1^2\to0. In the opposite, thin-gap limit R1→R2≡RR_1\to R_2\equiv R with gap d=R2−R1≪Rd=R_2-R_1\ll R, N→4πηω2R4/(2Rd)=2πηω2R3/dN\to4\pi\eta\omega_2R^4/(2Rd)=2\pi\eta\omega_2R^3/d, exactly the flat-plate estimate (the shear ηω2R/d\eta\omega_2R/d acting over the area 2πR2\pi R per unit length, and that force then carried on the moment arm RR) for a narrow Couette gap, confirming the curved result reduces correctly to the flat one.

Remark

The whole problem turns on one fact: with no torque source hiding anywhere in the fluid’s interior, the viscous torque crosses every concentric shell unchanged in steady state, even though the stress carrying it thins out as 1/r21/r^2. That single conserved quantity is what lets (1) set r3 dω/dr=constr^3\,d\omega/dr=\text{const} directly from the given σ\sigma, without ever writing down the fluid’s equation of motion; the same trick, applied to r dv/drr\,dv/dr, is exactly what makes 1.332’s straight coaxial sliding problem solvable by inspection too. In both problems the quantity that crosses every shell unchanged is whatever the steady balance conserves, the axial force in that straight film and the torque here, and in both the stress itself still falls off with radius.