alphaFeedback

WorksMathematics

Energy Solutions and Harnack Estimates for Nonlocal Drift–Diffusion

AI contributions
Not recorded for this paper
View details

Abstract

We consider nonlocal drift–diffusion equations of order 2s∈[1,2)2s\in[1,2) with bounded divergence-free drift. The symmetric measurable jump kernel satisfies local coercivity, a global upper form bound, and a scale-invariant bound on the intensity of long jumps. The transport form is bounded and skew-symmetric on H1/2×H1/2H^{1/2}\times H^{1/2}. In the fractional energy space, distributional and variational solutions coincide and satisfy nonlinear testing identities. Viscous approximations converge strongly in CtLx2∩Lt2HxsC_tL^2_x\cap L^2_tH^s_x to the unique solution of the Cauchy problem. For sources bounded in the cylinder, supersolutions nonnegative there satisfy a weak Harnack inequality involving the time integral of the negative exterior tail. Under this source bound, locally bounded solutions whose exterior tail belongs to LqL^q in time for some q>1q>1 are locally Hölder continuous. This also gives Hölder continuity at positive times for homogeneous global solutions with L2L^2 data. Under upper jumping smoothness, control of the positive tail on an earlier time interval yields the full parabolic Harnack inequality for solutions nonnegative on the full spatial slab. A kernel supported on alternating dyadic annuli satisfies all the kernel assumptions. The local estimates depend on the drift through R2s−1∥b∥∞R^{2s-1}\|b\|_\infty.

Introduction

Consider the equation

∂tu+b(t,x)⋅∇u+Ltu=f,(1)\partial_t u+b(t,x)\cdot\nabla u+\mathcal{L}_t u=f, \tag*{(1)}

where bb is a bounded divergence-free vector field and Lt\mathcal{L}_t is a nonlocal diffusion operator of order α=2s∈[1,2)\alpha=2s\in[1,2). Its symmetric jump kernel may depend measurably on time and need not have a pointwise lower bound of fractional Laplacian type. At the critical order α=1\alpha=1, the drift and diffusion have the same scaling, while the energy controls only half a spatial derivative. To compare a solution on an earlier cylinder with its infimum on a later one, one must also account for jumps from outside the region where the equation holds.

In the absence of drift, Felsinger and Kassmann [10] proved a weak Harnack inequality and Hölder estimates under quadratic-form assumptions. They also discussed first-order terms for α>1\alpha>1. Strömqvist [14] obtained Harnack estimates that allow negative values outside the cylinder. Kassmann and Weidner [11] proved a full parabolic Harnack inequality using time-integrated tails and distinguished the roles of local coercivity and upper jumping smoothness. In the probabilistic setting, Bass and Levin [1] proved a Harnack inequality for jump processes. Chen, Kumagai, and Wang [5] identified upper jumping smoothness as a necessary condition for the parabolic Harnack inequality for autonomous symmetric jump forms.

For fractional diffusion with drift, Caffarelli and Vasseur [4] proved Hölder regularity at the critical order in their study of the quasi-geostrophic equation. Chen and Zhang [6] obtained Harnack estimates for the fractional Laplacian with critical and supercritical divergence-free drifts. Liao [12] proved weak and full Harnack inequalities and Hölder estimates for pointwise comparable kernels and divergence-free drifts in scale-invariant classes.

The definition of a local solution there requires a locally square-integrable spatial gradient [12 Definition 1.1]. Ciani and Nakamura [7] derived regularity estimates for functions in nonlocal parabolic De Giorgi classes.

For bounded incompressible drift, interpolation gives a continuous transport form on H1/2×H1/2H^{1/2} \times H^{1/2}. Hence the form acts on the diffusion energy space even at s=1/2s = 1/2. This bound justifies nonlinear testing in § 2. The energy identity then gives strong convergence of the viscous approximations, while a fractional Nash inequality gives positive-time smoothing (§ 3).

For the local estimates, we add the time integral of the negative exterior tail and the negative source bound to the solution before taking powers. The derivative of this shift absorbs both terms in the energy inequality. The localized energy bounds in § 4 and a differential inequality for the weighted mean of the logarithm give the weak Harnack inequality in § 5. In § 6 we use it to reduce the oscillation, estimating the tail on successive annuli. For the full Harnack inequality, upper jumping smoothness bounds the earlier positive tail in terms of the later infimum, the source, and an arbitrarily small multiple of a supremum on a larger past cylinder. Iteration on nested cylinders absorbs the supremum (§ 7).

The examples in § 8 distinguish the kernel assumptions and illustrate the time dependence of the exterior tail. In Example 8.1 the kernel is supported on a union of separated dyadic annuli. Each missing jump splits into two allowed jumps inside the same ball, giving local coercivity. This kernel also satisfies upper jumping smoothness, even after multiplication by a symmetric measurable factor bounded above and bounded away from zero. In Example 8.3, concentrated jumps preserve the form and jump-intensity bounds but violate upper jumping smoothness.

Kernel assumptions and notation

For a ball D⊂RnD \subset\mathbb{R}^{n}, write

[v]s,D2=∬D×D∣v(x)−v(y)∣2∣x−y∣n+2s dx dy,[v]s=[v]s,Rn.[v]_{s,D}^{2} = \iint_{D \times D} \frac{\lvert v(x) - v(y)\rvert^{2}}{\lvert x-y\rvert^{n+2s}}\,\mathrm{d}x\,\mathrm{d}y, \qquad[v]_{s} = [v]_{s,\mathbb{R}^{n}}.

The norm on Hs(Rn)H^{s}(\mathbb{R}^{n}) is (∥v∥22+[v]s2)1/2(\lVert v\rVert_{2}^{2} + [v]_{s}^{2})^{1/2}. Set

Et,D(v,w)=12∬D×D(v(x)−v(y))(w(x)−w(y))K(t;x,y) dx dy,Et=Et,Rn.(2)\mathcal{E}_{t,D}(v,w) = \frac{1}{2}\iint_{D \times D} (v(x)-v(y))(w(x)-w(y))K(t;x,y)\,\mathrm{d}x\,\mathrm{d}y, \qquad\mathcal{E}_{t} = \mathcal{E}_{t,\mathbb{R}^{n}}. \tag*{(2)}

The kernel KK is nonnegative, jointly measurable, and symmetric in x,yx,y. Fix 0<λ≤Λ0 < \lambda\leq\Lambda. On a common set of times of full measure assume

Et,D(v,v)≥λ[v]s,D2for every ball D and v∈Hs(D),(3)\mathcal{E}_{t,D}(v,v) \geq\lambda[v]_{s,D}^{2} \qquad\text{for every ball } D \text{ and } v \in H^{s}(D), \tag*{(3)}
Et(v,v)≤Λ[v]s2for every v∈Hs(Rn),(4)\mathcal{E}_{t}(v,v) \leq\Lambda[v]_{s}^{2} \qquad\text{for every } v \in H^{s}(\mathbb{R}^{n}), \tag*{(4)}
ess supx∫∣x−y∣≥rK(t;x,y) dy≤Λr−αfor every r>0.(5)\mathop{\mathrm{ess\,sup}}_{x}\int_{\lvert x-y\rvert\geq r} K(t;x,y)\,\mathrm{d}y \leq\Lambda r^{-\alpha} \qquad\text{for every } r > 0. \tag*{(5)}

In (1.5), it suffices to take positive rational radii on a common set of spatial points of full measure. Monotone convergence then gives the estimate at every radius. We assume

b∈L∞((0,T)×Rn;Rn),div⁡xb=0in distributions.(6)b \in L^{\infty}((0,T) \times\mathbb{R}^{n};\mathbb{R}^{n}), \qquad\operatorname{div}_{x} b = 0 \quad\text{in distributions}. \tag*{(6)}

All spatial suprema and infima of measurable functions are essential extrema.

For z0=(t0,x0)z_{0} = (t_{0},x_{0}), define the ambient cylinder CRC_{R}, the past and future cylinders PRP_{R} and FRF_{R}, and the smaller past cylinder HRH_{R} by

CR(z0)=(t0−4Rα,t0+Rα)×B4R(x0),PR(z0)=(t0−3Rα,t0−Rα)×B2R(x0),FR(z0)=(t0,t0+Rα)×BR(x0),HR(z0)=(t0−52Rα,t0−2Rα)×BR(x0).(7)\begin{aligned} C_{R}(z_{0}) &= (t_{0}-4R^{\alpha},t_{0}+R^{\alpha}) \times B_{4R}(x_{0}),\\ P_{R}(z_{0}) &= (t_{0}-3R^{\alpha},t_{0}-R^{\alpha}) \times B_{2R}(x_{0}),\\ F_{R}(z_{0}) &= (t_{0},t_{0}+R^{\alpha}) \times B_{R}(x_{0}),\\ H_{R}(z_{0}) &= (t_{0}-\tfrac{5}{2}R^{\alpha},t_{0}-2R^{\alpha}) \times B_{R}(x_{0}). \tag*{(7)} \end{aligned}

We take cylinders whose closures lie in the solution domain, except possibly for the terminal time face. The rescaled drift has norm

MR=Rα−1∥b∥L∞(CR).(8)M_R = R^{\alpha-1}\|b\|_{L^\infty(C_R)}. \tag*{(8)}

Structural constants depend only on nn, ss, λ\lambda, Λ\Lambda and, in local estimates, on an upper bound for the rescaled drift norm. On CRC_R this norm is MRM_R, and on a backward cylinder QRQ_R it is Rα−1∥b∥L∞(QR)R^{\alpha-1}\|b\|_{L^\infty(Q_R)}. We indicate any further dependence.

For concentric balls with r<ρr<\rho, define

Θ(v;t,r,ρ)=ess sup⁡x∈Br(x0)∫Rn∖Bρ(x0)∣v(t,y)∣K(t;x,y) dy.(9)\Theta(v;t,r,\rho)=\operatorname*{ess\,sup}_{x\in B_r(x_0)}\int_{\mathbb{R}^n\setminus B_\rho(x_0)}|v(t,y)|K(t;x,y)\,dy. \tag*{(9)}

Here and below we omit the common center when it is clear. The weak Harnack inequality involves the negative tail

TR−(u)=∫t0−4Rαt0+RαΘ(u−;t,3R,4R) dt.(10)T_R^-(u)=\int_{t_0-4R^\alpha}^{t_0+R^\alpha}\Theta(u_-;t,3R,4R)\,dt. \tag*{(10)}

Under parabolic rescaling, TR−(u)T_R^-(u) has the same scaling as uu. For a bounded function vv, the row bound (5) gives

Θ(v;t,r,ρ)≤Λ(ρ−r)−α∥v(t)∥∞.(11)\Theta(v;t,r,\rho)\leq\Lambda(\rho-r)^{-\alpha}\|v(t)\|_\infty. \tag*{(11)}

Principal results

Write

XT=L∞(0,T;L2(Rn))∩L2(0,T;Hs(Rn)).\mathcal{X}_T=L^\infty(0,T;L^2(\mathbb{R}^n))\cap L^2(0,T;H^s(\mathbb{R}^n)).

For u0∈L2u_0\in L^2 and f∈L2(0,T;H−s(Rn))f\in L^2(0,T;H^{-s}(\mathbb{R}^n)), a distributional energy solution is a function u∈XTu\in\mathcal{X}_T such that

−∫0T∫u ∂tϕ dx dt−∫0T∫u b⋅∇ϕ dx dt+∫0TEt(u,ϕ) dt=∫0T⟨f,ϕ⟩ dt+∫u0ϕ(0) dx(12)-\int_0^T\int u\,\partial_t\phi\,dx\,dt-\int_0^T\int u\,b\cdot\nabla\phi\,dx\,dt+\int_0^T\mathcal{E}_t(u,\phi)\,dt =\int_0^T\langle f,\phi\rangle\,dt+\int u_0\phi(0)\,dx \tag*{(12)}

for every ϕ∈Cc∞([0,T)×Rn)\phi\in C^\infty_c([0,T)\times\mathbb{R}^n). Definition 4.2 specifies the local variational inequalities. For distributional energy solutions the equation implies the local time regularity required there.

References

  1. [1]R. F. Bass and D. A. Levin, Harnack inequalities for jump processes, Potential Anal. 17 (2002), no. 4, 375–388. doi:10.1023/A:1016378210944.DOI
  2. [2]J. Bergh and J. Löfström, Interpolation Spaces. An Introduction, Grundlehren der mathematischen Wissenschaften, vol. 223, Springer, Berlin, 1976. doi:10.1007/978-3-642-66451-9.
  3. [3]E. Bombieri and E. Giusti, Harnack’s inequality for elliptic differential equations on minimal surfaces, Invent. Math. 15 (1972), 24–46. doi:10.1007/BF01418921.
  4. [4]L. A. Caffarelli and A. Vasseur, Drift diffusion equations with fractional diffusion and the quasi-geostrophic equation, Ann. of Math. (2) 171 (2010), no. 3, 1903–1930. doi:10.4007/annals.2010.171.1903.DOI
  5. [5]Z.-Q. Chen, T. Kumagai, and J. Wang, Stability of parabolic Harnack inequalities for symmetric non-local Dirichlet forms, J. Eur. Math. Soc. 22 (2020), no. 11, 3747–3803. doi:10.4171/JEMS/996.DOI
  6. [6]Z.-Q. Chen and X. Zhang, Harnack inequalities for nonlocal operators with supercritical drifts and their applications, preprint, arXiv:2511.13012, 2025.arxiv.org/abs/2511.13012
  7. [7]S. Ciani and K. Nakamura, Nonlocal parabolic De Giorgi classes, preprint, arXiv:2508.16247v4, 2026.arxiv.org/abs/2508.16247
  8. [8]E. Di Nezza, G. Palatucci, and E. Valdinoci, Hitchhiker’s guide to the fractional Sobolev spaces, Bull. Sci. Math. 136 (2012), no. 5, 521–573. doi:10.1016/j.bulsci.2012.11.004.
  9. [9]B. Dyda and M. Kassmann, On weighted Poincaré inequalities, Ann. Acad. Sci. Fenn. Math. 38 (2013), no. 2, 721–726. doi:10.5186/aasfm.2013.3834.
  10. [10]M. Felsinger and M. Kassmann, Local regularity for parabolic nonlocal operators, Comm. Partial Differential Equations 38 (2013), no. 9, 1539–1573. doi:10.1080/03605302.2013.808211.DOI
  11. [11]M. Kassmann and M. Weidner, The parabolic Harnack inequality for nonlocal equations, Duke Math. J. 173 (2024), no. 17, 3413–3451. doi:10.1215/00127094-2024-0008.DOI
  12. [12]N. Liao, Harnack estimates for nonlocal drift-diffusion equations, preprint, arXiv:2402.11986v2, 2025.arxiv.org/abs/2402.11986
  13. [13]J.-L. Lions, Équations différentielles opérationnelles et problèmes aux limites, Grundlehren der mathematischen Wissenschaften, vol. 111, Springer, Berlin, 1961.
  14. [14]M. Strömqvist, Harnack’s inequality for parabolic nonlocal equations, Ann. Inst. H. Poincaré C Anal. Non Linéaire 36 (2019), no. 6, 1709–1745. doi:10.1016/j.anihpc.2019.03.003.DOI

Paper details

vlnc:W_fu2jaqgd
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents