A3.16

Transport Processes
Part II, 2001

(i) Incompressible fluid of kinematic viscosity ν\nu occupies a parallel-sided channel 0yh0,<x<0 \leqslant y \leqslant h_{0},-\infty<x<\infty. The wall y=0y=0 is moving parallel to itself, in the xx direction, with velocity Re{Ueiωt}\operatorname{Re}\left\{U e^{i \omega t}\right\}, where tt is time and U,ωU, \omega are real constants. The fluid velocity u(y,t)u(y, t) satisfies the equation

ut=νuyyu_{t}=\nu u_{y y}

write down the boundary conditions satisfied by uu.

Assuming that

u=Re{asinh[b(1η)]eiωt}u=\operatorname{Re}\left\{a \sinh [b(1-\eta)] e^{i \omega t}\right\}

where η=y/h0\eta=y / h_{0}, find the complex constants a,ba, b. Calculate the velocity (in real, not complex, form) in the limit h0(ω/ν)1/20h_{0}(\omega / \nu)^{1 / 2} \rightarrow 0.

(ii) Incompressible fluid of viscosity μ\mu fills the narrow gap between the rigid plane y=0y=0, which moves parallel to itself in the xx-direction with constant speed UU, and the rigid wavy wall y=h(x)y=h(x), which is at rest. The length-scale, LL, over which hh varies is much larger than a typical value, h0h_{0}, of hh.

Assume that inertia is negligible, and therefore that the governing equations for the velocity field (u,v)(u, v) and the pressure pp are

ux+vy=0,px=μ(uxx+uyy),py=μ(vxx+vyy)u_{x}+v_{y}=0, p_{x}=\mu\left(u_{x x}+u_{y y}\right), p_{y}=\mu\left(v_{x x}+v_{y y}\right)

Use scaling arguments to show that these equations reduce approximately to

px=μuyy,py=0p_{x}=\mu u_{y y}, \quad p_{y}=0

Hence calculate the velocity u(x,y)u(x, y), the flow rate

Q=0hudyQ=\int_{0}^{h} u d y

and the viscous shear stress exerted by the fluid on the plane wall,

τ=μuyy=0\tau=-\left.\mu u_{y}\right|_{y=0}

in terms of px,h,Up_{x}, h, U and μ\mu.

Now assume that h=h0(1+ϵsinkx)h=h_{0}(1+\epsilon \sin k x), where ϵ1\epsilon \ll 1 and kh01k h_{0} \ll 1, and that pp is periodic in xx with wavelength 2π/k2 \pi / k. Show that

Q=h0U2(132ϵ2+O(ϵ4))Q=\frac{h_{0} U}{2}\left(1-\frac{3}{2} \epsilon^{2}+O\left(\epsilon^{4}\right)\right)

and calculate τ\tau correct to O(ϵ2)O\left(\epsilon^{2}\right). Does increasing the amplitude ϵ\epsilon of the corrugation cause an increase or a decrease in the force required to move the plane y=0y=0 at the chosen speed U?U ?