CC BY
90
УДК 629.7.036
ТЕРМОДИНАМИЧЕСКИЙ РАСЧЕТ КЭУ НА БАЗЕ ГТУ И ТВЕРДООКСИДНЫХ ТОПЛИВНЫХ ЭЛЕМЕНТОВ SOFC
© А.А. Лоскутников1, И.М. Горюнов2, Н.С. Сенюшкин3, Р.Ф. Султанов4
Уфимский государственный авиационный технический университет, 450000, Россия, г. Уфа, ул. Карла Маркса, 12.
Рассмотрен алгоритм расчета твердооксидных топливных элементов (ТОТЭ) для расширения возможностей программного комплекса DVIGwT. Проведена проверка адекватности результатов расчетных исследований ТОТЭ с экспериментальными данными. Ил. 9. Библиогр. 12 назв.
Ключевые слова: твердооксидные топливные элементы; комбинированная энергоустановка; пароводяная конверсия; структурно-параметрический синтез и анализ; математическая модель; газотурбинный двигатель.
THERMODYNAMIC CALCULATION OF COMBINED POWER PLANT BASED ON GAS TURBINE UNIT AND SOLID OXIDE FUEL CELLS
A.A. Loskutnikov, I.M. Goryunov, N.S. Senyushkin, R.F. Sultanov
Ufa State Aviation Technical University, 12 Karl Marx St., Ufa, 450000, Russia.
The article examines an algorithm to calculate solid oxide fuel cells (SOFC) in order to extend the features of the software package DVIGwT. The validation of the results of SOFC computational modeling with experimental data is performed.
9 figures. 12 sources.
Key words: solid oxide fuel cells; combined power plant; steam water conversion; structural and parametric synthesis and analysis; mathematical model; gas turbine engine.
Увеличивающийся дефицит электрической и тепловой энергии покрывается внедрением газотурбинных установок (ГТУ) на базе конвертированных авиационных газотурбинных двигателей (ГТД). При этом совершенствуются старые и создаются новые перспективные методы конвертирования для уменьшения потерь тепла с уходящими газами, с повышением эффективного КПД установок [1, 2].
Конвертируемые ГТД чаще всего создаются на базе авиационных двигателей, принадлежащих к предыдущим поколениям, следовательно, обладающих умеренными параметры рабочего процесса, вы-
сокой температурой уходящих газов и достаточно низким эффективным КПД (порядка 25-30%) [3].
Значительно повысить эффективность энергоустановок на органическом топливе возможно с применением комбинированных энергоустановок (КЭУ) на основе ГТУ и топливных элементов (ТЭ). Отличительной способностью ТЭ является возможность производить как электроэнергию, так и теплоту, которую можно использовать в циклах ГТУ, паротурбинных установок (ПТУ) или парогазовых установок (ПГУ).
Над созданием ТЭ и КЭУ на их базе работают основные промышленные компании как в США, так и в
1Лоскутников Александр Александрович, кандидат технических наук, доцент кафедры авиационной теплотехники и теплоэнергетики, тел.: (347) 2737954, e-mail: [email protected]
Loskutnikov Alexander, Candidate of technical sciences, Associate Professor of the Department of Aviation Combustion and Heat-Power Engineering, tel.: (347) 2737954, e-mail: [email protected]
2Горюнов Иван Михайлович, доктор технических наук, профессор кафедры, тел.: (347) 2737954, e-mail: [email protected] Goryunov Ivan, Doctor of technical sciences, Professor of the Department, tel.: (347) 2737954, e-mail: [email protected] 3Сенюшкин Николай Сергеевич, кандидат технических наук, доцент кафедры авиационной теплотехники и теплоэнергетики, тел.: (347) 2737954,e-mail: [email protected]
Senyushkin Nikolai, Candidate of technical sciences, Associate Professor of the Department of Aviation Combustion and Heat-Power Engineering, tel.: (347) 2737954, e-mail: [email protected]
4Султанов Рузиль Фаилевич, аспирант, младший научный сотрудник, тел.: 89373566519, e-mail: [email protected] Sultanov Ruzil, Postgraduate, Junior Researcher, tel.: 89373566519, e-mail: [email protected]
136
ВЕСТНИК ИрГТУ №7 (78) 2013
России [1, 2]. По данным, представленным в работах [1, 2, 4], КПД по выработке электроэнергии на КЭУ, созданных на базе ГТУ и ТЭ, может достигать более 70%. Диапазон мощностей таких КЭУ рассматривается от 200 кВт до 10 МВт.
Существуют три основные схемы КЭУ, созданные на базе ГТД и ТЭ [1, 2]:
- схема 1 - КЭУ, где ТЭ замещает камеру сгорания ГТУ;
- схема 2 - КЭУ на основе ГТУ и ТЭ с внутренней конверсией природного газа, работающая при атмосферном давлении;
- схема 3 - КЭУ на основе ТЭ, работающая на продуктах газификации углей (мощность достигает 300 и более МВт).
Рассмотрим устройство топливных элементов. ТЭ - это объект, в котором химическая энергия топлива (восстановителя) и окислителя, непрерывно и раздельно подводимых к электродам, непосредственно превращается в электрическую энергию [1, 2, 5].
В твердооксидном ТЭ (ТОТЭ), по международной классификации SOFC, электрическая проводимость электролитов, обычно на основе оксидов циркония и иттрия, становится приемлемой при температуре выше 1173 К, следовательно, рабочая температура
ются универсальными и предоставляют большие функциональные возможности, обеспечивают естественную среду для создания и модификации модели.
В перечисленных программных комплексах отсутствует модуль для расчета ТЭ, что делает невозможным выполнение системных исследований рабочих процессов ТЭ, а также созданных на их основе КЭУ.
Анализ особенностей работы ТЭ и систем математического моделирования определяет необходимость разработки алгоритма расчета ТОТЭ, который позволит достаточно точно и оперативно воспроизводить различные схемы и режимы работы КЭУ, созданных на базе ГТУ и ТОТЭ с учетом основных значимых факторов, а также позволит осуществлять учет влияния изменения физических свойств рабочего тела на характеристики элементов КЭУ.
Рассмотрим основные процессы, происходящие в твердооксидном топливном элементе (рис. 1). Для системы топливоподготовки разумно использовать пароводяную конверсию (ПВК) углеводородного топлива (метана) в среде водяного пара, в результате которой выделяется водород и монооксид углерода. В результате протекания основных электрохимических реакций окисления водорода и СО выделяются водяной пар и СО2.
Окислитель
Горючее
Водяной пар
Пароводяная Основные
конверсия реакции
горючего в ТЭ
Продукты реакции ТЭ
Рис. 1. Схема расчета ТОТЭ
ТОТЭ лежит в диапазоне 1173-1273 К. Основным видом топлива в ТОТЭ служит водород, а окислителем, как правило, - кислород воздуха.
С повышением температуры в ТОТЭ увеличивается электрическая проводимость электролита и полупроводниковых межэлементных соединений, поляризация уменьшается, что приводит к улучшению вольт-амперной характеристики (ВАХ) [1]. Благодаря этим параметрам растет интерес к ТОТЭ со стороны энергетики.
Для выполнения расчетных исследований повышения эффективности схем КЭУ необходим математический комплекс расчета ТЭ. Остро стоит вопрос о создании эффективного, надежного и простого в эксплуатации модуля, способного рассчитывать параметры ТЭ.
В настоящее время наиболее известные математические модели и программные средства моделирования ГТД [6, 7] - это комплексные математические модели ЦИАМ, ГРАД, DVIGW7, GasTurb, GSP, GECAT, Ebsilon Professional, United cycle, ОГРА, АСТРА и др. Некоторые программы, например, GasTurb, АСТРА, ОГРА, предназначены для моделирования определённых схем ГТД, другие (ГРАД, комплексные ММ ЦИАМ, DVIGwТ GSP, Ebsilon Professional и др.) явля-
Для осуществления пароводяной конверсии необходимо подводить теплоту к стенкам камеры, в которой происходит ПВК. Энергия химической реакции окисления топлива непосредственно переходит в электрическую и в тепловую энергии.
Для проведения расчета и исследования КЭУ различных схем на базе ГТД и ТОТЭ создан программный модуль SOFC [8].
Оценка работоспособности и адекватности разработанного алгоритма, реализованного в программном модуле БОРО в системе DVIGwT, выполнена на примерах моделирования ячеек ТОТЭ как планарной конструкции [9], так и трубчатой конструкции в составе электрохимического генератора (ЭХГ) мощностью 1 кВт, описанного в работе [10].
Отличия в оценках экспериментальных параметров и параметров, полученных расчетом в модуле БОРО, не превышают 3%.
Выполнено моделирование ЭХГ модуля ТОТЭ для макета энергоустановки (ЭУ) мощностью 1 кВт (максимальное отклонение рассчитанных характеристик не превышает « 3,35% в сравнении с экспериментальными данными). Полученное расхождение экспериментальных и рассчитанных величин объясняется тем, что в программном продукте не удается учесть
влияние катализатора на внутренние процессы в ТОТЭ. Естественно, и экспериментальные данные также получены с определённой погрешностью измеренных параметров.
На основании рассмотрения преимуществ авиационных ГТД, в результате которых они могут быть конвертированы в ЭУ различного назначения, с учетом необходимости повышения их эффективности проведены расчетные исследования повышения эффективности ГТЭ-10/95 [11] путем модернизации в КЭУ следующих схем:
- КЭУ с ЭХГ на базе ТОТЭ, работающем на отборе воздуха за компрессором низкого давления (КНД);
- КЭУ с замененной камерой сгорания (КС) в ГТЭ-10/95 на ЭХГ;
- КЭУ на базе каскада низкого давления (НД) ГТЭ-10/95 и ЭХГ на базе ТОТЭ, питающегося воздухом от
КНД (0,3 МПа).
Рассмотрим способ повышения эффективности ГТЭ-10/95 заменой камеры сгорания на ЭХГ, состоящий из ТОТЭ. Принципиальная схема такой модернизированной КЭУ приведена на рис. 2.
Топливо подается только в ЭХГ. Температуру продуктов реакции ЭХГ, которая частично понижена за счет подогрева входящих в ЭХГ воздуха и топлива, необходимо повысить до температуры газов на выходе из камеры сгорания в ГТЭ на расчетном режиме. Суммарные затраты топлива в ЭХГ и камере дожигания позволяют существенно повысить мощность всей КЭУ.
На основании принципиальной схемы разработана математическая модель, структурная схема которой приведена на рис. 3.
Рис. 2. Принципиальная схема КЭУ на базе ГТЭ-10/95, в которой камера сгорания заменена на ЭХГ на базе ТОТЭ (схема 1): 1 - комплексное воздухоочистительное устройство (КВОУ); 2 - КНД; 3 - компрессор высокого давления (КВД); 4 - ЭХГ; 5 - турбина высокого давления (ТВД); 6 - турбина низкого давления (ТНД); 7 - свободная турбина (СТ); 8 - электрогенератор; 9 - КУ; 10 - насос питательной воды; 11 - бойлер; 12 - насос сетевой воды; 13 - потребитель тепловой энергии; 14 - выхлопное устройство; 15 - камера дожигания; э - экономайзер;
и - испаритель; пп - пароперегреватель
Структурная схема модели состоит из газотурбинного привода, модернизированного ЭХГ на базе ТОТЭ, элементов системы генерации водяного пара: котла-утилизатора, состоящего из экономайзера, испарителя, пароперегревателя и др. Устанавливаются рекурсивные связи, обеспечивающие передачу параметров от последующих элементов к предыдущим элементам структурной схемы (показаны пунктирными линиями).
Расходы направляемых в ЭХГ воздуха, топлива и пара были смоделированы по расходу воздуха в ГТЭ-10/95 на базе работ [1, 12]. Таким же образом, путем аналогии, по расходу воздуха в ГТУ подобрано количество параллельно и последовательно соединенных ячеек ТОТЭ.
Температура ПВК и температура основной реакции в ЭХГ равны [1, 12, 13] и поддерживаются за счет камеры дожигания. Значения данных параметров взяты с учетом соответствующих значений из работы [13], а также рекомендаций по ТОТЭ, изложенных в [1, 13]. Поддерживаемые в процессе расчета схемы температурных напоров на холодных концах различных теп-лообменных аппаратов выбраны по рекомендациям [12].
Для удобства анализа результаты расчетов базовой ГТЭ-10/95 и КЭУ на базе ГТЭ и ЭХГ из ТОТЭ (далее по тексту схема 1) на номинальном режиме приведены в таблице.
Схема с замененной камерой сгорания на ЭХГ в ГТЭ-10/95 требует значительных переделок, однако, экономический эффект от такого вложения компенсируется повышением электрической суммарной мощности КЭУ до 20630,4 кВт (в 2,56 раза больше мощности базовой ГТЭ) и повышением КПД по выработке электроэнергии до 30,82% (на 8,19% абсолютных по отношению к базовой ГТЭ).
ТОТЭ и ЭХГ на их основе сложно поддаются регулировке, так как работают в основном на изобарно-изотермическом режиме. Изменяемые параметры ГТУ при дросселировании КЭУ с ЭХГ непосредственно влияют на протекание основной химической реакции внутри ТОТЭ [13]. В ТЭ параметры воздуха для окисления и пара для ПВК изменяются на переменных режимах работы ГТУ в составе КЭУ.
По результатам расчетов построены дроссельные характеристики исходной ГТЭ и схемы (схема 1) КЭУ (при условии обеспечения регулирования ЭХГ в составе ГТУ).
Проанализированы зависимости расхода топлива, электрической мощности, , расхода воздуха на входе в КНД, приведенной частоты вращения КВД, вырабатываемой теплоты потребителю, п*кнд, п*квд, температуры и давления воздуха на выходе из КВД, температуры газов на выходе из основной камеры или камеры дожигания, температуры и давления газов на выходе из СТ, плотности тока ЭХГ, плотности мощности, температуры газов на выходе из ЭХГ от приведенной частоты вращения ротора КНД (лпр кнд).
Для представленных схем ЭУ параметры получе-
ны в диапазоне приведенных частот вращения ротора КНД, соответствующих мощности от 8 до 20 МВт.
Суммарная электрическая мощность базовой ГТУ и схемы 1 с падением лпр кнд уменьшается (рис. 4). При лпр кнд = 9627 об/мин мощность базовой ГТЭ равна 8000 кВт, в схеме 1 электрическая мощность возросла в 2,58 раза по сравнению с мощностью базовой ГТЭ.
Результаты расчета в системе DVIGwT базовой ГТЭ-10/95 и схемы 1
Параметр ГТЭ-10/95 Схема 1
N, кВт 8000 8000
Пнд> об/мин 9636,25 9621,64
Пвд ' 0б/мин 10385,70 10350,70
GB, кг/с 55,27 54,98
* Анд> кПа 280,7 280,8
ТКнд > К 403,35 403,37
Т*вд. К 550,95 550,59
Р*с > кПа 660,3 659,8
* Ткс , К 1116,45 1060,62
GT, кг/с 0,718 0,390
^кс 4,361 2,289
ТТвд, К 969,20 924,56
кг/с 55,07 56,40
ТТНд, К 866,18 827,98
Р*т > кПа 105,3 105,3
Т*т, к 733,05 702,37
ЫЭХГ, кВт - 12630,4
Вт ПЭХГ' 2 м2 - 870,15
. А гЭХГ , ~ м - 2088,46
* Тосн р, К - 1293,02
ТП ВК, К - 1293,15
I, А - 2631,46
U, В - 0,417
GTЭХГ' кг/с - 0,971
Q, кВт 21433,2 17522,0
КЭУ' кВТ 8000 20630,4
, %** 22,63 30,82
Примечание: подчеркнутые параметры даны для камеры дожигания, так как камера сгорания в схеме отсутствует; ** - параметр рассчитан по отношению электрической мощности к теплоте сжигания топлива.
кВт
25000 -|-
20000
15000
10000
5000
0 и------
7000 7500 8000 8500 9000 9500 10000
ппрКНД
ПЭ 10/95 КЭУ с ЭХГ вместо ОКС
Рис. 4. Зависимость общей электрической мощности ЭУ от ппр кнд
Такое резкое увеличение электрической мощности в схеме 1 объясняется более высокими термодинамическими параметрами входных потоков веществ в ЭХГ, что повлекло за собой резкое увеличение мощности ЭХГ и, следовательно, КЭУ в целом. С уменьшением лпр кнд суммарные расходы топлива для камеры сгорания в базовой ГТЭ для камеры дожигания в схеме 1 уменьшаются. При лпр КНд = 9627 об/мин, соответствующем номинальному режиму работы базовой ГТЭ, расход топлива в базовой схеме равен 0,718 кг/с, в схеме 1 - 1,36 кг/с, что на 52,8% больше расхода базовой схемы.
Анализ изменения электрического КПД рассматриваемых ЭУ (рис. 5) показ2ы5в0а0е0т, что при снижении
КПДэ
40,00
35,00 30,00 25,00 20,00 15,00 10,00
7000 7500 8000 8500 9000 9500 10000
ппрКНД
нд ПЭ 10/95 КЭУ с ЭХГ вмасто ОКС
Рсс. 5. Зависимость КПД ЭУ от пПр кпд
ппр пНд в базовор ГТЭ пэ уменьшается, а в схеме 1 цэ возрастает и при лпр КНд = 7100 об/мин достигает 35%, в то время эак у базовой ГТЭ пэ = 12%. Это объясняется тем, что в схеме 1 электрическая мощность снижается не так интенсивно, как уменьшается суммарный расход топлива на КЭУ при дроссельных режимах.
Повышение электрического КПД с уменьшением ппр КНд есть следствие распределения электрической мощности между ГТЭ и ЭХГ с учетом того, что ЭХГ работает практически на стационарном режиме.
Расход воздуха во входное устройство уменьшается со снижением ппр КНд во всех рассматриваемых схемах ЭУ.
Приведенная частота вращения КВД (ппр квд) при уменьшении ппр КНд снижается как для базовой ГТЭ, так и для схемы 1. Особенностью является наличие максимума при ппр КНд = 8300 об/мин.
Для рассматриваемых ЭУ переданная потребителю теплота Опот снижается с уменьшением ппр КНд (рис. 6). При ппр КНд = 7100 об/мин Опот в базовой ГТЭ равна 11000 кВт, в схеме 1 - 9300 кВт. Разница между значениями теплоты в исследуемых ЭУ во всем диапазоне ппр КНд сохраняется. Объясняется это тем, что в схеме 1 большая часть теплоты уходящих газов расходуется на генерацию пара для ПВК, а в базовой ГТЭ - на нужды потребителя. Уменьшение п*КНд для исследуемых ЭУ происходит с понижением ппр КНд. Значения п*КНд для схемы 2 и базовой совпадают во всем исследуемом диапазоне ппр КНд.
Зависимости изменения п*квд для двух ЭУ имеют прямопропорциональную зависимость от ппр КНд. Кривая зависимости пКВд для схемы 1 пересекает кривую зависимости п*квд для базовой ГТЭ при ппр КНд = 9700
об/мин.
Рассмотрим зависимость температуры продуктов сгорания на выходе из камеры сгорания (для базовой схемы ГТЭ) и для камеры дожигания схемы 1 от ппр КНд. В КЭУ исследуемая температура несколько ниже базовой, так как камера дожигания установлена для обеспечения поддержания пропускной способности ТВД. Параметры продуктов сгорания за СТ необходимы для подбора режимов работы котла-утилизатора. Давления за СТ для исследуемых схем равны, что было учтено при проведении расчетов. В исследуемых ЭУ значения температуры газов за СТ (Т*СТ) уменьшаются с уменьшением ппр КНд. Значения исследуемой температуры схемы 1 ниже значений базовой ГТЭ и при ппр кнд = 9600 об/мин достигают: для схемы 1 - 700 К, для базовой ГТЭ - 733 К. При рассмотрении характеристик ЭХГ на дроссельных режимах в составе КЭУ необходимо отметить тенденцию уменьшения плотности тока и плотности мощности (рис. 7, 8) с уменьшением приведенной частоты вращения КНД.
Q, кВт 24000
19000
14000
9000
4000
7000 7500 8000 8500 9000 9500 10000
ппр КНД
0Т" 1 //95 —■— КЭУ к "Х0 вместо СКК Рис. 6. Зависимость производимой теплоты от Ппр кнд
¡, А/м2
2500
2000
1500
1000
500
7000 7500 8000 8500 9000 9500 10000
ппр КНД
Рис. 7. Зависимость плотности тока ЭХГ от ппр кнд (схема 1)
Рис. 8. Зависимость плотности мощности ЭХГ от пПр кнд (схема 1)
Рис. 9. Зависимость температуры газов на выходе из ЭХГ от Ппр кнд (схема 1)
С уменьшением лпр кнд температура на выходе из ЭХГ в КЭУ практически не изменяется (порядка 15001520 К) (рис. 9).
Для обработки результатов моделирования были применены методы системного подхода и сравнения с экспериментальными данными. По итогам проведенных исследований можно сформулировать следующие выводы:
1. КЭУ с замененной камерой сгорания на ЭХГ в ГТЭ-10/95 требует значительных переделок базовой ГТЭ, но при этом достигается повышение электрической суммарной мощности КЭУ до 20630,4 кВт при КПД по выработке электроэнергии 30,82% на номинальном режиме.
2. Анализ дроссельных характеристик показывает, что со снижением приведенной частоты вращения КНД в базовой ГТЭ цэ уменьшается, а в КЭУ цэ возрастает и при лпр кнд = 7100 об/мин достигает 35%.
3. Электрическая мощность в КЭУ падает при
снижении приведенной частоты вращения КНД с 9621 до 7146, при этом мощность всей ЭУ уменьшается с 20630 до 12436 кВт (на 39,72% от номинала), в то время как в базовой схеме при таком же диапазоне изменения приведенной частоты вращения КНД электрическая мощность снижается с 8000 до 2000 кВт (75% от номинала).
Таким образом, результаты моделирования КЭУ на базе ГТД и ТОТЭ, приведенные на примере ГТЭ-10/95 и КЭУ с замененной камерой сгорания на ЭХГ из ТОТЭ при различных схемах совмещения, подтверждают функциональные возможности разработанного авторами программного комплекса. На его основе возможно создание новых схем КЭУ, выполнение многовариантных предпроектных расчетов, включая расчеты характеристик КЭУ, что в свою очередь сократит сроки и затраты на их разработку.
Работа выполнена при финансовой поддержке Минобрнауки РФ._
Библиографический список
1. Коровин Н.В. Топливные элементы и электрохимические энергоустановки. М.: Изд-во МЭИ, 2005. 280 с.
2. Липилин А. С. Состояние и будущее индивидуальной энергетики // Альтернативная энергетика и экология. 2009. № 9. С. 139-152.
3. Цанев С.В., Буров В.Д., Ремезов А. Н. Газотурбинные и парогазовые установки тепловых электростанций: учеб. пособие для вузов / под ред. С.В. Цанева. М.: Изд-во МЭИ, 2002. 584 с.
4. Захаренков Е.А. Исследование и оптимизация схем и параметров гибридных электростанций на основе топливных элементов и газотурбинных установок: автореф. дис. ... канд. техн. наук: 05.14.01. М.: МЭИ (ТУ), 2009. 20 с.
5. Разработка модуля термодинамического расчета твердо-оксидных топливных элементов БОРО / А.А. Лоскутников [и др.] // Вестник Воронежского государственного технического университета. 2010. Т. 6. № 10. С. 186-190.
6. Горюнов И.М. Структурно-параметрический синтез и анализ авиационных гТД и энергетических установок на их основе: дис. ... д-ра технических наук: 05.07.05. Уфа: УГАТУ, 2007. 330 с.
7. Горюнов И.М. Термогазодинамические расчеты ГТД и теплоэнергетических установок с использованием системы DVIGwT // Вестник УГАТУ. 2006. Т. 7. № 1. С. 61 -70.
8. Лоскутников А.А., Горюнов И.М., Бакиров Ф.Г. Свидетельство в Российском агентстве по патентам и товарным знакам (РОСПАТЕНТ). № 2009613945 от 24.07.20099. Исследование электрических характеристик планарного твердоок-сидного топливного элемента / Д.В. Белов [и др.] // Твердо-оксидные топливные элементы: сб. науч.-техн. ст. Снежинск: Изд-во РФЯЦ - ВНИИТФ, 2003. С. 334-339.
10. Электрохимический модуль твердооксидного топливного элемента для макета энергоустановки мощностью 1 кВт /
B.Ф. Чухарев [и др.] // Твердооксидные топливные элементы: сб. науч.-техн. ст. Снежинск: Изд-во РФЯЦ - ВНИИТФ, 2003.
C. 352-357.
11. Газотурбинная энергоустановка ГТЭ-10/953-002. Руководство по эксплуатации и техническому обслуживанию. ГТП 953.00.0.0000РЭ. Уфа: ФГУП «НПП»Мотор», 2004. 458 с.
12. High Efficiency PSOFC / ATS Gas Turbine Power System. Final report. Nat. Energy Tech. Laboratory. Siemens Westing-house Power Corp, 2001. 180 p.
13. Борисов В.Н., Лукашенко И.Г., Ахлюстин М.А. Введение в термодинамику топливного элемента // Твердооксидные топливные элементы: сб. науч.-техн. ст. Снежинск: Изд-во РФЯЦ - ВНИИТФ, 2003. С. 9-15.
