KPP-Gleichung

Die Kolmogorov-Petrovsky-Piscounov-Gleichung (KPP-Gleichung nach Andrei Kolmogorow, Iwan Petrowski und Nikolai Piskunow 1937) ist eine nichtlineare partielle Differentialgleichung von der Form einer Reaktions-Diffusions-Gleichung.

Ein Spezialfall ist Fishers-Gleichung (nach Ronald Aylmer Fisher 1937) der Populationsdynamik, eine stetige Variante der Logistischen Gleichung (siehe auch Logistische Funktion).

Hauptteil

Die KPP-Gleichung hat die Form:

$\partial _{t}u=\partial _{x}^{2}u+g(u)$

mit einer nichtlinearen Funktion $g(u)$, die erfüllt: $g(1)=g(0)=0$, $g(u')>0$ und ${\tfrac {dg}{du}}(u')<{\tfrac {dg}{du}}(0)$ für $0<u'<1$ (Das Intervall [0,1] ist häufig auch das Definitionsintervall der Variablen $u$, wenn diese eine Konzentration angibt).

Fishers Gleichung ist ein Spezialfall der Form:

$\partial _{t}u=\partial _{x}^{2}u+u\cdot (1-u)=\partial _{x}^{2}u+u-u^{2}$

Manchmal wird statt des Reaktionsterms $g(u)=u\cdot (1-u)$ auch ein Term $g(u)=r\cdot u\cdot (1-u)$ angegeben, ähnlich wie bei der Logistischen Gleichung.

Dies ist eine semilineare parabolische Gleichung zweiter Ordnung. Sie wird verwendet, um verschiedene Vorgänge in der Natur zu modellieren, beispielsweise die Populationsdynamik oder chemische Reaktionen.

Die Differentialgleichung besteht aus einem Diffusionsterm $\partial _{x}^{2}u$ und einem nichtlinearen Reaktionsterm $u-u^{2}$.

Verwendet man eine ortsunabhängige Funktion $f(t)=u(x,t)$, so erhält man die gewöhnliche Differentialgleichung

$\partial _{t}f=f-f^{2}$.

An dieser kann man erkennen, dass mit dem Modell ein exponentielles Wachstum $\partial _{t}f=f$ modelliert wird, das jedoch einen Sättigungsterm $-f^{2}$ enthält. Dieser steht z. B. bei der Populationsdynamik für die begrenzte Nahrungsversorgung oder bei chemischen Reaktionen für die Sättigung der Konzentration.

Reaktionsfronten

Verwendet man die Gleichung zur Modellierung einer örtlich lokalisiert startenden Reaktion, so ist klar, dass sich eine Reaktionsfront ausbildet. Diese besitzt, wie man zeigen kann, eine minimale Ausbreitungsgeschwindigkeit.

Verwendet man den für Wellen üblichen Ansatz

$u(x,t)=f(x-vt)=f(w)$,

so erhält man nach Einsetzen die gewöhnliche Differentialgleichung zweiter Ordnung

$\partial _{w}^{2}f+v\partial _{w}f+f-f^{2}=0$.

Nach Linearisierung und unter der Annahme, dass die "Konzentration" f nur Werte zwischen 0 und 1 annehmen kann, erhält man die Gleichung für die Eigenwerte

$\lambda _{1,2}={\frac {-v\pm {\sqrt {v^{2}-4}}}{2}}$.

Da diese für stabile Wellen reell sein müssen, muss $v\geq 2$ gelten.

Verallgemeinerungen

Die Fisher-Gleichung kann verallgemeinert werden zu:

$\partial _{t}u=\partial _{x}^{2}u+(1-u^{m})u$

mit einer positiven ganzen Zahl $m$.

Im Fall der Fisher-Gleichung gilt dann $m=1$.

Siehe auch

• Nernst-Planck-Gleichung