固体氧化物燃料电池/燃气轮机混合动力系统建模仿真研究进展
0 引 言
国家发展改革委在《“十四五”可再生能源发展规划》中提出,实现能源低碳化、清洁化、高效化利用是现阶段主要发展目标[1]。分布式发电具有投资少、安装地点灵活、建设周期短、能源利用率高及环境污染小等优势[2],同时可结合余热回收利用进一步提高能源利用效率[3]。常见的分布式发电装置主要有微型燃气轮机(Micro Gas Turbine,MGT)、内燃机、燃料电池(Fuel Cell,FC)等。其中,燃料电池是一种直接将燃料中化学能转化为电能的发电装置,由于发电效率不受卡诺循环效率限制,FC与传统热机相比能量利用效率可提高10%~20%,且污染物排放明显低于传统火力发电[4],因此燃料电池分布式发电技术是适应未来能源低碳化、清洁化、高效化发展趋势的重要方向。根据电解质不同,燃料电池可分为质子交换膜燃料电池(Proton Exchange Membrane Fuel Cell,PEMFC)、固体氧化物燃料电池、熔融碳酸盐燃料电池(Molten Carbonate Fuel Cell,MCFC)、磷酸燃料电池(Phosphoric Acid Fuel Cell,PAFC)和碱性燃料电池(Alkaline Fuel Cell,AFC)5类。SOFC、MCFC、AFC反应产物在阳极生成,PAFC、PEMFC反应产物在阴极生成。同时,PEMFC、PAFC、AFC均需使用H2作为燃料,如果使用碳氢化合物如CH4、乙醇等作为燃料则需在系统中额外引入外部燃料重整器,不仅增加系统复杂性和成本,还会降低整体效率[5]。与之相比,SOFC、MCFC作为高温燃料电池,电池内部可发生燃料重整,可直接使用碳氢化合物作为燃料,燃料灵活性更高。其中,SOFC可使用H2、CO、CH4、NH3等作为燃料,主要用于小型分布式发电站[6]。考虑到SOFC尾气温度高,具有高品位废热,因此可与采用有机朗肯循环(ORC)[7]、布雷顿循环[8-9]或底层蒸汽循环[10]等热机耦合构建热电联产(CHP)装置,用来额外发电或供热。研究表明,SOFC混合发电系统最大效率可达90%[11]。
常见热机有燃气轮机、蒸汽机等。燃气轮机是21世纪动力设备的核心[12]。SOFC/GT混合动力系统中,GT可回收利用SOFC尾气中的废热从而提高系统效率。关于SOFC/GT混合动力系统研究较多,考虑到SOFC/GT测试平台的建造成本,数学建模仿真为SOFC/GT混合系统提供了一种经济高效的方法。学者已基于gPROMs、Aspen Plus、Matlab/Simulink等商业化软件,针对SOFC/GT系统开展了前瞻性热力学建模和仿真分析工作,从而优化混合系统性能。但目前仍局限于利用已有商业化软件对SOFC/GT系统开展模块化建模工作,缺少相关开源代码,而基于商业化软件搭建的模型功能有限、针对性受限。目前亟需针对SOFC/GT系统建模和仿真技术开展进一步研究。
针对SOFC/GT混合动力系统的建模与仿真研究,笔者介绍了SOFC、GT基本结构、工作原理与数学建模研究现状;论述了SOFC/GT混合发电系统布置方式及SOFC/GT数学仿真模型;综述了现阶段SOFC/GT相关研究工作、SOFC/GT混合发电系统示范工程与相关建模仿真工作,其中SOFC/GT混合动力系统建模仿真工作按照SOFC/GT建模方法展开;最后对未来SOFC/GT混合动力系统建模仿真工作发展方向进行展望。
1 固体氧化物燃料电池及燃气轮机
1.1 固体氧化物燃料电池简介
1.1.1 SOFC基本结构及工作原理
SOFC作为高温燃料电池,主要部件为阳极、阴极、电解质构成的膜电极结构(Membrane-electrode Assembly,MEA),工作原理[13]如图1所示。其中,O2-为电化学反应中的离子导体,燃料(H2、CO)在阳极发生氧化反应并向外电路释放电子,氧化剂(O2)在阴极接受来自外部电路的电子并进行还原反应[14]。
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/bbac716bd42cee7b48fadfc1be519f60.jpg)
图1 SOFC工作原理[13]
Fig.1 Schematic principle of SOFC[13]
根据物理结构不同,SOFC可分为管式(TSOFC)和平板式(PSOFC)2种[15]。Westinghouse公司于1980年成功研制出具有高耐久性的TSOFC,不存在高温密封问题,但输出功率较低且成本较高,与之相比PSOFC成本相对低廉,且输出功率密度和性能较好,是目前主流商用化SOFC类型。
由于SOFC工作温度高,阳极可发生燃料内重整反应,具有较高的燃料灵活性,CH4[16-17]、生物燃料[18]等烃类燃料均可经燃料重整反应后作为SOFC反应物[19],Bloom Energy公司推出的Bloom Box固体氧化物燃料电池微型电站几乎可以使用所有烃类燃料(如乙醇、生物机油、沼气、天然气)发电,且在使用石油、煤气层等化石能源发电时,碳排量仅为传统火电厂的40%[20]。由于与H2相比,CH4价格低廉且更易管理,SOFC/GT混合系统常采用CH4作为燃料[15]。CH4作为燃料的SOFC/GT系统首先需进行燃料重整,根据重整反应发生位置不同,SOFC系统可分为内部燃料重整(Internal Reforming,IR)和引入外部燃料重整器2种。IRSOFC系统由于不需引入额外的重整器,成本低,且燃料重整反应作为吸热过程可对SOFC电堆进行冷却,因此更常用[21-24],较多SOFC/GT混合动力系统建模仿真研究工作涉及SOFC内部重整过程建模[25-29]。在CH4作为燃料的IRSOFC系统中,阳极除发生H2电化学反应外,还发生CH4蒸汽重整反应(Methane Steam Reforming Reaction,MSRR)和水煤气变换反应(Water Gas Shift Reaction,WGSR),电池内部发生的反应方程见表1。
表1 CH4作为燃料的SOFC系统内部反应方程式[17]
Table 1 Internal reaction equation of SOFC system with CH4 as fuel[17]
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/7b140aab6283c77a7abc0b095b9d3b8e.jpg)
1.1.2 SOFC建模工作现状
SOFC建模工作主要围绕多孔电极中的气体流动、传热、多组分传质、化学和电化学反应展开,并对电池极化曲线(Ⅴ-Ⅰ曲线)、稳态运行状况、瞬态行为进行分析。按照建模对象不同,SOFC模型可分为SOFC电极机理模型、SOFC单电池模型、SOFC电堆模型、系统级SOFC模型等[14]。
近年来,计算机技术进步使SOFC仿真可实现从量子(亚原子)到原子再到连续尺度的多尺度建模[30],建模方法分为微观尺度方法(如密度泛函理论、分子动力学、量子化学等)、中尺度方法(玻尔兹曼法、格点法、蒙特卡罗法)、宏观尺度方法(如有限元法、有限差分法、有限体积法)[14]。
对于宏观SOFC建模而言,按照建模维度划分可将SOFC分为零维(0D)、一维(1D)、二维(2D)、三维(3D)[31],按照建模维度,列举了一些已有的SOFC仿真模型[25-28,30-33](表2)。以3D建模为例,SOFC三维模型主要利用已有商业CFD软件(如Fluent、COMSOL、CFD-ACE 、gPROMs、STAR-CD[32]等)对燃料电池全空间长度进行建模,进而实现CFD仿真计算。如BESSETTE等[26]较早开展TSOFC三维建模工作,对SOFC建立单管稳态三维电化学、热力学数学模型;ZHU等[33]利用STAR-CD软件开展PSOFC电堆稳态建模,研究了SOFC电堆中燃料浓度、温度、局部电流浓度等参数的分布。RAMREZ-MINGUELA等[27]对SOFC采用COMSOL软件建立了三维CFD模型,分析了H2、CO的电化学反应,燃料重整反应及水煤气变换反应对SOFC系统熵产率和性能的影响,突出燃料电池内部发生的主要不可逆性。
表2 SOFC建模举例
Table 2 Examples of SOFC models
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/058764486a6d84aa59d1e1555ad9910c.jpg)
高维度模型(如2D、3D模型)需要使用CFD软件,并基于SOFC中详细的传热过程及质量传输方程展开计算。高维度模型的优点是可用于捕捉局部物质浓度、局部电流密度、SOFC热梯度分布等重要参数[13,31],但计算量较大,成本较高,因此SOFC/GT混合动力系统多采用SOFC 0D、1D模型进行热力学建模。0D模型是将SOFC单电池或电堆视为集总模型,即不考虑相关参数(物质浓度、电流密度、电压等)沿空间坐标的变化;1D模型仅考虑相关参数沿主要维度的参数分布;因此0D、1D CFD模型计算量较小,适合于SOFC系统性能优化,尤其是控制方面的研究。
同时,不同维度SOFC模型按照是否考虑各参数随时间变化分为稳态模型及瞬态模型2种。在SOFC/GT混合动力系统中引入SOFC动态模型工作量大且计算成本较高,因此现有模型中SOFC稳态模型较多[13]。相关研究常基于Matlab、Aspen Plus等软件建立系统级SOFC稳态模型,即不考虑温度、物质浓度等参数随时间变化。因此第3.2节主要基于SOFC的0D稳态模型对SOFC/GT混合动力系统数学模型展开介绍。
1.2 燃气轮机简介
1.2.1 燃气轮机基本结构与工作原理
燃气轮机是目前技术较成熟的分布式发电装置,系统循环以布雷顿循环为基础,包括压缩机、燃烧室、透平、回热器、发电机等部件[34],常见单轴式燃气轮机结构如图2所示。空气经压缩机压缩后送至燃烧室燃烧,废气送至涡轮膨胀驱动发电机和压缩机,热交换器可预热进入燃烧室的气体[35]。换热器部件有助于进一步提高燃气轮机效率,实现节能减排。常见燃料包括天然气、汽油、柴油等。
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/1fd07034acc45c93960c714de0d56c9d.jpg)
图2 单轴式燃气轮机(含换热器)
Fig.2 Single-shaft gas turbine (including heat exchanger)
燃气轮机系统可分为重型燃气轮机发电系统、小型燃气轮机发电系统、微型燃气轮机发电系统。其中,MGT常见功率为25~300 kW,具有结构紧凑、维修成本低、振动小等优点;换热器部件有助于进一步提高MGT效率,实现节能减排。基于以上特点,MGT应用范围广,适用于分布式发电、CHP系统和FC/MGT联合发电系统。
1.2.2 燃气轮机建模现状
燃气轮机整体性能主要取决于压缩机、燃烧室、透平等重要部件[36],因此建模工作主要围绕GT系统中的压缩机、透平、燃烧室等展开。常见的燃气轮机建模仿真按照侧重点不同分为系统级建模与针对某一部件开展多维CFD仿真。
燃气轮机系统级建模包括机理建模与黑箱建模2种,黑箱建模是基于试验数据库建立,不需要机理,建模更常用[37]。机理建模是基于系统机理得出的数学模型,如林向雨[37]基于PAnysimu仿真平台采用模块化机理建模方法建立了SGT5-4000F型重型燃气轮机的仿真模型。张元哲等[38]在Matlab/Simulink中基于能量守恒、质量守恒及热力学分析采用模块化建模方法建立了单轴重型燃气轮机的数学模型。目前针对燃气轮机系统级的研究多从热动力学角度建立数学模型,一般采用状态方程描述[39]。考虑到系统级GT建模中未表征燃烧室部件的燃烧行为,相关研究针对燃烧室部件建立CFD模型进行分析。如CHIARIELLO等[40]使用商业软件ANSYS对GT系统燃烧室部件在不同负载下的热力学和流体动力学行为进行三维CFD仿真分析。GOBBATO等[41]对H2作为燃料的单轴燃气轮机系统燃烧室部件内的空气-氢反应流动进行CFD分析,并评估了所提出的CFD模型预测燃烧室内部温度场的能力。考虑到CFD模型计算量大等问题,一般不适用于SOFC/GT混合动力系统的仿真。
在GT系统级建模仿真中,根据研究需要可分为稳态模型和动态模型2种,等熵效率模型[35]是常见稳态模型,性能曲线拟合模型[42]是常见动态模型。
压缩机、透平作为燃气轮机系统中关键部件,压缩机、透平部件的稳态模型基于等熵过程热力学公式建立[14],考虑到压缩机和透平实际工作过程并非等熵,因此常引入等熵效率、多变量效率表示循环过程的不可逆性和热损失,虽然多变量效率模型相比等熵效率模型更精确,但由于计算过程较复杂,因此等熵效率模型更适用于系统级建模[43]。压缩机、透平性能曲线拟合模型通过引入压缩机、透平的性能曲线进行系统级GT瞬态建模。性能曲线图是在试验台测试中获得的压缩机、透平详细信息[42],通常包括压缩比、无量纲质量流量、转速之间关系图及等熵效率与无量纲质量流量之间效率图[44]。以压缩机为例,压缩机性能图中,4个基本性能参数压缩比ε、效率η、质量流量或校正后的质量流量
为进入GT燃料/空气的温度,p0为大气压强,pin为GT进气压强)、转速N(或校正后的转速
已知任意2个变量,即可计算其余变量。压缩机进行数学建模时,需要特定格式的映射数据[45],因此需对压缩机性能图进行预处理生成相关映射数据,然后作为数学建模的半经验模型输入。由于压缩机性能图中,喘振区附近的速度轮廓几乎垂直于压力比轴(y轴),缺乏高速与低速数据,因此需借助相关数据处理方法对性能图进行预处理,进而对数据合理外推。KURZKE[45]引入辅助坐标(β线)对性能曲线图进行处理,β线模型是性能曲线图处理中应用最广泛的方法。此外还有Jensen &Kristensen法[46]、穆勒法(Mueller method)、神经网络法等[47]。
换热器、燃烧室等部件常近似为等压开口系统。GT系统级建模中,常忽略燃烧室的燃烧过程,将燃烧室近似作为与外界有热交换的控制体处理,因此常基于质量守恒、能量守恒方程采用简单热力学公式对燃烧室进行零维简单建模分析[13],按相关进气参数(空气、燃料的质量流量、比热容等)计算燃烧室出口温度。换热器部件常基于对数平均温差法[48]、ε-NTU法[49-50]进行设计分析,如薛提微[51]基于改进的传热有效度-传热单元法(ε-φ NTU法)对换热器进行设计分析。但为简化计算,燃气轮机进行系统级建模仿真时一般假设换热器无热损失,近似认为冷热流热交换量相同。
2 SOFC/GT系统布置方式及控制方程
2.1 SOFC/GT混合动力系统布置方式
SOFC/GT混合动力系统根据SOFC与GT布置位置不同,可分为底层循环、顶层循环2种。顶层循环是在SOFC/GT混合系统中将燃料电池置于燃气轮机系统中燃烧室前,SOFC反应后尾气进入GT燃烧室燃烧后进入透平膨胀做功,通过发电机对外发电并为压缩机提供动力[52-53]。底层循环是燃气轮机位于燃料电池上游,透平排气作为燃料电池阴极空气来源。燃料电池与燃气轮机之间通过换热器耦合,即空气经压缩机压缩后,先经由换热器与燃料电池高温尾气发生热交换,回收利用SOFC尾气中废热后再进入透平膨胀做功[52]。
由于燃料电池工作压力不同,顶层循环又称为增压型SOFC/GT系统,如图3所示,底层循环又称为常压型SOFC/GT系统,如图4所示。其中,底层循环效率要低于顶层循环,且底层循环换热器部件成本比较高,SOFC/GT混合动力系统常采用顶层循环。
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/15be805fa97673d0a053555088478751.jpg)
图3 增压式SOFC/GT系统
Fig.3 Pressurized SOFC/GT hybrid system
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/229886b6b49c9239c09a497eb852012f.jpg)
图4 常压式SOFC/GT系统
Fig.4 Non-pressurized SOFC/GT hybrid system
2.2 SOFC/GT混合动力系统数学模型
常见SOFC/GT系统建模按照研究目的可分为稳态模型与动态模型2种,稳态模型主要研究设计工况下的SOFC/GT系统,有时为分析非设计工况下系统行为,会对比分析不同负载下SOFC/GT系统稳态特性。动态模型与稳态模型的区别在于动态模型考虑了相关参数随时间变化,但仍以稳态模型相关数学公式作为基础。本节模型主要计算SOFC输出电压、排气温度、输出功率与GT各部件的排气温度、输出功率,最后计算系统效率,以下分别从建模的基本假设、建模方法和数学模型公式展开介绍。
2.2.1 SOFC建模及基本假设
主要论述SOFC/GT混合动力系统中SOFC稳态模型的电化学反应及热力学相关数学模型公式。该模型以SOFC进气温度、各物质分压、电流密度作为输入,计算SOFC工作电压、工作电流,SOFC对外输出电功的功率可根据工作电流、工作电压计算,通过反复迭代确定SOFC排气温度。
为简化计算,SOFC建模时常做如下假设:① SOFC单电池绝热[4,17];② 阴极出口气体温度等于阳极气体出口温度;③ 阴极气体由O2(21%)、N2(79%)组成,阳极气体由CH4、CO、CO2、H2、H2O组成[54],均为不可压缩理想气体;④ CH4作为燃料的SOFC中,假设阳极只有H2参与电化学反应,CH4只参与蒸汽重整反应,燃料中所有CO均通过水煤气变换反应消耗[17],且上述反应均处于化学平衡状态[54];⑤ 燃料电池内部气体流动均为层流状态,为不可压缩牛顿理想气体混合物,且具有恒定努塞尔数[17]。
2.2.2 SOFC电化学反应数学模型
在CH4作为燃料的SOFC/GT混合动力系统中,SOFC只考虑H2参与的电化学反应(假设④),因此基于能斯特方程可求得SOFC在某一特定工作温度(800~1 000 ℃)下的开路电压E:
(1)
式中,E0为SOFC在某一工作温度下的可逆电动势,常基于经验公式计算[55-56],也可基于某一特定温度下该反应的吉布斯自由能ΔG计算;F为法拉第常数,取96 500 C/mol;R为理想气体常数,取8.314 J/(mol·K);p(i)为i物质的分压,i表示H2、O2、H2O;TSOFC为SOFC工作温度。
E0=1.272 3-2.764 5×10-4TSOFC,
(2)
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/7f8bfeee16c230245c21709598bebbe0.jpg)
(3)
式中,TSOFC计算时可使用阴极空气的进气温度作为初始值进行迭代计算。
但SOFC实际输出工作电压需考虑各种损耗(活化损耗ηact、欧姆损耗ηohmic、浓度损耗ηcon[57])引起的电压下降,因此实际工作电压Vcell表示为
Vcell=E-ηact-ηcon-ηohmic。
(4)
活化损耗ηact可借助巴特勒-沃尔默(B-V)方程确定:
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/05a50739a62f8602fe98292d961fc5ac.jpg)
(5)
式中,j为SOFC工作电流密度;j0为交换电流密度;α为传输系数;ne为电化学反应转移的电子数;TSOFC 为燃料电池工作温度。
在高电流密度下,可用塔菲尔公式近似简化计算活化损耗ηact:
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/0e927d066837f429b02438e2b8408a65.jpg)
(6)
浓度损耗ηcon 计算公式为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/2d2db463ad5314a4eee183be90be1e4d.jpg)
(7)
式中,jL为极限电流密度,表示物质传输的极限情况下(催化层内的反应物浓度降至0时)对应的电流密度。
欧姆损耗ηohmic 为
ηohmic =jReff ,
(8)
式中,Reff为电堆的有效电阻。
SOFC对外所做电功的功率Pel可根据SOFC的输出电流、实际工作电压Vcell计算:
Pel=NcellVcellI。
(9)
其中,Ncell为电堆中单电池个数;I为SOFC的输出电流,可基于单电池有效工作面积Acell计算。建模过程中一般会给出SOFC交换电流密度j0、极限电流密度 jL、有效总电池电阻Reff、传输系数α等重要参数,可确定SOFC电堆的工作电压、工作电流,进而确定Pel。
I=jAcell。
(10)
2.2.3 SOFC稳态热力学模型
基于SOFC电堆能量守恒方程,可计算SOFC排气温度
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/258c456c1bd64650d8261c2502369138.jpg)
(11)
其中,ΔHk为SOFC内部反应焓变,可查表获得参与反应物质的标准生成焓结合盖斯定律计算获得;下标i为电堆中的化学物质,在IRSOFC系统中,i为空气流道中的O2、N2,燃料流道中的H2、CO、CH4等物质;下标k为电堆中发生的化学反应或电化学反应(表1);Rk为第k个反应的化学反应速率;cp,i为SOFC中i物质的定压比热容;上标in、out分别表示进气、排气过程;
为i物质的进气摩尔流量,kmol/s,作为输入量已知;
为i物质的排气摩尔流量,kmol/s,可根据SOFC系统的物料守恒计算:
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/d6715d40cde4fe0ce1de7ca0a9b947d3.jpg)
(12)
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/5a1be1bcc0a8c833ccd9bf4bf91a7159.jpg)
(13)
式中,vi,k为i物质在k反应(Ⅰ、Ⅱ、Ⅲ、Ⅳ)中的化学计量数。
IRSOFC内部k反应的速率Rk为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/b3e35746866367d5f4fddb28743e01e4.jpg)
(14)
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/5f0aefb1cf4c2762768d12e0766c71ed.jpg)
(15)
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/b77372fa4f61304afd340eca21cf3854.jpg)
(16)
其中,Ⅱ、Ⅱ、Ⅲ、Ⅳ表示具体反应方程(表1);kact为MSRR的速率常数;Eact为MSRR的活化能,取82×103 J/mol;kshift为WGSR反应平衡常数;kact、kWGSR分别为MSRR、WGSR的反应速率常数,与SOFC电池内部温度TSOFC有关,计算公式为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/16c0a969795d65a5d617c26df8157e4b.jpg)
(17)
其中,ki为kact、kWGSR;相关系数A、B、C、D、E取值见表3[21,58]。
表3 重整反应及水气变换反应的平衡常数计算公式参数
Table 3 Values of equilibrium constants of reforming and shifting processes
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/9eb4474dd22c485464738a1d7a8b35e6.jpg)
在增压式SOFC/GT混合动力系统中,SOFC排气温度即为GT中燃烧室的进气温度Tin,comb,因此可通过SOFC排气温度将2个系统耦合,进行SOFC/GT系统分析计算。同时,燃料电池功率已知,可进一步计算系统效率。
2.2.4 GT建模及基本假设
SOFC/GT混合动力系统中的GT数学模型主要是对压缩机、透平、换热器等部件进行热力学相关分析。燃气轮机系统每个部件都可看作开口系统。压缩机、透平作为燃气轮机系统中关键部件,排气压力、排气温度、等熵效率和质量流量取决于轴速、进气温度、进气压力和气体热力学特性。一般压缩机进气压力、进气温度作为输入参数已知,通过热力学分析,可计算压缩机、透平的排气压力,通过相关方程联立获得压缩机、燃烧室、换热器、透平各个部件的进气、排气温度。
为简化计算,进行GT建模时做出如下假设:① 压缩机压缩过程和透平膨胀过程近似认为绝热;② 忽略各部件入口出口工作流体动能、势能变化,忽略各部件连接管道的压力损失;③ 工作流体认为是具有恒定比热的理想气体;④ 假设换热器中冷热流热传递为绝热过程,因此冷测与热测的热交换量相同;⑤ 燃烧室燃烧过程近似为恒压热交换控制体处理。
2.2.5 GT控制方程
在SOFC/GT混合系统中,针对GT研究大多是建立热力学模型。其中,压缩机和透平常基于压缩比和等熵效率给出输入输出关系的热力学方程。
以压缩机为例,排气压力pout,c为
pout,c=εcpin,c,
(18)
式中,εc为压缩机压缩比,一般以表格形式列出;pin,c为压缩机进气压力,常取环境压力,作为已知输入条件。
压缩机排气温度Tout,c为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/72da489adca92ce7134fb8f81fbc6577.jpg)
(19)
式中,Tin,c为压缩机进气温度,常取环境温度,作为已知输入条件;ηc为压缩机效率,一般以表格形式列出;κ为气体比热容比,对于空气取κ=1.4。
压缩机对外做功的功率Pc计算为式(20),需要注意的是压缩机消耗系统中透平产能,对外做负功。
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/09b99d9d0b82fbcce7db73e71c8a6099.jpg)
(20)
其中,为压缩机中空气质量流量,透平压缩比、转速已知的情况下可通过透平图获得;cp为空气定压比热,根据工作温度确定。在整个GT系统中,为简化计算,一般近似认为cp为定值,不考虑随温度的变化。
在GT控制方程中,压缩机进气温度、进气压力为环境温度、环境大气压。相关研究通常会注明GT型号并给出GT在设计工况下的重要参数,如压缩比、压缩机等熵效率、透平等熵效率、空气质量流量、燃料质量流量、换热器效率等,这些参数常作为系统输入参数给出[59]。LI等[16]在工作中以表格形式给出了CAPSTONE公司的C30型单轴微型燃气轮机设计工况相关参数。ROSSI等[59]给出了Turbec T100燃气轮机设计工况下相关参数。常引入压缩机性能曲线图、透平性能曲线图对不同负载下非设计工况进行分析。
与压缩机类似,透平排气压力pout,t为
pout,t=εtpin,t,
(21)
式中,εt为表示压缩机压缩比;pin,t为压缩机进气压力。
透平排气温度Tout,t为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/c1f2c987343f06acacb3d5ce9f1b1c32.jpg)
(22)
式中,Tin,t为透平的进气温度,近似等于燃烧室排气温度;ηt为透平效率,相关研究一般会以表格形式列出。
透平对外做功的功率Pt为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/1cdc8569265c1ccda619bd4cda879a0f.jpg)
(23)
换热器可基于假设④对冷流、热流之间的热交换进行计算:
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/5ad8e8662bbc6ef4c1536c70c2b87777.jpg)
(24)
Ph=m2cp2(Tout,h-Tin,h),
(25)
Tin,l=Tout,c,
(26)
Tin,h=Tout,t,
(27)
式中,Pl为冷流的热交换功率;Ph为热流的热交换功率;Tin,l为换热器冷流的进气温度,近似等于压缩机排气温度;Tout,l为冷流排气温度;Tin,h为热流的进气温度,近似等于透平的排气温度;Tout,h为热流的排气温度;cp1、cp2分别为冷流和热流工作温度下对应的定压比热容;分别为冷流、热流的质量流量。
甲烷作为燃料的燃烧室部件排气温度Tout,comb为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/9cae431d0fc990a3e3d5b621ec3303e9.jpg)
(28)
式中,Tin,comb近似等于SOFC的排气温度;at为完全燃烧1 kg燃料所需理论空气质量,kg;n为过量空气系数;LHV(CH4)为甲烷的低热值,MJ/kg。
Tin,comb=,
(29)
Tin,t=Tout,comb,
(30)
at=11.56xC 34.48xH,
(31)
LHV(CH4)=37.966 MJ/kg。
(32)
式中,xi为甲烷中i元素的质量分数(i为C、H)。
联立式(19)(22)(26)(27)(28)(29)(30)可计算各部件的进气、排气温度,并与试验数据进行对比,进而验证所建模型的准确性,有助于进一步分析SOFC/GT系统性能。
2.2.6 SOFC/GT混合动力系统效率计算
常见SOFC/GT系统性能分析主要用于优化SOFC/GT系统效率,因此需对系统各部件效率进行分析计算,以确定造成系统效率低下的部件及原因。因此,可基于SOFC燃料利用率、实际工作电压、电化学反应的吉布斯自由能变、焓变等参数计算SOFC效率,结合GT的效率计算SOFC/GT混合动力系统的总体效率。
SOFC效率εSOFC为
εSOFC=εthermoεvoltageεfuel。
(33)
式中,εthermo为燃料电池的可逆热效率;εvoltage为燃料电池的电压效率;εfuel为燃料电池的燃料利用率。
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/e4f0947f83eee3c2f9a5e8547c667168.jpg)
(34)
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/899aaaeaead67bc88fcff3e49e971662.jpg)
(35)
εfuel=Uf,
(36)
式中,ΔH为SOFC工作温度下内部电化学反应的焓变;E0为SOFC热力学可逆电压;Uf为燃料利用率,一般设为固定值。
其中,对于氢作为反应物的电化学反应,在标准状况下,反应方程式对应的吉布斯自由能变ΔGθ与焓变ΔHθ如式(37)所示,显然SOFC实际工作温度远高于标准态(298.15 K),因此对于SOFC在某一特定工作温度(800~1 000 ℃)、工作压力下效率的计算可考虑查找水蒸气的焓熵图确定焓、熵值[60],进而计算该条件下的吉布斯自由能变,或通过引入NASA多项式计算某一特定工作温度下的反应焓变[61]。
ΔHθ=-241.818 kJ/mol。
(37)
对于燃气轮机而言,效率表示为实际做功与通入燃料化学能的比值,其中燃料化学能用LHV表示,因此,燃气轮机效率εGT为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/c7bdf362d33c83af3936d4ce6c4c5a98.jpg)
(38)
其中,为进入SOFC燃料的质量流量,等于燃料中各组分的进气摩尔流量
之和,kmol/s,具体为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/4a43204483a8328bee7738b749d53f84.jpg)
(39)
CH4作为燃料的SOFC/GT混合动力系统效率εSOFC/GT为
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/f36edaae0c7c84802c6ca2e752e1d078.jpg)
(40)
式中,Pel为GT对外输出电功的功率。
基于上述公式可以初步分析SOFC、GT及混合系统效率。
3 SOFC/GT系统示范性工程及相关建模仿真
3.1 SOFC/GT混合动力系统示范性工程
Westinghouse Electric公司较早开始SOFC研究,并在此基础上研发了100 kW级SOFC/CHP系统[62]。同时,Siemens Westinghouse Power公司与加州大学合作在加州大学国家燃料电池中心安装了第1台220 kW TSOFC/GT混合发电系统,该系统为典型的顶层循环SOFC/GT混合动力系统,效率为53%,其中SOFC功率为180 kW,燃气轮机功率为40 kW[44],该混合动力系统也是目前应用最广泛的示范模型,广泛用于数学仿真模型验证[44,63-65]。作为瑞典高温燃料电池国家计划的一部分,瑞典隆德大学热能与动力工程系开展了SOFC/GT联合循环的理论研究[66]。Ceres Power公司开发了一种微型热电联产SOFC发电装置,其中采用金属支撑的SOFC燃料电池电堆,工作温度为450~600 ℃[67]。日本三菱重工公司于2012年开发了1种三联合循环发电系统[68],该系统经燃料电池、燃气轮机和蒸汽轮机3个子系统发电,由此产生的燃料电池联合循环系统在天然气作为燃料时发电效率预计超过60%;于2015年向九州大学交付了SOFC-MGT混合循环系统[69],该系统与九州大学校园电力系统连接,发电功率为250 kW,效率为55%。总体而言,SOFC/GT混合动力系统发电效率高,可作为小型分布式发电系统普及推广。
3.2 SOFC/GT系统相关建模仿真
考虑到SOFC/GT混合动力系统示范性工程的建造成本,开展SOFC/GT数学建模仿真是一种高效、低成本的研究方法[44,65]。鉴于BUONOMANO等[15]详细综述2015年及以前SOFC/GT混合动力系统的布局、耦合方式,且BAO等[13]论述了2018年及以前SOFC/GT混合动力系统布局[70-71]、安全运行工况及约束参数[72-75]、部分负载非设计工况[22,54,76]、控制策略等[77-79]研究工作,因此本文重点关注2019年及以后的SOFC/GT混合动力系统数学建模仿真工作,同时围绕SOFC/GT混合动力系统不同建模方法展开综述。
大部分SOFC/GT混合动力系统建模仿真工作围绕稳态热力学建模仿真展开。如WILSON等[80]基于ProMax仿真软件建立并验证了用于电动航空的SOFC/GT混合动力系统稳态热力学模型,分析SOFC设计和运行对混合动力系统性能的影响。DING等[6]基于Matlab/simulink建立了SOFC/GT热力学模型,对以不同类型生物质气体为燃料的ITSOFC/GT混合动力系统进行热力学分析。NI等[9]对SOFC/GT各子系统(包括SOFC、后燃烧室、压缩机、透平和换热器)分别进行热力学分析,设计了一种顶部循环为布雷顿超临界CO2循环,底部循环为ORC的SOFC/GT系统,进一步提高系统效率。NI等[9]基于Matlab软件对SOFC/GT进行热力学建模,分析了燃料利用率、最大压力、最高温度等参数对系统性能的影响。同时,由于该模型对SOFC进行了简化,还与其他文献仿真结果进行对比以确保模型的可靠性。BADUR等[81]基于非开源代码COM-GAS建立IR-TSOFC零维鲁棒模型对SOFC电流、电压、电压损失、燃烧成分、功率输出、效率等基本参数进行预测,并与现有试验数据进行比较。结果表明仿真差异不超过5%;SOFC/GT顶层循环混合动力系统可显著提高效率(提升至55%~58%)。
为进一步提升SOFC/GT混合动力系统效率,需关注造成系统效率低下的部件,同时确定该部件造成系统效率低下的程度、原因。相关研究常在SOFC/GT稳态热力学模型中引入分析,即确定各部件能的品位,进而评估效率提升的潜力与可行性[82]。GHOLAMIAN等[23]基于热力学第一和第二定律建立了SOFC/GT系统热力学模型,并使用工程方程求解器EES进行求解,研究SOFC/GT系统耦合2种循环(ORC、卡丽娜循环)对系统效率的影响。结果表明,SOFC/GT-ORC复合体系的
效率更高(达62.35%)。EISAVI等[24]在EES中对建立的SOFC/GT热力学模型进行求解,通过参数化研究模拟循环性能,从能源、
、环境和经济等方面对3种不同配置的SOFC/GT内重整混合动力系统进行对比分析。结果表明,相较单电堆的SOFC/GT系统,使用双电堆(每个电堆的单电池数量减半)的SOFC/GT系统经济性更好;同时采用空气串联进给的双电堆SOFC/GT系统的系统效率相较采用平行进给系统更高,能量效率和
效率分别达63.22%和60.81%,比传统单电堆SOFC/GT系统分别提高20.30%和8.89%。CALISE等[83]利用Matlab软件对IRSOFC/GT混合动力系统进行热力学建模仿真,在不同操作压力、燃料利用率、燃料空气比、蒸汽燃料比和电流密度下对系统及各部件进行能量和
分析,确定系统效率低下的来源。结果表明,由于化学和电化学反应(蒸汽重整和电化学氧化)的发生导致IRSOFC电堆成为系统
破坏的最重要来源,电堆也是不可逆性最高组分,因此SOFC/GT混合动力系统设计需重点关注SOFC电堆部分。BEIGZADEH等[84]利用Matlab软件基于能量平衡方程对SOFC/GT混合动力系统各部件进行热力学建模,考虑蒸汽重整情况下对系统进行能量和
分析。从热力学角度比较天然气和其他生物燃料(包括污水沼气、农业和工业废物沼气、合成气、生物燃料和气化生物质)对单个设备和整个系统性能的影响,发现以天然气为燃料的系统电力效率最高,而生物燃料的总电力效率最低。
由于SOFC集总模型(0D)中电池温度、电流密度和燃料利用率等参数均采用集总平均值[85],而SOFC工作温度较高,沿电池长度方向的温度梯度对电池正常工作十分重要,此时热力学平衡和集总模型的预测结果会存在较大偏差,且无法考虑顺流和逆流等流动形式的影响[65],因此无法准确描述燃料电池的局部特性。为进一步提高建模精度以获得更真实的温度分布等性能参数,对SOFC/GT混合动力系统进行性能分析时,通常引入SOFC一维模型[86]。因此CHEN等[86]基于Matlab软件建立SOFC、后燃烧室、重整器集总模型,同时基于Ebsilon软件建立GT模型,并与NETL开发的1D SOFC模型进行耦合,研究燃料利用率和燃料重整过程对100 kW SOFC/GT混合动力系统性能的影响。结果表明,燃料电池Uf越高,系统效率越高;同时SOFC/GT系统采用阳极重整效率高于外部重整。LIU等[87]通过建立SOFC一维模型,预测了O/C比、S/C比(原子比)和重整温度等重整条件对正癸烷重整SOFC/GT混合动力航空系统性能的影响。结果表明,高O/C比、S/C比和转化温度对碳的形成有抑制作用;改变转化温度或O/C比、S/C比是减小SOFC温度梯度的有效途径。HUANG等[88]基于平衡态重整器模型、一维SOFC模型及搭建于Ebsilon软件中的燃气轮机系统模型,分析了不同燃料预重整份额和SOFC燃料利用率在设计工况下对SOFC/GT混合发电系统的性能影响,仿真结果表明,通过合理配置SOFC燃料利用率和燃料重整过程有望在不依赖燃料内重整的情况下提高混合发电系统的发电效率。同时,HUANG等[89]还建立了包括基于Matlab软件的SOFC一维动态模型、蒸汽甲烷重整器集总模型及基于Ebsilon软件的GT循环模型的100 kW级 NG SOFC/GT混合动力系统模型,分析了SOFC工作温度和燃料利用率对SOFC/GT混合动力系统净发电效率的影响。结果表明,较高的工作温度可提高整个混合动力系统效率,但过高的燃料利用率可能会导致系统效率降低。CHEN等[90]基于Ebsilon软件建立了燃料重整器及顶层一维SOFC/GT混合循环系统模型,对回热式燃气轮机循环相关性能进行模拟。发现与阳极内部重整相比,燃料外部重整器可灵活实现低燃料利用率、高燃料预重整率,同时保持较高系统效率。
大多数研究工作集中于SOFC/GT混合动力装置的稳态建模仿真,以确定实现混合动力系统的最佳性能系统配置[9],与之相比,关于SOFC/GT动态建模文献较少[64]。SOFC/GT的动态建模仿真是系统进行负荷跟随控制的基础[65],混合动力系统动态行为的建模、优化与控制研究对于系统集成和商业化普及具有重要意义[91]。相关研究常基于Matlab[64]软件建立SOFC/GT动态模型,同时常以热那亚大学的热化学电力集团(TPG)[92]、德国航空航天中心[93]和美国国家能源技术实验室(NETL)[59,94]开发的SOFC/GT半实物动态仿真平台获得的试验数据作为参考,对所建立动态模型的准确性进行验证[44,59,95],进而分析系统动态特性。ROBERTS等[44]基于Matlab/Simulink建立了顶层循环TSOFC/GT混合动力系统模型对Siemens Westinghouse Power公司开发的220 kW TSOFC/GT系统进行启动过程的动态仿真,并与NETL开发的220 kW SOFC/GT动态模型及运行测试过程中收集的试验数据进行比较。结果表明,该动态模型可以很好捕捉启动期间系统的动态性能,但未捕捉到旁通阀相关系统动力学。FERRARI等[95]展示了NETL基于Matlab/Simulink开发的瞬态模型,并与TPG基于TRANSEO工具开发的瞬态模型进行对比,结果表明,NETL开发的瞬态模型整体上可以对系统非设计工况下的瞬态响应进行实时预测,与之相比,TPG开发的瞬态模型精确度更高,但计算时间较长。在完成SOFC/GT动态模型的搭建后常对SOFC/GT系统瞬态过程(如启停工况[44,59])的动态特性及非设计工况进行分析、优化。BARELLI等[64]基于Matlab/Simulink建立SOFC/GT动态模型进行建模仿真进而对整个混合动力系统进行动力学分析。考虑到SOFC与GT之间的相互作用以优化系统组件,同时对系统效率、时间响应进行评估,特别关注燃气轮机和热交换器对燃料电池的惯性效应。开发的动态模型可用于分析SOFC/GT混合动力系统的整体性能,在燃气轮机部分负载下优化设备组件及控制逻辑,并填补SOFC/GT混合系统组件之间动态相互作用分析。
除基于Matlab软件开发SOFC/GT动态模型外,ZHANG等[96]基于质量守恒方程、能量守恒方程和动量守恒方程在仿真工具Aspen Custom Modeler中建立了SOFC/GT动态模型。其中,SOFC模型基于指数衰减函数和指数关联函数建立。利用所建动态模型进行动态仿真时,选取电流密度作为扰动变量,计算结果表明,该动态模型能较好拟合系统的动态特性。包成等[65]基于gPROMS平台建立PSOFC/GT动态模型,并将仿真结果与SWPC 220 kW TSOFC/GT示范系统额定工况试验数据、美国燃料电池工程中心APSAT工具预测结果进行对比。结果表明,该文献中模型相比于APSAT具有更好预测精度,可用于系统动态特性分析,为后续控制器设计奠定基础。同时,考虑到Matlab等软件对技术人员要求较高,难以商业化推广,ROSSI等[59]基于经验参数在Microsoft Excel/Visual Basic中建立了SOFC/GT混合动力系统动态模型的简化框架,分析系统启动工况,同时将该模型的响应与NETL瞬态模型的数据进行比较。结果表明,该模型能实现快速响应的同时具有良好的预测精确度,因此更适于商业化推广。
4 结语与展望
1)燃料电池分布式发电技术是适应未来能源低碳化、清洁化、高效化发展趋势的重要应用方向。其中SOFC作为高温燃料电池,可发生燃料内部重整反应,燃料适应性好;同时,SOFC可与GT等装置构成混合发电系统,GT回收利用SOFC高温尾气中的能量进一步提高发电效率。
2)考虑到SOFC/GT示范性工程有限且建造成本较高,一般采用数学建模仿真方法开展相关工作。数学建模仿真对SOFC/GT混合发电系统发展有重要影响。与单独SOFC或GT模型不同,SOFC/GT混合动力系统常基于宏观尺度进行热力学分析。SOFC/GT集总模型通过计算温度、压力等主要参数对系统稳态性能进行分析优化。由于集总模型无法准确描述燃料电池的局部特性,常引入一维SOFC的CFD模型对燃料重整过程进行分析,但现阶段SOFC/GT混合系统建模中引入一维SOFC模型研究较少,后续工作中可考虑在SOFC/GT系统中引入一维甚至更高维度的SOFC模型以提高建模精度。
3)SOFC/GT动态建模仿真是系统进行负荷跟随控制的基础,对于SOFC/GT混合发电系统商业化普及具有重要意义。现阶段进行瞬态及非设计工况下的建模仿真工作较少,在后续工作中可考虑进行瞬态/非设计工况的相关建模工作。
4)现阶段SOFC/GT混合动力系统建模主要基于Matlab/Simulink、Aspen Plus等商业化软件开展模块化建模,因此功能有限、不易拓展,建立的模型针对性有限,很难对SOFC/GT中某一特定部件进行参数调整,难以实现某些特定功能的仿真。同时,基于商业化软件建立的SOFC/GT混合系统模型较难耦合SOFC、GT的高维CFD仿真模型,计算精度很难提高。后续工作中可尝试对SOFC/GT混合动力系统开发专用仿真软件,提高模型的拓展性和针对性;同时耦合高维度SOFC、GT CFD仿真模型进一步提高计算精度。
[1] 国家发展改革委、国家能源局等.“十四五”可再生能源发展规划[EB/OL].(2022-06-01)[2023-02-10]https://www.ndrc.gov.cn/xxgk/zcfb/ghwb/202206/P020220601501054858882.pdf.
[2] 皇甫成, 郭金智, 贾非, 等. 一种考虑新能源出力相关性及不确定性的分布式发电规划策略[J].电网与清洁能源, 2022,38(3):105-110,118.
HUANG Pucheng, GUO Jinzhi, JIA Fei, et al. A distributed generation planning strategy considering the correlation and uncertainty of renewable energy output[J]. Advances of Power System &Hydroelectric Engineering, 2022,38(3):105-110,118.
[3] 钟财富. 国内外分布式燃料电池发电应用现状及前景分析[J]. 中国能源,2021,43(2):34-37,73.
ZHONG Caifu. Application status and prospect analysis of distributed fuel cell power generation at home and abroad[J]. Energy of China,2021,43(2):34-37,73.
[4] 李卓言, 段丽平, 李少华, 等. 基于AspenPlus的兆瓦级燃料电池分布式发电系统建模及仿真分析 [J]. 南方能源建设, 2022, 9(4): 78-86.
LI Zhuoyan, DUAN Liping, LI Shaohua, et al. Modeling and simulation analysis of MW-level fuel cell distributed generation system based on AspenPlus [J].Southern Energy Construction, 2022, 9 (4):78-86.
[5] STEELE B C H, HEINZEL A. Materials for fuel-cell technologies[J]. Nature, 2001, 414: 345-352.
[6] DING X, LV X, WENG Y. Fuel-adaptability analysis of intermediate-temperature-SOFC/gas turbine hybrid system with biomass gas[J]. Journal of Energy Resources Technology, 2021, 143(2): 022104-1-022104-10.
[7] DUONG P A, RYU B, JUNG J, et al. Thermal evaluation of a novel integrated system based on solid oxide fuel cells and combined heat and power production using ammonia as fuel [J]. Applied Sciences, 2022, 12(12): 6287.
[8] TOONSSEN R, SOLLAI S, ARAVIND P V, et al. Alternative system designs of biomass gasification SOFC/GT hybrid systems[J]. International Journal of Hydrogen Energy, 2011, 36(16): 10414-10425.
[9] NI T M, SI J W, GONG X H, et al. Thermodynamic and economic analysis of a novel cascade waste heat recovery system for solid oxide fuel cell[J]. Energy Conversion and Management, 2022, 259:115562-1-115562-14.
[10] LI Z, ZHANG X, HE X, et al. Comparative analysis of thermal economy of two SOFC-GT-ST triple hybrid power systems with carbon capture and LNG cold energy utilization [J]. Energy conversion and management, 2022, 256: 115385-1-115385-25.
[11] CHOUDHURY A, CHANDRA H, ARORA A. Application of so-lid oxide fuel cell technology for power generation:A review [J]. Renewable and Sustainable Energy Reviews, 2013, 20: 430-442.
[12] SARAVANAMUTTOO H I, ROGERS G F C, COHEN H. Gas turbine theory[M].London: Pearson Education, 2001.
[13] BAO C, WANG Y, FENG D, et al. Macroscopic modeling of solid oxide fuel cell (SOFC) and model-based control of SOFC and gas turbine hybrid system [J]. Progress in Energy and Combustion Science, 2018, 66: 83-140.
[14] ZHOU X, PRAMUANJAROENKIJ A, KAKAÇ S. A review on miniaturization of solid oxide fuel cell power sources-II: From system to material [J]. Mini-Micro Fuel Cells, 2008: 319-347. https://doi.org/10.1007/978-1-4020-8295-5_22.
[15] BUONOMANO A, CALISE F, D′ACCADIA M D, et al. Hybrid solid oxide fuel cells-gas turbine systems for combined heat and power: A review [J]. Applied Energy, 2015, 156: 32-85.
[16] LI Y, WENG Y. Performance study of a solid oxide fuel cell and gas turbine hybrid system designed for methane operating with non-designed fuels [J]. Journal of Power Sources, 2011, 196(8): 3824-3835.
[17] WANG L, ZHANG H, WENG S. Modeling and simulation of solid oxide fuel cell based on the volume-resistance characteristic modeling technique [J]. Journal of Power Sources, 2008, 177(2): 579-589.
[18] KANEKO T, BROUWER J, SAMUELSEN G S. Power and temperature control of fluctuating biomass gas fueled solid oxide fuel cell and micro gas turbine hybrid system [J]. Journal of Power Sources, 2006, 160(1): 316-325.
[19] SERVICES E G G,RALPH M. Fuel Cell Handbook[M].[S.l.]:DIANE Publishing, 2000.
[20] 贾旭平.美国Bloom energy公司Bloomenergy推出固体氧化物燃料电池微型电站[EB/OL].(2021-01-29)[2023-02-10].https://xueshu.baidu.com/usercenter/paper/show?paperid=559bea83d19c7cb37ff2b874a46bd90b.
[21] CHAN S, HO H, TIAN Y. Modelling of simple hybrid solid oxide fuel cell and gas turbine power plant [J]. Journal of Power Sources, 2002, 109(1): 111-120.
[22] CHAN S H, HO H K, TIAN Y. Modelling for part-load oper-ation of solid oxide fuel cell-gas turbine hybrid power plant [J]. Journal of Power Sources, 2003, 114(2): 213-227.
[23] GHOLAMIAN E, ZARE V. A comparative thermodynamic investigation with environmental analysis of SOFC waste heat to power conversion employing Kalina and Organic Rankine Cycles [J]. Energy Conversion and Management, 2016, 117: 150-161.
[24] EISAVI B, CHITSAZ A, HOSSEINPOUR J, et al. Thermo-environmental and economic comparison of three different arrangements of solid oxide fuel cell-gas turbine (SOFC-GT) hybrid systems [J]. Energy Conversion and Management, 2018, 168: 343-356.
[25] YAKABE H, OGIWARA T, HISHINUMA M, et al. 3-D model calculation for planar SOFC [J]. Journal of Power Sources, 2001, 102(1): 144-154.
[26] BESSETTE N F, WEPFER W J, WINNICK J. A mathematical model of a solid oxide fuel cell [J]. Journal of The Electrochemical Society, 1995, 142(11): 3792-3800.
[27] RAMREZ-MINGUELA J J, RANGEL-HERN
NDEZ V H, ALFARO-AYALA J A, et al. Energy and entropy study of a SOFC using biogas from different sources considering internal reforming of methane [J]. International Journal of Heat and Mass Transfer, 2018, 120: 1044-1054.
[28] JANARDHANAN V M, DEUTSCHMANN O. CFD analysis of a solid oxide fuel cell with internal reforming: Coupled interactions of transport, heterogeneous catalysis and electrochemical processes [J]. Journal of Power Sources, 2006, 162(2): 1192-1202.
[29] CHAN S H, HO H K, TIAN Y. Multi-level modeling of SOFC-gas turbine hybrid system [J]. International Journal of Hydrogen Energy, 2003, 28(8): 889-900.
[30] RECKNAGLE K P, WILLIFORD R E, CHICK L A, et al. Three-dimensional thermo-fluid electrochemical modeling of planar SOFC stacks[J]. Journal of Power Sources, 2003, 113(1): 109-114.
[31] SUWANWARANGKUL R, CROISET E, PRITZKER M D, et al. Mechanistic modelling of a cathode-supported tubular solid oxide fuel cell [J]. Journal of Power Sources, 2006, 154(1): 74-85.
[32] ACKMANN T, DE HAART L G J, LEHNERT W, et al. Modeling of mass and heat transport in planar substrate type SOFCs [J]. Journal of the Electrochemical Society, 2003, 150(6): A783-A789.
[33] ZHU H, KEE R J. Modeling electrochemical impedance spectra in SOFC button cells with internal methane reforming[J]. Journal of The Electrochemical Society, 2006, 153(9): A1765.
[34] SARAVANAMUTTOO H I H, ROGERS G F C, COHEN H. Gas turbine theory[M]. London: Pearson Education, 2001.
[35] CHINDA P, BRAULT P. The hybrid solid oxide fuel cell (SO-FC) and gas turbine (GT) systems steady state modeling [J]. International Journal of Hydrogen Energy, 2012, 37(11): 9237-9248.
[36] SONG T, KIM T, KIM J, et al. Performance prediction of axial flow compressors using stage characteristics and simultaneous calculation of interstage parameters[J]. Journal of Power and Energy, 2001,215(1): 89-98.
[37] 林向雨. SGT5-4000F燃气轮机本体建模与仿真[D].北京:华北电力大学,2015.
[38] 张元哲, 刘培, 李政. 燃气轮机机理建模与性能监测 [J]. 工程热物理学报,2020,41(3): 548-553.
ZHANG Yuanzhe, LIU Pei, LI Zheng. Mechanism modeling and performance monitoring of gas turbine. [J]. Journal of Engineering Thermophysics, 2020,41(3): 548-553.
[39] ROSSI I, TRAVERSO A, HOHLOCH M, et al. Physics-based dynamic models of three SOFC/GT emulator test rigs[J]. Journal of Engineering for Gas Turbines and Power, 2018, 140(5): 051702-1-051702-10.
[40] CHIARIELLO F, REALE F, CALABRIA R, et al. Off design behavior of a 100 kW turbec T100P micro gas turbine[J]. Applied Mechanics and Materials, 2013, 390: 275-280.
[41] GOBBATO P, MASI M, TOFFOLO A, et al. Numerical simulation of a hydrogen fuelled gas turbine combustor [J]. International Journal of Hydrogen Energy, 2011, 36(13): 7993-8002.
[42] TSOUTSANIS E, MESKIN N, BENAMMAR M, et al. An efficient component map generation method for prediction of gas turbine performance[C]//Proceedings of the Turbo expo: Power for Land, Sea, and Air.Düsseldorf:American Society of Mechanical Engineers,2014.
[43] CUNNEL C, PANGALIS M, MARTINEZ-BOTAS R. Integration of solid oxide fuel cells into gas turbine power generation cycles. Part 2: Hybrid model for various integration schemes[J]. Journal of Power and Energy, 2002, 216(2): 145-154.
[44] ROBERTS R A, BROUWER J. Dynamic simulation of a pressurized 220 kW solid oxide fuel-cell-gas-turbine hybrid system: Modeled performance compared to measured results [J]. Journal of Fuel Cell Science and Technology, 2006, 3(1):18-25.
[45] KURZKE J. How to get component maps for aircraft gas turbine performance calculations[C]//Proceedings of the Turbo Expo: Power for Land, Sea, and Air. Düsseldorf:American Society of Mechanical Engineers,1996.
[46] JENSEN J P, KRISTENSEN A, SORENSON S C, et al. Mean value modeling of a small turbocharged diesel engine [R]. Detroit, Michigan: SAE Technical Paper, 1991.
[47] MORAAL P, KOLMANOVSKY I. Turbocharger modeling for automotive control applications [R].Detroit, Michigan: SAE Technical Paper, 1999.
[48] CARTAXO S J, FERNANDES F A. Counterflow logarithmic me-an temperature difference is actually the upper bound: A demonstration [J]. Applied Thermal Engineering, 2011, 31(6/7): 1172-1175.
[49] NAVARRO H, CABEZAS-GMEZ L. Effectiveness-NTU computation with a mathematical model for cross-flow heat exchangers [J]. Brazilian Journal of Chemical Engineering, 2007, 24: 509-521.
[50] NAVARRO H A, CABEZAS-GMEZ L. A new approach for thermal performance calculation of cross-flow heat exchangers [J]. International Journal of Heat and Mass Transfer, 2005, 48(18): 3880-3888.
[51] 薛提微. 换热器的热力设计方法与性能评价[D].济南:山东大学, 2017.
[52] SAMUELSEN S, BROUWER J. Fuel cell/gas turbine hybrid[J]. Encyclopedia of Electrochemical Power Sources, 2009:DOI:10.1016/B978-044452745-5.00911-4.
[53] 翁一武, 苏明, 翁史烈. 先进微型燃气轮机的特点与应用前景 [J]. 热能动力工程, 2003, 6(2): 111-116,212.
WENG Yiwu, SU Ming, WENG Shilie. Specific features of advanced micro gas turbines and their application prospects [J]. Journal of Engineering for Thermal Energy and Power, 2003, 6(2): 111-116,212.
[54] COSTAMAGNA P, MAGISTRI L, MASSARDO A F. Design and part-load performance of a hybrid system based on a solid oxide fuel cell reactor and a micro gas turbine [J]. Journal of Power Sources, 2001, 96(2): 352-368.
[55] APPLEBY A J, FOULKES F R. Fuel cell handbook[M]. New York: Van Nostrand Reinhold, 1988.
[56] CAMPANARI S, IORA P. Definition and sensitivity analysis of a finite volume SOFC model for a tubular cell geometry [J]. Journal of Power Sources, 2004, 132(1/2): 113-126.
[57] O′HAYRE Ryan,车硕源,COLELLA Whitney,等. 燃料电池基础[M].北京:电子工业出版社, 2007.
[58] MASSARDO A F, LUBELLI F. Internal reforming solid oxide fu-el cell-gas turbine combined cycles (IRSOFC-GT): Part A:Cell model and cycle thermodynamic analysis[J]. Journal of Engineering for Gas Turbines and Power, 2000, 122(1): 27-35.
[59] ROSSI I, TRAVERSO A, TUCKER D. SOFC/Gas Turbine Hybrid System: A simplified framework for dynamic simulation [J]. Applied Energy, 2019, 238: 1543-1550.
[60] 金康, 付照中, 吴滨, 等. 基于Soave-Redlich-Kwong气体物态方程的CH4焓、熵图数值计算及关联应用程序开发 [J]. 云南大学学报(自然科学版), 2020, 42(3): 485-491.
[61] SHAFAGH I, HUGHES K J, POURKASHANIAN M. Modified enthalpies of formation for hydrocarbons from DFT and ab initio thermal energies [J]. Computational and Theoretical Chemistry, 2011, 964(1/3): 100-107.
[62] HUANG K, SINGHAL S C. Cathode-supported tubular solid ox-ide fuel cell technology: A critical review [J]. Journal of Power Sources, 2013, 237: 84-97.
[63] YI Y. Simulation of a 220 kW hybrid SOFC gas turbine system and data comparison [J]. ECS Proceedings Volumes, 2003, 7(1): 1442-1454.
[64] BARELLI L, BIDINI G, OTTAVIANO A. Part load operation of a SOFC/GT hybrid system: Dynamic analysis [J]. Applied Energy, 2013, 110: 173-189.
[65] 包成, 史翊翔, 李晨, 等. SOFC-GT混合系统的参数分布与动态特性数值模拟 [J]. 工程热物理学报, 2010, 31(3): 469-472.
BAO Cheng, SHI Yixiang, LI Chen, et al. Numeric simulation of parameter distribution and transient behaviors for SOFC-GT hybrid system [J]. Journal of Engineering Thermophysics, 2010, 31(3): 469-472.
[66] PALSSON J, SELIMOVIC A, SJUNNESSON L. Combined solid oxide fuel cell and gas turbine systems for efficient power and heat generation [J]. Journal of Power Sources, 2000, 86(1/2): 442-448.
[67] LEAH R T, BONE A, SELCUK A, et al. Development of highly robust, volume-manufacturable metal-supported SOFCs for operation below 600 ℃[J]. ECS Transactions, 2011, 35(1): 351.
[68] None. Mitsubishi to develop SOFC-turbine triple combined cycle system [J]. Fuel Cells Bulletin, 2012, 2012(7): 5-6.
[69] ANDO Y, OOZAWA H, MIHARA M, et al. Demonstration of SOFC-micro gas turbine (MGT) hybrid systems for commercialization[J]. Mitsubishi Heavy Industries Technical Review, 2015, 52(4): 47-52.
[70] MASSARDO A F, LUBELLI F. Internal reforming solid oxide fu-el cell-gas turbine combined cycles (IRSOFC-GT): Part A. Cell model and cycle thermodynamic analysis [J]. Journal of Engineering for Gas Turbines and Power, 2000, 122(1): 27-35.
[71] RAO A D, SAMUELSEN G S. A thermodynamic analysis of tubular solid oxide fuel cell based hybrid systems [J]. Journal of Engineering for Gas Turbines and Power, 2003, 125(1): 59-66.
[72] STILLER C, THORUD B, BOLLAND O. Safe dynamic operation of a simple SOFC/GT hybrid system [J]. Journal of Engineering for Gas Turbines and Power, 2006, 128(3): 551-559.
[73] LV X, LIU X, GU C, et al. Determination of safe operation zone for an intermediate-temperature solid oxide fuel cell and gas turbine hybrid system [J]. Energy, 2016, 99: 91-102.
[74] STILLER C. Design, operation and control modelling of SOFC/GT hybrid systems [M]. Tromso:Norges Arktiske Universitet,2006.
[75] COSTAMAGNA P, HONEGGER K. Modeling of solid oxide heat exchanger integrated stacks and simulation at high fuel utilization [J]. Journal of The Electrochemical Society, 1998, 145(11): 3995-4007.
[76] BAVARSAD P G. Energy and exergy analysis of internal reforming solid oxide fuel cell-gas turbine hybrid system [J]. International Journal of Hydrogen Energy, 2007, 32(17): 4591-4599.
[77] STILLER C, THORUD B, BOLLAND O, et al. Control strategy for a solid oxide fuel cell and gas turbine hybrid system [J]. Journal of Power Sources, 2006, 158(1): 303-315.
[78] ROBERTS R, BROUWER J, JABBARI F, et al. Control design of an atmospheric solid oxide fuel cell/gas turbine hybrid system: Variable versus fixed speed gas turbine operation [J]. Journal of Power Sources, 2006, 161(1): 484-491.
[79] ROBERTS R A, BROUWER J, LIESE E, et al. Dynamic simulation of carbonate fuel cell-gas turbine hybrid systems [J]. Journal of Engineering for Gas Turbines and Power, 2006, 128(2): 294-301.
[80] WILSON J A, WANG Y, CARROLL J, et al. Hybrid solid oxide fuel cell/gas turbine model development for electric aviation [J]. Energies, 2022, 15(8): 2885.
[81] BADUR J, M, KOWALCZYK T, et al. Zero-dimensional robust model of an SOFC with internal reforming for hybrid energy cycles [J]. Energy, 2018, 158: 128-138.
[82] KOTAS T J. The exergy method of thermal plant analysis [M]. New York:Paragon Publishing, 2012.
[83] CALISE F, D′ACCADIA M D, PALOMBO A, et al. Simulation and exergy analysis of a hybrid solid oxide fuel cell (SOFC)-gas turbine system [J]. Energy, 2006, 31(15): 3278-3299.
[84] BEIGZADEH M, POURFAYAZ F, GHAZVINI M, et al. Energy and exergy analyses of solid oxide fuel cell-gas turbine hybrid systems fed by different renewable biofuels: A comparative study [J]. Journal of Cleaner Production, 2021, 280: 124383.
[85] ORYSHCHYN D, HARUN N F, TUCKER D, et al. Fuel utilization effects on system efficiency in solid oxide fuel cell gas turbine hybrid systems [J]. Applied Energy, 2018, 228: 1953-1965.
[86] CHEN H, YANG C, ZHANG B, et al. Coupling effects of fuel reforming process and fuel utilization on the system performance of a natural gas solid oxide fuel cell/gas turbine hybrid system [J]. International Journal of Energy Research, 2021, 45(12): 17664-17690.
[87] LIU H, QIN J, LI C J, et al. Influence of fuel reforming with low water to carbon ratio on thermodynamic performance of aviation solid oxide fuel cell-gas turbine hybrid system[J]. Journal of Power Sources, 2022, 546: 231978-1-231978-12.
[88] HUANG S, YANG C, CHEN H, et al. Effects of fuel pre-reforming ratio on system design performance of a natural gas solid oxide fuel cell/gas turbine hybrid system[J]. Proceedings of the Chinese Society of Electrical Engineering, 2022, 42(10): 3691-3700.
[89] HUANG S, YANG C, CHEN H, et al. Coupling impacts of SO-FC operating temperature and fuel utilization on system net efficiency in natural gas hybrid SOFC/GT system[J]. Case Studies in Thermal Engineering, 2022, 31: 101868.
[90] CHEN H, YANG C, ZHOU N, et al. High efficiencies with low fuel utilization and thermally integrated fuel reforming in a hybrid solid oxide fuel cell gas turbine system [J]. Applied Energy, 2020, 272: 115160-1-115160-16.
[91] 包成, 蔡宁生. 固体氧化物燃料电池:燃气轮机混合发电系统建模与控制的研究现状与进展 [J]. 机械工程学报, 2008,44(2): 1-7.
BAO Cheng, CAI Ningsheng. Research status and advances in modeling and control of solid oxide fuel cell gas turbine hybrid generation system[J]. Journal of Mechanical Engineering 2008,44(2): 1-7.
[92] ROBERTS R A, ROSSI I, TRAVERSO A. Dynamic simulation of energy systems: Comparison of a physics-based against time constant based approach applied to a microturbine test rig[C]//Turbo Expo: Power for Land, Sea, and Air.[S.l.]:American Society of Mechanical Engineers, 2015:23-31.
[93] HOHLOCH M, HUBER A, AIGNER M. Experimental investigation of a SOFC/MGT hybrid power plant test rig: Impact and characterization of a fuel cell emulator [C]//Turbo Expo: Power for Land, Sea, and Air.[S.l.]:American Society of Mechanical Engineers,2016:6-18.
[94] ZHOU N N, YANG C, TUCKER D, et al. Transfer function development for control of cathode airflow transients in fuel cell gas turbine hybrid systems[J]. International Journal of Hydrogen Energy, 2015, 40(4): 1967-1979.
[95] FERRARI M L, LIESE E, TUCKER D, et al. Transient mod-eling of the NETL hybrid fuel cell/gas turbine facility and experimental validation [J]. Journal of Engineering for Gas Turbines and Power, 2007, 129(4): 1012-1019.
[96] ZHANG X, LI J, LI G, et al. Dynamic modeling of a hybrid system of the solid oxide fuel cell and recuperative gas turbine [J]. Journal of Power Sources, 2006, 163(1): 523-531.
Recent research progresses of mathematical modeling and simulation of solid oxide fuel cell/gas turbine (SOFC/GT) hybrid systems
![width=auto,height=auto,dpi=110](http://qn1.jjmjs.com.cn/upload/art/20230612-2/bf2d88bfc6392f3486c6c390c9e1e14a.png)
LIAN Yanke,MING Pingwen,CAI Liming.Recent research progresses of mathematical modeling and simulation of solid oxide fuel cell/gas turbine (SOFC/GT) hybrid systems[J].Clean Coal Technology,2023,29(3):26-39.