advanced29 min studyLesson 178 in Physics

Couette and Poiseuille solutions

Solve fully developed flow between parallel plates, distinguish wall-driven and pressure-driven motion, and check stress and flow-rate balances.

On this page

What you’ll understand

  • Derive the linear Couette and parabolic plane Poiseuille profiles from explicit boundary conditions.
  • Compute flux per unit width and wall stresses for combined pressure and wall driving.
  • Distinguish plane-channel formulas from cylindrical pipe-flow formulas and from developing flow.

Exact solutions expose their assumptions. Consider a long channel between parallel plates, with one plate possibly moving and a pressure difference possibly imposed. Newtonian momentum balance gives a straight profile for pure wall driving and a parabola for pure pressure driving. We will solve their combined problem and calculate transport and wall stresses. Geometry, development length, constitutive behavior, and stability remain part of the model.

Reduce the equations without dropping the driving

Let the plates lie at y = 0 and y = h. The lower plate is stationary and the upper plate moves in the positive x direction at speed U. Assume steady incompressible Newtonian flow with constant viscosity, no streamwise body force, and a fully developed velocity u(y) in the x direction. There is no cross-stream velocity. Under these assumptions the advective acceleration vanishes exactly, because the velocity has no x dependence; this is a geometric cancellation rather than an appeal to small Reynolds number.

The remaining streamwise equation balances the pressure gradient against the second derivative of velocity. Define G = −dp/dx, positive when pressure falls downstream. Integrate twice and impose no-slip at both plates. The result contains a linear wall-driven contribution and a parabolic pressure-driven contribution. Their superposition is valid here because the reduced equation and boundary conditions are linear; it does not imply that arbitrary solutions of the full nonlinear Navier–Stokes equation can be added.

μd2udy2=−G,u(0)=0,u(h)=U,u(y)=Uyh+G2μy(h−y)\mu\frac{d^2u}{dy^2}=-G,\qquad u(0)=0,\quad u(h)=U,\qquad u(y)=\frac{Uy}{h}+\frac{G}{2\mu}y(h-y)
Here h is the full plate separation. G has units of pressure per length. The linear contribution is plane Couette flow; the quadratic contribution is plane Poiseuille flow.↔ Scroll wide equations horizontally.

Integrate transport and differentiate stress

The volumetric flux per unit width is the integral of u across the gap, with units of square meters per second. Multiplying by a physical channel width gives volumetric flow rate only when side-wall effects are negligible. The average speed is this per-width flux divided by h. For pure plane Poiseuille flow, the maximum occurs at the channel midpoint and is one and a half times the average. For pure Couette flow, the average is one-half of the moving-wall speed.

Differentiate the profile to find the shear-stress component tau_xy. It varies linearly with y when G is nonzero, supplying the viscous-force density needed to balance the constant pressure gradient. Remember that traction also depends on surface orientation. On the fluid domain, the outward normals at the upper and lower boundaries point in opposite directions. Consequently the same positive stress component would correspond to opposite streamwise tractions on those two surfaces. This sign check prevents confusing stress components with forces on walls.

q=∫0hu(y) dy=Uh2+Gh312μ,τxy(y)=μUh+G2(h−2y)q=\int_0^h u(y)\,dy=\frac{Uh}{2}+\frac{Gh^3}{12\mu},\qquad \tau_{xy}(y)=\mu\frac{U}{h}+\frac{G}{2}(h-2y)
The symbol q is flux per unit width, not total volumetric flow rate. The stress is a tensor component; boundary traction requires multiplication by the outward normal.↔ Scroll wide equations horizontally.

Keep geometry and validity visible

A circular pipe has a different transverse differential operator and a different profile. For stationary walls and radius R, cylindrical Poiseuille flow has u(r) = G(R² − r²)/(4mu), total flow rate pi G R⁴/(8mu), and maximum speed twice the cross-sectional mean. Do not insert a pipe diameter into the plane-channel formula. The numerical factors encode geometry and boundary conditions. In a finite-width rectangular duct, both transverse directions matter unless the width is large compared with the plate separation.

Fully developed means the streamwise velocity profile has stopped changing along the channel, not that the fluid stopped moving or that pressure stopped decreasing. Near an inlet, acceleration and development invalidate the simple reduction. At sufficiently strong disturbances the assumed laminar state may also be inappropriate, even though its algebra remains an exact solution of the ideal equations. Temperature-dependent viscosity, non-Newtonian rheology, slip, and significant edge effects require revisions. Use the exact profile as a benchmark for a stated regime, then measure where the real system departs.

Combine moving-wall transport with a pressure gradient

An ideal plane channel has full gap h = 0.01 m, upper-wall speed U = 0.2 m/s, dynamic viscosity mu = 0.05 Pa s, and downstream pressure-gradient magnitude G = 100 Pa/m. Calculate the midpoint speed, the flux per unit width, and the two boundary values of tau_xy.

  1. At y = h/2 = 0.005 m, the wall-driven contribution is U/2 = 0.1 m/s. The pressure contribution is G h²/(8mu) = 100 × 0.0001 / 0.4 = 0.025 m/s. Their sum is 0.125 m/s.

  2. The wall-driven flux is Uh/2 = 0.001 m²/s. The pressure-driven flux is Gh³/(12mu) = 0.0001/0.6 = 0.0001666667 m²/s. The combined flux is approximately 0.0011666667 m²/s, giving mean speed approximately 0.1166667 m/s.

  3. The uniform Couette stress contribution is mu U/h = 1 Pa. The pressure contribution is plus 0.5 Pa at the lower wall and minus 0.5 Pa at the upper wall. Thus tau_xy(0) = 1.5 Pa and tau_xy(h) = 0.5 Pa; their difference is −1 Pa.

The result

The midpoint speed is 0.125 m/s and the per-width flux is approximately 0.00116667 m²/s. Boundary stress components are 1.5 Pa and 0.5 Pa. For a unit-width, unit-length control volume, pressure supplies 1 N downstream and the wall tractions supply a net −1 N.

Explore the model

Interactive lab

Interactive lab · function

Inspect the curvature of a plane Poiseuille profile

Use the centered dimensionless coordinate x = 2y/h − 1, so the channel walls are x = −1 and x = 1. For pure plane Poiseuille flow, normalized velocity is 1 − x². Set Coefficient a to −1: the lab shows the velocity relative to its centerline value, namely −x². Predict a zero slope at the center and slopes of opposite sign at the walls. Change Position x within the wall interval to check.

Inspect the curvature of a plane Poiseuille profileThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forAdd the missing constant 1 mentally to recover normalized velocity. The lab plots only a quadratic, its tangent, and its signed integral; it does not add a Couette term, solve fluid equations, or determine the dimensional pressure gradient. Keep the interpretation inside −1 ≤ x ≤ 1.

Check your understanding

Two stationary parallel plates support steady fully developed Newtonian pressure-driven flow with no slip. For this plane-channel geometry, how does the maximum streamwise speed compare with the mean speed obtained by averaging over the plate separation?

Try it yourself

Take stationary plates separated by 0.02 m, with viscosity 0.1 Pa s and pressure falling at 40 Pa/m. Compute the centerline speed and average speed for fully developed plane Poiseuille flow. Check that the lower and upper wall stress components have equal magnitudes and opposite signs.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Set the wall speed U to zero in the combined formulas. Use the full gap h, not the half-gap, and distinguish stress components from tractions on oppositely oriented surfaces.

Compare with the solution

The centerline speed is G h²/(8mu) = 40 × 0.0004 / 0.8 = 0.02 m/s. The mean speed is G h²/(12mu) = 0.016/1.2, approximately 0.0133333 m/s. Their ratio is 1.5 as expected. At y = 0 the stress component is G h/2 = 0.4 Pa; at y = h it is −0.4 Pa. The corresponding wall tractions on the fluid both oppose the downstream motion because the outward normals have opposite signs.

Take it with you

  • Plane Couette and Poiseuille profiles follow from the same reduced momentum balance with different driving.
  • Flow rate is an integral of velocity, while shear stress is a derivative; both provide independent checks on a solution.
  • Plane channels and circular pipes have different geometric factors, and fully developed flow can still have a pressure gradient.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.