Завихреність рідини та нестійкість Кельвіна–Гельмгольца

Вихровий шар Біо–Савар · 2D ідеальна рідина · Скручування вихорів · Зсувна нестійкість

t = 0.000 | крок = 0 | N = 32 вихори | Γ = ±0.006 | ε = 0.008 | dt = 0.0050

🌊 Нестійкість Кельвіна–Гельмгольца

Якщо два шари рідини рухаються з різними швидкостями, будь-яке збурення на їхній межі зростає експоненційно — це нестійкість Кельвіна–Гольмгольца (НКГ). Вона є одним із найважливіших механізмів виникнення турбулентності в природі: хмарні хвилі-«білі гребені», внутрішні хвилі океану, Велика Червона Пляма Юпітера, сонячні корональні джети та магнітопауза Землі — всюди можна побачити характерні КГ-рулони.

🔬 Фізика

Неперервний вихровий шар (межа двох рівномірних потоків) апроксимується N дискретними точковими вихорами з циркуляцією Γ. Кожен вихор переміщується зі швидкістю, яку індукують усі інші, за законом Біо–Савара (регуляризованим скінченним ядром ε). Малі синусоїдальні збурення зростають зі швидкістю σ = |Γ|k/4π (k — хвильове число) — скручування Кейдена–Бірхгофа. Врешті-решт виникають окремі спіральні вихрові структури та злиття пар вихорів.

🎮 Як користуватись

Натисніть Пауза, щоб зупинити. Скинути — перезапустити з поточними параметрами. Вихори N — задати кількість дискретних вихорів. Збільшіть Збурення для швидшого скручування. Два шари додає контр-обертальний ряд на 0,08 вище — це дає вихрову вулицю Кармана. Випадкові — хаотичні початкові умови.

💡 Приклади в природі

НКГ формує характерні хмари-«хвилі» (хмари Кельвіна–Гельмгольца), океанічні внутрішні хвилі (видно з супутника), перемішувальні шари в зоряних надрах, плазмові нестійкості на магнітопаузі Землі, пальці Релея–Тейлора в залишках наднових і зсувні нестійкості в акреційних дисках чорних дір.

Про нестійкість Кельвіна–Гельмгольца

Ця симуляція моделює двовимірний вихровий шар без в'язкості на межі двох потоків рідини, що рухаються з різними швидкостями. Неперервний шар дискретизується у N точкових вихорів із циркуляцією Γ, кожен із яких переміщується під дією швидкісного поля, яке індукують усі інші вихори через регуляризований закон Біо–Савара: vx = −(Γ/2π)·dy/r², vy = +(Γ/2π)·dx/r², де скінченне ядро ε запобігає сингулярностям. Periodicity у напрямку x забезпечується зображеними стовпцями.

Повзунки задають кількість вихорів N (8–64), циркуляцію Γ кожного вихору, радіус ядра ε, явний крок Ейлера dt та амплітуду початкового синусоїдального збурення. Більше збурення прискорює скручування. Кнопки «Один шар», «Два шари» і «Випадкові» змінюють початкову конфігурацію. Ця зсувна нестійкість зумовлює турбулентність у хмарних хвилях, океанічних течіях, атмосфері Юпітера і сонячному вітрі — і є фундаментальною для гідродинаміки.

Поширені запитання

Що таке нестійкість Кельвіна–Гельмгольца?

Це нестійкість, що виникає, коли два шари рідини ковзають один відносно одного з різними швидкостями. Будь-яке мале збурення на спільній межі зростає експоненційно, закручуючи її у ряд спіральних вихорів. Це один із найпоширеніших механізмів переходу до турбулентності в природі.

Як симуляція обчислює рух?

Вихровий шар представлено N дискретними точковими вихорами. На кожному кроці кожен вихор переноситься сумарною швидкістю, яку індукують усі інші через закон Біо–Савара без в'язкості, і інтегрується явною схемою Ейлера з кроком dt.

Навіщо потрібен радіус ядра ε?

Чисто точковий вихор індукує швидкість, що прямує до нескінченності як 1/r при зближенні двох вихорів, — це спричиняє числові вибухи. Додавання ε² до знаменника регуляризує ядро, надаючи кожному вихору плавне скінченне ядро, що зберігає стабільність симуляції, майже не впливаючи на великомасштабне скручування.

Що роблять повзунки «Вихори N» і «Циркуляція Γ»?

N (8–64) задає, наскільки детально неперервний вихровий шар апроксимується дискретними вихорами; більше вихорів — чіткіші спіральні ядра. Γ — циркуляція кожного вихору, що масштабує силу індукованого швидкісного поля і, відповідно, швидкість розвитку нестійкості.

Яка швидкість зростання нестійкості?

Для вихрового шару збурення з хвильовим числом k зростає експоненційно зі швидкістю σ = |Γ|k/4π в ідеальному наближенні. Коротші хвилі зростають найшвидше, тому дискретизований шар швидко утворює дрібномасштабні рулони, які згодом зливаються у більші структури через спарювання вихорів.

Що показує режим «Два шари»?

«Два шари» додає другий ряд вихорів із протилежним обертанням приблизно на 0,08 одиниці вище першого. Протилежні циркуляції утворюють конфігурацію, подібну до вулиці вихорів Кармана за тупим тілом, і ілюструють взаємодію двох зсувних шарів.

Чи є симуляція фізично точною?

Це точна двовимірна модель точкових вихорів без в'язкості, що вірно відтворює раннє скручування і спарювання вихорів, передбачені класичною теорією. Однак вона не враховує в'язкість, тривимірні ефекти та стисливість, тому пізні хаотичні стадії мають якісний, а не кількісний характер.

Чому течія з часом стає хаотичною?

Системи точкових вихорів із більш ніж трьома вихорами загалом неінтегровні і демонструють чутливу залежність від початкових умов. У міру скручування і спарювання вихорів найменші відмінності посилюються, породжуючи хаотичний стан інтенсивного перемішування, характерний для розвиненої зсувної турбулентності.

Чому вихори пофарбовані по-різному?

Сині вихори мають додатну циркуляцію (Γ > 0), червоні — від'ємну (Γ < 0), що вказує на протилежні напрямки обертання. В режимі «Один шар» усі вихори мають один знак, тоді як «Два шари» і «Випадкові» змішують обидва знаки — це видно по слідах, які вони залишають.

Де нестійкість Кельвіна–Гельмгольца зустрічається в реальному світі?

Вона породжує характерні хмари-хвилі Кельвіна–Гельмгольца в небі, океанічні внутрішні хвилі, перемішувальні шари у зорях, плазмові нестійкості на магнітопаузі Землі та смугасті зсувні потоки Юпітера. Та сама математика описує зрив вихорів за перешкодами і перемішування в струменях.