Fluid mechanics

Fluid mechanics studies the motion and forces of liquids and gases. Its central simplification is the continuum approximation: replace individual molecules by bulk properties defined at every position and time. The averaging region must contain many molecules while remaining small compared with the distances over which those properties vary. A fluid element is such a small material parcel, rather than a single molecule. This connects fluid mechanics to continuum mechanics and makes calculus useful for describing flow. ^[ref]

Three basic fields are the velocity field $\mathbf{u}(\mathbf{x},t)$, the mass density $\rho (\mathbf{x},t)$ and the pressure $p(\mathbf{x},t)$. The first gives a parcel's instantaneous velocity; the second gives mass per unit volume. The pressure contributes an inward normal traction $-p\mathbf{n}$ across a surface with outward unit normal $\mathbf{n}$. Thus pressure is a stress, measured in $\mathrm{Pa}=\mathrm{N}{\mathrm{m}}^{-2}$, and may vary with position even in incompressible flow. Additional viscous stresses can act tangentially. ^[ref]

The guiding questions are what must be conserved and which forces matter. Mass conservation constrains the velocity field; momentum conservation determines its response to pressure gradients, viscosity and body forces. An incompressible flow model, a Reynolds number estimate and Bernoulli's principle answer different parts of a problem. Their appeal is that a few carefully stated assumptions can determine useful quantities before the entire velocity field is known.

A flow containing two or more distinguishable material phases, such as two immiscible liquids or a gas-liquid mixture. Each phase can have its own volume fraction, density and velocity. Cross-sectional or ensemble averaging produces coupled phase balance equations whose mass, momentum and energy exchange terms must be specified consistently. A common pressure alone does not automatically make an averaged model well posed.

Take incompressible phase mass equations for fractions $\alpha ,1-\alpha$ and momentum equations ${\rho }_i[(u_i{)}_t+u_i(u_i{)}_x]+p_x=0$ with a common pressure. At a uniform state with $0<\alpha <1$, elimination of the mass and momentum perturbation amplitudes givesIts roots are the displayed speeds, whose imaginary parts are $\pm s(u_1-u_2)/(1+s^2)$. A nonzero slip velocity therefore gives a Fourier branch growing as $e^{|k||\mathrm{Im}c|t}$. No finite-order Sobolev estimate can control all high wavenumbers, showing ill-posedness of this unregularized closure. At zero slip the roots coincide and this nonreal-root argument no longer settles the degenerate case. Phase-fraction endpoints also lie outside the root formula's derivation.

Here $⟨\cdot ⟩$ means integration over a fixed pipe cross-section divided by its area, and ${\chi }_1$ is the indicator of one immiscible fluid. The volume fraction is $\alpha$; intrinsic density averages over the phase-occupied area rather than over the whole section. The velocity must be density-weighted for its product with mean phase density to represent mass flux exactly.

If there is no interphase mass transfer, integrating ${\partial }_t(\rho {\chi }_1)+\nabla \cdot (\rho v{\chi }_1)=0$ over a section with impermeable outer wall gives ${\partial }_t({\rho }_1\alpha )+{\partial }_x({\rho }_1\alpha u_1)=0$. The second phase uses $1-\alpha$ and its own intrinsic averages. For constant phase density the velocity reduces to the ordinary phase-area average. These definitions apply where the corresponding phase fraction is positive; zero-fraction intrinsic quantities need a separate convention.

For an incompressible phase, this coefficient compares its actual axial momentum flux with the uniform-profile flux at the same intrinsic mean velocity. Cauchy-Schwarz gives $D_1\ge 1$ where the mean is nonzero. A spatially varying profile can have a larger second moment even when its mean speed is unchanged.

With constant $D_1$, expansion of ${\partial }_t(\alpha u_1)+{\partial }_x(\alpha D_1{u}_1^2)$ and use of ${\alpha }_t+(\alpha u_1{)}_x=0$ gives $\alpha (u_1{)}_t+(2D_1-1)\alpha u_1(u_1{)}_x+(D_1-1)u_1^2{\alpha }_x$. Thus a profile correction changes both convection and phase-fraction-gradient coupling in the averaged momentum equation. If $D_1$ varies, an extra $\alpha u_1^2(D_1{)}_x$ remains; a zero mean should be treated through the undivided second moment rather than an undefined coefficient.

The source ${\Gamma }_i$ is the mass gained by phase $i$ per mixture volume and time, through processes such as evaporation or condensation. Internal exchange conserves total mass, so its phase sources sum to zero. In an immiscible material-interface model without phase change, every ${\Gamma }_i$ is zero. The interface may still move: absence of transfer means that mass does not cross it relative to its own motion.

For a pipe of area $A$, let a phase indicator function be $\chi$. Define $\alpha =A^{-1}\int \chi dA$, $\bar{\rho }=(A\alpha {)}^{-1}\int \chi \rho dA$ and the mass-weighted axial velocity $\bar{u}=(A\alpha \bar{\rho }{)}^{-1}\int \chi \rho u_{x}dA$. Without interfacial mass transfer, integrate ${\partial }_t(\chi \rho )+\nabla \cdot (\chi \rho \mathbf{u})=0$ over a cross-section with impermeable lateral wall. The resulting averaged mass and flux are exactly $\alpha \bar{\rho }$ and $\alpha \bar{\rho }\bar{u}$, proving the displayed law. The complementary phase uses $1-\chi$ and $1-\alpha$. An unweighted phase velocity would not generally factor the mass flux correctly.

The Chézy law relates mean open-channel speed to hydraulic radius and energy slope through a friction coefficient. Equivalently one can use bed shear ${\tau }_b=\rho fu|u|$; for positive uniform flow, balance $\rho gAI={\tau }_b{P}_b$ gives $u=\sqrt{gI/f}\sqrt{A/P_b}$. The absolute value makes drag oppose motion if reverse flows are admitted.

The hydraulic radius is cross-sectional flow area divided by wetted perimeter. A perimeter-averaged boundary shear stress $\tau$ exerts drag per unit fluid volume $\tau \ell /A=\tau /R$. It differs from both actual depth and the radius of a noncircular channel.

A quadratic boundary-drag model relates shear stress to mean flow velocity, with dimensionless resistance coefficient $c_d>0$. Uniform steady flow down slope $S$ balances $\rho gAS=\tau \ell$, giving $|u|=\sqrt{gRS/c_d}$. Its sign is chosen so drag opposes the motion.

For a phase occupying area $\alpha A$ of a constant-area pipe, define $\bar{\rho }=(\alpha A{)}^{-1}\int \rho dA$, $\bar{v}=\int \rho vdA/(\alpha A\bar{\rho })$ and $\overline{v^2}=\int \rho v^{2}dA/(\alpha A\bar{\rho })$. The momentum per length is $A\alpha \bar{\rho }\bar{v}$, the advected momentum flux is $A\alpha \bar{\rho }\overline{v^2}$, and end pressures contribute $-A{\partial }_x(\alpha \bar{p})$. Applying Newton's second law to a slice proves the displayed balance, with $M$ all other phase forces per length. Interphase pressure terms can cancel the derivative of the area fraction. Equality of the second moment and squared mean requires a separate closure.

A floating density instrument. At equilibrium, Archimedes' principle gives $\rho V_{\mathrm{displaced}}=m$ for fixed instrument mass $m$, so it displaces more volume in a less dense liquid. A graduated narrow stem converts this volume change into a visible immersion reading. Calibration specifies a temperature and a reading convention for the meniscus; a sample's other dissolved substances can affect a sugar-calibrated scale.

A hydrometer used to measure the specific gravity of urine; readings depend on calibration conditions and dissolved constituents.

A submerged lifting surface that develops lift when water moves relative to it. A thin two-dimensional model relates its motion-induced lift per unit span to bound circulation; buoyancy is a separate contribution to the upward fluid force.

An anelastic model filters acoustic density fluctuations while retaining a prescribed stratified reference mass density. For ${\rho }_0\propto e^{-z/H}$, continuity is $u_x+v_y+w_z-w/H=0$. Momentum and buoyancy conventions must be consistent: with reduced pressure $\Phi =p^{\prime }/{\rho }_0$, a pressure acceleration $-\nabla \Phi$ contains a vertical density-gradient term that is absent from $-{\rho }_0^{-1}\nabla p^{\prime }$. The two expressions cannot be interchanged when deriving waves in a stratified atmosphere.

A device that transfers air into or out of a vessel. Historical vacuum experiments used a receiver, valves and a mechanical pump; leakage limited the attainable pressure.

A coherent or fluctuating swirling region of a velocity field. Eddies transport momentum, heat and material constituents, and their averaged fluxes can modify a mean current. They need not be nonlinear: a linear wave packet can also have nonzero averaged quadratic fluxes.

Relative tangential motion between material and its solid boundary. In rheometry it can produce an apparent shear rate greater than the actual bulk rate and distort inferred viscosity or modulus.

A liquid in steady rigid rotation has a free surface $z-z_0={\Omega }^2{r}^2/(2g)$. Comparing with a paraboloid $z-z_0=r^2/(4f)$ gives $f=g/(2{\Omega }^2)$. The result follows from hydrostatic equilibrium in the rotating frame and is independent of density.

A device that supplies mechanical energy to a fluid to drive a flow. A specified pump installation does not by itself determine the resulting mass flow rate; the connected system's resistance also matters.

A rotating-impeller pump that transfers mechanical energy to a fluid. At a fixed speed its head varies with flow; its characteristic curve and the connected system determine the pump operating point.

The flow and head at which a pump's characteristic curve intersects the hydraulic system's required-head curve. Increasing system resistance at unchanged speed normally reduces flow on a conventional falling pump curve.^[ref]

Four times a flow passage's cross-sectional area divided by its wetted perimeter. It provides a characteristic length for many internal-flow correlations, with applicability depending on geometry and flow regime.

Numerical solution of fluid-flow models. In fossil functional analysis, reconstructed geometry, material assumptions and boundary conditions constrain the interpretation.

Microfluidics controls and studies small fluid volumes in confined channels. In cellular assays it can impose reproducible flow and chemical fields, but channel geometry, wall interactions and flow-induced stresses can themselves alter behaviour.

A nanopore is a channel or aperture with a nanometre-scale constriction. A solid-state realization can connect conducting liquid reservoirs, permitting ionic current to report particle passage through resistive-pulse sensing. A geometrical pore diameter alone does not determine its conductance: axial shape, reservoir access, surface charge, solution conductivity and interfacial effects can all matter. This definition concerns a physical aperture, without a biological application.

In a one-dimensional Ohmic approximation, a linear conical frustum of length $L$, end radii $a,b>0$ and ionic resistivity $\rho$ has $r(z)=a+(b-a)z/L$. Conservation of electric current gives $d\mathcal{R}=\rho dz/[\pi r(z{)}^2]$. Integration yields the displayed resistance. Two identical opposed frusta, each of length $l$, have resistance $2\rho l/(\pi ab)$, whether their central junction is the wider or narrower end. Their local field profiles nevertheless differ. The approximation neglects surface conduction and resistance outside the defined pore, and is not an exact three-dimensional field solution for arbitrary wall slope.

For a point-like object with constant electrophoretic mobility, deterministic traversal is the integral of inverse local drift speed. In a conical frustum of length $L$ at pore voltage $V$, the one-dimensional field is $E(z)=Vab/[Lr(z{)}^2]$. Integrating gives $t=L^2(b/a+1+a/b)/(3\mu V)$ for positive $\mu V$. For two symmetric frusta of total length $2l$ and the same total pore voltage, replace $L$ by $2l$. The inequality $b/a+a/b\ge 2$ shows that this nonuniform field cannot produce a shorter deterministic point-particle transit than a uniform field at fixed total length and voltage. Reservoir forces, diffusion, adsorption and finite-object conformation are additional physics.

Electrical detection of particles passing through a conducting-liquid aperture. A particle changes the effective pore conductance transiently; pulse amplitude and duration depend on geometry, transport and interfacial charge, so calibration is essential.

A small insulating object of excluded volume $\Omega$ perturbs an Ohmic pore most strongly where the pore cross-sectional area $A$ is small. In a local weak-blockage model, replacing a small conducting volume gives the displayed resistance scaling, with a shape-dependent coefficient. At fixed voltage, $\Delta I\simeq -I_0\Delta \mathcal{R}/{\mathcal{R}}_0$. A symmetric pore with narrow ends and a wide middle can therefore resolve two end features in a passage pulse. This is a conditional volume-exclusion model: charged-particle ion transport, surface interactions and bandwidth can change the observed sign, amplitude and shape.

The Magnus effect is a transverse force on a spinning body in relative motion through a fluid. Rotation changes the surrounding circulation, pressure distribution and wake. In a two-dimensional cylinder model, the Kutta–Joukowski theorem gives force per unit span ${\mathbf{F}}^{\prime }=\rho \mathbf{U}\times (\Gamma \hat{\mathbf{z}})$, with relative fluid velocity $\mathbf{U}$ and circulation $\Gamma$ positive anticlockwise viewed from $+z$. Real viscosity and boundary-layer behavior help determine the circulation; impermeability alone does not determine it from spin rate.

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.

In his 1883 paper, Osborne Reynolds investigated how the character of pipe flow depends jointly on velocity, pipe diameter and viscosity. His final dye apparatus drew water from a large glass tank through a horizontal glass tube with a smooth trumpet-shaped entrance. A separate reservoir supplied strongly colored water just upstream of the entrance; a valve below the tank controlled the outflow. The glass tube was immersed in the tank, so its internal dye band could be watched. Reynolds repeated the trials with different tube diameters and water temperatures. ^[ref]

The volumetric flow rate was measured independently of the appearance of the dye. Reynolds timed a one-inch fall of the tank's water level, corresponding to a calibrated discharged volume. In modern notation, if this volume is $V$ and the elapsed time is $\Delta t$, then $Q=V/\Delta t$ and $U=4V/(\pi D^2\Delta t)$. He determined the glass tubes' internal diameters by filling them with water. For a uniform bore, a known filled length $\ell$ and measured water volume $V_{\mathrm{tube}}$ determine $D=\sqrt{4V_{\mathrm{tube}}/(\pi \ell )}$; checking bore uniformity matters because the same volume would only give an average cross-section in a tapered tube. Thermometers at different tank levels checked both the temperature and its uniformity. Reynolds compared the viscosity at different temperatures using Poiseuille's results and checked his own water by separate experiments of that type. Thus the Reynolds number estimate did not come from judging the dye band or from assuming that transition occurred at a preset value. ^[ref]

A practical repetition can instead collect a known outflow volume in a measuring vessel and time its collection, or weigh the collected water and divide its mass by mass density. Measure the internal diameter, rather than the external diameter, and obtain kinematic viscosity at the measured water temperature from independently established properties or a separate viscosity measurement. Record the dye appearance and observation position alongside $Q$, $D$, temperature and the resulting Reynolds number. These are alternative measurement procedures following $Q=V/\Delta t$, not additional historical claims about Reynolds' dye apparatus.

Begin with settled water and a small dye supply, then open the outlet valve in small steps. Reynolds saw a thin straight band at sufficiently low speeds. With a tank that had not fully settled, the band could wander while retaining its coherence. At higher discharge it could abruptly spread into the surrounding water downstream of the entrance. Short illumination by an electric spark resolved curls within what otherwise looked like a colored cloud. This connected the visible breakup with internal irregular motion rather than merely the appearance of an opaque pipe. ^[ref]

The controls are part of the result. Reynolds used a smooth entrance and let the tank rest to reduce inlet disturbances. He found that a heated tank could retain slow circulation as it cooled, so measuring one temperature was insufficient: thermal uniformity and residual motion also mattered. His wire-coil trial broke up the band below the previous onset speed, but he explicitly cautioned that once the dye was mixed he could not tell from that dye alone whether the motion became orderly again farther along the tube. In a repetition, keep the dye supply small enough to act approximately as a passive scalar, avoid a strong injection jet, and record the inlet condition and downstream observation distance. If also testing the Hagen-Poiseuille equation, allow the profile to develop and measure pressure drop between positions within the developed region; do not treat the entrance and outlet losses as part of its uniform pressure gradient. The dye experiment by itself does not establish fully developed flow. ^[ref]

The experiment's usefulness is that it compares a visible change in the motion with independently measured scales. Repeating it with several diameters and temperatures, while controlling disturbances, tests the organizing role of the Reynolds number. Its limitations are equally instructive: dye is a tracer of transport, diffusion still acts in laminar flow, and a local dye observation cannot by itself prove that turbulence persists throughout an arbitrarily long pipe.

The void fraction is the fraction of a mixture's local volume occupied by gas. For a horizontally uniform bubbly liquid it is the transported concentration in a one-dimensional gas-volume continuity equation. It is distinct from gas mass fraction because gas and liquid have different mass densities.

For steady upward bubbles in a stationary liquid, assume fixed bubble mass, quadratic drag force, negligible gas weight and inertia, and negligible wall force. Then $a\propto {\rho }_g^{-1/3}$. Hydrostatic mixture support and terminal drag force balance give $v^2\propto a(1-\alpha )$. The gas mass flux $J=\alpha {\rho }_{g}v$ is constant, hence $J^2\propto {\alpha }^2(1-\alpha ){\rho }_g^{5/3}$. Comparing at a reference point proves the displayed relation. The approximation requires a small gas/liquid mass density ratio and macroscopic variation scales; it need not hold near a gas-dominated endpoint or rapid acceleration region.

If the vessel radius is proportional to $e^{-\beta z}$, its area is proportional to $e^{-2\beta z}$. Gas-volume conservation gives the displayed equation. Along concentration characteristic curves, $\dot{z}=V_s(1-2\varphi )$ and $\dot{\varphi }=2\beta V_s\varphi (1-\varphi )$. Uniform initial concentration stays uniform in the uninterrupted bubbly region and follows this logistic law. Individual bubbles move at $V_s(1-\varphi )$ instead: their rise distance there is $\eta =-(2\beta {)}^{-1}\log [{\varphi }_0+(1-{\varphi }_0)e^{-2\beta V_{s}t}]$. Boundary interactions and dense-mixture backflow can invalidate the simplified closure.

For initially uniform bubbles with upward speed $V_s(1-\varphi )$, no new gas production and free escape at the top, gas conservation gives flux $V_s\varphi (1-\varphi )$. A clearing shock rises from the impermeable bottom with speed $V_s(1-{\varphi }_0)$ by the Rankine-Hugoniot condition. The displayed mean void fraction averages the remaining uniform region over the whole liquid depth $H$. This assumes the dilute regime and no retained foam layer.

For ${\alpha }_t+[V_s\alpha (1-\alpha ){]}_x=0$ in a finite upward pipe, bubbles move upward at $V_s(1-\alpha )$, but dense-mixture phase-fraction characteristics move downward. Hence an open upper end may be an incoming boundary for the scalar conservation law even though the bubbles themselves leave through it. Atmospheric pressure alone does not prescribe that incoming phase-fraction trace.

The bottom-clearing shock from zero to a uniform initial ${\alpha }_0$ has speed $V_s(1-{\alpha }_0)$. Continuing that constant right state to the outlet gives clearance time $L/[V_s(1-{\alpha }_0)]$, but at ${\alpha }_0>1/2$ it requires a compatible maintained incoming trace. For ${\alpha }_0=3/4$, the shock is $x=V_{s}t/4$, and a right-state characteristic traced backward from it has initial position $3V_{s}t/4$. It reaches the upper boundary at time $4L/(3V_s)$, before the formal clearance time $4L/V_s$. Later evolution therefore depends on the top closure. This boundary issue is distinct from the validity of the local slip and conservation derivation.

The Strouhal number compares the advective time $h/U$ with a prescribed forcing time $1/\omega$. Some conventions use ordinary frequency instead of angular frequency, introducing a factor $2\pi$. For a longitudinal disturbance extending over length $\lambda h$, the relevant ratio is $\lambda St$ rather than $St$ alone.

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 signed lift divided by dynamic pressure times a stated reference area. For a two-dimensional airfoil, use $c_l=L^{\prime }/(q_{\infty }c)$ instead. Reference area and experimental flow conditions are part of the definition of a reported coefficient; see lift for a worked scale and a pressure/balance experiment.

Dynamic pressure is the bulk kinetic energy per unit volume associated with fluid motion at speed $U$ in a specified reference frame:Here $\rho =m/V$ is mass density. The word "pressure" identifies its units and its use as a pressure scale; $q$ is not an extra isotropic stress to add to the local static pressure in every surface traction. A stationary fluid can have large static pressure with $q=0$. Conversely, knowing $U$ and $\rho$ fixes $q$ without fixing static pressure. The definition remains useful in viscous fluid flow and compressible flow, even when the simplest Bernoulli equation is inapplicable. ^[ref]

For stipulated mass density ${\rho }_{\infty }=1.20\mathrm{kg}{\mathrm{m}}^{-3}$ and reference speed $U_{\infty }=20\mathrm{m}{\mathrm{s}}^{-1}$, $q_{\infty }=240\mathrm{Pa}$. At $40\mathrm{m}{\mathrm{s}}^{-1}$ with the same mass density, $q_{\infty }=960\mathrm{Pa}$. This quadratic scaling makes dynamic pressure a useful normalization for the lift coefficient, drag coefficient and pressure coefficient; it predicts four times a particular force only if its corresponding coefficient remains unchanged.

The connection with stagnation pressure requires hypotheses. For steady flow, constant mass density, inviscid flow, and gravity as the only body force, the Bernoulli equation gives, along one streamline,Here $p$ is static pressure, $g$ is gravitational acceleration, and $z$ is height. Bring that flow to rest without losses at the same height and the limiting static pressure is the stagnation pressure $p_0$:There must be no intervening pump supplying work or other energy input. Constancy along one streamline does not assert equal $B$ on all streamlines; irrotational flow in a connected region supplies that additional conclusion. Changing height requires retaining $\rho gz$. In compressible flow the stopping pressure difference generally differs from $\rho U^2/2$, and viscous dissipation can reduce recovered stagnation pressure. ^[ref]

A Pitot-static tube turns this relation into a measurement: its upstream-facing opening senses stagnation pressure, while separate side openings sample static pressure. Connect these to the positive and negative ports of a differential pressure transducer or manometer. Under the hypotheses above, the measured difference $\Delta p$ givesThus the stipulated $240\mathrm{Pa}$ and $1.20\mathrm{kg}{\mathrm{m}}^{-3}$ imply $20\mathrm{m}{\mathrm{s}}^{-1}$, not an independently measured speed. Align the probe with the incoming flow and keep its passages clear: blockage and pinched tubing invalidate the comparison. At very low speed, the difference can be smaller than instrument error; at high Mach number, compressibility and possibly a shock wave require a different reduction. ^[ref]

A useful control is differential-pressure calibration. For a digital manometer, zeroing establishes the reading at equal port pressures; it does not establish the gain away from zero. The manufacturer's Series 475 Handheld Digital Manometer: Specifications - Installation and Operating Instructions by Dwyer Instruments, Inc. directs the user to vent both ports to the atmosphere before zeroing. Its stated accuracy is a fraction of full scale, not of the reading, and depends on instrument temperature. Select a suitable range and record the instrument variant and units; the manual covers several ranges.^[ref]

An independent physical check can use a liquid-column manometer: apply a stable pressure difference to both instruments and compare the digital reading with the hydrostatic pressure inferred from the vertical level difference. For equal-height connections carrying a gas of mass density ${\rho }_g$ above a manometric liquid of mass density ${\rho }_m$,This follows by balancing hydrostatic pressure down and up the two connected columns. Use vertical height, even if a tube is inclined. Known mass densities, measured height and gravitational acceleration set the comparison scale. Check several differences across the working range and reverse the port connections to check the sign. This is a proposed independent check, not a claim that the cited teaching experiment performed a full instrument calibration. Recheck zero after reconnecting the sealed tubing; a zero adjustment alone cannot detect a blocked passage or an incorrect scale.

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$.

A small flush opening communicating the local surface static pressure through tubing to a manometer or pressure transducer. An upstream-facing Pitot tube instead brings incident flow to rest to measure stagnation pressure.

The local fluid pressure, distinguished from dynamic pressure and stagnation pressure. "Static" does not require the fluid to be stationary.

Pressure referred to an ideal vacuum, rather than to atmospheric pressure or another chosen reference. Gauge pressure instead subtracts atmospheric pressure.

The pressure obtained by bringing a flow to rest without losses; in compressible flow the usual definition specifies an isentropic process for stopping. For the constant-mass density, same-height Bernoulli equation discussed under dynamic pressure, $p_0=p+q$.

An upstream-facing probe that senses stagnation pressure. Determining speed also requires a separate static pressure and mass density.

A combined Pitot tube and static pressure probe with separate passages for the two pressures.

An instrument for measuring pressure relative to another pressure. A liquid-column manometer uses hydrostatic pressure; a digital manometer uses a pressure transducer.

The appeal of a liquid-column manometer is that a pressure measurement becomes a visible length measurement: mass density and gravitational acceleration supply the conversion. It can compare two pressure tappings in a flowing fluid, even though the measuring liquid itself must settle into hydrostatic equilibrium. The reference connection determines which pressure difference is being measured.

A U-tube differential manometer contains a connected column of liquid with one gas connection above each liquid meniscus. Write $p_L,p_R$ for the left and right pressures, and $z_L,z_R$ for the meniscus heights on an upward vertical coordinate. Assume constant mass density $\rho$ and downward gravitational acceleration of magnitude $g$; neglect the gas column's mass density and differences in capillary pressure. The signed vertical length difference $\Delta h=z_R-z_L$ then givesLarger left pressure means a lower left liquid meniscus: $p_L>p_R$ implies $z_L<z_R$. Equal pressures give equal meniscus heights.

The balance follows directly from hydrostatic pressure. At a common height $z_0$ below both surfaces, the two points in the connected homogeneous liquid have equal pressure:Subtracting cancels the common bottom depth and leaves $p_L-p_R=\rho g(z_R-z_L)$. The same-height rule follows from hydrostatic equilibrium, $\nabla p=-\rho g{\mathbf{e}}_z$: a horizontal displacement changes no pressure, and for constant $\rho$ the quantity $p+\rho gz$ is constant throughout the connected liquid. Equal height alone is insufficient for unrelated vessels or for fluid flow with significant acceleration.