DSMC i Gaz Rzadki: Symulacja Gazu, Którego Nie Da się Przewidzieć

Kiedy gaz przestaje zachowywać się jak ciecz

Równania Naviera-Stokesa opisujące zwykły przepływ płynów opierają się na ukrytej założeniu: że gaz jest wystarczająco gęsty, aby cząsteczka mogła zderzyć się z sąsiadami częściej niż zmieniają się w tym samym obszarze własne cechy przepływu. Rozważ niezwykle rzadki gaz – atmosferę górną podczas ponownego wejścia statków kosmicznych, próżnię wewnątrz komory próżniowej lub gaz wokół urządzenia MEMS o wymiarach zaledwie mikrometry – i to założenie ulega zmianie, a równania Naviera-Stokesa przestają dawać poprawne wyniki, niezależnie od tego, jak starannie je rozwiązujesz.

Licencja liczby Knudsena decyduje o reżimie

Niezmienna liczba bezwymiarowa, która mówi Ci, w jakim reżimie się znajdujesz, to Licencja Knudsena, Kn = lambda / L, gdzie lambda to średni rozstryż molekularny (średnia odległość, jaką pokonuje cząstka między zderzeniami), a L to charakterystyczna długość przepływu, np. rozmiar obiektu lub szerokość kanału.

Kn < 0.01 : przepływ ciągły - Navier-Stokes jest dokładny (cząsteczki uderzają w siebie stale w stosunku do L) 0.01 < Kn < 0.1 : przepływ z tarciem - Navier-Stokes działa, jeśli dodasz warunek tarcia/temperatury 0.1 < Kn < 10 : przepływ przejściowy - ani równania ciągłości PDE, ani granica cząsteczkowa nie mają zastosowania; potrzebujesz pełnej równania Boltzmanna Kn > 10 : przepływ cząsteczkowy - cząsteczki zasadniczo nigdy nie zderzają się ze sobą, tylko z powierzchniami W reżimie przejściowym dokładnym równaniem rządzącym jest równanie Boltzmanna, które śledzi pełną rozkład prawdopodobieństwa pozycji i prędkości molekularnych, zamiast średnich wartości (gęstość, prędkość obrotowa, temperatura), które wykorzystuje Navier-Stokes. Równanie Boltzmanna jest znany z trudności w bezpośrednim rozwiązywaniu dla realistycznych trójwymiarowych geometrii ze względu na wysokowymiarowy integral kolizji, który dokładnie wypełniał lukę Symulacja Monte Carlo Bezpośredniej.

Kn <  0.01   :  continuum flow      - Navier-Stokes is accurate (molecules collide constantly relative to L)
0.01 < Kn < 0.1  :  slip flow      - Navier-Stokes works if you add a velocity/temperature slip boundary condition
0.1  < Kn < 10  :  transition flow  - neither continuum PDEs nor free-molecular limits apply; need the full Boltzmann eq.
Kn >  10   :  free molecular flow  - molecules essentially never collide with each other, only with surfaces

Jak działa DSMC

DSMC, opracowany przez Graeme Bird’ego począwszy od lat 60. XX wieku, w ogóle nie rozwiązuje równania Boltzmanna jako równania różniczkowego – symuluje je statystycznie. Zamiast śledzić każdą z astronomicznie wielu rzeczywistych cząstek gazowych, DSMC śledzi mniejszą liczbę symulowanych cząstek, z których każda reprezentuje dużą liczbę rzeczywistych cząstek, a następnie aktualizuje je w powtarzalnych cyklach dwóch oddzielonych kroków:

dla każdego kroku czasowego: 1. PRZESUNIĘCIE: każda symulowana cząstka porusza się prostoliniowo na dystans dt (ruch bez kolizji, uwzględniając wszelkie interakcje z powierzchniami/granicami) oraz 2. KOLIZJA: w każdej komórce przestrzennej losuje się pary cząstek i pozwala się im uderzyć prawdopodobieństwie (prawdopodobieństwo akceptacji proporcjonalne do względnej prędkości x przekroje czynnego zderzenia – popularnym i wydajnym wyborem jest schemat „No-Time-Counter”) Podział ruchu i kolizji na oddzielne kroki jest możliwy tylko wtedy, gdy krok czasowy jest mniejszy niż średni czas kolizji lokalnej oraz komórki przestrzenne są mniejsze niż środek swobodny, w przeciwnym razie cząstki, które powinny uderzyć się nawzajem w trakcie podróży, mijają się bez przeszkód. Ponieważ kolizje są obsługiwane statystycznie w komórce zamiast sprawdzania każdego możliwego przypadku zderzenia cząstek, DSMC skaluje się znacznie lepiej niż dosłowne symulacje N-ciała molekularnego, a jednocześnie odzyskuje prawidłowe zachowanie gazu na makroskali (gęstość, prędkość, temperatura, nawet niehomogeniczne efekty takie jak dwumodelowa lub skośna dystrybucja prędkości) jako emergentny wynik statystyczny zamiast założonego.

for each time step:
    1. MOVE:    every simulated particle streams in a straight line for dt (collision-free motion, apply any boundary/surface interactions)
    2. COLLIDE: within each spatial cell, sample pairs of particles and let them collide probabilistically
               (acceptance probability proportional to relative speed x collision cross-section - the "No-Time-Counter" scheme is a common efficient choice)

Co uchwycają DSMC, czego nie potrafią równania Navier-Stokesa

Ponieważ DSMC pracuje bezpośrednio z mechaniką na poziomie cząsteczkowym, naturalnie wychwytuje efekty niestosunkowości – sytuacje, w których rozkład prędkości gazu nie jest gładką dystrybucją Maxwell-Boltzmanna, którą zakładają równania Navier-Stokesa – takie jak silne skoki temperatury i prędkości na froncie rozszczepienia się pęcherza, lub sposób, w jaki cząsteczki gazu odbijają się od powierzchni stałej z różną efektywną temperaturą (nawet opuszczając ją z preferowanym kierunkiem) zamiast idealnie oddewajania się po kontakcie. Te efekty dokładnie pojawiają się w obliczeniach nagrzewania podczas powrotu kapsuły kosmicznej, ponieważ warstwa rozszczepienia się pęcherza wokół wracającej kapsuły na dużych wysokościach znajduje się głęboko w obszarze rzadkiego/przejściowym przed tym, gdy atmosfera staje się gęsta i można zastosować równania ciągłe.

Gdzie liczą się dynamika rozrzedzonych gazów

Poza aerodynamiką reentry, projektowaniem układów próżniowych, analizą przepływów gazów w systemach MEMS oraz modelowaniem kolizyjnych obszarów atmosfer planet i ohonów, analiza DSMC (Direct Simulation Monte Carlo) i analiza rozrzedzonych gazów są niezbędne wszędzie tam, gdzie skale długości i ciśnienia stają się bardzo małe: projektowanie dysz jonotrzaskowych i plazmowych, ekspansja obłoków strumieniowych jonotrzaskowych w próżni bliskim otoczeniu kosmicznym, przepływ gazów w mikroelektromechanicznych systemach (MEMS) o szerokości kanałów rzędu kilku mikrometrów oraz zewnętrzne, kolizyjnie wolne obszary atmosfer planet i ohonów, gdzie gazy uciekają w przestrzeń w trybie swobodnego molekuł, a nie jako ciągły płyn.

Często zadawane pytania

Wypróbuj na żywo

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