Задача о нагреве цилиндрической емкости

[Кадников 0]

Помогите разобраться, откуда берутся такие результаты.

Рассматривается плоская нестационарная задача о нагреве цилиндрической емкости.

В начальный момент температура воздуха в емкости T0. Температура стенок поддерживается T1>T0. Рассматривается режим естественной конвекции.

Вроде бы начальные и граничные условия задал правильно. Задача абсолютно симметрична, а на выходе система выдает, что правая половина емкости прогревается медленне чем левая и картина температурного поля ассиметрична. Линии тока по форме вроде похожи на правду.

Подскажите, пожалуйста, в чем может быть дело

Примечание

1. Результат сравниваю с расчетами и экспериментами МВТУ

2. При решении пришлось применить метод Modified Inertial Relaxation(не настолько силен в языке, чтоб точно перевести :) ) по температуре с коэффициентом 0.1. Насколько я понимаю этот метод искусственно делает положительными диагональные коэффициенты. Которые иначе получались отрицательными.

Заранее благодарен

[takedo 1]

Как то продавцы показывали возможности пятой серии Ansys(толи настран толи насран) в области газодинамики. Так вот предложил я им решить симметричную задачу продува трубы... 1) в пятых верстиях точно не было сииметричного решения!! 2) было явно видно, что результаты расчета ложаться по линиям сетки, что говорило о первом порядке кон. разн. схем.

Исходя из изложенного могу лишь посоветовать решить задачку в последнем Ansys, ибо те же продавцы говорили, что уже в версии 6 с этими делами получше стало

А можно и во FloWorks е или FlowVision е, CFX(тоже Ansys), FLUENT сделать тоже самое

[lks 0]

Задача абсолютно симметрична, а на выходе система выдает, что правая половина емкости прогревается медленне чем левая и картина температурного поля ассиметрична

Задача абсолютно симметрична, а симметрична ли сетка?

[Кадников 0]

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

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

А вот пробовал решить ту же задачу для кольцевого зазора - получаются совсем бредовые результаты с отрицательными температурами

[Mocus 0]

Отрицательные температуры ? Измельчай сетку !!!