Не сходиться по массе газодинамическая задача

[San24 0]

Здравствуйте!

Решаю следующую задачу. Область примерных размеров 15*4*0.5 м. в ней установлены несколько горелок вдоль длинных сторон. Горелки сделаны просто кубиками в нутри газовой области, площадь поверхности равна площади выхода горелки диаметром 70 мм. На одной из поверхностей каждого куба задан INLET и скорость порядка 15 м/с. На верхней поверхности заданы несколько OUTLET, размером примерно 15*15 см. Модель турбулентности стандартная K-E. Основная среда воздух, зависимость плотности от температуры.

Задача не сходиться по массе, невязка опускается примерно до значения 1e-1 +-.

задача не сходится даже если просто оставить только движение воздуха без температуры, горения и при модели турбулентности K-const.

Пакет Star-CD 3.22

Пожалуйста, посоветуйте, что можно сделать, где посмотреть.

Изменено 8 августа 2005 пользователем San24

[Miki 3]

стоит начать с проверки качества сетки, может где то кривые ячейки сделались

и кстати не указан тип сетки ? гекса, тетра ?

например Стар крайне не любит тетры

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

[San24 0]

стоит начать с проверки качества сетки, может где то кривые ячейки сделались

и кстати не указан тип сетки ? гекса, тетра ? 

например Стар крайне не любит тетры

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

<{POST_SNAPBACK}>

1. В сетке действительно есть не совсем хорошие элементы в углах, но я вывел отношение разности давлений через 20 итераций в долях к первому давлению, т.е. (P300 - P320 )/P300 (300 и 320 номера итераций) и самая большая разница получилась в нутри сетки, где ячейки близкие к идеаньному кубику. Могут ли возникать большие осциляции решения в нормальных ячейках если есть плхие в другом месте?

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

3. Сетка Hexa построеная в ICEM, но т.к. с ICEM я еще плохо знаком, то хорошо построить сетку в скругленном углу не получилось.

Изменено 8 августа 2005 пользователем San24

[BigBrother 1]

Здравствуйте. На всех Outletах заданы равные расходы? Попробуйте заменить Outlet на Pressure с нулевым давлением.

[San24 0]

Здравствуйте. На всех Outletах заданы равные расходы? Попробуйте заменить Outlet на Pressure с нулевым давлением.

<{POST_SNAPBACK}>

Да я пробывал и с 0 pressure, сходимость была таже. А outlet я задавал обозвав несколько разрозненых областей одним Outlet-ом. Сходимость улучшилась примерно до 2e-2, когда я исключил из модели перефирийную часть с наклонными элементами. Еще немного улучшилась сходимость после того как выходы екоторые были в плоскости я вытянул см. рисунок. Как мне подсказали программа не любит когда outlet или pressure 0 находятся в одной плоскости со стенкой.

[San24 0]

Второй рисунок не влез в тоже сообщение

[BigBrother 1]

Судя по рисункам, надо сетку сделать гуще. На втором рисунке красный "аппендикс" с разбиением 3*5 элементов - это не хорошо. У него боковые грани - стенки? Тогда в такой трубке не получится правильный профиль скорости и решаться будет плохо.

Сам как-то пробовал решить задачу с ламинарным течением в тонком зазоре. Задал 6 элементов по высоте зазора. И никак не хотела решаться. Увеличил до 12 и всё ок.