Codex @codex  0

Consider planar horizontal layers with positive scalar intrinsic intrinsic permeabilities $k_i$, thicknesses $b_i$, and fractions $f_i=b_i/{\sum }_j{b}_j$. Assume each layer is uniform and isotropic, fully saturated by the same Newtonian fluid of constant mass density and dynamic viscosity, with perfect hydraulic contact and steady incompressible flow. Use boundaries that impose pure flow normal to the layers or a common gradient of hydraulic head along them; additional interface resistance and finite-edge effects are excluded.

For vertical flow normal to the layers, every layer carries the same normal Darcy flux by mass conservation. The hydraulic head drops add. Writing $B={\sum }_i{b}_i$ and taking positive $\Delta h_i$ as the drop along the flow direction givesThis is the thickness-weighted harmonic mean. It is hydraulic head, rather than fluid pressure alone, that accounts correctly for gravity in the vertical calculation. A thin layer of very small permeability of a porous medium can control the whole discharge because all flow must pass through it.

For horizontal flow parallel to the layers, each layer experiences the same gradient of hydraulic head. For width $W$, each area is $A_i=W{b}_i$, and the individual volume fluxes add:This is the thickness-weighted arithmetic mean. A highly permeable layer can carry most of the horizontal flow because it supplies a continuous route from inlet to outlet. Woessner and Poeter derive both cases for hydraulic conductivity; identical fluid properties allow the common factor $\rho g/\mu$ to cancel.^[ref]

For two equally thick layers, take $k_1=k_0$ and $k_2=9k_0$, where $k_0=10^{-12}{\mathrm{m}}^2$. The values 1 and 9 are therefore permeabilities in units of $k_0$, rather than dimensionless intrinsic intrinsic permeabilities. ThenAcross the layers, the less permeable half takes $90\%$ of the total hydraulic head drop and the more permeable half only $10\%$. Along the layers, the more permeable half carries $90\%$ of the total volume flux. At the same magnitude of imposed gradient of hydraulic head, horizontal Darcy flux is $5/1.8=25/9\simeq 2.78$ times the vertical Darcy flux.

The arrangement thus creates physical anisotropy from locally isotropic ingredients. In coordinates with $x,y$ along the layers and $z$ normal to them, the effective permeability of a porous medium tensor is ${\mathbf{k}}_{\mathrm{eff}}=\mathrm{diag}(5,5,1.8)k_0$. A driving gradient of hydraulic head with both horizontal and vertical components generally produces a discharge in a different direction. This distinction matters when deciding whether a geological bed conducts aquifer flow along itself or acts as a barrier to flow across it.

Regenerate the diagram by running python ../layered-porous-permeability.py from _media; it writes its PNG into the current directory.

Lift is the component of the fluid force perpendicular to the far-field flow relative to the body. It need not point vertically: a tilted flow changes the lift direction, and a banked wing can produce a sideways component. In a planar problem choose a unit vector ${\mathbf{e}}_D$ along the incoming relative velocity ${\mathbf{U}}_{\infty }$ and a perpendicular unit vector ${\mathbf{e}}_L$ pointing towards the chosen positive-lift side. ThenHere $D$ is drag and $L$ is signed lift; $L<0$ means a force towards the opposite side. In three dimensions the plane perpendicular to ${\mathbf{U}}_{\infty }$ also admits a separate side-force direction, so a lift axis must be specified. Lift is also distinct from buoyancy, which can exist without a relative flow. ^[ref]

The force comes from adding local surface tractions, rather than assigning a single pressure to each side. Let $\Sigma$ be the wetted body surface and $\mathbf{n}$ its unit normal vector pointing out of the body into the fluid. The pressure contribution pushes inwards; the viscous stress tensor also contributes, including tangential shear stress:Using a constant reference pressure $p_{\infty }$ is legitimate on a closed surface because ${\int }_{\Sigma }\mathbf{n}dA=0$. This surface integral makes the direction explicit: a surface-normal force is generally neither pure lift nor pure drag. Pressure and skin shear stress both enter the resultant, although the pressure term usually dominates lift on an airfoil. ^[ref]

The lift coefficient isolates a useful scale. For a finite wing of reference area $S$ (normally its projected planform area), and for a two-dimensional airfoil section of chord length $c$, respectively,Here $q_{\infty }$ is dynamic pressure, ${\rho }_{\infty }$ is far-field mass density, and $L^{\prime }$ is lift per unit span, with units $\mathrm{N}{\mathrm{m}}^{-1}$. The signs of $C_L$ and $c_l$ follow the chosen lift axis. For stipulated values ${\rho }_{\infty }=1.20\mathrm{kg}{\mathrm{m}}^{-3}$, $U_{\infty }=20\mathrm{m}{\mathrm{s}}^{-1}$, $S=0.50{\mathrm{m}}^2$ and $C_L=0.60$, $q_{\infty }=240\mathrm{Pa}$ and $L=72\mathrm{N}$. Doubling $U_{\infty }$ gives four times the lift only if $C_L$ stays fixed: changes in Reynolds number, Mach number, surface condition or angle of attack can change it. ^[ref]

A motivating experiment is to compare two independently observed quantities: the lift registered by a force balance, and the lift inferred from surface pressures. A low-speed wind tunnel holds the model fixed while air passes it at adjustable speed; a calibrated balance measures forces, while small flush holes on both surfaces connect by tubing to differential-pressure sensors referenced to $p_{\infty }$. Iowa State's GA(W)-1 teaching model has $43$ taps around a $101\mathrm{mm}$ chord; its instructions give their coordinates and close the contour by repeating the first tap. A classroom comparison can combine this pressure-panel method with a balance on the same model; this is an experiment design, not a claim that those instructions report a combined measurement. ^{[ref] [ref]}

Before turning on the wind, apply known positive and negative loads to the balance through its resolving centre; use $mg$ for the known load from a mass. Record all balance channels during loading and unloading, including response to drag and torque loads, to detect cross coupling, hysteresis and zero drift. Record a wind-off tare at each model orientation so its weight does not become an aerodynamic lift reading. Map tunnel flow direction and uniformity, rather than using the fan setting as a speed measurement. Covert's MIT guide describes these calibration, tare and flow-survey requirements. ^[ref]

For each pressure channel, apply zero and several known pressure differences spanning the expected range; determine the voltage-to-pressure relation and recheck it after a run. Inspect tubing for leaks and blocked taps, and wait for readings to settle before averaging: tubing can delay a pressure signal. This calibration principle is used in NASA's electronically scanned pressure systems, where several applied calibration pressures determine each sensor's response and repeated calibration tracks drift. A simple pressure standard can instead use a liquid-column height difference and the hydrostatic pressure law. ^[ref]

Measure far-field static and total pressure outside the model's fluid wake; in low-speed incompressible flow their difference gives $q_{\infty }$ by the Bernoulli equation, so $U_{\infty }=\sqrt{2q_{\infty }/{\rho }_{\infty }}$. Record temperature and ambient pressure to determine mass density using the ideal gas law. At fixed $q_{\infty }$, step the angle of attack, repeat selected settings, and sweep back down to check reproducibility. Record balance lift, each tap's mean pressure, and fluctuation sizes. NASA describes pressure probes connected to differential transducers and emphasizes recording their reference pressure. ^[ref]

Convert each tap to a pressure coefficient $C_{p,i}=(p_i-p_{\infty })/q_{\infty }$. For consecutive taps bounding a panel of length $\Delta s_i$ with outward normal vector ${\mathbf{n}}_i$, a trapezoidal rule estimate of the pressure force per unit span isCompare this with balance lift divided by span only when the spanwise flow is approximately uniform. Agreement tests the pressure-resultant picture; disagreement calls for checking unresolved leading-edge pressure variations, skin shear stress, support loads and three-dimensional end effects before rejecting it. A force balance senses the total load, while taps alone reconstruct the pressure part. Tunnel walls and model blockage alter the effective flow and require corrections; end plates or a model extending between walls approximate two-dimensional flow but do not guarantee it. Record the Reynolds number and Mach number with the result, and do not interpret tunnel lift coefficients as universal wing data. ^{[ref] [ref]}

The Navier slip boundary condition allows a fluid to slide tangentially past a solid wall while resisting that motion through wall friction in the Navier slip boundary condition. For parallel flow past a flat, impermeable wall translating uniformly with tangential velocity $U_w$, write the fluid's tangential velocity as $u$ and its relative tangential velocity as $u_{\mathrm{rel}}=u-U_w$. Thenwhere $b$ is the slip length, also written ${\ell }_s$. The unit normal $\mathbf{n}$ points from the wall into the fluid, and ${\partial }_{n}u=\mathbf{n}\cdot \nabla u$ is the normal derivative. The separate no-penetration boundary condition requires $(\mathbf{u}-{\mathbf{U}}_w)\cdot \mathbf{n}=0$: tangential sliding does not allow flow through the wall.

The sign convention matters. In a channel $0<z<H$, the lower wall has inward unit normal $+{\mathbf{e}}_z$ and the upper wall has inward unit normal $-{\mathbf{e}}_z$. If both walls obey the Navier slip boundary condition, their tangential velocities satisfyHere the prime means the derivative with respect to $z$. Using $+b_H{u}^{\prime }(H)$ at the upper wall would reverse the friction law. The examples Couette flow with a slipping lower wall and planar Poiseuille flow with equal wall slip use different upper-wall boundary conditions.

This flat-wall formula is a local linear model. On a curved wall the Newtonian fluid stress tensor must supply the tangential traction; the normal derivative alone generally misses curvature terms. The author-written review by Lauga, Brenner and Stone states the full strain-based condition with the unit normal directed into the liquid. ^[ref]

The same review credits Navier's 1823 treatise with the linear slip law. Its history is a reminder that the no-slip boundary condition is an experimentally successful model, whose validity depends on the interface and scale. The history here follows that review rather than a new reading of Navier's original memoir. ^[ref]

Neutral buoyancy means a fully immersed body's upward buoyancy equals its downward force from gravity. For a body of fixed volume in a constant-mass density fluid under a uniform gravitational field, this occurs when its mean mass density equals the surrounding mass density. A $0.10\mathrm{m}$ cube of mass density $1000\mathrm{kg}{\mathrm{m}}^{-3}$ in the stipulated water has mass $1.000\mathrm{kg}$, so both forces are $9.81\mathrm{N}$.

Ideally it can remain at rest at any depth for which it stays fully immersed and clear of boundaries. Raising or lowering it changes the absolute pressure, but neither its displaced volume nor the difference between the vertical pressure forces. There is therefore no vertical restoring force toward a preferred depth in this model. Neutral buoyancy is a vertical force balance, not a claim about resistance to tilting or about a compressible body whose volume changes with pressure.

A pathline records the motion of one identified fluid element: keep its starting position $\mathbf{a}$ and starting time $t_0$ fixed, and let physical time $t$ advance. In a velocity field $\mathbf{u}(\mathbf{x},t)$ its Lagrangian trajectory $\mathbf{X}(t;t_0,\mathbf{a})$ solves the initial value problemA streamline instead freezes the time argument of the velocity field; a streakline fixes the observation time and varies the release time of many fluid elements. A pathline is useful for answering where a particular parcel goes and how long it takes, rather than merely which way the fluid is moving now. ^[ref]

For a concrete analytical model, take spatial coordinates $x,y$ in metres, time $t\ge 0$ in seconds, and the spatially uniform but time-dependent velocity fieldIts numerical components in these units are $(1,t)$. A fluid element starting at the origin at $t_0=0$ has the pathlineAt $t=1\mathrm{s}$ it is at $(1,0.5)\mathrm{m}$; at $t=2\mathrm{s}$ it is at $(2,2)\mathrm{m}$. The parabola records earlier, shallower velocity directions as well as the final direction. Differentiating gives $(\dot{x},\dot{y})=(U,at)$ and acceleration $(0,a)$, so this curved pathline needs no spatial variation of the velocity field. More generally the chain rule along a pathline gives $df(\mathbf{X}(t),t)/dt={\partial }_{t}f+\mathbf{u}\cdot \nabla f$, the material derivative of a transported scalar field.

A motivating tracer demonstration appears in S. J. Kline's film notes: water passes through a contracting channel, and hydrogen bubbles mark a small element at a wire at time $t_0$. The illustrated element is shown again $4.7\mathrm{s}$ later, with its accumulated pathline superposed. Following that same marked element across successive film frames records a history; the entire bubble pattern in one frame has a different meaning. ^[ref]

For a modern version of the readout, use sparse small markers in a transparent water channel, illuminate the viewed plane, and keep an optical camera fixed with a ruler in that plane. Record timestamps and track the same isolated marker between frames; successive positions estimate its pathline and $[\mathbf{X}(t+\Delta t)-\mathbf{X}(t)]/\Delta t$ estimates its velocity. A long exposure records a pathline segment but loses the timing within it, so its total length divided by exposure time gives only mean speed. These are proposed measurement steps based on Jean Hertzberg's teaching guidance, rather than additional measurements from Kline's film. Markers leaving the viewed plane, overlapping tracks and perspective can break particle identification; repeat with a thinner viewed slice and a lower marker concentration. ^[ref]

Permeability of a porous medium is the coefficient $k$ that measures how readily its connected pore geometry transmits a fluid in Darcy law. Its units are ${\mathrm{m}}^2$, since $k/\mu$ converts a pressure gradient in $\mathrm{Pa}{\mathrm{m}}^{-1}$ into a Darcy velocity in $\mathrm{m}{\mathrm{s}}^{-1}$. Large open, well-connected passages generally give larger permeability of a porous medium than narrow or winding passages. Porosity alone cannot determine it: two samples can contain the same void volume but have different passage sizes and connectivity.

For an unchanged framework and the basic no-slip Newtonian fluid model, intrinsic permeability of a porous medium does not contain dynamic viscosity or mass density. These properties instead enter hydraulic conductivity $K=k\rho g/\mu$. At the same horizontal pressure gradient, replacing a liquid with another of twice the dynamic viscosity halves Darcy flux even though $k$ stays fixed. Structural changes, such as swelling or clogging, can change $k$; the independence claim assumes the pore geometry stays fixed.^[ref]

An isotropic sample has scalar permeability of a porous medium. An anisotropic sample requires a tensor $\mathbf{k}$ in $\mathbf{q}=-\mathbf{k}(\nabla p-\rho \mathbf{g})/\mu$, so discharge need not be parallel to the net driving force. Layered porous permeability gives a concrete example: arranging individually isotropic materials in horizontal layers produces different vertical and horizontal effective permeabilities.

Now both plates at $z=0$ and $z=H$ are stationary and both obey the Navier slip boundary condition with the same slip length $b$. Drive steady flow of a Newtonian fluid in incompressible flow by the constant pressure gradient $G=-dp/dx>0$, with constant dynamic viscosity $\mu >0$ and no other streamwise body force. The planar Poiseuille flow profile and its volumetric flow rate per unit width areThe full volumetric flow rate for channel width $W$ is $Q=Wq$ when side-wall effects are negligible. Thus $q$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$, while $Q$ has units ${\mathrm{m}}^3{\mathrm{s}}^{-1}$. Denote the value at $b=0$ by $q_0=G{H}^3/(12\mu )$.

The two wall checks expose the normal signs:Both Navier slip boundary conditions therefore hold exactly. The two walls exert backward tangential tractions of magnitude $GH/2$ on the fluid. Increasing $b$ at fixed $G$ adds a uniform velocity offset, leaving the velocity gradient and wall shear stress unchanged. This differs from Couette flow with a slipping lower wall, where upper-wall velocity is fixed and shear stress decreases.

For $H=10\mu \mathrm{m}$, $b=2\mu \mathrm{m}$, $G=10^3\mathrm{Pa}{\mathrm{m}}^{-1}$ and $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$, the wall velocity is $10\mu \mathrm{m}{\mathrm{s}}^{-1}$, the centre velocity is $22.5\mu \mathrm{m}{\mathrm{s}}^{-1}$, and $q=1.83\times 10^{-10}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. The enhancement is $q/q_0=2.2$. Equivalently, a prescribed volumetric flow rate needs only $1/2.2$ of the original pressure gradient. This makes slip length useful for estimating pumping requirements in small channels.

At fixed $H$, $G>0$ and $\mu$, $b\to 0$ recovers ordinary planar Poiseuille flow. As $b\to \infty$, the velocity and $q$ grow without bound in this model, although the varying part of the profile stays the same. At exactly zero tangential wall traction, no steady solution can balance this nonzero pressure gradient: $u^{\prime }(0)=u^{\prime }(H)=0$ conflicts with $\mu [u^{\prime }(H)-u^{\prime }(0)]=-GH$. Unlimited enhancement is consequently not a prediction for a real channel; finite wall friction, finite driving capacity and the validity of the assumed steady flow must be assessed. At fixed finite $q$ instead, $G\to 0$ and the profile approaches uniform velocity $q/H$.

A short derivation integrates the reduced Navier-Stokes equation $\mu u^{\prime \prime }=-G$ twice. Equal slip lengths and the opposite inward unit normals give $u=G[z(H-z)+bH]/(2\mu )$. Its integral over the gap supplies the enhancement factor $1+6b/H$.

Pore velocity is the fluid velocity averaged over pore volume relative to the solid framework. For uniform porosity $\varphi$, full saturation and all pores participating, $\mathbf{v}=\mathbf{q}/\varphi$, with Darcy velocity $\mathbf{q}$ defined per bulk area. The identity follows from volume averaging: extending the relative fluid velocity by zero in the solid gives a bulk average equal to the pore fraction times the pore-volume average. It does not assert that every local passage has open area $\varphi A$ or identical speed.

In a uniform advective model, travel over length $L$ takes $L/v=\varphi L/q$. If some pore volume is inaccessible to transport, the corresponding estimate uses effective porosity ${\varphi }_e$ and $v=q/{\varphi }_e$. A connected but nearly stagnant passage can exchange solute slowly, so even this mean does not describe all arrival times. Flux-weighted residence time in a porous layer explains one precise averaging convention.^[ref]

Porosity is the fraction of a representative bulk material volume occupied by voids:It is dimensionless and measures storage space. A sample with bulk volume $1.0\mathrm{L}$ and pore volume $0.25\mathrm{L}$ has porosity $\varphi =0.25=25\%$. If all its pores are saturated by one liquid, the liquid occupies $0.25\mathrm{L}$; partial fluid saturation $s$ instead gives liquid volume $s\varphi V_{\mathrm{bulk}}$.

Porosity and permeability of a porous medium answer different questions: how much fluid fits, and how easily it flows. Isolated bubbles or cavities can contribute to total porosity while providing no throughgoing path. Even connected passages may be narrow enough to transmit very little discharge. In the simple fully participating saturated model, Darcy velocity is $\varphi$ times pore velocity; changing porosity at a stipulated fixed Darcy flux changes ideal advective travel time, but porosity alone cannot predict the Darcy flux.

Porous-media flow is fluid motion through the interconnected voids of a porous medium, such as a sand bed, a rock or a filter. Individual passages can be tiny and tortuous, yet their combined discharge can be predicted without resolving every passage. Darcy law supplies this bulk description when viscous resistance dominates pore-scale inertia; porosity measures how much space stores fluid, while permeability of a porous medium measures how readily that space transmits it.

The basic fully saturated model assumes one Newtonian fluid fills the connected pore space of a stationary solid framework, with incompressible flow, a small pore-scale Reynolds number, and averaging volumes large compared with individual pores but small compared with macroscopic variations. These are assumptions of the model developed below, rather than a claim that Darcy law describes every porous-media flow. Two-phase porous-media flow and poroelastic aquifers require additional constitutive and storage information.

For a first calculation, see Darcy flow through a homogeneous slab: a pressure drop determines the discharge, but porosity determines the corresponding pore velocity. Layered porous permeability then shows how geometry alone can make the bulk response anisotropic.

A porous medium is a solid framework containing connected voids in which a fluid can reside and flow. Its porosity is the void fraction; its permeability of a porous medium measures the ease of pressure-driven flow through the connected pores. The averaged Darcy flux is measured per total cross-sectional area, whereas the pore velocity is that flux divided by porosity. Not every material with voids has connected permeable pores.

The pressure coefficient compares a local static pressure with a specified reference static pressure, using a fixed reference dynamic pressure:The reference mass density, speed, pressure and location must accompany a reported pressure coefficient. Keep the same $q_{\infty }$ throughout a surface distribution; replacing it by local $\rho U^2/2$ gives a different quantity and fails at a stagnation point. The definition itself requires neither Bernoulli's principle nor inviscid flow, so measured pressure coefficients remain useful when those idealizations fail.

Positive $C_p$ means $p>p_{\infty }$; zero means equality; negative $C_p$ means $p<p_{\infty }$. This sign does not specify the sign of absolute pressure. For stipulated $p_{\infty }=101325\mathrm{Pa}$ absolute and $q_{\infty }=240\mathrm{Pa}$, $C_p=-3$ gives $p=100605\mathrm{Pa}$ absolute: a $720\mathrm{Pa}$ deficit, not a negative absolute pressure. If $p_{\infty }$ equals atmospheric pressure, the numerator is gauge pressure; a wind tunnel's reference static pressure need not equal room atmospheric pressure. These are reference choices, not different physical pressures.

Under the same-height, constant-mass density $\rho ={\rho }_{\infty }$, steady flow and inviscid flow hypotheses of the Bernoulli equation, and with the same Bernoulli invariant as the reference flow,The first equation immediately gives the second after subtracting $p_{\infty }$ and dividing by $q_{\infty }$. Hence a lossless stagnation point has $C_p=1$, $U=U_{\infty }$ gives $C_p=0$, and $U=2U_{\infty }$ gives $C_p=-3$. For a real stationary wall satisfying the no-slip boundary condition, the velocity is zero right at the wall; inserting that wall velocity in an inviscid flow formula would wrongly predict $C_p=1$ everywhere. A thin boundary layer approximation instead relates the wall pressure to the outer-flow pressure and uses the outer velocity where the Bernoulli equation applies.

For potential flow around a circular cylinder of radius $a$ with zero circulation, the surface tangential velocity is $u_{\theta }=-2U_{\infty }\sin \theta$. Thus the ideal pressure coefficient isWith $\theta$ measured from downstream, the front is at $\pi$ and the rear at zero. Both ideal stagnation points have $C_p=1$, while the sides at $\pi /2$ and $3\pi /2$ have $C_p=-3$. This exact model is a useful baseline for an experiment; its complete rearward pressure recovery is not a prediction for a separated real fluid wake.

To obtain a force, integrate the local normal contributions with their directions. Let $\mathbf{n}$ be the unit normal vector from the body into the fluid and $\Sigma$ its closed wetted surface:The second equality removes a spatially constant reference pressure over the closed surface; it does not remove a nonuniform hydrostatic pressure or justify dropping the reference term on an isolated open patch. A negative $C_p$ contributes outwards in this reference-subtracted calculation, even though a positive absolute pressure still pushes inwards locally. Project ${\mathbf{F}}_p$ onto the chosen lift and drag directions. The full fluid force also includes ${\int }_{\Sigma }{\mathbf{\sigma }}_{\mathrm{visc}}\mathbf{n}dA$ from the viscous stress tensor, including tangential shear stress; pressure tapping does not measure that contribution.

For the two-dimensional cylinder, force per unit span follows from $dA=ad\theta$ per unit span and $\mathbf{n}=(\cos \theta ,\sin \theta )$:Substituting $C_p=1-4{\sin }^2\theta$ gives $(0,0)$: opposite points have identical pressure coefficients and opposite normal vectors. Zero drag despite a strongly nonuniform pressure is the instructive limitation associated with D'Alembert's paradox. Real viscosity, boundary layers and flow separation can prevent rearward pressure recovery and produce pressure drag. Thus the measured curve can test the model and estimate a pressure contribution without proving that total drag vanishes.

A motivating measurement is Virginia Tech's Experiment 3 - Flow Past a Circular Cylinder. Its open-jet wind tunnel holds a rotatable $140\mathrm{mm}$ cylinder with end plates and 36 one-millimetre pressure tappings at ten-degree intervals. A reference Pitot-static tube supplies $p_{0\infty }$ and $p_{\infty }$ to one digital manometer; a second connects a surface tap to its positive port and reference static pressure to its negative port. Record the simultaneous readings:The last equality uses the low-speed stopping relation, not an additional speed measurement. The manual requires sealed tubing and zero checks with the tunnel off. It also provides a water-column manometer: read each column's initial level, and subtract the independently measured reference-to-atmosphere pressure from the tap-to-atmosphere readings. Record tap positions after rotation. Compare the two instruments and check spanwise uniformity; end plates do not guarantee two-dimensional flow. ^[ref]

The inferential distinction matters: the manometers directly measure pressure differences, the ratio estimates pressure coefficients, and a surface integral estimates the pressure part of the force. Measure room absolute pressure and temperature if converting the reference difference to speed through the ideal gas law $\rho =p_{\mathrm{abs}}/(RT)$; use thermodynamic temperature in that equation. For force reconstruction, cover the complete contour and use numerical integration with the tap geometry. Finite tap spacing, unknown end effects, probe disturbance and fluctuations limit that reconstruction; tubes and instruments also average or delay changing pressures. Retain the signed raw readings, repeated readings and reference readings so that these limits can be assessed. A small reference difference magnifies any zero error in the ratio, rather than producing a reliably large $C_p$.

Quadratic drag has magnitude proportional to the square of the relative speed and points opposite to the relative velocity.

For an isotropic model of a body of fixed shape, let $\mathbf{U}={\mathbf{v}}_{\mathrm{body}}-{\mathbf{u}}_{\infty }$ denote its relative velocity through the distant fluid, with $U=|\mathbf{U}|$. A constant positive resistance $k$ givesHere density $\rho$, reference area $A$ and drag coefficient $C_D$ are held fixed over the modeled range of relative speeds. The vector formula makes the direction explicit: the drag force on the body opposes its motion relative to the fluid, and vanishes at $U=0$. For a general body, orientation and other flow conditions can affect the drag coefficient; treating it as a constant requires a regime in which that is an adequate approximation. ^[ref]

Compare two distinct models at a chosen relative speed $U_0>0$. A linear drag model has $F_D(U)/F_D(U_0)=U/U_0$, while a constant-coefficient quadratic drag model has $F_D(U)/F_D(U_0)=(U/U_0{)}^2$. At $U=2U_0$ their drag force ratios are two and four respectively.

This distinction matters when estimating terminal velocity or trajectories such as descent from rest under quadratic drag: doubling relative speed increases the resistance much more strongly in the quadratic drag model. An approximately constant drag coefficient can occur in another Reynolds number regime; its use must be justified there. Conversely, substituting the creeping-flow sphere coefficient $C_D=24/{\mathrm{Re}}_d$ into $F_D=\frac{1}{2}\rho U^{2}A{C}_D$ recovers linear drag, so writing this general formula does not make the resistance quadratic. Neither approximation should be extended across all Reynolds numbers.

The distinction already appears in the original investigation underlying the Stokes drag law, read on 9 December 1850: article 42 identifies resistance proportional to relative speed and distinguishes it from a contribution proportional to the square of relative speed. ^[ref]

The Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous effects. Here $U$ is a characteristic velocity, $L$ a characteristic length, $\mu$ the dynamic viscosity, $\rho$ the mass density, and $\nu =\mu /\rho$ the kinematic viscosity. It is a dimensionless quantity, not a property of the fluid alone: the same material has different Reynolds numbers at different speeds and sizes. Specify the chosen scales when reporting it.

The comparison comes from the Navier-Stokes equation. For velocity variations on scale $L$, convective acceleration is of order $U^2/L$, while the viscous term is of order $\nu U/L^2$. Their ratio is $UL/\nu$. Small Reynolds number favours Stokes flow; large Reynolds number permits an inviscid core, although viscous boundary layers can remain important. Matching Reynolds numbers is useful when comparing geometrically similar flows, provided other relevant dimensionless quantities and boundary conditions also match. ^[ref]

In the constriction example under incompressible flow, use the mean axial velocity and pipe diameter, $L=D=2R$. Taking $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and a stipulated dynamic viscosity $\mu =0.001\mathrm{Pa}\mathrm{s}$ givesHalving the radius at fixed volumetric flow rate doubles this diameter-based Reynolds number, because mean velocity increases fourfold. These values alone neither establish laminar flow nor prove that viscosity can be omitted from the pressure calculation. Increasing $\mu$ to $0.10\mathrm{Pa}\mathrm{s}$ with the same imposed flow gives $40$ and $80$, but strengthens viscous effects; small Reynolds number does not justify Bernoulli's principle.

The name recalls Osborne Reynolds's 1883 experimental investigation of water motion and resistance in pipes. His dye experiments distinguished smooth laminar flow from mixed, irregular turbulence. ^{[ref] [ref]} There is no universal transition value applicable to every geometry and disturbance environment. ^[ref]