Fronty Reakcji: Jak Fala Chemiczna Wybiera Swoją Prędkość

Fala, która się wzmacnia

Wypuść iskru do mieszaniny reaktywów, które katalizują ich własne wytwarzanie, i coś bardziej uporządkowanego niż losowy rozkład nastąpi: powstaje ostry front, który porusza się z jednego, dobrze zdefiniowanego tempa, pozostawiając za sobą całkowicie zareagowane materiały i niezmienione reaktywy z przodu.

Zachowanie to – dyfuzja połączona z lokalnym, autocatalitycznym wzrostem, tworząc samorządną falę podróżującą – opisuje się jednym z najbardziej rozpowszechnionych równań w matematyce stosowanej, równanie Fishera-KPP, niezależnie zaproponowane w 1937 roku przez Ronald Gibsona (badającego przestrzenną rozprzestrzenianie się korzystnego genu w populacji) i Andrieja Kolmogorov, Iwana Petrowskiego i Nikołaja Piskunowa (badających tę samą matematykę bardziej rygorystycznie i ogólnie).

Równanie

∂u/∂t = D · ∇²u + r · u · (1 - u) u(x,t) = lokalna koncentracja "zareagowanego" materiału, znormalizowana do 0…1 D = współczynnik dyfuzji — jak szybko rozprzestrzenia się produkt reakcji przestrzennie r = wskaźnik wzrostu (reakcja) wewnętrzny D·∇²u = TERMIN DYFUZYJNY — rozprzestrzenia reakcję na zewnątrz od obszaru, gdzie jest ona wysoka r·u·(1-u) = TERMIN REAKCYJNY — wzrost logistyczny: szybki, gdy u jest małe, saturując się, gdy u→1 Termin reakcyjny jest dokładnie modelem wzrostu logistycznego: produkcja jest proporcjonalna do ilości obecnej już reakwenta (u, autocataliza) i do ilości pozostałego niezapewnionego materiału (1-u, pojemność nośna saturująca). Bez dyfuzji ten termin sam w sobie sprawiłby, że u rośnie płynnie od 0 do 1 wszędzie naraz. Dyfuzja to właśnie przekształca proces o charakterze jednorodnym przestrzennie w prawdziwą falę podróżną: pozwala niewielkiej ilości produktu tuż przed obszarem zareagowanym „posiać” kolejny obszar niezapewnionego materiału, który następnie rośnie autocatalitycznie i dalej się dyfunduje, rozprzestrzeniając front do przodu.

∂u/∂t = D · ∇²u  +  r · u · (1 - u)

u(x,t) = local concentration of "reacted" material, normalised to 0…1
D      = diffusion coefficient — how fast the reaction product spreads spatially
r      = intrinsic growth (reaction) rate
D·∇²u  = DIFFUSION term — spreads the reaction outward from wherever it's high
r·u·(1-u) = REACTION term — logistic growth: fast when u is small, saturating as u→1

Dlaczego prędkość wynosi dokładnie v = 2√(Dr)

Rozwiązanie Fisher-KPP dopuszcza ruchome fale reakcji w postaci u(x,t) = U(x − vt) dla całego zakresu prędkości v ≥ v_min, ale przy założeniu lokalnego punktu reakcji. Jednak front, który faktycznie powstaje – udowodniono to w oryginalnym artykule KPP – zawsze wybiera najmniejszą dopuszczalną prędkość:

v = 2·√(D·r) (prędkość wybrana przez każdą ostro zlokalizowaną warstwę początkową) Intuicja jest taka, że daleko przed frontem, u jest bardzo małe, a w tym obszarze wyrażenie reakcji r·u(1-u) upraszcza się do r·u – czysty wzrost wykładniczy, bez saturacji. Równanie różniczkowe liniowe z tego obszaru prowadzącej krawędzi daje ∂u/∂t = D·∇²u + r·u, którego najszybciej rosnąca, nieujemna i ograniczona rozwiązanie porusza się dokładnie z 2√(Dr). Nieliniowa saturacja za frontem nie zmienia tej prędkości prowadzącej krawędzi – wpływa jedynie na kształt profilu, który pojawia się za nią, "odnaczając" to, co już zostało zasiane przez liniową krawędź prowadzącą.

v = 2·√(D·r)          (speed selected by any sharply-localized initial condition)

Ta sama zasada, wiele przebiegłości

Ten wybór prędkości dyfuzji i wzrostu nie jest specyficzny dla chemii. Fisher wywodził go pierwotnie do opisu rozprzestrzeniania się korzystnych genów w populacji o zróżnicowanej przestrzennej dystrybucji – szybkość wzrostu populacji odgrywa rolę r, a losowe rozprowadzanie (migracja) odgrywa rolę D. W procesach spalania analogiczna prędkość płomienia mieszanek wstępnych zależy od współczynnika dyfuzji termicznej i szybkości reakcji paliwa w sposób podobny. Ekologowie wykorzystują tę samą równanie do opisu forsowania gatunków – gatunek wprowadzony o lokalnej szybkości wzrostu r i losowym rozprowadzaniu D posuwa się własną granicę zasięgu o wartości 2√(Dr), co jest przewidywaniem, które dobrze oddawało obserwowane fronty inwazyjne dla gatunków od muskraków w Europie po jolowe kłębaki w Australii jako pierwsza aproksymacja.

Symulowanie na siatce

Prosty schemat różnic skończonych z eksplizitnym krokiem Eulera odtwarza front czysto, o ile przestrzega się granicy stabilności dyfuzji:

Schemat opiera się na kroku Eulera (explicit Euler step) dla 1D modelu Fishera-KPP, z rozdzielczością siatki dx i krokiem czasowym dt. Następnie, dla każdego punktu w siatce (od indeksu 1 do N-2), obliczana jest pochodna drugiego rzędu (Laplacian) przy pomocy wzoru: `const laplacian = (u[i+1] - 2*u[i] + u[i-1]) / (dx * dx);` Wartość `uNext[i]` jest obliczana jako następny krok w czasie, uwzględniając dyfuzję (D) i ruch pod wpływem gradientu koncentracji (r): `uNext[i] = u[i] + dt * (D * laplacian + r * u[i] * (1 - u[i]));`

Stabilność symulacji wymaga, aby dt ≤ dx² / (2D) – standardowego limitu dyfuzji CFL. Pomiar pozycji frontu (np. gdzie u przekracza 0,5) w czasie i dopasowanie do niego prostej pozwala na sprawdzenie symulacji względem prognozy 2√(Dr), po początkowym okresie przejścia, gdy profil frontu ustala się w kształcie samopowtarzalnym. Dopasowanie powinno zbiegać się do analitycznej prędkości w granicach 1-2% dla siatki o wystarczającej rozdzielczości.

// explicit Euler step for 1D Fisher-KPP, grid spacing dx, step dt
for (let i = 1; i < N - 1; i++) {
  const laplacian = (u[i+1] - 2*u[i] + u[i-1]) / (dx * dx);
  uNext[i] = u[i] + dt * (D * laplacian + r * u[i] * (1 - u[i]));
}
// stability requires dt <= dx² / (2D)  — the standard diffusion CFL limit

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Reaction Front i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.