Gradient higher integrability for degenerate parabolic double phase systems with two modulating coefficients

1. Problemstellung

Das Paper untersucht die Regularität schwacher Lösungen von entarteten parabolischen Double-Phase-Systemen, die zwei modulierende Koeffizienten beinhalten. Das System ist gegeben durch:
$$u_t - \text{div} \left( a(z)|Du|^{p-2}Du + b(z)|Du|^{q-2}Du \right) = -\text{div} \left( a(z)|F|^{p-2}F + b(z)|F|^{q-2}F \right)$$
in einem Raum-Zeit-Zylinder $\Omega_T = \Omega \times (0, T)$, wobei:

• $2 \le p \le q < \infty$.
• $a(z)$ und $b(z)$ nicht-negative modulierende Koeffizienten sind.
• Hauptannahmen:
  • $a(z)$ ist gleichmäßig stetig.
  • $b(z)$ ist Hölder-stetig ($C^{\alpha, \alpha/2}$).
  • Die Summe der Koeffizienten ist nach unten beschränkt: $a(z) + b(z) \ge \delta > 0$ (verhindert gleichzeitiges Verschwinden).
  • Die Exponenten erfüllen die Standard-Wachstumsbedingung für Double-Phase-Probleme: $q \le p + \frac{2\alpha}{n+2}$.

Das primäre Ziel ist die Etablierung einer lokalen höheren Integrierbarkeit des räumlichen Gradienten $Du$. Konkret streben die Autoren an zu beweisen, dass, wenn der Quellterm $F$ und der Gradient $Du$ eine Integrierbarkeit aufweisen, die mit der Energiedichte $H(z, \xi) = a(z)|\xi|^p + b(z)|\xi|^q$ zusammenhängt, $Du$ tatsächlich eine höhere Integrierbarkeit besitzt (d. h. $H(z, |Du|) \in L^{1+\varepsilon}_{loc}$).

2. Methodik

Die Autoren verwenden einen ausgefeilten Ansatz der intrinsischen Geometrie, kombiniert mit einer Phasentrennungsstrategie. Dies stellt eine signifikante Abkehr von früheren Arbeiten zu Systemen mit nur einem modulierenden Koeffizienten dar.

A. Intrinsische Geometrie

Das System zeigt unterschiedliches Wachstumsverhalten, abhängig von der relativen Größe der Koeffizienten $a(z)$ und $b(z)$. Um dies zu handhaben, definieren die Autoren drei Arten von intrinsischen Zylindern, skaliert mit einem Parameter $\lambda \ge 1$:

1. $p$-intrinsische Zylinder ($Q_{\rho, \lambda}$): Wo der $p$-Wachstumsterm dominiert ($a(z)\lambda^p \gtrsim b(z)\lambda^q$). Die Zeitskalierung ist $\sim \lambda^{2-p}$.
2. $q$-intrinsische Zylinder ($J_{\rho, \lambda}$): Wo der $q$-Wachstumsterm dominiert ($b(z)\lambda^q \gtrsim a(z)\lambda^p$). Die Zeitskalierung ist $\sim \lambda^{2-q}$.
3. $(p, q)$-intrinsische Zylinder ($G_{\rho, \lambda}$): Wo beide Terme vergleichbar sind. Die Zeitskalierung wird durch die Summe $\Lambda = a(z)\lambda^p + b(z)\lambda^q$ bestimmt.

B. Phasentrennung und Vergleichsschema

Die zentrale Neuerung liegt im delikaten Vergleichsschema, das zur Trennung der Phasen verwendet wird. Im Gegensatz zu früheren elliptischen oder parabolischen Studien mit nur einem Koeffizienten erzeugen das Vorhandensein zweier Koeffizienten und der parabolische Charakter neun verschiedene Fälle, basierend auf der Größe von $a(z)$ und $b(z)$ innerhalb der intrinsischen Zylinder:

• $p$-Phase: $a(z)$ ist groß, $b(z)$ ist klein.
• $q$-Phase: $b(z)$ ist groß, $a(z)$ ist klein.
• $(p, q)$-Phase: Beide Koeffizienten sind signifikant.

Die Autoren müssen den Übergang zwischen diesen Phasen sorgfältig handhaben und sicherstellen, dass die intrinsischen Skalierungsparameter ($\lambda$) und die Radien der Zylinder konsistent gewählt werden, um die notwendigen umgekehrten Hölder-Ungleichungen zu erfüllen.

C. Wichtige technische Schritte

1. Energieabschätzungen: Herleitung von Caccioppoli-artigen Ungleichungen, die an die Double-Phase-Struktur angepasst sind.
2. Umgekehrte Hölder-Ungleichungen: Beweis dieser Ungleichungen separat für jede der drei Phasen ($p$, $q$ und gemischt). Dies beinhaltet komplexe Abschätzungen unter Verwendung der Stetigkeit von $a(z)$ und der Hölder-Stetigkeit von $b(z)$, um die „Kreuzterme" zu kontrollieren, bei denen ein Koeffizient klein, aber nicht null ist.
3. Stop-Zeit-Argument: Ein verfeinertes Argument wird verwendet, um den geeigneten intrinsischen Radius für jeden Punkt in der Superlevel-Menge des Gradienten auszuwählen. Dies stellt sicher, dass der Punkt in eine der drei definierten Phasen fällt.
4. Vitali-Überdeckungslemma: Ein verallgemeinertes Überdeckungslemma wird konstruiert, um die neun möglichen Kombinationen intrinsischer Zylinder zu handhaben (z. B. ein $p$-Zylinder, der einen $q$-Zylinder überdeckt, usw.). Die Autoren beweisen, dass die Skalierungsfaktoren $\lambda$ über sich schneidenden Zylindern vergleichbar sind, was ein einheitliches Überdeckungsargument ermöglicht.
5. Gehring-Lemma: Der letzte Schritt kombiniert die umgekehrten Hölder-Ungleichungen mit dem Überdeckungslemma und dem Satz von Fubini, um die Integrierbarkeit von $L^1$ auf $L^{1+\varepsilon}$ zu steigern.

3. Hauptbeiträge

• Erstes parabolisches Ergebnis für zwei Koeffizienten: Dies ist das erste Regularitätsergebnis für parabolische Double-Phase-Systeme mit zwei modulierenden Koeffizienten. Frühere Ergebnisse beschränkten sich auf den elliptischen Fall oder parabolische Fälle mit nur einem Koeffizienten (z. B. $a(z) \equiv 1$).
• Einheitliche intrinsische Geometrie: Das Paper erweitert die Methode der intrinsischen Geometrie erfolgreich auf ein System, bei dem die dominante Phase zwischen $p$ und $q$ wechseln kann (oder gemischt ist), abhängig von den lokalen Werten zweier verschiedener Koeffizienten.
• Handhabung der gleichmäßigen Stetigkeit: Die Autoren lockern die Regularitätsanforderung an $a(z)$ auf gleichmäßige Stetigkeit (anstatt Hölder-Stetigkeit), was im Kontext von Double-Phase-Problemen eine allgemeinere und technisch anspruchsvollere Annahme darstellt.
• Neun-Fälle-Analyse: Die rigorose Behandlung der neun verschiedenen geometrischen Fälle, die aus der Interaktion der beiden Koeffizienten und der drei intrinsischen Phasen entstehen, stellt einen signifikanten technischen Fortschritt gegenüber der bestehenden Literatur dar.

4. Hauptergebnisse

Der Hauptsatz (Satz 1.1) etabliert die lokale höhere Integrierbarkeit des Gradienten. Konkret existieren für eine schwache Lösung $u$ Konstanten $\varepsilon_0 > 0$ und $c \ge 1$, sodass für jeden Zylinder $Q_{2r}(z_0) \subset \Omega_T$:

$$\iint_{Q_r(z_0)} H_1(z, |Du|)^{1+\varepsilon} \, dz \le c \left( \iint_{Q_{2r}(z_0)} H_1(z, |Du|) \, dz \right)^{1 + \frac{\varepsilon q}{2}} + c \left( \iint_{Q_{2r}(z_0)} H_1(z, |F|)^{1+\varepsilon} \, dz \right)^{\frac{q}{2}}$$

wobei $H_1(z, \kappa) = a(z)\kappa^p + b(z)\kappa^q + 1$ und $\varepsilon \in (0, \varepsilon_0)$.

Dieses Ergebnis impliziert, dass der Gradient $Du$ zu einem höheren Lebesgue-Raum gehört als ursprünglich angenommen, was ein entscheidender Schritt zur Beweisung weiterer Regularitätseigenschaften wie Hölder-Stetigkeit oder Calderón-Zygmund-Abschätzungen ist.

5. Bedeutung

• Modellierung anisotroper Materialien: Double-Phase-Systeme modellieren Materialien mit stark anisotropen Eigenschaften (z. B. Verbundwerkstoffe mit unterschiedlichen Verfestigungsverhalten). Die Einbeziehung zweier modulierender Koeffizienten bietet einen realistischeren Rahmen für solche Materialien als Modelle mit nur einem Koeffizienten.
• Theoretische Grundlage: Durch die Etablierung einer höheren Integrierbarkeit legt diese Arbeit das notwendige Fundament für zukünftige Forschung zur partiellen Regularität (Hölder-Stetigkeit) und zu Calderón-Zygmund-Abschätzungen für parabolische Double-Phase-Systeme mit zwei Koeffizienten.
• Methodischer Fortschritt: Die entwickelten Techniken, insbesondere die Handhabung der „gemischten" Phase und des Neun-Fälle-Vitali-Überdeckungsarguments, bieten ein robustes Werkzeugset zur Analyse anderer entarteter parabolischer Systeme mit nicht-standardisierten Wachstumsbedingungen.

Zusammenfassend löst dieses Paper ein großes offenes Problem in der Regularitätstheorie parabolischer partieller Differentialgleichungen, indem es Methoden der intrinsischen Skalierung erfolgreich an die komplexe, mehrphasige Umgebung anpasst, die durch zwei modulierende Koeffizienten erzeugt wird.