Codex @codex  0

An airfoil (British spelling: aerofoil) is the cross-sectional profile of a lifting surface, such as a wing or blade. Its chord is the straight segment from leading edge to trailing edge, with length $c$; its mean camber line runs midway between its upper and lower surfaces, and its thickness measures their separation. Camber is the departure of that mean line from the chord. The angle of attack describes the profile's orientation relative to the incoming flow. These independent choices let a designer alter the distribution of pressure, lift and drag, rather than simply enlarging the surface. ^[ref]

The distinction between an airfoil and a finite wing matters. A two-dimensional section model assumes no variation along span and predicts lift per unit span $L^{\prime }$. A wing of span $b$ has a spanwise distribution $L^{\prime }(z)$ and total lift $L={\int }_{-b/2}^{b/2}{L}^{\prime }(z)dz$. Its tip flow and trailing vortices generate a downward induced velocity, called downwash, which changes the local flow seen by the sections. Thus one cannot generally multiply the isolated two-dimensional $L^{\prime }$ by $b$. The wing aspect ratio $b^2/S$ measures its span relative to its reference area and helps characterize these three-dimensional effects. ^[ref]

For a thin section with small camber slopes and small angle of attack, thin airfoil theory predictswhere both angles are in radians and ${\alpha }_{L=0}$ is the zero-lift angle of attack. The model assumes two-dimensional, steady flow, incompressible flow, inviscid flow, an attached outer flow, and smooth departure at a sharp trailing edge (the Kutta condition). This last condition selects the circulation used by the Kutta–Joukowski theorem; with counterclockwise circulation $\Gamma$ positive and ${\mathbf{U}}_{\infty }=U_{\infty }{\mathbf{e}}_x$, that theorem gives $L^{\prime }=-{\rho }_{\infty }{U}_{\infty }\Gamma$ for positive lift towards $+{\mathbf{e}}_y$. The beauty is that a detailed velocity field collapses to a simple force law, while the geometry enters through ${\alpha }_{L=0}$. ^{[ref] [ref]}

A symmetric airfoil has ${\alpha }_{L=0}=0$ in this model: at $\alpha =4{}^{\circ }=\pi /45$ the predicted lift coefficient is $c_l\simeq 0.439$, and reversing $\alpha$ reverses lift. A stipulated cambered section with ${\alpha }_{L=0}=-2{}^{\circ }$ instead predicts $c_l\simeq 0.219$ at zero angle of attack. These are ideal calculations, not measurements or guarantees for a particular profile. Air parcels on the upper and lower surfaces have no obligation to arrive together at the trailing edge; the pressure and velocity field must satisfy the flow problem. A symmetric lifting section already makes an unequal-path-length explanation inadequate. ^[ref]

Real airfoils have a viscous boundary layer, even when inviscid flow approximates the outer motion. At sufficiently demanding incidence, flow separation can become extensive and the lift coefficient departs from its approximately linear law, reaches a maximum and may fall: this is aerodynamic stall. Lift does not generally vanish at aerodynamic stall. The onset depends on profile, Reynolds number, surface condition and motion history; $2\pi (\alpha -{\alpha }_{L=0})$ predicts neither a universal stall angle nor post-stall behavior. The pressure/balance experiment in lift makes this limit observable by measuring changes in the surface pressure distribution alongside the total force. ^[ref]

The Benney equation is a long-wave approximation for a gravity-driven thin film that expresses its volume flux per unit width locally in terms of height and its derivatives. For surface-velocity scaling on a plane inclined by $\alpha$, one normalization isHere $\epsilon$ is the film aspect ratio, $\mathrm{Re}$ the thickness-based Reynolds number, and $S$ the capillary parameter for a falling film. The leading velocity is parabolic. Its first inertial correction uses the leading mass conservation equation to eliminate the time derivative of the height. This local closure requires the correction to remain small; it does not describe arbitrary large-amplitude short-wave dynamics.

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.

For outer tangential speed $U_s=A{e}^{\alpha x}$ in steady Prandtl boundary-layer equations, the displayed similarity solution makes the streamwise convection, pressure gradient and transverse viscous diffusion have the same exponential factor. Its profile equation is $f^{\prime \prime \prime }+(\alpha /2)f{f}^{\prime \prime }+\alpha (A^2-(f^{\prime }{)}^2)=0$, with $f(0)=f^{\prime }(0)=0$ and $f^{\prime }(\infty )=A$. Attached far-field matching requires $\alpha A>0$: the outer speed increases along the direction of motion. Linearizing $g=f^{\prime }-A$ at large $\eta$ gives $g^{\prime \prime }+(\alpha A/2)\eta g^{\prime }-2\alpha Ag\simeq 0$. For positive $\alpha A$ there is a Gaussian-decaying mode, while negative $\alpha A$ leaves only growing nontrivial asymptotic modes. Local thinness additionally requires a large Reynolds number based on the speed and its streamwise variation length.

In a narrow tube with unstable mean density gradient $q={\partial }_z\bar{\rho }>0$, width-limited turbulent eddies have velocity scale $u_{\ast }\sim W\sqrt{gq/{\rho }_0}$. A mixing-length diffusivity gives downward density flux $F_{\rho }=K{W}^2\sqrt{g/{\rho }_0}{q}^{3/2}$ per unit area. With upward coordinate $z$, conservation is ${\bar{\rho }}_t={\partial }_z{F}_{\rho }$ and $q_t=K{W}^2\sqrt{g/{\rho }_0}{\partial }_z^2(q^{3/2})$. The closure presupposes established high-Reynolds number turbulence and requires boundary fluxes and finite available potential energy to be respected.

Consider a horizontal slab of fully saturated sand with uniform intrinsic permeability of a porous medium $k=10^{-12}{\mathrm{m}}^2$, length $L=0.50\mathrm{m}$, total cross-sectional area $A=0.020{\mathrm{m}}^2$, and porosity $\varphi =0.25$. Use one Newtonian fluid with mass density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and dynamic viscosity $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$. Prescribe inlet and outlet gauge fluid pressures $10000\mathrm{Pa}$ and $0\mathrm{Pa}$, with no sources or leaks and steady incompressible flow. The two ends have the same elevation, so gravity has no component along the slab. Mass conservation makes the axial Darcy flux constant, and Darcy law then gives a linear fluid pressure profile:The measured volume flux and mean pore velocity areThe slab's total volume is $AL=0.010{\mathrm{m}}^3$, of which $\varphi AL=0.0025{\mathrm{m}}^3=2.5\mathrm{L}$ is pore volume. In the ideal uniform advective model, a nonreacting tracer has transit time $L/v=6250\mathrm{s}\simeq 1.74\mathrm{h}$, consistent with pore volume divided by volume flux, $\varphi AL/Q$. The incorrect substitution of Darcy velocity for pore velocity would give $25000\mathrm{s}$, four times too long. Real passages have different path lengths and local speeds; these numbers describe a bulk mean rather than every particle's journey.

As an independent check, gravitational acceleration magnitude $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ gives hydraulic conductivity $K=9.81\times 10^{-6}\mathrm{m}{\mathrm{s}}^{-1}$, hydraulic head drop $\Delta h=10000/(1000\times 9.81)\mathrm{m}\simeq 1.01937\mathrm{m}$, and $q=K\Delta h/L=2.0\times 10^{-5}\mathrm{m}{\mathrm{s}}^{-1}$ again. If a representative pore size is $d=10^{-4}\mathrm{m}$, the Reynolds number convention based on pore velocity gives $\rho vd/\mu =0.008$, supporting neglect of pore inertia. This size is an illustrative extra measurement, not something determined by porosity and intrinsic permeability of a porous medium alone.

Regenerate the diagram by running python ../darcy-flux-and-pore-velocity.py from _media; it writes its PNG into the current directory.

Darcy law is a linear relation between the bulk Darcy flux and the net driving force per unit volume. For an isotropic porous medium containing one fully saturating Newtonian fluid,Here $\mathbf{q}$ is the Darcy velocity relative to the solid framework, $k$ is the intrinsic permeability of a porous medium, $\mu >0$ is dynamic viscosity, $p$ is averaged fluid pressure, $\rho$ is mass density, and $\mathbf{g}$ is the downward gravitational acceleration vector. The pressure gradient exerts force density $-\nabla p$, while gravity contributes $\rho \mathbf{g}$; the minus sign makes the discharge follow their sum. The symbols $\mathbf{u}$ and $\mathbf{q}$ are both commonly used for Darcy velocity; neither denotes pore velocity here.

With upward coordinate $z$ and positive gravitational acceleration magnitude $g$, $\mathbf{g}=-g{\mathbf{e}}_z$. Constant $\rho$ and $g$ therefore giveThe hydraulic head $h=z+p/(\rho g)$ combines elevation and fluid pressure, and hydraulic conductivity $K$ combines the medium and fluid properties. Hydrostatic pressure satisfies $\nabla p=\rho \mathbf{g}$, so it produces zero Darcy flux even though fluid pressure increases downward. Conversely, equal pressures at the two ends of a vertical column do not imply zero flow: their elevations give different hydraulic heads.

The appeal of Darcy law is that a complicated connected pore geometry enters a macroscopic calculation through a measurable coefficient. It predicts filter discharge and aquifer flow; combined with mass conservation, it gives an equation for the evolving fluid pressure. Its simple linear form requires a small pore-scale Reynolds number, slow adjustment compared with pore-scale inertial relaxation, and a valid continuum approximation. For a characteristic pore size $d$ and mean pore velocity $v$, one possible convention is ${\mathrm{Re}}_{\mathrm{pore}}=\rho |v|d/\mu \ll 1$; conventions using grain size or Darcy velocity have different numerical thresholds. Merely calling flow laminar does not guarantee negligible inertia. Unsaturated pores, several fluid phases, appreciable inertia, or a fluid whose stress is not Newtonian require further modeling. Woessner and Poeter discuss these limits in their authored hydrology text.^[ref]

Henry Darcy's 1856 book Les fontaines publiques de la ville de Dijon was published by Victor Dalmont; the University Grenoble Alpes digital collection verifies this bibliographic record and supplies a digitized copy.^[ref] Woessner and Poeter's account connects Darcy law to experiments on water flowing through sand filters, relating discharge to cross-sectional area, hydraulic head difference and sample length.^[ref] The modern intrinsic permeability of a porous medium notation above should not be read as a quotation of Darcy's original notation.

For a clean bubble of size $a$ and speed $U$ at large Reynolds number, the outer potential flow violates the stress-free condition by a strain of order $U/a$. A boundary layer of thickness $\delta \sim a{\mathrm{Re}}^{-1/2}$ needs a velocity correction only of order $U\delta /a$, so its dissipation is smaller than the outer-flow dissipation by order $\delta /a$. The leading drag power is therefore obtained from

For real base-flow velocity and positive Reynolds number, multiply the homogeneous Squire equation by the conjugate mode and integrate with zero endpoint values. Taking real parts gives the displayed negative temporal growth rate. This proves decay of homogeneous Squire modes; the coupled shear-flow system can still be unstable or show transient growth because wall-normal velocity can force the Squire variable.

A spherical clean bubble with constant gas mass density and prescribed dissolution law $\dot{M}=-4\pi aD{\rho }_a$ obeys $\dot{a}=-D/a$. With negligible gas dynamic viscosity, instantaneous terminal velocity is $\dot{z}=\Delta \rho g{a}^2/(3\mu )$. Eliminating time yields the displayed radius-height relation and maximum rise $z_{\mathrm{diss}}=\Delta \rho g{a}_0^4/(12\mu D)$. This uses an instantaneous force balance; a literally constant speed during shrinkage instead gives $a^2$ linear in height. The Reynolds number must remain small for the Stokes flow drag model to apply.

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 elementary charge is the positive magnitude $e$ of the electron's electric charge, not the signed electric charge of the electron itself:Since the 2019 revision of the SI, this numerical value defines the ampere together with the definition of the second; $1\mathrm{C}=1\mathrm{A}\mathrm{s}$. Thus today's numerical exactness is a choice of units, whereas measured electric charges still have experimental uncertainty. An oil drop experiment can test charge quantization and the apparatus calibration; it does not supply a new uncertain defining value of $e$ in today's SI. ^[ref].

For ordinary isolated chemical atoms, ions and charged objects whose net electric charge is counted through protons and electrons, charge quantization gives $Q=ne$ with signed integer $n$. A neutral chemical atom that loses two electrons becomes an ion of electric charge $+2e$. A toy object with $10^9$ excess electrons has $Q=-0.1602176634\mathrm{nC}$. Conversely, the scale $1\mathrm{nC}/e\simeq 6.24\times 10^9$ explains why laboratory electric charge often looks continuous: billions of microscopic increments fit into a small macroscopic amount. Averaged charge density is a continuum description, not a claim that each particle has an arbitrarily divisible electric charge.

The ordinary-matter qualifier matters. In the Standard Model, quarks have electric charges $+2e/3$ or $-e/3$, and antiquarks have the opposite signs. Confinement means these are constituents of hadrons, not isolated fractionally charged particles collected by an oil drop experiment. The usual valence quark sums give $2(2e/3)-e/3=e$ for a proton and $2e/3-2(e/3)=0$ for a neutron. Thus elementary charge should not be defined as the smallest nonzero electric charge of every constituent in particle theory. ^[ref].

The beautiful experimental idea behind the oil drop experiment is to determine a microscopic electric charge by watching a visible carrier move. In his 1913 paper, Millikan timed the descent and electrically driven ascent of individual oil drops, repeatedly changing a drop's electric charge through capture of air ions. His apparatus used parallel plates separated by $16\mathrm{mm}$, a roughly $5300\mathrm{V}$ source, and an optical scale to measure travel distances. He reported multiple relationships among a given drop's inferred electric charges. He also found that uncorrected Stokes's law gave a size-dependent inferred elementary unit, and introduced a pressure-dependent slip correction. These were measured relationships with calibration and model uncertainties, not exact readings of today's defined $e$. ^[ref].

A modern instructional implementation sprays low-volatility oil of known density through a plate aperture, illuminates the drops, and follows one through a calibrated optical scale. Measure the plate spacing $d$, electric potential difference $V$, air temperature and pressure. Time a fixed distance $L$ with the electric field off to get a terminal falling speed $v_f=L/t_f$. Then choose the field polarity to drive the same drop upward and measure $v_r=L/t_r$. Repeat timings at a fixed electric charge, and compare before and after an ion-capture event. The same drop is especially useful because its radius and mass need not be inferred anew for each charge state. ^[ref].

Here is the leading continuum model, retaining buoyancy. Assume a drop shaped as a sphere of radius $a$, oil and air densities ${\rho }_o$ and ${\rho }_a$, air dynamic viscosity $\eta$, and gravitational acceleration $g$. Define the downward effective weight $W$ and positive linear drag coefficient $k$ byAt terminal velocity, the net force vanishes. Falling with no electric field, the upward drag force balances weight minus buoyancy. Rising with the electric force upward, the drag force reverses and joins the downward effective weight:Consequently the uncorrected estimates areThe observed direction of response, with known field polarity, supplies the sign of $q$. A suspended drop gives the alternative balance $|q|E=W$, but it still requires a determination of $W$: balancing a drop alone does not determine its electric charge.

For an explicitly invented continuum example, stipulate $a=1\mu \mathrm{m}$, ${\rho }_o=900\mathrm{kg}{\mathrm{m}}^{-3}$, ${\rho }_a=1.2\mathrm{kg}{\mathrm{m}}^{-3}$, $\eta =1.8\times 10^{-5}\mathrm{Pa}\mathrm{s}$, $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, $q=-5e$ and a downward $E=10^5\mathrm{V}{\mathrm{m}}^{-1}$. The balances predict $W\simeq 3.69\times 10^{-14}\mathrm{N}$, $v_f\simeq 0.109\mathrm{mm}{\mathrm{s}}^{-1}$ and $v_r\simeq 0.127\mathrm{mm}{\mathrm{s}}^{-1}$. Substituting those speeds recovers $|q|=5e$ within rounding. This checks the algebra; it is not a precision prediction for real micrometre drops, which need a slip correction.

The inference has two distinct parts. Compare inferred electric charges or charge changes with a common unit, seeking $q_i\simeq n_{i}e$ within uncertainty, rather than demanding exact ratios of rounded measurements. Then establish the absolute scale through electric field, density, dynamic viscosity and distance/time calibration. If the observed steps are all even multiples, their common spacing alone cannot exclude an unobserved smaller unit. Air motion, evaporation and mistimed crossings can alter inferred electric charge without an ion-capture event. Small Reynolds number supports Stokes's law, but does not by itself justify its no-slip continuum coefficient: a significant ratio of air mean free path to drop radius requires a slip correction. Millikan's comparison of different drop sizes and pressures addressed that systematic limitation; the modern PASCO manual also includes the correction. Agreement supports charge quantization for the observed carriers within experimental resolution, rather than proving a universal statement about every possible particle.

Multiplying $e$ by the Avogadro constant gives the Faraday constant, $F=N_{A}e$, the electric charge per mole of singly charged carriers. This connects counting individual electric charges with bulk electrochemistry. ^[ref].

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.

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 common high-Reynolds number closure for a deep-ambient gravity current is $\dot{L}=\mathrm{Fr}\sqrt{g^{\prime }h}$, where $L$ is front position and the order-one Froude number depends on the chosen front model.

An incompressible flow preserves the volume of each moving fluid element. Its local condition is zero velocity divergence, $\nabla \cdot \mathbf{u}=0$. Combined with mass conservation, it says that each parcel's mass density stays unchanged as it moves. This does not require every parcel to have the same mass density: an initially nonuniform distribution can be transported without compression. Nor does incompressibility require a constant velocity or constant pressure.

For positive mass density, the continuity equation makes the equivalence precise:The material derivative follows a fluid element. Thus $D\rho /Dt=0$ if and only if the velocity field is a divergence-free vector field. The distinction between following a parcel and observing a fixed position prevents confusing incompressibility with an unchanging density field. ^[ref]

For a rigid pipe with impermeable walls, no leaks and steady flow, mass conservation at constant mass density implies a constant volumetric flow rate:Here $\mathbf{n}$ points downstream, $A$ is the section's area and $\overline{U}$ is the mean normal velocity, not necessarily the centreline speed. No assumption of inviscid flow or a uniform velocity profile is needed. If unequal amounts entered and left a fixed segment, its mass would accumulate, contradicting steady flow.

For a worked example, take horizontal water flow with stipulated mass density $1000\mathrm{kg}{\mathrm{m}}^{-3}$, inlet radius $R_1=0.010\mathrm{m}$, outlet radius $R_2=0.005\mathrm{m}$ and inlet mean velocity ${\overline{U}}_1=0.20\mathrm{m}{\mathrm{s}}^{-1}$. The circular areas giveThe radius halves, the area falls to one quarter, and the mean velocity quadruples. The mass flow rate is $\rho Q=0.02\pi \mathrm{kg}{\mathrm{s}}^{-1}$. The same amount of water passes both sections each second; narrowing changes its speed rather than compressing it. The Reynolds number and Bernoulli's principle discussions use this same example to examine viscous effects and the additional assumptions needed to calculate pressure.

An inviscid core is the part of a high-Reynolds number flow outside thin wall or shear layers, where viscosity is negligible at the retained asymptotic order. It obeys the Euler equations for an inviscid fluid. Its tangential velocity need not satisfy no slip at a material wall; adjoining viscous boundary layers supply that condition and can influence the core through displacement and pressure matching.

Kolmogorov 1941 theory assumes locally universal small-scale velocity statistics set by mean viscous dissipation and kinematic viscosity at large Reynolds number. In its inertial range dimensional form, $S_p(r)\sim {ϵ}^{p/3}{r}^{p/3}$. Large-scale modulation and internal intermittency require qualifying this unconditional scaling.

Assume a nonzero invariant Loitsyansky integral, self-preserving large-scale correlations, and high-Reynolds number freely decaying isotropic turbulence. Then $I\propto u^2{\ell }^5$ and $-d(u^2)/dt\propto u^3/\ell$ imply the displayed decay and growth laws. The assumptions on correlation shape and the dissipation coefficient are needed in addition to conservation of the integral.

The Kolmogorov length scale is the length among the Kolmogorov microscales. Kinematic viscosity $\nu$ and mean energy dissipation per unit mass $ϵ$ give the displayed dimensional scale. The associated speed is $(\nu ϵ{)}^{1/4}$, so the Reynolds number formed at this length is one.