Учебная программа · вычислительная гидродинамика
Расчёт течения в круглой трубе
Ламинарное и турбулентное течение в трубе круглого сечения. Расчёт ведётся в переменных U (радиальная), V (азимутальная), W (осевая) и P (давление), методом конечных объёмов.
Параметры расчёта
—
На этой геометрии проверены на сервере эти два режима.
—
Картины течения и коэффициенты
| Величина | Из расчёта | Теория | Расхождение |
|---|
Как читать результат. Коэффициент Дарси берётся по развитому участку трубы
(последние 40 %): на входе течение ещё формируется, и градиент давления падает.
Дискретный максимум скорости достигается не на оси, а в первой ячейке, поэтому
u_max/u_сред ≈ 2(1 − 0.25/n_z²) — при 12 ячейках по радиусу это около 1.997.
Теория и история задачи
Раздел загружается…
Расчёты на больших компьютерах
Раздел загружается…
Инструкция по пользованию и презентация
Инструкция по пользованию
- Задайте число Рейнольдса. Ползунком или числом. Подсказка под ним покажет режим течения и опорные значения коэффициента трения.
- Выберите модель турбулентности. Для Re < 2300 — ламинарный режим, иначе realizable k-эпсилон. ν подставляется автоматически как D / Re (скорость на входе в модели равна 1 м/с, D = 1 м).
- Задайте сетку. Осевое число — длина трубы (100 ячеек на 100 диаметров), радиальное — разрешение у стенки. Осевое направление меняется быстрее всего.
- Нажмите «Рассчитать». Сервер считает 15–60 секунд; вместо крутилки показывается понятный текст. Пока расчёт идёт, кнопка заблокирована.
- Смотрите результат. Цветное сечение, профиль скорости с наложенной параболой Пуазейля, давление вдоль трубы и точка расчёта на кривых Пуазейля и Блазиуса.
- Если сервер недоступен — переключитесь на демо-режим: он считает аналитическое решение мгновенно и работает без сети.
Демо-режим — это аналитика, а не CFD. Он подставляет готовое решение уравнения Пуазейля или степенной закон 1/7 и показывает ту же картину. Это удобно, чтобы посмотреть на форму течения, но расхождение с реальным расчётом на нём измерить нельзя.
Презентация программы — сценарий показа
- Постановка (2 мин). Открыть «Теория», раздел «Постановка задачи». Показать уравнения импульса в цилиндрических координатах и объяснить, что при U = V = 0 они сводятся к одному уравнению Лапласа — именно это и решает программа.
- Историческая справка (1 мин). Раздел «История»: Пуазейль 1840, Рейнольдс 1883, Блазиус 1913. Главный вывод: для ламинарного режима есть точное аналитическое решение, для турбулентного — нет, и это породило проблему замыкания.
- Ламинарный контрольный пример (3 мин). Задать Re = 1000, выбрать ламинарный режим, рассчитать. Главный вывод: коэффициент трения из расчёта совпадает с формулой 64/Re в пределах 1 %. Программа здесь точная.
- Турбулентный режим (3 мин). Re = 100000, realizable k-эпсилон. Профиль становится плоским в ядре, коэффициент приближается к закону Блазиуса с расхождением в несколько процентов. Объяснить, откуда берётся остаток: сетка грубая у стенки и используются пристенные функции.
- Чувствительность к сетке (2 мин). Увеличить радиальное число ячеек и показать, что результат приближается к теории. Это честная демонстрация того, как сетка влияет на ответ.
- Модели турбулентности (2 мин). Раздел «Турбулентность»: объяснить, почему вязкость турбулентного течения эффективная и откуда берётся потребность в замыкании.
- ИИ в моделях (2 мин). Раздел «ИИ»: снижение сопротивления на 43 % при Re_т = 180 и честные ограничения методов.
- Масштабирование (2 мин). Раздел про большие компьютеры: почему наш расчёт на 7200 ячейках занимает секунды, а DNS требует миллиардов ячеек, и куда такое отправляют.
Частые вопросы и типичные ошибки
- Коэффициент больше теоретического на 20–100 %. Почти всегда труба слишком короткая или сетка слишком грубая по радиусу.
- u_max/u_сред ≈ 2. Это признак ламинарного режима. В турбулентном отношение около 1.2.
- Расчёт не сходится. Уменьшите число ячеек или увеличьте число итераций. Сервер один и очень загруженный — параллельный запуск не помогает.
- Профиль скорости «плоский» у стенки. Первые ячейки у стенки разрешены грубо — это работа пристенных функций.