Gaz Latticowy FHP: Navier-Stokes Z Prawdziwych/Fałszywych Bitów na Heksagonie

Dynamika płynu z siatki bitów true/false

Model FHP – nazwany na cześć Uriela Frisha, Brosla Hasslachera i Yves Pomeau, którzy go wprowadzili w 1986 roku – jest jednym z najbardziej niezwykłych podejść do dynamiki płynów: odzyskuje równania Navier-Stokes o nieskonsecjowanej fazie z siatki cząstek logicznych, które przestrzegają absurdalnie prostych lokalnych zasad, bez jakiegokolwiek pola ciśnienia lub prędkości zmiennoprzecinkowej w modelu. Każde pole na dwuwymiarowej, sześciokątnej siatce może zawierać maksymalnie jeden cząstek poruszającą się w jednym z sześciu możliwych kierunków (sześć krawędzi heksagonu), więc stan całej układu jest reprezentowany jedynie przez kilka bitów na pole.

Dwie kroki, powtarzane w nieskończoność: strumień i zderzenie

W każdym kroku czasowym następują dwie operacje lokalne. W kroku propagacji (strumienia) każdy cząstek przesuwa się o jedną komórkę siatki w kierunku, w którym już się porusza, lądując na odpowiadającym mu sąsiednim elemencie. W kroku zderzenia cząstki, które właśnie dotarły do tego samego elementu, są przekierowywane zgodnie z ustaloną tabelą reguł, która zachowuje zarówno liczbę cząstek, jak i pęd całkowity w tym elemencie – na przykład dwa cząstki uderzające frontalnie (poruszające się w przeciwnych kierunkach) stanowią klasyczny przypadek kolizji FHP, a reguła rozprasza je w dwóch kierunkach prostopadłych do ich pierwotnej linii ruchu, wybierane losowo monetą, aby uniknąć nieprawdziowego uprzedzenia.

// one FHP time step, conceptually
propagate();   // every particle moves 1 site in its direction of travel
collide();     // at each site: apply the local rule table
               //   - conserves particle count at that site
               //   - conserves total momentum at that site
               //   - the specific outgoing directions are looked up from
               //     a small table indexed by which incoming directions
               //     were occupied

Dlaczego siatka sześciokrotna, a nie kwadratowa?

Wcześniejsza siatka gazu na kwadratowej planszy (model HPP, 1976) posiada poważny błąd: jej równowagowe rozkład jest niewystarczająco izotropiczny, a czwartorzędowy tensor występujący w wyprowadzeniu równania momentu pędu makroskopowego nie redukuje się do postaci symetrycznej względem obrotów, której potrzebuje Navier-Stokes – czwórkowa symetria siatki kwadratowej przecieka do zachowania płynu jako zależność od kierunku. Siatka sześciokrotnia FHP ma symetrię sześciokrotną, która okazuje się być dokładnie wystarczająca dyskreta rotacyjna symetria dla danego tensora, aby stał się izotropiczny w granicy ciągłej. To był konkretny wgląd, który sprawił, że praca FHP odniosła sukces, gdzie model HPP zawiódł – wybór geometrii siatki nie jest czysto kosmetyczny, a determinuje poprawność emergentnej fizyki.

Przybliżone granulowanie do równań Naviera-Stocksa

Żaden pojedynczy cząstek nie przestrzega niczym zbliżnie równania płynu – fizyka pojawia się dopiero po przybliżonym granulowaniu: uśrednienie zajętości cząstek na wielu węzłach siatki i przez wiele kroków czasowych, aby uzyskać gładko zmieniające się pole gęstości i prędkości, a to uśrednione pole przestrzega równań Naviera-Stocksa niepłynnych do korekty. Ta korekta jest subtelnością, którą warto znać: pochodzone przez FHP równanie zachowania pędu nie jest idealnie Galileuszowskie, ponieważ rozkład równowagi gazu w węzłach siatki nie jest dokładnie Maxwellowskiem. Wyraża się to jako współczynnik mnożący g(ρ) zależny od gęstości przed wyrazem nieliniowym adwersyjnym, który należy skorygować (analitycznie lub numerycznie), aby odzyskać dokładną formę równań Naviera-Stocksa zamiast przybliżenia o smaku gazu w węzłach siatki.

Od gazu na siatce do metody Boltzmanna dyskretnej

Liczby zajętości FHP oparte na wartościach logicznych powodują w modelu inherentny szum — każda z pozyskiwanych zmiennych makroskopowych jest w rzeczywistości średnią statystyczną obliczoną na podstawie liczby cząstek, a uzyskanie gładkiego pola prędkości wymaga uśredniania przestrzennego lub czasowego po wielu węzłach lub po wielu uruchomieniach, co jest kosztowne. Metoda Boltzmanna dyskretnej (LBM), opracowana w latach następujących po FHP, zastępuje liczby zajętości logicznych oparte na wartościach logicznych przez rzeczywiste funkcje rozkładu pojedynczych cząstek — zasadniczo wstępnie uśrednia populację logiczną bezpośrednio zamiast uśredniać ją po fakcie nad zespołem symulacji logicznych. LBM zachowuje podstawową ideę FHP (stała siatka, lokalne przepływy i zasady kolizji, odzyskiwanie Naviera-Stokes w kontinuum), ale eliminuje szum statystyczny, a właśnie LBM, a nie oryginalny model FHP oparty na wartościach logicznych, znalazł zastosowanie w poważnych obliczeniach płynowych dla złożonych geometrii, mediów porywistych i przepływów wielofazowych.

Frequently asked questions

Wypróbuj na żywo

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