采用常规条形流场的质子交换膜燃料电池阴极数值模拟(Ⅰ)模型建立

对采用常规条形流场的H₂-Air PEMFCs阴极建立了二维数学模型,模型的控制方程耦合了连续性方程、Darcy方程、电传导方程以及O₂和H₂O的对流-扩散方程,对氧的电化学还原反应过程采用Butler-Volmer方程描述.利用模型计算了阴极扩散层中电流密度、O₂和H₂O浓度、催化层界面上局部电流密度的分布,分析了采用常规条形流场时气体在阴极扩散层中的传递机制及各组分浓度分布的特点.     

脉冲放电NO脱除过程模拟

建立了数学模型,把所考虑的NO/N₂/O₂/H₂O脉冲放电体系中的主要化学反应,归结为各物种浓度的刚性常微分方程组的初值问题来求解,计算方法选用Rosenbrock法,得到了物种浓度随停留时间的演化规律。结果表明：NO还原和氧化在其脱除过程中同时共存,当NO,N₂,O₂,H₂O和CO₂的体积浓度分别为200×10^-6,80%, 5％,6％和9％时,NO氧化所占的比例比NO还原的大很多;脉冲频率增大导致NO还原率和氧化率均增大;H₂O浓度增大导致HNO₃浓度增大,表明NO氧化所占的比例随H₂O浓度增大而增大。     

气体与煤表面吸附作用的量子化学研究

选用褐煤、次烟煤、高挥发分烟煤、低挥发分烟煤和无烟煤5种煤表面结构模型,采用量子化学半经验方法INDO,从分子水平描述了CO、O₂、H₂O(g)、CO₂、CH₄和H₂等6种气体在煤表面的吸附作用,计算了气体在煤表面的吸附能、吸附距离、吸附作用键级和净电荷变化等微观参数,用Morse函数拟合了气体与煤表面的结合能曲线,得到了气体吸附作用强弱次序为：CO和O₂最强,H₂O和CO₂次之,CH₄和H₂最弱。     

采用常规条形流场的质子交换膜燃料电池阴极数值模拟(Ⅱ)电池性能影响因素的分析

利用已建立的数学模型考察了阴极扩散特性参数、阴极流场几何参数、电极电阻、催化剂活性、流道上O₂浓度变化对H₂-Air PEMFCs性能的影响.计算表明,增大氧有效扩散系数、减小扩散层厚度可增大电池的工作电流密度;提高氧还原反应交换电流密度、减小电极电阻、优化流道尺寸可改善电池性能;沿反应气体流动方向逐渐增加MEA的催化剂负载量可提高电流密度分布的均匀性.     

烟气气相组分及Ca(OH)_2对KMnO_4氧化NO的影响机理

在固定床反应器中考察了KMnO₄氧化烟气中NO的过程,分析了烟气组分H₂O、O₂及SO₂对NO氧化过程的影响规律,得到了Ca（OH）₂对KMnO₄氧化NO的影响机理。实验结果表明,H₂O是KMnO₄氧化NO的必要条件;在含H₂O条件下,O₂可以提高NO氧化率。SO₂与氧化剂反应生成无水钾镁钒类复盐K₂Mn₂（SO₄）₃对NO氧化具有负面作用;Ca（OH）₂的加入提高了氧化剂表面的固体碱度从而促进氧化过程进行;通过添加Ca（OH）₂可以降低SO₂对NO氧化过程的负面影响。根据气体成分和产物分析可知,KMnO₄在钙基吸收剂表面氧化烟气中NO的机理可能是KMnO₄以离子态将吸附在氧化剂表面的NO和SO₂氧化为NO₂和SO₃,生成的NO₂、SO₃再传递到氧化位临近的碱性位被吸收。     

简化PDF模型对Texaco气化炉的三维数值模拟

应用商业CFD软件Fluent建立气化过程热态模型,对某化肥厂Texaco水煤浆气化炉进行了三维数值模拟。计算中采用简化PDF模型描述炉内的化学反应,将水煤浆看作燃料流,氧气看作氧化剂流;根据对冷态流场的计算,采用六面体结构网格为主的网格划分;Realizable k-ε湍流模型封闭湍流方程;dpm模型考察气体和颗粒相耦合;随机轨道模型对颗粒相进行追踪,P-1辐射模型计算炉内辐射特性;同时编制UDF函数模拟焦炭和O₂、H₂O、CO₂以及H₂的颗粒异相反应。通过与工业数据的对比,证明该模型能够真实反应气化炉内的物理特性,同时表明工业炉内的同相反应基本达到化学平衡。     

烟气气相组分及Ca（OH）2对KMnO4氧化NO的影响机理

在固定床反应器中考察了KMnO₄氧化烟气中NO的过程,分析了烟气组分H₂O、O₂及SO₂对NO氧化过程的影响规律,得到了Ca（OH）₂对KMnO₄氧化NO的影响机理。实验结果表明,H₂O是KMnO₄氧化NO的必要条件;在含H₂O条件下,O₂可以提高NO氧化率。SO₂与氧化剂反应生成无水钾镁钒类复盐K₂Mn₂（SO₄）₃对NO氧化具有负面作用;Ca（OH）₂的加入提高了氧化剂表面的固体碱度从而促进氧化过程进行;通过添加Ca（OH）₂可以降低SO₂对NO氧化过程的负面影响。根据气体成分和产物分析可知,KMnO₄在钙基吸收剂表面氧化烟气中NO的机理可能是KMnO₄以离子态将吸附在氧化剂表面的NO和SO₂氧化为NO₂和SO₃,生成的NO₂、SO₃再传递到氧化位临近的碱性位被吸收。     

天然气-煤共气化制备合成气热态模拟

天然气-煤共气化制备通用合成气技术是基于天然气蒸汽转化法和煤气化过程开发的新工艺，通过技术原理分析，可以直接制备出H₂/CO在1.0～2.0之间可调的粗合成气.应用移动床反应器进行热态模拟实验，主要研究了不同操作参数对火焰区温度及合成气有效成分(H₂＋CO)和H₂/CO的影响.结果表明：天然气与氧气在同一位置喷入反应器，控制喷吹参数H₂O/CH₄/O₂以及流量，在气化炉炉温不低于1000 ℃的条件下，煤或焦炭中挥发分基本裂解，可以直接制备出H₂/CO在1.0～2.0之间，有效成分大于90.0%，残留的CH₄小于2.0%的粗合成气.     

工艺条件对Ni/Al2O3和Ni/CeO2-Al2O3催化剂在甲烷自热重整制氢反应中催化性能的影响

采用共沉淀法制备γ-Al₂O₃载体和不同Ce添加量的CeO₂-Al₂O₃载体,然后用浸渍法制备Ni负载质量分数10%的Ni/γ-Al₂O₃和Ni/CeO₂-Al₂O₃催化剂。在固定床微反装置中考察了反应温度、原料气配比和CH₄空速等工艺条件对Ni/γ-Al₂O₃和Ni/Ce₃₀Al₇₀O_δ催化剂在甲烷自热重整制氢反应中催化性能的影响。结果表明,添加Ce的催化剂催化性能有较大提高,在Ni/Ce₃₀Al₇₀O_δ催化剂上,反应温度750 ℃时, CH₄转化率94.3％,与Ni/Al₂O₃催化剂相比,提高20％。Ni/γ-Al₂O₃和Ni/CeO₂-Al₂O₃催化剂的CH₄转化率均随反应温度的升高而增大。原料气中n(O₂)∶n(CH₄)和n(H₂O)∶n(CH₄)的增加均能提高各催化剂的CH₄转化率。但n(O₂)∶n(CH₄)和n(H₂O)∶n(CH₄)的变化对各催化剂的催化性能的影响不同。随着n(O₂)∶n(CH₄)的增大,产物中n（H₂）∶n（CO）降低,n（CO₂）∶n(CO＋CO₂)升高;而n(H₂O)∶n(CH₄)增大时,产物中n（H₂）∶n（CO）和n（CO₂）∶n(CO＋CO₂)均升高。随着CH₄空速的增加,Ni/Al₂O₃催化剂上CH₄转化率、n（H₂）∶n（CO）和n（CO₂）∶n(CO＋CO₂)均较大程度下降;而在Ni/Ce₃₀Al₇₀O_δ催化剂上,随着CH₄空速的增加,CH₄转化率、n（H₂）∶n（CO）和n(CO₂）∶n(CO＋CO₂)变化不大。     

Cu1Zr1Ce9Oδ催化剂选择性氧化CO性能的研究

用共沉淀法制备了Cu₁Zr₁Ce₉O_δ催化剂,考察了反应温度和反应气体中各组分对Cu₁Zr₁Ce₉O_δ催化剂上选择性氧化CO反应的影响。结果表明,降温的过程中Cu₁Zr₁Ce₉O_δ催化剂的活性滞后。H₂的存在有利于CO的脱附,促进了低温下选择性氧化CO的反应;而温度较高时,H₂氧化副反应的发生降低了CO的转化率,反应气中H₂O和CO₂降低了催化剂的活性和选择性,最佳反应温度为（160~200） ℃,O₂的进入量取3为宜。     

质子交换膜燃料电池两维、两相流动模型

提出了考虑电池内部两相流动的质子交换膜燃料电池数学模型,模拟了阳极、阴极两侧的流道和扩散层中同时发生两相流动时电池内部的各种传递特性,并用实验数据验证了该模型的准确性。模拟结果显示,当电池阴极扩散层中有液态水存在时会大大降低膜中的局部电流密度;质子交换膜中水的净通量方向可正、可负,因此电池的增湿策略应根据不同的运行工况而不断变化。     

NH_3-H_3PO_4-H_2O体系汽液平衡的测定与热力学平衡计算

采用静态法测定了NH_3-H_3PO_4-H_2O体系在NH_3/H_3PO_4分子比为1.0-1.4,相应温度范围为75-175℃时的汽液平衡数据,并把实验数据关联成数学式.根据热力学理论,对NH_3-H_3PO_4-H_2O体系的热力学平衡计算的模型作了初步探讨.液相活度及活度系数采用扩充Pitzer方程,汽相逸度系数采用根据极性修改的RK方程.根据Edwards提出的弱电解质溶液的分子热力学,基于质量平衡、化学平衡、电荷平衡及汽液平衡理论,提出了NH_3-H_3PO_4-H_2O体系汽液平衡计算的热力学模型.计算结果表明:摩尔浓度低于25mol/kg H_2O时,模型计算结果与实验数据吻合得较好.     

固定床反应器中气液分配器的流体力学性能

设计了一种底部带有碎流板的新型管式气液分配器。以水和空气代替工业上的原油和氢气进行冷模实验。实验过程为:水由水箱经水泵抽出,经液体流量计计量后进入实验塔上部的气液扩散器,气液两相流经急冷箱冷却后同时向下通过分配器。进入到接液装置的液体通过橡胶管导入到放置在地面上的17个标有编号的量筒中,未进入接液装置的液体流入水箱,气体则排放到大气中。最后用U形管差压计测量分配器的压力降。实验研究了该分配器的分配性能、压力降损失、气液相操作弹性和分布不均度,优化并确定该分配器的结构形式和结构参数。结果表明,该分配器的最佳工作条件为:液相量为0.3m³/h左右,气相量为20~30m³/h。     

Modelling Approach for Planar Self‐Breathing PEMFC and Comparison with Experimental Results

This paper presents a model‐based analysis of a proton exchange membrane fuel cell (PEMFC) with a planar design as the power supply for portable applications. The cell is operated with hydrogen and consists of an open cathode side allowing for passive, self‐breathing, operation. This planar fuel cell is fabricated using printed circuit board (PCB) technology. Long‐term stability of this type of fuel cell has been demonstrated. A stationary, two‐dimensional, isothermal, mathematical model of the planar fuel cell is developed. Fickian diffusion of the gaseous components (O₂, H₂, H₂O) in the gas diffusion layers and the catalyst layers is accounted for. The transport of water is considered in the gaseous phase only. The electrochemical reactions are described by the Tafel equation. The potential and current balance equations are solved separately for protons and electrons. The resulting system of partial differential equations is solved by a finite element method using FEMLAB (COMSOL Inc.) software. Three different cathode opening ratios are realized and the corresponding polarization curves are measured. The measurements are compared to numerical simulation results. The model reproduces the shape of the measured polarization curves and comparable limiting current density values, due to mass transport limitation, are obtained. The simulated distribution of gaseous water shows that an increase of the water concentration under the rib occurs. It is concluded that liquid water may condense under the rib leading to a reduction of the open pore space accessible for gas transport. Thus, a broad rib not only hinders the oxygen supply itself, but may also cause additional mass transport problems due to the condensation of water.     

THE MODIFIED PGR EQUATION OF STATE: PURE-FLUID PREDICTIONS

The Park-Gasem-Robinson (PGR) equation of state (EOS) has been modified to improve its volumetric and equilibrium predictions. Specifically, the attractive term of the PGR equation was modified to enhance the flexibility of the model, and a new expression was developed for the temperature dependence of the attractive term in this segment-segment interaction model.

The predictive capability of the modified PGR EOS for vapor pressure and saturated liquid and vapor densities was evaluated for selected normal paraffins, normal alkenes, cyclo-paraffins, light aromatics, argon, carbon dioxide, and water. The generalized EOS constants and substance-specific characteristic parameters in the modified PGR EOS were obtained from pure component vapor pressures and saturated liquid and vapor molar volumes. The calculated phase properties were compared with those of the Peng-Robinson (PR), the simplified-perturbed-hard-chain theory (SPHCT), and the original PGR equations. Generally, the performance of the proposed EOS (%AAD of 1.3, 2.8, and 3.7 for vapor pressure, saturated liquid, and vapor densities, respectively) was better than the SPHCT and original PGR equations in predicting the pure fluid properties.     

Vapor-liquid equilibria for the binary monoethanolamine+ water and monoethanolamine+ethanol systems

Isobaric vapor-liquid equilibrium data for monocthanolamine (MEA)+water and monoethanolamine (MEA)+ ethanol were measured in an equilibrium cell circulating both the vapor and liquid phases. The measurements of the boiling temperatures and corresponding equilibrium compositions for two binary mixtures were carried out at the atmospheric pressure. The results were checked for thermodynamic consistency by the point test method. The activity coefficients for nonideality of the liquid phase were calculated by using the Wilson, NRTL, and UNIQUAC equations. The fugacity coefficients for nonideality of the vapor phase were calculated from the virial equation of state with the second virial coefficient. The experimental data were verified to be thermodynamically consistent and showed a good agreement with the calculated values.     

On the representation of liquid properties of pure polar compounds by means of a modified redlich-kwong equation of state

The temperature-dependent parameters ω_a and ω_b of the Redlich-Kwong equation of state were evaluated from saturated liquid volumes and vapor pressures for three polar compounds, water, ammonia and sulfur dioxide, and then correlated in terms of T_r. To demonstrate the applicability of the resulting correlation, the values of ω_a and ω_b computed from these equations were used to evaluate the liquid properties of these compounds, such as vapor pressures and liquid volumes at saturation, and the heat of vaporization at the normal boiling point. In addition, the enthalpy departures of liquid water as well as the liquid volumes of compressed water were calculated. The results obtained from the proposed correlations are generally superior than those obtained from the equations of state recently proposed by Fuller, and Peng and Robinson.     

PEM燃料电池阴极中的液态水及其影响因素

阴极多孔介质中液态水的含量对PEM燃料电池阴极中的传质及其性能具有极其重要的影响。提出了一个二维、两相、稳态数学模型,研究PEM燃料电池阴极中两相水的传递及其对电池性能的影响。模型耦合了连续方程、动量方程和组分守恒方程,并将质子膜中的净水迁移通量作为边界条件之一来处理。通过实验的方法和数值模拟的方法,研究了电池操作压力和温度对电池性能的影响,同时验证了模型的有效性。模拟发现:提高操作压力和升高阴极加湿温度使电池阴极催化剂层(CTL)和扩散层(GDL)界面上的液态水含量大幅提高;升高阳极加湿温度,电池阴极CTL和GDL界面上的液态水含量变化不明显;而升高燃料电池的操作温度,阴极CTL和GDL界面上液态水的含量降低。     

Numerical simulation of enhancement of mass transfer in the cathode electrode of a PEM fuel cell by magnet particles deposited in the cathode-side catalyst layer

Recent experimental results show that performance of a proton exchange membrane (PEM) fuel cell is improved when small particles of permanent magnet are deposited in the catalyst layer on the cathode side. In this study, the mechanism of this phenomenon was clarified by using numerical simulation. Permanent magnet particles induce a Kelvin (magnetic) repulsive force against liquid water and an attractive force towards oxygen gas. To precisely study the behavior of oxygen and liquid water flows near the interface between the catalyst layer and gas diffuser, the simulation domain included the catalyst layer and gas diffuser. The simulation results revealed the following. The Kelvin repulsive force against liquid water mainly “manages” the liquid water flow and improves the performance of the fuel cell, especially in the current-limited region. In the presence of Kelvin forces, the liquid water saturation near the interface between the catalyst layer and gas diffuser decreases, thus making more pore space available for transport of oxygen gas. Furthermore, the water velocity moving from the interface increases in the upstream region. The gas velocity near the interface also increases, and thus more oxygen is supplied to the reaction sites. In summary, the Kelvin force promotes removal of liquid water from the catalyst layer, thus providing more oxygen to the catalyst and improving the performance of the fuel cell.