Турбулентность DES

[Bonusfrag 104]

Шаг интегрирования?

Критерий сходимости? Невязки? До какого уровня? И до какой точности по умолчанию считаются уравнения? 1%, больше, меньше?

я через личку всё уточню)

чтоб тут не получилось чата.

Потом подробно всё выложу, с картинками.

[Гость suvorov.as]

я через личку всё уточню)

Сетка чисто 3d. В cfx нет 2d.

На боковые грани накладывается условие симметрии и задача становится как бы 2D. Не совсем конечно, но для этого теста неважно.

Шаг интегрирования?

Стационарная задача. Код автоматически подбирал.

Критерий сходимости? Невязки? До какого уровня? И до какой точности по умолчанию считаются уравнения? 1%, больше, меньше?

Я там написал, что эти настройки по умолчанию.

В CFX это невязки RMS по скорости и давлению 1e-4.

Во флюенте х.з. там не поймешь. Написано для всего 0.001, но течение там глобально продолжало чуть-чуть меняться и после достижения этих критериев. Так, что предлагаю ориентироваться на CFX или на глазок .

[Bonusfrag 104]

Стационарная задача. Код автоматически подбирал.

В FV нестационарная решалка: надо указывать либо ШПВ, либо CFL.

И менять их можно только вручную.

А в ANSYS он вроде бы эти параметры меняет сам, ускоряя решение.

Так что сравнивать на постоянном шаге надо.

Вот картиночка, но это на 2D - 380 тыс.ячеек.

2D сетка 150*800, в доп.объемах адаптация 1 и 2 уровня, в погране - 3го уровня.

[a_schelyaev 367]

Сетка чисто 3d. В cfx нет 2d.

В генераторе сетки 1 ячейку поперек и размер получить нужный именно начальной сеткой, чтобы потом никаких адаптаций не делать.

Стационарная задача. Код автоматически подбирал.

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

Также в выдаче он указываться должен.

Во флюенте х.з. там не поймешь. Написано для всего 0.001, но течение там глобально продолжало чуть-чуть меняться и после достижения этих критериев. Так, что предлагаю ориентироваться на CFX или на глазок .

Просто в CFD по невязкам задачи решать бессмысленно, особенно если вскрывается нестационарность.

Нужно смотреть сходимость по физическим характеристикам, например, по продольной силе (сопротивление).

[Bonusfrag 104]

невязки

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

[Гость suvorov.as]

В генераторе сетки 1 ячейку поперек и размер получить нужный именно начальной сеткой, чтобы потом никаких адаптаций не делать.

Так никаких автоадаптаций нет. Сетка такая, какая она сгенерирована изначально. Каринка сечения выше приведена.

Так что сравнивать на постоянном шаге надо.

Так я не против постоянного шага. Ансис автоматом выбирает этот постоянный шаг. Вечером могу его указать.

Вот картиночка, но это на 2D - 380 тыс.ячеек.

Нет, 2D не интересно сравнивать. Алгоритм другой и другие нюансы.

Нужно смотреть сходимость по физическим характеристикам, например, по продольной силе (сопротивление).

Я не спорю. Можно ряд контрольных точек заложить для скорости, для давления, допустим.

[Bonusfrag 104]

Нет, 2D не интересно сравнивать. Алгоритм другой и другие нюансы.

алгоритм и уравнения - 3D. 2D только сетка.

Сетки сделаем одинаковыми, 3х-мерными.

Сейчас немного раскидаю текучку - распишу подробно.

Например, в FV только параллелепипеды в сетке, т.е. перейти с 4 ячеек на 10 я никак не могу. Могу только на 8 перейти.

[Гость suvorov.as]

Могу только на 8 перейти.

Хорошо. Поскольку, по моему мнению, сгущение сетки 1 к 8 в 3D это не то же самое, что 1 к 4 в 2D.

Сейчас немного раскидаю текучку - распишу подробно.

ОК. Не горит.

[FSA-FSA 249]

На фоне дохода от Fluent батька-CFX казался сопливым мальчуганом.

ОБС много разработчиков из флюента в стар ссм перешло. Верно говорят? И что Ансис сам решателей уже давно сильно не менял?

[Гость suvorov.as]

ОБС много разработчиков из флюента в стар ссм перешло.

стар ссм очень похож на флюент. При этом, в отличае от флюента, имеет user friendly interface.

[a_schelyaev 367]

ОБС много разработчиков из флюента в стар ссм перешло. Верно говорят? И что Ансис сам решателей уже давно сильно не менял?

Тоже слышал такую байку и что офис чуть ли не по другую сторону улицы был сделан на первых порах.

Фиг знает.

Ансис решатели не меняет, он их дорабатывает.

Если что-то и изменят серьезно, то только под GPU или MIC.

[vl 327]

Если что-то и изменят серьезно, то только под GPU или MIC.

Что есть MIC?

[a_schelyaev 367]

Зеон Фи, он же Майк, в девичестве Лараби

Мэни Интегрейтед Корес

[Гость suvorov.as]

Шаг 0,033с

[Bonusfrag 104]

И так:

опорная сетка 1х60х320. Внутри первого дополнительного тела - адаптация 1го уровня (делим ячейку на 8 маленьких), внутри второго тела - адаптация 2го уровня. Погран-слой на всей нижней стенке дополнительно проадаптирован - толщина 12 ячеек 3го уровня адаптации.

Т.е. поперек слоя сверху вниз - 1 / 2 /4 / 8 ячеек.

Итого в сетке 375504 ячейки

сверху и по бокам симметрия, снизу стенка с прилипанием и пристенной функцией.

При пересчете из k/e у меня получились пульсации 0.001721 и масштаб турбулентности 0.00284.

Начальная скорость 0 м/с, граничная 15 м/с (задана нормальная массовая скорость 15000 кг/м-2 с-2)

3D URANS + k-e модель, неявная схема 2го порядка точности, шаг по времени 0,033 с

осталось разобраться с точностью решения внутренних уравнений и с критерием сходимости

я никогда ранее автоостановкой не пользовался.

чего я забыл?

[Гость suvorov.as]

И так:

опорная сетка 1х60х320. Внутри первого дополнительного тела - адаптация 1го уровня (делим ячейку на 8 маленьких), внутри второго тела - адаптация 2го уровня. Погран-слой на всей нижней стенке дополнительно проадаптирован - толщина 12 ячеек 3го уровня адаптации.

Т.е. поперек слоя сверху вниз - 1 / 2 /4 / 8 ячеек.

Итого в сетке 375504 ячейки

сверху и по бокам симметрия, снизу стенка с прилипанием и пристенной функцией.

При пересчете из k/e у меня получились пульсации 0.001721 и масштаб турбулентности 0.00284.

Начальная скорость 0 м/с, граничная 15 м/с (задана нормальная массовая скорость 15000 кг/м-2 с-2)

3D URANS + k-e модель, неявная схема 2го порядка точности, шаг по времени 0,033 с

осталось разобраться с точностью решения внутренних уравнений и с критерием сходимости

я никогда ранее автоостановкой не пользовался.

чего я забыл?

Я правильно понял, это тридэ?

Вроде все так .

[a_schelyaev 367]

Если будешь по невязкам стоперы ставить, то закажи невязки в закладке Солвер, а в каждой из них пропиши нужный уровень 1e-4 например.

Но лучше по физическим характеристикам смотрите.

Или максимальное давление в области или силу/интеграл на стенке.

Количество итераций для уравнений задается в доп.параметрах.

По умолчанию стоит 2000, но это на случай ядерной войны, поэтому можешь не трогать. Количество решенных итераций для уравнений увидишь в окошке Мониторинга. По идее только для давления будет больше итераций.

Точность задается по умолчанию 1%.

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

Дома у меня на стареньком Коредуо (еще первом) с двумя ядрами я получал 130 секунд на итерацию, а потом 90.

Также стоит обратить внимание на тип решателя.

В CFX стоит по умолчанию совместная решалка, если я правильно понимаю, то это аналог Pressure Based Coupled в Fluent.

В Струегляде, скорее всего, аналог Сегрегейта Fluenta.

По замерам, выполненных в Сарове, CFX делал всех на простых задачках из-за алгебраического мультигрида (AMG), особенно, если все градиентные зоны разрешены сеткой.

Поэтому, сперва можно сравнить в лоб при схожих параметрах, а потом показать кто на что способен - залупить шаг интегрирования побольше или задействоать в FV Совместную схему с двумя внутренними итерациями.

Поглядим какого цвета у него кишки.

;)

Да, и процессоры выставьте в одинаковый режим относительно HT.

[Гость suvorov.as]

По умолчанию стоит 2000, но это на случай ядерной войны, поэтому можешь не трогать.

[a_schelyaev 367]

чего я забыл?

Слои не строй. И характеристику для мониторинга сходимости бери по объему, а не по стенке - быстрее будет.

Сечение расчетной сетки тоже убери.

Адаптацию выключи после построения сетки.

;)

Я правильно понял, это тридэ?

Саш, сетку поперек покажи, чтобы вопросов не было.

[Гость suvorov.as]

Да, и процессоры выставьте в одинаковый режим относительно HT.

Это что значит?

чтобы вопросов не было.

Да ладно, я наслово поверю.