---
title: Time-Fractional Porous Medium Equation
url: https://www.emergentmind.com/topics/time-fractional-porous-medium-equation
type: topic
---

# Time-Fractional Porous Medium Equation

Searching arXiv for recent and foundational papers on the time-fractional porous medium equation.
The time-fractional porous medium equation denotes a family of nonlinear diffusion equations in which the classical first-order time derivative is replaced by a fractional derivative, typically of Caputo or Riemann–Liouville type. In its prototypical form, it combines porous-medium degeneracy with temporal nonlocality, so that the instantaneous rate of change depends on the full prior history of the solution. Several distinct but related models appear in the literature, including the local-pressure equation
$$
\partial_t^\alpha u=\Delta(u^m),
$$
the divergence-form equation
$$
\partial_t^\alpha u=\partial_x\!\bigl(u^m\,\partial_x u\bigr),
$$
and nonlocal-pressure variants such as
$$
D_t^{\alpha} u - \operatorname{div}\!\bigl(u\,\nabla(-\Delta)^{-\sigma}u\bigr)=f,
$$
with \(0<\alpha<1\) and suitable ranges of \(m\) or \(\sigma\) [1509.06325]. Across these formulations, the central analytical themes are anomalous subdiffusion, self-similarity, free boundaries, compact support versus heavy tails, weak solvability, and regularity.

## 1. Model classes and fractional-time operators

A standard starting point is the replacement of the integer-order derivative in the porous medium equation by a Caputo derivative of order \(0<\alpha<1\). In one widely studied formulation,
$$
{}^{C}\partial_t^\alpha u(x,t)=\partial_x\!\bigl(u^m(x,t)\,\partial_xu(x,t)\bigr),
$$
with \(u(x,0)=0\), where
$$
{}^{C}\partial_t^\alpha u(x,t)=\frac{1}{\Gamma(1-\alpha)}\int_0^t (t-s)^{-\alpha}\,\partial_su(x,s)\,ds
$$
[1409.1825]. Because \(u(x,0)=0\) in that setting, Caputo and Riemann–Liouville derivatives coincide [1409.1825].

A second common local-pressure model is written as
$$
\partial_t^\alpha u(x,t)=\Delta\bigl(u^m(x,t)\bigr)
$$
in \(\mathbb R^d\), with a Dirac-mass initial trace and conserved total mass [2604.09281]. In bounded domains, a more general porous-medium-type problem is
$$
D_t^\alpha u(t,x)=-\mathcal A\bigl(u^m(t,x)\bigr),
$$
where \(\mathcal A\) is a local or nonlocal self-adjoint diffusion operator and zero Dirichlet boundary conditions are imposed [2401.03234].

The literature also contains nonlocal-pressure equations that couple temporal memory with long-range spatial interaction. The model
$$
D_t^{\alpha} u - \operatorname{div}\!\bigl(u\,\nabla(-\Delta)^{-\sigma}u\bigr)=f,
\qquad 0<\sigma<\tfrac12,
$$
uses a Caputo-type time derivative together with the inverse fractional Laplacian
$$
(-\Delta)^{-\sigma}w(x)=C_{n,\sigma}\int_{\mathbb R^n}\frac{w(y)}{|x-y|^{n-2\sigma}}\,dy
$$
[1509.06325]. Another nonlocal-pressure formulation is
$$
\partial_t^\alpha u=\nabla\!\cdot\!\bigl(u^\beta\nabla p\bigr),\qquad p=(-\Delta)^{-s}u,
$$
with \(m=\beta+1>1\) [2409.08441].

The time-fractional porous medium equation is also treated in stochastic and variational settings. On a bounded domain \(\Omega\subset\mathbb R^d\), the stochastic equation
$$
{}^{RL}\!D_t^\beta\bigl(X(t)-x\bigr)=\Delta\Psi(X(t))
+\frac{1}{\Gamma(1+\beta-\gamma)}\int_0^t(t-s)^{\beta-\gamma}B(s)\,dW(s)
$$
includes time-fractional porous media equations as special cases when \(B\equiv0\) and \(\Psi\) is a porous-medium nonlinearity [1708.05649].

These formulations are linked by a common mechanism: the time-fractional derivative introduces memory, while the porous-medium nonlinearity creates degeneracy or singularity according to whether \(m>1\) or \(0<m<1\). This suggests that the term “time-fractional porous medium equation” is best understood as a class of models rather than a single canonical PDE.

## 2. Physical derivation and interpretation

One deterministic derivation begins from Darcy’s law \(q=-D(u)\nabla u\) and mass conservation \(u_t+\nabla\!\cdot q=0\), then replaces the local continuity law by a nonlocal balance
$$
u_t=-\int_0^t w(t-s)\,\nabla\!\cdot q(x,s)\,ds.
$$
Choosing the kernel
$$
w_\alpha(s)=\frac{s^{\alpha-1}}{\Gamma(\alpha)}
$$
leads to
$$
{}^{C}\partial_t^\alpha u=-\nabla\!\cdot q,
$$
and, with \(q=-D(u)\nabla u\), to the time-fractional Richards equation and its porous-medium specialization [1409.1825]. That derivation is presented as dual to the stochastic continuous-time random walk framework and is motivated by waiting-time or trapping effects in porous media [1409.1825].

The interpretation of the Caputo derivative as a memory term is explicit in the nonlocal-pressure model of Allen, Caffarelli, and Vasseur, where the local rate at time \(t\) depends on the full history \(u(s)\) for \(s<t\) [1509.06325]. In that paper, the one-sided integration-by-parts structure of the Caputo derivative is central to the weak formulation and regularity theory [1509.06325].

Applications cited in the literature include subdiffusive moisture uptake in building materials such as brick, cement, and zeolite, hydrocarbons in tight rocks, and contaminant transport in heterogeneous soil [1409.1825]. Experimental support for subdiffusive anomalies is discussed in connection with porous zeolite data showing anomalous scaling \(x/t^{\alpha/2}\) with \(\alpha\approx0.36\), fitted by a three-term approximate profile with parameters
$$
\alpha=0.36,\qquad m=1,\qquad D_0\simeq0.4568\,\mathrm{cm}^2/\mathrm{s}^{0.36}
$$
[1409.1825].

For nonlocal-pressure variants, the inverse fractional Laplacian implies that the pressure at a point depends on remote density values, so the model incorporates both temporal memory and spatial long-range interaction [1509.06325]. A plausible implication is that such equations interpolate between anomalous transport theory and nonlocal aggregation–diffusion dynamics.

## 3. Self-similarity, scaling, and profile equations

Self-similar solutions are a principal analytical tool throughout the subject. For the source-type local equation
$$
\partial_t^\alpha u=\Delta(u^m)
$$
in \(\mathbb R^d\), one seeks
$$
u(x,t)=t^{-\kappa_1}F\!\left(\frac{|x|}{t^{\kappa_2}}\right),
$$
and mass conservation yields
$$
\kappa_1=\frac{\alpha d}{2+d(m-1)},\qquad
\kappa_2=\frac{\alpha}{2+d(m-1)}.
$$
These exponents are positive precisely when
$$
m>m_c:=\frac{(d-2)_+}{d},
$$
which is the optimal range for existence of finite-mass self-similar profiles [2604.09281].

In one dimension for the equation
$$
\partial_t^\alpha u=\partial_x(u^m\partial_xu),
$$
the source-type ansatz
$$
u(x,t)=t^{-a}U(z),\qquad z=x\,t^{-b},
$$
combined with mass conservation gives \(a=b\), and substitution yields
$$
a=b=\frac{\alpha}{m+2}
$$
[2507.19217]. The corresponding classical limit as \(\alpha\to1\) is the Barenblatt–Zel’dovich–Kompaneets solution
$$
B(x,t)=t^{-1/(m+2)}
\Bigl(D-\frac{x^2}{2(m+2)t^{2/(m+2)}}\Bigr)_+^{1/m},
$$
with \(D\) determined by the mass normalization [2507.19217].

Boundary-driven self-similar problems on the half-line produce different exponents because the normalization is set by Dirichlet, Neumann, or Robin data rather than by mass conservation. Under constant-flux boundary conditions,
$$
-u^m(0,t)u_x(0,t)=1,\qquad u(x,0)=0,
$$
Płociniczak obtains
$$
a=\frac{\alpha}{m+2},\qquad b=\frac{\alpha(m+1)}{m+2}
$$
for
$$
u(x,t)=t^aU(\eta),\qquad \eta=x/t^b
$$
[1409.1825]. More generally, for half-line free-boundary problems with Dirichlet, Neumann, or Robin data, the self-similar reduction leads to an ordinary integro-differential equation involving the Erdélyi–Kober operator \(I^{\beta,\gamma}_\delta\) [2106.05138].

The Erdélyi–Kober operator is a recurring structure in the self-similar reduction of time-fractional porous medium equations. In the one-dimensional half-line setting with \(u(x,t)=U(\eta)\), \(\eta=x\,t^{-\alpha/2}\), one obtains equations of the form
$$
(D(U)U')'(\eta)=[A-B\,\eta\,d/d\eta]\,F_\alpha U(\eta),
$$
where
$$
F_\alpha U(\eta)=\frac{1}{\Gamma(1-\alpha)}
\int_0^1(1-s)^{-\alpha}U(s^{-B}\eta)\,ds
$$
[2303.01725]. Equivalent Volterra-type integral equations then follow after integrating twice and using compact support.

For the nonlocal-pressure model, Allen, Caffarelli, and Vasseur prove that if
$$
D_t^\alpha u-\operatorname{div}\!\bigl(u\nabla(-\Delta)^{-\sigma}u\bigr)=f,
$$
then the rescaled function
$$
v(t,x)=A\,u(Bt,Cx)
$$
is again a solution provided
$$
A=B^\alpha C^{2-2\sigma}
$$
[1509.06325]. Their Hölder-regularity argument iterates oscillation decay on nested cylinders at scales synchronized by this relation [1509.06325].

## 4. Existence, uniqueness, and weak-solution frameworks

Existence theory depends strongly on the chosen model.

For the nonlocal-pressure equation
$$
D_t^\alpha u-\operatorname{div}\!\bigl(u\nabla(-\Delta)^{-\sigma}u\bigr)=f
$$
with \(0<\sigma<1/2\), nonnegative data, and exponential decay
$$
0\le u_0(x),\,f(t,x)\le A e^{-|x|},
$$
Allen, Caffarelli, and Vasseur prove existence of weak solutions [1509.06325]. Their construction regularizes \((-\Delta)^{-\sigma}\) by a smooth kernel \(K_\zeta\), adds a viscosity term \(\delta\,\operatorname{div}((u+d)\nabla u)\), discretizes time with step \(\varepsilon\), and passes to the limit using a priori bounds, discrete Aubin–Lions compactness in time, and Besov/Riesz-potential lifting in space [1509.06325].

For the one-dimensional equation
$$
\partial_t^\alpha u=(D(u)u_x)_x,\qquad u(x,0)=0,\qquad u(0,t)=M>0,
$$
with \(D\in C^1(\mathbb R_+)\), \(D(0)=0\), \(D(u)>0\), and \(D'(u)>0\) for \(u>0\), existence and uniqueness of the general initial-boundary-value weak solution are proved in [2303.01725]. The proof uses integration by parts for the Caputo derivative via the right-side fractional integral \(J_t^{1-\alpha}\) and the Ladyženskaja trick applied to the difference of two solutions [2303.01725]. In the self-similar setting, the profile equation is recast as
$$
K(U(\eta))=\int_\eta^{\eta^*}G(\eta,z)\,F_\alpha U(z)\,dz,
$$
with \(K(u)=\int_0^uD(s)\,ds\), and existence follows from a Leray–Schauder fixed-point argument [2303.01725].

A related half-line Dirichlet problem with \(m>1\),
$$
D_t^\alpha u=\partial_x\!\bigl(u^m\partial_xu\bigr),\qquad
u(x,0)=0,\qquad u(0,t)=1,
$$
together with a moving free boundary and no-flux condition, is reduced to a nonlinear Volterra integral equation in \(Y=U^{m+1}\) [1803.03016]. A shooting parameter \(B=-U'(0)\) is introduced; Schauder’s fixed-point theorem yields a compactly supported candidate for each \(B\), and an intermediate-value argument identifies exactly one \(B^*\) producing the no-flux condition [1803.03016].

For broader operator classes on bounded domains, the equation
$$
D_t^\alpha u=-\mathcal A(u^m)
$$
admits a unique nonnegative \(H^*\)-solution for every \(u_0\in L^{1+m}(\Omega)\cap H^*(\Omega)\), assuming the inverse kernel \(G(x,y)\) satisfies
$$
0\le G(x,y)\le C|x-y|^{-(N-2s)}
$$
and the weak Kato inequality holds [2401.03234]. The comparison principle is established in the same general setting [2401.03234].

In the stochastic and variational framework of Gelfand triples, Liu, Röckner, and da Silva prove existence and uniqueness for
$$
{}^{RL}\!D_t^\beta(X(t)-x)=\Delta\Psi(X(t))+\cdots
$$
under assumptions \((\Psi1)\)–\((\Psi4)\) on the monotone nonlinearity \(\Psi\), with the deterministic time-fractional porous-medium equation recovered when the noise vanishes [1708.05649]. The proof combines pseudomonotonicity of the fractional derivative operator with Yosida approximation and Brezis-type surjectivity [1708.05649].

A point of contrast is that uniqueness is open for the fully nonlocal time-and-pressure model of [1509.06325], whereas uniqueness is obtained for several local-pressure one-dimensional and bounded-domain problems [2303.01725].

## 5. Regularity, decay, and qualitative behavior

Regularity theory for time-fractional porous medium equations must address both degeneracy in \(u\) and one-sided temporal nonlocality.

For the equation with fractional potential pressure,
$$
D_t^\alpha u-\operatorname{div}\!\bigl(u\nabla(-\Delta)^{-\sigma}u\bigr)=f,
$$
Allen, Caffarelli, and Vasseur prove that every nonnegative weak solution constructed under the exponential-decay assumptions is locally Hölder continuous in \((t,x)\) [1509.06325]. The proof adapts the De Giorgi–Nash–Moser method to a nonlocal-in-time, nonlocal-in-space, degenerate setting. It uses “pull-up” and “pull-down” lemmas, one-sided integration by parts for the Caputo derivative, mixed space-time energy inequalities, Sobolev embeddings in time and space, and a nonlinear \(L^p\) iteration [1509.06325].

For the nonlocal time porous medium equation
$$
\partial_t^y w+(-\Delta)^s(w^m)=f,
$$
with Marchaud time derivative and \(m>1\), existence and uniqueness of weak solutions are established under
$$
m>m^*=\frac{yN-2s}{yN},\qquad N>\frac{2s}{y},
$$
and bounded weak solutions are shown to be locally Hölder continuous for positive time by a De Giorgi–Nash–Moser argument [1803.03413].

In bounded domains, sharp decay and regularization results are available for
$$
D_t^\alpha u=-\mathcal A(u^m).
$$
Under the assumptions of [2401.03234], there is an \(L^p\to L^\infty\) smoothing estimate
$$
\|u(t)\|_{L^\infty(\Omega)}
\le
K\,\frac{\|u_0\|_{L^p(\Omega)}^{1/m}}{t^{\alpha/m}},
\qquad t>0,
$$
for \(m\ne1\) and suitable \(p\) depending on \(m,s,N\) [2401.03234]. The same paper proves Benilan–Crandall-type monotonicity:
- if \(0<m<1\), then \(t^{1-m}u(t,x)\) is nonincreasing;
- if \(m>1\), then \(t^{m-1}u(t,x)\) is nondecreasing [2401.03234].

The long-time behavior in that bounded-domain framework is particularly notable. For every \(m>0\), there is no finite-time extinction, and every \(L^p\)-norm decays like
$$
\|u(t)\|_{L^p(\Omega)}\sim C(1+t)^{-\alpha/m}
\qquad\text{as }t\to\infty
$$
[2401.03234]. The paper emphasizes that the decay exponent \(\alpha/m\) is independent of \(p\), and that temporal memory mitigates the distinction between slow and fast diffusion regimes [2401.03234].

For one-dimensional Barenblatt profiles of
$$
\partial_t^\alpha u=\partial_x(u^m\partial_xu),
$$
the profile \(U\) is strictly increasing on \((-z_0,0]\), real-analytic on \((-z_0,0)\), globally Hölder continuous with exponent \((1-\alpha)/(m+1)\), and satisfies explicit upper and lower bounds together with free-boundary asymptotics
$$
U(z)\sim C\,(z_0+z)^{(2-\alpha)/m}
\qquad\text{as }z\to -z_0^+
$$
[2507.19217].

These results indicate that the regularity theory bifurcates according to geometry and operator structure: bounded domains favor comparison-based decay estimates, while whole-space self-similar and nonlocal-pressure problems emphasize Hölder continuity and free-boundary structure.

## 6. Barenblatt solutions, compact support, and fast-diffusion profiles

The analogue of the classical Barenblatt solution is central to the subject. For the one-dimensional equation
$$
\partial_t^\alpha u=\partial_x(u^m\partial_xu),
$$
a 2025 construction proves existence of a compactly supported, symmetric mass-one self-similar profile \(U\) with support \([ -z_0^*, z_0^*]\), giving
$$
u(x,t)=t^{-\alpha/(m+2)}U\!\bigl(|x|t^{-\alpha/(m+2)};z_0^*\bigr)
$$
inside the support and zero outside [2507.19217]. The profile is obtained from a nonlinear Volterra integral equation with an Erdélyi–Kober-type kernel expressed באמצעות incomplete Beta functions [2507.19217]. The half-mass condition
$$
\int_{-z_0}^0U(z)\,dz=\frac12
$$
selects a unique \(z_0^*(\alpha,m)\) [2507.19217].

A broader whole-space classification is given for
$$
\partial_t^\alpha u=\Delta(u^m)
$$
in all dimensions \(d\ge1\). There exist self-similar finite-mass solutions for all
$$
m>m_c=\frac{(d-2)_+}{d},
$$
and this range is optimal [2604.09281]. Two regimes are identified:
- for \(m>1\), there is a unique, nonincreasing, compactly supported profile with support \([0,R]\);
- for \(m_c<m<1\), there is a family of positive, strictly decreasing profiles with heavy tails [2604.09281].

In the slow-diffusion range \(m>1\), the profile vanishes at the free boundary according to
$$
F(\xi)\sim C\,(R-\xi)^{\frac{2-\alpha}{m-1}}
\qquad\text{as }\xi\to R^-,
$$
with an explicit coefficient involving Gamma functions [2604.09281]. In the mildly fast-diffusion range \(m_c<m<1\), the asymptotics are
$$
F_M(\xi)=F^*(\xi)\Bigl[1-(L\xi)^{-\gamma^*}+o(\xi^{-\gamma^*})\Bigr],
\qquad
F^*(\xi)=c^*\xi^{-2/(1-m)},
$$
with \(\gamma^*\) defined implicitly by an integral equation involving \(\mathsf Q(s)\) [2604.09281].

The behavior near the similarity origin also depends on dimension. As \(\xi\to0^+\),
- in \(d\ge3\), \(U(\xi)\to\infty\) with power-law blow-up;
- in \(d=2\), \(U(\xi)\) grows like a logarithm;
- in \(d=1\), \(U(\xi)\) tends to a constant
[2604.09281]. The authors describe this as reflecting a “Newtonian-potential-type” memory at the origin [2604.09281].

These results clarify a frequent misconception inherited from the classical theory. In the time-fractional setting, compact support is not universal: it holds in the slow-diffusion regime \(m>1\), but for \(m_c<m<1\) the finite-mass self-similar solutions have heavy tails rather than compact support [2604.09281].

## 7. Symmetry analysis and numerical methodologies

Lie-symmetry methods provide explicit reductions for the one-dimensional fractional porous medium equation
$$
\partial_t^\alpha u=(u^m)_{xx},
\qquad 0<\alpha<1,\quad m\ne1,
$$
when the time derivative is taken in the Riemann–Liouville sense [1901.03666]. The admitted Lie algebra is generated by
$$
V_{11}=t\partial_t-\frac{\alpha}{r-1}u\partial_u,\qquad
V_{12}=x\partial_x+\frac{2}{r-1}u\partial_u,\qquad
V_{13}=\partial_x,
$$
with \([V_{12},V_{13}]=-V_{13}\) [1901.03666]. An optimal one-dimensional subalgebra system is computed, similarity reductions are derived, and several invariant solutions are exhibited, including
$$
u(t)=C\,t^{\alpha-1}
$$
for the translation-invariant case \(V_{13}\) [1901.03666].

Numerically, much of the literature exploits the self-similar reduction to convert the PDE into a Volterra integral equation or an EK-based integral equation. For the half-line equation with Dirichlet data, the nonlinear Abel-type Volterra equation
$$
y(z)^{m+1}=\int_0^zK(z,u)\,y(u)\,du
$$
admits a family of explicit discretizations, including second-order trapezoidal schemes [2106.05138]. The main convergence theorem states that if
$$
\mu_m=\frac{4K_+}{(m+1)V_-}<3,
$$
then the scheme converges with order
$$
\min\{2,\,3-\mu_m\}
$$
[2106.05138].

For the EK operator itself, rectangle and trapezoid discretizations satisfy error bounds
$$
\|F_\alpha U-\widehat F_{\alpha,N}U\|_\infty
\le
\frac{C_1h}{\Gamma(2-\alpha)}
\quad\text{or}\quad
\frac{C_2h^2}{\Gamma(2-\alpha)},
$$
with optimal truncations
$$
N\approx h^{-B}\quad\text{or}\quad N\approx h^{-2B}
$$
respectively [2303.01725]. In the same work, the rectangle backward-shooting scheme is explicit and yields an \(O(N^2)\) algorithm [2303.01725].

For the Barenblatt integral equation in one dimension, product-integration on a uniform mesh yields a discrete fixed-point system
$$
U_n=
\Bigl[h\sum_{i=1}^{n-1}w_{n,i}U_i\Bigr]^{1/(m+1)},
$$
with the initial step selected from the free-boundary asymptotics [2507.19217]. If the quadrature error is \(O(\delta(h))\), then the global error is
$$
O\!\bigl(\ln\ln(1/h)\,\ln(1/h)\,\delta(h)\bigr)
$$
[2507.19217].

A distinct numerical and variational direction is the JKO formulation for the nonlocal-pressure equation
$$
\partial_t^\alpha u=\nabla\!\cdot(u^\beta\nabla(-\Delta)^{-s}u),
$$
where the discrete dynamics minimize
$$
\rho^k\in\operatorname{Arg\,min}_{\rho\in P_2(\mathbb R^d)}
\left\{
\frac{1}{2\tau}W_m(\rho,\rho^{k-1})^2+\mathcal F(\rho)
\right\},
$$
with modified Wasserstein distance \(W_m\) and energy
$$
\mathcal F(\rho)=\tfrac12\|\rho\|_{\dot H^{-s}(\mathbb R^d)}^2
$$
[2409.08441]. Under \(u_0\in P_2(\mathbb R^d)\cap\dot H^{-s}(\mathbb R^d)\), \(d\ge2\), \(\beta\in(0,1)\), \(\alpha\in(0,\min\{1,2\})\), and \(s\in(0,\min\{1,d/2\})\), the scheme converges to a global weak solution [2409.08441].

The numerical literature therefore reflects two dominant strategies. One reduces to self-similar integral equations and exploits their one-dimensional weakly singular structure; the other builds variational time discretizations directly at the PDE level. This suggests that the computational treatment of time-fractional porous medium equations is shaped as much by their memory structure as by their nonlinear degeneracy.

Source: https://www.emergentmind.com/topics/time-fractional-porous-medium-equation