Diffusive Stabilität

Diffusive Stabilität ist eine Eigenschaft von Lösungen von Reaktionsdiffusionsgleichungen. Unter bestimmten Bedingungen haben diese, die für Diffusionen typische Eigenschaft, dass ihre $L^{\mathrm{\infty }}$-Norm für $t\to \mathrm{\infty }$ gegen 0 konvergiert, während die $L^1$-Norm gleichzeitig beschränkt bleibt. Insbesondere ist damit die konstante Nulllösung asymptotisch stabil. Anschaulich bedeutet dies, für eine Lösung, die die Konzentration eines Stoffes im Raum in Abhängigkeit von der Zeit modelliert, dass sich die Konzentration des Stoffes mit der Zeit gleichmäßig im Raum verteilt, während die gesamte Konzentration selbst allerdings immer beschränkt bleibt. Diese Eigenschaft von Lösungen wird auch als diffusive Stabilität bezeichnet.

Verhalten bei linearer Diffusionsgleichung

Für die lineare Diffusionsgleichung in einer Raumdimension

Lösung der linearen Diffusionsgleichung bei $t=1$ mit Anfangsbedingung $u_0(x)=e^{-x^2}$.

$\frac{\partial }{\partial t}u=\frac{{\partial }^2}{\partial x^2}u$ mit $x\in ℝ,t\in [0,\mathrm{\infty })$

mit Anfangsbedingung $u(x,0)=u_0(x)$ ist die Lösung gegeben durch Faltung der Anfangsbedingung mit der Fundamentallösung

$H(x,t)=\frac{1}{\sqrt{4\pi t}}{e}^{-\frac{x^2}{4t}}$.

Sei nun $u_0(\cdot )\in L^1$.

Dann lässt sich die Lösung allgemein schreiben als

$u(x,t)=(H(\cdot ,t)*u_0(\cdot ))(x)=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}H(x-y,t)u_0(y)dy$.

Nach der Young-Ungleichung für Faltungen folgt somit einerseits

$\Vert u(\cdot ,t){\Vert }_{L^{\mathrm{\infty }}}={\Vert \underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}H(\cdot -y,t)u_0(y)dy\Vert }_{L^{\mathrm{\infty }}}\le \Vert H(\cdot ,t){\Vert }_{L^{\mathrm{\infty }}}\Vert u_0(\cdot ){\Vert }_{L^1}$

Lösung der linearen Diffusionsgleichung bei $t=5$ mit Anfangsbedingung $u_0(x)=e^{-x^2}$.

andererseits auch

$\Vert u(\cdot ,t){\Vert }_{L^1}={\Vert \underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}H(\cdot -y,t)u_0(y)dy\Vert }_{L^1}\le \Vert H(\cdot ,t){\Vert }_{L^1}\Vert u_0(\cdot ){\Vert }_{L^1}$.

Außerdem können die Normen der Fundamentallösung explizit berechnet werden (siehe auch Normalverteilung und Fehlerintegral).

Es gilt

$\Vert H(\cdot ,t){\Vert }_{L^{\mathrm{\infty }}}=\frac{1}{\sqrt{2t\pi }}$

und

$\Vert H(\cdot ,t){\Vert }_{L^1}=\frac{1}{\sqrt{4\pi t}}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{-\frac{x^2}{4t}}=1$.

Somit gelten insgesamt die Normabschätzungen

Lösung der linearen Diffusionsgleichung bei $t=10$ mit Anfangsbedingung $u_0(x)=e^{-x^2}$.

$\Vert u(\cdot ,t){\Vert }_{L^{\mathrm{\infty }}}\le \frac{1}{\sqrt{4\pi t}}\Vert u_0(\cdot ){\Vert }_{L^1}\stackrel{t\to \mathrm{\infty }}{\to }0$

$\Vert u(\cdot ,t){\Vert }_{L^1}\le \Vert u_0(\cdot ){\Vert }_{L^1}$.

Diese Eigenschaften entsprechen genau dem, was man von einer Diffusion erwarten würde, nämlich, dass sich die Stoffkonzentration $u(x,t)$ für wachsendes $t$ immer weiter im Raum verteilt, ohne dass sich dabei die gesamte Stoffmenge erhöht.

Verallgemeinerung auf Reaktionsdiffusionsgleichungen

Unter bestimmten Bedingungen kann die Eigenschaft der diffusiven Stabilität auch auf allgemeine Reaktionsdiffusionsgleichungen der Form

$\frac{\partial }{\partial t}u=\frac{{\partial }^2}{\partial x^2}u+f(u)$ mit $u(0,x)=u_0(x)$

verallgemeinert werden.

Man betrachte beispielsweise die Gleichung

$\frac{\partial }{\partial t}u=\frac{{\partial }^2}{\partial x^2}u+u^p$

für $p\ge 1$. Die Lösung dieser Gleichung lässt sich schreiben als

$u(x,t)=(H(\cdot ,t)*u_0(\cdot ))(x)+\underset{0}{\overset{t}{\int }}(H(\cdot ,t-s)*u(\cdot ,s))(x)ds$.

Für Anfangsbedingungen der Form $u_0\in L^1\cap L^{\mathrm{\infty }}$ und $p>3$ lässt sich damit zeigen, dass für alle $C\ge 0$ ein $\delta >0$ existiert, sodass

aus

$\Vert u_0(\cdot ){\Vert }_{L^1}+\Vert u_0(\cdot ){\Vert }_{L^{\mathrm{\infty }}}\le \delta$

folgt, dass

$\Vert u(\cdot ,t){\Vert }_{L^{\mathrm{\infty }}}\le C(1+t{)}^{-\frac{1}{2}}$

$\Vert u(\cdot ,t){\Vert }_{L^1}\le C$.

Die Bedingung $p>3$ ist dabei wichtig, da für kleinere $p$ der Reaktionsterm $f(u)=u^p$ stärker als Diffusionsterm $\frac{\partial }{\partial x^2}$ins Gewicht fällt und damit Lösungen der Gleichung gegen unendlich konvergieren würden.

Für die Betrachtung allgemeiner Reaktionsdiffusionsgleichungen spielt daher der Reaktionsterm und dessen Verhalten für Lösungen für große Werte von $t$ im Vergleich zum Verhalten der Lösungen der linearen Diffusionsgleichung eine große Rolle.