CFD · Нестисливий потік · Числовий аналіз
📅 Липень 2026 ⏱ ≈ 10 хв читання 🎯 Просунутий

Паразитні моди у нестисливих рівняннях Нав'є-Стокса

Нестисливі рівняння Нав'є-Стокса пов'язують швидкість і тиск через обмеження — ∇·u = 0 (дивергенція поля швидкості дорівнює нулю всюди, тобто об'єм рідини ніде не створюється й не зникає) — а не через рівняння еволюції тиску. Дискретизуйте швидкість і тиск у тих самих вузлах сітки без обережності — і дискретна система допускає нефізичні паттерни тиску, які точно задовольняють дискретне обмеження дивергенції, залишаючись цілком неправильними: шахові коливання, що безконтрольно ростуть і ніколи не з'являються у справжньому неперервному розв'язку. Розуміння того, чому вони виникають — і як розрізнені сітки та стабілізовані методи їх усувають — необхідне, перш ніж довіряти будь-якому полю тиску, яке видає розв'язувач.

Коротко: Якщо швидкість і тиск зберігати в одних і тих самих вузлах сітки, розв'язувачі нестисливого потоку можуть породжувати фіктивні «шахові» коливання тиску, яких рівняння не бачать і нічим не гасять. Виправлення — рознести місця зберігання тиску та швидкості (розрізнені MAC-сітки), точніше відновлювати швидкості на гранях (інтерполяція Рі-Чоу) або додати стабілізаційні члени в скінченно-елементних методах.

Автор: Команда MySimulator · Редакційна перевірка: Редакція MySimulator

1. Чому в тиску немає рівняння еволюції

Система нестисливого Нав'є-Стокса:

∂u/∂t + (u·∇)u = −∇p/ρ + ν∇²u (імпульс)
∇·u = 0 (нерозривність / нестисливість)

Ніде немає члена ∂p/∂t — тиск не є незалежним фізичним полем, що еволюціонує самостійно; він — множник Лагранжа, який забезпечує бездивергентність у кожен момент часу. Взявши дивергенцію рівняння імпульсу і використавши ∇·u = 0, отримуємо рівняння Пуассона для тиску ∇²p = −ρ∇·((u·∇)u), яке потрібно розв'язувати узгоджено з дискретним полем швидкості на кожному кроці часу — саме тут наївна дискретизація може дати збій.

2. Шахова мода тиску

На колокованій сітці — швидкість і тиск обидва зберігаються в центрах комірок — стандартна центральна різниця другого порядку для градієнта тиску у вузлі i використовує лише p_{i-1} та p_{i+1}:

(∂p/∂x)_i ≈ (p_{i+1} − p_{i-1}) / (2Δx)

Цей шаблон повністю сліпий до самого p_i. Розгляньмо поле тиску, що чергується: p_i = (−1)^i · A — «шахова» картина. Її центрально-різницевий градієнт дорівнює нулю в кожному вузлі — так само, як градієнт рівномірного (сталого) поля тиску. Дискретне рівняння імпульсу не може відрізнити справжню, фізично беззмістовну шахову моду від нешкідливого сталого зсуву — обидві невидимі для дискретного оператора градієнта, і ніщо в дискретній системі їх не гасить.

+A
−A
+A
−A
−A
+A
−A
+A
+A
−A
+A
−A

У 2D/3D той самий збій виникає одночасно вздовж обох напрямків, і нуль-простір дискретного оператора градієнта може бути великим. Похибки округлення або трохи асиметрична гранична умова достатні, щоб посіяти цю моду, і оскільки ніщо в рівняннях її не штрафує, вона може рости, доки не домінуватиме у видимому полі тиску — симптом, з яким рано чи пізно стикається кожен практик CFD.

3. Умова сумісності LBB / inf-sup

Строге формулювання того, «які пари дискретизацій швидкість/тиск безпечні», — це умова Ладиженської-Бабушки-Брецці (LBB), також звана умовою inf-sup. Для задачі сідлової точки з простором швидкостей Vh та простором тисків Qh стійкість вимагає існування β > 0, незалежного від розміру сітки h, такого, що:

infq∈Q_h supv∈V_h (∇·v, q) / (‖v‖ ‖q‖) ≥ β > 0

Інтуїтивно: для кожного нетривіального поля тиску q (включно з кандидатом-шахівницею) простір швидкостей повинен містити поле швидкості, дивергенція якого "бачить" q із силою, незалежною від сітки. Колоковані швидкість/тиск однакового порядку (елементи Q1-Q1 чи P1-P1, або наївна скінченно-різницева сітка з обома полями в тих самих точках) як правило порушують LBB — шахова мода і є саме тим полем тиску, яке дискретний оператор дивергенції не бачить за жодної швидкості.

4. Розрізнені (MAC) сітки

Сітка "маркер-і-комірка" (MAC) Гарлоу та Велча (1965) — класичне, елегантне вирішення: зберігати тиск у центрах комірок, але зберігати горизонтальну компоненту швидкості u на вертикальних гранях комірки, а вертикальну компоненту v — на горизонтальних гранях:

(∂p/∂x)_{i+1/2,j} = (p_{i+1,j} − p_{i,j}) / Δx ← тепер використовує сусідні p, жоден вузол не пропущено

(∇·u)_{i,j} = (u_{i+1/2,j} − u_{i-1/2,j})/Δx + (v_{i,j+1/2} − v_{i,j-1/2})/Δy

На розрізненій сітці оператори градієнта тиску та дивергенції — точні дискретні спряжені один одного і залучають кожне сусіднє значення тиску — шахова мода тепер дає великий, правильно спрямований градієнт і гаситься дискретним розв'язанням Пуассона. MAC-сітки задовольняють LBB за побудовою і залишаються стандартним вибором для скінченно-різницевих/ скінченно-об'ємних нестисливих розв'язувачів, включно з розв'язувачем методу проекції, використаним в інших місцях цього сайту.

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

5. Інтерполяція Рі-Чоу для колокованих сіток

Рі та Чоу (1983) показали, що колоковану сітку можна зробити стійкою без фізичного розрізнення зберігання, шляхом інтерполяції за імпульсом швидкості на грані, що використовується в рівнянні збереження маси (тиску), замість простого усереднення швидкостей сусідніх центрів комірок:

u_грань = ū_грань − D_грань · [ (∂p/∂x)_грань − (∂p/∂x)‾_грань ]

ū_грань: лінійне усереднення швидкостей центрів комірок
D_грань: коефіцієнт рівняння імпульсу (≈ Δt/ρ або 1/A_P)
(∂p/∂x)_грань: справжній градієнт через грань — чутливий до шахівниці
(∂p/∂x)‾_грань: інтерпольований усереднений градієнт — нечутливий

Коригувальний член швидкості грані явно повертає чутливість до різниці тиску через грань, яку наївне усереднення пропустило б, фактично додаючи контрольований згладжувальний член тиску четвертого порядку. Це домінантна техніка в сучасних колокованих скінченно-об'ємних кодах (OpenFOAM, більшість комерційних пакетів CFD), тому що вона зберігає всі поля в одній точці — спрощуючи неструктуровані сітки, мультисітку та обробку границь — водночас задовольняючи дискретний аналог LBB.

6. Паразитні моди у скінченних елементах — стабілізація

Те саме явище з'являється у скінченно-елементних розв'язувачах Стокса/Нав'є-Стокса. Елементи однакового порядку P1-P1 (лінійна швидкість, лінійний тиск) порушують LBB. Стандартними є два сімейства виправлень:

Стабілізаційний член PSPG (доданий до слабкого рівняння нерозривності):
+ Σ_K τ_K ∫_K ∇q · [ −∇p + ρ((u·∇)u) − ν∇²u ] dΩ

τ_K: параметр стабілізації, що залежить від розміру елемента
Зникає, коли нев'язка сильної форми дорівнює нулю → узгоджено

7. Виявлення паразитних мод на практиці

8. Практичні наслідки для розв'язувачів реального часу

Розв'язувачі рідини в реальному часі на WebGL ("Stable Fluids" Йоса Стама та його нащадки, використані в інших місцях цього сайту) зазвичай уникають проблеми шахівниці, зберігаючи швидкість на розрізненій сітці MAC-типу всередині ping-pong буферів фрагментного шейдера, або навмисно надмірно згладжуючи розв'язання тиску (кілька ітерацій Якобі замість збіжного розв'язку) — прагматичний компроміс, що ховає низькорівневі артефакти шахівниці за числовою дисипацією, прийнятний для візуальної симуляції, але не для інженерно-точного CFD.

🌊 Рідина Нав'є-Стокса (WebGL)

Подивіться, як працює розрізнений розв'язувач з проекцією тиску наживо — без шахівниці, розв'язання тиску Якобі в реальному часі.

Відкрити симуляцію →