Second order upwind

[flinstone 0]

задача истечение газа из сопла,истечение струи большой нерасчетности в вакуум. Если считать задачу,сначала используя first order upwind все отлично сходится задача считается,если переключить/начать считать second order upwind то через некоторое время взависимости от числа куранта появляются какие то точки в коорых задча расходится,будь то мах или температура.

выглядят они вот так

Перепробовал все что можно,наверно,разную сетку,разные начальные давления,понижал его,понижать under-relax factor,считать при очень малых числах куратна порядка 0.0001,все равно через некоторое время задача расходится,появляются точки с громадными числами Маха около 120,или температуры свяше 5000 К, после чего задача вылетает. При этом при first order все нормально считается,но результат по точности неудовлетворяет.

Сетка выглядит так

ГУ

pressuare inlet задается мах 2.5

presuuare outlet 1па

symmetry

operating conditions 0

решается density based implicit

ideal gas

В чем может быть ошибка,в сетки?может что то не так задаю?

[flinstone 0]

Неужели никто даже предположить не может в какую сторону копать? Никак не получается свести подобные задачи (истечение газа из сопла\сопел в вакуум) на второй точности.Для задач внешнего обтекания все нормально получается на second upwind.Не понимаю что в этих задачах не так делаю. В чем математически отличается first order от second order?Почему на first order когда задача начинает расходится (time step reduce) продолжается расчет и через некоторое время все нормально,а на second order задача почти сразу вылетает как появляется time step reduce?

Изменено 21 марта 2011 пользователем flinstone

[Yuri69 1]

Схема первого порядка привносит большую искусственную вязкость => алгоритм на ее основе более устойчив.

Если с 1-м порядком считается без проблем, то можно через некоторое число итераций попробовать переключиться на схему 2-го порядка.

А вы в области струи делаете сгущение сетки?

[flinstone 0]

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

Нет,сгущение в области струи не делал,только около стенки,срез сопла и рядом с ним. Сейчас как раз этим занимаюсь пытаюсь построить более мелкую сетку в области струи.

[dimaZ 1]

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

Нет,сгущение в области струи не делал,только около стенки,срез сопла и рядом с ним. Сейчас как раз этим занимаюсь пытаюсь построить более мелкую сетку в области струи.

какой у вас решатель? pressure-based или density-based?

[flinstone 0]

ГУ

pressuare inlet задается мах 2.5

presuuare outlet 1па

symmetry

operating conditions 0

решается density based implicit

ideal gas