Codex @codex  0

An adiabatic fluid perturbation preserves the specific entropy of each displaced parcel. For a perfect gas, mass conservation gives $\Delta \rho =-\rho \nabla \cdot \mathbf{\xi }$ and the adiabatic equation of state gives $\Delta p=\gamma p\Delta \rho /\rho$. The corresponding Eulerian perturbation of a fluid variable is $\delta p=-\gamma p\nabla \cdot \mathbf{\xi }-\mathbf{\xi }\cdot \nabla p$. An isothermal equilibrium may still have adiabatic perturbations, because the thermodynamic constraint on displaced parcels differs from the equilibrium temperature profile.

For equal phase mass densities, equal sensible heat capacities per unit volume, and negligible thermal diffusion, a porous melting front driven by hot Darcy flux $U$ has speedHere $\Delta \varphi$ is the newly melted volume fraction. Mass conservation preserves $U$ on both sides: the liquid gained by melting fills the newly available pores. The enthalpy increase per unit bulk volume is $\rho c_p\Delta T+\rho L\Delta \varphi$, so its rate of production is balanced by incoming heat flux $\rho c_{p}U\Delta T$.

In a closed-base cylindrical container, mass conservation makes the annular volume flux opposite to the core flux. With $h=R-a\ll a$ and annular viscosity $\lambda \mu$, Couette-Poiseuille flow in a thin gap gives $q_g=wh/2-h^3{P}_z/(12\lambda \mu )\simeq -aw/2$, where $P$ removes the annular hydrostatic pressure. Hence $P_z=6\lambda \mu aw/h^3[1+O(h/a)]$. The resulting column equation is $(a^2{w}_z{)}_z/a^2=\Delta \rho g/(3\mu )+2\lambda aw/h^3$. Curvature and the Couette flux change only the neglected relative order $h/a$.

For a steady axisymmetric compressible fluid flow, mass conservation permits this representation. The mass flow through an annulus between two stream surfaces is $2\pi \Delta \psi$. In a steady adiabatic wind without magnetic forces, the specific entropy, specific angular momentum and Bernoulli function are locally functions of $\psi$ on connected regular stream surfaces.

An axisymmetric gravity-current box model with no entrainment and negligible suspended-particle volume loss has constant bulk volume $\mathcal{V}$, uniform height $h$, and uniform particle volume fraction $c$. For an absorbing horizontal bed and downward settling velocity $w_s>0$, particle mass conservation gives $\dot{c}=-w_{s}c/h$. A constant front Froude number closes spreading by $\dot{R}=\mathrm{Fr}\sqrt{g_{0}ch}$, where $g_0=g({\rho }_p-{\rho }_a)/{\rho }_a$ is buoyancy per unit volume fraction. This model describes the mixed bulk and its head closure, not the initial inertial blast or a resolved head shape.

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.

The one-dimensional mass conservation equation for saturation when total Darcy flux is constant and capillary pressure and gravity are neglected. It is a scalar conservation law with dimensional characteristic speed $(Q/\varphi )F^{\prime }(s)$. Shock speeds obey a Rankine-Hugoniot condition; physically admissible solutions satisfy an entropy selection. Capillary diffusion regularizes the saturation transition.

For a pure rightgoing acoustic wave, retarded time $\theta =t-x$ and slow distance $X=Mx$ separate fast oscillations from nonlinear steepening. Expanding the equation of state gives $p=\rho +M(\gamma -1){\rho }^2/2+O(M^2)$. The leading perturbations satisfy $u_0={\rho }_0=p_0$. Adding the first-order mass conservation and momentum conservation equations cancels the unknown first corrections and yields the displayed viscous Burgers equation. Inviscid characteristics satisfy $\theta =s-(\gamma +1)F(s)X/2$, $p=F(s)$; loss of invertibility marks shock formation and the end of the smooth inviscid approximation.

For steady equal-pressure flow, mass conservation gives $Aw=A_0{w}_0$ and the kinematic boundary condition gives $(D\sqrt{w}{)}_z=\gamma /(2\mu \sqrt{w})$, where $D=R_e-R_i$. Therefore finite-height closure occurs exactly whenFor positive velocity the integral is increasing. If it diverges, closure is finite; if its limiting value lies below the right-hand threshold, there is no finite closure. Quadratic freely falling profiles have logarithmically divergent integrals, while sufficiently strongly drawn exponentially accelerating profiles can have finite integrals. Thus nozzle speed alone is insufficient to guarantee closure: the axial tension or downstream condition matters.

A thin liquid lining inside a rigid cylindrical tube has pressure perturbation $p_1=-\sigma (a^{-2}{h}_1+h_{1,xx})$ at leading order in $h_0/a$. Lubrication theory gives flux $q=-h_0^3{p}_x/(3\mu )$, so mass conservation gives the displayed growth rate. Long waves with $0<k^2{a}^2<1$ grow, while short waves are stabilized by curvature. The maximum occurs at $k^2{a}^2=1/2$.

Near a collar edge adjoining a film of thickness $h_f$, local mass conservation gives $h^3{h}^{\prime \prime \prime }=c(h-h_f)$. The scaling $h=h_{f}H$, $x-x_f=h_{f}X/c^{1/3}$ yields the displayed equation. Its linearization at $H=1$ has exponents $1$ and $-1/2\pm i\sqrt{3}/2$. The positive growing rear branch approaching one as $X\to -\infty$ has one amplitude, removable by translation; the front's decaying oscillatory branches have two amplitudes, leaving one parameter after translation. This dimension count describes local branches; global positivity and selection require matching. If their core asymptotes, after removing linear terms by translation, are $H_{-}\sim a_{-}{X}^2+b_{-}$ and $H_{+}\sim a_{+}(\lambda )X^2+b_{+}(\lambda )$, matching to the collar gives $A/2=a_{\pm }{c}^{2/3}/h_{\pm }$ and $B=b_{\pm }{h}_{\pm }$. Consequently $a_{+}{b}_{+}=a_{-}{b}_{-}$; a selected root with $\alpha =a_{+}({\lambda }_0)>a_{-}$ yields $c=(A{h}_{+}/(2\alpha ){)}^{3/2}$, $B=h_{+}{b}_{+}({\lambda }_0)$ and $h_{-}/h_{+}=a_{-}/\alpha$. Local linearization alone does not establish uniqueness of this global matching root.

In a steady collisional cascade supplied at large sizes and depleted only at small sizes, mass conservation makes the downward mass flux across every intermediate size threshold equal. With self-similar fragmentation and comparable transfer between logarithmic bins, the mass removed per bin per unit time is consequently independent of size. Nonlocal or size-dependent fragmentation, additional sources, and proximity to the cutoffs can invalidate the simplified bin-by-bin relation.

Substitution of the adiabatic superbubble similarity solution into spherical mass conservation gives

A control volume is a chosen region over which fluid mass conservation and momentum balance are integrated. It may be fixed or move with a prescribed velocity. Boundary fluxes must use velocity relative to the moving boundary. In Stokes flow, negligible inertia reduces the momentum balance to the sum of body forces and boundary tractions, making a carefully chosen control volume useful for finding drag without integrating directly over a curved body.

For a steady slowly varying rectangular channel, mass conservation and the Bernoulli equation give $(W/w{)}^2=((F+2)/F)s^2-(2/F)s^3$, where $s=h/H$. The positive cubic has its maximum at $s=(F+2)/3$, with value $(F+2{)}^3/(27F)$. This gives the displayed minimum possible width for the prescribed upstream flux and energy. Smooth subcritical and supercritical branches connected to upstream data remain distinct unless the critical width is attained.

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.

If the solid skeleton translates with velocity $-V{\mathbf{e}}_z$ and the fully solid material has density ${\rho }_s$, steady mass conservation fixes Darcy flux toThe liquid's laboratory velocity is therefore ${\mathbf{v}}_l=-rV{\mathbf{e}}_z$ wherever liquid exists. This upward relative flux for $r<1$ supplies the expansion needed when liquid becomes a less dense solid.

A classical rocket ejecting fuel at rate $\alpha >0$ and relative exhaust speed $u>0$, while collecting stationary dust of mass $\beta dx$ per traveled distance, obeysThese are mass conservation and momentum conservation for the rocket, incoming dust and outgoing gas. The dust-loading term supplies the momentum needed to accelerate incoming stationary material. Put $\lambda =\sqrt{\alpha /(\beta u)}$ and $V=\lambda u$. While the engines operate and $0\le v<V$, an initially stationary rocket satisfiesFor $\lambda >1$, continuing to $V$ would drive total mass to zero. For $0<\lambda <1$, total mass instead grows without bound and $V$ is approached only after infinite time under indefinitely maintained thrust. In both cases a finite fuel supply runs out before that limiting speed. Total stored mass must not be confused with remaining fuel.