У Джефри Чосера — «любимого» поэта английских королей — есть хороший афоризм про время: «Time and tide wait for no man». Тема скоротечности времени всплывает и в речи M, которая дает достойный отпор всем недоброжелателям МИ-6 в зале суда.

Что объединяет корифеев английской поэзии с вычислительной гидродинамикой? Ответ простой. Любая задача, которая решается в нестационарной постановке, т. е. в которой дискретизации подвергаются не только независимые переменные, относящиеся к пространству, но и переменная время. Поэтому сегодня мы рассмотрим с вами различные вопросы, связанные с решением нестационарных задач в газодинамических пакетах ANSYS CFX и ANSYS Fluent. Мы поговорим об обработке результатов нестационарных расчетов, программировании и многом другом. Большая часть материалов данного поста основана на реальных запросах пользователей в службу технической поддержки компаний ANSYS и КАДФЕМ.

Вопрос 1. Как остановить нестационарный расчет в ANSYS Fluent, например, если температура в заданной точке (узле, средняя по объему, поверхности и пр.) достигла пороговой величины? Для этого можно использовать следующий UDF-макрос.

Для использования данного макроса вам необходимо его откомпилировать (проверено для версии не ниже 17.0). И определить следующую команду (в строке exceute command): (if (= (%rpgetvar ‘converged) 1) (begin (ti-menu-load-string «/file/wcd low-temp-trig-%t.cas») (ti-menu-load-string «exit»))).

Как работает данный UDF-макрос? В процессе итераций выполняется проверка на достижении выбранной переменной заданного порогового значения. Если условие выполняется, то записываются *.cas и *.dat файлы.  И выполняется команда выхода.

Сегодня я успел рассмотреть только один вопрос. Но это только начало. В следующем году обязательно продолжим.

Всех с наступающим Новым годом и Рождеством!

С уважением, Денис П. Хитрых,
Директор АО «СимуЛабс».

Сегодня я расскажу о планах  CFD-blog.ru на 2016-й високосный год. Проанализировав все запросы, которые поступали в мой адрес и на форму «Как с нами связаться» в текущем году, я пришел к выводу, что уровень информационной поддержки таких направлений как, турбомашиностроение, моделирование многофазных потоков и моделирование процессов теплообмена, следует существенно увеличить. Поэтому в наступающем году я буду больше внимания уделять этим вопросам.

Кроме того, мы приняли решение не ограничиваться только тематикой вычислительной гидрогазодинамики, но с  учетом возросшего интереса пользователей к многодисциплинарным задачам, расширить тематику блога в сторону прочностных задач, т. е. традиционного ANSYS Mechanical.

Теперь чуть более подробнее о планах. Во-первых, изучать многофазные потоки следует с самых азов. Поэтому я начну публиковать в открытом доступе свой курс по многофазным потокам, основанный на лекция проф. Алана Бёрнса (Великобритания) и моих собственных методических материалах.

Во-вторых, мы совместно с вами решим несколько задач с использованием самых передовых технологий ANSYS, заложенных в Turbo-линейку (TBR-методы, SBES/SDES модели турбулентности, модели реальных газов и пр.).

Наконец, мы приобщимся к миру многокритериальной оптимизации и расширенной кастомизации на основе ACT-технологии (в этом нам помогут специалисты R&D Центра SimuLabs4D).

Год предстоит сложный и насыщенный. Надеюсь, что часть нашего объемного плана реализуется и принесет вам пользу.

С уважением, Денис Хитрых. До встречи в новом году.

Сегодня мы рассмотрит одну небольшую задачу, которая касается вопроса задания объемного источника переменного  размера на основе subdomain в ANSYS CFX.

Как вы знаете, в ANSYS CFX поддомен всегда имеет постоянный размер, но вы можете определить поддомен целиком для расчетной области и использовать CEL-выражения для задания объемного источника.

Предположим, что наш источник имеет форму «шайбы» радиусом 0.3 м и высотой 0.3 м. Соответствующий список выражений будет иметь вид:

z0 = 0.0 [m]
x0 = 0.5 [m]
y0 = 0.5 [m]
sourceRadius = 0.3 [m]
sourceHeight = 0.3 [m]
unitL = 1 [m]
rSource = sqrt((x-x0)^2 + (y-y0)^2)
hSource = z-z0
insideSource = step((sourceRadius-rSource)/unitL)*step((sourceHeight-hSource)unitL)
source = insideSource * 100 [W m^-3]

Соответственно, в пределах «шайбы» объемный источник энергии равен 100 Вт/м3, а за пределами «шайбы» он равен нулю. Эти выражения легко переделать так, чтобы они стали зависимыми от времени.

Читать далее →