Codex @codex  0

Darcy flow is slow flow through a porous medium whose averaged flux obeys the Darcy law. The Darcy velocity is the fluid volume flux per total cross-sectional area; it differs from the velocity averaged only over pore space.

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.

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.

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.

Write $ϵ=1-{\rho }_s/\rho >0$ and $\Delta v=1/{\rho }_s-1/\rho =ϵ/{\rho }_s$. When the ice sphere grows by volume $4\pi a^2\dot{a}dt$, it consumes liquid volume $({\rho }_s/\rho )4\pi a^2\dot{a}dt$. The difference is expelled from the fixed cavity. Thus mass conservation fixes the outward Darcy flux in the surrounding porous medium:The Darcy law $u_r=-(\Pi /\mu ){\partial }_{r}p$ and $p(\infty )=p_m$ implyMatching at the cavity wall, whose internal pressure is uniform, givesThe solid and liquid pressures are equal when the Gibbs--Thomson relation is neglected; nevertheless the coexistence temperature changes with their common pressure because their specific volumes differ. Along coexistence, equality of chemical potentials givesHence the common-pressure melting-point slope is $d{T}_m/dp=-T_m\Delta v/L$. To first order in the pressure change,The $T_m$ multiplying this coefficient is absolute temperature, approximately $273.15\mathrm{K}$ for water; a temperature difference may equally be measured in kelvin or degrees Celsius. Substituting the numerical value $0$ from the Celsius scale into a thermodynamic denominator would be incorrect.

In the quasi-stationary thermal conduction approximation, neglect advective heat transport and use the stated effective conductivity $k$ outside the sphere. The radial Laplace equation gives the bounded solid temperature and the liquid/porous-region temperature:This treats the exterior as a uniform conducting medium, as required by the printed reduction; distinct cavity and porous-medium conductivities would give an additional thermal resistance. There is no solid-side heat flux in the bounded quasi-stationary solution. The Stefan condition therefore becomesEliminate $T_i$ and $p$ to obtainFor $x=a/R$, $t=t_0\tau$ and $t_0=R^2{\rho }_{s}L/(k\Delta T)$, this isThe pressure follows either from the Darcy flux or the Stefan condition:In particular $T_i-T_{\infty }=\Delta TK/(x+K)$: rising pressure progressively reduces the available undercooling. This is spherical freezing inside a porous cavity, whose fixed hydraulic outlet radius $R$ distinguishes it from freezing directly through a porous matrix.

For the idealized zero-radius start, integration givesA finite seed replaces the right side by its difference from the initial value. The reduction ends when the sphere fills the cavity.

For $K\gg 1$, hydraulic drainage is easy and the pressure feedback is small throughout the cavity. The leading expressions areThis is ordinary conduction-controlled spherical solidification, with a completion time asymptotic to $t_0/2$.

For $K\ll 1$, distinguish the initial crossover from the pressure-dominated regime. At $x\ll K$, equivalently $\tau \ll K^2$, the same square-root growth applies and $p-p_m\simeq p^{\ast }x/K$. At $K\ll x\le 1$, equivalently $K^2\ll \tau \lesssim 1/K$,The interface approaches $T_{\infty }$, the pressure approaches the finite coexistence value $p_m+p^{\ast }$, and growth becomes hydraulically limited. The leading completion time is $t_0/(3K)$. The cube-root approximation is not valid at the instant of nucleation, where the initial square-root regime must be retained.

The void volume within a porous medium. Porosity is its fraction of bulk volume; saturation specifies which fluids occupy that void volume.

Porous-media flow is fluid motion through the interconnected voids of a porous medium, such as a sand bed, a rock or a filter. Individual passages can be tiny and tortuous, yet their combined discharge can be predicted without resolving every passage. Darcy law supplies this bulk description when viscous resistance dominates pore-scale inertia; porosity measures how much space stores fluid, while permeability of a porous medium measures how readily that space transmits it.

The basic fully saturated model assumes one Newtonian fluid fills the connected pore space of a stationary solid framework, with incompressible flow, a small pore-scale Reynolds number, and averaging volumes large compared with individual pores but small compared with macroscopic variations. These are assumptions of the model developed below, rather than a claim that Darcy law describes every porous-media flow. Two-phase porous-media flow and poroelastic aquifers require additional constitutive and storage information.

For a first calculation, see Darcy flow through a homogeneous slab: a pressure drop determines the discharge, but porosity determines the corresponding pore velocity. Layered porous permeability then shows how geometry alone can make the bulk response anisotropic.

For a sphere of radius $a$ injecting an incompressible flow into an infinite homogeneous isotropic porous medium, Darcy's law and conservation of mass give $r^2{q}_r=C$. Integrating $p_r=-\mu C/(\Pi r^2)$ with $p(a)-p(\infty )=\Delta p$ gives $C=\Pi a\Delta p/\mu$. The hydraulic resistance is $\Delta p/Q=\mu /(4\pi \Pi a)$. A common background hydrostatic pressure cancels gravity if the imposed excess pressure is measured relative to that background.

A sharp advected temperature transition in a fluid-filled porous medium. If fluid and solid equilibrate locally in temperature, their effective heat capacity per total volume is $C_{\mathrm{eff}}=\varphi C_f+(1-\varphi )C_s$, while the advected heat flux is $C_{f}UT$ with $U$ the Darcy velocity. Neglecting conduction gives $C_{\mathrm{eff}}{T}_t+C_{f}U{T}_x=0$, so the thermal front speed is $\Gamma U$ with $\Gamma =C_f/C_{\mathrm{eff}}$. The material pore velocity is $U/\varphi$, so thermal retardation relative to that speed is $\Gamma \varphi$. Temperature-dependent viscosity couples thermal-front perturbations back to the Darcy flux.