Сходимость

[Apologist 0]

3. Второй порядок, ибо первым порядком все размажем и решение будет нефизичным.

Первый порядок нормальное решение дает. Потом проходим вторым порядком чтобы отлакировать.

Методика:

1. Считаем с постоянной плотностью, первый порядок, получаем очень грубое (и местами нефизичное) решение.

2. Добавляем сжимаемость и энергию. Считаем, всё норм.

3. Включаем второй порядок. Считаем, сходится.

a_schelyaev

Массовый расход какой на входе/выходе? Выход не удлинняли? В моей геометрии он вынесен на 10 метров. Сетка на сколько ячеек?

Изменено 5 августа 2010 пользователем Apologist

[a_schelyaev 366]

1. Вход не удлинял.

2. На выходе прикрутил 5 метровый короб.

3. Сходиться начинает уже на грубой сетке. Если мельчить, то сходимость продолжается. Дошел до 80000 ячеек, как минимум.

4. Расход не смотрел. Задачу уже остановил, т.к. счет идет нормальным. А если считать, то нужно смотреть сходимость по сетке, но это уже на более подробной сетке значит и тут сходимость тоже будет идти нормально.

5. Осциляции в невязках могут появиться на более подробной сетке из-за появившихся срывных вихрей.

6. Никаких допущений не делал: ни первого порядка, ни несжимаемой жидкости.

P.S. И всем советую первые несколько шагов считать строго на Куранте 1, чтобы посмотреть что да как.

Первый порядок нормальное решение дает. Потом проходим вторым порядком чтобы отлакировать.

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

2. Если задача стационарная, то такая методика Upwind + High Order вполне подойдет. Но если задача нестационарная, то нужно уже с первой итерации уже иметь получать решение. В этой ситуации такой подход не прокатит.

[Apologist 0]

Если задача нстационарная - то надо решать нестационар, а не париться :)

У меня на данный момент пара проблем:

1. На выходе поток немного "гуляет" по коробу. Из-за этого массовый расход на выходе нестабилен + нестабильные обратные вихри. Хотя продолжает нормально сходиться.

2. Не нравится распределение температур.

Изменено 5 августа 2010 пользователем Apologist

[konstantinm 0]

Вот такую картину получил при скорости 50 м/с. Видимо имеется какая-то взаимосвязь со скоростью?

Изменено 5 августа 2010 пользователем konstantinm

[konstantinm 0]

Еще вот что обнаружилось: оказывается, имеет значение то, как создается таблица с данными. Я создавал в Excel, а потом сохранял в *.csv. А потом посмотрел в tutorial; они там создают таблицу с помощью блокнота. Создал табличку с помощью блокнота.

Также изменил модель течения на k-Omega Turbulence + SST (Menter) K-Omega.

Задача сходится если в табличке заданы распределения скоростей только по четырем радиусам (вверх, вниз, влево, вправо). Если дополнительно задать распределения скоростей по радиусам, которые под углом 45 град, то задача - не сходится.

Изменено 12 августа 2010 пользователем konstantinm

[a_schelyaev 366]

И Константин берет гитарууу...

А ты попробуй вывести распределение того, что ты выводишь. Вдруг какие-то скачки появятся - мало-ли.

Ну заодно выглади профиль свой - тоже не помешает.