---
title: Very Thin Porous Medium (VTPM)
url: https://www.emergentmind.com/topics/very-thin-porous-medium-vtpm
type: topic
---

# Very Thin Porous Medium (VTPM)

Searching arXiv for recent and foundational papers on very thin porous media, thin porous layers, and effective interface models.
Very Thin Porous Medium (VTPM) denotes, in the cited literature, a porous layer whose thickness is small relative to its in-plane dimensions and whose effective behavior is determined by a combined thin-domain and homogenization limit. Depending on the geometric scaling, a VTPM may appear as a three-dimensional perforated layer of thickness \(\varepsilon\), as a thin porous medium with thickness \(h_\varepsilon\) and pore scale \(\varepsilon \ll h_\varepsilon\), or, in a singular limit, as a zero-thickness interface or screen carrying transmission conditions instead of a volumetric porous law [2508.03214, 2508.04629, 2105.02063, 2106.05994].

## 1. Geometric archetypes and scale regimes

A standard geometric realization of a VTPM is a domain of the form
\[
\Omega_\varepsilon=\omega_\varepsilon\times(0,\varepsilon),
\]
where \(\omega_\varepsilon\subset\mathbb{R}^2\) is obtained by periodically removing solid cylinders from a bounded planar set \(\omega\). In the Carreau formulation of VTPM, the cylinders are periodic at scale \(\varepsilon^\ell\) with \(0<\ell<1\), so that \(\varepsilon\ll \varepsilon^\ell\); this regime is explicitly identified as “very thin porous medium (VTPM)” and contrasted with homogeneously thin porous media (HTPM), where \(\varepsilon^\ell\ll \varepsilon\), and proportionally thin porous media (PTPM), where \(\varepsilon^\ell=\varepsilon\) [2508.03214]. Closely related micropolar and Bingham analyses use the ratio \(a_\varepsilon/\varepsilon\), with the VTPM regime corresponding to \(a_\varepsilon\gg \varepsilon\) [2004.11756, 1903.02769].

A second archetype is a thin porous medium whose thickness \(h_\varepsilon\) is small at the macroscopic scale but still much larger than the pore size, with
\[
\lim_{\varepsilon\to 0}\frac{\varepsilon}{h_\varepsilon}=0.
\]
In this setting there are many pore layers across the thickness, and the thin-domain limit produces a lower-dimensional Darcy law while preserving three-dimensional cell problems in the periodic pore geometry [2508.04629]. A coupled thin-film/thin-porous-medium configuration extends this geometry by adjoining a thin free-fluid film of thickness scale \(\eta_\varepsilon\) across an interface \(\Sigma\), again under the assumptions \(\varepsilon\ll h_\varepsilon\) and \(\varepsilon\ll \eta_\varepsilon\) [2512.14275].

A third archetype replaces the thin porous layer by a collapsed interface. In the tumour-invasion model, a membrane of thickness \(\varepsilon\) separating two bulk domains shrinks to the plane \(\tilde\Gamma_{1,3}\), and the asymptotic problem is posed directly on the two bulk domains plus a transmission law on the limiting interface [2105.02063]. In hemodynamics, a flow diverter is modeled not as a volumetric porous block but as a zero-thickness internal surface, or “screen,” carrying a pressure-jump law whose coefficients vary along the device [2106.05994]. Taken together, these configurations show that VTPM is not tied to a single geometry: the shared feature is the singular role of a porous structure whose thickness is negligible at the macroscopic scale.

## 2. Governing equations and constitutive descriptions

The governing equations associated with VTPM span several constitutive classes. They include porous-medium transport equations for cell populations, Newtonian Stokes flow with slip, micropolar Stokes systems, generalized Newtonian Carreau laws, viscoplastic Bingham laws, and thin-screen momentum jumps.

| Model class | Representative relation or structure | Representative papers |
|---|---|---|
| Porous-medium transport | \(\partial_t u-\nabla\cdot(\mu u\nabla p)=uG(p)\), \(p=u^\gamma\) | [2105.02063] |
| Newtonian Stokes with slip | Stokes system in a thin porous domain with non-homogeneous slip depending on \(\gamma\) | [1712.05998] |
| Micropolar Stokes | Velocity \(u_\varepsilon\), pressure \(p_\varepsilon\), microrotation \(w_\varepsilon\) | [2004.11756], [2508.04629] |
| Carreau fluid | \(\eta_r(\mathbb{D}[u])=(\eta_0-\eta_\infty)(1+\lambda|\mathbb{D}[u]|^2)^{\frac r2-1}+\eta_\infty\) | [2312.01844], [2508.03214] |
| Bingham fluid | \(\sigma=-pI+2\mu D(u)+\sqrt{2}\,g\,D(u)/|D(u)|\) | [1903.02769] |
| Thin screen law | \(\Delta p=-\alpha_v u-\alpha_i u^2\) | [2106.05994] |

In the porous-medium membrane problem, the unknown is a cell density \(u\ge 0\) with pressure \(p(u)=u^\gamma\), \(\gamma>1\), and Darcy velocity \(\mathbf v=-\mu\nabla p\). Mass conservation yields
\[
\partial_t u+\nabla\cdot(u\mathbf v)=uG(p),
\]
equivalently
\[
\partial_t u-\nabla\cdot(\mu u\nabla p)=uG(p),
\]
so the diffusion operator is of porous-medium type and degenerates where \(u=0\) [2105.02063]. In the thin-screen hemodynamic model, by contrast, the governing object is not a volumetric permeability field but an interfacial pressure drop
\[
\Delta p=-\left(\frac{\mu}{d_h}u+\frac{\rho}{2}\Big[1.3(1-\beta)+\Big(\frac1\beta-1\Big)^2\Big]u^2\right),
\]
with local coefficients determined from porosity \(\beta\) and hydraulic diameter \(d_h\) [2106.05994].

The Newtonian slip problem studies the Stokes system in a thin porous medium perforated by periodically distributed solid cylinders, with a Robin-type non-homogeneous slip law on the obstacle boundaries involving a parameter \(\gamma\). The body force is horizontal, the surface force is periodic, and the limit law depends sensitively on the slip scaling \(\alpha\varepsilon^\gamma\) [1712.05998]. The micropolar formulations augment velocity and pressure with a microrotation field \(w_\varepsilon\), coupled through curl terms and controlled by a coupling number \(N^2\) and a characteristic micropolar length \(R_M\) [2004.11756, 2508.04629].

For generalized Newtonian fluids, the Carreau viscosity produces several distinct VTPM limits. In one thin porous setting the asymptotics depend on both the flow index \(r\) and a viscosity scaling \(\varepsilon^\gamma\), yielding linear Darcy laws, nonlinear Carreau-type Darcy laws, or nonlinear power-law Darcy laws according to the regime [2312.01844]. In a VTPM with pore scale \(\varepsilon^\ell\), the same rheology leads to reduced Hele–Shaw-type cell problems and explicit filtration laws for pseudoplastic and dilatant fluids [2508.03214]. The Bingham case replaces shear-dependent viscosity by a yield-stress inequality, and the effective laws remain nonlinear in all three geometric regimes of \(a_\varepsilon/\varepsilon\) [1903.02769].

## 3. Asymptotic machinery: rescaling, unfolding, extension, and compactness

A common analytic step is vertical rescaling. For layers of thickness \(\varepsilon\), one sets \(y_3=x_3/\varepsilon\); for layers of thickness \(h_\varepsilon\), one sets \(z_3=x_3/h_\varepsilon\). This transforms the thin domain into a fixed-height set and introduces anisotropic differential operators such as
\[
\nabla_\varepsilon\psi=(\nabla_{x'}\psi,\varepsilon^{-1}\partial_{z_3}\psi)^t
\]
or
\[
\nabla_{h_\varepsilon}\widetilde\psi=(\nabla_{z'}\widetilde\psi,h_\varepsilon^{-1}\partial_{z_3}\widetilde\psi)^t,
\]
so that vertical derivatives explicitly carry the thinness scale [1712.05998, 2312.01844, 2508.04629].

Periodic unfolding is then used to separate slow and fast variables. In the Newtonian slip problem the unfolding acts only in the horizontal variables and converts bulk and boundary integrals in the thin perforated domain into integrals over \(\omega\times Y_f\) and \(\omega\times\partial T\) [1712.05998]. In the micropolar thin porous medium with \(\varepsilon\ll h_\varepsilon\), unfolding maps the rescaled porous set into \(\mathbb R^3\times Y_f\) and yields uniform control of unfolded velocity, microrotation, and pressure-corrector fields [2508.04629]. In the VTPM Carreau problem, the unfolding is adapted to the mixed scales \(\varepsilon\) and \(\varepsilon^\ell\), converting functions on \(\widetilde\Omega_\varepsilon\) into functions on \(\omega\times Z_f\) and allowing the limit equations to be written on a fixed product domain [2508.03214]. In the thin-film/VTPM coupling, an unfolding tailored to the porous layer is combined with separate rescaling in the film region [2512.14275].

When the coefficients degenerate inside the thin layer, compactness may fail unless the solutions are extended. In the membrane problem, the mobility in the membrane satisfies \(\mu_{2,\varepsilon}\to0\), so the natural \(L^2\)-bound controls only \(\mu_{2,\varepsilon}\int|\nabla p_{2,\varepsilon}|^2\), not \(\|\nabla p_{2,\varepsilon}\|_{L^2}\). The analysis therefore constructs an explicit extension operator \(P_\varepsilon\) by reflecting functions from the bulk domains into the thin layer; this yields strong convergence of \(P_\varepsilon u_\varepsilon\) and \(P_\varepsilon p_\varepsilon\) and weak convergence of their gradients on \(\Omega\setminus\tilde\Gamma_{1,3}\) [2105.02063]. In the Carreau VTPM, pressure decomposition \(p_\varepsilon=p_\varepsilon^0+p_\varepsilon^1\) and extension operators provide compactness for the macro-pressure while showing that the micro-pressure vanishes in the limit [2508.03214]. The power-law thin-film/VTPM coupling uses a restriction operator preserving divergence-free fields in the porous layer and a De Rham-type pressure reconstruction [2512.14275].

For nonlinear constitutive laws, limit passage relies on monotonicity. The Carreau analyses explicitly use Minty’s method to identify the nonlinear limit operator in the critical scaling \(\gamma=1\) and in the dilatant power-law regime \(\gamma>1\) [2312.01844, 2508.03214]. The Bingham analysis passes to the limit in convex plastic dissipation functionals and obtains variational inequalities of Bingham type in the homogenized problems [1903.02769]. A plausible implication is that VTPM analysis is less a single theorem than a recurring technical pattern: anisotropic rescaling identifies the dominant gradients, unfolding resolves the pore geometry, and compactness or monotonicity arguments decide which nonlinearities survive.

## 4. Effective laws and lower-dimensional closures

The effective models produced by VTPM asymptotics are uniformly lower-dimensional, but they are not uniform in form. Linear Darcy laws, generalized Darcy laws with microrotation, nonlinear power-law and Carreau filtrations, Bingham-type constitutive maps, Kedem–Katchalsky transmission conditions, and quadratic screen jumps all arise in the cited literature.

For generalized Newtonian Carreau flow in a VTPM with \(0<\ell<1\), the pseudoplastic problem with \(\gamma\neq1\) and the dilatant problem with \(\gamma<1\) reduce to a linear Darcy law
\[
\widetilde V'(x')=-\frac{1}{6\eta}A\bigl(f'(x')-\nabla_{x'}\widetilde p(x')\bigr),\qquad \widetilde V_3\equiv0,
\]
where \(A\) is a permeability tensor defined by 2D Hele–Shaw cell problems on \(Z_f'\), and \(\eta\) is \(\eta_0\) or \(\eta_\infty\) depending on the scaling regime [2508.03214]. In the proportionally thin Carreau medium, the corresponding linear regimes yield
\[
\tilde V'(x')=\frac1\eta\,\mathcal A\bigl(f'(x')-\nabla_{x'}\tilde P(x')\bigr),
\]
with \(\mathcal A\) computed from 3D cell Stokes problems in \(Z_f\) [2312.01844]. These are Darcy laws in the strict sense: the mean in-plane flux is linearly related to the in-plane pressure gradient and forcing.

Micropolar VTPM limits retain coupled translational and rotational response. In the thin porous medium with \(\varepsilon\ll h_\varepsilon\), the averaged fields satisfy
\[
\mathcal{\widetilde U}'(z')=K^{(1)}\Big(f'(z')-\nabla_{z'}\widetilde p(z')\Big)+K^{(2)}g'(z'),
\qquad
\mathcal{\widetilde W}'(z')=L^{(1)}\Big(f'(z')-\nabla_{z'}\widetilde p(z')\Big)+L^{(2)}g'(z'),
\]
together with
\[
{\rm div}_{z'}\left(K^{(1)}\Big(f'(z')-\nabla_{z'}\widetilde p(z')\Big)+K^{(2)}g'(z')\right)=0
\]
in \(\omega\) [2508.04629]. In the VTPM regime \(a_\varepsilon\gg\varepsilon\), the earlier micropolar model yields
\[
U'(x')=K^{(1)}\big(f'(x')-\nabla_{x'}P(x')\big)+K^{(2)}g'(x'),
\qquad
W'(x')=L^{(2)}g'(x'),
\]
with a generalized Darcy equation for \(P\) and with \(W_3=U_3=0\) [2004.11756].

Nonlinear rheologies survive in critical scalings. In the VTPM Carreau model, the case \(\gamma=1\) produces a nonlinear Darcy law of Carreau type, with \(\widetilde V'\) given by an explicit nonlinear function of \(\nabla_{x'}\widetilde p-f'\) through a 2D nonlinear cell problem [2508.03214]. In the proportionally thin Carreau model, the same critical scaling yields
\[
\tilde V'(x')=\mathcal U\bigl(f'(x')-\nabla_{x'}\tilde P(x')\bigr),
\]
where \(\mathcal U\) is defined through 3D Carreau cell problems in \(Z_f\) [2312.01844]. For dilatant fluids with \(\gamma>1\), the VTPM limit becomes a nonlinear power-law Darcy law governed by a monotone coercive map \(\mathcal U:\mathbb R^2\to\mathbb R^2\) [2508.03214].

For Bingham fluids the limit law remains nonlinear in all three geometric regimes \(a_\varepsilon/\varepsilon\to0\), \(a_\varepsilon/\varepsilon\to\lambda\in(0,\infty)\), and \(a_\varepsilon/\varepsilon\to+\infty\). The macroscopic filtration velocity takes the form
\[
\tilde U'(x')=K^\lambda\bigl(f'(x')-\nabla_{x'}\hat P(x')\bigr),\qquad {\rm div}_{x'}\tilde U'(x')=0,
\]
where \(K^\lambda\) is a nonlinear map determined by a Bingham cell problem that is 3D in the critical regime, 2D in the subcritical regime, and 1D in the supercritical regime [1903.02769].

The interface and screen formulations are lower-dimensional in a different sense. In the porous-medium membrane problem, the limit bulk equations are posed on the two sides of \(\tilde\Gamma_{1,3}\), while the vanished membrane is encoded through
\[
\tilde\mu_{1,3}\,\llbracket \Pi \rrbracket
=
\tilde\mu_1\,\tilde u_1\nabla\tilde p_1\cdot\tilde{\boldsymbol n}_{1,3}
=
\tilde\mu_3\,\tilde u_3\nabla\tilde p_3\cdot\tilde{\boldsymbol n}_{1,3},
\qquad
\Pi(u)=\frac{\gamma}{\gamma+1}u^{\gamma+1},
\]
which is identified as a nonlinear generalized Kedem–Katchalsky condition [2105.02063]. In the flow-diverter screen model, the lower-dimensional closure is
\[
[p]=\Delta p=-\alpha_v u_n-\alpha_i u_n^2,
\]
with \(\alpha_v=\mu/d_h\) and \(\alpha_i=\frac{\rho}{2}\left[1.3(1-\beta)+\left(\frac1\beta-1\right)^2\right]\) [2106.05994]. The literature therefore does not support a single universal “VTPM law”; it supports a family of asymptotically consistent closures whose form is dictated by geometry, rheology, and scale separation.

## 5. Role of scale ratios and critical parameters

The central structural variable in VTPM theory is the ratio between thickness and pore scale. One formulation uses \(a_\varepsilon/\varepsilon\): in micropolar and Bingham models, \(a_\varepsilon\sim\varepsilon\) yields a critical regime, \(a_\varepsilon\ll\varepsilon\) a homogeneously thin regime, and \(a_\varepsilon\gg\varepsilon\) the VTPM regime proper [2004.11756, 1903.02769]. Another formulation uses \(\varepsilon^\ell\) for the pore scale in a layer of thickness \(\varepsilon\); then \(0<\ell<1\) implies \(\varepsilon\ll\varepsilon^\ell\) and defines VTPM, while \(\ell=1\) and \(\ell>1\) correspond to PTPM and HTPM respectively [2508.03214]. A third formulation uses a thin porous thickness \(h_\varepsilon\) with \(\varepsilon\ll h_\varepsilon\), so the medium is thin at the macro-scale yet contains many pore layers across its thickness [2508.04629, 2512.14275].

These ratios determine the local problem dimension. In the micropolar thin porous medium, the VTPM regime leads to a micropolar Reynolds-type system with derivatives only in the vertical coordinate \(y_3\) and 2D cross-sectional cell problems, whereas the PTPM regime requires fully 3D local Stokes–micropolar problems and the HTPM regime yields 2D micropolar local problems [2004.11756]. In the Bingham case, the same trichotomy leads to 3D, 2D, and 1D Bingham-type cell problems respectively [1903.02769]. In the Carreau VTPM, the authors state that compared to proportionally thin media, the microproblems simplify to Hele–Shaw-type problems involving only pressure, not a full velocity–pressure Stokes system [2508.03214].

A second class of critical parameters concerns constitutive scaling. In the slip Stokes problem, the Robin coefficient scales like \(\alpha\varepsilon^\gamma\), and three limit regimes are derived for \(\gamma<-1\), \(-1\le\gamma<1\), and \(\gamma\ge1\). These regimes determine whether the limit Darcy law is driven only by pressure, by pressure plus body and surface forces, or by a tensorial permeability coming from 2D cell Stokes problems [1712.05998]. In Carreau flow, the viscosity prefactor \(\varepsilon^\gamma\) plays the analogous role: \(\gamma<1\) yields a linear Darcy law with viscosity \(\eta_0\), \(\gamma=1\) preserves the full Carreau nonlinearity, and \(\gamma>1\) produces either a linear law with viscosity \(\eta_\infty\) or a nonlinear power-law Darcy law depending on whether \(1<r<2\) or \(r>2\) [2312.01844, 2508.03214].

The coupled thin-film/VTPM problem adds a third scale and a genuine critical balance. For a power-law fluid of flow index \(r\in(1,+\infty)\), the critical regime is
\[
h_\varepsilon\approx \eta_\varepsilon^{\frac{2r-1}{r-1}}\varepsilon^{-\frac{r}{r-1}},
\]
under which the effective model becomes a coupled \(1\)D Darcy law in the porous layer and \(1\)D Reynolds law in the film [2512.14275]. This suggests that “critical” in VTPM theory has a precise meaning: it designates scale balances under which multiple substructures remain visible in the macroscopic law instead of one asymptotic mechanism dominating the others.

## 6. Applications, modeling choices, and limitations

The motivating applications are diverse. The porous-medium membrane problem is explicitly motivated by biological applications on tumour invasion through thin membranes, with cell-density transport under Darcy’s law and an effective nonlinear generalized Kedem–Katchalsky transmission condition on the collapsed membrane [2105.02063]. The Newtonian slip problem lists manufacturing of fibre-reinforced polymer composites via liquid moulding processes, passive mixing in microfluidic devices, paper-making processes, block copolymer thin-film self-assembly on nanometer scales, and heat and mass transfer equipment using fibrous or porous layers [1712.05998]. The micropolar thin porous formulations point to filters, coatings, gaskets, microfluidic devices, lubrication with structured fluids, and biological tissues [2508.04629]. The coupled thin-film/VTPM model is motivated by flows in the oil industry, such as non-Newtonian drilling muds or fracturing fluids infiltrating porous rock around a wellbore [2512.14275].

The hemodynamic application is the most quantitatively developed. The thin inhomogeneous porous medium (iPM) model represents a cerebral aneurysm flow diverter as a zero-thickness internal surface whose resistance is assigned pore by pore from local device geometry. In three patient-specific aneurysms, iPM predictions of aneurysm-averaged velocity, inflow rate, shear rate, and turnover time differed by only \(1\%\)–\(6\%\) from explicit-device CFD, summarized as \(94\%\)–\(99\%\) accuracy, while “iPM CFD ran 500% faster” than explicit CFD; the reported mesh reduction is about \(78\%\) and runtime reduction about \(80\%\) [2106.05994]. This does not imply that every VTPM model is a screen model, but it does show that interface-type representations can be both asymptotically appropriate and computationally decisive when the porous structure is extremely thin.

A recurrent modeling distinction concerns whether the porous structure should be resolved as a volume or collapsed to a surface. The cited works show both possibilities. Very thin perforated domains with no-slip, no-spin, slip, Carreau, micropolar, or Bingham constitutive laws are treated by volume homogenization and dimension reduction [1712.05998, 2004.11756, 1903.02769, 2312.01844, 2508.03214, 2508.04629]. By contrast, the tumour membrane and the flow-diverter screen are replaced in the limit by interfacial laws on zero-thickness sets [2105.02063, 2106.05994]. A common misconception is therefore that a VTPM must always be represented by a bulk Darcy–Forchheimer block. The literature instead supports a more conditional statement: when the thickness becomes asymptotically negligible, an interface law or screen law may be more faithful than a volumetric porous approximation.

The limitations are likewise systematic. The thin porous analyses assume periodic microstructure, smooth obstacles, rigid solid matrices, and steady or quasi-steady low-Reynolds-number settings in the Stokes regime [1712.05998, 2004.11756, 2312.01844, 2508.04629]. The aneurysm iPM model assumes incompressible Newtonian blood, rigid no-slip walls, and steady-state flow, and it was validated only for a single type of flow diverter, the Pipeline Embolization Device, with extension to other devices or overlapping layers left open [2106.05994]. The micropolar and Carreau studies fix specific constitutive scalings such as \(R_M=\varepsilon^2R_c\) or \(\varepsilon^\gamma\), and the resulting effective laws are rigorously valid within those regimes rather than universally [2004.11756, 2508.04629, 2312.01844]. The Bingham work similarly ties the persistence of nonlinear yield behavior to carefully chosen scalings of \(g(\varepsilon)\) with \(a_\varepsilon\) or \(\varepsilon\) [1903.02769].

Across these applications and limitations, VTPM emerges not as a single equation but as a modeling framework for singularly thin porous structures. Its defining mathematical content is the replacement of a three-dimensional microscopic problem by an effective lower-dimensional description whose coefficients, constitutive maps, or transmission laws still encode pore geometry, rheology, and dominant scale ratios.

Source: https://www.emergentmind.com/topics/very-thin-porous-medium-vtpm