Strona główna Chemia Wzory reakcji-dyfuzji 3D

🎨 Wzory reakcji-dyfuzji 3D

Model reakcji-dyfuzji Graya-Scotta w trójwymiarowej siatce wokseli renderowanej w Three.js. Dwa morfogeny u i v tworzą plamy, paski i spirale w zależności od parametrów zasilania i wygaszania.

Chemia2DŚredni60 FPS
reaction-diffusion-3d ↗ Otwórz osobno
Interfejs samej symulacji jest w języku angielskim.

Jak to działa

Model Graya-Scotta opisuje dwie reagujące ze sobą substancje chemiczne, U i V, na dyskretnej siatce. U dyfunduje szybciej i jest zasilane z zewnętrznego rezerwuaru; V dyfunduje wolniej i powstaje w wyniku reakcji U + 2V → 3V. Rywalizacja między lokalną aktywacją (V katalizuje własne powstawanie) a globalną inhibicją (U jest zużywane szybciej, niż dyfunduje z powrotem) tworzy samoorganizujące się wzory przestrzenne.

W każdym kroku: obliczany jest laplasjan każdej komórki za pomocą ważonego szablonu 9-punktowego, a następnie U i V są aktualizowane jednocześnie metodą Eulera w przód. Na klatkę animacji przypada wiele kroków, co przyspiesza formowanie się wzoru.

∂u/∂t = Du∙∇²u - u∙v² + f∙(1-u)
∂v/∂t = Dv∙∇²v + u∙v² - (f+k)∙v
∇²u[i,j] = 0.2∙(u[i±1,j]+u[i,j±1]) + 0.05∙diag - u[i,j]

Kolor każdego piksela koduje stężenie substancji V: niskie V jest ciemnoniebieskie, wysokie V jest jasnobursztynowe lub białe, w zależności od wybranej mapy kolorów. Kliknięcie na obszarze symulacji zasiewa mały fragment o wysokim stężeniu V, wywołując nowy wzrost wzoru z tego punktu.

Najczęściej zadawane pytania

Czym jest model reakcji-dyfuzji Graya-Scotta?

Model Graya-Scotta opisuje dwie substancje chemiczne U i V, które reagują i dyfundują. U przekształca się w V, gdy się spotkają: U + 2V → 3V, a V zanika z szybkością k. U jest uzupełniane z szybkością zasilania f. W zależności od f i k układ tworzy plamki, paski, spirale lub wzory labiryntowe.

Czym są wzory Turinga?

Wzory Turinga to wzory przestrzenne wyłaniające się z układu reakcji-dyfuzji, w którym aktywator i inhibitor dyfundują z różną szybkością. Alan Turing zaproponował w 1952 roku, że ten mechanizm może wyjaśniać biologiczne powstawanie wzorów, takich jak umaszczenie zwierząt, wzory pasków u ryb czy rozwój narządów.

Co kontrolują szybkość zasilania f i szybkość zabijania k?

Szybkość zasilania f kontroluje, jak szybko substancja U jest uzupełniana z zewnętrznego rezerwuaru. Szybkość zabijania k kontroluje, jak szybko substancja V jest usuwana. Różne pary (f,k) tworzą dramatycznie różne wzory: niskie f/k dają plamki, pośrednie dają paski, wysokie dają robaki lub labirynty.

Czym są równania różniczkowe cząstkowe Graya-Scotta?

∂u/∂t = Du∙∇²u - u∙v² + f∙(1-u); ∂v/∂t = Dv∙∇²v + u∙v² - (f+k)∙v. Tutaj Du i Dv to współczynniki dyfuzji (Du > Dv), ∇² to laplasjan, f to szybkość zasilania, a k to szybkość zabijania.

Jak obliczany jest laplasjan na siatce dyskretnej?

Na siatce 2D stosowany jest szablon 9-punktowy z wagami przekątnymi 0,05 i wagami kardynalnymi 0,2: ∇²u ≈ 0,2∙(N+S+E+W) + 0,05∙(NE+NW+SE+SW) - u[i,j]. Daje to gładsze wyniki niż podstawowy szablon 4-sąsiedzki.

Dlaczego reakcja-dyfuzja tworzy wzory umaszczenia zwierząt?

Podczas rozwoju zarodkowego morfogeny dyfundują przez tkankę i reagują zgodnie z kinetyką aktywator-inhibitor. Zasada lokalnej aktywacji i dalekozasięgowej inhibicji tworzy regularne rozmieszczenie plamek lub pasków, które odzwierciedlają rzeczywiste umaszczenie zwierząt, jak plamki geparda czy paski zebry.

Czym jest morfogeneza?

Morfogeneza to biologiczny proces, w którym organizm rozwija swój kształt. Układy reakcji-dyfuzji to jeden z matematycznych modeli tego, jak proste gradienty chemiczne mogą tworzyć złożone, powtarzalne wzory przestrzenne bez żadnego odgórnego planu.

Jak współczynniki dyfuzji wpływają na wzór?

Kluczowy jest stosunek Du/Dv. Jeśli obie substancje dyfundują jednakowo, wzory nie mogą powstać. Aktywator (V) musi dyfundować wolniej niż inhibitor (U). Zwykle Du ≈ 2∙Dv. Większy stosunek daje drobniejsze wzory; mniejszy stosunek czyni wzory grubszymi.

Czy układy reakcji-dyfuzji mogą wytwarzać chaos?

Tak. W pobliżu punktów bifurkacji układy reakcji-dyfuzji mogą wykazywać chaos czasoprzestrzenny - wzory nieustannie się rozpadają i odtwarzają w nieprzewidywalny sposób. Jest to powiązane z chaosem chemicznym w reakcji Biełousowa-Żabotyńskiego.

Czym jest reakcja Biełousowa-Żabotyńskiego?

Reakcja Biełousowa-Żabotyńskiego (BZ) to prawdziwy oscylator chemiczny, w którym stężenia oscylują między stanami utlenionymi i zredukowanymi. W cienkiej warstwie cieczy wytwarza rozszerzające się fale spiralne - żywy fizyczny pokaz dynamiki reakcji-dyfuzji.

O tej symulacji

Ważony szablon 9-punktowy oblicza laplasjan dwóch substancji chemicznych U i V w każdej komórce siatki, pięć kroków Eulera jest wykonywanych na klatkę animacji, a wynikowe stężenie V jest malowane za pomocą jednej z trzech map kolorów. Ponieważ pięć wbudowanych gotowych ustawień zmienia tylko dwie liczby, szybkość zasilania f i szybkość zabijania k, dokładnie to samo równanie reakcji-dyfuzji może ustalić się w plamki, paski, spirale, robaki lub labirynty wyłącznie w zależności od tego, gdzie te dwa pokrętła znajdują się względem siebie.

🔬 Co pokazuje

Wzór reakcji-dyfuzji Graya-Scotta rosnący na zewnątrz od zasianej centralnej plamki, samoorganizujący się w plamki, paski lub labirynty, w miarę jak substancja aktywująca V i inhibitor U rywalizują o dyfuzję i reakcję na całej siatce.

🎮 Jak korzystać

Wybierz gotowe ustawienie lub ręcznie dostosuj szybkość zasilania f, szybkość zabijania k, Du i Dv, przełącz mapę kolorów między Heat/Cool/Plasma, kliknij na obszarze, by zasiać nowe stężenie V w dowolnym miejscu, lub naciśnij P/R, by wstrzymać/zresetować.

💡 Czy wiesz, że?

Alan Turing zaproponował dokładnie ten mechanizm aktywator-inhibitor w 1952 roku, by wyjaśnić, jak jednorodny zarodek może samoistnie rozwinąć niejednorodne wzory - ta sama matematyka stojąca za plamami lamparta i paskami zebry działa na żywo w twojej przeglądarce.

Często zadawane pytania

Dlaczego różne ustawienia dają plamki, paski lub labirynty?

applyPreset() zmienia tylko wartości szybkości zasilania f i szybkości zabijania k wprowadzane do tych samych równań Graya-Scotta, a przestrzeń parametrów Graya-Scotta jest znana z wielkiej wrażliwości - niewielkie przesunięcia w tej parze (f,k) przesuwają tę samą reakcję w zupełnie inne rodziny wzorów stanu ustalonego.

Co się dzieje, gdy kliknę na obszarze symulacji?

Kliknięcie ustawia u=0,5 i v=0,25 w niewielkim kołowym obszarze komórek wokół kursora, sztucznie zasiewając nowe miejsce reakcji, które następnie rośnie i wchodzi w interakcję z istniejącym wzorem.

Dlaczego symulacja wykonuje 5 kroków na klatkę animacji?

stepRD(5) jest wywoływane raz na wywołanie requestAnimationFrame, uruchamiając aktualizację reakcji-dyfuzji pięć razy przed każdym przerysowaniem, dzięki czemu wzór widocznie ewoluuje i stabilizuje się znacznie szybciej, niż pozwoliłby jeden krok na klatkę.

Dlaczego Du musi być większe niż Dv, aby wzory mogły powstać?

U musi dyfundować szybciej niż V, aby lokalna aktywacja przez V mogła wyprzedzić dalekozasięgową inhibicję przez U; gdyby oba współczynniki dyfuzji były równe, obie substancje zlałyby się w płaską, pozbawioną wzoru równowagę.

Co dokładnie oblicza szablon 9-punktowy w funkcji laplacian()?

Waży cztery kardynalne sąsiednie komórki wagą 0,2, a cztery przekątne wagą 0,05, dając gładsze, bardziej izotropowe przybliżenie ciągłego laplasjanu niż prosta średnia 4-sąsiedzka.

Podobne symulacje