🌊 Розпад Вихору Тейлора-Ґріна

Симулюйте вихор Тейлора-Ґріна — еталон турбулентності. Спостерігайте, як початкове синусоїдальне поле переносить енергію до менших масштабів, а енстрофія зростає.

ГідродинамікаІнтерактивний
Колір = завихреність · P пауза · R скидання

Як це працює

Вихор Тейлора-Ґріна починається з periodичної синусоїдальної початкової умови у квадраті [0, 2π]². Поле швидкостей u(x,y,0) = U·sin(kx)cos(ky) і v(x,y,0) = −U·cos(kx)sin(ky) точно задовольняє умову нестисливості. Ми вирішуємо формулювання завихреність-функція течії 2D рівнянь Нав'є-Стокса на сітці методом скінчених різниць.

Завихреність ω еволюціонує як Dω/Dt = ν∇²ω. Енергія, спочатку сконцентрована на хвильовому числі k, каскадує до вищих хвильових чисел у міру того, як вихрові структури розтягуються і складаються, тоді як в'язкість дисипує енергію на найменших масштабах.

u(x,y,0) = U sin(kx) cos(ky)
v(x,y,0) = −U cos(kx) sin(ky)
ω = ∂v/∂x − ∂u/∂y
∂ω/∂t + u·∇ω = ν∇²ω

Часті запитання

Що таке вихор Тейлора-Ґріна?

Вихор Тейлора-Ґріна — точний розв'язок нестисливих рівнянь Нав'є-Стокса при t=0, заданий синусоїдальним полем швидкостей. Це канонічний тестовий приклад для моделей турбулентності.

Що таке каскад енергії у турбулентності?

Каскад енергії — процес передачі кінетичної енергії від великих вихорів до дедалі менших, поки в'язкісна дисипація не перетворить її на теплоту. Теорія Колмогорова описує спектр E(k) ~ k-5/3.

Що таке енстрофія?

Енстрофія — об'ємний інтеграл квадрата завихреності: E_n = 0.5·∫ω²dV. У 2D-турбулентності вона зберігається; у 3D — зростає під час каскаду і досягає максимуму перед дисипацією.

Що таке рівняння Нав'є-Стокса?

Рівняння Нав'є-Стокса описують рух в'язкої рідини: ρ(Du/Dt) = −∇p + μ∇²u + f. Вони містять закон збереження імпульсу та, разом з рівнянням нерозривності, закон збереження маси.

Чому вихор Тейлора-Ґріна є еталонним тестом?

Аналітичні початкові умови, periodичні граничні умови та відомі темпи затухання роблять його ідеальним тестом для DNS-кодів, LES-моделей і спектральних методів.

Що таке число Рейнольдса і як воно впливає на симуляцію?

Число Рейнольдса Re = UL/ν вимірює відношення інерційних сил до в'язкісних. При більших Re каскад енергії охоплює менші масштаби, а енстрофія досягає максимуму пізніше.

Що таке пряме числове моделювання (DNS)?

DNS вирішує всі масштаби турбулентності від великих вихорів до масштабу Колмогорова без моделювання. Для вихору Тейлора-Ґріна при Re=1600 потрібно ~512³ вузлів сітки.

Що таке масштаб Колмогорова?

Масштаб Колмогорова η = (ν³/ε)1/4 — найменший розмір, де кінетична енергія дисипується в'язкістю. Нижче цього масштабу течія гладка і ламінарна.

Як розтяжка вихорів впливає на турбулентність?

Розтяжка вихорів виникає, коли вихрова нитка вирівнюється з власним вектором тензора швидкостей деформацій. Розтяжка підсилює завихреність і переносить енергію на менші масштаби.

Що відбувається з кінетичною енергією з часом?

Кінетична енергія монотонно спадає від початкового значення. Енстрофія спочатку зростає в міру формування і розтяжки вихрових структур, потім спадає разом із дисипацією.

Про цю симуляцію

Симуляція розв'язує рівняння переносу завихреності у 2D на періодичній сітці 64×64, стартуючи з точного початкового поля Тейлора-Гріна u=U·sin(kx)cos(ky), v=−U·cos(kx)sin(ky). Кожен кадр оновлює завихреність через скінченнорізницеву адвекцію та в'язку дифузію (ν=1/Re), після чого перераховує поле швидкості. Спостерігайте, як кінетична енергія монотонно спадає, а ензострофія спершу зростає, поки вихрові структури розтягуються, а потім падає, коли дисипація бере гору.

🔬 Що показано

Масив 2D вихрових комірок, що чергують напрямок обертання, розтягуються, згортаються та згасають з часом, ілюструючи турбулентний каскад енергії від великих до дрібних масштабів у двічі-періодичному боксі.

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

Перетягуйте Число Рейнольдса Re, щоб змінити співвідношення інерції та в'язкості, Хвильове число k — щоб задати початкову кількість вихорів, і Амплітуду U — щоб масштабувати стартову швидкість. Перемикайте випадаючий список Display між Vorticity, Speed і Streamlines, або натискайте P/R для паузи/скидання.

💡 Чи знали ви?

Завдяки відомому аналітичному початковому умовному полю та відсутності стінок, вихор Тейлора-Гріна — один із найпопулярніших бенчмарків для тестування DNS та LES кодів турбулентності на роздільностях до 512³.

Часті запитання

Чому ензострофія спершу зростає, а потім падає?

Розтягування вихорів посилює градієнти завихреності на початковому етапі, штовхаючи енергію до дрібніших масштабів і піднімаючи ензострофію, поки в'язка дифузія (масштабована як 1/Re тут) не почне домінувати і не розсіє її.

Що робить збільшення Re із симуляцією?

Вище число Рейнольдса зменшує в'язкість nu=1/Re, дозволяючи каскаду енергії досягати тонших структур перш ніж дисипація візьме гору, тож завихрені візерунки залишаються чіткими довше.

Чому сітка періодична, а не обмежена стінками?

Періодичні межі (реалізовані тут через модульну індексацію у функції idx()) прибирають крайові ефекти, відповідаючи класичній постановці Тейлора-Гріна і дозволяючи вихорам безшовно «загортатись» через домен.

Яка різниця між трьома режимами відображення?

Vorticity фарбує кожну комірку за напрямком та величиною обертання, Speed показує локальну величину швидкості у відтінках блакитного, а Streamlines малює траєкторії часток крізь миттєве поле швидкості.

Як обчислюється кінетична енергія на екрані?

Симуляція підсумовує u²+v² по кожній комірці сітки та нормалізує на 2N², даючи поточну статистику «Kinetic Energy» поряд із ензострофією та швидкістю дисипації.