Странный глюк Gambit-а

[Yuri69 1]

При создании и присоединении пограничного слоя к любой edge, принадлежащей некоторой face, помимо собственно слоя возникает еще сбоку побочный слой ячеек - на одной из соседних edge из той же face (показано черными стрелками на рисунке):

<noindex></noindex>

И ничем сей артефакт не удается убрать!

[Iatsouk 3]

Что-то странное нарисовано.

Если выложите модель и лог можно будет посмотреть.

[Yuri69 1]

Запускается следующий скрипт. А на рисунке - одна из подобластей в препарированном виде...

geom1.zip

Изменено 29 сентября 2008 пользователем Yuri69

[Yuri69 1]

German_psv: спасибо за помощь. Суть действительно в генерации всех ПС до мешинга edges и faces. Такой порядок можно бы было узнать и из расположения гамбитовых кнопок =)) Единственно, немного не понял, как одновременно сгенерить ВСЕ ПС ОДНОЙ командой, на что вы вроде намекаете?. Что до хранения модели во внешнем CAD-формате, то для моего уровня владения пакетами такой подход будет неоправданным усложнением. (Кстати скажу, что сам скрипт я создаю с помощью программы на Фортран, этаком пре-препроцессоре, ибо пустыми руками здесь не навоюешь, важна именно автоматизация перестроения контура при внесении в него исправлений)...

Еще пару вопросов:

1) Создание ПС на одной из edge (верхнем контуре выходной области "dif_throat") приводит при генерации сетки для "dif_throat" к ухудшению ее качества, в виде островков самопересечений (где отрицательный Cell Volume, как потом глаголет Fluent). Не помогает ни изменение типа элементов, ни схемы мешинга. Пример приведу чуть позже. А пока просто констатация факта. Как бороться с сим явлением, не отказываясь от ПС на данной стенке?

2) У меня стоит Gambit 1.3 (вроде бы) и кондовость интерфейса с частыми вылетаниями из-за ошибок в "hpp"-файлах просто поражают. Стоит ли в этом смысле переходить на более новую версию программы (скажем, 2.3.16)?

Изменено 29 сентября 2008 пользователем Yuri69

[Yuri69 1]

Вроде разобрался с проблемой отрицательного объема некоторых ячеек сетки: перешел в проблемной области на метод triangle...

Кстати, вопросов не становится меньше =)) Моделирую стационарное турбулентное течение в приведенном выше аэродинамическом тракте на основе k-eps модели с использованием Enhanced Wall Treatment. Величину y+ поддерживаю на всех стенках равной примерно 3 (что соотносится с рекомендациями y+<4-5). Воздух прокачивается по каналу вентилятором (fan) и истекает в окружающую среду. Имеются местные потери, которые представлены пористым слоем (porous zone) и скачком пористости (porous jump). В граничных условиях на входе (pressure inlet) и выходе (pressure outlet) задаю интенсивность турбулентости I=1% и соответствующий гидравлический диаметр. Расчет турбулентного течения осуществляю на основе стационарного решения, полученного в рамках ламинарной модели.

В результате получаю аномально высокий уровень турбулентности: I=100% и выше в следе за втулкой вентилятора и в струе, не говоря уже о пристенных областях! Что-то где-то не так делаю? Или трудности в разрешении физики? (Течение в области выхлопа безотрывное, а за центральным телом - отрывное, но ведь со следовыми возмущениями борется сетка, расположенная вниз по потоку)

Изменено 3 октября 2008 пользователем Yuri69

[Yuri69 1]

German_psv:

На входе задаю избыточное давление скажем (-100)Па, а на выходе 0Па и соответствующие турбулентные параметры. НО, насколько я понял, при нахождении стационарного решения эти величины используются только КАК ЗАТРАВОЧНЫЕ, т.е. чтобы стартануть итерационный процесс. В итоге, какое бы избыточное давление на входе не задать, все равно получишь ОДНО-единственное стационарное решение, определяемое только характеристикой вентилятора dp=f1(v) (зависимостью повышения давления за вентилятором от средней скорости перед ним) и гидравлическими потерями по тракту. Ведь, как известно из гидравлики, рабочий (установившийся) режим трубы будет определяться точкой пересечения двух графиков - характеристики нагнетателя dp=f1(v) (f1-функция убывающая) и кривой потерь dp=f2(v) (f2-возрастающая функция скорости v).

Насчет фейсов с 4-мя ребрами - это не так для "edge_throat_2", если я правильно понял ваше замечание.

А второе соображение уже приводил: для выходного устройства использование именно triangle позволяет избежать проблем с отрицательным cell volume некоторых ячеек.

Изменено 8 октября 2008 пользователем Yuri69

[Yuri69 1]

"...

Нет, граничные условия действуют постоянно..."

--> тем не менее в задаче с вентилятором на входе вырабатывается одно и то же давление, вне зависимости от того значения, что ставится на inlet. Причем, чем дальше эта граница отнесена от оси трубы, тем меньше статическое давление на ней (в среднем), p->-0 (слева к нулю), т.к. условие о покое окружающей среды выполняется все лучше...

"...

Теперь можно, увеличивая размер картинки..."

--> эта опция у меня почему-то не работает (пытаюсь приближать и правой и центральной кнопкой мышки)

[Yuri69 1]

German_psv:

Спасибо про Mouse бутоны =)) Бум пробовать.