Skip to main content
Sandbox Physics

M096 · Fluid mechanics

Couette–Poiseuille Viscous Flow

Set a plate speed, pressure gradient and gap. Move a probe through the velocity and stress profiles, then infer viscosity and pressure drive from acquired readings.

Interactive modelCouette–Poiseuille Viscous Flow
Probe velocityReady\text{Ready}
Probe signed shearReady\text{Ready}
Signed flow per widthReady\text{Ready}
Upper wall power per footprintReady\text{Ready}
Pressure power per footprintReady\text{Ready}
Viscous heating per footprintReady\text{Ready}
Peak-speed Reynolds numberReady\text{Ready}
Power closure residualReady\text{Ready}
Viscosity inferred from raw readingsReady\text{Ready}
Pressure drive inferred from raw readingsReady\text{Ready}
Acquired stationsReady\text{Ready}
Raw-coordinate fit residualReady\text{Ready}
Inference statusReady\text{Ready}
Model domainReady\text{Ready}

Physics tutorial

Velocity, stress and dissipation

BackgroundTwo independent measurement benches share one acquisition workflow; each keeps its own physical law, geometry and validity screen.

Why it mattersConnect physical geometry, reference profiles and a raw sensor record before interpreting an estimate.

Start with the essentials

Focus question
Can shear reveal a hidden viscosity?
One-sentence intuition
Infer parameters from independent raw channels and check the approximation domain first.

Core mathematical model

Signed velocity and shear

η=y/H,G=−dp/dxu=Uη+GH22μη(1−η)τ=μU/H+GH(1/2−η)\begin{aligned}\eta&=y/H,\quad G=-dp/dx\\u&=U\eta+\frac{GH^2}{2\mu}\eta(1-\eta)\\\tau&=\mu U/H+GH(1/2-\eta)\end{aligned}

The lower wall is fixed; shear is the stress on an upward-facing fluid plane.

Flux per unit width

q=∫0Hu dy=UH2+GH312μq=\int_0^H u\,dy=\frac{UH}{2}+\frac{GH^3}{12\mu}

Positive and negative layers can cancel even while the fluid dissipates energy.

Power closure per footprint

Gq+Uτ(H)=DD=∫0Hμ(uy)2 dy\begin{aligned}Gq+U\tau(H)&=\mathcal D\\\mathcal D&=\int_0^H\mu(u_y)^2\,dy\end{aligned}

The wall can extract mechanical energy; pressure plus wall work equals heating.

Infer from acquired profiles

u^=c+aη+bη(1−η)μ^=∑iu^y,iτi∑iu^y,i2G^=2μ^b/H2\begin{aligned}\widehat u&=c+a\eta+b\eta(1-\eta)\\\widehat\mu&=\frac{\sum_i\widehat u_{y,i}\tau_i}{\sum_i\widehat u_{y,i}^2}\\\widehat G&=2\widehat\mu b/H^2\end{aligned}

A free velocity offset is fitted. The stress fit assumes zero sensor bias; the derivative must be nonzero.

Domain diagnostic

Repeak=ρHmax⁡y∣u∣/μ\mathrm{Re}_{\mathrm{peak}}=\rho H\max_y|u|/\mu

A teaching threshold of 200 blocks acquisition outside the selected low-inertia domain, without proving stability.

Common difficulties

Reference versus inference

Typical misconceptionThe estimate can be read without recording data.

Better mental modelClear the record and the estimate disappears. Sensor perturbations can produce nonphysical fits; estimates are not clipped or independently calibrated.

Run the experiment

  1. 01

    Separate the two drives

    Compare Wall only, Pressure only and Combined drive.

    What to observe: Linear and parabolic contributions add; both meet no-slip wall values.
  2. 02

    Find counterflow

    Select Counterflow and move the gold probe across the gap.

    What to observe: A moving upper wall can coexist with negative interior velocity.
  3. 03

    Close the energy ledger

    Change wall speed and inspect signed wall power.

    What to observe: Pressure and wall work sum to nonnegative viscous heating.
  4. 04

    Acquire and infer

    Acquire the gap profile, increase perturbations and export CSV.

    What to observe: The fit uses recorded velocity and shear, without reference viscosity or gradient.