Hölder regularity of parabolic equations with Dirichlet boundary conditions and application to reaction-diffusion and reaction-cross-diffusion systems

基本信息

• 论文ID: 2511.18872
• 标题: Hölder regularity of parabolic equations with Dirichlet boundary conditions and application to reaction-diffusion and reaction-cross-diffusion systems
• 作者: Hector Bouton, Laurent Desvillettes, Helge Dietert (Université Paris Cité and Sorbonne Université, CNRS, IMJ-PRG)
• 分类: math.AP (偏微分方程分析)
• 提交时间: 2025年11月24日
• 论文链接: https://arxiv.org/abs/2511.18872v1

摘要

本文将作者最近关于Neumann边界条件的工作推广到Dirichlet边界条件情形。核心研究对象是具有粗糙系数$a$的抛物方程$a\partial_t w - \Delta w = f$，配以齐次Dirichlet边界条件和特殊假设$\partial_t w \geq 0$。基于此结果，证明了三维空间中带Lotka-Volterra反应项的三角Shigesada-Kawasaki-Teramoto (SKT)交叉扩散系统全局强解的存在性，以及可逆化学反应-扩散系统解的估计，均在Dirichlet边界条件下成立。

研究背景与动机

问题背景

1. 核心问题：研究具有非常数粗糙系数的抛物方程在Dirichlet边界条件下的Hölder正则性。具体考虑方程： $a(t,x)\partial_t w(t,x) - \Delta w(t,x) = f(t,x), \quad \partial_t w \geq 0$ 其中$a$满足椭圆性条件$0 < a_0 \leq a \leq c_0 a_0 < \infty$。
2. 重要性：
   • Hölder正则性是偏微分方程理论的基本问题，对理解解的性质至关重要
   • 该结果可应用于生物数学、化学反应动力学等领域的非线性系统
   • 为交叉扩散系统和反应-扩散系统的全局存在性提供关键技术工具
3. 现有方法局限性：
   • Krylov-Safonov (1981)的经典结果虽然涵盖更一般情形，但常数不显式且需要更强的空间假设
   • 作者之前的工作boutondesvillettes2025仅处理Neumann边界条件
   • Dirichlet边界条件下，解在边界附近的行为更复杂，需要新的技术处理
4. 研究动机：
   • 利用单调性假设$\partial_t w \geq 0$简化证明，获得显式常数
   • 仅需临界Lebesgue空间中对$f$的控制
   • 完整处理边界条件，为应用提供明确的理论支撑

核心贡献

1. 主要理论结果（定理1）：证明了在条件$\gamma := 2 - \frac{2}{p} - \frac{d}{q} > 0$下，方程(1)的非负解$w$属于$C^{0,\alpha}([0,T]\times\Omega)$，且有估计： $\|w\|_{C^{0,\alpha}([0,T]\times\Omega)} \leq C_*\left(\|f\|_{L^p((0,T];L^q(\Omega))} + \|w^{in}\|_{C^1(\Omega)}\right)$
2. 应用于可逆化学系统（命题1）：证明了四维及以下空间中，可逆化学反应系统 $\partial_t u_i - d_i\Delta u_i = (-1)^i(u_1u_3 - u_2u_4), \quad i=1,\ldots,4$ 在Dirichlet边界条件下强解的全局存在性
3. 应用于SKT交叉扩散系统（命题2）：证明了四维及以下空间中，三角SKT系统\partial_t u - \Delta[(d_1+\sigma v)u] = f_u(u,v)\\ \partial_t v - d_2\Delta v = f_v(u,v) \end{cases}$$ 在Dirichlet边界条件下强解的全局存在性
4. 技术创新：发展了专门适用于Dirichlet边界条件的振荡衰减技术，与Neumann情形有本质区别

方法详解

任务定义

输入：

• 有界$C^2$区域$\Omega \subset \mathbb{R}^d$
• 满足椭圆性条件的系数$a(t,x)$
• 强迫项$f \in L^p((0,T];L^q(\Omega))$，其中$\gamma = 2-\frac{2}{p}-\frac{d}{q} > 0$
• 初值$w^{in} \in C^1(\Omega)$满足$w^{in}|_{\partial\Omega} = 0$

输出：

• 解$w$的Hölder连续性：$w \in C^{0,\alpha}([0,T]\times\Omega)$
• 显式的Hölder范数估计

约束条件：

• 齐次Dirichlet边界条件：$w = 0$ on $(0,T]\times\partial\Omega$
• 单调性：$\partial_t w \geq 0$
• 非负性：$w \geq 0$

核心技术框架

1. Dirichlet热核估计（引理1-2）

引理2（关键下界估计）：对于球$B(0,R)$中的Dirichlet热核$\Gamma_{B(0,R),D}$，取$T_R := a_0\frac{9R^2}{32d}$，则 $\inf_{x,y\in B(0,R/4)} \inf_{t\in[\frac{T_R}{2c_0a_0}, \frac{T_R}{a_0}]} \Gamma_{B(0,R),D}(t,x,y) \gtrsim_{c_0,d} R^{-d}$

证明思路：

• 构造比较函数$\Psi(t,x,y) := p(t,|x-y|) - \sup_{0\leq s\leq t}p(s,3R/4)$
• 利用最大值原理得到$\Gamma_{B(0,R),D}(t,x,y) > \Psi(t,x,y)$
• 对于$x,y\in B(0,R/4)$和适当的时间范围，精确计算得到正的下界

2. 振荡衰减（命题3）

核心思想：在远离边界的区域证明振荡的几何衰减。

命题3：对于$R>0$满足$\beta R^2 \leq T$（其中$\beta = \frac{9a_0}{32d}$），解$w:[−\beta R^2, 0]\times B(0,R)\to[0,1]$满足 $\text{osc}_{(-\beta R^2/16,0]\times B(0,R/4)} w \leq 1-\delta + C_f R^\gamma \|f\|_{L^p((−\beta R^2,0);L^q(B(0,R)))}$

证明策略：

• 分两种情况：$|\{w(-T_R,\cdot)\geq\frac{1}{2}\}\cap B(0,R/4)| \geq \frac{1}{2}|B(0,R/4)|$或其补集
• 构造比较函数$v$满足$(a_0c_0\partial_t - \Delta)v = f$
• 利用比较原理和引理2的下界估计，得到$w(t,x) > \delta - C_2R^\gamma\|f\|$
• 通过迭代获得几何衰减率

3. 边界附近估计（引理3）

关键引理3：对于$z\in\Omega$，令$d_z := d(z,\partial\Omega)$为到边界的距离，则 $w(t,x) \leq C(\|w^{in}\|_{\text{Lip}(\Omega)} + \|f\|_{L^p([0,T],L^q(\Omega))}) d_x^{\tilde{\gamma}}$ 其中$\tilde{\gamma} := \min(\gamma, 1-2\varepsilon)$

证明要点：

• 情况1（小时间$t^{1/2}\leq d_x$）：利用热核的积分表示和引理1
• 情况2（大时间$t > d_x^2$）：使用Gaussian估计 $\Gamma_{\Omega,D}(t,x,y) \lesssim \frac{d_x d_y}{t^{\frac{d+2}{2}}}e^{-c\frac{|x-y|^2}{t}}$ 将积分分解为三部分精细估计

4. 小时间估计（引理4）

引理4： $|w(t,x) - w(0,x)| \leq C(\|w^{in}\|_{C^1(\Omega)} + \|f\|_{L^p([0,T],L^q(\Omega))}) t^{\frac{\min(1,\gamma)}{2}}$

证明方法：

• 利用比较原理：$v_{a_0c_0} \leq w \leq v_{a_0}$
• 将$v_{a_0}$分解为$v_{a_0,f} + \tilde{v}_{a_0}$
• 对$\tilde{v}_{a_0}$使用Schauder估计得到抛物Hölder连续性
• 对$v_{a_0,f}$使用引理1

定理1的完整证明架构

三种情况的统一处理：

1. 边界附近（$\min(d_x,d_y) \leq 2|x-y|$）：
   • 直接应用引理3
   • 得到$|w(t,x)-w(t',y)| \lesssim |x-y|^{\tilde{\gamma}}$
2. 小时间（$t'^{1/2}, t^{1/2} \leq |x-y|$）：
   • 应用引理4和初值的Lipschitz性
   • 得到$|w(t,x)-w(t',y)| \lesssim |x-y|^{\tilde{\gamma}}$
3. 远离边界且大时间（$|x-y| \leq \min(t'^{1/2}, t^{1/2}, \frac{\min(d_x,d_y)}{2})$）：
   • 取尺度$R := \min(t^{1/2}, \min(d_x,d_y))$
   • 如果$\min(d_x,d_y) \leq t^{1/2}$：使用引理3
   • 如果$\min(d_x,d_y) > t^{1/2}$：使用引理4
   • 应用命题3的振荡衰减，通过迭代得到： $o_k \lesssim R^{\tilde{\gamma}}(\|f\| + \|w^{in}\|)\Lambda^k$ 其中$\Lambda > \max(\frac{1}{4^\gamma}, 1-\delta)$
   • 选择$k_0$使得$R4^{-k_0} \leq \max(|x-y|, (\frac{|t-t'|}{\beta})^{1/2})$
   • 得到$|w(t,x)-w(t',y)| \lesssim \max(|x-y|, |t-t'|^{1/2})^{\alpha_\Lambda}$

应用系统详解

应用1：可逆化学系统（命题1）

系统模型： $\partial_t u_i - d_i\Delta u_i = (-1)^i(u_1u_3 - u_2u_4), \quad i=1,2,3,4$ 表示可逆反应$A_1 + A_3 \rightleftharpoons A_2 + A_4$

证明策略：

1. 标准先验估计（式12）： $\sum_{i=1}^4 \int_\Omega \frac{u_i^{p+1}}{p+1}(T) + \frac{4p}{(p+1)^2}\sum_{i=1}^4 d_i\int_0^T\int_\Omega |\nabla(u_i^{\frac{p+1}{2}})|^2 \leq \text{初值} + C_p\sum_{i=1}^4\int_0^T\int_\Omega u_i^{p+2}$
2. 定义辅助函数：$w := \int_0^t(\sum_{i=1}^4 d_i u_i)$
   • 验证：$w\geq 0$, $\partial_t w \geq 0$, $w|_{\partial\Omega}=0$
   • 计算：$\Delta w = \sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}$
   • 定义$a := \frac{\sum_{i=1}^4 u_i}{\sum_{i=1}^4 d_i u_i}$满足$\frac{1}{\max d_i} \leq a \leq \frac{1}{\min d_i}$
   • 得到：$a\partial_t w - \Delta w = \sum_{i=1}^4 u_i^{in}$
3. 应用定理1：得到$\sum_{i=1}^4 u_i \in C^{0,\alpha}([0,T]\times\Omega)$
4. 插值不等式（式13）： $\left\|\sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}\right\|_{L^{\frac{2(3-\alpha)}{2-\alpha}}(\Omega)}^3 \leq C\|w\|_{C^{0,\alpha}(\Omega)}^{\frac{3}{3-\alpha}} \left\|\nabla\left[\sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}\right]\right\|_{L^2(\Omega)}^{\frac{3(2-\alpha)}{3-\alpha}}$
5. Bootstrap论证：结合估计(12)和(13)，得到$u_i \in L^{3+\delta}([0,T]\times\Omega)$，进而得到所有$L^p$估计

应用2：SKT交叉扩散系统（命题2）

系统模型：