Розпад Торріцелових Воротів Taylor–Green

Розв’язок вічної задачі

Рівняння Навіє-Стокса, що описують рух рідин, майже ніколи не мають чистого аналітичного розв’язку, коли присутні в’язкість та нелінійність – саме тому турбульозність залишається однією з найбільших невирішених проблем класичної фізики. Вісце-віяльний віж, опублікований Дж. І. Тейлором і А. Е. Гріном у 1937 році, є рідкісною та цінною винятковою ситуацією: повністю аналітичне розв’язання рівнянь Навіє-Стокса для незмінного (інкомпресоване) потоку в двох вимірах, яке описує масив протистоячих віж, що плавно розпадаються під дією лише в’язкості.

жива демонстрація · пов'язана симуляція● LIVE

2D розв’язник

Класичне двовимірне потоки Taylor-Green починається з простого синусоїдального поля швидкості та експоненційно зменшується в часі, точна залежність від часу випливає з підстановки ансац у рівняння Нав’є-Стоксів і підтверджує її тотожність, включаючи вічну терм:

u(x,y,t) = cos(x) sin(y) * exp(-2*nu*t)

v(x,y,t) = -sin(x) cos(y) * exp(-2*nu*t)

p(x,y,t) = -(1/4)(cos(2x) + cos(2y)) * exp(-4*nu*t)

// nu = кінематична в’язкість; це точно задовольняє рівняння Нав’є-Стоксів для незмірюваного потоку в кожний момент часу t

Все зменшується з простим експоненційним законом: амплітуда швидкості падає як exp(-2*nu*t), а кінетична енергія, пропорційна квадрату швидкості, падає як exp(-4*nu*t). Немає турбулентності, не відбувається каскаду - лише впорядкована марина віртів, що розсіюються за допомогою в’язкості, безперервно зберігаючи однакову просторову структуру. Саме ця простота робить її цінною: будь-який чисельний розв'язок може бути спрямований на початкові умови потоку Taylor-Green, і його вихід порівнюється, пункт за пунктом та цифра за цифрою, з відомим аналітичним рішенням.

u(x,y,t) =  cos(x) sin(y) * exp(-2*nu*t)
v(x,y,t) = -sin(x) cos(y) * exp(-2*nu*t)
p(x,y,t) = -(1/4)(cos(2x) + cos(2y)) * exp(-4*nu*t)
// nu = kinematic viscosity; this satisfies incompressible Navier-Stokes EXACTLY at every t

Версія 3D: бенчмарк турбулентності

Тривимірна торнадо Тейлора-Гріна використовує подібний синусоїдальний початковий стан, але без аналітичного розв’язку – і ця відсутність точного розв’язку є причиною, чому випадок 3D став відомим. Починаючи з гладких великих масштабів синусоїдальних торнадо, нелінійний член потоку розтягує та нахиляє лінії торнадо, руйнуючи великі впорядковані структури на менші через розтягнення торнадо – фундаментальний механізм, за допомогою якого тривимірна турбулентність передає енергію від великих масштабів до малих (механізм, який просто не існує у двох вимірах, тому що випадок 2D залишається ламинарним і впорядкованим).

Два діагностування відстежують цей перехід. Кінетична енергія монотонно та плавно зменшується протягом усього процесу, розсіюючись зрештою через в’язкість. Ентропія – об’ємний інтеграл кутової швидкості в квадраті, пряме вимірювання того, наскільки багато обертальної, мікромасштабної структури присутня – робить протилежне: вона зростає, коли потік розвиває все більш дрібні ворсинки, досягає піку близько безрозмірної часу приблизно 9 (для стандартного початкового числа Рейнольдса, що використовується у випробувальних дослідженнях), і лише тоді починає падати через в’язке розсіювання, яке масштабується з ентропією, зрештою перевершуючи виробництво нових малих масштабів каскадом.

Чому кожен CFD код тестується на ньому

3D Taylor-Green vortex став де-факто стандартом першого бенчмарка для будь-якого нового обчислювального потоку динаміки (CFD) розв’язання, яке стверджує, що воно може обробляти перехідні або турбулентні потоки, оскільки воно пропонує рідкісну річ: повністю специфічну, відтворювану геометрію без геометрії (періодний куб, без стін, без умовних меж входу/виходу, щоб помилитися), яка все ж таки розвиває повну фізику розтягування воронки, перехід до турбулентності та енергетичний каскад вниз до найменших розділених масштабів, які вирішуються в прямій чисельній симуляції (DNS). Оскільки налаштування таке чисте, результати від різних кодів і різних числових методів - кінцевий об’єм, спектральний, латте-болотний – можуть бути порівняні один з одним та з високороздільними посиланнями DNS майже пункт за пунктом, ізолюючи чисельні помилки методу від будь-якої невизначеності в граничних умовах або геометрії.

Як саме час пікової антропії вимірює

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

Frequently asked questions

Чи є віялок Тейлора-Гріна реальним фізичним потоком, який можна було б побачити в природі?

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

Чому 2D віялок Тейлора-Гріна має точне розв’язання, а 3D – ні?

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

Що показує пікове значення енстропії для інженера CFD?

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

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Taylor-Green Vortex і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.

▶ Відкрити симуляцію Taylor-Green Vortex

Розпад Торріцелових Воротів Taylor–Green: Улюблений Бенчмарк Турбулентності