Про симуляцію
Нестабільність Кельвіна-Гельмгольца (КГ) у плазмі виникає щоразу, коли два намагнічені шари плазми рухаються повз один одного з різними швидкостями. Розрив швидкості на їхній межі діє як вихровий лист: будь-яке мале збурення посилюється дисбалансом тиску Бернуллі між швидкою та повільною сторонами, змушуючи межу згортатися у характерні вихори типу «котяче око». Лінійна швидкість зростання становить γ = k·ΔV·√(ρ₁ρ₂)/(ρ₁+ρ₂), тому моди з більшим хвильовим числом (коротші довжини хвилі) зростають найшвидше — аж до нелінійного насичення та злиття вихорів, які переносять енергію до більших масштабів.
Ця симуляція моделює межу плазми за допомогою регуляризованого вихрового листа (рівняння Біркгофа-Ротта з десингуляризацією вихрових сгустків). N маркерних точок, розподілених уздовж межі, несуть циркуляцію, пропорційну місцевому стрибку швидкості ΔV. Їхні позиції розвиваються за дискретизованим ядром Біо-Савара, який просувається в часі за допомогою інтегратора Рунге-Кутта. Панель дисперсії під основним полотном показує теоретичну криву швидкості зростання γ(k) з маркером для поточної вибраної моди. Керуйте зсувом швидкості, відношенням густин між двома шарами плазми та початковим хвильовим числом k за допомогою повзунків на панелі HUD.
Часті запитання
Що таке нестабільність Кельвіна-Гельмгольца в плазмі?
Нестабільність Кельвіна-Гельмгольца (КГ) виникає на межі між двома шарами плазми, що рухаються з різними швидкостями. Зсув швидкості на кордоні посилює малі збурення: тиск падає там, де течія швидша (ефект Бернуллі), що ще більше викривляє межу і породжує вихори. У намагніченій плазмі магнітне поле може стабілізувати або дестабілізувати межу залежно від орієнтації поля відносно напрямку течії.
Що означає формула швидкості зростання γ = k·ΔV/2?
Для двох нестискуваних шарів однакової густини з різницею швидкостей ΔV = V₁ − V₂ швидкість зростання КГ становить γ = k·ΔV/2, де k = 2π/λ — хвильове число поверхневого збурення. Це означає, що коротші довжини хвилі (більше k) зростають швидше в межі рівних густин. При різних густинах: γ = k·ΔV·√(ρ₁ρ₂)/(ρ₁+ρ₂), що зменшує зростання при великих контрастах густин.
Де нестабільність КГ плазми зустрічається в природі?
Нестабільність КГ плазми поширена в космічних плазмах. На магнітопаузі Землі сонячний вітер, проходячи повз магнітосферу, породжує КГ-вихори, які переносять сонячну плазму всередину. На хмарних поясах Юпітера видно КГ-хвилі вздовж меж струменевих течій. У сонячній короні безпосередньо спостерігали КГ-вихори у пір'ях і коронному дощі. Вітрові туманності пульсарів та астрофізичні струмені також демонструють характерні КГ-структури.
Як відношення густин впливає на нестабільність КГ?
Відношення густин ρ₂/ρ₁ змінює швидкість зростання через множник √(ρ₁ρ₂)/(ρ₁+ρ₂). При рівних густинах (відношення = 1) множник дорівнює 1/2, забезпечуючи максимальну швидкість зростання k·ΔV/2. Зі збільшенням контрасту густин ефективна швидкість зростання зменшується, і важкий шар чинить опір відхиленню. Дуже великий контраст густин пригнічує нестабільність.
Що таке вихрове згортання і як воно розвивається?
Вихрове згортання — це нелінійне насичення нестабільності КГ. Спочатку синусоїдальні збурення зростають експоненціально (лінійна фаза). Коли амплітуди стають порівнянними з довжиною хвилі, завихреність концентрується в дискретних ядрах: межа згортається у вихори типу «котяче oko». Сусідні вихори врешті-решт зливаються (злиття вихорів), утворюючи більші структури — інверсний каскад енергії, що забезпечує турбулентне перемішування.
Як мода збурення k впливає на симуляцію?
Ціла мода k задає, скільки повних синусоїдальних довжин хвилі вміщується у домені симуляції. Мода k=1 показує одне велике вихрове згортання; k=4 — чотири одночасні вихори. Вищі моди зростають швидше (γ ∝ k) у лінійній фазі, але й насичуються при меншій амплітуді відносно розміру домену. Найбільш видовищне згортання зазвичай відбувається при середніх модах k=2–4.
Що стабілізує нестабільність КГ у намагніченій плазмі?
Компонента магнітного поля, паралельна течії, забезпечує магнітний натяг, що чинить опір вигину ліній поля під час утворення вихорів. Нестабільність стабілізується, коли альфвенівська швидкість V_A = B/√(μ₀ρ) перевищує ΔV/2, тобто коли магнітна енергія перевищує кінетичну енергію зсуву. Перпендикулярне магнітне поле не має стабілізуючого ефекту.
Яка швидкість зростання шару перемішування у нелінійній фазі?
Після початкового експоненціального зростання ширина шару перемішування δ зростає лінійно з часом: δ = α·ΔV·t, де α ≈ 0,06–0,1 — коефіцієнт турбулентного перемішування (константа Брауна-Рошко). Ця лінійна фаза зростання зберігається до заповнення шаром зсуву всього домену. Товщина імпульсу шару перемішування θ зростає як dθ/dt ≈ 0,017·ΔV.
Як нестабільність КГ плазми пов'язана з термоядерною енергетикою?
У крайовій плазмі токамаків нестабільність КГ спричиняє турбулентний перенос через сепаратрису, змішуючи гарячу плазму ядра з холоднішою крайовою плазмою. В інерційному термоядерному синтезі (ІТС) нестабільності КГ на межі паливо-аблятор погіршують симетрію імплозії та змішують холодний матеріал аблятора з гарячим паливом, гасячи запалення. Контроль нестабільності КГ є ключовим для досягнення термоядерної енергії.
Який чисельний метод використовує ця симуляція?
Симуляція використовує 2D метод вихрового листа на переодичному домені. Межа між двома шарами плазми представлена N маркерними точками, кожна з яких несе завихреність, пропорційну місцевому стрибку швидкості. Їхні позиції розвиваються за дискретизованим ядром Біо-Савара (десингуляризація вихрових сгустків). Інтегратор Рунге-Кутта 4-го порядку просуває позиції листа, захоплюючи згортання аж до фази насичення.