Parabolic Double Phase Equation

Updated 9 July 2026

The parabolic double phase equation is a nonlinear parabolic PDE featuring two competing growth regimes (p-phase and q-phase) modulated by a coefficient a(x,t).
Methodologies involve phase analysis, intrinsic scaling, and energy estimates like Caccioppoli inequalities to establish local boundedness and regularity.
Implications include refined weak and strong solution frameworks, enhanced interior and boundary regularity, and extended applications in variable-growth problems.

The parabolic double phase equation is a nonlinear parabolic PDE whose prototype is

$u_t-\operatorname{div}\!\left(|\nabla u|^{p-2}\nabla u+a(x,t)|\nabla u|^{q-2}\nabla u\right)=0,$

with energy density

$|\nabla u|^p+a(x,t)|\nabla u|^q.$

It is called “double phase” because the diffusion has two competing power-growth regimes: a $p$ -phase and a $q$ -phase activated by the nonnegative coefficient $a(x,t)$ . Where $a(x,t)=0$ , the equation reduces to the $p$ -parabolic equation; where $a(x,t)$ is positive, the $q$ -growth part becomes active. In this sense the medium is heterogeneous in both the size of $\nabla u$ and the space-time location $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 0 (Adimurthi et al., 2020).

1. Canonical form and structural regimes

The standard model is

$|\nabla u|^p+a(x,t)|\nabla u|^q.$ 1

A recurrent structural theme is that the equation is locally similar either to a $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 8-Laplacian flow or to a $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 9-Laplacian flow. In the Hölder-regular coefficient theory this is formalized by a phase analysis. In the $p$ 0-phase, the $p$ 1-term dominates in the intrinsic scaling; in the $p$ 2-phase, the $p$ 3-term is subordinate and the estimates reduce to a $p$ 4-Laplacian type argument (Kim et al., 2024). This local phase classification is not merely heuristic: it determines the intrinsic cylinders, the Sobolev–Poincaré inequalities, and the admissible gap conditions between $p$ 5 and $p$ 6.

The same double-phase mechanism appears in variable-growth problems. A representative variable-exponent model is

$p$ 7

with Lipschitz-continuous exponents and coefficients, together with the balance conditions

$p$ 8

(Arora et al., 2021). In that setting, coercivity is preserved because at each point at least one phase is active.

2. Weak and strong solution frameworks

For the prototype constant-exponent problem, weak solutions are the natural solution concept. One formulation requires

$p$ 9

and the integral identity

$q$ 0

for all $q$ 1, where

$q$ 2

(Kim et al., 17 Nov 2025). In more general inhomogeneous problems one replaces the right-hand side by divergence-form data involving a field $q$ 3 and the same double-phase structure (Kim et al., 4 Jan 2026, Kim et al., 2023).

Several works use stronger ambient spaces, for example

$q$ 4

especially when proving Hölder continuity or boundedness under $q$ 5-growth assumptions (Kim et al., 2024, Kim et al., 16 Apr 2026). A separate technical point is that parabolic Lipschitz truncation can recover the more natural energy space without assuming $q$ 6 a priori (Kim et al., 2023).

Strong solutions arise in the variable-growth theory. For the Dirichlet problem with

$q$ 7

one obtains existence and uniqueness of strong solutions with

$q$ 8

and global higher integrability

$q$ 9

under a small-gap condition $a(x,t)$ 0 with $a(x,t)$ 1 (Arora et al., 2020). Related strong-solution theories with global higher integrability and second-order spatial regularity have also been proved for variable-growth double phase fluxes (Arora et al., 2021, Arora et al., 5 Jul 2026).

3. Interior regularity theory

The basic interior regularity problem is local boundedness. For the prototype equation with measurable $a(x,t)$ 2, one optimal boundedness result proves that weak solutions are locally bounded when

$a(x,t)$ 3

(Adimurthi et al., 2020). A later bounded-coefficient result establishes local boundedness under

$a(x,t)$ 4

using a Caccioppoli inequality, a parabolic embedding theorem, and an iteration method (Kim et al., 16 Apr 2026). These theorems address different structural regimes and should not be conflated into a single universal threshold.

For bounded weak solutions with Hölder-continuous coefficients, local Hölder continuity is now available in the degenerate case. If

$a(x,t)$ 5

then bounded weak solutions are locally Hölder continuous, and the proof is based on phase analysis and methods for the $a(x,t)$ 6-Laplace equation (Kim et al., 2024). In that framework, the phase analysis determines whether the equation is locally similar to the $a(x,t)$ 7-Laplace or the $a(x,t)$ 8-Laplace equation.

A further step is gradient higher integrability. For bounded solutions to

$a(x,t)$ 9

with $a(x,t)=0$ 0, one obtains local higher integrability of

$a(x,t)=0$ 1

under the gap condition

$a(x,t)=0$ 2

More generally, if

$a(x,t)=0$ 3

then higher integrability holds under

$a(x,t)=0$ 4

which interpolates between the $a(x,t)=0$ 5 threshold and the bounded-solution threshold (Kim et al., 17 Nov 2025). In the singular regime $a(x,t)=0$ 6, a parallel interpolation theory gives higher integrability under

$a(x,t)=0$ 7

for bounded solutions, or

$a(x,t)=0$ 8

for $a(x,t)=0$ 9 (Kim et al., 4 Jan 2026).

Regularity beyond Hölder continuity has also been obtained. For bounded continuous weak solutions with $p$ 0, bounded $p$ 1, local spatial Lipschitz continuity of $p$ 2, and continuity in time, local spatial Lipschitz regularity and time Hölder regularity follow for the equation with gradient-dependent forcing

$p$ 3

In the degenerate regime $p$ 4, the time exponent improves to the sharp $p$ 5 (Sen et al., 22 Aug 2025).

4. Boundary behavior, Harnack theory, and pointwise continuity estimates

Boundary regularity in the parabolic double-phase setting is phase-dependent. For bounded weak solutions of

$p$ 6

with prototype

$p$ 7

a Wiener-type sufficient criterion describes continuity up to the lateral boundary. If $p$ 8, regularity is governed by the $p$ 9-capacity; if $a(x,t)$ 0, it is governed by the $a(x,t)$ 1-capacity (Ciani et al., 2024). This phase-sensitive criterion leads, under uniform fatness or density assumptions, to boundary Hölder continuity.

Weak Harnack theory also reflects the double-phase threshold phenomena. For non-negative super-solutions of

$a(x,t)$ 2

with $a(x,t)$ 3, a weak Harnack inequality is proved in the case

$a(x,t)$ 4

under the additional assumption

$a(x,t)$ 5

(Savchenko et al., 2023). The need for extra integrability is a structural feature of the supercritical $a(x,t)$ 6-regime.

A different continuity theory controls the oscillation of bounded weak solutions by nonlinear potentials of the data. For

$a(x,t)$ 7

with

$a(x,t)$ 8

and assumptions

$a(x,t)$ 9

one obtains local continuity estimates in terms of elliptic Riesz potentials $q$ 0 and $q$ 1 (Li, 3 Feb 2025). This provides a potential-theoretic continuity criterion adapted to the degenerate double-phase parabolic structure.

5. Main variants and generalizations

The parabolic double phase literature now includes several distinct model classes, each with its own analytic thresholds.

Model class	Representative form	Representative result
Variable-growth local equations	$q$ 2	Strong solutions, global higher integrability, second-order regularity (Arora et al., 2021)
Variable-exponent irregular evolution	$q$ 3	Global Calderón–Zygmund transfer, higher integrability, second-order space regularity (Arora et al., 5 Jul 2026)
Nonlocal double phase equations	Fractional $q$ 4- and $q$ 5-phase integral operator with bounded measurable $q$ 6	Local boundedness of variational solutions (Prasad et al., 2021)
Mixed local/nonlocal equations	$q$ 7	Local boundedness; semicontinuity of super- and subsolutions (Shang et al., 2023)
Normalized double phase flow	$q$ 8	Asymptotic mean value characterization in the viscosity sense (Meng et al., 2022)
Parabolic systems of double phase type	$q$ 9	Partial regularity: $\nabla u$ 0 locally Hölder continuous except on a set of measure zero (Ok et al., 4 Oct 2025)

For systems, a recent partial regularity theorem proves that if

$\nabla u$ 1

then the spatial gradient of any weak solution is locally Hölder continuous except on a set of measure zero (Ok et al., 4 Oct 2025). This is a genuinely parabolic system result, not merely a scalar extension.

The existence, uniqueness, and energy theory for inhomogeneous Dirichlet problems has also been strengthened by a parabolic double-phase Lipschitz truncation method based on a Whitney-type covering result and a related partition of unity in the intrinsic geometry for the double-phase problem (Kim et al., 2023).

6. Analytic mechanisms, gap conditions, and terminological distinctions

The dominant analytic tools are Caccioppoli inequalities, intrinsic cylinders, reverse Hölder inequalities, De Giorgi iteration, Whitney-type coverings, and covering arguments of Vitali type (Kim et al., 2023, Kim et al., 17 Nov 2025). In the singular higher-integrability theory, the proof uses intrinsic cylinders adapted to the singular scaling, a stopping-time argument splitting into $\nabla u$ 2-phase and $\nabla u$ 3-phase cylinders, Caccioppoli inequalities and Poincaré/Gagliardo–Nirenberg inequalities, a proof that the “bad” third case cannot occur thanks to the gap condition, a reverse Hölder inequality for $\nabla u$ 4, and a Vitali covering argument plus Fubini to deduce higher integrability (Kim et al., 4 Jan 2026).

Gap conditions are a central organizing principle, but the literature does not support a single universal “optimal” formula valid across all settings. Bounded measurable coefficients yield local boundedness under one class of thresholds (Adimurthi et al., 2020, Kim et al., 16 Apr 2026); Hölder-continuous coefficients yield Hölder continuity and higher integrability under $\nabla u$ 5 or its interpolative refinements (Kim et al., 2024, Kim et al., 17 Nov 2025); variable-exponent strong-solution theories impose balance conditions such as $\nabla u$ 6 (Arora et al., 2021); and the coefficient class $\nabla u$ 7 leads to the parabolic gap bound

$\nabla u$ 8

for Hölder continuous weak solutions (Sen et al., 9 Jun 2026). That last condition is explicitly described as purely parabolic in nature and stricter than the optimal gap relation associated with the Lavrentiev phenomenon for the elliptic double phase functional (Sen et al., 9 Jun 2026).

A further conceptual distinction concerns terminology. In the standard nonlinear-parabolic literature, “double phase” refers to $\nabla u$ 9- and $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 00-growth modulated by a nonnegative coefficient $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 01. A different line of work uses “two-phase” for one-dimensional forward-backward diffusion laws with two forward-diffusion phases separated by an unstable interval. There the model equation is

$|\nabla u|^p+a(x,t)|\nabla u|^q.$ 02

or, in the simpler entropy framework,

$|\nabla u|^p+a(x,t)|\nabla u|^q.$ 03

and the emphasis is on phase transitions, hysteresis, nonuniqueness, and generalized Abel equations rather than on $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 04-growth regularity (Kim et al., 2016, Terracina, 2013). The shared vocabulary reflects phase coexistence, but the PDE structures are distinct.

Taken together, these developments show that the parabolic double phase equation is not a single theorem or a fixed regularity paradigm. It is a class of nonuniformly elliptic or parabolic evolution equations in which the interaction between the growth gap $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 05, the regularity or geometry of the coefficient $|\nabla u|^p+a(x,t)|\nabla u|^q.$ 06, and the intrinsic parabolic scaling determines whether one can prove local boundedness, Hölder or Lipschitz continuity, higher integrability, strong solvability, boundary regularity, or only partial regularity.