Darcy law ($\mathbf{q}=-k(\nabla p-\rho \mathbf{g})/\mu$)

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.

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]

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.

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 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.

In a thin fluid-filled crack, lubrication theory gives flux per unit depth $q=-h^3{p}_x/(12\mu )$. The surrounding Darcy flow has normal velocity $-k{p}_y/\mu$. Conservation in the crack gives $q_x=-[u_y{]}_{-}^{+}$, yielding the displayed condition. Pressure is continuous across the crack. This reduced model neglects transverse pressure variation and assumes wall slip lengths are small compared with $h$.

For $h=H_0(1-x^2/a^2{)}^{1/6}$, elliptic coordinates make $h^3/\sin \eta$ constant on the upper crack face. Pressure continuity and harmonicity reduce the boundary condition to $F^{\prime }(0)=\alpha F(0)$ for the cosine mode, while uniform far flow gives $F\sim -e^{\xi }/2$. Solving $F^{\prime \prime }=F$ yields the displayed pressure. The midpoint crack flux is $2aU\alpha /(1+\alpha )$; small $\alpha$ leaves the background pressure gradient nearly unchanged, while large $\alpha$ makes the crack nearly equipotential.

Subtracting the uniform-flow pressure $-\mu Ux/k$ leaves $(\mu Ua/k)\alpha e^{-\xi }\cos \eta /(1+\alpha )$. The coordinate map has $r\sim a{e}^{\xi }/2$, so $e^{-\xi }\sim a/(2r)$, giving the displayed dipole coefficient. The one-half is essential. In the perfect-conductor limit the independent expression $p=-(\mu U/k)\mathrm{Re}\sqrt{z^2-a^2}$ has the same far-field expansion.

For dipole pressure $\delta p\sim Cx/(x^2+y^2)$, choose a rectangular control volume whose transverse size tends to infinity faster than its streamwise size. Each end face contributes $\pi C$ to $\int \delta p{n}_{x}ds$, so fixed-throughflow boundary work gives $\Delta \mathcal{D}=-2\pi UC$, the displayed result. This control-surface choice removes the integrated far-field velocity correction. A circular cutoff of the pressure-only functional gives $-\pi UC$ instead and does not implement the same throughflow normalization.

For constant wall heat-flux density $F>0$, effective thermal conductivity $k_e$, buoyancy mobility $B={\rho }_{0}g\alpha K/\mu$ and thermal transport coefficient $\kappa =k_e/({\rho }_0{c}_f)$, a slender Darcy flow boundary layer has width $\delta =(\kappa k_{e}z/(BF){)}^{1/3}$. Writing $\eta =y/\delta$, $\psi =\kappa zf(\eta )/\delta$ and $T-T_a=(F\delta /k_e)f^{\prime }(\eta )$ gives the displayed equation with $f(0)=0$, $f^{\prime \prime }(0)=-1$ and $f^{\prime }(\infty )=0$. The exact integral is ${\int }_0^{\infty }(f^{\prime }{)}^{2}d\eta =1$. Consequently its kinematic Darcy momentum flux is $\int w^{2}dy=B\kappa Fz/k_e$ and its buoyancy flux is $g\alpha \kappa Fz/k_e$. These are per span on one heated face.

Assume an exponentially decreasing temperature profile while imposing wall heat flux, impermeability and the exact integrated heat budget. For $f=C(1-e^{-\eta /d})$, the wall condition gives $C=d^2$ and $\int (f^{\prime }{)}^{2}d\eta =1$ gives $d^3=2$. This yields the displayed integral approximation and discharge $\int wdy\simeq 2^{2/3}\kappa z/\delta$. It satisfies the boundary conditions and exact flux budget, but not the differential equation pointwise; its mass-flux coefficient is approximate.

For two incompressible fluids separated by a planar material front $X(t)$, the Darcy flux is uniform and the pressure drop sums the resistances of the two layers: $U=k\Delta P/[\mu X+{\mu }_h(L-X)]$. The front moves at the pore velocity $X^{\prime }=U/\varphi$. Integrating from $X(0)=0$ yields the displayed relation. The formula holds up to breakthrough $X=L$, with either sign of the viscosity difference; equal viscosities give constant speed.

If a thermal front has speed $Y^{\prime }=\Gamma U$ with $U$ the Darcy velocity, while a material front has speed $X^{\prime }=U/\varphi$, their position ratio from coincident initial locations is $Y/X=\Gamma \varphi$. Calling a retardation factor the ratio to pore velocity instead would give $Y/X=\Gamma$; these conventions must not be mixed. With cold viscosity $\mu$ and heated viscosity $b\mu$, the displacement's effective viscosity is ${\mu }_{\mathrm{eff}}=\mu [b-(b-1)\Gamma \varphi ]$. In the pressure resistance the cold layer has length $Y$, the hot injected layer has length $X-Y$, and the displaced layer has length $L-X$.

A slender planar plume in a saturated porous medium has buoyancy velocity $w=\beta (T-T_0)$, with $\beta ={\rho }_{0}g\alpha \Pi /\mu$. The leading heat equation retains both advection components and lateral thermal conduction: $u{T}_x+w{T}_z=\kappa T_{xx}$. Its convective heat input per span is conserved, $F_h=C_h\int w(T-T_0)dx$. Defining $Q=\beta F_h/C_h$ gives width proportional to $z^{2/3}$ and vertical speed/temperature excess proportional to $z^{-1/3}$.

For a slender porous thermal plume, Darcy law gives $w=\beta (T-T_0)$, with $\beta =Kg\alpha /\nu$. If $C_h={\rho }_0{c}_p$ and $F_h=C_h\int w(T-T_0)dx$ is the conserved excess convective heat flux per span, then the kinematic momentum flux is $M=\int w^{2}dx=\beta F_h/C_h$, and the buoyancy flux is $B=\int g\alpha (T-T_0)wdx=g\alpha F_h/C_h$. The hyperbolic-secant porous plume profile has $Q=(36\kappa Mz{)}^{1/3}$. Multiplication of $Q$ and $M$ by ${\rho }_0$ gives physical mass flux and the density-weighted Darcy second moment. Intrinsic pore-velocity momentum flux instead has an additional inverse-porosity factor. Momentum conservation here follows from the heat/velocity proportionality, not from an inertia balance: porous drag balances buoyancy locally.

The symmetric steady porous thermal plume has $w=W{\mathrm{sech}}^2(x/b)$, $T-T_0=w/\beta$ and $\psi =-Wb\tanh (x/b)$. Its exact normalizations are $W{b}^2=6\kappa z$ and $(4/3)W^{2}b=Q$, so $b=(48{\kappa }^2{z}^2/Q{)}^{1/3}$ and $W=(3Q^2/(32\kappa z){)}^{1/3}$. The streamfunction supplies lateral entrainment required by continuity.

Truncating the steady porous thermal plume profile below a head height and conserving supplied heat gives ${\dot{Z}}_h\simeq \int w^{2}dx/\int wdx=(2/3)W(Z_h)$. Thus $Z_h\propto t^{3/4}$. This is an integral bulk-head estimate, while fastest centreline parcels have characteristic speed $W$; the approximation does not resolve the transient nose.

A Darcy velocity $\mathbf{u}$ driven by pressure in a stationary porous skeleton dissipates mechanical power $\mu |\mathbf{u}{|}^2/\Pi$ per bulk volume. It follows by dotting Darcy law with the discharge: $-\mathbf{u}\cdot \nabla p=\mu |\mathbf{u}{|}^2/\Pi$. Its contribution to a heat equation can be neglected only when it is small relative to the relevant thermal-energy transfer.

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]

Pore velocity is the fluid velocity averaged over pore volume relative to the solid framework. For uniform porosity $\varphi$, full saturation and all pores participating, $\mathbf{v}=\mathbf{q}/\varphi$, with Darcy velocity $\mathbf{q}$ defined per bulk area. The identity follows from volume averaging: extending the relative fluid velocity by zero in the solid gives a bulk average equal to the pore fraction times the pore-volume average. It does not assert that every local passage has open area $\varphi A$ or identical speed.

In a uniform advective model, travel over length $L$ takes $L/v=\varphi L/q$. If some pore volume is inaccessible to transport, the corresponding estimate uses effective porosity ${\varphi }_e$ and $v=q/{\varphi }_e$. A connected but nearly stagnant passage can exchange solute slowly, so even this mean does not describe all arrival times. Flux-weighted residence time in a porous layer explains one precise averaging convention.^[ref]

Darcy-Bénard convection is thermal convection in a porous layer heated from below, with velocity governed by Darcy's law. With fixed temperatures and impermeable boundaries, a normalized two-dimensional model uses streamfunction $\psi$ and temperature $\theta$:The control parameter is a porous-medium Rayleigh number. On $0<z<1$ with temperature values $0,-1$, the conductive state is $\psi =0$, $\theta =-z$.

A tall porous column with square width $a$, insulated impermeable sidewalls and a maintained upward temperature decrease $\Gamma$ admits nearly vertical, zero-net-flux overturning modes. Let $B={\rho }_{0}g\alpha K/\mu$ be buoyancy mobility and $\kappa =k_e/({\rho }_0{c}_f)$ the thermal transport coefficient accompanying Darcy flux. Neumann transverse modes have $k_{\perp }^2={\pi }^2(m^2+n^2)/a^2$. In the long-column limit the lowest nonconstant modes have growth rate proportional to $B\Gamma -\kappa {\pi }^2/a^2$. They are warm ascending and cold descending columns with return flow at distant ends. The two lowest modes, varying along either side direction, are degenerate; linear theory alone does not select their nonlinear superposition.

In a unit square with prescribed conductive boundary temperatures and impermeable boundaries, the stream function convention $\mathbf{u}=(-{\psi }_z,0,{\psi }_x)$ gives linear equations $\Delta \psi =R{\theta }_x$, ${\theta }_t={\psi }_x+\Delta \theta$, with $\psi =\theta =0$ on all sides. The velocity equation follows by taking the vertical-plane curl of Darcy law, and the temperature equation by linearizing the heat equation about $1-z$. Unlike a laterally unbounded Darcy layer, the conducting side boundaries produce a two-dimensional real marginal eigenspace at $R_c=8{\pi }^2$.

For $q=q_c+ϵq_1$, ${\theta }_t=ϵ\lambda \theta$, the first-order equations are $\Delta {\psi }_1-q_c^2{\theta }_{1x}=2q_c{q}_1{\theta }_{0x}$ and ${\psi }_{1x}+\Delta {\theta }_1=\lambda {\theta }_0$. Multiply them by ${\psi }_0,q_c^2{\theta }_0$ and integrate. Integration by parts and the leading equations cancel the homogeneous correction terms, giving $2q_c{q}_1⟨{\psi }_0{\theta }_{0x}⟩+q_c^2\lambda ⟨{\theta }_0^2⟩=0$. Since $⟨{\psi }_0{\theta }_{0x}⟩=-⟨|\nabla {\theta }_0{|}^2⟩$, the displayed title formula follows. The reciprocal norm ratio is incompatible with this solvability condition. Reflection parity diagonalizes the first-order splitting of the two marginal modes.

At zero growth rate put $P=\psi +iq\theta$, $q=\sqrt{R}$. The coupled Darcy equations become $\Delta P+iq{P}_x=0$. Substituting $P=e^{-iq(x-1/2)/2}F$ gives $\Delta F+(q^2/4)F=0$ with homogeneous Dirichlet boundary conditions. The square Dirichlet Laplacian eigenvalues give $q^2=4{\pi }^2(m^2+n^2)$, $m,n\ge 1$, whose minimum is $8{\pi }^2$. For $S=\sin \pi x\sin \pi z$, $a=q_c/2$, the complex multiples $F=S,iS$ give real pairs $(S\cos (a\xi ),-S\sin (a\xi )/q_c)$ and $(S\sin (a\xi ),S\cos (a\xi )/q_c)$, $\xi =x-1/2$. Their reflection parities are opposite, yielding two independent physical modes.

With permeability $K$, depth $d$, kinematic viscosity $\nu$ and the relevant effective thermal diffusivity ${\kappa }_T$, the Darcy Rayleigh number compares buoyancy-driven Darcy advection to thermal diffusion. It differs from the clear-fluid Rayleigh number by replacing $d^3$ with $Kd$. Effective heat capacity and porosity factors may be absorbed into the definition of ${\kappa }_T$ and the thermal-time scale, and must be kept consistent with the chosen Darcy law convention.

A perturbation proportional to $\sin (n\pi z)e^{st+ikx}$ has growth rateMinimizing the neutral curve $(k^2+n^2{\pi }^2{)}^2/k^2$ gives $n=1$, $k=\pi$ and $R_c=4{\pi }^2$. Thus the first convection rolls have horizontal wavelength twice the layer depth.

For $R=4{\pi }^2+{ϵ}^2$, fixed roll phase, $T={ϵ}^{2}t$ and normalization $\psi =ϵA(T)\cos (\pi x)\sin (\pi z)+\cdots$, the Landau amplitude equation isThe mean-temperature correction in weakly nonlinear Darcy convection supplies the cubic negative feedback. Projection of the third-order equations onto the critical eigenfunction via the Fredholm solvability condition for a self-adjoint operator fixes both coefficients. The stable nonzero amplitudes are $\pm 2/\pi$ within the chosen phase, giving a supercritical pitchfork bifurcation in this real-amplitude reduction.

For ${\psi }_1=\cos (\pi x)\sin (\pi z)$, the first temperature mode is ${\theta }_1=(2\pi {)}^{-1}\sin (\pi x)\sin (\pi z)$. The streamfunction advection bracket is $J({\psi }_1,{\theta }_1)=(\pi /4)\sin (2\pi z)$, independent of $x$. The second-order equations therefore have ${\psi }_2=0$ andThis modifies the vertical temperature gradient. Its interaction with the critical roll produces the cubic term in the amplitude equation.

The Saffman–Taylor instability occurs when a more mobile fluid displaces a less mobile fluid in a porous medium or Hele–Shaw cell. Viscosity contrast amplifies interface corrugations, while surface tension suppresses sufficiently short wavelengths.

For two fluids with dynamic viscosities ${\mu }_1<{\mu }_2$, uniform permeability of a porous medium $k$, porosity $\varphi$, and imposed Darcy velocity $U$, the Saffman–Taylor instability of a planar interface has growth rate $\sigma =\alpha U({\mu }_2-{\mu }_1)/[\varphi ({\mu }_2+{\mu }_1)]$ without capillarity. With constant interfacial tension $\gamma$ and the stabilizing pressure-jump convention, a perturbation of transverse wavenumber $\alpha >0$ instead hasThe fastest-growing mode satisfies ${\alpha }_m^2=({\mu }_2-{\mu }_1)U/(3k\gamma )$.

For a sharp thermal front moving at fraction $\Gamma$ of the Darcy velocity $U$, a cold fluid of viscosity $\mu$ displacing its hot counterpart of viscosity $b\mu$ has the displayed growth rate in the local, unbounded planar approximation without thermal conduction. The harmonic pressure perturbations decay on either side. Pressure continuity and continuity of perturbed Darcy flux give the velocity response $\delta U=U|a|(b-1)\eta /(b+1)$; the thermal front's kinematic law is ${\eta }_t=\Gamma \delta U$. The prefactor is the thermal-front speed, not the material pore velocity. Increasing viscosity with heating therefore makes the trailing thermal interface unstable even when the leading displacement of the original fluid is stable.

Let fluid 1 lie above fluid 2, with mass densities ${\rho }_1<{\rho }_2$ and dynamic viscosities ${\mu }_1<{\mu }_2$. Define $\Delta \rho ={\rho }_2-{\rho }_1>0$, $\Delta \mu ={\mu }_2-{\mu }_1>0$, porosity $\varphi$, permeability of a porous medium $k$, and downward interface speed $W$. With Darcy flux $\varphi W$, no capillarity, and two unbounded layers, Darcy's law givesFor a Fourier mode of interface displacement, the two perturbation pressures decay exponentially into the respective layers. The kinematic boundary condition makes their amplitudes proportional to $\varphi {\mu }_1\sigma /(k|\alpha |)$ and $-\varphi {\mu }_2\sigma /(k|\alpha |)$. Linearized pressure continuity adds the base-gradient difference $\varphi W\Delta \mu /k-g\Delta \rho$. Eliminating the amplitudes proves the formula. The critical interface speed is $kg\Delta \rho /(\varphi \Delta \mu )$; the critical Darcy velocity is $kg\Delta \rho /\Delta \mu$. Without surface tension this model has unbounded short-wave growth rate on its unstable side.

Suppose the apparent capillary pressure is $p_1-p_2=\gamma (1+\beta x)({\kappa }_0+\nabla \cdot \mathbf{n})$, with normal towards the displaced fluid and positive apparent tension at the flat front $x=X$. Linearizing the Darcy law interface problem gives the local growth rateThe pressure increase experienced by a forward protrusion is stabilizing when ${\kappa }_0\beta >0$. It suppresses every positive wavenumber when $U\le k\gamma {\kappa }_0\beta /({\mu }_2-{\mu }_1)$. For a fixed spatial gradient, $X=Ut/\varphi$ makes the rate time dependent; constant exponential growth is a frozen-position approximation.

Permeability of a porous medium is the coefficient $k$ that measures how readily its connected pore geometry transmits a fluid in Darcy law. Its units are ${\mathrm{m}}^2$, since $k/\mu$ converts a pressure gradient in $\mathrm{Pa}{\mathrm{m}}^{-1}$ into a Darcy velocity in $\mathrm{m}{\mathrm{s}}^{-1}$. Large open, well-connected passages generally give larger permeability of a porous medium than narrow or winding passages. Porosity alone cannot determine it: two samples can contain the same void volume but have different passage sizes and connectivity.

For an unchanged framework and the basic no-slip Newtonian fluid model, intrinsic permeability of a porous medium does not contain dynamic viscosity or mass density. These properties instead enter hydraulic conductivity $K=k\rho g/\mu$. At the same horizontal pressure gradient, replacing a liquid with another of twice the dynamic viscosity halves Darcy flux even though $k$ stays fixed. Structural changes, such as swelling or clogging, can change $k$; the independence claim assumes the pore geometry stays fixed.^[ref]

An isotropic sample has scalar permeability of a porous medium. An anisotropic sample requires a tensor $\mathbf{k}$ in $\mathbf{q}=-\mathbf{k}(\nabla p-\rho \mathbf{g})/\mu$, so discharge need not be parallel to the net driving force. Layered porous permeability gives a concrete example: arranging individually isotropic materials in horizontal layers produces different vertical and horizontal effective permeabilities.

A dimensionless reduction of a porous formation's intrinsic permeability of a porous medium for one fluid phase in multiphase flow. Its effective phase permeability is $K{k}_{ri}$, and its phase mobility is $k_{ri}/{\mu }_i$. It depends on saturation and pore-scale phase configuration; it is not the magnetic use of relative permeability.

Relative permeability divided by the phase's dynamic viscosity. In the phase Darcy law, $q_i=-K{\lambda }_i\nabla p_i$ without body forces. Ratios of mobilities determine fractional flow when the phases have a common pressure gradient.

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.

Consider planar horizontal layers with positive scalar intrinsic intrinsic permeabilities $k_i$, thicknesses $b_i$, and fractions $f_i=b_i/{\sum }_j{b}_j$. Assume each layer is uniform and isotropic, fully saturated by the same Newtonian fluid of constant mass density and dynamic viscosity, with perfect hydraulic contact and steady incompressible flow. Use boundaries that impose pure flow normal to the layers or a common gradient of hydraulic head along them; additional interface resistance and finite-edge effects are excluded.

For vertical flow normal to the layers, every layer carries the same normal Darcy flux by mass conservation. The hydraulic head drops add. Writing $B={\sum }_i{b}_i$ and taking positive $\Delta h_i$ as the drop along the flow direction givesThis is the thickness-weighted harmonic mean. It is hydraulic head, rather than fluid pressure alone, that accounts correctly for gravity in the vertical calculation. A thin layer of very small permeability of a porous medium can control the whole discharge because all flow must pass through it.

For horizontal flow parallel to the layers, each layer experiences the same gradient of hydraulic head. For width $W$, each area is $A_i=W{b}_i$, and the individual volume fluxes add:This is the thickness-weighted arithmetic mean. A highly permeable layer can carry most of the horizontal flow because it supplies a continuous route from inlet to outlet. Woessner and Poeter derive both cases for hydraulic conductivity; identical fluid properties allow the common factor $\rho g/\mu$ to cancel.^[ref]

For two equally thick layers, take $k_1=k_0$ and $k_2=9k_0$, where $k_0=10^{-12}{\mathrm{m}}^2$. The values 1 and 9 are therefore permeabilities in units of $k_0$, rather than dimensionless intrinsic intrinsic permeabilities. ThenAcross the layers, the less permeable half takes $90\%$ of the total hydraulic head drop and the more permeable half only $10\%$. Along the layers, the more permeable half carries $90\%$ of the total volume flux. At the same magnitude of imposed gradient of hydraulic head, horizontal Darcy flux is $5/1.8=25/9\simeq 2.78$ times the vertical Darcy flux.

The arrangement thus creates physical anisotropy from locally isotropic ingredients. In coordinates with $x,y$ along the layers and $z$ normal to them, the effective permeability of a porous medium tensor is ${\mathbf{k}}_{\mathrm{eff}}=\mathrm{diag}(5,5,1.8)k_0$. A driving gradient of hydraulic head with both horizontal and vertical components generally produces a discharge in a different direction. This distinction matters when deciding whether a geological bed conducts aquifer flow along itself or acts as a barrier to flow across it.

Regenerate the diagram by running python ../layered-porous-permeability.py from _media; it writes its PNG into the current directory.

At small areal number density $\varphi$, identical aligned cracks contribute independently the fixed-throughflow dissipation reduction of a conductive crack. Equating the mean dissipation density to $\mu U^2/k_{xx}^{\ast }$ gives the displayed first-order enhancement. Writing $k/[1-\pi \varphi a^2\alpha /(1+\alpha )]$ is an algebraic noninteraction estimate, not a claim of accuracy at higher density. The enhancement saturates as crack conductance becomes large compared with transport through the surrounding porous matrix.

In a slender layer, two parallel sublayers of fractions $x/L$ and $1-x/L$ have local effective permeability $K(x)=k_{1}x/L+k_2(1-x/L)$. Integrating their successive hydraulic resistance gives the logarithmic mean displayed above. A leading pressure gradient uniform across the thickness does not make parcel paths horizontal; a small transverse velocity preserves mass conservation as the cross-sectional conductance changes.

Porosity is the fraction of a representative bulk material volume occupied by voids:It is dimensionless and measures storage space. A sample with bulk volume $1.0\mathrm{L}$ and pore volume $0.25\mathrm{L}$ has porosity $\varphi =0.25=25\%$. If all its pores are saturated by one liquid, the liquid occupies $0.25\mathrm{L}$; partial fluid saturation $s$ instead gives liquid volume $s\varphi V_{\mathrm{bulk}}$.

Porosity and permeability of a porous medium answer different questions: how much fluid fits, and how easily it flows. Isolated bubbles or cavities can contribute to total porosity while providing no throughgoing path. Even connected passages may be narrow enough to transmit very little discharge. In the simple fully participating saturated model, Darcy velocity is $\varphi$ times pore velocity; changing porosity at a stipulated fixed Darcy flux changes ideal advective travel time, but porosity alone cannot predict the Darcy flux.

Effective porosity is the pore-volume fraction accessible to connected transport under the specified measurement or model. It excludes isolated pore volume and may also exclude poorly exchanging regions on the transport timescale, so $0\le {\varphi }_e\le \varphi$ for total porosity $\varphi$. In a simple mobile-volume model the mean pore velocity is $q/{\varphi }_e$; the assumptions behind that model matter when interpreting tracer arrival times. A sample with total porosity $0.25$ but mobile pore fraction $0.20$ and Darcy flux $2.0\times 10^{-5}\mathrm{m}{\mathrm{s}}^{-1}$ has mobile mean pore velocity $1.0\times 10^{-4}\mathrm{m}{\mathrm{s}}^{-1}$. Woessner and Poeter discuss the role of connectivity and slowly exchanging pores.^[ref]

Suppose a porous layer has thickness $h(p)$ and permeability of a porous medium $k(p)$, and the pressure drop $p$ acts across its thickness. Darcy's law gives the drainage volume fluxIf $k/h=K{p}^{-\beta }$, then $w=(K/\mu )p^{1-\beta }$.