Анализ энергоэффективности системы кондиционирования воздуха пассажирского самолета энтропийно-статистическим методом

Полный текст

Условные обозначения

P - давление (кПа)

Rв - газовая постоянна для воздуха (287 Дж×кг^-1×K^-1)

T - температура (K)

k - оказатель адиабаты для воздуха (1,41)

i - Энтальпия (Дж×кг^-1)

S - Энтропия (Дж×кг–1×K^-1)

G - Расход воздуха (кг×с^-1)  

η - КПД узла

ОБОСНОВАНИЕ

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

Существует большое количество возможных вариантов схемных решений СКВ. В данной статье анализируется схема с трехколесным турбохолодильным агрегатом (турбохолодильником) с петлевым влагоотделением на высоком давлении. Влага имеет большое влияние на эффективность работы СКВ [1]. Но, помимо этого, она так же влияет на конструктивные особенности и варианты режимов работы СКВ. Наличие влаги при работе СКВ на отрицательных температурах приводит к ухудшению эффективности теплообменных аппаратов вследствие их обмерзания, что вынуждает вводить дополнительные агрегаты для решения этой проблемы, а также ограничивает диапазон рабочих режимов. Схема СКВ с петлей влагоотделения на высоком давлении является одним из вариантов борьбы с негативным влиянием влаги на СКВ. Ниже приведено описание математической модели подобной схемы для пассажирского самолета. Математические модели часто применяют для холодильных и криогенных систем [2, 3]. Рассматриваемая в настоящей работе модель позволяет проводить статический расчет СКВ на всех эксплуатационных режимах. На основе этой модели был проведен энтропийно-статистический анализ системы на выбранном режиме полета самолета. Данный метод был разработан в МГТУ им. Н.Э.Баумана [4]. На основе указанного метода активно проводят анализы различных криогенных систем [5–7]. Энтропийно-статистический метод анализа системы позволяет определить влияние эффективности отдельных агрегатов из состава СКВ на эффективность СКВ в целом.

ЦЕЛЬ РАБОТЫ

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

ОПИСАНИЕ РАБОТЫ СИСТЕМЫ

СКВ пассажирского самолета — это система, предназначенная для регулирования температуры, расхода и влагосодержания воздуха, подаваемого в самолет. В данной статье рассматривается работа СКВ с трехколёсным турбохолодильным агрегатом (турбохолодильником) и петлевой схемой влагоотделения на высоком давлении. Схема СКВ, анализ которой проводится в статье, представлена на рис. 1.

 

Рис. 1. Схема СКВ. В — ступень вентилятора турбохолодильника; ВО — влагоотделитель; К — компрессорная ступень турбохолодильника; ОК –— обратный клапан; РУ — регулирующее устройство для перепуска рабочего воздуха; Т — детандерная ступень турбохолодильника; ТОА 1 — первичный теплообменный аппарат; ТОА 2 — Вторичный теплообменный аппарат; ТОА 3 — перегреватель; ТОА 4 — конденсатор.

Fig. 1. ACS scheme. V — air cycle machine fan stage; VO — water collector; K — compressor stage of the air cycle machine; OK — check valve; RU — control device for bypassing working air; T — expander stage of the air cycle machine; TOA 1 — primary heat exchanger; TOA 2 — secondary heat exchanger; TOA 3 — reheater; ТОА 4 — condenser.

 

Источником сжатого воздуха в СКВ является компрессор газотурбинного двигателя (ГТД) самолета.

По рабочему тракту воздух поступает в первичный теплообменный аппарат (ТОА1), где происходит его охлаждение забортным воздухом. Пройдя его, рабочий воздух поступает в компрессорную ступень турбохолодильника (К), где происходит его дополнительное сжатие для улучшения, в последствии, отделения влаги. После этого, воздух поступает во вторичные теплообменный аппарат (ТОА2), где так же происходит теплообмен с забортным воздухом. Продувка первичного и вторичного теплообменных аппаратов (ТОА1 и ТОА2) обеспечивается за счет работы ступени вентилятора (В) турбохолодильника или набегающим потоком воздуха.

Далее воздух поступает в так называемую «петлю», состоящую из двух теплообменных аппаратов (ТОА3 и ТОА4), детандерной ступени турбохолодильника (Т) и влагоотделителя (ВО). В данной «петле» охлаждение рабочего воздуха по горячему тракту происходит за счет того же потока воздуха, после его охлаждения в детандерной ступени турбохолодильника. Рабочий воздух поступает в перегреватель (ТОА3), а затем в конденсатор (ТОА4). В данных теплообменных аппаратах происходит конденсация влаги, которая затем отделяется во влагоотделителе (ВО).

Дабы избежать попадания оставшейся влаги в ступень детандера, а также предотвратить преждевременную конденсацию последней, рабочий воздух поступает в перегреватель по холодному тракту, где подогревается, обеспечивая отвод теплоты от горячего тракта, и поступает в детандерную ступень (Т) турбохолодильного агрегата. В детандерной ступени (Т) происходит расширение и охлаждение рабочего потока воздуха, а совершаемая им при этом работа по валу турбохолодильного агрегата передается на компрессорную ступень (К) и ступень вентилятора (В). После этого основная часть потока воздуха поступает по холодному тракту в конденсатор (ТОА4) обеспечивая тем самым отвод теплоты от горячего тракта. Отведённая влага во влагоотделителе (ВО) поступает на впрыск в продувочный тракт вторичного теплообменника СКВ, увеличивая тем самым эффективность охлаждения в первичном и вторичном теплообмен- никах (ТОА1 и ТОА2).

Система кондиционирования воздуха может работать в двух режимах: охлаждения и обогрева. В режиме обогрева происходит регулирование температуры воздуха на выходе путем отвода горячего потока через регулирующее устройство (РУ) и подачи его на смешивание с потоком, выходящим из детандерной ступени (Т) турбохолодильного агрегата.К преимуществам данной схемы можно отнести небольшое количество агрегатов, необходимых для создания «петли», простоту, надежность, высокую технологичность агрегатов «петли», а так же применения в качестве хладагента в петле рабочего воздуха.

ОПИСАНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ

Для анализа работы СКВ была разработана математическая модель системы в среде Matlab Simulink. В модели производится термодинамический расчет системы. Аналогичные расчеты систем и агрегатов проводились в работах [8–10]. Модель рассчитывает параметры узловых точек СКВ и некоторых параметров её узлов, в частности. Применение данной среды было обусловлено возможностью легкой замены составляющих агрегатов СКВ и тем самым унификации модели под аналогичные схемы. В данной схеме используются математические модели аппаратов с учетом характеристик работы их реальных аналогов. В модели использован метод последовательных итераций для расчета всей схемы в целом. Указанная модель учитывает влияние фазового перехода влаги на температуру рабочего воздуха (в частности, в ней учтено влияние фазового перехода при впрыске влаги в продувочный контур вторичного теплообменного аппарата, в иных теплообменных аппаратах в рабочем тракте и ступени детандера). Работа модели разделена на два блока. Блок № 1 рассчитывает работу первичного и вторичного теплообменных аппаратов, ступень компрессора, а также продувочный контур со ступенью вентилятора. Блок № 2 рассчитывает работу «петли» системы, включающей в себя перегреватель, конденсатор, влагоотделитель и ступень турбины. Каждый блок имеет свои итерационные параметры. Блок № 1 последовательно рассчитывает параметры основных узловых точек и подбирает расход продувочного воздуха. Блок № 2 определяет параметры узловых точек, количество отведенной влаги и мощность, отведённую от детандерной ступени турбохолодильника.

 

Рис. 2. Расчетная схема с точками.

Fig. 2. Design diagram with points.

 

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

АНАЛИЗ РАБОТЫ СИСТЕМЫ

При разработке и производстве СКВ важно иметь подробное представление о работе такой системы. Для каждой СКВ необходимо проводить свой термодинамический расчет и анализ работы такой системы. В большинстве случаев эффективность работы СКВ анализируют по её энергопотреблению (а точнее, в основном по потреблению горячего воздуха). В МГТУ им. Н.Э. Баумана активно используется метод энтропийно-статистического анализа, с помощью которого можно подробно проанализировать работу СКВ на режимах охлаждения. Данный метод позволяет определить узлы, потери в которых имеют наибольшее влияние на эффективность СКВ, а также позволяет оценить перспективу улучшения этих узлов. В данной статье этот метод применен для анализа наземного режима охлаждения самолета, как наиболее нагруженного для СКВ. Методика расчета взята на основе работы [4]. Результаты расчета наземного режима СКВ по математической модели представлены в таблице 1.

 

Таблица 1. Расчетные параметры СКВ

Table 1. Calculated Parameters of ECS

№	Давление, кПа	Температура, К	Расход, кг/ч	Влагосодержание в виде пара, г/кг сух. возд.	Влагосодержание в виде капель, г/кг сух. возд.	Энтальпия, кДж/кг	Энтропия, кДж/кг*К
Параметры Атмосферы
А	101,301	310	1185	18	-	310,303	6,898
Линия рабочего потока
0	304,65	473	1185	18	0	475,463	7,010
1	302,425	473	1185	18	0	475,465	7,012
2	280,602	378,994	1185	18	0	379,681	6,808
4	375,503	424,192	1183,8	18	0	425,504	6,838
6	348,652	314,583	1183,8	14,764	3,236	314,413	6,557
12	338,705	311,004	1183,8	12,522	5,478	310,818	6,553
13	329,388	298,716	1183,8	6,351	11,658	298,425	6,521
16	316,235	299,136	1163,8	6,931	0,816	298,88	6,534
17	315,514	305,103	1163,8	7,681	0,066	304,907	6,554
19	103,193	256,906	1144,1	0,916	6,832	256,887	6,704
20	101,301	275,423	1144,1	4,516	3,232	275,513	6,779
Продувочная линия
7	101,301	299,561	2788	22,45	0	299,796	6,864
8	99,639	348,887	2788	22,45	0	349,498	7,022
9	98,791	389,104	2788	22,45	0	390,154	7,135
10	101,3	396,283	2788	22,45	0	397,425	7,146

 

Расчёт показателей цикла по рабочему тракту осуществлялся по следующим формулам:

Удельная холодопроизводительность цикла:

q_ц = $\left({\mathrm{i}}_{\mathrm{A}}-{\mathrm{i}}_2\right) · \frac{{\mathrm{G}}_2}{{\mathrm{G}}_0}+\left({\mathrm{i}}_2-{\mathrm{i}}_6\right) · \frac{{\mathrm{G}}_6}{{\mathrm{G}}_0}+\left({\mathrm{i}}_{17}-{\mathrm{i}}_{19}\right) · \frac{{\mathrm{G}}_{19}}{{\mathrm{G}}_0}$ = 42,189 $\frac{\text{кДж}}{\text{кг}}$;

Минимальная удельная работа цикла:

l_min = q_ц $· \left(\frac{T_A}{T_A-T_{20}} · ln\left(\frac{T_A}{T_{20}}\right)-1\right)+q_{\text{ц}} · \left(\frac{T_A}{T_{20}-T_{19}} · ln\left(\frac{T_{20}}{T_{19}}\right)-1\right)$= 9,512 $\frac{\text{кДж}}{\text{кг}}$;

Реальная удельная работа двигателя:

l_реал = $\frac{k}{k-1}$$\frac{{\left(\frac{P_0}{P_A}\right)}^{\frac{k-1}{k}}-1}{{\eta }_{\text{дв}}}$ · R_B · T_A   = 160,881$\frac{\text{кДж}}{\text{кг}}$;

Действительное значение холодильного коэффициента

ε = $\frac{q_{\text{ц}}}{l_{\text{реал}}}$ = 0,262.

Удельные работы, необходимые для компенсации производства энтропий:

• для двигателя

$\mathrm{\Delta }{l}_{\text{дв}}=$$\frac{\mathrm{k}}{\mathrm{k}-1}$R · T_A ·$[{\left(\frac{{\mathrm{P}}_0}{{\mathrm{P}}_{\mathrm{A}}}\right)}^{\frac{\mathrm{k}-1}{\mathrm{k}}}-1]$ · $(\frac{1}{{\mathrm{\eta }}_{\text{дв}}}-1)$ = 45,419$\frac{\text{кДж}}{\text{кг}}$;

• для первичного теплообменного аппарата

$∆l_{\text{тоа}1}=$${\mathrm{T}}_{\mathrm{A}} ·\left[\left({\mathrm{S}}_9-{\mathrm{S}}_8\right) · \frac{{\mathrm{G}}_9}{{\mathrm{G}}_0}-\left({\mathrm{S}}_1-{\mathrm{S}}_2\right) · \frac{{\mathrm{G}}_2}{{\mathrm{G}}_0}\right]$= 18,923$\frac{\text{кДж}}{\text{кг}}$;

• для вторичного теплообменного аппарата

$∆l_{\text{тоа2}}=$${\mathrm{T}}_{\mathrm{A}} · \left[\left({\mathrm{S}}_8-{\mathrm{S}}_7\right) · \frac{{\mathrm{G}}_8}{{\mathrm{G}}_0}-\left({\mathrm{S}}_4-{\mathrm{S}}_6\right) · \frac{{\mathrm{G}}_6}{{\mathrm{G}}_0}\right]$= 28,268$\frac{\text{кДж}}{\text{кг}}$;

• для компрессора

$\Delta l_{\text{комп}}=$$\frac{\mathrm{k}}{\mathrm{k}-1} · {\mathrm{R}}_{\text{в}} · {\mathrm{T}}_2 · [{\left(\frac{{\mathrm{P}}_4}{{\mathrm{P}}_2}\right)}^{\frac{\mathrm{k}-1}{\mathrm{k}}}-1] · \left(\frac{1}{{\mathrm{\eta }}_{\text{комп}}}-1\right)$ = 12,306 $\frac{\text{кДж}}{\text{кг}}$;

• для влагоотделителя

$∆l_{\text{влаг}}={\mathrm{T}}_{\mathrm{A}} · \left({\mathrm{S}}_{16}-{\mathrm{S}}_{13}\right) ·$$\frac{{\mathrm{G}}_{16}}{{\mathrm{G}}_0}$= 4,021$\frac{\text{кДж}}{\text{кг}}$;

• для учета теплообмена с окружающей средой

 $∆l_{\text{мним}}=[{\mathrm{T}}_{\mathrm{A}} · \left({\mathrm{S}}_{\mathrm{A}}-{\mathrm{S}}_{20}\right)-\left({\mathrm{i}}_{\mathrm{a}}-{\mathrm{i}}_{20}\right)] ·$$\frac{{\mathrm{G}}_{20}}{{\mathrm{G}}_0}$= 2,025$\frac{\text{кДж}}{\text{кг}}$ ;

• для перегревателя

$∆l_{\text{тоа}3}={\mathrm{T}}_{\mathrm{A}}·$$\left[\left({\mathrm{S}}_6-{\mathrm{S}}_{12}\right) · \frac{{\mathrm{G}}_{12}}{{\mathrm{G}}_0}-\left({\mathrm{S}}_{16}-{\mathrm{S}}_{17}\right) · \frac{{\mathrm{G}}_{17}}{{\mathrm{G}}_0}\right]$= 7,261 $\frac{\text{кДж}}{\text{кг}}$;

• для конденсатора

$∆l_{\text{тоа}4}={\mathrm{T}}_{\mathrm{A}} ·$$\left[\left(S_{20}-S_{19}\right) · \frac{G_{19}}{G_0}-\left(S_{12}-S_{13}\right) · \frac{G_{13}}{G_0}\right]$= 12,341 $\frac{\text{кДж}}{\text{кг}}$;

• для детандера

$\Delta l_{\text{дет}}=$$\frac{k}{k-1}$$· R_{\text{в}} · T_{17}·$$\left[1-{\left(\frac{P_{19}}{P_{17}}\right)}^{\frac{k-1}{k}}\right]$$· \left(1-{\eta }_{\text{турб}}\right)$ = 17,659$\frac{\text{кДж}}{\text{кг}}$;

• для учёта гидравлических потерь в трубопроводе перед первичным теплообменным аппаратом

$∆l_{\text{гидр}}=T_A · \left(S_1-S_0\right) ·$$\frac{G_1}{G_0}$= 0,654$\frac{\text{кДж}}{\text{кг}}$;

Суммарная работа цикла

$∆l_{c\text{ум}}=l_{min}+ ∆l_{\text{тоа}1}+ ∆l_{\text{тоа}2}+\mathrm{\Delta }{l}_{\text{комп}}+∆l_{\text{влаг}}+∆l_{\text{мним}}+\Delta l_{\text{дет}}+∆l_{\text{тоа}3}+∆l_{\text{тоа}4}+∆l_{\text{гидр}}$,

$∆l_{c\text{ум}}=$158,388 $\frac{\text{кДж}}{\text{кг}}$.

Погрешность расчета

$\frac{l_{\text{реал}}-l_{\text{сум}}}{l_{\text{реал}}}$ = 1,25%

Рабочий цикл представлен на T-S диаграмме на рис. 3, а соотношение работ для компенсации прироста энтропии на рис. 4.

 

Рис. 3. Рабочий цикл в T-S диаграмме.

Fig. 3. Duty cycle in the T–S diagram.

 

Рис. 4. Соотношение работ для компенсации прироста энтропии по узлам.

Fig. 4. Ratio of work required to compensate for the increase in entropy across units.

 

ВЫВОДЫ

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

ДОПОЛНИТЕЛЬНО

Вклад авторов. Все авторы внесли существенный вклад в разработку концепции и подготовку статьи, прочли и одобрили финальную версию перед публикацией.

Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов, связанного с подготовкой и публикацией статьи.

Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования и подготовке публикации.

ADDITIONAL INFORMATION

Authors’ contributions. All authors made a substantial contribution to the conceptual development and preparation of this article and read and approved the final version before publication.

Competing interests. The authors declare that they have no competing interests.

Funding source. This study was not supported by external sources of funding.

Об авторах

Игорь Валерьевич Тищенко

Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет); ПАО НПО «Наука»

Email: iv.tischenko@bmstu.ru
ORCID iD: 0000-0001-6094-8723
SPIN-код: 5630-4301

кандидат технических нак, доцент

Россия, 105005, Москва, Лефортовская наб., д. 1; 105005, 125124, Москва, 3-я ул. Ямского поля, вл. 2

Сергей Александрович Абалакин

ПАО НПО «Наука»

Email: sserviuss@gmail.com
ORCID iD: 0009-0001-9193-1147
SPIN-код: 7580-9545

кандидат технических нак

Россия, 125124, Москва, 3-я ул. Ямского поля, вл. 2

Владислав Иванович Меркулов

Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет); ПАО НПО «Наука»

Email: merkulov@bmstu.ru
ORCID iD: 0000-0002-3455-3709
SPIN-код: 8409-8754

доктор техничеких наук, профессор

Россия, 105005, Москва, Лефортовская наб., д. 1; 105005, 125124, Москва, 3-я ул. Ямского поля, вл. 2

Антон Андреевич Жаров

Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет)

Email: zharov_a@bmstu.ru
ORCID iD: 0000-0001-9945-0850
SPIN-код: 8581-1809

кандидат технических нак

Россия, 105005, Москва, 2-я Бауманская ул., д. 5

Артем Витальевич Борисенко

Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет)

Автор, ответственный за переписку.
Email: borart@bmstu.ru
ORCID iD: 0000-0002-4818-3702
SPIN-код: 2859-5006

кандидат технических нак

Россия, 105005, Москва, 2-я Бауманская ул., д. 5

Список литературы

Дополнительные файлы