Plazma to gaz zjonizowany na tyle, że jego zachowanie zdeterminowane jest przez dalekozasięgowe siły elektromagnetyczne, a nie krótkozasięgowe zderzenia cząstek obojętnych. Zerwij elektrony z atomów, a zostanie ci mieszanina swobodnych elektronów i dodatnich jonów, która reaguje zbiorowo: szturchnij jeden obszar, a całe medium odpowie poprzez pola elektryczne. Właśnie ta zbiorowa natura sprawia, że plazma jest zarazem fascynująca i trudna do wyobrażenia sobie na podstawie samych równań. Poniższe modele poświęcają część fizycznej kompletności na rzecz interaktywności, dzięki czemu możesz przeciągnąć suwak i obejrzeć, jak naładowany gaz reorganizuje się w czasie rzeczywistym.
Ekranowanie Debye'a i skala długości plazmy
Pierwszą wielkością, jakiej potrzebuje każdy model plazmy, jest długość
Debye'a, zapisywana jako λ_D. Umieść pojedynczy dodatni
ładunek próbny w plazmie, a ruchliwe elektrony napłyną ku niemu rojem, podczas gdy jony
odpłyną. Na niewielkiej odległości pole ładunku próbnego zostaje niemal całkowicie
zniesione — plazma je „wyekranowała". Długość Debye'a to charakterystyczny promień,
w obrębie którego zachodzi to ekranowanie, wyznaczony przez równowagę między ruchem
termicznym, który dąży do randomizacji położeń, a przyciąganiem elektrostatycznym, które
dąży do ich uporządkowania. Formalnie
λ_D = sqrt(ε_0 k_B T_e / (n_e e²)): gorętsza plazma
ekranuje na większych odległościach, gęstsza — na mniejszych.
Długość Debye'a to nie akademicka ciekawostka. Wyznacza ona granicę między dwoma
reżimami. W skalach mniejszych niż λ_D cząstka odczuwa nagie pola
kulombowskie swoich sąsiadów i zachowuje się niemal jak zbiór pojedynczych ładunków. W
skalach większych niż λ_D te pola są wyekranowane, a medium zachowuje
się jak gładki, quasi-neutralny płyn. Użytecznym testem poprawności dla dowolnego gazu
zjonizowanego jest liczba cząstek wewnątrz sfery Debye'a — dopiero gdy ta liczba
jest duża, obowiązuje statystyczny, zbiorowy obraz i system w ogóle zasługuje na nazwę
„plazma".
Dlaczego to ważne dla symulacji: każdy siatkowy solver plazmy musi
rozdzielić długość Debye'a, inaczej fizyka ekranowania zostaje utracona i pojawia się
pozorne nagrzewanie. Modele przeglądarkowe idą na skróty, pracując w jednostkach
przeskalowanych, gdzie jeden piksel to ustalona część λ_D, co
utrzymuje stabilność dynamiki nawet na laptopie.
Drgania plazmowe: puls naładowanego gazu
Przesuń nieco populację elektronów w plazmie w jedną stronę, a pojawi się siła
przywracająca: przesunięte elektrony zostawiają za sobą warstwę dodatnich jonów, powstałe
pole elektryczne ściąga je z powrotem, przeskakują położenie równowagi, a cała chmura
zaczyna „dzwonić". To drganie plazmowe, a jego częstotliwość —
elektronowa częstotliwość plazmowa ω_p — zależy wyłącznie od
gęstości elektronów: ω_p = sqrt(n_e e² / (ε_0 m_e)).
Gęstsza plazma dzwoni szybciej. Ponieważ jony są tysiące razy cięższe od elektronów, w
tej skali czasowej praktycznie się nie poruszają, więc drganie jest zasadniczo zjawiskiem
czysto elektronowym.
Częstotliwość plazmowa wyznacza też granicę odcięcia dla fal elektromagnetycznych.
Nadchodząca fala o częstotliwości poniżej ω_p nie może się propagować
— elektrony reagują wystarczająco szybko, by ją znieść, i fala zostaje odbita.
Dlatego jonosfera odbija fale radiowe AM wokół krzywizny Ziemi i dlatego korona Słońca
jest nieprzezroczysta dla pewnych częstotliwości radiowych. To samo odcięcie tłumaczy
błyszczący, lustrzany wygląd metali, których elektrony przewodnictwa tworzą zdegenerowaną
plazmę. Nasz model drgań pozwala podnosić gęstość i obserwować, jak jednocześnie rośnie
częstotliwość „dzwonienia" i poszerza się pasmo odbicia.
Niestabilności: gdy małe zaburzenia wymykają się spod kontroli
Plazma jest znana ze swoich niestabilności — konfiguracji, w których drobne zaburzenie żeruje na wolnej energii układu i rośnie wykładniczo. Niestabilność dwóch strumieni to podręcznikowy przykład: przepuść dwie wiązki naładowanych cząstek jedną przez drugą, a każda mała fala gęstości grupuje cząstki w skupiska, co wzmacnia pole, które grupuje je jeszcze mocniej. Energia przepływa z ruchu strumieniowego do wzrostu fal, aż wiązki rozpadają się na turbulentne wiry. To plazmowy kuzyn sprzężenia zwrotnego (feedbacku) w mikrofonie.
Niestabilność Kelvina-Helmholtza pojawia się wszędzie tam, gdzie dwa płyny przesuwają się względem siebie ścinaniem, a w namagnesowanej plazmie granica między szybkim a wolnym przepływem zwija się w łańcuch wirów. Widać ją na krawędzi magnetosfery Ziemi, gdzie wiatr słoneczny prześlizguje się obok, w pasiastych chmurach Jowisza i w laboratoryjnych strugach plazmowych. Fizyka jest identyczna z zawirowaniami, które powstają, gdy wiatr wieje nad wodą, dlatego właśnie nasza symulacja fali kilwaterowej Kelvina jest tak przydatnym towarzyszem: wersja z płynem obojętnym czyni geometrię oczywistą, zanim dodasz siły elektromagnetyczne. Dla przypadku ścinania namagnesowanego symulacja niestabilności plazmy zasiewa warstwę ścinania prędkości i pozwala obserwować, jak wiry rosną, łączą się i nasycają.
Niestabilności to również miejsce, w którym fizyka plazmy spotyka się z chaosem. Gdy kilka niestabilnych modów zacznie ze sobą oddziaływać, trajektorie stają się wrażliwe na warunki początkowe, a przepływ staje się turbulentny. Te same matematyczne sygnatury — rozciąganie i składanie w przestrzeni fazowej — pojawiają się w symulacji atraktora Lorenza, która jest okrojonym modelem turbulencji konwekcyjnej i łagodnym wprowadzeniem w dynamikę, przez którą utrzymanie plazmy w reaktorach fuzyjnych jest tak trudne do kontrolowania.
Energia, zderzenia i nagrzewanie
Prawdziwa plazma rzadko znajduje się w porządnej równowadze termicznej. Energia wprowadzona przez wiązkę, falę uderzeniową lub niestabilność musi rozłożyć się między cząstkami, a tempo, w jakim to następuje, zależy od zderzeń. W gorącej, rozrzedzonej plazmie zderzenia są rzadkie i gaz może długo utrzymywać bardzo różne temperatury elektronów i jonów. W gęstszym, ziarnistym reżimie energia rozprzestrzenia się szybko, a rozkład zbliża się do jednej temperatury. Dynamikę rozprzestrzeniania się energii przez gęsty zbiór zderzających się cząstek bezpośrednio bada symulacja nagrzewania granularnego, która — choć modeluje obojętne ziarna — oddaje tę samą intuicję statystyczno-mechaniczną: energia kinetyczna kaskadowo przechodzi z ruchu uporządkowanego w ciepło poprzez niezliczone drobne zderzenia.
W przypadku plazmy właśnie to tłumaczy, dlaczego nagrzewanie urządzenia fuzyjnego jest trudne. Wpompuj energię w elektrony, a mogą wypromieniować ją, zanim dotrze do jonów, które faktycznie musisz doprowadzić do syntezy. Każdy praktyczny schemat nagrzewania — omowy, radiowy, wiązką neutralną — to strategia sprzęgania energii z właściwą populacją cząstek, zanim plazma straci ją na rzecz ścianek lub promieniowania.
Przeniesienie do przeglądarki
Pełna symulacja plazmy wykorzystuje metody cząstka-w-komórce: śledzi się miliony makrocząstek, ich ładunek osadza się na siatce, pola rozwiązuje się solverem Poissona lub Maxwella, a następnie popycha cząstki wynikowymi siłami. To zdecydowanie zbyt wiele dla telefonu, więc modele przeglądarkowe wprowadzają świadome uproszczenia — zredukowaną liczbę cząstek, dwa wymiary, przeskalowane jednostki i jawne całkowanie w czasie z krokiem wystarczająco małym, by pozostać poniżej częstotliwości plazmowej. Kompromis jest uczciwy: tracisz dokładność ilościową, ale zachowujesz fizykę jakościową, a to właśnie ona buduje intuicję.
Fizyka plazmy nagradza podejście interaktywne bardziej niż niemal każda inna gałąź fizyki, ponieważ tak wiele w jej zachowaniu jest zbiorowe i emergentne. Żadna pojedyncza cząstka nie „decyduje" o utworzeniu powłoki Debye'a ani zwinięciu się w wir; struktura powstaje z interakcji milionów cząstek. Zwolnij tę interakcję, pokoloruj pola i pozwól czytelnikowi poruszać parametrami, a pojęcia zajmujące cały rozdział gęstej algebry stają się czymś, co po prostu widać. Otwórz dowolną z powyższych symulacji, zmieniaj po jednym suwaku naraz, a abstrakcje czwartego stanu materii zaczną wydawać się konkretne.