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 Bernoulli principle relates pressure, kinetic energy per unit volume and Newtonian gravitational potential along a streamline of steady flow with inviscid flow. For constant mass density, the Bernoulli equation gives $p/\rho +u^2/2+\Phi$ constant along that streamline. Constancy across different streamlines requires further hypotheses, such as irrotational flow.

Precisely, suppose a fluid of constant positive mass density obeys the Euler equations for an inviscid fluid, its velocity field and pressure are steady, and the body force per unit mass is $-\nabla \Phi$. ThenEvery term has units of energy per unit mass. For a uniform gravitational field, Newtonian gravitational potential is $\Phi =gz$, up to an additive constant. The pressure term represents work by surrounding fluid, so acceleration can exchange pressure work with kinetic energy. Dotting the steady Euler equations with $\mathbf{u}$ gives the displayed directional derivative: this is why constancy follows a streamline. ^[ref]

Continue the horizontal water constriction example from incompressible flow. Assume negligible viscous losses, steady flow and approximately uniform flow at both measurement sections, so the streamline speeds there equal the mean speeds $0.20$ and $0.80\mathrm{m}{\mathrm{s}}^{-1}$. Equal height cancels the Newtonian gravitational potential terms. The ideal pressure drop isThis is the Venturi effect in its ideal limit: the gain in kinetic energy density is supplied by pressure work. Cambridge's treatment uses this relation to connect a measured pressure difference to volumetric flow rate. ^[ref]

The $300\mathrm{Pa}$ result is conditional, rather than a prediction of the actual viscous pipe. Mass conservation fixes $Q$ and the ratio of mean speeds even with viscosity, but nonuniform velocity profiles require an integrated energy balance, and wall friction dissipates mechanical energy. A pump would supply additional work. In particular, constant mean speed in a straight pipe does not imply constant pressure: Poiseuille flow needs a continuing pressure gradient to overcome viscosity. ^[ref]

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.

Convective acceleration is the spatial part $(\mathbf{u}\cdot \nabla )\mathbf{u}$ of the material derivative of fluid velocity. It can remain nonzero in steady flow, because a moving particle passes through positions with different velocities. The identity $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (|\mathbf{u}{|}^2/2)-\mathbf{u}\times (\nabla \times \mathbf{u})$ relates it to vorticity.

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

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

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

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

Ideal steady flow with incompressible flow, inviscid flow and irrotational flow can predict zero drag on a body that has nonzero drag in a real fluid. The zero-circulation cylinder pressure coefficient distribution gives an explicit cancellation example; viscosity and flow separation are absent from that model.

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.

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

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

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

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

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

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.

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.

Consider an original model of a cold section enclosing a throttle and collecting region, with the precooling heat exchanger outside its boundary. The inlet has mass flow rate $\dot{m}>0$ and specific enthalpy $h_{\mathrm{in}}$. The outlets are equilibrium liquid and vapor at one fixed pressure, with specific enthalpies $h_l$ and $h_g$. Let $y$ be the outlet liquid mass fraction, and let $q=\dot{Q}/\dot{m}\ge 0$ be incoming unwanted heat per unit incoming mass. Assume steady flow, no shaft work, and negligible bulk kinetic energy and gravitational potential energy changes. The steady-flow energy balance givesHere $L$ is the latent heat per unit mass at the chosen outlet pressure. The mass flow rates of the outlets are $y\dot{m}$ and $(1-y)\dot{m}$. The formula describes a flowing cold section; a filling vessel with appreciable energy or mass accumulation requires storage terms.

For a deliberately dimensionless example, choosePer unit incoming mass, one quarter leaves as liquid and three quarters as vapor. With zero incoming heat, the same inlet would give $y=0.30$. Increasing $q/L$ to $0.20$ reduces $y$ to $0.10$, and $q/L=0.30$ reaches the no-liquid boundary. These are chosen model values, not helium property data or a historical yield. In normalized enthalpy units with $h_l=0$ and $h_g=1$, conservation checks directly: $0.70+0.05=0.25\times 0+0.75\times 1$.

The admissible interval is $h_l\le h_{\mathrm{in}}+q\le h_g$, or $0\le y\le 1$. Outside it, the assumed equilibrium two-phase outlet fails: an enthalpy above $h_g$ requires a hotter single-vapor outlet under these assumptions, rather than a negative amount of liquid. Below $h_l$, a colder liquid or a different thermodynamic phase must be considered. The balance does not predict its own outlet pressure, temperature or helium property functions.

At fixed outlet pressure, hence fixed $L$, the derivatives are $\partial y/\partial h_{\mathrm{in}}=-1/L$ and $\partial y/\partial q=-1/L$. Lowering inlet specific enthalpy by regeneration and reducing incoming heat therefore buy the same increase in liquid mass fraction per unit energy per unit mass in this model. Adding any constant to all three specific enthalpies leaves $y$ unchanged. This independence of the arbitrary enthalpy reference makes the balance useful even when property tables choose different zero points.

The momentum thickness weights the streamwise velocity deficit by the local velocity carrying it:Thus $\rho U^2\theta$ measures a momentum flux deficit per unit span relative to transporting the same local mass flow rate at $U$. It is not $\rho \int (U^2-u^2)dy$, which compares a different reference flux. A Blasius boundary layer has $\theta \simeq 0.664115\sqrt{\nu x/U}$; the ratio of displacement thickness to momentum thickness is about $2.591$.^[ref]

For the steady flow over an impermeable flat plate at constant $U$, integrating the boundary-layer equation and using mass conservation givesThis connects growth of the momentum thickness to local wall shear stress. Consequently the drag per unit span on one face up to $x=L$ is $D^{\prime }=\rho U^2\theta (L)$ in the ideal Blasius boundary layer with $\theta (0)=0$. Measured velocity profiles can therefore estimate integrated drag without differentiating noisy near-wall readings. The identity assumes no wall suction and zero outer pressure gradient.^[ref]

Following a parcel. Let $\mathbf{X}(t)$ be a fluid pathline, with $\dot{\mathbf{X}}=\mathbf{u}(\mathbf{X},t)$. The multivariable chain rule gives the material derivative,The first term measures change at fixed position; the advection term accounts for the parcel moving through spatial gradients. The derivative is also called the Lagrangian derivative.

The forces on a parcel. A material volume $dV$ has mass $\rho dV$. The pressure forces on its boundary sum to $-\nabla pdV$, while the gravitational acceleration is $-\nabla \psi$. Newton's second law therefore gives the Euler equations for an inviscid fluid,This neglects viscosity and other body forces. The left side is the parcel's acceleration, not just the time derivative at a fixed spatial location.

Bernoulli conservation. The claimed streamline result needs steady flow and a barotropic fluid, $p=p(\rho )$, or an equivalent enthalpy relation along each isentropic streamline. These hypotheses are necessary; the unsteady Euler equation alone does not imply the stated conservation law. Define the barotropic enthalpy functionDot the steady Euler equation with $\mathbf{u}$, using $\mathbf{u}\cdot [(\mathbf{u}\cdot \nabla )\mathbf{u}]=\mathbf{u}\cdot \nabla (u^2/2)$. This givesThus the Bernoulli function isEvery term has units of energy per mass, ${\mathrm{m}}^2{\mathrm{s}}^{-2}$. Different streamlines can have different constants.

The spherical accretion equation. Write the inward speed as $u(r)>0$, so the velocity vector is $-u(r)\hat{\mathbf{r}}$ and the positive mass accretion rate isAssume a prescribed central potential $\psi =-GM/r$ and negligible gas self-gravity. The steady radial Euler equation isFor a barotropic equation of state, $dp/dr=c_s^{2}d\rho /dr$, with $c_s^2=dp/d\rho$. Eliminate the density derivative to obtainA smooth sonic point requires the right side to vanish where $u=c_s$, henceFor a nonisothermal gas this is an implicit relation involving the local sound speed. In the isothermal case $c_s$ is constant and this is the Bondi sonic point.

The isothermal integral. For $p=c_s^2\rho$, choose the enthalpy reference so that $h_b=c_s^2\ln (\rho /{\rho }_{\infty })$. The reservoir is at rest at infinity and $\psi (\infty )=0$, making its Bernoulli constant zero. ThusorThe logarithm has a dimensionless argument. Constant temperature during compression generally requires heat removal, so this is Isothermal Bondi accretion, rather than adiabatic accretion with $\gamma >1$.

The critical accretion rate. Select the smooth reservoir-connected transonic branch. Evaluate the last relation at $u_s=c_s$ and $GM/r_s=2c_s^2$:Mass conservation then gives the isothermal Bondi accretion rate,The units are $\mathrm{kg}{\mathrm{s}}^{-1}$. Equivalently $\dot{M}=4\pi \lambda G^2{M}^2{\rho }_{\infty }/c_s^3$, with $\lambda =e^{3/2}/4\simeq 1.12$.

The inner-boundary selection of isothermal Bondi flow is important: a standard absorbing accretor inside the sonic radius selects this critical solution. The steady equations and ${\rho }_{\infty }$ alone also admit subcritical subsonic branches supported by a different inner boundary; they do not by themselves fix a unique rate. The usual answer presumes transonic capture.

Write ${\omega }_h=\dot{\alpha }$, positive on opening. The no-slip boundary condition equates the fluid velocity with each plate's velocity:The upper condition is a moving-wall no-penetration boundary condition, not zero normal velocity in the laboratory frame. At the hinge the velocity tends to zero. With the arbitrary constant of the streamfunction chosen so that $\psi (r,0)=0$, these conditions becomeA changing angle cannot give a literally steady flow in fixed coordinates. We interpret the stated steadiness as the quasi-steady Stokes approximation at each instantaneous $\alpha$.

The printed expression is not biharmonic. To test it without reproducing the source formula, write its bracket as $G_p(\theta )$ and its streamfunction as ${\psi }_p=-\Omega r^2{G}_p/2$. For any angular function $G$, the Laplacian givesAt the concrete angle $\alpha =\pi /4$, differentiation of the printed bracket givesFor example, at the interior point $\theta =\pi /8$,The printed bracket does satisfy $G_p(0)=G_p^{\prime }(0)=G_p^{\prime }(\alpha )=0$ and $G_p(\alpha )=1$, so its boundary velocities would require $\Omega ={\omega }_h$. Satisfying these boundary conditions does not cure its nonzero biharmonic operator. The stated series is likewise not a general biharmonic basis: its $r^2\theta \sin 2\theta$ term already fails the displayed test. The separated polar solutions of the biharmonic equation contain additional radial powers; for the required $r^2$ dependence, the angular basis is $1,\theta ,\sin 2\theta ,\cos 2\theta$.

The corrected slowly opening Stokes wedge similarity solution follows by setting $\psi =-{\omega }_h{r}^{2}G(\theta )/2$. The biharmonic equation gives $G^{(4)}+4G^{\prime \prime }=0$. The two lower-wall conditions reduce its general solution toThe upper-wall condition $G^{\prime }(\alpha )=0$ fixes the ratio of $A$ and $B$, and $G(\alpha )=1$ fixes their scale. An expression that also remains finite at $\alpha =\pi /2$ isFor $0<\alpha <\pi$, its denominator is positive: the derivative of $\sin \alpha -\alpha \cos \alpha$ is $\alpha \sin \alpha >0$. Direct differentiation verifies all four boundary conditions and the biharmonic equation. The resulting velocity components are $v_r=-\dot{\alpha }r{G}^{\prime }/2$ and $v_{\theta }=\dot{\alpha }rG$.

For a fixed angle of order one, $U(r)\sim |\dot{\alpha }|r$, so the local Reynolds number is ${\mathrm{Re}}_r\sim |\dot{\alpha }|r^2/\nu$. One breakdown length is thereforeA second is the geometric length $l_2\sim \min (L,W)$, where $L$ is the distance to the radial ends and $W$ the spanwise width. At $r\sim l_2$, matching to the surrounding bath or three-dimensional end flow matters. Thus $r\ll \min (l_1,l_2)$ is required for this local two-dimensional Stokes flow, with angle-dependent factors in the inertial estimate. Ideal infinite plates remove the geometric cutoff but not the inertial one.

There is also a temporal restriction: the relevant viscous diffusion time must be short compared with the time $T_s$ over which the shape and its velocity field change appreciably. Using $r$ as a conservative diffusion length gives $r\ll \sqrt{\nu T_s}$. For order-one angles and constant $\dot{\alpha }$, $T_s\sim |\dot{\alpha }{|}^{-1}$ gives the same scale as $l_1$; narrow angles require using their actual gap and angular factors rather than extrapolating this order-one-angle estimate.

The coefficient of the growing term is fixed by the far-field complex potential: $A=\beta e^{i\alpha }$. On $z=a{e}^{i\theta }$, the no-penetration boundary condition makes the cylinder a streamline, so the imaginary part of $A{a}^2{e}^{2i\theta }+B{a}^{-2}{e}^{-2i\theta }$ must be constant. Choosing the two terms to be complex conjugates makes that constant zero, givingThis is the zero-circulation solution, also obtained by the Milne-Thomson circle theorem. Impermeability and the leading far field alone permit an additional circulation term $-i\Gamma \log z/(2\pi )$; the requested algebraic form implicitly selects $\Gamma =0$. No no-slip boundary condition is imposed on an inviscid flow.

The streamfunction in polar coordinates isIt gives $v_r(a,\theta )=0$ andThe liquid still exerts normal pressure on each surface element even though tangential viscous traction is absent. For this steady flow with zero vorticity, the Bernoulli equation has one constant $C$ throughout the connected fluid region:There is no finite uniform far-field pressure to substitute for $C$, because this straining velocity field grows with radius. Taking $\mathbf{n}=(\cos \theta ,\sin \theta )$ outward from the cylinder, the fluid force per unit axial width isIndeed $p(\theta +\pi )=p(\theta )$ while $\mathbf{n}(\theta +\pi )=-\mathbf{n}(\theta )$: diametrically opposite elements cancel exactly. Zero resultant force does not mean zero local pressure or zero deformation tendency.

Take $x$ along a horizontal gap and $s$ perpendicular to its plates, with $0<s<b$. Assume a Newtonian fluid, stationary plates, a steady flow, a laminar flow and a fully developed flow away from the entrances and side edges. The no-slip boundary condition holds on both plates. Any transverse gravity is absorbed in the hydrostatic pressure; there is no imposed transverse through-flow. Translational symmetry and the no-penetration boundary condition then give $\mathbf{u}=(u(s),0,0)$. The Navier-Stokes equation reduces toTwice integrating gives the plane Poiseuille flow:The volumetric flow rate through one gap isThere are $n-1$ internal gaps. Their total cross-sectional area is $(n-1)wb$, because the plate thickness is negligible, so the Darcy velocity isThus the permeability of a porous medium has the dimensions of squared length. If finite plate thickness were included, the Darcy velocity averaged over the whole manifold would also contain its porosity.

Both the local acceleration and the convective acceleration vanish: ${\partial }_t\mathbf{u}=0$ and $(\mathbf{u}\cdot \nabla )\mathbf{u}=u_x{\partial }_x\mathbf{u}=0$. The dominant balance is therefore between the pressure gradient and viscous stress, with no inertia contribution in this fully developed solution. This exact cancellation does not guarantee that arbitrary large-Reynolds number realizations remain laminar flows; entrances, disturbances and turbulence can invalidate the assumptions.

A small Reynolds number, $\mathrm{Re}=\rho UL/\eta \ll 1$, means that inertia is small compared with viscous forces. It can arise from small speed, small length or large dynamic viscosity. Examples are the locomotion of a microswimmer and the slow sedimentation of a small particle in a viscous liquid.

To neglect local acceleration as well, an externally imposed time scale $T$ must satisfy $\rho L^2/(\eta T)=L^2/(\nu T)\ll 1$. This holds for a steady flow or a sufficiently slowly varying quasi-steady Stokes approximation. Small Reynolds number alone is insufficient for that extra step. For an incompressible flow with constant material coefficients and uniform gravity, the Navier-Stokes equation then reduces to the Stokes equation:Taking the divergence and using commuting spatial derivatives givesTaking the curl, which annihilates the pressure gradient and uniform gravitational force, givesThese are the Harmonic pressure and vorticity in Stokes flow results.

A concrete counterexample to inferring both solely from small Reynolds number is Stokes' second problem: an oscillating plate with arbitrarily small velocity amplitude has negligible inertia, but ${\partial }_t\mathbf{\omega }=\nu {\nabla }^2\mathbf{\omega }$ is nonzero. Its vorticity is not a harmonic function. The quasi-steady qualification is therefore necessary.

Let $\mathbf{U}=(u,v)$, $\zeta =v_x-u_y$ and $D/Dt={\partial }_t+\mathbf{U}\cdot \nabla$. The inviscid shallow water equations areTaking the vertical curl of the momentum equation, with constant $f$, gives $D\zeta /Dt=-(f+\zeta )\nabla \cdot \mathbf{U}$. Combining this with the depth equation provesThus shallow-water potential vorticity is materially conserved: changes in absolute vorticity are exactly compensated by stretching or compression of the layer depth.

The relation to a Bernoulli function presupposes steady flow. In that case define $B=g(h+H)+\frac{1}{2}|\mathbf{U}{|}^2$. The vector identity $(\mathbf{U}\cdot \nabla )\mathbf{U}=\nabla (\frac{1}{2}|\mathbf{U}{|}^2)+\zeta \hat{\mathbf{z}}\times \mathbf{U}$ givesSteady volume conservation allows a transport stream function with $h\mathbf{U}=({\psi }_y,-{\psi }_x)$. Hence $\hat{\mathbf{z}}\times (h\mathbf{U})=\nabla \psi$, soIn particular, $\mathbf{U}\cdot \nabla B=0$: the Bernoulli function is constant along each streamline. Where the flow permits a single-valued dependence on the stream function, the displayed derivative describes how that constant changes between neighbouring streamlines. For an unsteady flow, shallow-water potential vorticity remains conserved but the stated steady Bernoulli function relation does not generally hold.