Добрый день.
Моделирую центробежный (радиальный) вентилятор, Модель состоит из трех доменов: входной участок (статор), рабочее колесо (ротор), улитка (статор).
Решаю задачу в стационарной постановке. Связываю домены интерфейсом Stage (mixing-plane). На входе в вентилятор задаю постоянное полное давление равное 0 (Pt=0), на выходе ГУ - outlet, с постоянным массовым расходом. Опорное давление Pref=1 атм. Модель турбулентности k-e.
Настройки решателя приложены к теме. Решатель запуска по схеме двойной точности.
Для построения характеристики вентилятора параметризую массовый расход воздуха на выходе из вентилятора. Расчет проходит без ошибок для расходов от 0,6 до 0,84 кг/с, невязки опускаются до -4 степени. Для малых расходов от 0,12 до 0,6 после нескольких итераций выскакивает ошибка. Текст ошибки приложен к теме.
Также ошибка выскакивает в окне workbench, также прикладываю ее к теме.
Подскажите с чем связана ошибка, как исправить ее?
ГУ outlet пробовал изменить на Opening для лучшей сходимости, но ошибка также выскакивает.
Вложение | Размер |
---|---|
solver_control.png | 84.24 КБ |
display_monitor.png | 106.64 КБ |
error_cfx_2.png | 9.93 КБ |
error_workbench.png | 8.44 КБ |
Здравствуйте!
Скорее всего, вы выходите за границы газодинамической устойчивости. Попробуйте уменьшить timescale factor, если не поможет - пробуйте нестационар. Как вариант можно попробовать плавно наращивать обороты.
С пониженным timescale factor решение сходится очень долго, вариант с плавным наростанием оборотов интересно попробовать. Возможно написать выражение меняющее количество оборотов в зависимости от шагов итерации? Например каждый 200 итераций менять количество оборотов в минуту с 1000 до 4500 с шагом 500 оборотов в мин.
Написал выражение для изменения количества оборотов ротора во время расчета, буду пробовать:
if (citern >= 800, 4500, if (citern >= 600, 3500, if (citern >= 400, 2500, if (citern >= 200, 1500, 1000)))) [rev min^-1]
или так
if (aitern > 800, 4500, (4500/800*aitern))[rev min^-1]
Добавить комментарий