В Star-CD версии 4.06, а также в ранних релизах версии 4.08, неправильно вычисляется полное давление при пост-процессинге. Макрос вычисляет полное давление по изэнтропическому соотношению.
Согласно изэнтропическим соотношениям полное давление (или давление торможения) может быть рассчитано следующим образом:
где T - термодинамическая температура, T0 - полная температура (температура торможения), ? - показатель адиабаты, pt - термодинамическое давление.
Показатель адиабаты - отношение теплоемкостей при постоянном давлении cp и при постоянном объеме cv, с учетом соотношения Майера может быть записан как:
где R - универсальная газовая постоянная (8314,3 Дж/кмоль•К), М - молярный вес (в макросе задано значение для воздуха - 28,96 кг/кмоль).
Полная температура (температура торможения) и термодинамическое давление рассчитываются самим постпроцессором. Они находятся в разделе Calculated Scalar Data, называются соответственно Total Temperature (имя параметра для использования в командах GETC и GETV - TTOT) и Therm. Pressure (имя параметра для использования в командах GETC и GETV - PTER).
Термодинамическое давление, если кому интересно, отличается от «обычного» (пьезометрического) давления, определяемого в ходе решения уравнений Навье-Стокса, и рассчитывается по формуле:
где ?0 - опорное значение плотности (Reference Density), g - ускорение свободного падения, z - превышение над тем уровнем, для которого задано опорное значение плотности, ? - местное значение плотности:
,
k - удельная кинетическая энергия турбулентности (соответственно, для ламинарного течения последнее слагаемое обращается в ноль).
В макросе итоговое значение полного давления вычисляется как избыточное по отношению к атмосферному (атмосферное берется равным 100000 Па). Если нужно абсолютное значение давления - закомментируйте последние две строчки.
Полное давление - стандартный расчет через постпроцессор.
Полное давление - работа макроса:
Скачать Star-CD макрос вычисление полного давления