Сообщения при расчёте в Fluent

[Railnolds 2]

При расчёте на каждой итерации Флюент выдаёт такие сообщения:

reversed flow in 20 faces on pressure-inlet 4.

reversed flow in 3 faces on pressure-outlet 5.

absolute pressure limited to 1.000000e+000 in 1 cells on zone 2 ->1.->2.->3.->4.<<<<

turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 138 cells

Первые 2 я так понимаю это не страшно, но 2 вторые меня смущают. При их появлении как правило в расчёте рано или поздно будет ошибка.

В чём может быть причина этих 2-х последних сообщений и как исправить, чтоб их небыло?

[Dim555 0]

По моему опыту:

reversed flow - обратное направление течения (реверс потока), т.е. поток выходит из модели на границе входа и входит в неё на границе выхода. Это не критично для расчёта, но говорит о неверности граничных условий, возникновении вихрей и т.д.

Сообщения о высокой турбулентной вязкости (turbulent viscosity) появляются при турбулентных потоках на грубой сетке. Это не критично для расчёта но нежелательно.

Сообщение absolute pressure limited to 1.000000e+000 - достигнут предел давление в 1 Па наиболее серьёзно и говорит о плохой сетке или неверных граничных условиях.

Можно: переразбить модель, провести адаптацию сетки в зонах с высокими градиентами и значениями давления и вязкости, снизить коэффициенты релоксации, изменить точность расчёта на First Order Upwing и т.д.

[Railnolds 2]

По моему опыту:

reversed flow - обратное направление течения (реверс потока), т.е. поток выходит из модели на границе входа и входит в неё на границе выхода. Это не критично для расчёта, но говорит о неверности граничных условий, возникновении вихрей и т.д.

Сообщения о высокой турбулентной вязкости (turbulent viscosity) появляются при турбулентных потоках на грубой сетке. Это не критично для расчёта но нежелательно.

Сообщение absolute pressure limited to 1.000000e+000 - достигнут предел давление в 1 Па наиболее серьёзно и говорит о плохой сетке или неверных граничных условиях.

Можно: переразбить модель, провести адаптацию сетки в зонах с высокими градиентами и значениями давления и вязкости, снизить коэффициенты релоксации, изменить точность расчёта на First Order Upwing и т.д.

Спасибо.

[chupal 0]

absolute pressure limited to 1.000000e+000 in 1 cells on zone 2 ->1.->2.->3.->4.<<<<

turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 138 cells

В чём может быть причина этих 2-х последних сообщений и как исправить, чтоб их небыло?

Лимиты по параметрам устанавливаются в Solve->Controls->Limits... Нужны они для того, чтобы помочь сходимости на начальном этапе расчета, чтобы решение не развалилось. Для улучшения сходимости, пределы можно подстраивать под свою задачу. Например, если Вы знаете примерный диапазон давлений и температур в вашей задаче (ну хотя бы порядки величин), то их желательно установить. Значения, стоящие по-умолчанию, обеспечивают слишком широкий диапазон изменения параметров и часто не могут предотвратить падения расчета в конкретной задаче. А сами сообщения говорят о том, что просто сработали какие-то из этих ограничений там-то и там-то. Но, к моменту, когда решение сойдется, они исчезнут (если пределы выбраны правильно).

[Victoria 104]

chupal

Railnolds

Подскажите, пожалуйста, порядок вычисления интегральных сил, действующих на тело при обтекании его внешним потоком.

Догадываюсь, что нужно проинтегрировать по поверхности давление, но как не промахнуться при этом.

Спасибо!

[chupal 0]

chupal

Railnolds

Подскажите, пожалуйста, порядок вычисления интегральных сил, действующих на тело при обтекании его внешним потоком.

Догадываюсь, что нужно проинтегрировать по поверхности давление, но как не промахнуться при этом.

Спасибо!

Report->Forces... Выбираем зоны, по которым нужно проинтегрировать, выбираем направление (Force vector), на которое спроектировать силу (например, если интересует Fx, нужно ввести x=1, y=0, z=0) и жмем Print.

[Railnolds 2]

Лимиты по параметрам устанавливаются в Solve->Controls->Limits... Нужны они для того, чтобы помочь сходимости на начальном этапе расчета, чтобы решение не развалилось. Для улучшения сходимости, пределы можно подстраивать под свою задачу. Например, если Вы знаете примерный диапазон давлений и температур в вашей задаче (ну хотя бы порядки величин), то их желательно установить. Значения, стоящие по-умолчанию, обеспечивают слишком широкий диапазон изменения параметров и часто не могут предотвратить падения расчета в конкретной задаче. А сами сообщения говорят о том, что просто сработали какие-то из этих ограничений там-то и там-то. Но, к моменту, когда решение сойдется, они исчезнут (если пределы выбраны правильно).

У меня симметричная 2Д модель. Когда решаю её как осисимметричную расчёт нормально идёт даже на стандартных лимитах. А если считать её как просто 2Д модель (соответственно профиль полный рисую) то выпадают эти сообщения и расчёт завершается с ошибкой Error: Floating point error: invalid number. Error Object: (). Пытаюсь разобраться в чём дело, уже и сетку переделывал - никак не выходит. Да и результаты промежуточные не такие, как в осисимметричной модели. Никак не могу понять.

chupal

Railnolds

Подскажите, пожалуйста, порядок вычисления интегральных сил, действующих на тело при обтекании его внешним потоком.

Догадываюсь, что нужно проинтегрировать по поверхности давление, но как не промахнуться при этом.

Спасибо!

Я немного в другой области :) вряд ли что-то подскажу.

[chupal 0]

У меня симметричная 2Д модель. Когда решаю её как осисимметричную расчёт нормально идёт даже на стандартных лимитах. А если считать её как просто 2Д модель (соответственно профиль полный рисую) то выпадают эти сообщения и расчёт завершается с ошибкой Error: Floating point error: invalid number. Error Object: (). Пытаюсь разобраться в чём дело, уже и сетку переделывал - никак не выходит. Да и результаты промежуточные не такие, как в осисимметричной модели. Никак не могу понять.

1. Grid->Check проходит без ошибок (имеется в виду та сетка, на которой расчет падает)?

2. Какую модель среды используете (жидкость, идеальный газ или другую)?

3. Какой солвер? (Pressure based, density based или Pressure based coupled)

4. На какой итерации происходит падение?

[Railnolds 2]

1. Grid->Check проходит без ошибок (имеется в виду та сетка, на которой расчет падает)?

2. Какую модель среды используете (жидкость, идеальный газ или другую)?

3. Какой солвер? (Pressure based, density based или Pressure based coupled)

4. На какой итерации происходит падение?

1. Да, без ошибок

2. Идеальный газ

3. Солвер Coupled, Explicit, Steady

4. На разных, в зависимости от плотности сетки.

Вот ссылка на тему http://fsapr2000.ru/index.php?showtopic=30256

Изменено 11 сентября 2009 пользователем Railnolds

[chupal 0]

1. Да, без ошибок

2. Идеальный газ

3. Солвер Coupled, Explicit, Steady

4. На разных, в зависимости от плотности сетки.

Вот ссылка на тему http://fsapr2000.ru/index.php?showtopic=30256

1. А какие граничные условия на входах и выходе?

2. Пробовали ли Вы уменьшать число Куранта (Courant number в Solve->Controls->Solution...)? Заданное по-умолчанию часто оказывается слишком большим. Поэтому на старте можно уменьшить до 0.1 ... 1.

PS: Не сразу заметил, что Вы используете Explicit Solver. А почему не Implicit? У Explicit более серьезные ограничения на Курант. Вообще, чтобы быстрей добиться сходимости, думаю, что лучше начать считать по Pressure Based solver. При этом, начать расчет можно с постоянной плотностью (задать примерно предполагаемую среднюю плотность по расчетной области). Расчет с постоянной плотностью более устойчив. Когда он начнет сходиться, можно переключиться на идеальный газ.

[Railnolds 2]

1. А какие граничные условия на входах и выходе?

2. Пробовали ли Вы уменьшать число Куранта (Courant number в Solve->Controls->Solution...)? Заданное по-умолчанию часто оказывается слишком большим. Поэтому на старте можно уменьшить до 0.1 ... 1.

PS: Не сразу заметил, что Вы используете Explicit Solver. А почему не Implicit? У Explicit более серьезные ограничения на Курант. Вообще, чтобы быстрей добиться сходимости, думаю, что лучше начать считать по Pressure Based solver. При этом, начать расчет можно с постоянной плотностью (задать примерно предполагаемую среднюю плотность по расчетной области). Расчет с постоянной плотностью более устойчив. Когда он начнет сходиться, можно переключиться на идеальный газ.

1. На входе - Prassure Inlet, на выходе Pressure Outlet. Значения давлений берутся из холодильного цикла, поэтому они известны.

2. Нет, не пробовал уменьшать. Стоит 1 по умолчанию. Для осисимметричной модели работало. Попробую уменьшить.

С Implicit осисимметричная модель толи не работала, толи выдавала не те результаты. Не помню уже точно. Да и потом при выборе модели я ориентировался на предыдущие работы по моделированию таких задач.

Переключиться на идеальный газ - Вы имеете ввиду, что в начале расчёта задать плотность постоянной по области, а через некоторое число итераций просто переключить на идеальный газ?

Pressure Based solver - что за зверь? У меня такого нет в солверах. Или я что-то не так понял?

[chupal 0]

Переключиться на идеальный газ - Вы имеете ввиду, что в начале расчёта задать плотность постоянной по области, а через некоторое число итераций просто переключить на идеальный газ?

Да, я это и имею ввиду.

Pressure Based solver - что за зверь? У меня такого нет в солверах. Или я что-то не так понял?

У Вас, видимо, старая версия. Раньше Pressure Based назывался Segregated solver.

[Railnolds 2]

Да, я это и имею ввиду.

У Вас, видимо, старая версия. Раньше Pressure Based назывался Segregated solver.

Попробую этот вариант, если ничего не поможет.

У меня 6.2.16 версия, Segregated solver там есть :)

[Victoria 104]

Report->Forces... Выбираем зоны, по которым нужно проинтегрировать, выбираем направление (Force vector), на которое спроектировать силу (например, если интересует Fx, нужно ввести x=1, y=0, z=0) и жмем Print.

Спасибо!!!

[Railnolds 2]

Ничего пока не получается. Число Куранту уменьшил, расчёт теперь идёт без ошибок, но результат тот же - бредовый. Наверное надо кардинально что-то делать с сеткой.

[chupal 0]

Ничего пока не получается. Число Куранту уменьшил, расчёт теперь идёт без ошибок, но результат тот же - бредовый. Наверное надо кардинально что-то делать с сеткой.

Попробуйте еще посчитать по Segregated Solver. Только не забудьте включить Viscous Heating в Define->Models->Viscous...