Fully Discrete Pointwise Smoothing Error Estimates for Measure Valued Initial Data

Technische Zusammenfassung: Voll diskrete punktweise Glättungsfehlerabschätzungen für maßwertige Anfangsdaten

Problemstellung
Dieser Artikel analysiert die numerische Approximation des homogenen parabolischen Problems:
$$\partial_t v - \Delta v = 0 \quad \text{in } I \times \Omega, \quad v = 0 \quad \text{auf } I \times \partial\Omega, \quad v(0) = v_0 \quad \text{in } \Omega,$$
wobei $\Omega \subset \mathbb{R}^N$ ($N=2,3$) ein beschränktes, konvexes, polygonales/polyedrisches Gebiet ist. Der Schwerpunkt liegt auf dem Fall, bei dem die Anfangsdaten $v_0$ ein reguläres Borel-Maß ($v_0 \in \mathcal{M}(\Omega)$) sind, wie beispielsweise eine Linearkombination von Dirac-Delta-Funktionen, und nicht eine Standard-$L^2$-Funktion. Die Autoren untersuchen die voll diskrete Approximation unter Verwendung kontinuierlicher finiter Elemente der Ordnung $s \in \{1, 2\}$ im Raum und Diskontinuierlicher-Galerkin-(dG)-Methoden beliebigen Grades $r \geq 0$ in der Zeit.

Das zentrale Ziel besteht darin, punktweise ($L^\infty$) Fehlerabschätzungen zum Endzeitpunkt $T$ für den Fehler $(v - v_{kh})(T)$ zu etablieren, speziell in einem Teilgebiet $\Omega_0 \subset \Omega$, in dem der Träger des Anfangsmaßes enthalten ist. Die gewünschte Abschätzung hat die Form:
$$\|(v - v_{kh})(T)\|_{L^\infty(\Omega_0)} \leq C(\Omega_0, T) \left( k^{2r+1} + \ell_{kh} h^{s+1} \right) \|v_0\|_{\mathcal{M}(\Omega)},$$
wobei $k$ der maximale Zeitschritt, $h$ die räumliche Gittergröße und $\ell_{kh}$ ein logarithmischer Term ist.

Methodik
Die Analyse verläuft durch mehrere distincte Stufen:

1. Schwache Formulierungen und Regularität: Die Autoren definieren „sehr schwache Lösungen" für das kontinuierliche Problem mit maßwertigen Anfangsdaten, wobei Testfunktionen in Räumen verwendet werden, die sicherstellen, dass die Dualitätspaarung mit Maßen wohldefiniert ist. Sie etablieren Regularitätsergebnisse, die zeigen, dass Lösungen der Wärmeleitungsgleichung mit Maßdaten für $t > 0$ sofort $L^2$-regulär werden, mit spezifischen zeitabhängigen Glättungsabschätzungen der Form $\|v(t)\|_{L^2} \leq C t^{-N/4} \|v_0\|_{\mathcal{M}}$.

2. Diskretisierungsschemata:

   • Zeit: Das Problem wird unter Verwendung eines dG($r$)-Schemas auf potenziell nicht-uniformen Zeiteinteilungen diskretisiert, die spezifische Gitterbedingungen erfüllen (quasi-Uniformität und beschränkte Variation). Für Maßdaten werden die Trial- und Testräume zu $W^{1,p}_0$ und $W^{1,p'}_0$ modifiziert, um Kontinuitätsanforderungen für die Anfangsauswertung sicherzustellen.
   • Raum: Es werden kontinuierliche Lagrange-finite Elemente vom Grad $s$ verwendet. Die Autoren definieren den diskreten Laplace-Operator $\Delta_h$ und die Ritz-Projektion $R_h$.

3. Glättungsabschätzungen: Ein kritischer Bestandteil der Methodik ist die Herleitung von Glättungsabschätzungen sowohl für die kontinuierlichen als auch für die diskreten (halb-diskreten und voll-diskreten) Lösungen. Die Autoren beweisen, dass diskrete Lösungen die Glättungseigenschaften des kontinuierlichen Problems erben, insbesondere dass $\|v_{kh}(t)\|_{L^\infty}$ relativ zur Norm der Anfangsdaten wie $t^{-N/4}$ abfällt. Diese Abschätzungen werden von $L^2$-Anfangsdaten auf Maßdaten unter Verwendung von Dualitätsargumenten erweitert.

4. Fehlerzerlegung: Der Gesamtfehler wird in drei Komponenten aufgeteilt:
   $$(v - v_{kh})(T) = (v - v_k)(T) + (v_k - R_h v_k)(T) + (R_h v_k - v_{kh})(T).$$

   • Der erste Term stellt den Zeitdiskretisierungsfehler dar.
   • Der zweite Term ist der räumliche Approximationsfehler der halb-diskreten Lösung (Ritz-Projektionsfehler).
   • Der dritte Term ist der voll-diskrete Fehler (Unterschied zwischen der halb-diskreten und der voll-diskreten Lösung).

5. Umgang mit Maßdaten: Um $L^\infty$-Abschätzungen für Maßdaten zu erhalten, wenden die Autoren eine „Cut-off"-Technik an. Sie teilen das Zeitintervall bei $T/2$ auf und behandeln die Lösung im ersten Teil als Quelle für den zweiten Teil. Dies ermöglicht es ihnen, den Glättungseffekt zu nutzen, um von den rauen Maßdaten zu $L^2$-Daten überzugehen, bevor sie auf der zweiten Hälfte des Zeitintervalls Standard-Fehlerabschätzungen anwenden.

Hauptbeiträge und Ergebnisse

• Hauptergebnis für lineare Elemente ($s=1$): Der Artikel etabliert die voll diskrete punktweise Fehlerabschätzung (Gleichung 2) für stückweise lineare finite Elemente auf konvexen Gebieten. Entscheidend ist, dass dieses Ergebnis keine zusätzlichen Regularitätsannahmen über die auf konvexen Gebieten verfügbare Standardregularität hinaus erfordert. Die Fehlerschranke lautet:
  $$\|(v - v_{kh})(T)\|_{L^\infty(\Omega_0)} \leq C(T, \Omega_0) \left( k^{2r+1} + \ell_{kh} h^2 \right) \|v_0\|_{\mathcal{M}(\Omega)}.$$

• Erweiterung auf quadratische Elemente ($s=2$): Die Autoren erweitern die Analyse auf quadratische finite Elemente. Diese Erweiterung erfordert Annahme 1, die $H^3$-Regularität für die Lösung des elliptischen Problems mit $H^1_0$-Daten postuliert (eine Bedingung, die beispielsweise auf rechteckigen Gebieten erfüllt ist). Unter dieser Annahme beweisen sie:
  $$\|(v - v_{kh})(T)\|_{L^\infty(\Omega_0)} \leq C(T, \Omega_0) \left( k^{2r+1} + \ell_{kh} h^3 \right) \|v_0\|_{\mathcal{M}(\Omega)}.$$
  Dies verbessert frühere Ergebnisse der Autoren, die auf lineare Elemente beschränkt waren.

• Interior-Abschätzungen: Die Abschätzungen sind im Wesentlichen „interior", was bedeutet, dass der $L^\infty$-Fehler auf einem Teilgebiet $\Omega_0$ ausgewertet wird, das den Träger des Anfangsmaßes enthält, und nicht global auf $\Omega$. Diese Unterscheidung ist entscheidend für den Umgang mit den Singularitäten von Maßdaten.

• Allgemeine Zeiteinteilungen: Die Ergebnisse gelten für allgemeine Zeiteinteilungen, die milde Bedingungen erfüllen, und nicht nur für uniforme Gitter, wodurch die Anwendbarkeit früherer $L^2$- und $L^\infty$-Stabilitätsergebnisse erweitert wird.

Bedeutung und Behauptungen
Die Autoren stellen fest, dass ihr Hauptbeitrag die Etablierung voll diskreter Fehlerabschätzungen für Galerkin-Methoden mit maßwertigen Anfangsdaten auf potenziell stark variierenden Zeiteinteilungen und quasi-uniformen Gittern ist.

Der Artikel hebt hervor, dass diese Ergebnisse für Anwendungen erforderlich sind, die eine Identifikation von Anfangsdaten aus Beobachtungen zum Endzeitpunkt betreffen, insbesondere wenn Beobachtungen an einer endlichen Anzahl von Punkten durchgeführt werden (was zu einer maßwertigen Randbedingung für das adjungierte Problem führt). Frühere Arbeiten der Autoren [15] lieferten ähnliche Abschätzungen für $L^2$-Anfangsdaten, doch ist die Erweiterung auf Maßdaten (Gleichung 2) erforderlich, um die spezifische Struktur adjungierter Probleme zu behandeln, die in Szenarien mit sparsamer optimaler Steuerung und punktueller Beobachtung auftreten.

Der Artikel vermerkt bescheiden, dass, während der Fall stückweise linearer Elemente ($s=1$) keine zusätzlichen Regularitäten erfordert, die Konvergenz höherer Ordnung für quadratische Elemente ($s=2$) von spezifischen Gebietsgeometrien (wie Rechtecken) abhängt, die die notwendige $H^3$-Regularität garantieren. Die Arbeit schließt eine Lücke in der Literatur bezüglich der punktweisen Fehleranalyse für parabolische Probleme mit singulären Anfangsdaten im voll-diskreten Setting.