Hölder-Regularität parabolischer Gleichungen mit Dirichlet-Randbedingungen und Anwendung auf Reaktions-Diffusions- und Reaktions-Kreuziffusionssysteme

Grundinformationen

• Paper-ID: 2511.18872
• Titel: Hölder regularity of parabolic equations with Dirichlet boundary conditions and application to reaction-diffusion and reaction-cross-diffusion systems
• Autoren: Hector Bouton, Laurent Desvillettes, Helge Dietert (Université Paris Cité and Sorbonne Université, CNRS, IMJ-PRG)
• Klassifikation: math.AP (Analyse partieller Differentialgleichungen)
• Einreichungsdatum: 24. November 2025
• Paper-Link: https://arxiv.org/abs/2511.18872v1

Zusammenfassung

Der vorliegende Artikel verallgemeinert die jüngsten Arbeiten der Autoren zu Neumann-Randbedingungen auf den Fall von Dirichlet-Randbedingungen. Das Kernforschungsobjekt ist die parabolische Gleichung mit rauen Koeffizienten $a$ der Form $a\partial_t w - \Delta w = f$, versehen mit homogenen Dirichlet-Randbedingungen und der speziellen Annahme $\partial_t w \geq 0$. Basierend auf diesem Ergebnis wird die Existenz globaler starker Lösungen für das dreidimensionale Shigesada-Kawasaki-Teramoto (SKT)-Kreuziffusionssystem mit Lotka-Volterra-Reaktionstermen sowie Schätzungen für Lösungen reversibler chemischer Reaktions-Diffusions-Systeme unter Dirichlet-Randbedingungen nachgewiesen.

Forschungshintergrund und Motivation

Problemstellung

1. Kernproblem: Untersuchung der Hölder-Regularität parabolischer Gleichungen mit nichtkonstanten rauen Koeffizienten unter Dirichlet-Randbedingungen. Konkret wird die Gleichung betrachtet: $a(t,x)\partial_t w(t,x) - \Delta w(t,x) = f(t,x), \quad \partial_t w \geq 0$ wobei $a$ die Elliptizitätsbedingung $0 < a_0 \leq a \leq c_0 a_0 < \infty$ erfüllt.
2. Bedeutung:
   • Hölder-Regularität ist ein fundamentales Problem der Theorie partieller Differentialgleichungen und entscheidend für das Verständnis von Lösungseigenschaften
   • Das Ergebnis findet Anwendung in der Biomathematik, chemischen Reaktionsdynamik und anderen Bereichen nichtlinearer Systeme
   • Es liefert ein wichtiges technisches Werkzeug für die globale Existenz von Kreuziffusions- und Reaktions-Diffusions-Systemen
3. Beschränkungen bestehender Methoden:
   • Das klassische Ergebnis von Krylov-Safonov (1981) deckt zwar allgemeinere Fälle ab, aber die Konstanten sind nicht explizit und erfordern stärkere räumliche Annahmen
   • Die vorherige Arbeit der Autoren boutondesvillettes2025 behandelt nur Neumann-Randbedingungen
   • Bei Dirichlet-Randbedingungen ist das Verhalten der Lösung in der Nähe des Randes komplexer und erfordert neue technische Ansätze
4. Forschungsmotivation:
   • Nutzung der Monotonieannahme $\partial_t w \geq 0$ zur Vereinfachung des Beweises und Gewinnung expliziter Konstanten
   • Nur Kontrolle von $f$ in kritischen Lebesgue-Räumen erforderlich
   • Vollständige Behandlung der Randbedingungen mit expliziter theoretischer Unterstützung für Anwendungen

Kernbeiträge

1. Haupttheoretisches Ergebnis (Satz 1): Unter der Bedingung $\gamma := 2 - \frac{2}{p} - \frac{d}{q} > 0$ wird nachgewiesen, dass die nichtnegative Lösung $w$ der Gleichung (1) zu $C^{0,\alpha}([0,T]\times\Omega)$ gehört mit der Schätzung: $\|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. Anwendung auf reversible chemische Systeme (Proposition 1): Nachweis der globalen Existenz starker Lösungen für das reversible chemische Reaktionssystem $\partial_t u_i - d_i\Delta u_i = (-1)^i(u_1u_3 - u_2u_4), \quad i=1,\ldots,4$ unter Dirichlet-Randbedingungen in vier Dimensionen und darunter
3. Anwendung auf SKT-Kreuziffusionssystem (Proposition 2): Nachweis der globalen Existenz starker Lösungen für das dreieckige SKT-System\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}$$ unter Dirichlet-Randbedingungen in vier Dimensionen und darunter
4. Technische Innovationen: Entwicklung spezialisierter Oszillationsdämpfungstechniken für Dirichlet-Randbedingungen, die sich wesentlich vom Neumann-Fall unterscheiden

Methodische Darstellung

Aufgabendefinition

Eingaben:

• Beschränktes $C^2$-Gebiet $\Omega \subset \mathbb{R}^d$
• Koeffizient $a(t,x)$ erfüllend die Elliptizitätsbedingung
• Quellterm $f \in L^p((0,T];L^q(\Omega))$, wobei $\gamma = 2-\frac{2}{p}-\frac{d}{q} > 0$
• Anfangswert $w^{in} \in C^1(\Omega)$ mit $w^{in}|_{\partial\Omega} = 0$

Ausgaben:

• Hölder-Stetigkeit der Lösung: $w \in C^{0,\alpha}([0,T]\times\Omega)$
• Explizite Schätzung der Hölder-Norm

Nebenbedingungen:

• Homogene Dirichlet-Randbedingung: $w = 0$ auf $(0,T]\times\partial\Omega$
• Monotonie: $\partial_t w \geq 0$
• Nichtnegativität: $w \geq 0$

Technischer Kernrahmen

1. Dirichlet-Wärmekernschätzungen (Lemmata 1-2)

Lemma 2 (Kritische Unterschranken-Schätzung): Für den Dirichlet-Wärmekern $\Gamma_{B(0,R),D}$ in der Kugel $B(0,R)$, setze $T_R := a_0\frac{9R^2}{32d}$, dann gilt $\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}$

Beweisidee:

• Konstruktion einer Vergleichsfunktion $\Psi(t,x,y) := p(t,|x-y|) - \sup_{0\leq s\leq t}p(s,3R/4)$
• Anwendung des Maximumprinzips zur Gewinnung von $\Gamma_{B(0,R),D}(t,x,y) > \Psi(t,x,y)$
• Für $x,y\in B(0,R/4)$ und angemessene Zeitbereiche: präzise Berechnung einer positiven Unterschranke

2. Oszillationsdämpfung (Proposition 3)

Kernidee: Nachweis geometrischen Zerfalls der Oszillation in Bereichen weit entfernt vom Rand.

Proposition 3: Für $R>0$ mit $\beta R^2 \leq T$ (wobei $\beta = \frac{9a_0}{32d}$) erfüllt die Lösung $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)))}$

Beweisstrategien:

• Fallunterscheidung: $|\{w(-T_R,\cdot)\geq\frac{1}{2}\}\cap B(0,R/4)| \geq \frac{1}{2}|B(0,R/4)|$ oder Komplement
• Konstruktion einer Vergleichsfunktion $v$ erfüllend $(a_0c_0\partial_t - \Delta)v = f$
• Anwendung des Vergleichsprinzips und der Unterschranke aus Lemma 2, um $w(t,x) > \delta - C_2R^\gamma\|f\|$ zu erhalten
• Iteration zur Gewinnung geometrischer Zerfallsraten

3. Schätzungen in Randnähe (Lemma 3)

Schlüssellemma 3: Für $z\in\Omega$ sei $d_z := d(z,\partial\Omega)$ der Abstand zum Rand, dann gilt $w(t,x) \leq C(\|w^{in}\|_{\text{Lip}(\Omega)} + \|f\|_{L^p([0,T],L^q(\Omega))}) d_x^{\tilde{\gamma}}$ wobei $\tilde{\gamma} := \min(\gamma, 1-2\varepsilon)$

Beweishauptpunkte:

• Fall 1 (kleine Zeit $t^{1/2}\leq d_x$): Nutzung der Integraldarstellung des Wärmekerns und Lemma 1
• Fall 2 (große Zeit $t > d_x^2$): Verwendung von Gaussian-Schätzungen $\Gamma_{\Omega,D}(t,x,y) \lesssim \frac{d_x d_y}{t^{\frac{d+2}{2}}}e^{-c\frac{|x-y|^2}{t}}$ Zerlegung des Integrals in drei Teile mit feinen Schätzungen

4. Schätzungen für kleine Zeit (Lemma 4)

Lemma 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}}$

Beweismethode:

• Anwendung des Vergleichsprinzips: $v_{a_0c_0} \leq w \leq v_{a_0}$
• Zerlegung von $v_{a_0}$ in $v_{a_0,f} + \tilde{v}_{a_0}$
• Für $\tilde{v}_{a_0}$ Verwendung von Schauder-Schätzungen zur Gewinnung parabolischer Hölder-Stetigkeit
• Für $v_{a_0,f}$ Anwendung von Lemma 1

Vollständige Beweisstruktur von Satz 1

Einheitliche Behandlung von drei Fällen:

1. Randnähe ($\min(d_x,d_y) \leq 2|x-y|$):
   • Direkte Anwendung von Lemma 3
   • Ergebnis: $|w(t,x)-w(t',y)| \lesssim |x-y|^{\tilde{\gamma}}$
2. Kleine Zeit ($t'^{1/2}, t^{1/2} \leq |x-y|$):
   • Anwendung von Lemma 4 und Lipschitz-Eigenschaft des Anfangswertes
   • Ergebnis: $|w(t,x)-w(t',y)| \lesssim |x-y|^{\tilde{\gamma}}$
3. Fern vom Rand und große Zeit ($|x-y| \leq \min(t'^{1/2}, t^{1/2}, \frac{\min(d_x,d_y)}{2})$):
   • Wahl der Skala $R := \min(t^{1/2}, \min(d_x,d_y))$
   • Falls $\min(d_x,d_y) \leq t^{1/2}$: Verwendung von Lemma 3
   • Falls $\min(d_x,d_y) > t^{1/2}$: Verwendung von Lemma 4
   • Anwendung von Proposition 3 zur Oszillationsdämpfung, durch Iteration: $o_k \lesssim R^{\tilde{\gamma}}(\|f\| + \|w^{in}\|)\Lambda^k$ wobei $\Lambda > \max(\frac{1}{4^\gamma}, 1-\delta)$
   • Wahl von $k_0$ so dass $R4^{-k_0} \leq \max(|x-y|, (\frac{|t-t'|}{\beta})^{1/2})$
   • Ergebnis: $|w(t,x)-w(t',y)| \lesssim \max(|x-y|, |t-t'|^{1/2})^{\alpha_\Lambda}$

Detaillierte Darstellung der Anwendungssysteme

Anwendung 1: Reversible chemische Systeme (Proposition 1)

Systemmodell: $\partial_t u_i - d_i\Delta u_i = (-1)^i(u_1u_3 - u_2u_4), \quad i=1,2,3,4$ repräsentiert die reversible Reaktion $A_1 + A_3 \rightleftharpoons A_2 + A_4$

Beweisstrategien:

1. Standardmäßige A-priori-Schätzungen (Gleichung 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{Anfangswert} + C_p\sum_{i=1}^4\int_0^T\int_\Omega u_i^{p+2}$
2. Definition einer Hilfsfunktion: $w := \int_0^t(\sum_{i=1}^4 d_i u_i)$
   • Verifikation: $w\geq 0$, $\partial_t w \geq 0$, $w|_{\partial\Omega}=0$
   • Berechnung: $\Delta w = \sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}$
   • Definition $a := \frac{\sum_{i=1}^4 u_i}{\sum_{i=1}^4 d_i u_i}$ erfüllend $\frac{1}{\max d_i} \leq a \leq \frac{1}{\min d_i}$
   • Ergebnis: $a\partial_t w - \Delta w = \sum_{i=1}^4 u_i^{in}$
3. Anwendung von Satz 1: Gewinnung von $\sum_{i=1}^4 u_i \in C^{0,\alpha}([0,T]\times\Omega)$
4. Interpolationsungleichung (Gleichung 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-Argument: Kombination der Schätzungen (12) und (13) ergibt $u_i \in L^{3+\delta}([0,T]\times\Omega)$, woraus alle $L^p$-Schätzungen folgen

Anwendung 2: SKT-Kreuziffusionssystem (Proposition 2)

Systemmodell: