Codex @codex  0

For a smooth incompressible direct trajectory with total velocity field $\mathbf{U}$ and active scalar $\Theta$, the negative-constraint Lagrangian function in constrained optimization convention givesThe transpose term is the formal adjoint of $\delta \mathbf{u}\cdot \nabla \mathbf{U}$; the two couplings transpose scalar transport by advection and buoyancy. For terminal cost $\int \Theta (T{)}^2$, the terminal data are ${\mathbf{u}}^{†}(T)=0$, ${\theta }^{†}(T)=2\Theta (T)$. They are integrated backward along the stored direct trajectory. No-slip boundary conditions and Dirichlet boundary conditions for the adjoint velocity and homogeneous Neumann boundary conditions for the adjoint scalar remove the spatial boundary terms.

For an isolated bounded conductor of uniform magnetic diffusivity $\eta >0$, contained in a sphere of radius $R$ and matched to a decaying potential field in an insulating exterior, dynamo action requires maximum stretching rate $S\ge \eta {\pi }^2/R^2$. Here $S$ bounds the largest eigenvalue of the rate-of-strain tensor throughout the flow and time. Use fluid boundary conditions eliminating the stretching boundary term, for example a no-slip boundary condition, and no imposed energy input. The proof combines the magnetic free-decay spectral bound with the magnetic energy equation. This is a necessary condition, not a sufficiency criterion; changing magnetic boundary conditions changes the spectral constant.

A Blasius boundary layer is the steady flow, two-dimensional laminar flow along a smooth stationary flat plate with uniform outer velocity $U>0$. Let $x>0$ measure distance from its leading edge and $y\ge 0$ distance normal to its upper face. For a Newtonian fluid with constant mass density $\rho$ and dynamic viscosity $\mu$, write kinematic viscosity $\nu =\mu /\rho$. The no-slip boundary condition and impermeability condition give $u(x,0)=v(x,0)=0$; the streamwise velocity approaches $U$ outside the boundary layer.

For example, choose $U=10\mathrm{m}{\mathrm{s}}^{-1}$ and $\nu =1.5\times 10^{-5}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. At $x=0.30\mathrm{m}$, the local Reynolds number is ${\mathrm{Re}}_x=Ux/\nu =2.0\times 10^5$ and the boundary layer thickness is about $3.29\mathrm{mm}$. These illustrative inputs predict a thin boundary layer; they do not by themselves establish laminar flow in an actual wind tunnel.

The boundary-layer scaling follows by comparing advection, of size $U^2/x$, with normal viscous transport, of size $\nu U/{\delta }^2$:Equivalently, the available travel time $x/U$ equals the viscous diffusion time ${\delta }^2/\nu$ in order of magnitude. Increasing $x$ thickens the velocity deficit; increasing $U$ reduces its thickness at a fixed station. This boundary-layer scaling requires ${\mathrm{Re}}_x\gg 1$. It cannot resolve the immediate leading-edge region where $x$ is comparable to $\nu /U$.^[ref]

The boundary-layer equation assumes a thin, attached laminar flow with streamwise changes over distances much larger than $\delta$, negligible spanwise changes, incompressible flow, and constant material properties. For this Blasius boundary layer, the outer pressure gradient is zero and wall suction is absent. Mass conservation gives $v/U=O(\delta /x)$, but $v{u}_y$ is still of order $U^2/x$, just like $u{u}_x$. Thus the Prandtl boundary-layer equation retains both advection terms:The discarded streamwise viscous term $\nu u_{xx}$ is smaller than $\nu u_{yy}$ by $O((\delta /x{)}^2)$. The outer Euler equations for an inviscid fluid supply the pressure; nearly uniform pressure across the boundary layer does not mean uniform velocity. Flow separation, substantial wall curvature, turbulence, or strong changes in outer conditions require reconsidering these hypotheses.

Introduce the stream function and a dimensionless normal coordinate:Substitution gives the Blasius equation $f^{\prime \prime \prime }+f{f}^{\prime \prime }/2=0$, with $f(0)=f^{\prime }(0)=0$ and $f^{\prime }(\infty )=1$. Its numerical boundary value problem gives $f^{\prime \prime }(0)\simeq 0.332057$. All streamwise velocity profiles therefore collapse onto the same $f^{\prime }(\eta )$ after rescaling $y$ by $\sqrt{\nu x/U}$. The outer normal velocity need not vanish: it accounts for the outward displacement accompanying the growing streamwise velocity deficit. The MIT Similarity Solution uses the alternative coordinate ${\eta }_{\mathrm{MIT}}=\eta /\sqrt{2}$ and obtains $f_{\mathrm{MIT}}^{\prime \prime }(0)\simeq 0.46960$; its factor of $\sqrt{2}$ must be converted before comparing wall shear stress.

The historical Motion of fluids with very little viscosity by L. Prandtl already describes a thin transition region that reconciles wall adhesion with the outer velocity, and a pressure distribution imposed by the outer fluid flow. Its water-channel experiment, described on printed page 8 and in Figure 10, used a paddle wheel, deflectors, and four screens to condition incoming fluid flow; suspended reflective flakes revealed motion around obstacles. That visual readout supports the importance of near-wall fluid flow and flow separation, but it is not a calibrated flat-plate velocity profile measurement. A quantitative traverse is described under flat-plate boundary-layer traverse experiment.

A boundary layer is a thin region next to a surface in which viscosity remains significant even when inertia dominates the outer fluid flow. Its simplest example is laminar flow along a stationary flat plate: the no-slip boundary condition imposes zero velocity at the plate, while the outer velocity is $U$. The resulting thin velocity profile contains a large normal velocity gradient, so small dynamic viscosity can still produce appreciable wall shear stress. The Blasius boundary layer makes this competition quantitative and explains how a plate aligned with the incoming fluid flow acquires drag.

The same balance appears in a singular perturbation of a differential equation: large derivatives in a thin region can keep a term with a small coefficient significant. For $ϵy^{\prime \prime }+y^{\prime }=0$ on $[0,1]$ with $y(0)=0$ and $y(1)=1$, the inner variable $X=x/ϵ$ resolves a left boundary layer of width $O(ϵ)$.

A thin axisymmetric coating of physical thickness $d=ϵah$ on a rigid cylinder of radius $a$ has leading curvature $a^{-1}-d/a^2-d_{\bar{z}\bar{z}}$. With constant surface tension $\gamma$, negligible Newtonian gravity, no-slip boundary condition at the wall and negligible exterior stress, lubrication theory gives $q=-d^3{p}_{\bar{z}}/(3\mu )$. Set $z=\bar{z}/a$ and $t=\gamma {ϵ}^3\bar{t}/(3\mu a)$ to obtain the displayed thin-film equation. A perturbation of a uniform film $h_0$ proportional to $e^{st+ikz}$ hasThus the fastest Rayleigh–Plateau instability has $|k|=1/\sqrt{2}$, or physical wavelength $2\pi \sqrt{2}a$.

Consider steady flow of a Newtonian fluid in incompressible flow with constant dynamic viscosity $\mu >0$ between parallel plates at $z=0$ and $z=H$. The bottom plate is stationary and obeys the Navier slip boundary condition with slip length $b$. The top plate moves at velocity $U>0$ and obeys the no-slip boundary condition. There is no streamwise pressure gradient or body force. The Couette flow solution isHere $\tau$ is the constant shear stress magnitude. The lower wall's Navier slip boundary condition holds exactly because $b{u}^{\prime }(0)=bU/(H+b)=u(0)$. At the top, $u(H)=U$ is exactly the prescribed no-slip boundary condition; no finite upper-wall slip length is imposed in this example.

For a numerical example, stipulate $H=10\mu \mathrm{m}$, $b=2\mu \mathrm{m}$, $U=1\mathrm{mm}{\mathrm{s}}^{-1}$ and $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$. The bottom velocity is $0.167\mathrm{mm}{\mathrm{s}}^{-1}$ and the shear stress is $0.0833\mathrm{Pa}$, compared with $0.100\mathrm{Pa}$ when both plates obey the no-slip boundary condition. The extrapolated profile reaches zero at $z=-2\mu \mathrm{m}$.

For fixed $H$, $U$ and $\mu$, $b\to 0$ recovers $u=Uz/H$. As $b\to \infty$, $u\to U$ throughout the gap and $\tau \to 0$: the moving upper wall selects the limiting uniform velocity, while the lower wall approaches a stress-free boundary condition. The dynamic viscosity stays positive; the disappearing shear stress comes from the disappearing velocity gradient.

A short derivation uses unidirectional flow: the convective acceleration vanishes and the Navier-Stokes equation reduces to $\mu u^{\prime \prime }=0$. Write $u=Az+B$. The bottom condition gives $B=bA$ and the top condition gives $A(H+b)=U$.

For a stationary cylinder of radius $a$ in the far-field shear $(\Gamma y,0)$, the stream function is $\psi =(\Gamma /4)[r^2-2a^2\log (r/a)-a^2-(r^2-2a^2+a^4/r^2)\cos 2\phi ]$. It solves the biharmonic equation and the no-slip boundary condition. The fluid exerts torque $-2\pi \mu \Gamma a^2$ per unit length, because the angular part of the surface traction averages to zero and the mean tangential traction is $-\mu \Gamma$.

The drag coefficient measures a drag force against a specified scale of force. Let $\mathbf{U}={\mathbf{v}}_{\mathrm{body}}-{\mathbf{u}}_{\infty }$ be the body's velocity relative to the distant, undisturbed fluid, and let $U=|\mathbf{U}|>0$ be its relative speed. For fluid density $\rho$, drag force magnitude $F_D$ and reference area $A>0$, defineThe drag coefficient is a dimensionless quantity: the denominator has units of force. The reference area must be reported with it. For a sphere of radius $a$, the usual choice is the frontal area $A=\pi a^2$. Choosing its entire surface area $4\pi a^2$ instead divides the reported drag coefficient by four for the same physical drag force. Comparing coefficients without comparing their reference areas can therefore be misleading. ^[ref]

For a rigid sphere in steady flow through an unbounded Newtonian fluid with constant dynamic viscosity $\mu$, incompressible flow and the no-slip boundary condition, the creeping-flow limit requires the diameter-based particle Reynolds number ${\mathrm{Re}}_d=2a\rho U/\mu \ll 1$. The Stokes drag law gives $F_D=6\pi \mu aU$. Substitution into the definition with $A=\pi a^2$ gives the short derivationThus the drag coefficient is neither necessarily constant nor bounded above by one. The formula containing $U^2$ still describes linear drag here, because the drag coefficient varies as $1/U$. The relation $C_D=24/{\mathrm{Re}}_d$ is the leading small-Reynolds number approximation, not a measured curve valid at every Reynolds number. ^[ref]

For a worked example, stipulate dynamic viscosity $\mu =0.10\mathrm{Pa}\mathrm{s}$, density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, radius $a=10^{-4}\mathrm{m}$ and relative speed $U=10^{-3}\mathrm{m}{\mathrm{s}}^{-1}$. Within the stated Stokes drag law assumptions,Doubling the relative speed with all other parameters fixed gives ${\mathrm{Re}}_d=0.004$ and $C_D=6000$: the predicted drag force doubles to $3.76991118\times 10^{-7}\mathrm{N}$. Falsely holding the drag coefficient at $12000$ would instead quadruple the drag force. Both particle Reynolds numbers remain small, so the comparison is within the creeping-flow approximation.

The useful idea is to separate the force scale from the dependence on geometry and flow conditions. A drag coefficient can help compare experiments, but changing Reynolds number can change the coefficient itself; NASA Glenn explicitly stresses matching the relevant flow conditions. The elegant cancellation above also exposes a common mistake: a drag coefficient definition does not by itself impose a quadratic drag law.

The dynamic viscosity $\mu$ is the coefficient relating shear stress to shear rate in a Newtonian fluid. For a simple shear flow $\mathbf{u}=(u(y),0,0)$,Here ${\tau }_{xy}$ is the $x$ component of viscous traction on a plane with normal vector $+{\mathbf{e}}_y$. The shear rate $\dot{\gamma }=du/dy$ has units ${\mathrm{s}}^{-1}$. Thus dynamic viscosity measures the stress required to sustain a given rate of deformation. It is the shear coefficient; the separate coefficient for volume change is volume viscosity. ^[ref]

For a worked Couette flow example, place a Newtonian fluid of constant dynamic viscosity $\mu =0.10\mathrm{Pa}\mathrm{s}$ between parallel plates separated by $h=0.0010\mathrm{m}$. Hold the lower plate fixed and move the upper plate at $U=0.020\mathrm{m}{\mathrm{s}}^{-1}$. Assume steady flow, no streamwise pressure gradient, and the no-slip boundary condition. ThenThe fluid pulls backward on the moving plate with force magnitude $F={\tau }_{xy}A$, where $A$ is its wetted area. For $A=0.010{\mathrm{m}}^2$ this is $0.020\mathrm{N}$. Measuring that force, the gap and the plate velocity determines $\mu =Fh/(AU)$: a direct way to measure dynamic viscosity. The linear profile and opposing wall tractions are the standard plane Couette flow solution. ^[ref]

The velocity gradient matters, rather than the velocity itself. Adding the same constant velocity to both plates and the fluid leaves $du/dy$, and hence shear stress, unchanged. In uniform flow, every layer translates together and $du/dy=0$. More generally, even a nonzero velocity gradient need not generate viscous stress: solid-body rotation has zero strain-rate tensor, as the viscous stress tensor example shows. Thin gaps amplify shear stress for a fixed plate-speed difference, making dynamic viscosity useful for predicting the force needed to move lubricated surfaces.

A rotating inner cylinder of radius $R$ and height $H$ has a uniform gap $h$ from a fixed housing at its side and both ends. For $h\ll R,H$, no-slip boundary condition and laminar Couette flow give side torque $2\pi \eta \omega R^{3}H/h$. Each end contributes ${\int }_0^R(\eta \omega s/h)(2\pi sds)s=\pi \eta \omega R^4/(2h)$. Summing both end faces proves the displayed expression. Shaft area, corner regions, gas slip and viscous heating are neglected.

A force-driven rigid particle with surface $r=a+\epsilon f(\mathbf{n})$ has a leading translating sphere in Stokes flow with ${\mathbf{U}}_0=\mathbf{F}/(6\pi \mu a)$. Expanding the no-slip boundary condition gives the first-order surface velocity ${\mathbf{u}}_1={\mathbf{U}}_1+{\mathbf{\Omega }}_1\times \mathbf{x}+\mathbf{g}$, where $\mathbf{g}=3f(\mathbf{I}-\mathbf{nn}){\mathbf{U}}_0/(2a)$. The first-order resultant force and torque vanish when the applied load is fixed. The Lorentz reciprocal theorem for Stokes flow, using the uniform surface traction in translating-sphere Stokes flow, therefore gives ${\mathbf{U}}_1=-\int \mathbf{g}dS/(4\pi a^2)$. For $f=a\mathbf{n}\cdot \mathbf{Dn}$ with symmetric traceless $\mathbf{D}$, the isotropic fourth moment gives ${\mathbf{U}}_1=\mathbf{D}{\mathbf{U}}_0/5$. A rotating-sphere test similarly gives ${\mathbf{\Omega }}_1=-3\int \mathbf{n}\times \mathbf{g}dS/(8\pi a^3)$, which vanishes for an inversion-symmetric deformation.

For a dimensionless Navier-slip Taylor swimming sheet, the decaying first-harmonic streamfunction is $(1+Y)e^{-Y}\sin (X-\tau )$. Its tangential surface velocity and surface shear both vanish at $Y=0$. Hence it obeys the first-order Navier slip boundary condition for every nonnegative slip length, including the no-slip boundary condition.

The Louisiana Tech LAB 8 - Velocity Profile above a Flat Plate describes a 14-inch plate in a wind tunnel. A Pitot-static tube of 0.10-inch outside diameter traverses normally about 12 inches downstream; its screw scale has 0.001-inch increments. A water-column manometer reads the difference between stagnation pressure and static pressure. The procedure establishes a wall position, zeros the tunnel indicator, sets fan speed, and records at least ten positions inside the boundary layer, together with room temperature and barometric pressure. The instrument list includes a flow-stream thermometer implemented with a temperature-sensitive resistor. The manual requests comparison of outer velocity with the tunnel indicator and of the recovered velocity profile with laminar or turbulent predictions; it provides no observed dataset. These are apparatus specifications and instructions, not a reported confirmation of laminar flow.

Under incompressible flow conditions, an ideal aligned Pitot tube gives $q(y)=p_0(y)-p(y)=\rho u(y{)}^2/2$. Convert the manometer head $h$ to $q=({\rho }_{\mathrm{water}}-\rho )gh$, then infer $u=\sqrt{2q/\rho }$. Here $g$ is gravitational acceleration; mass density $\rho$ can be estimated from temperature and absolute pressure using the ideal gas law. A normalized velocity profile rises from the unresolved near-wall region to a plateau, whose value defines $U$; downstream profiles should broaden and collapse when plotted against $\eta$ if the Blasius boundary layer assumptions hold. Integrating a recovered velocity profile gives displacement thickness and momentum thickness.

Additional controls for interpreting this experiment follow from the measurement model: check differential-pressure calibration with equal port pressures and known liquid heads; leak-check tubing; align the Pitot tube with the incoming velocity; convert wall contact to the inlet-centre coordinate rather than treating outside contact as $y=0$; repeat traverses and outer velocity readings to detect drift. Keep the plate smooth and aligned, and check outer static pressure at multiple downstream stations before assuming zero pressure gradient. Repeating at another station or outer velocity tests the predicted boundary-layer scaling more directly than fitting one curve. These controls are proposed checks, not actions documented as completed in the manual.

Finite probe size, wall interference, flow fluctuations, and weak near-wall dynamic pressure limit inference. Even the manual's probe-based Reynolds number criterion, ${\mathrm{Re}}_d=ud/\nu >1000$ for its ideal Bernoulli equation conversion, must be checked with local $u$, not just $U$. In the illustrative model at $x=0.3048\mathrm{m}$, ${\delta }_{99}\simeq 3.32\mathrm{mm}$; the quoted outside diameter is $2.54\mathrm{mm}$, a substantial fraction of that thickness. A real inlet's sensing footprint and correction cannot be deduced from outside diameter alone. For an ideal footprint that averages dynamic pressure, the inferred velocity is $\sqrt{⟨u^2⟩}$ rather than $⟨u⟩$. Differentiating such sparse readings to obtain wall shear stress magnifies measurement error; the momentum thickness route is less sensitive to this differentiation but still needs the unresolved near-wall contribution. The no-slip boundary condition is a model constraint here, not a direct reading from a probe centred on the wall.

A fluid is a material that flows under sustained shear stress and cannot maintain a static shear stress in its fluid state. Liquids and gases are the usual examples. A Newtonian fluid has deviatoric stress proportional to its rate-of-strain tensor, whereas more general fluid constitutive equations can be nonlinear or have memory. In a porous medium, the fluid occupies connected voids within a load-bearing solid framework. The important distinction is the response to a continuing tangential force: an ordinary elastic solid can settle into a fixed deformation, while an ordinary fluid continues to deform. A fluid can nevertheless be at rest when no unbalanced shear stress drives it. ^[ref]

For a worked example, place a Newtonian fluid between parallel plates separated by $h=1\mathrm{mm}$, keep the lower plate fixed and move the upper one at $U=0.10\mathrm{m}{\mathrm{s}}^{-1}$. With no-slip boundary conditions, steady flow, no imposed pressure gradient and negligible edge effects, Couette flow has $u(y)=Uy/h$. For dynamic viscosity $\mu =0.001\mathrm{Pa}\mathrm{s}$, the shear stress magnitude isThe velocity gradient describes neighbouring layers sliding at different speeds. Doubling $U$ doubles the required shear stress; maintaining the same velocity keeps producing deformation, rather than fixing a final shape. This gives viscosity a concrete meaning: resistance to a rate of motion, not a prohibition on motion.

The second-order mean no-slip boundary condition correction contains products of first-order deformation modes and velocity gradients. Distinct integer wavenumbers have zero averaged mixed products. Their separate mean longitudinal and transverse contributions therefore add, even though mixed products create nonzero oscillatory flow modes.

For a unidirectional incompressible flow between stationary plates a distance $H>0$ apart, removing the driving pressure gradient leaves $u_t=\nu u_{yy}$ when the kinematic viscosity $\nu >0$ is constant. The no-slip boundary condition gives $u(0,t)=u(H,t)=0$. For an initial profile $g$ in L2 space, separation of variables and the Fourier sine series give the displayed solution, withEach normal mode decays at a rate proportional to $n^2$. The lowest mode present dominates at late times; when $b_1\ne 0$, its time constant is $H^2/(\nu {\pi }^2)$, consistent with the viscous diffusion time. The solution approaches the initial data in L2 space; compatible smooth profiles also give the familiar pointwise initial condition. This is a heat equation for momentum diffusion, rather than temperature.

A Hartmann layer is the thin viscous-magnetic boundary layer next to a wall crossed by the applied magnetic field. For channel half-width $L$ and large Hartmann number $H$, its thickness is $\delta \sim L/H=\sqrt{{\mu }_0\rho \nu \eta }/|B_0|$. It restores the no-slip boundary condition and the specified magnetic wall condition outside the nearly uniform Hartmann flow core.

For $f>0$, a stationary lower no-slip boundary condition below uniform exterior geostrophic flow $W_g$ produces $W_E=-W_g{e}^{-(1+i)z/{\delta }_E}$. The upward shear component at the bottom is $T_b=\rho \nu W_E^{\prime }(0)=\rho \sqrt{f\nu /2}(1+i)W_g$; the bottom traction on the fluid is $-T_b$. This law cannot be assigned to an arbitrary prescribed surface wind stress without an additional boundary condition.

Laminar flow is a regime of fluid flow in which the motion remains organized, with neighboring fluid parcels following smooth, coherent paths rather than the irregular interacting fluctuations characteristic of turbulence. In a straight pipe, the simplest example consists of coaxial layers sliding past one another. This describes the character of the motion, not a requirement that every parcel have the same velocity. Nor does it require steady flow: a smooth profile can evolve in time as viscosity transmits motion from a newly moving wall into a resting fluid. The unsteady wall-driven solution provides an explicit example. ^[ref]

The Reynolds pipe-and-dye experiment makes this distinction visible. A thin injected dye band can remain coherent in laminar flow, whereas irregular motion can fold it and carry dye across the pipe. However, molecular diffusion can spread dye even in laminar flow, and an inlet disturbance can move a coherent band without producing turbulence. A dye band is a streakline, not a measurement of the entire velocity field; in steady flow an ideal passive tracer's streakline coincides with a streamline. Reynolds explicitly distinguished a wandering band from one broken into eddies. ^[ref]

For a circular pipe of internal diameter $D$, define the mean axial velocity $U$ using the volumetric flow rate $Q$ and cross-sectional area $A$:Here $\rho$ is mass density, $\mu$ is dynamic viscosity and $\nu =\mu /\rho$ is kinematic viscosity. This particular Reynolds number uses the diameter and mean velocity, not the radius or centreline velocity. Small Reynolds number generally favors laminar flow, but a single universal critical Reynolds number does not distinguish it from turbulence in all geometries and experiments. Even in a given pipe, inlet conditions and the magnitude of disturbances affect the observed onset. Reynolds demonstrated this dependence by deliberately inserting a wire coil. ^[ref]

The quantitative example above assumes steady flow, incompressible flow, a Newtonian fluid with constant dynamic viscosity, a straight horizontal circular pipe of radius $R=D/2$, stationary walls satisfying the no-slip boundary condition, and fully developed flow. The last assumption means that the axial velocity profile no longer changes downstream. A positive pressure drop $\Delta p=p_{\mathrm{in}}-p_{\mathrm{out}}$ over a developed length $L$ gives the Hagen-Poiseuille equation:The velocity is zero at the wall and maximal at the centre. Averaging over the circular area, whose annuli have area $2\pi rdr$, gives $U=u(0)/2$; averaging along one diameter would give a different answer. These assumptions produce a parabola, rather than every laminar flow having a parabolic profile. ^[ref]

As a stipulated calculation, take $D=0.010\mathrm{m}$, $U=0.10\mathrm{m}{\mathrm{s}}^{-1}$, $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, $\mu =0.0010\mathrm{Pa}\mathrm{s}$ and $L=1.0\mathrm{m}$, and assume the developed laminar flow above. Then $\nu =10^{-6}{\mathrm{m}}^2{\mathrm{s}}^{-1}$, ${\mathrm{Re}}_D=1000$, $u(0)=0.20\mathrm{m}{\mathrm{s}}^{-1}$, $Q=7.85\times 10^{-6}{\mathrm{m}}^3{\mathrm{s}}^{-1}$ (about $0.471\mathrm{L}\mathrm{mi}{\mathrm{n}}^{-1}$), and $\Delta p=32\mu LU/D^2=32\mathrm{Pa}$. These are predictions from chosen parameters, not measurements or a claim that ${\mathrm{Re}}_D=1000$ guarantees laminar flow in every apparatus.

The beauty of this example is that a spatially varying velocity field still gives an exact relation between two accessible measurements, pressure drop and volumetric flow rate. It predicts the hydraulic resistance $\Delta p/Q=8\mu L/(\pi R^4)$. Doubling the radius therefore multiplies $Q$ by sixteen at fixed pressure drop, length and dynamic viscosity, provided the model remains applicable. This makes the laminar flow model useful for sizing narrow tubes or determining dynamic viscosity. It also explains why a persistent narrow dye band can be informative without implying that the fluid is moving as a uniform plug.

For a thin viscous gravity current with hydrostatic pressure $p=\rho g(h-z)$, a no-slip boundary condition at its base and a stress-free boundary condition at its top, horizontal momentum balance gives $u=(g/\nu )h_x(z^2/2-hz)$. Integrating across the depth gives the volume flux $q=-(g/3\nu )h^3{h}_x$. Horizontal spreading follows $h_t+q_x=0$ if the substrate is impermeable.