PL 
EnglishDomyślny
Українська🔥 Front reakcji — biegnąca fala chemiczna
Symuluj biegnący front falowy z równania Fishera-KPP. Patrz, jak reakcje autokatalityczne rozchodzą się z prędkością v = 2√(Dr).
Podobne symulacje
O symulacji frontu reakcji
Równanie Fishera-KPP - nazwane od nazwiska Ronalda Fishera (1937) oraz niezależnie od Kołmogorowa, Pietrowskiego i Piskunowa - to jedno z najbardziej fundamentalnych równań różniczkowych cząstkowych w biologii matematycznej i chemii. Łączy ono dwa procesy fizyczne: dyfuzję, która rozprzestrzenia substancję w przestrzeni, oraz wzrost logistyczny, który modeluje autokatalityczne wzmocnienie, nasycające się, gdy substrat się wyczerpuje. Efektem jest wędrująca fala, która zachowuje swój kształt w nieskończoność, propagując się ze stałą prędkością określoną wyłącznie przez dwa parametry D i r.
Ta symulacja rozwiązuje 1D równanie różniczkowe cząstkowe Fishera-KPP numerycznie za pomocą różnic skończonych z jawnym schematem czasowym Eulera. Domena przestrzenna ma N=600 komórek z warunkami brzegowymi Neumanna (zerowy strumień). Krok czasowy jest dobierany adaptacyjnie, by spełnić kryterium stabilności dt ≤ dx²/(2D), zapewniając ograniczoność rozwiązania numerycznego. Na klatkę animacji przypadają trzy podkroki, co poprawia rozdzielczość czasową przy zachowaniu renderowania 60 kl./s.
Główny obszar wyświetla bieżący profil stężenia u(x) jako wypełniony wykres powierzchniowy z nałożoną linią. Wstawiony w prawym górnym rogu kymograf pokazuje diagram czasoprzestrzenny: każdy poziomy wiersz to migawka u(x) uchwycona w kolejnych krokach czasowych, przewijająca się w dół. Ukośna smuga frontu fali na kymografie bezpośrednio koduje prędkość fali - bardziej stromy ukos oznacza szybszą propagację. Porównaj zmierzoną prędkość fali z przewidywaniem teoretycznym v = 2√(D·r), by zweryfikować wzór Fishera-KPP.
Najczęściej zadawane pytania
Czym jest równanie Fishera-KPP?
Równanie Fishera-KPP (Kołmogorowa-Pietrowskiego-Piskunowa) ∂u/∂t = D∂²u/∂x² + ru(1-u) opisuje rozprzestrzenianie się korzystnego genu (Fisher, 1937) lub autokatalitycznej reakcji chemicznej. Wyraz D∂²u/∂x² modeluje dyfuzję, natomiast ru(1-u) to człon wzrostu logistycznego, który nasyca się przy u=1.
Czym jest wędrujący front falowy?
Wędrujący front falowy to rozwiązanie u(x,t) = f(x - vt), które zachowuje swój kształt, poruszając się ze stałą prędkością v. W równaniu Fishera-KPP dowolny warunek początkowy o zwartym nośniku ewoluuje w falę poruszającą się z minimalną prędkością v = 2√(Dr), wybraną zgodnie z zasadą najstromszego zejścia (frontu ciągnionego).
Co oznacza wzór na prędkość fali v=2√(Dr)?
Wzór na prędkość fali v = 2√(Dr) pokazuje, że zarówno szybsza dyfuzja (większe D), jak i szybszy wzrost (większe r) zwiększają prędkość propagacji. Ten wzór wynika z linearyzacji równania wokół u=0: front jest „ciągnięty” przez swoją czołową krawędź, gdzie u jest małe, a wzrost jest w przybliżeniu wykładniczy.
Czym jest autokataliza?
Autokataliza to reakcja chemiczna, w której produkt katalizuje własne powstawanie. Jeśli istnieje niewielka ilość produktu, przyspiesza ona dalszą produkcję, prowadząc do wzrostu wykładniczego, aż substraty się wyczerpią - dokładnie takie zachowanie logistyczne ru(1-u) jest tu modelowane.
Czym różni się fala bistabilna od Fishera-KPP?
Fala bistabilna wykorzystuje sześcienny człon reakcji r·u·(u-a)·(1-u) z progiem a. W przeciwieństwie do Fishera-KPP, która rozprzestrzenia się z dowolnego dodatniego zarodka, fale bistabilne rozprzestrzeniają się tylko wtedy, gdy początkowe zaburzenie przekracza próg a - co ma znaczenie dla modeli impulsów nerwowych i przełączników chemicznych.
Czym jest analogia do impulsu nerwowego?
Impuls nerwowy Hodgkina-Huxleya to wędrująca fala w układzie reakcji-dyfuzji, w którym napięcie błonowe u dyfunduje wzdłuż aksonu, a bramki jonowe zapewniają nieliniowość bistabilną. Model Fishera-KPP jest uproszczoną wersją oddającą istotną fizykę propagacji fali.
Czym jest zasada najstromszego zejścia przy wyborze prędkości fali?
Gdy warunek początkowy zanika szybciej niż e^(-λx) dla λ = √(r/D), fala Fishera-KPP ustala się na minimalnej prędkości v_min = 2√(Dr). Bardziej strome warunki początkowe (szybszy zanik) dają szybsze fale; prędkość 2√(Dr) to najwolniejsza możliwa stabilna fala, wybierana przy t→∞.
Czym jest diagram czasoprzestrzenny (kymograf)?
Wstawka pokazuje kymograf: obraz 2D, w którym oś pozioma to przestrzeń x, a oś pionowa to czas t (przewijający się w dół). Każdy poziomy pasek to migawka u(x) w jednym momencie. Ukośna smuga to wędrujący front falowy, którego nachylenie bezpośrednio daje prędkość fali.
Dlaczego front reakcji zachowuje swój kształt?
Kształt jest zachowywany, ponieważ dyfuzja ma tendencję do rozszerzania frontu (czyniąc go szerszym), podczas gdy człon reakcji przeciwdziała temu poszerzaniu (szybszy wzrost na froncie go stroma). Te dwa efekty równoważą się przy charakterystycznej szerokości ~√(D/r), nadając fali stabilny profil.
Gdzie w naturze pojawiają się fale reakcji-dyfuzji?
Wędrujące fale chemiczne pojawiają się w reakcji Biełousowa-Żabotyńskiego (spiralne fale w naczyniach), sygnalizacji wapniowej w komórkach (fale zapłodnienia), rozprzestrzenianiu się gatunków inwazyjnych (pierwotna motywacja Fishera), potencjałach czynnościowych w neuronach oraz rozprzestrzenianiu się epidemii w populacjach.