Как повысить устойчивость работы AMG метода

[German_psv 20]

Добрый день, коллеги.

Прошу Вашей консультации по AMG-методу.

Немного о сеточной модели.

Назначение – вычисление характеристик гребного винта с насадкой в реверсе при наличии потока, имитирующего движение судна.

Создана цилиндрическая расчетная область в виде сектора 90 град (радиуса 5м и длиной 25м) с периодическими границами. В средине, соосно с расчетной областью, выделен небольшой объем (тоже сектор 90 град), в который помещена лопасть винта и насадка. В Гамбите сгенерирована тетра-сетка (первый слой 0.5мм), во Флуенте – преобразована в полиэдральную. Сетка получилась небольшая – всего 220тыс. ячеек.

Солвер – cтационарный (параллельный – 8 процессов на i7 – Vista Business),

Модель вязкости – k-e,RNG, Swirl….

Solution control: все релакс. коэффициенты уменьшил более чем вдвое; PRESTO,

PISO (только neighbor correction), в Multigrid’e включил BCGSTAB.

Задал скорость потока жидкости 1м/с,

вращение винта (MRF для объема в котором помещен винт) тоже небольшое 1рад/с.

При запуске на вычисление (бывает и в ходе вычислений) вылетает Error:

“divergence detected in AMG solver: pressure correction”.

В Multigrid’e по давлению установлен V-Cycle (defolt) – поэтому нечем уменьшить коррекцию по давлению.

В нестационарном солвере на шаге 0.01с тоже самое.

Перечитал еще раз все об AMG. Вроде все учел.

P.S.

1. Укорочение на 50 процентов длины канала улучшает стабильность работы AMG.

Но нужен длинный канал, на больших оборотах винт на реверсе

пробивает струей до самого входа и разрушает решение.

2. Проверка работы на тетра-сетке (без преобразования в полихедру)

показала аналогичную нестабильную работу.

Вопрос: как добиться стабильности работы AMG солвера

С уважением, Герман.

Изменено 22 июня 2009 пользователем German_psv

[DenizeN 0]

обычно AMG трогать не нужно. его параметры оптимальны для большинства задач.

в Вашем случае проблема заключается скорее в расчетной сетке, ГУ или параметрах решателя.

Почему вы уменьшили релаксации в 2 раза по всем уравнениям?

С какой целью Вы используете PISO в стационаре?

Какие параметры качества расчетной сетки?

Неплохо бы посмотреть на скриншот расчетной сетки.

[German_psv 20]

обычно AMG трогать не нужно. его параметры оптимальны для большинства задач.

в Вашем случае проблема заключается скорее в расчетной сетке, ГУ или параметрах решателя.

Почему вы уменьшили релаксации в 2 раза по всем уравнениям?

С какой целью Вы используете PISO в стационаре?

Какие параметры качества расчетной сетки?

Неплохо бы посмотреть на скриншот расчетной сетки.

Добрый день, Денис.

1. AMG я ничего не трогал, лишь включил BCGSTAB.

Улучшает стабильность, меньше срывов.

2. Опытным путем уменьшал коэфф. релаксации.

Улучшается стабильность и не срывается вычисление.

3. Понятно, что PISO применяется в нестационаре,

но я попробовал включить - считает.

Все это позволило увеличить длину канала до 20м.

4. Качество сетки: тетра-сетка построена в Гамбите,

во Флуенте преобразована в полиэдральную.

На винте по периметру размер шага от 1мм до 10мм.

На поверхности винта - до 5см. Первый слой сетки 0.4мм

В объеме MRF (там где винт) максимальный размер сетки 5см.

В канале максимальный размер ячейки 30см.

Все переходы от мелкой сетки к крупной сделаны с помощью

Size Function.

Картинки прилагаю.

[DenizeN 0]

Здравствуйте, Герман.

Меня больше интересуют параметры качаетсва сетки - скошенность, максимальное соотношение сторон элементов (больше относится к призматическим элементам у стенки).

какие схемы вы используете для других уравнений?

[German_psv 20]

Тетра-сетку делал в Гамбите. Предупреждений о некачественной сетке не было.

Но потом я ее преобразовал в полихедру во Флуенте, поэтому в основании призм

на стенках получились многоугольники, т.е. скошенность там близка к идеальной.

В личку мне сообщили об ошибке в cas-файле.

Дело в том, что периодические границе я устанавливал в Гамбите,

выполнив предварительно их линковку.

Во Флуенте я с ними ничего не делал.

Мне подсказали, что надо было их отредактировать, т.е. заменить

Translational (то что было по умолчанию) на Rotational.

По-видимому это и явилось причиной нестабильности.