Skip to main content

Couette flow with a pressure gradient

Short answer

Velocity u = (Uy/h) - (1/(2μ))(dp/dx)(hy - y^2); q = Uh/2 - (h^3/(12μ))(dp/dx). Pressure gradient alters profile: zero gives linear, favorable bulges outward, adverse may reverse flow. Shear stress on moving plate τ_h = μU/h + (h/2)(dp/dx): 75 Pa for dp/dx=0, 55 Pa for dp/dx=-10 kPa/m.

(a) Derive the velocity distribution for steady laminar flow between two parallel plates when the upper plate moves with velocity UU and the lower plate is fixed, with a constant pressure gradient dp/dxdp/dx along the flow. Obtain the discharge per unit width. (b) Discuss the effect of the pressure gradient on the velocity profile. For U=3 m/sU = 3\ \text{m/s}, gap h=4 mmh = 4\ \text{mm} and oil of viscosity 0.1 Pa s, find the shear stress on the moving plate when dp/dx=0dp/dx = 0 and when dp/dx=−10 kPa/mdp/dx = -10\ \text{kPa/m}.

(a) Velocity distribution

For steady, fully developed laminar flow, the x-momentum (Navier-Stokes) equation reduces to

μd2udy2=dpdx\mu\frac{d^2u}{dy^2} = \frac{dp}{dx}

Integrating twice: u=12μdpdxy2+C1y+C2u = \dfrac{1}{2\mu}\dfrac{dp}{dx}y^2 + C_1y + C_2.

Boundary conditions: u=0u = 0 at y=0y = 0, and u=Uu = U at y=hy = h. Then C2=0C_2 = 0 and C1=Uh−h2μdpdxC_1 = \dfrac{U}{h} - \dfrac{h}{2\mu}\dfrac{dp}{dx}.

u=Uyh−12μdpdx(hy−y2)u = \frac{U y}{h} - \frac{1}{2\mu}\frac{dp}{dx}\left(hy - y^2\right)

The first term is the Couette (linear) part due to the moving plate, and the second is the Poiseuille (parabolic) part due to the pressure gradient.

Discharge per unit width:

q=∫0hu dy=Uh2−h312μdpdxq = \int_0^h u\,dy = \frac{Uh}{2} - \frac{h^3}{12\mu}\frac{dp}{dx}

(b) Effect of pressure gradient

Let P=−h22μUdpdxP = -\dfrac{h^2}{2\mu U}\dfrac{dp}{dx}.

Casedp/dxdp/dxProfile
Simple Couette flow00Straight line from 0 to UU
Favourable gradient<0< 0 (P>0P > 0)Profile bulges outward, velocity above the linear value
Adverse gradient>0> 0 (P<0P < 0)Profile bends inward; if dp/dx>2μU/h2dp/dx > 2\mu U/h^2 the velocity near the fixed plate reverses (back flow)
  dp/dx = 0     dp/dx < 0     dp/dx > 0 (large)
      /|           _.|             |
     / |         .'  |            _|
    /  |        /    |          <-'/
   /___|       /_____|         /___|

Shear stress: τ=μdudy=μUh−(h2−y)dpdx\tau = \mu\dfrac{du}{dy} = \mu\dfrac{U}{h} - \left(\dfrac{h}{2} - y\right)\dfrac{dp}{dx}

At the moving plate (y=hy = h): τh=μUh+h2dpdx\tau_h = \dfrac{\mu U}{h} + \dfrac{h}{2}\dfrac{dp}{dx}

Numerical

μU/h=0.1×3/0.004=75 Pa\mu U/h = 0.1 \times 3/0.004 = 75\ \text{Pa}.

  • For dp/dx=0dp/dx = 0: τh=75 Pa\tau_h = 75\ \text{Pa}.
  • For dp/dx=−10 000 Pa/mdp/dx = -10\,000\ \text{Pa/m}:
τh=75+0.0042(−10 000)=75−20=55 Pa\tau_h = 75 + \frac{0.004}{2}(-10\,000) = 75 - 20 = 55\ \text{Pa}

(The shear stress on the fixed plate would be 75+20=95 Pa75 + 20 = 95\ \text{Pa}.) The discharge per unit width rises from Uh/2=0.006 m2/sUh/2 = 0.006\ \text{m}^2/\text{s} to 0.006+(0.004)3×10 00012×0.1=0.00653 m2/s0.006 + \dfrac{(0.004)^3 \times 10\,000}{12 \times 0.1} = 0.00653\ \text{m}^2/\text{s}.

Answer: τh=75 Pa\tau_h = 75\ \text{Pa} for dp/dx=0dp/dx = 0, and 55 Pa55\ \text{Pa} for dp/dx=−10 kPa/mdp/dx = -10\ \text{kPa/m}

✓ Answer checked by Claude Opus 5.5 (high effort)

All Fluid Mechanics chapter 6 questions →