Capillarity-Screened Darcy Flow

• Capillarity-screened Darcy flow is a class of models where capillarity modifies classical Darcy behavior through curvature, screening, and dynamic mechanisms.
• The formulations include interface regularization, modified Helmholtz reaction terms, and pore-scale pressure thresholds to filter or limit the Darcy response.
• These approaches enhance the understanding of free-boundary dynamics, nanoporous drying, and phase drift in multiphase systems, offering practical insights for fluid transport analysis.

to=arxiv_search.search 天天中彩票无法 手机天天中彩票{"query":"all:capillarity-screened Darcy flow Hele-Shaw surface tension modified Helmholtz porous media", "max_results": 10} to=arxiv_search.search สล็อตโैनिशत्कार្យություն to=arxiv_search.search 】【。】【”】【json {"query":"id:(Matioc et al., 21 Aug 2025) OR id:(Scrobogna, 2020) OR id:(Tantardini et al., 19 Sep 2025) OR id:(Tantardini et al., 9 Nov 2025)", "max_results": 10} Capillarity-screened Darcy flow denotes a family of Darcy-type transport and free-boundary models in which capillarity changes the effective evolution by imposing curvature-driven stresses, dynamic capillary relaxations, explicit screened balance laws, or capillary thresholds that limit transmission of pressure into the flow. The phrase is not used uniformly in the literature: several relevant papers analyze the mechanism without adopting the term itself. Across these works, capillarity may enter through the Young–Laplace boundary condition in Hele–Shaw flow, through a third-order nonlocal operator in asymptotic Muskat models, through a modified-Helmholtz reaction term $K\lambda^{-2}p$ in random media, or through pore-scale Laplace ceilings that clip the Darcy pressure drop during drying (Matioc et al., 21 Aug 2025, Scrobogna, 2020, Tantardini et al., 19 Sep 2025, Vincent et al., 2015).

1. Meanings of screening in Darcy settings

In the sources considered here, “screening” does not denote a single constitutive law. At least four mathematically distinct mechanisms appear. The first is interface-level regularization: the bulk flow remains classical Darcy flow, but capillarity enters through curvature in the free-boundary condition, and the reduced interface dynamics acquires a higher-order dissipative symbol. The second is bulk screened transport in the strict modified-Helmholtz sense, where the balance law itself becomes

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$

The third is capillary thresholding or clipping, where capillary entry or Laplace pressure limits the pressure drop that can actually be transmitted into Darcy flow. The fourth is dynamic-capillarity regularization, where a capillary law of the form $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$ generates pseudo-parabolic transport (Tantardini et al., 19 Sep 2025, Vincent et al., 2015, Tantardini et al., 9 Nov 2025).

This diversity is already visible in the governing equations. In the free-boundary graph formulation of Darcy or Muskat flow, the interface velocity is written as

$h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$

so capillarity acts through the curvature $\kappa$ and the Dirichlet-to-Neumann map $G[h]$. In the pure-capillarity regime $g=0$, the linearized operator is of order three, $A^3=|D|^3$, rather than first order. By contrast, in the random-granular screened model the constitutive flux remains Darcy-like, $\mathbf q=-K\nabla p$, but the governing balance adds a zero-order screening term. In nanoporous drying, the bulk Darcy law

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$

persists only until the Kelvin-imposed tension reaches the pore Laplace threshold, after which the flux saturates at

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$0

These are different notions of “capillarity-screened Darcy flow,” even though each places capillarity in the role of suppressing or filtering the unscreened Darcy response (Scrobogna, 2020, Tantardini et al., 19 Sep 2025, Vincent et al., 2015).

A central conceptual distinction therefore separates screening of the bulk operator from screening of the effective dynamics. In many free-boundary Hele–Shaw and Muskat formulations, capillarity does not replace the harmonic pressure problem by a Helmholtz or Yukawa law. Instead, it modifies the interface law through curvature and thereby suppresses short wavelengths. This suggests that, in much of the current literature, “screening” is often best understood as a higher-order damping or regularizing effect rather than as finite-range screening in the strict elliptic sense (Matioc et al., 21 Aug 2025, Nguyen, 2023).

2. Free-boundary Darcy flow: curvature, nonlocality, and cubic damping

A mathematically precise instance of capillarity-modified Darcy flow is the two-dimensional quasistationary Hele–Shaw problem with surface tension. In a bounded planar domain $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$1 with free boundary $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$2, the bulk equations are

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$3

while the boundary conditions are

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$4

All physical parameters are normalized to $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$5. In this formulation, capillarity is entirely encoded by the Young–Laplace-type boundary condition $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$6; the bulk Darcy law itself is unchanged. For star-shaped bounded domains parameterized by a radius function $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$7, the problem is reduced to a quasilinear nonlocal evolution

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$8

with $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$9 built from double-layer and singular integral operators that eliminate the bulk harmonic pressure exactly (Matioc et al., 21 Aug 2025).

The capillary regularization mechanism is explicit in the curvature decomposition. For $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$0,

$$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$1

so curvature contributes two interface derivatives. After solving the bulk harmonic problem and taking the normal derivative, the resulting interface law has principal order three. At the unit circle, the leading operator is

$$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$2

whose Fourier symbol is $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$3. This is the sharpest spectral expression of interface-level screening in the Hele–Shaw setting: high Fourier modes are strongly damped, and the free-boundary problem becomes quasilinear parabolic in Bessel potential spaces $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$4 (Matioc et al., 21 Aug 2025).

The same cubic mechanism appears in asymptotic porous-medium models. In the quadratic approximation of one-phase Muskat or Hele–Shaw flow with gravity and capillarity,

$$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$5

and in the pure-capillarity case $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$6,

$$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$7

Here $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$8 is the Calderón operator with Fourier multiplier $$p_c^{\mathrm{dyn}}=p_c^{\mathrm{AM}}-\tau\,\partial_t S$$9, so $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$0. The paper proves local well-posedness for arbitrary $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$1 data and global well-posedness for sufficiently small $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$2 data, identifying $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$3 as the critical Sobolev space for the problem. The linear term is therefore a third-order nonlocal dissipative operator generated entirely by capillarity (Scrobogna, 2020).

Traveling-wave theory for Darcy flow with surface tension shows the same order separation. For a graph $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$4 in finite or infinite depth, the free-boundary equation is

$h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$5

and under a traveling ansatz $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$6,

$h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$7

Since $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$8, the deep-water principal operator becomes $h_t=G[h]\big((\rho_+-\rho_-)gh-\gamma \kappa\big),$9: gravity contributes a lower-order restoring term, while capillarity contributes the higher-order regularizer. The paper constructs a unique local curve of small periodic traveling waves and a connected continuum $\kappa$0 containing either arbitrarily large $\kappa$1 waves or, in finite depth, waves approaching the bottom (Nguyen, 2023).

An asymptotically flat-domain reduction yields a related screening effect in the Brinkman regime. The Brinkman Vertical Equilibrium model takes the form

$\kappa$2

with $\kappa$3 and $\kappa$4 determined nonlocally by vertically integrated mobility. The pseudo-parabolic mixed time-space terms act as a higher-order relaxation of the reduced Darcy transport, and the numerical model exhibits overshoot behavior similar to that observed for dynamical capillarity laws (Armiti-Juber et al., 2017).

3. Dynamic capillarity, global pressure, and higher-order Darcy regularizations

In multicomponent, multiphase flow, capillarity can be screened out of the total Darcy flux while remaining active in the phase-relative drifts. The isothermal Global Buckley–Leverett framework defines phase momentum by

$\kappa$5

with

$\kappa$6

The total Darcy flux is then driven by a single global pressure,

$\kappa$7

while each phase velocity is recovered by an exact fractional-flow split consisting of an advective share of total flux, a buoyancy drift, and a capillary drift. Static capillarity contributes a second-order smoothing, and dynamic capillarity contributes a pseudo-parabolic term involving $\kappa$8. Together with Maxwell–Stefan diffusion, this regularization resolves the loss of strict hyperbolicity that plagues classical three-phase Buckley–Leverett formulations (Tantardini et al., 9 Nov 2025).

A more general capillarity-regularized Darcy/Fick transport law appears in large-strain poromechanics. On a reference domain $\kappa$9, with deformation $G[h]$0, deformation gradient $G[h]$1, and concentration $G[h]$2, the actual-configuration capillarity energy is

$G[h]$3

The corresponding chemical potential is

$G[h]$4

and the species balance is

$G[h]$5

In the fixed-domain limit $G[h]$6, this reduces to a Cahn–Hilliard form with fourth-order regularization. The same actual-configuration capillarity generates a Korteweg-like stress in the momentum balance. This is not a screened Darcy law in pressure form, but it is a capillarity-regularized Darcy/Fick system with an intrinsic length scale and higher-order smoothing (Roubíček, 2019).

Diffuse-interface porous-media reduction supplies a third route. Starting from Navier–Stokes–Cahn–Hilliard dynamics for a binary fluid in the void space of a rigid porous solid, the momentum equation is intrinsically phase-averaged in the creeping-flow regime. For isotropic permeability,

$G[h]$7

where $G[h]$8 contains gravity, averaged capillary stress divergence, and wall-wetting traction contributions. A harmonic potential field is then reconstructed to satisfy impenetrability while preserving the averaged transport velocity. This suggests a coarse-grained form of capillary screening: the local momentum response is not resolved in full, but the net capillary and wetting forces are transmitted through a Darcy-like mean closure (Roudbari et al., 2016).

4. Pore-scale thresholds, nanoporous transport, and capillary limitation of Darcy flux

In nanoporous drying, capillarity screens Darcy flow not by adding derivatives to an interface equation but by imposing a pressure ceiling. For a porous silicon layer of thickness $G[h]$9, porosity $g=0$0, and nitrogen-porosimetry pore radius $g=0$1 nm, the drying flux initially obeys

$g=0$2

with Kelvin equilibrium relating external humidity to liquid pressure. The pore capillary pressure is

$g=0$3

and once $g=0$4, the flux saturates at

$g=0$5

Experimentally, the linear regime remains robust down to $g=0$6 MPa, the overall drying stresses reach approximately $g=0$7 MPa of tension, the capillary pressure inferred from the crossover is $g=0$8 MPa, and the permeability from the linear drying slope is $g=0$9. The discrepancy between the hydraulic effective radius $A^3=|D|^3$0 nm and the capillary pore radius $A^3=|D|^3$1 nm is resolved by a negative slip length $A^3=|D|^3$2 nm, interpreted as an immobile molecular layer (Vincent et al., 2015).

Capillarity-driven imbibition in mesoporous Vycor exhibits a related Darcy-scale description. For V5 and V10 matrices with mean pore radii $A^3=|D|^3$3 and $A^3=|D|^3$4, spontaneous imbibition of silicone oils obeys Lucas–Washburn square-root kinetics from seconds to $A^3=|D|^3$5 days. The height and mass laws are

$A^3=|D|^3$6

with the Lucas–Washburn–Darcy imbibition ability

$A^3=|D|^3$7

The experimentally inferred slip lengths are consistently negative, with mean $A^3=|D|^3$8, corresponding to an immobile boundary layer that reduces the hydraulically active pore radius while leaving the continuum Darcy description intact (Gruener et al., 2018).

Capillarity can also localize Darcy-driven deformation rather than merely limit flux. In the Multiphase Darcy-Brinkman-Biot framework for drainage in soft porous media, the fluid and solid momentum balances contain the capillary-force density $A^3=|D|^3$9, and the capillary entry pressure is

$\mathbf q=-K\nabla p$0

The capillary fracturing number

$\mathbf q=-K\nabla p$1

measures capillary stress against structural resistance. When $\mathbf q=-K\nabla p$2 and $\mathbf q=-K\nabla p$3, the paper identifies a capillary fracturing transition in which fractures are preceded by an invasion front, the saturation front is more uniform, and deformation is driven by more evenly distributed capillarity-induced stresses localized at the advancing invasion front. This is another screening-like mechanism: capillarity prevents broad volumetric transmission of the viscous pressure drop and concentrates forcing near the moving interface (Carrillo et al., 2020).

5. Homogenization, random media, and explicit screened-Darcy operators

Among the formulations summarized here, the most literal screened-Darcy operator appears in stationary random, polydisperse granular media. On a periodic supercell $\mathbf q=-K\nabla p$4, with phase indicator $\mathbf q=-K\nabla p$5 and transport coefficient

$\mathbf q=-K\nabla p$6

capillarity is modeled by the screened equation

$\mathbf q=-K\nabla p$7

or, in interface-height form,

$\mathbf q=-K\nabla p$8

Here

$\mathbf q=-K\nabla p$9

is the capillary length. Homogenization yields a macroscopic screened law

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$0

with

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$1

In Fourier space, the screening weight

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$2

acts as a spatial low-pass filter, and the Green function behaves like $J=-\kappa(P_{\mathrm{edge}}-P_0)/L$3 (Tantardini et al., 19 Sep 2025).

The same paper links screening to statistical representativity. The microstructure is characterized by the covariance $J=-\kappa(P_{\mathrm{edge}}-P_0)/L$4, spectral density $J=-\kappa(P_{\mathrm{edge}}-P_0)/L$5, and integral range

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$6

Screening modifies these descriptors through the capillarity-weighted volume fraction

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$7

and the screened integral range

$J=-\kappa(P_{\mathrm{edge}}-P_0)/L$8

These recover the unscreened descriptors as $J=-\kappa(P_{\mathrm{edge}}-P_0)/L$9. The resulting supercell criteria combine a shortest-edge rule,

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$00

with a volume criterion such as

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$01

This is a fully distribution-aware, quantitative formulation of capillarity-screened Darcy flow in the modified-Helmholtz sense (Tantardini et al., 19 Sep 2025).

A different form of screening emerges in Darcy-scale upscaling of multiphase flow. A recent review argues that Darcy behavior in capillary-dominated steady states is not obtained by spatial averaging alone but by space-time averaging. The paper’s key statement is that the collective energy contribution of capillary fluctuations can vanish in space-time average,

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$02

so that a linear Darcy-type constitutive law survives in steady-state connected-pathway flow. Outside that regime, capillary barriers, ganglion dynamics, and topology changes generate nonlinear rheology, intermittency, disconnected flow, and memory effects. This suggests a statistical-thermodynamic notion of screening: pore-scale capillary fluctuations remain active, but their net energetic contribution can disappear from the leading-order Darcy law after suitable averaging (Berg et al., 22 Oct 2025).

6. Stability, regime structure, and recurring misconceptions

A recurrent misconception is that capillarity-screened Darcy flow must refer to a bulk constitutive replacement of Darcy’s law. In the present literature, that is true only for the modified-Helmholtz class. In the Hele–Shaw and Muskat problems, the bulk equations remain harmonic or Darcy-like, and capillarity acts through the interface condition. The two-dimensional Hele–Shaw analysis makes this explicit: capillarity does not add a term such as $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$03 in the bulk, but it still turns the free-boundary problem into a quasilinear parabolic evolution with principal symbol $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$04, instantaneous smoothing, and exponential stability of circles modulo area and translation invariances. The linearization at the unit circle has spectrum

$$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$05

with $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$06 for $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$07 and neutral modes $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$08 arising from geometry (Matioc et al., 21 Aug 2025).

A second misconception is that screening necessarily removes nonlocality. Several papers show the opposite. In traveling Darcy waves the Dirichlet-to-Neumann operator $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$09 remains nonlocal, with symbol $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$10 in infinite depth or $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$11 in finite depth, while capillarity changes the interface operator through the higher-order factor $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$12. In steady-state multiphase upscaling, capillary topology and interfacial geometry remain essential even when the net fluctuation energetics average out. In the fractured-media Global Buckley–Leverett formulation, capillarity is screened from the total flux by the global-pressure decomposition, but it reappears as a capillary drift in each phase flux and as a pseudo-parabolic regularizer in transport (Nguyen, 2023, Berg et al., 22 Oct 2025, Tantardini et al., 9 Nov 2025).

A third misconception is that capillarity always stabilizes in the same way. The sources show several distinct outcomes. In free-boundary Hele–Shaw and asymptotic Muskat models, capillarity is a restoring mechanism that damps high wavenumbers and yields smoothing or exponential decay. In nanoporous drying, it imposes a finite pressure ceiling and produces a humidity-independent flux plateau. In drainage of soft porous media, it can localize deformation and trigger capillary fracturing when $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$13. In Brinkman-type asymptotic models and dynamic-capillarity Buckley–Leverett systems, it is associated with pseudo-parabolic regularization and overshoot rather than simple monotone decay (Vincent et al., 2015, Carrillo et al., 2020, Armiti-Juber et al., 2017).

The most precise general statement supported by these works is therefore conditional. Capillarity-screened Darcy flow is not a single equation but a class of Darcy-type descriptions in which capillarity suppresses, filters, clips, or regularizes the unscreened response. Depending on the model, the operative mechanism may be a curvature-induced $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$14 interface symbol, a dynamic-capillary $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$15 term, a modified-Helmholtz reaction $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$16, a Laplace pressure ceiling $$-\nabla\cdot\big(K(\chi)\nabla p\big)+K(\chi)\lambda^{-2}p=\rho_0.$$17, or a space-time averaging principle under which capillary fluctuations disappear from the leading-order macroscopic balance while still controlling the admissible flow regime (Matioc et al., 21 Aug 2025, Tantardini et al., 19 Sep 2025, Vincent et al., 2015, Berg et al., 22 Oct 2025).