Редактируем RULES в CFX для задания специфического ГУ при расчете течения в турбомашине

Аэроупругий расчет лопатки высоконагруженного компрессора (изображение с сайта http://www.exploreenergy.eu).

Аэроупругий расчет лопатки высоконагруженного компрессора (изображение с сайта http://www.exploreenergy.eu).

Можно ли в ANSYS CFX задать зависимость полной температуры в неподвижной системе координат («Total Temperature in Stn Frame») от полной температуры во вращающейся системе координат («Total Temperature in Rel Frame»), и использовать это выражение для определения граничного условия на входе во вращающийся домен?

По-умолчанию такая возможность отсутствует в препроцессоре ANSYS CFX. Но вы можете отредактировать RULES-файл, или запустить CFX через командную строчку, указав ссылку на CCL-файл с соответствующим содержанием:

RULES:

  PARAMETER: Stationary Frame Total Temperature

  Dependency List = XYZT,rNoDim,aitern,citern,atstep, ctstep,\

  Ttot, Ttotstn,Ttotrel,acplgstep,cstagger

  END

END

Команда для запуска решателя CFX с ключом CCL выглядит так:

cfx5solve -def DEFFILE -ccl CCLFILE,

здесь CCLFILE — это обычный текстовый файл с командами и опциями, которые мы описали выше.

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

Вычисляем значение переменной Wall Distance в CFD-Post

соваСегодня мы рассмотрим с вами вопрос, который регулярно всплывает на разных форумах и которым обычно бывают озадачены все начинающие пользователи. Вы почти верно угадали, что этот вопрос связан с Y+ или чем-то подобным. Дело в том, что некоторые модели турбулентности напрямую не вычисляют расстояние от стенки до первого узла (эта переменная носит в ANSYS CFX название Wall Distance) и, соответственно, мы не можем «вытащить» значение этой переменной в постпроцессоре CFD-Post. Что делать в таких случаях?

Решение будет следующим. Создайте Additional Variable в препроцессоре CFX со следующими параметрами:

Variable Type = Unspecified
Units = [m]
Tensor = Scalar

В настройках расчетного домена выберите для Additional Variable опцию ‘Algebraic Equation’ и далее укажите Wall Distance.

Можно поступить еще проще. Создайте любое выражение, в котором будет содержаться ссылка на переменную Wall Distance. Далее используйте это выражение в качестве «монитора» (закладка Output Control). После этого решатель начнет записывать значения переменой Wall Distance в res-файл.

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

Визуализация пограничного слоя в ANSYS CFD-Post

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

Для этого обычно используют градиенты основных или производных скалярных полей. Часто используется критерий Lambda 2 – это критерий вихревого усиления и Q-критерий – инвариант тензора градиентов скорости. Отметим, что в постпроцессоре CFD-Post эти критерии присутствуют по-умолчанию и доступны через закладку Vortex Core Region.

Для анализа сжимаемых течений (высокие числа Маха) «рисуют» численные шлирен-фотографии, подобные тем, что показаны на рисунке слева.

Сегодня мы не будем подробно рассматривать вопрос визуализации структуры турбулентного потока и подходы, которые при этом используются. Ограничимся только одним вопросом, связанным с визуализацией пограничного слоя в CFD-Post. Далее немного рассуждений общего характера.

Возмущения, которые вносятся в поток обтекаемым телом (например, профилем крыла), при числах Re >>1 локализуются в относительно тонких пристеночных слоях, что позволяет использовать при инженерных расчетах модель пограничного слоя конечной толщины delat. За пределами этого слоя течение можно считать невозмущенным. Соответственно, основное влияние вязкости сосредоточено в пределах пограничного слоя, где скорость потока по нормали к стенке монотонно изменяется от нуля на стенке (условие прилипания) до скорости невозмущенного течения.

При этом за толщины слоев delatпринимают значения координат y по нормали к стенке, при которых скорость, температура и пр. в них отличаются от соответствующих значений этих величин во внешнем потоке на 1% (см. рисунок ниже, на котором показана схема развития пограничного слоя).

Laminar boundary layer

Этим «удобным» допущением обычно и пользуются для визуализации пограничного слоя в ANSYS CFD-Post. Безусловно, в данный подход заложены определенные неточности, которые могут быть сняты при использовании интегральных характеристик толщин пограничных слоев (толщины вытеснения, толщины потери импульса и пр.).

Читать далее

Как удалить «лишнюю» O-сетку в ANSYS ICEM CFD

сетки

Cегодня я покажу вам как за один раз удалить всю O-сетку, построенную в сеточном препроцессоре ANSYS ICEM CFD.

Для этого существует специальная опция, которая называется Merge Vertices. Последовательность действий следующая: открываете одноименную закладку, нажимаете на иконку, которая находится слева. Затем включаете опцию Propagate merge (иначе все ваши потуги ограничатся только двумя узлами).

Наконец, выбираете два любых «диагональных» узла О-сетки и нажимаете на кнопку [Confirm].

На этом все!о_сетка С уважением, Денис Хитрых.

Расчет углов падения/отражения частиц в Fluent

линейкаРассмотрим типичную ситуацию, характерную для большинства коммерческих и открытых CFD-кодов: пакет выполняет вычисления определенных переменных, градиентов и пр., но при этом не сохраняет их значения для последующей обработки в постпроцессоре.

Подобную ситуацию мы можем наблюдать, например, при расчете углов падения/отражения частиц при взаимодействии со стенкой в ANSYS Fluent:  углы вычисляются, но информация об их значениях не сохраняется в файлы результатов.

Доступ к этой информации можно получить с помощью соответствующих предустановленных UDF-макросов: DEFINE_DPM_EROSION и DEFINE_DPM_BC (подробнее см. документацию).

С другой стороны, косинус угла падения можно вычислить по формуле  cos_a = NV_DOT(A,V)/(NV_MAG(A)*NV_MAG(V)). Тогда угол падения Alpha равен acos(cos_a). Здесь A — площадь поверхности, с которой взаимодействует частица;  V  — вектор скорости частицы.

На этом всё!
С уважением, Денис Хитрых.

Варианты задания объемного источника с переменным размером в ANSYS CFX

Сегодня мы рассмотрит одну небольшую задачу, которая касается вопроса задания объемного источника переменного  размера на основе 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, а за пределами «шайбы» он равен нулю. Эти выражения легко переделать так, чтобы они стали зависимыми от времени.

Читать далее

Моделирование процесса адсорбции в ANSYS Fluent

Сегодня мы рассмотрим задачу моделирования процесса адсорбции c использованием решателя ANSYS Fluent. Напомним, что адсорбция – это процесс поглощения поверхностно-активных веществ из раствора или газовой смеси поверхностным слоем жидкости или твердого тела – специальными адсорбентами. В качестве адсорбентов могут выступать различные материалы с высокой удельной поверхностью, например, пористый углерод (активированный уголь) или силикагели.

Сgeometryледует отличать адсорбцию от абсорбции. В последнем случае вещество диффундирует в объем жидкости или твердого тела и образует раствор или гель.

Адсорбция – достаточно распространенное явление в химической и нефтегазовой отраслях. Оно возникает повсеместно, где присутствует поверхность раздела между фазами. Устройство для проведения адсорбции называется адсорбером.

Для моделирования адсорбции мы будем использовать предопределенные шаблоны для задания реакций на поверхностности и гетерогенных реакций по всему объёму в ANSYS Fluent.

Количественно процесс адсорбции (без учета межмолекулярного взаимодействия адсорбата) описывается уравнением Ленгмюра:

2018-07-10_12-37-57где Θ – доля площади поверхности адсорбента, занятая адсорбатом, α – адсорбционный коэффициент Ленгмюра, P – концентрация адсорбтива.

Уравнение Ленгмюра представляет собой одну из форм уравнения изотермы адсорбции, которая определяет зависимость равновесной величины адсорбции от концентрации адсорбтива при постоянной температуре. Концентрация адсорбтива при адсорбции из жидкости (растворов) выражается в мольных долях. При адсорбции из газовой фазы концентрация выражается в единицах абсолютного давления. В случае адсорбции на твёрдых адсорбентах, используют отношение массы поглощённого вещества к массе самого адсорбента.

Рассмотрим двумерную задачу адсорбции газовой смеси (воздух+углекислый газ) твердым телом (слоем насадки адсорбера).

Поглощение углекислого газа адсорбентом описывается следующим уравнением:
2018-07-10_12-50-14

Далее опишем основные шаги, связанные с постановкой задачи адсорбции в ANSYS Fluent.

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

Включите модель переноса компонентов Enable Species transport и в открывшейся панели Species Model активируйте в разделе Reactions опции Volumetric (реакция в объеме) и Wall Surface (реакция на поверхности).

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