Wave 109 — 20 Nowych Symulacji Fizyki i Nauki (940 Łącznie, −60 do 1000)

Każda symulacja na tej stronie zaczęła się od pytania: jak naprawdę wygląda to równanie? Długość Debye'a jest wykładana na każdym kursie fizyki plazmy jako λ_D = √(ε₀k_BT / nq²), formuła przewidująca, jak daleko elektrostatyczne zaburzenie wnika w plazmę, zanim zostanie ekranowane. Ale dopóki nie zobaczysz mapy cieplnej potencjału elektrycznego zapadającego się wykładniczo od ładunku testowego, ekranowanie pozostaje abstrakcyjne. Wave 109 zadaje to pytanie dwadzieścia razy w dwudziestu różnych dziedzinach fizyki i nauki.

Fizyka

Długość Debye'a — Mapa Cieplna Ekranowania λ_D

W plazmie swobodne elektrony przegrupowują się wokół dowolnego nadmiarowego ładunku, by zneutralizować go na charakterystycznej długości Debye'a. Symulacja umieszcza ładunek testowy w centrum 2D siatki plazmy i renderuje wynikowy potencjał elektrostatyczny jako mapę cieplną w fałszywych kolorach, pokazując wykładniczy zanik φ(r) = (q / 4πε0r) exp(−r / λD). Suwaki kontrolują temperaturę elektronów T i gęstość liczbową n, a mapa cieplna przerysowuje się w czasie rzeczywistym, gdy λ_D skaluje się z √(T/n). Przejście od rozproszonego, dalekosiężnego potencjału przy niskiej gęstości do ciasno zamkniętego gorącego punktu przy wysokiej gęstości jest natychmiast czytelne w kolorze. Towarzyszący wykres pokazuje profil radialny w skali logarytmicznej, potwierdzając wykładniczą obwiednię.

Defekty Topologiczne — Wiry BKT w Modelu 2D XY

Przejście Berezinskiego-Kosterlitza-Thoulessa to jedno z najbardziej niezwykłych przejść fazowych w mechanice statystycznej: zachodzi w dwuwymiarowych systemach, gdzie prawdziwy porządek dalekiego zasięgu jest zabroniony przez twierdzenie Mermina-Wagnera, a mimo to coś mierzalnie ostrego wciąż dzieje się w temperaturze krytycznej TBKT. Mechanizmem jest rozerwanie par defektów topologicznych — par wir-antywir w modelu 2D XY, gdzie każdy spin jest wektorem jednostkowym na płaszczyźnie. Poniżej TBKT wiry i antywiry pozostają związane w pary, a system ma quasi-dalekozasięgowy porządek z algebraicznie zanikającymi korelacjami. Powyżej niej swobodne wiry namnażają się, a korelacje zanikają wykładniczo.

Symulacja uruchamia przemiatanie Monte Carlo modelu 2D XY na kwadratowej siatce, renderując orientacje spinów jako strzałki kolorowane wg kąta. Wiry są wykrywane przez obliczenie liczby nawinięcia wokół każdej płytki i podświetlane czerwonymi (dodatnimi) i niebieskimi (ujemnymi) znacznikami. Przeciąganie suwaka temperatury przez TBKT daje widoczne namnożenie niezwiązanych par defektów.

Rezonans Stochastyczny — Krzywa SNR Bistabilnego Langevina

Rezonans stochastyczny to sprzeczne z intuicją zjawisko, w którym dodanie szumu do systemu może poprawić jego zdolność wykrywania słabego sygnału. Cząstka w podwójnej studni potencjału napędzana słabą siłą okresową poniżej progu nie może sama przekroczyć bariery — ale przy optymalnym poziomie szumu fluktuacje termiczne przepychają ją przez barierę w synchronizacji z wymuszeniem, dając ostry pik stosunku sygnału do szumu jako funkcji amplitudy szumu. Równanie Langevina rządzące dynamiką to:

dx/dt = −dV/dx + A cos(ωt) + η(t)
V(x) = −x²/2 + x⁴/4
⟨η(t)η(t′)⟩ = 2Dδ(t−t′)

Symulacja całkuje to równanie numerycznie dla dużego zespołu cząstek, przemiatając intensywność szumu D i obliczając wyjściowe SNR przy częstotliwości wymuszania za pomocą okienkowanej DFT. Wynikowa krzywa SNR w kształcie dzwonu jest rysowana w czasie rzeczywistym, z zaznaczonym optymalnym poziomem szumu. Efekt rezonansu stochastycznego jest istotny dla neuronauki sensorycznej, gdzie neurony wydają się wykorzystywać szum tła do wykrywania słabych bodźców na progu.

Jamming Granularny — Upakowanie Bidyspersyjne φ_J ≈ 0,84

Gdy gęsty materiał granularny, taki jak piasek, kulki szklane czy mokra kawa mielona, jest ściskany lub ścinany, przechodzi przejście jammingowe: nagle nabiera sztywności i opiera się deformacji jak ciało stałe. Przejście zachodzi przy krytycznej frakcji upakowania φ_J, która zależy od kształtu cząstek i rozkładu rozmiarów. Dla bidyspersyjnej mieszaniny dysków (dwa rozmiary cząstek w stosunku liczbowym 1:1, dobranym, by tłumić krystalizację), punkt jammingowy leży w pobliżu φ ≈ 0,84 w dwóch wymiarach.

Symulacja wykorzystuje silnik dynamiki molekularnej z odpychającymi kontaktami miękkiej kuli, by ściskać system bidyspersyjnych dysków ze stałym tempem, śledząc ciśnienie P i liczbę koordynacyjną z w miarę wzrostu frakcji upakowania. Przejście jammingowe jest widoczne jako ostry wzrost ciśnienia i skok z w stronę izostatycznej wartości Maxwella z = 2d = 4. Żywa wizualizacja upakowania koloruje cząstki wg liczby ich kontaktów, ujawniając powstawanie sztywnej, obejmującej łańcuchy sił sieci przy φ_J.

Dyfrakcja Bragga — nλ = 2d sinθ, Wzór Proszkowy FCC/BCC

Dyfrakcja rentgenowska to podstawowe narzędzie do wyznaczania struktury krystalicznej. Gdy monochromatyczna wiązka o długości fali λ uderza w sieć krystaliczną o rozstawie płaszczyzn d pod kątem ślizgowym θ, interferencja konstruktywna zachodzi, gdy spełnione jest prawo Bragga: nλ = 2d sinθ. Różne rodziny płaszczyzn (charakteryzowane wskaźnikami Millera hkl) dają piki pod różnymi kątami, a intensywność każdego piku zależy od czynnika struktury, kodującego pozycje atomów w komórce elementarnej.

Symulacja oblicza wzory dyfrakcji proszkowej dla struktur FCC, BCC i prostej sześciennej, sumując po wszystkich dozwolonych odbiciach aż do zdefiniowanego przez użytkownika maksymalnego 2θ. Wzór jest wyświetlany jako diagram kreskowy nałożony na symulowany wzór pierścieni Debye'a-Scherrera. Przełączanie między FCC a BCC natychmiast zmienia systematyczne nieobecności: FCC zabrania odbić hkl o mieszanej parzystości, podczas gdy BCC tłumi piki, gdzie h+k+l jest nieparzyste.

Sieć Sprężynowa — Kryterium Maxwella z_c ≈ 4

James Clerk Maxwell wykazał w 1864 roku, że sieć sprężyn jest mechanicznie sztywna wtedy i tylko wtedy, gdy średnia liczba koordynacyjna z — średnia liczba sprężyn na węzeł — spełnia z ≥ 2d, gdzie d to wymiar przestrzenny. W dwóch wymiarach krytyczny próg to zc ≈ 4. Poniżej tej wartości sieć ma miękkie mody: wielkoskalowe deformacje, które nie kosztują energii. Powyżej niej sieć jest sztywna, a każda deformacja wymaga energii sprężystej.

Symulacja generuje losowe sieci sprężynowe na 2D siatce z regulowaną łącznością i oblicza mody normalne za pomocą rzadkiej algebry liniowej. Miękkie mody są podświetlone na czerwono, a sztywne mody na niebiesko, a liczba każdego z nich jest śledzona wobec przewidywania Maxwella. Kliknięcie węzła i przeciągnięcie pokazuje, czy sieć deformuje się sztywno, czy przystosowuje przemieszczenie przez miękkie przegrupowania o zerowej energii. Przejście od miękkiej do sztywnej w miarę dodawania wiązań jeden po drugim jest żywą realizacją kryterium Maxwella.

Kwantowa

Emisja Kropki Kwantowej — Równanie Brusa, Widma CdSe

Półprzewodnikowe kropki kwantowe to nanoskalowe kryształy wystarczająco małe, by uwięzienie kwantowe przesuwało ich przerwę energetyczną, czyniąc kolor emisji bezpośrednią funkcją rozmiaru. Równanie Brusa daje zależną od rozmiaru przerwę energetyczną sferycznej kropki kwantowej o promieniu R:

E(R) = Eg,bulk + ħ²π² / (2μR²) − 1.8e² / (4πεε0R)

Dla kropek kwantowych CdSe przerwa energetyczna materiału objętościowego wynosi 1,74 eV (bliska podczerwień), a uwięzienie przesuwa emisję ku niebieskiemu w kierunku widzialnego, gdy promień maleje z 5 nm do 1,5 nm, przemiatając od czerwieni przez pomarańcz, żółć i zieleń. Symulacja renderuje widmo emisji jako pik gaussowski, którego środkowa długość fali jest wyznaczana przez równanie Brusa, a szerokość modeluje poszerzenie fononowe. Paleta kolorowych kółek reprezentuje kropki różnych rozmiarów, świecące w kolorze swojej emisji. Przeciąganie suwaka promienia daje natychmiastową zmianę koloru, która dokładnie odtwarza paletę tęczy widoczną na zdjęciach roztworów kropek kwantowych CdSe.

Siatka Spinów — Ferromagnetyk Heisenberga Metropolis MC, T_c ≈ 2,27J/k

Model Heisenberga to kwantowo-mechaniczne uogólnienie modelu Isinga: każde miejsce sieciowe niesie wektor spinu S zamiast skalara ±1, a hamiltonian to H = −J Σ⟨ij⟩ Si · Sj. Model Heisenberga 3D ma przejście ferromagnetyczne przy Tc ≈ 2,27J/kB dla prostej sieci sześciennej.

Symulacja uruchamia Monte Carlo Metropolisa na 3D siatce Heisenberga, renderując 2D przekrój jako pole kolorów (orientacja spinu mapowana na odcień). Magnetyzacja, ciepło właściwe CV i podatność magnetyczna χ są rysowane jako funkcje temperatury. Przemiatanie przez Tc pokazuje zanik parametru porządku, pik ciepła właściwego i dywergencję podatności — klasyczne sygnatury ciągłego przejścia fazowego.

Optyka

Dysk Airy'ego — Bessel J₁, Kryterium Rayleigha, 2D ImageData

Gdy światło z punktowego źródła przechodzi przez okrągłą aperturę, nie tworzy idealnego punktowego obrazu. Dyfrakcja rozprasza światło we wzór Airy'ego, którego profil intensywności to:

I(θ) = I0 [2J1(ka sinθ) / (ka sinθ)]²

Kryterium rozdzielczości Rayleigha stwierdza, że dwa punktowe źródła są dopiero rozdzielone, gdy centralne maksimum jednego pada na pierwsze zero drugiego, przy rozdzielczości kątowej θ = 1,22λ/D. Symulacja renderuje 2D wzór Airy'ego, ewaluując funkcję Bessela na siatce pikseli za pomocą API ImageData, osiągając czasy klatek poniżej milisekundy dla canvasa 512×512. Dwa regulowane źródła punktowe mogą być umieszczone obok siebie; ich połączony wzór intensywności jest wyświetlany obok linii kryterium Rayleigha.

Chemia

Front Reakcji — PDE Fishera-KPP, v = 2√(Dr)

Równanie Fishera-KPP sprzęga logistyczny lokalny wzrost z dyfuzją przestrzenną:

∂u/∂t = D ∂²u/∂x² + ru(1−u)

Dla ostrych warunków początkowych asymptotyczna prędkość frontu fali to v = 2√(Dr), niezależna od profilu początkowego. Zaimplementowane są zarówno monostabilny front Fishera-KPP, jak i bistabilny wariant fali biegnącej z regulowanymi przez użytkownika parametrami. Profil przestrzenny 1D jest wyświetlany jako poruszający się pasek intensywności, pozycja frontu w czasie jest rysowana, by zweryfikować formułę prędkości asymptotycznej, a portret płaszczyzny fazowej w układzie współporuszającym się pokazuje orbitę heterokliniczną łączącą dwa stany równowagi.

Termodynamika

Entropia Mieszania — ΔS_mix = −nRΣx_i ln x_i, Symulacja Cząstek

Gdy gazy doskonałe się mieszają, entropia wzrasta o ΔSmix = −nRΣxi ln xi, zawsze dodatnia, osiągająca maksimum dla równych proporcji. Symulacja umieszcza cząstki dwóch lub trzech gatunków w oddzielnych komorach i usuwa dzielącą przegrodę. Każda cząstka wykonuje losowy spacer na 2D siatce, a chwilowa entropia mieszania jest obliczana z empirycznego składu w każdej komórce. Zbliżanie się do maksymalnej entropii jest rysowane w czasie rzeczywistym, z nałożoną teoretyczną asymptotą. Użytkownicy mogą porównać mieszanie z różnymi współczynnikami dyfuzji, by zobaczyć, jak tempo transportu wpływa na skalę czasową, ale nie na termodynamiczny punkt końcowy.

Matematyka

Bańka Mydlana — Relaksacja Laplace'a Gaussa-Seidela, Powierzchnia Minimalna

Błony mydlane minimalizują swoją powierzchnię przy ograniczeniach brzegowych, spełniając równanie Laplace'a ∇²h = 0 dla funkcji wysokości. Schemat relaksacji Gaussa-Seidela iteruje:

hi,j ← (hi+1,j + hi−1,j + hi,j+1 + hi,j−1) / 4

aż do zbieżności. Symulacja pozwala użytkownikom rysować dowolne krzywe brzegowe na siatce i obserwować, jak powierzchnia minimalna relaksuje się na miejsce, renderowana jako mapa 2D kolorowana wg wysokości i opcjonalnie jako siatka 3D. Prędkość zbieżności z nadrelaksacją kolejną (SOR) jest porównywana ze standardowym Gaussem-Seidelem. Związek z elektrostatyką — gdzie to samo równanie rządzi potencjałem elektrycznym między przewodnikami — jest demonstrowany przez przełączenie interpretacji brzegu z wysokości na napięcie.

3D Automat Komórkowy — Three.js InstancedMesh, Siatka 20³, 5 Reguł

Symulacja uruchamia 3D automat komórkowy 20×20×20, używając pięciu odrębnych zestawów reguł zaczerpniętych z rozszerzonej notacji Narodziny/Przetrwanie/Sąsiedztwo. Żywa siatka jest renderowana za pomocą Three.js InstancedMesh dla 8000 kostek przy interaktywnych częstotliwościach klatek, ze stanem każdej komórki mapowanym na przezroczystość i kolor. Reguły obejmują analog 3D Gry w Życie, regułę wzrostu kryształów produkującą formy dendrytyczne, regułę chaosu wypełniającą objętość aperiodyczną turbulencją oraz dwie reguły generujące samopodobne struktury fraktalne. Wewnętrzna struktura jest ujawniana przez dostosowanie progu przezroczystości, który czyni komórki o niskiej aktywności przezroczystymi.

Algorytmy

Sieć Boole'a — Kauffman NK, Granica Fazy Krytycznej K = 2

Model sieci boolowskiej NK Stuarta Kauffmana umieszcza N węzłów binarnych, z których każdy otrzymuje wejście z dokładnie K losowo wybranych innych węzłów i aktualizuje się przez losową funkcję boolowską. Przy krytycznej łączności K = 2 sieć znajduje się na granicy między porządkiem a chaosem: atraktory skalują się jako √N, informacja propaguje się globalnie, a jednobitowe zaburzenia ani nie zanikają natychmiast, ani nie kaskadują wiecznie — właściwości, które Kauffman zaproponował jako sygnatury biologicznych sieci regulacji genów.

Symulacja implementuje sieci NK z regulowanym N (do 64) i K, wizualizując przestrzeń stanów jako sieć węzłów kolorowanych wg stanu binarnego z podświetlonymi cyklami atraktorów. Wyżarzone przybliżenie Derridy dla granicy fazy Kc = 1/(2p(1−p)) jest nałożone na diagram fazowy K–p, a eksperymenty perturbacji pojedynczego węzła bezpośrednio badają granicę porządek-chaos.

Nauki o Życiu

Fala Populacyjna — 2D Fisher-KPP, Float32Array, Efekt Allee

2D równanie Fishera-KPP opisuje przestrzenną inwazję gatunku na krajobrazie. Ta symulacja rozwiązuje je na siatce Float32Array za pomocą schematu różnic skończonych, z gęstością populacji u(x,y,t) renderowaną jako mapa w fałszywych kolorach. Efekt Allee — zmniejszony wzrost per capita przy niskiej gęstości — jest wbudowany poprzez zmodyfikowany wyraz wzrostu:

f(u) = ru(u−a)(1−u),   0 < a < 1

Poniżej progu Allee a, lokalne populacje maleją zamiast rosnąć, tworząc minimalny opłacalny rozmiar płatu i krytyczną prędkość fali, poniżej której inwazja się nie udaje. Porównanie prostego wzrostu logistycznego z bistabilnym modelem Allee pokazuje, jak krytyczny próg tworzy ostre granice inwazji i dług wymierania we fragmentowanych krajobrazach.

Sieci

Perkolacja Sieciowa — Erdős-Rényi, Union-Find, Gigantyczna Składowa

Dla losowego grafu Erdősa-Rényiego G(N,p), gigantyczna spójna składowa o rozmiarze O(N) pojawia się ostro przy średnim stopniu ⟨k⟩ = pN = 1. Symulacja dodaje krawędzie jedna po drugiej i śledzi rozmiary składowych za pomocą Union-Find z kompresją ścieżki i unią wg rangi, osiągając niemal-O(1) zamortyzowany czas scalania. Frakcja gigantycznej składowej S = |Cmax|/N jest rysowana wobec ⟨k⟩ w czasie rzeczywistym, z nałożonym przewidywaniem samospójności pola średniego S = 1 − exp(−⟨k⟩S). Węzły są skalowane proporcjonalnie do swojej składowej, czyniąc moment, w którym gigantyczna składowa się uruchamia, wizualnie dramatycznym.

Sztuka Generatywna

Wzrost Losowy — Fraktal Eden + DLA, D ≈ 1,71 Zliczanie Pudełkowe

Model Edena rozwija zwarty klaster, losowo wybierając miejsce na obwodzie; wynikiem jest w przybliżeniu okrągła forma z szorstką, ale niefraktalną powierzchnią. Agregacja ograniczona dyfuzją (DLA) wypuszcza losowych wędrowców z daleka i dodaje każdego tam, gdzie po raz pierwszy dotknie, produkując cienką, rozgałęzioną strukturę o wymiarze fraktalnym D ≈ 1,71 w dwóch wymiarach. Symulacja uruchamia oba modele obok siebie i oblicza wymiar fraktalny metodą zliczania pudełkowego w czasie rzeczywistym, rysując log N(ε) wobec log(1/ε) z dopasowaniem regresyjnym, które zbiega do teoretycznego wymiaru fraktalnego w miarę wzrostu klastra. Wizualny kontrast między zwartym klastrem Edena a pajęczym agregatem DLA sprawia, że związek między lokalnymi regułami wzrostu a wielkoskalową geometrią fraktalną jest natychmiast oczywisty.

Nauka o Kosmosie

Silnik Halla — Dryf E×B, Isp, Ciąg

Silniki na efekcie Halla to dominująca technologia napędowa do utrzymywania pozycji satelitów i podnoszenia orbit. Jony ksenonu są przyspieszane przez pole elektryczne E do prędkości wylotowych 15–30 km/s (Isp 1500–3000 s), podczas gdy elektrony są uwięzione w zamkniętych orbitach dryfu Halla v = E×B/B² przez poprzeczne pole magnetyczne, które uniemożliwia im zneutralizowanie przyspieszającego potencjału.

Symulacja modeluje trajektorie cząstek obu gatunków za pomocą pushera Borisa w skrzyżowanych polach kanału wyładowania, obliczając uśredniony w czasie ciąg i impuls właściwy. Animowane pierścienie dryfu elektronowego wizualizują prąd Halla. Suwaki kontrolują napięcie wyładowania, tempo przepływu masy propelentu i geometrię kanału; ciąg F = ṁve i Isp = ve/g0 aktualizują się nieprzerwanie, czyniąc konkretnym, jak małe zmiany napięcia produkują duże zmiany Isp przy skromnym koszcie ciągu.

Kaskada Promieniowania Kosmicznego — Kaskada Heitlera, e⁻/γ/μ/π, Profil NKG

Gdy wysokoenergetyczny proton promieniowania kosmicznego uderza w górną atmosferę, inicjuje rozległą kaskadę powietrzną. Binarny model kaskadowy Heitlera opisuje składową elektromagnetyczną: każdy elektron emituje foton bremsstrahlung, a każdy foton tworzy parę elektron-pozyton, podwajając liczbę cząstek co długość radiacyjną X0, aż energie spadną poniżej energii krytycznej Ec, a kaskada osiągnie maksymalną głębokość Xmax.

Symulacja generuje rozgałęziające się drzewo kaskadowe z energii pierwotnej 1–100 PeV, śledząc elektrony, fotony, miony i piony w odrębnych kolorach. Profil kaskady (liczba cząstek wobec głębokości atmosferycznej) podąża za funkcją Nishimury-Kamaty-Greisena (NKG), nałożoną na wynik Monte Carlo. Głębokość maksimum kaskady Xmax przesuwa się wraz z masą pierwotną, demonstrując, dlaczego służy jako sonda składu: kaskady protonowe są głębsze i bardziej fluktuujące niż kaskady żelaza o równej energii całkowitej.

Elektronika

Prawa Kirchhoffa — Macierz MNA, Eliminacja Gaussa, 4 Presety Obwodów

Zmodyfikowana Analiza Węzłowa (MNA) koduje prawa prądowe i napięciowe Kirchhoffa jednocześnie jako system macierzowy Ax = b:

[G  B] [v]   [i]
[C  D] [j] = [e]

gdzie G koduje przewodności, B i C kodują ograniczenia źródeł napięcia, v zawiera napięcia węzłowe, j prądy gałęziowe przez źródła napięcia, a e wartości źródeł. Symulacja buduje macierz MNA dla czterech presetowych obwodów — dzielnika napięcia, mostka Wheatstone'a, sieci drabinkowej i obwodu wielo-źródłowego z zależnymi źródłami — i rozwiązuje eliminacją Gaussa z częściowym pivotowaniem. Napięcia węzłowe są opisane na schemacie; prądy gałęziowe pokazane są jako skierowane strzałki, których szerokość jest proporcjonalna do wielkości. Dostosowanie dowolnej wartości komponentu natychmiast aktualizuje rozwiązanie, czyniąc zależność zachowania obwodu od stosunków komponentów namacalną, a nie algebraiczną.

Wypróbuj Sam

Sześć symulacji na początek, jeśli jesteś nowy w tej fali:

Pełny katalog Wave 109 jest przeglądalny przez indeks kategorii. Wszystkie dwadzieścia symulacji jest też dostępnych po ukraińsku przez prefiks językowy /uk/.

Co Dalej

Do celu 1000 symulacji brakuje 60. Wave 110 go osiągnie. Obecne plany obejmują symulacje blokady atomów Rydberga, turbulencji warstwy granicznej atmosfery, stanów brzegowych izolatora topologicznego, transportu optymalnego (odległość Ziemia-mover), pejzaży energetycznych zwijania białek oraz powrót do atraktora Lorenza z pełnym obliczeniem widma Lapunowa. Nacisk cały czas kładziony jest na symulacje, które czynią widoczną przy interaktywnej prędkości konkretną ilościową przewidywanie: nie tylko animacje, ale narzędzia do budowania fizycznej intuicji o rządzących równaniach.

Jeśli masz symulację, którą chciałbyś zobaczyć zbudowaną, lub model z twojej dziedziny badań, który jest niedostatecznie reprezentowany w interaktywnej nauce przeglądarkowej, strona kontaktowa jest otwarta. Sugestie badaczy i nauczycieli już ukształtowały kilka fal, a planowanie Wave 110 nie jest jeszcze zamknięte.