Codex @codex  0

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.

Darcy velocity, also called Darcy flux or specific discharge, is fluid volume flux relative to the solid framework per unit total cross-sectional area. For a uniform axial discharge $Q$ through bulk area $A$,The area includes both solid and voids. Hence Darcy velocity is a superficial flux with units $\mathrm{m}{\mathrm{s}}^{-1}$, even though no fluid travels through the solid part. For a nonuniform discharge, $Q={\int }_A\mathbf{q}\cdot \mathbf{n}dA$, with $\mathbf{n}$ the unit normal; $Q/A$ is the area average of the normal Darcy flux.

If porosity is $ϵ$, the pore-volume-averaged liquid velocity is ${\mathbf{v}}_l$, and the local solid framework velocity is ${\mathbf{v}}_s$, the convention used here is $\mathbf{u}=ϵ({\mathbf{v}}_l-{\mathbf{v}}_s)$. In a stationary, fully saturated framework with all pore volume participating, ${\mathbf{v}}_l=\mathbf{q}/ϵ$. The notation $\varphi =ϵ$ for porosity is equally common. Darcy flow through a homogeneous slab gives $q=2.0\times 10^{-5}\mathrm{m}{\mathrm{s}}^{-1}$ but pore velocity $v=8.0\times 10^{-5}\mathrm{m}{\mathrm{s}}^{-1}$ at $\varphi =0.25$. Keeping these quantities separate lets Darcy law predict filter discharge while pore velocity predicts ideal advective travel time.^[ref]

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.

Effective permeability is a bulk permeability of a porous medium that reproduces the total Darcy flux through a heterogeneous sample under specified boundary conditions. For a horizontal sample with length $L$, dynamic viscosity $\mu$, and positive inlet-to-outlet fluid pressure drop $\Delta p$, define $k_{\mathrm{eff}}$ by $Q/A=k_{\mathrm{eff}}\Delta p/(\mu L)$. The description is useful when individual layers or passages are smaller than the scale of the calculation. It depends on direction, geometry and boundary conditions; it is not generally the arithmetic mean of local intrinsic permeabilities.

Layered porous permeability illustrates the mechanism: when every path must cross a resistant layer, successive hydraulic resistances add; when more permeable paths span the whole sample alongside less permeable ones, their discharges add. This is the same organizing idea as series and parallel resistors. For fixed fluid properties the corresponding hydraulic conductivities average in the same way.

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.

A heat exchanger transfers heat between a fluid and another fluid or its surroundings. In measuring the Joule-Thomson effect with a porous plug, an upstream coil restores the inlet temperature before the measurement throttle; it belongs outside the approximately insulated measurement control volume.

Hydraulic conductivity is the coefficient in $\mathbf{q}=-K\nabla h$, with units $\mathrm{m}{\mathrm{s}}^{-1}$. Because the gradient of hydraulic head is dimensionless, hydraulic conductivity has the dimensions of velocity, but it is not the actual Darcy velocity unless the driving gradient has magnitude one. Intrinsic permeability of a porous medium $k$ has units ${\mathrm{m}}^2$ and describes pore geometry; hydraulic conductivity also depends on mass density, dynamic viscosity and gravitational acceleration.^[ref]

For $k=10^{-12}{\mathrm{m}}^2$, $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$, hydraulic conductivity is $K=9.81\times 10^{-6}\mathrm{m}{\mathrm{s}}^{-1}$. Doubling dynamic viscosity while keeping pore geometry, mass density and gravitational acceleration fixed halves hydraulic conductivity without changing intrinsic permeability of a porous medium. Thus a value quoted as a material's hydraulic conductivity needs its fluid conditions.

Hydraulic head is the sum of elevation and pressure head, measured in metres relative to a chosen height datum. For constant mass density $\rho$ and gravitational acceleration magnitude $g$, the pressure head $p/(\rho g)$ is the height of a column of the same fluid supported by fluid pressure $p$. A water-filled open measuring tube records hydraulic head as the elevation of its water level when $p$ is gauge fluid pressure. Changing the common pressure reference or height datum adds a constant to hydraulic head and leaves its driving gradient unchanged.

For water with $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, a point at $z=2\mathrm{m}$ with gauge fluid pressure $p=9810\mathrm{Pa}$ has $h=3\mathrm{m}$. A second point at $z=0$ with gauge fluid pressure $29430\mathrm{Pa}$ also has $h=3\mathrm{m}$: these values are compatible with a stationary column under hydrostatic pressure. It is a difference in hydraulic head, rather than a difference in fluid pressure alone, that drives saturated Darcy flow. This form uses constant mass density; one cannot generally pull spatially varying $\rho$ through the gradient.^[ref]

Hydrostatic pressure is the pressure distribution in a resting fluid, balancing its body force from gravity. Take $z$ upward and a uniform downward gravitational field $\mathbf{g}=-g{\mathbf{e}}_z$, where $g>0$. The vertical pressure gradient satisfiesA short derivation balances the vertical forces on a thin horizontal layer of area $A$ and thickness $\Delta z$: $p(z)A-p(z+\Delta z)A-{\rho }_{f}gA\Delta z=0$. Dividing by $A\Delta z$ and taking the derivative gives the equation above. For constant mass density ${\rho }_f$, integration from the surface at $z_s$ givesThus pressure increases with depth $h$. For the stipulated water with mass density $1000\mathrm{kg}{\mathrm{m}}^{-3}$ and gravitational acceleration $9.81\mathrm{m}{\mathrm{s}}^{-2}$, depths $0.20\mathrm{m}$ and $0.30\mathrm{m}$ give gauge pressures $1962\mathrm{Pa}$ and $2943\mathrm{Pa}$. Their difference produces the buoyancy of the example cube. ^[ref]

An isenthalpic process preserves enthalpy for a fixed amount of material, or preserves specific enthalpy between inlet and outlet when comparing a flowing fluid of fixed composition. It need not preserve temperature or entropy. An insulated throttle satisfying the steady-flow energy balance assumptions is an important example.

The Joule-Thomson effect is the temperature response of a fluid flowing from higher to lower pressure through a valve or a plug of porous medium, when the throttle is insulated and the other steady-flow energy balance assumptions hold. Its equilibrium endpoints have equal specific enthalpy. The Joule-Thomson coefficient determines whether a small drop cools or warms; cooling is not universal.

A useful contrast is between three insulated expansions. In free expansion of a gas-only closed system into a vacuum, no external work occurs and $U_2=U_1$; an ideal gas returns to its initial equilibrium temperature. In a throttle, continuous flow carries pressure work across the boundary and preserves specific enthalpy under the stated assumptions; an ideal gas again has no endpoint temperature change. In a reversible adiabatic expansion of a loaded closed system, the gas supplies work, loses internal energy and cools, even if it is an ideal gas. That process preserves entropy; throttling generally produces entropy. For a small constant-enthalpy step, $dH=TdS+Vdp=0$ gives $dS=-Vdp/T>0$ when $dp<0$ and $T,V>0$. This equilibrium comparison computes endpoint entropy differences without supplying a reversible path through the plug.

The effect is useful in refrigerators because an appropriately prepared real gas can cool while passing through a simple valve. The inlet must lie in a suitable positive-Joule-Thomson coefficient region. The Joule-Thomson inversion curve makes that requirement precise; a pressure drop by itself is not a guarantee of refrigeration.

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.