Spotlight #52 – Geometria, inteligencja rojowa i dynamika płynów

Trzy zbadane kategorie

mysimulator.uk śledzi głębię kategorii jako jeden z mierników jakości. Fala 55 wzmacnia trzy odrębne dziedziny — od ponadczasowej matematyki wielościanów po fizykę obliczeniową płynów opartych na cząstkach.

🔷 Geometria dyskretna — wielościany i wzór Eulera

Badanie wielościanów — brył ograniczonych płaskimi ścianami wielokątnymi — sięga od Platońskiego „Timajosa” (ok. 360 r. p.n.e.), który utożsamił cztery bryły z klasycznymi żywiołami, przez zaskakujące twierdzenie topologiczne Eulera, aż po współczesne zastosowania w krystalografii i grafice komputerowej.

Dlaczego istnieje tylko pięć brył platońskich?

Bryła platońska wymaga, aby (a) wszystkie ściany były tym samym wielokątem foremnym oraz (b) w każdym wierzchołku spotykała się taka sama liczba ścian. Aby bryła mogła się zamknąć, suma kątów w każdym wierzchołku musi być mniejsza niż 360°, w przeciwnym razie ściany leżą płasko lub zawijają się na siebie. Istnieje dokładnie pięć rozwiązań:

• Czworościan — 3 trójkąty równoboczne na wierzchołek, suma kątów 180°.
• Sześcian (heksaedr) — 3 kwadraty na wierzchołek, suma kątów 270°.
• Ośmiościan — 4 trójkąty równoboczne na wierzchołek, suma kątów 240°.
• Dwunastościan — 3 pięciokąty foremne na wierzchołek, suma kątów 324°.
• Dwudziestościan — 5 trójkątów równobocznych na wierzchołek, suma kątów 300°.

Sześć trójkątów równobocznych wypełniłoby płaską płaszczyznę (360°) i nie zamknęłoby się; pięć kwadratów przekroczyłoby 360°. Ograniczenie jest ścisłe — nie ma miejsca na szóstą bryłę platońską.

Charakterystyka Eulera: V − E + F = 2

W 1752 roku Leonhard Euler zauważył, że dla każdego wielościanu wypukłego wartość V−E+F wynosi dokładnie 2. Było to zaskakujące, ponieważ zależność ta zachodziła niezależnie od kształtu: czworościan (4−6+4=2), sześcian (8−12+6=2), dwudziestościan (12−30+20=2), a nawet dwudziestościan ścięty (60−90+32=2).

Bryły archimedesowe i ścinanie

13 brył archimedesowych to wypukłe, wierzchołkowo-przechodnie wielościany o dwóch lub więcej rodzajach ścian będących wielokątami foremnymi. Powstają naturalnie z brył platońskich poprzez ścinanie (odcinanie wierzchołków), rektyfikację (cięcie do środków krawędzi), łamanie (skręcanie) oraz kantelację (rozsuwanie ścian na zewnątrz). Dwudziestościan ścięty — 60 wierzchołków, 90 krawędzi, 32 ściany w postaci pięciokątów i sześciokątów — jest identyczny ze wzorem na standardowej piłce nożnej oraz z cząsteczką fulerenu C₆₀ (buckminsterfullerenu).

🐟 Inteligencja rojowa — emergencja z lokalnych reguł

Inteligencja rojowa opisuje zbiorowe, samoorganizujące się zachowanie, które wyłania się z lokalnych interakcji między prostymi agentami. Żaden centralny kontroler nie kieruje rojem — każda jednostka reaguje wyłącznie na informacje z najbliższego otoczenia. Mimo to wzorce na poziomie grupy — spójne stada, rozgałęziające się sieci szlaków, skoordynowane reakcje paniki — mogą dorównać efektom działania dowolnego systemu scentralizowanego.

Boids Craiga Reynoldsa (1986)

Oryginalna praca Reynoldsa o Boids (SIGGRAPH '87) pokazała, że trzy reguły wystarczą, aby uzyskać realistyczne flokowanie:

Trzy promienie — r_sep < r_ali < r_coh — tworzą zagnieżdżone strefy wpływu. Szkoła ryb 3D wykorzystuje r_sep=28, r_ali=55, r_coh=80 px z trzema warstwami głębi, dającymi pseudo-3D wygląd przy koszcie Canvas 2D.

Dynamika drapieżnik–ofiara w stadach

Nowa symulacja Szkoła ryb 3D dodaje polującego rekina, który wywołuje reakcję paniki u ryb znajdujących się w promieniu jego patrolu. Odtwarza to kilka udokumentowanych zjawisk obserwowanych w terenie:

• Efekt fontanny: ławica rozdziela się wokół drapieżnika i łączy ponownie za nim, tworząc ciągły ruch, który dezorientuje śledzenie sensoryczne drapieżnika.
• Błyskawiczna ekspansja: gdy drapieżnik atakuje, najbliższe ryby przyspieszają gwałtownie na zewnątrz — natychmiastowy spadek gęstości zmusza drapieżnika do skupienia się na jednym celu.
• Efekt rozcieńczenia: prawdopodobieństwo upolowania danej jednostki wynosi 1/N dla ławicy liczącej N ryb, co sprawia, że większe ławice są bezpieczniejsze w przeliczeniu na rybę, mimo że są bardziej rzucające się w oczy.
• Efekt dezorientacji: szybki ruch wielu identycznych ofiar przeciąża system śledzenia obiektów drapieżnika. Eksperymenty z gupikami pokazują, że drapieżniki łapią mniej ryb z większych ławic, nawet gdy częstość ataków pozostaje stała.

💧 Obliczeniowa dynamika płynów — SPH

Obliczeniowa dynamika płynów (CFD) to numeryczna symulacja równań Naviera-Stokesa. Dominują dwie rodziny metod: eulerowskie (oparte na siatce: objętości skończonych, elementów skończonych, różnic skończonych) oraz lagranżowskie (oparte na cząstkach: SPH, DEM, MPM). SPH należy do rodziny lagranżowskiej i posiada właściwości, których brakuje metodom siatkowym.

Dlaczego podejście lagranżowskie?

Metody siatkowe mają trudności z problemami obejmującymi duże deformacje, powierzchnie swobodne i przepływy wielofazowe: siatka musi być regenerowana lub deformowana na każdym kroku. W SPH cząstki niosą przepływ ze sobą — powierzchnie swobodne są śledzone automatycznie, siatka nie jest potrzebna, a fragmentacja jest trywialna. Cena za to jest taka, że SPH w naiwnej implementacji ma złożoność O(N²) (choć haszowanie przestrzenne redukuje ją do O(N log N)) i jest mniej dokładne niż metody elementów skończonych wysokiego rzędu dla przepływów gładkich.

Układ ciśnienie–lepkość–grawitacja

Symulator całkuje trzy składowe siły w każdym podkroku:

• Ciśnienie — wynikające z lokalnych nadwyżek gęstości ponad gęstość spoczynkową ρ₀. Wykorzystuje symetryzowany gradient jądra typu spiky: F_p = −m (pᵢ+pⱼ)/2ρⱼ ∇W_spiky
• Lepkość — dyfuzja pędu, wygładzająca różnice prędkości: F_ν = νm(vⱼ−vᵢ)/ρⱼ ∇²W_visc
• Grawitacja — stała siła masowa skierowana w dół, działająca na wszystkie cząstki.

Sprężyste odbicie od granic ze współczynnikiem tłumienia 0.4 modeluje pochłanianie energii przez rzeczywiste betonowe ściany. Trzy podkroki na klatkę animacji (dt=0,004 s każdy) utrzymują spełnienie warunku stabilności CFL przy 60 kl./s.