CFD-моделирование насосов позволяет инженерам анализировать гидродинамические процессы внутри проточной части без создания физических прототипов. Уравнения Навье-Стокса, решённые численно, дают распределение давления, скоростей и турбулентности с детализацией, недоступной упрощённым гидравлическим расчётам. Показательный пример: при неудачной форме лопаток потеря напора достигает 15 % и напрямую бьёт по КПД.
Коротко главное
- Для стационарных расчётов характеристик используйте MRF, для анализа пульсаций, Sliding Mesh.
- k-omega SST предпочтительнее k-epsilon из-за лучшей работы с отрывными течениями и пристеночными слоями.
- Критерий сходимости: невязки упали на 3–4 порядка, интегральные параметры (напор — момент) стабилизировались.
- Требование к сетке: y+ ≈ 1 для моделей с интегрированием пристеночной области, y+ > 30 для wall functions.
- Симметричные колеса считают в режиме single passage, экономия до 90 % ресурсов.
Фундаментальные принципы и уравнения
Расчёт начинается с дискретизации области и численного решения уравнений сохранения массы, импульса и энергии. Для несжимаемой жидкости, воды или масел, система упрощается:
| Уравнение | Физический смысл | Математическая форма |
|---|---|---|
| Уравнение неразрывности | Сохранение массы | \nabla \cdot \vec{v} = 0 |
| Уравнение Навье-Стокса | Сохранение импульса | \rho \left( \frac{\partial \vec{v}}{\partial t} + \vec{v} \cdot \nabla \vec{v} \right) = -\nabla p + \nabla \cdot \bar{\bar{\tau}} + \rho \vec{g} |
Здесь \rho, плотность, \vec{v}, вектор скорости, p, давление, \bar{\bar{\tau}}, тензор вязких напряжений, \vec{g}, ускорение свободного падения. Турбулентность замыкают дополнительными уравнениями. Стандартная пара, k (кинетическая энергия турбулентности) и \epsilon (скорость её диссипации).
Подготовка геометрии и расчётной сетки
Точность расчёта упирается в качество сетки. Геометрию насоса, рабочее колесо, отвод, входной патрубок, перед разбивкой очищают от мелочи, на гидродинамику не влияющей: резьбы, логотипы, мелкие фаски.
Типы ячеек и гибридные сетки
Для насосов в CFD обычно применяют гибридные сетки:
- Гексаэдры (Hex) — основная область потока. Меньшая численная диссипация, выше точность.
- Призмы (Prisms) — пристеночные слои. Разрешают градиенты скорости в пограничном слое при меньшем числе ячеек, чем тетраэдры.
- Тетраэдры / тетра-призмы (Tet/Prism) — сложные зоны: углы отвода, переходы между участками.
Требования к y+ и пристеночным функциям
Параметр y+, безразмерное расстояние от стенки до центра первой ячейки, определяет, как разрешается пристеночный слой. Выбор зависит от модели турбулентности:
Low-Re / Интеграция до стенки
Целевой y+: 1 ± 0.5
Обязательно для k-\omega SST и LES. Сетка должна разрешать вязкий подслой. Ячеек у стенки много, но без этого модель не сойдётся.
Wall Functions
Целевой y+: 30–300
Применяется в k-\epsilon со стандартными пристеночными функциями. Первая ячейка попадает в логарифмическую область. Ячеек меньше, но метод плохо работает при отрыве потока.
Single Passage (Расчёт одного канала)
Если рабочее колесо симметричное (Z лопастей), можно считать только один канал между лопатками с периодическими граничными условиями. Количество ячеек в зоне вращения падает в Z раз, экономия до 90 % вычислительных ресурсов. Метод не учитывает ротор-статорное взаимодействие, неравномерность потока от предыдущих лопаток и асимметрию отвода.
Моделирование вращения и турбулентности
Главная сложность расчёта насоса, взаимодействие вращающегося рабочего колеса с неподвижными элементами: отводом, корпусом. Выбор метода зависит от требуемой точности и доступных ресурсов.
Методы взаимодействия зон
| Метод | Как работает | Когда применять |
|---|---|---|
| MRF (Multiple Reference Frame) | Стационарный подход. Вращающаяся зона считается стационарной, но в уравнения входят центробежные и кориолисовы силы. Интерфейс между зонами, «Frozen Rotor» (фиксированное взаимное расположение). | Быстрый расчёт Q-H, первичная оптимизация. Пульсации и ротор-статорное взаимодействие не учитывает. |
| Sliding Mesh | Нестационарный. Сетка вращающейся зоны поворачивается относительно неподвижной. Данные передаются через скользящую границу. | Точные расчёты пульсаций давления, акустики, кавитационных флуктуаций. Требует малого шага по времени. |
| Transient Rotor-Stator | Передача данных между сетками через преобразование Фурье. Позволяет использовать крупные шаги по времени, сохраняя точность взаимодействия. | Компромисс между MRF и Sliding Mesh, когда нужны нестационарные эффекты. |
Модели турбулентности
В инженерной практике чаще всего используют RANS. Конкретная модель зависит от характера течения:
| Модель | Особенности | Для насосов |
|---|---|---|
| k-\epsilon (Realizable) | Устойчивая, быстрая. Хромает в зонах отрыва и при больших градиентах давления. | Грубые оценки на больших расходах, где поток пристеночный и не оторван. |
| k-\omega SST | Лучше описывает пограничный слой и отрыв. Объединяет k-\epsilon в ядре потока и k-\omega у стенки. | Стандартный выбор для оптимизации рабочих колёс. Подходит для большинства гидродинамических задач. |
| LES / DES | Крупные вихри разрешаются, мелкие моделируются. Точность высокая, цена, вычислительная. | Научные исследования, акустика, сложная кавитация. В серийном проектировании, редко. |
Граничные условия и настройка расчёта
Чтобы расчёт сошёлся, граничные условия задают так:
- Входной патрубок: Mass Flow Inlet или Velocity Inlet. Лучше задать профиль скорости и турбулентные параметры (интенсивность ~5 %, гидравлический диаметр), иначе поток искусственно сгладится.
- Выходной патрубок: Pressure Outlet. Статическое давление, обычно 0 Pa gauge. Выведите выход за пределы отвода, чтобы обратная волна давления не влияла на область расчёта.
- Стенки: No-slip. Для вращающихся поверхностей, заданная угловая скорость \omega.
- Интерфейсы: Interface между вращающейся и неподвижной зонами. Тип зависит от метода (MRF, Sliding Mesh).
Расчёт ключевых гидродинамических параметров
По завершении расчёта инженеры извлекают интегральные параметры, по которым оценивают работу насоса.
Гидравлический КПД
Отношение полезной мощности к затраченной на валу:
Q, расход, H, напор, T, крутящий момент на рабочем колесе, \omega, угловая скорость. Момент берётся только с колеса: механические потери в подшипниках и уплотнениях в формулу не входят.
Кавитационный запас (NPSHr)
Требуемый кавитационный запас, минимальное избыточное давление на входе, при котором кавитация ещё не начинается. В CFD его находят снижением давления на входе до момента, когда напор начинает падать на 3 %:
p_{inlet, abs}, абсолютное давление на входе, p_v, давление насыщенных паров жидкости. Зоны кавитации визуализируют через многофазные модели, например, модель Зайцева или гомогенную двухфазную, где объёмная доля пара считается по уравнению Рэлея–Плассеца.
Валидация и верификация модели
Результаты CFD не имеют ценности без проверки. Процесс состоит из верификации (правильность решения уравнений) и валидации (соответствие физическому миру).
Верификация сеточной независимости
Делают серию расчётов на сетках разного разрешения, 1 млн, 3 млн, 6 млн ячеек. Если изменение КПД или напора между последней и предпоследней сеткой укладывается в 1–2 %, решение считают независимым от сетки.
Сравнение с экспериментом
При сопоставлении CFD со стендовыми данными учитывают:
- Погрешность датчиков: Расходомеры и манометры дают погрешность 0,5–1 %.
- Условия тестирования: Температура жидкости влияет на вязкость и плотность; шероховатость стенок в CFD часто идеализируется как гладкая.
- Потери в присоединительных трубопроводах: CFD считает только проточную часть насоса. Потери во входном и выходном патрубках стенда нужно либо вычитать, либо моделировать удлинённые участки.
- Стандарты: Ориентир, ISO 9906. Погрешность ±2–3 % по напору соответствует Grade 2B и считается хорошим результатом для CFD.
Обзор программных решений
Рынок CFD-софта для насосостроения делится на универсальные коммерческие пакеты и open-source решения.
Commercial (ANSYS, Siemens)
ANSYS CFX/Fluent, STAR-CCM+. Пакеты с модулями для турбомашин (TurboSystem, Turbomachinery). Автоматизация сеткогенерации, поддержка многофазных течений.
Open Source (OpenFOAM)
OpenFOAM. Бесплатный, открытый код. Требует знаний C++ и ручной настройки. Популярен в академической среде и для кастомизации физических моделей. Сообщество активно развивает насосные модули.
Выбор зависит от бюджета, компетенций и сложности задачи. Для типовой оптимизации КПД коммерческие пакеты с визуальными интерфейсами и готовыми шаблонами экономят 4–6 часов на настройку расчёта.
