阙加雄，杨双桥，左国坤，王蔚国

(中国科学院宁波材料技术与工程研究所燃料电池与能源技术事业部，浙江 宁波 315201)

　　摘要：仿真技术是研究固体氧化物燃料电池并加快其商业化开发周期的重要方法之一。基于固体氧化物燃料电池的基本原理，即综合考虑电池内部的质量平衡、能量平衡和电化学反应过程，运用Matlab/Simulink建立了固体氧化物燃料电池的集总模型。运用此模型仿真了当负载发生变化时，电池的温度和输出功率的动态响应情况。并对单电池模型引入了基于电流的控制策略，用以保证在负载发生变化时单电池的安全运行。仿真结果表明此类简单的控制策略并不能保证电池的正常工作温度和良好的负载追随性，因此在商业化开发类似控制系统时，必须对系统加以更为完善的控制。

　　随着能源短缺和环境问题成为本世纪全球面临的最重要课题，固体氧化物燃料电池（SOFC）技术日益受到重视。

　　SOFC技术环境友好、高效，它通过电化学反应将化石燃料中的化学能直接转化为电能，几乎不排放任何氮氧化物和硫氧化物，同时能量转换方式不受卡诺循环限制，能量利用效率也大为提高。目前SOFC技术的实用化、商业化研究被广泛开展，其中对SOFC电池系统进行有效的控制是其中一个重要环节。但如果仅通过实际的SOFC系统进行研究需要大量资金和人力的投入，在实际操作中也有很多困难，因此基于SOFC的数学模型的研究方法被广泛采用，也被证明是相当经济有效。

　　在SOFC电池内部，主要包括以下三个动力学过程：

　　(1)与燃料及空气流动以及电化学反应放热相关联的热动力学过程；

　　(2)电池内化学组分分压变化所引起的电化学反应动力学过程；

　　(3)与燃料及空气流动相关的流体力学动力学过程。

　　本文基于以上动力学过程，利用Matlab/Simulink建立了SOFC单电池的动态模型，研究了在负载发生变化时，电池的温度响应和输出功率响应，并引入了与之相应的控制策略，以保证电池的安全运行。

　　1数学模型原理

　　详尽地计算流体力学（CFD）模型可以精确地反应电堆内部各个点上的温度及温度梯度^[1]，需要耗费大量的计算时间。本文对电池模型做了适当的简化^[1]，主要考虑电池内部四个部分的温度分布：连接件、阳极气体、PEN（包括电池阳极、电解质及阴极）和阴极气体，如图1所示，这样的模型即为集总模型。

　　相耦合的复杂过程，需要将整个模型分为质量平衡、能量平衡以及电化学三个子模型进行分析。由于燃料气体为氢气，故本文所述模型仅考虑与氢气相关的化学反应及传热、传质过程，而没有考虑甲烷的重整过程。

　　1.1质量平衡子模型

　　电池内部各组分的质量平衡主要考虑阴极侧空气流道中氧气的平衡和阳极侧氢气流道中氢气和水的平衡

　　2仿真

　　2.1模型建立方法

　　Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境，可方便地进行模型的交互和扩展。本文根据以上数学模型，利用Matlab/Simulink建立了10cm×10cm阳极支撑平板式SOFC单电池的动态模型。电池的结构参数如表1。

　　图2是用Matlab/Simulink建立的单电池数学模型。

　　通过式（1）～（13）相耦合，设定SOFC单电池的工作条件，给SOFC电池入口气体成分、流速，即可对SOFC的整个工作过程进行仿真。

　　2.2控制条件

　　对于固体氧化物燃料电池，随着负载的变化，燃料利用率会发生变化。当燃料利用率超过90%时，电池会因为内部燃料气体耗尽而永久性损坏^[5]。可采用基于电流的燃料控制策略^[6]，即通过保持固定的燃料利用率来保证电池内部的燃料不致耗尽，保障电池的安全、稳定运行。

　　2.3结果与分析

　　根据建立的模型仿真了在保持燃料利用率不变和阴极空气流量不变的条件下，电池由空载变为负载电流为50A时，系统的动态响应，过程如图4所示。

　　阴极空气入口流速为5.5276L/min，阳极H2入口流速随负载电流I按照式（14）变化，温度均为1073K。加上负载后，电池及出口处尾气的温度均上升，这是由电池内部电化学反应放热及欧姆放热所致。PEN板温度与连接件温度非常接近（在图中两条曲线重合），这是因为连接件与PEN板之间有很强的热传导及热辐射。另外由仿真结果可以看出，电池温度达到稳态所需时间较长，约103数量级，且电池温度上升幅度很大，约300K。SOFC在正常工作时，温度在稳定工作点附近的波动不宜过大（～10K），故其温度需要加以控制。电池输出功率达到稳态所需要的时间比较长，说明其负载追随性能还需要加强[7]。

　　保持以上工作条件不变，阴极入口空气流速增大到22L/min，系统的动态响应过程如图5所示。可见，电池各部分的温度均明显下降，这是因为增大空气流量，气体与PEN板和连接件之间的对流换热增强，有利于降低电池的温度。同时电池温度和电池输出功率达到稳态所需的时间也明显减少。

　　因此可以考虑通过调节空气流速来调节电池温度，并改善其负载追随性。

　　3结论与讨论

　　以上建立的模型能够很好地反应负载发生变化时，固体氧化物燃料电池的动态响应。所采用的基于电流的燃料控制策略在理想情况下，能够保证电池内部的燃料气体不致耗尽。

　　但实际燃料供应系统由于存在延迟，必须采取更为完善的控制方法，保证系统的安全、稳定运行。

　　另外，在实际的固体氧化物燃料电池系统中，为了防止电池固体材料因热应力过大而损坏，对其内部的温度梯度分布亦有严格的限制。本文所建立的集总模型不能反应电池内部的温度梯度，因此需要进一步完善，将电池进行一维或二维离散，得到更为完善的模型，便于控制系统的设计。

　　同时，本文也没有从系统的角度，考虑固体氧化物燃料电池系统中的电池堆前置燃料处理、后置尾气回收处理等子系统的相互影响因素，这部分内容将在后续工作中深入研究。

　　参考文献：

　　[1]RECKNAGLE K P,WILLIFORD R E,CHICK L A.Three-dimen-sionalt hermo-fluid electro- chemical modeling of planar SOFC stacks[J].Journal of Power Sources，2002，113:109-114.

　　[2]CAMPANARI S，IORA P.Definition and sensitivity anslysis of a finite volume SOFC model for a tubular cell geometry[J].Journal of Power Sources，2004，132:113-126.

　　[3]SELIMOVIC A.Modeling of Solid Oxide Fuel Cells Applied to the Analysis of Integrated Systems with Gas Turbines[D].weden:Lund University，2002.

　　[4]AGUIAR P，ADJIMAN C S,BRANDON N P.Anode-supported intermediate temperature direct internal reforming solid oxide fuel cell.I: model-based steady-state performance[J].Journal of Power Sources，2004，38:120-136.

　　[5]PADULLE'S J，AULT G W，MCDonALD J R.An integrated SOFC plant dynamic model for power systems simulation[J].Journal of Power Source，2000，86:495-500.

　　[6]MUELLER F，BROUWER J，JABBARI F.Dynamic simulation of an integrated solid oxide fuel cell system including current-based fuel flow control[J].Journal of Fuel Cell Science and Technology,2006，3:144-155.

　　[7]FERRIARI M L，MAGISTRI L A，MASSARDO F S.Proceedings of the control system for solid oxide fuel hybrid systems[J].ASME,Reno-Tahoe，2005，6: 1-9.