Codex @codex  0

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.

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]

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]

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.

For a stationary solid skeleton, the Darcy velocity $\mathbf{u}$ is liquid volume transport per unit total cross-sectional area, while the pore velocity $\mathbf{v}$ is the intrinsic mean fluid velocity within the connected pore space. ThusFor an incompressible liquid in a saturated region, conservation of mass is ${\varphi }_t+\nabla \cdot \mathbf{u}=0$. With time-independent but spatially varying porosity, it becomesnot generally $\nabla \cdot \mathbf{v}=0$. Neglecting dispersion, a fluid interface is advected by the pore velocity; its mean normal speed is $V_n=\mathbf{v}\cdot \mathbf{n}=\mathbf{u}\cdot \mathbf{n}/\varphi$.

For the thin layer on the inclined plane, the gravitational acceleration components are $g\sin \theta$ downslope and $-g\cos \theta$ normal to the plane. The hydrostatic approximation and constant pressure $p_a$ above the current giveWith permeability of a porous medium $k(z)$, Darcy's law therefore givesHere $u_x,u_y$ denote the components of the Darcy velocity. Integrate over depth, using the prescribed power laws:Both integrals are finite for the positive exponents in the model. Conservation of the stored fluid gives ${\mathcal{V}}_t+{\nabla }_{x,y}\cdot \mathbf{Q}=0$, orFor $0<\theta <\pi /2$, set $X=x\tan \theta$, $Y=y\tan \theta$ andSince $\cos \theta {\tan }^2\theta =\sin \theta \tan \theta$, the depth-dependent inclined porous-current equation isThe longitudinal flux carries fluid downslope, while the thickness-gradient terms spread it. The degenerate horizontal or vertical slope limits need their own coordinate scaling.

Choose $x$ upslope and $z$ normally downwards from the impermeable roof, so the buoyant layer occupies $0<z<h(x,t)$. Use the hydrostatic approximation and neglect the ambient return-flow pressure correction of relative order $h/H$. Pressure continuity at the interface gives the pressure inside the buoyant layer asThe along-slope Darcy velocity is thereforeDefine $u=kg\Delta \rho \sin \theta /\mu$. Integrating over depth gives the volume flux per unit width $F=uh-u\cot \theta h{h}_x$. Local fluid-volume conservation, including porosity, is $\varphi h_t+F_x=0$. Hence the inclined porous gravity current satisfiesThe first term is upslope advection driven by the mass density difference; the second is nonlinear hydrostatic spreading. The Darcy velocity $u$ and the corresponding pore velocity $u/\varphi$ must be distinguished.

Let $Y(t)$ be the porous thermal front and $X(t)$ the polymer front. With the literal stated convention for Darcy velocity,Therefore $Y=\beta X$ with $\beta =\Gamma \varphi$. This is thermal-front retardation in Darcy flow. The region $0<x<Y$ is cold injected fluid of dynamic viscosity $\mu$; $Y<x<X$ is heated injected fluid of dynamic viscosity $b\mu$; $X<x<L$ is oil of dynamic viscosity ${\mu }_h$. The pressure drop is consequentlywhere ${\mu }_{\mathrm{eff}}=\mu [\beta +b(1-\beta )]$. Apply the preceding fixed-pressure planar Darcy displacement calculation with $\mu$ replaced by this effective dynamic viscosity. The polymer-front speed isIts position solves ${\mu }_{h}LX+({\mu }_{\mathrm{eff}}-{\mu }_h)X^2/2=k\Delta Pt/\varphi$, and $Y=\Gamma \varphi X$. These expressions hold while both interfaces remain in the field, the porous thermal front is sharp and stable, and $0\le \Gamma \varphi \le 1$. For $b=1$ the original constant-viscosity answer is recovered.

Sometimes a thermal retardation factor is defined relative to the pore velocity instead. Under that alternative definition one substitutes $\beta =\Gamma$ throughout. The factor of $\varphi$ in the boxed answer follows specifically from the PDF's phrase “fraction of the Darcy speed.”

Let $X(t)$ denote the polymer–oil material interface, measured from the injection well, and let $U(t)$ be the Darcy velocity. Incompressibility and uniform cross-sectional area make $U$ independent of $x$ in this planar model. Integrating Darcy's law through the injected layer and the oil layer givesThe material interface moves at the pore velocity, so $X^{\prime }=U/\varphi$. Put $a=\mu -{\mu }_h$ and $A={\mu }_{h}L$. Then $(A+aX)X^{\prime }=k\Delta P/\varphi$, which integrates from $X(0)=0$ toThus the fixed-pressure planar Darcy displacement hasFor $a\ne 0$, $X=[\sqrt{A^2+2ak\Delta Pt/\varphi }-A]/a$; its continuous $a=0$ limit is $X=k\Delta Pt/(\varphi {\mu }_{h}L)$. Positive layer dynamic viscosities ensure the resistance stays positive until breakthrough, which occurs at $t_b=\varphi L^2(\mu +{\mu }_h)/(2k\Delta P)$. The formulas stop there, since the two-layer geometry then changes.

The imposed pressure gradient and Darcy's law give a linear velocity profile. Use $U$ as the mean pore velocity in the contaminant advection-diffusion equation; thenThe last conditions express zero transverse solute flux under the standard closed-layer interpretation. Transverse exchange through the layer boundaries would change this cell problem and requires additional boundary data. If the prescribed mean speed is a Darcy velocity, its corresponding pore velocity is that speed divided by porosity.

The given parameter is the ratio of transverse mixing time to travel time:Thus transverse dispersion mixes the contaminant across the layer many times during its passage along the long layer. Different fluid speeds still matter: repeated exchange between slow and fast portions produces Taylor dispersion, rather than transport solely at each parcel's original speed. The averaged description applies after $t\gg H^2/D_T$ and on longitudinal scales large compared with the distance travelled during that mixing time.

Write $C=\overline{c}$ and, to leading correction order, $c=C+\chi (y)C_x+\cdots$, with $\overline{\chi }=0$. Taking $C_t\simeq -U{C}_x$ in the local equation gives the cell problemWith $\eta =y/H$, its solution isAveraging the local advection-diffusion equation now yields $C_t+U{C}_x=[D_L-\overline{(u-U)\chi }]C_{xx}$. An integration by parts in the cell problem gives $-\overline{(u-U)\chi }=D_T\overline{({\chi }^{\prime }{)}^2}$. Since ${\int }_0^1({\eta }^2-\eta {)}^{2}d\eta =1/30$, the effective equation isThis is Taylor dispersion in a linear porous-layer velocity profile. For a localized pulse away from the inlet and outlet, its mean position advances at $U$ and its longitudinal variance grows as $2D_{\mathrm{eff}}t$.

For the inlet step, use the semi-infinite inlet approximation $x>0$, $C(0,t)=1$, $C(x,0)=0$ and $C\to 0$ as $x\to \infty$ at fixed $t$. The resulting constant-concentration inlet solution isOne derivation is to Laplace transform in time: $\tilde{C}(x,p)=p^{-1}\exp [(U-\sqrt{U^2+4D_{\mathrm{eff}}p})x/(2D_{\mathrm{eff}})]$. Inverting gives the two complementary error functions. Their sum is one at $x=0$; for fixed $x>0$ both terms vanish as $t\downarrow 0$, and substitution verifies the averaged equation. Retaining only the first term gives the familiar moving error-function front far from the inlet when longitudinal advection dominates diffusion, but does not satisfy the inlet condition exactly.

For a literal finite layer $0<x<L$, a downstream boundary condition is also needed once the outlet influences the solution. The PDF's step-inlet request specifies only the initially clean region $x>0$, so the formula above is the semi-infinite/long-layer solution, not a claim of a unique finite-interval solution with unspecified outlet data.

Let $K$ be intrinsic permeability of a porous medium, $\varphi$ the mobile pore fraction, and let $h$ measure thickness normal to the sloping cap. Define cold and reservoir-temperature properties byThe imposed inequality makes ${\delta }_c>0$, so the cold fluid is still buoyant. Assume $K>0$, $\sin \theta >0$, positive densities and viscosities, and the usual positive thermal-expansion and viscosity coefficients. Away from localized fronts, neglect the along-current thickness gradient. Hydrostatic water pressure and Darcy's law then give the upslope Darcy velocitiesTypically $u_h>u_c$: warming increases buoyancy and decreases dynamic viscosity.

Current depths and the moving moving thermal front in a porous current. For this slope-driven porous gravity current, the injected mass flux per unit well length fixes the cold plateau thickness:The moving thermal front in a porous current moves at $V_T=\Gamma u_c$. Since $u_c$ is explicitly a Darcy velocity, the fluid mass stored per unit horizontal area is $\varphi \rho h$. Apply mass conservation across the moving moving thermal front in a porous current:Writing $\gamma =\varphi \Gamma$ and $r=u_c/u_h$, this gives the distal warm depthFor a retarded moving thermal front in a porous current take $0\le \gamma <1$; otherwise this two-plateau geometry needs reconsideration. The factor $(1-\gamma )/(1-\gamma r)$ is important: the distal mass flux is not generally $Q$, because mass is being stored as the cold region replaces the warm region. Setting the two fluxes equal would silently assume a stationary moving thermal front in a porous current.

If the model absorbs porosity into storage and uses pore velocity as its transport speed, the same formulas use $\gamma =\Gamma$. The explicit Darcy convention in the PDF instead gives $\gamma =\varphi \Gamma$; porosity must be specified or absorbed consistently. In the stationary-front limit $\gamma =0$, the two steady depths reduce to $Q/({\rho }_c{u}_c)$ and $Q/({\rho }_h{u}_h)$. In the Boussinesq approximation, use a common reference density in the mass factors while retaining ${\delta }_c,{\delta }_h$ in the driving force.

Leakage thresholds. A normal thickness $h_i$ generates cap overpressure ${\Pi }_i={\delta }_{i}g{h}_i\cos \theta$. Before leakage modifies the current, the two plateau values areFor the caprock leakage threshold $\Delta p>0$, the candidate critical mass-injection rates are thereforeCold-only leakage occurs for $Q_c^{\ast }\le Q<Q_h^{\ast }$, provided $Q_c^{\ast }<Q_h^{\ast }$. Both plateaus can leak once $Q\ge \max (Q_c^{\ast },Q_h^{\ast })$. In the usual common-density approximation, ${\mu }_c>{\mu }_h$ and $r<1$ ensure $Q_c^{\ast }<Q_h^{\ast }$, giving the expected sequence: no leakage, cold-only leakage, then leakage from both cold and warm regions. Without that approximation the ordering must be checked; the given inequality $\Delta \rho >\alpha \Delta T$ alone does not establish it. These are onset criteria computed on the nonleaking current, not a post-leakage mass budget.

For small slopes $\cos \theta \simeq 1$ and $\tan \theta \simeq \sin \theta \simeq \theta$. If another thickness convention is used, its hydrostatic column and projected flux must be changed consistently; one should not mix a normal thickness with a vertical-pressure formula lacking the cosine.

Cross-current heat conduction. Heat transfer from the warm formation makes temperature vary across the carbon dioxide depth and introduces a warming time controlled by thickness, thermal diffusivity and the surrounding rock's heat capacity. It smooths the sharp thermal adjustment and causes gradual changes in buoyancy, viscosity and the velocity profile. The cold fluid generally warms, becomes more mobile and requires less depth to carry a prescribed mass flux; the enhanced cold-region overpressure and its distinct leakage zone tend to shrink. Warming may occur before a parcel reaches the idealized advective moving thermal front in a porous current, especially for a thin current. A single supplied $\Gamma$ no longer describes all heat transport. Quantitative depths and thresholds then require a coupled temperature equation and thermal boundary data; their exact changes cannot be inferred from the linear property laws alone.

Take $z$ vertically upward and the porous skeleton stationary. Darcy's law and incompressibility areHere $\mathbf{u}$ is the Darcy velocity, a volume flux per unit total cross-sectional area, including solid and void. If $\mathbf{v}$ is the liquid pore velocity and all the pore space is saturated and connected, thenThe porosity $\varphi$ therefore enters storage and the conversion to pore velocity; it does not multiply the Darcy velocity again when calculating discharge.

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]

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.

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.

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.

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.