CN106960970A - 质子交换膜燃料电池内部水分布预测方法 - Google Patents

质子交换膜燃料电池内部水分布预测方法 Download PDF

Info

Publication number
CN106960970A
CN106960970A CN201710098474.8A CN201710098474A CN106960970A CN 106960970 A CN106960970 A CN 106960970A CN 201710098474 A CN201710098474 A CN 201710098474A CN 106960970 A CN106960970 A CN 106960970A
Authority
CN
China
Prior art keywords
gas
diffusion layers
water
gas diffusion
represent
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710098474.8A
Other languages
English (en)
Other versions
CN106960970B (zh
Inventor
焦魁
张国宾
樊林浩
邓豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201710098474.8A priority Critical patent/CN106960970B/zh
Publication of CN106960970A publication Critical patent/CN106960970A/zh
Application granted granted Critical
Publication of CN106960970B publication Critical patent/CN106960970B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M8/00Fuel cells; Manufacture thereof
    • H01M8/04Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
    • H01M8/04298Processes for controlling fuel cells or fuel cell systems
    • H01M8/04313Processes for controlling fuel cells or fuel cell systems characterised by the detection or assessment of variables; characterised by the detection or assessment of failure or abnormal function
    • H01M8/04492Humidity; Ambient humidity; Water content
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M8/00Fuel cells; Manufacture thereof
    • H01M8/04Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
    • H01M8/04298Processes for controlling fuel cells or fuel cell systems
    • H01M8/04992Processes for controlling fuel cells or fuel cell systems characterised by the implementation of mathematical or computational algorithms, e.g. feedback control loops, fuzzy logic, neural networks or artificial intelligence
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/30Hydrogen technology
    • Y02E60/50Fuel cells

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Chemical & Material Sciences (AREA)
  • Electrochemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Chemical & Material Sciences (AREA)
  • Sustainable Energy (AREA)
  • Sustainable Development (AREA)
  • Manufacturing & Machinery (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Medical Informatics (AREA)
  • Fuzzy Systems (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Artificial Intelligence (AREA)
  • Fuel Cell (AREA)

Abstract

本发明公开了一种质子交换膜燃料电池内部水分布预测方法,用于构建三维数值模型进行求解计算,步骤包括5部分:构建计算质子交换膜及催化层中,膜态水λ分布情况的守恒方程;构建计算气体扩散层、微孔层及催化层中水分布情况的守恒方程;求解燃料电池流道内液态水分布;定义流道与气体扩散层交界面两侧质量守恒;在整体计算域内构建能量守恒方程,获得温度分布。本发明利用两相流模型对质子交换膜燃料电池流道内的气液两相流动进行模拟计算,考虑表面张力和壁面吸附作用,并通过特定方法将其与膜电极(气体扩散层、微孔层、催化层及质子交换膜)内的气液两相流动进行结合,提供了一种通过数值计算预测给定工况下单个质子交换膜燃料电池内部整体水分布的方法。

Description

质子交换膜燃料电池内部水分布预测方法
技术领域
本发明属于电化学燃料电池领域,具体涉及一种用于对燃料电池的性能进行预测计算的方法。
背景技术
质子交换膜燃料电池(PEMFC)是一种将燃料中的化学能直接转化为电能的电化学反应动力装置,具有零排放和高功率密度等优点,并被广泛认为是将来最有可能取代内燃机作为汽车动力源的动力机械。但是,目前质子交换膜燃料电池性能及寿命还存在很大的提升空间,其中水管理的优劣是影响其性能的关键因素之一。
之所以要对质子交换膜燃料电池进行水管理,是因其质子交换膜需要保持一定的水含量才能保持高性能,同时,由于质子交换膜燃料电池水仅在阴极产生,特定工况下(如:高电流密度)又有可能发生“水淹”现象,这样就会降低燃料电池性能。因此,通过燃料电池内部结构优化设计等方法保持内部动态水平衡是提升燃料电池性能及寿命的一个关键因素。
建模仿真是获取质子交换膜燃料电池内部水分布的一条重要途径,不仅可节约大量的实验成本,还可预测通过实验手段很难观测到的燃料电池实际工作时的实时水分布情况。然而,由于质子交换膜燃料电池各组成部件以及水在不同部件中存在形态的不同,一般很难同时获得电池内部整体水分布情况。目前常见的燃料电池内部水分布研究大多将燃料电池分为膜电极(其中又包括:气体扩散层(GDL)、微孔层(MPL)、催化层(CL)和质子交换膜(Membrane))和流道两部分,或忽略流道内液态水存在对膜电极水分布的影响,或仅研究流道内气液两相流动对燃料电池排水的影响,显然这两种研究方法均有其局限性,不能准确反映出质子交换膜燃料电池在特定工况下内部整体水分布情况。
本发明将提出一种基于数值计算的方法,对质子交换膜燃料电池内部电化学反应及导电过程、气液两相流动、膜吸放水和传热过程等进行综合分析,从而预测燃料电池内部水分布。
发明内容
本发明的目的是提供一种质子交换膜燃料电池内部水分布预测方法,用于对质子交换膜燃料电池内部电化学反应及导电过程、气液两相流动、膜吸放水和传热过程等进行综合分析,从而预测燃料电池内部水的分布。
质子交换膜燃料电池内部水分布预测方法,用于构建三维数值模型进行求解计算,所述方法包括5个步骤:
(1)构建计算质子交换膜及催化层中,膜态水λ分布情况的守恒方程
其中:ρmem为质子交换膜干态密度,EW为当量质量;ω为Nafion电解质体积分数;F为法拉第常数;Jm为矢量离子电流密度;为有效扩散系数;Smw为源项,由膜态水和液态水之间相变Sd-l和质子交换膜两侧压力不同导致的渗透Sp两部分组成。Smw的具体表达式为:
(2)构建计算气体扩散层、微孔层及催化层中水分布情况的守恒方程
(2.1)通过构建液压Pl和气压Pg守恒方程来获得毛细压力Pc,然后根据Leverett方程提供的多孔介质中毛细压力与液态水体积分数s的关系式获得气体扩散层、微孔层及催化层中液态水分布情况,在此过程中,考虑气体扩散层纤维束平面方向和垂直纤维束平面方向的各向异性。
液压和气压守恒方程分别为:
其中下标l和g分别代表液态水和气体;ε表示孔隙率;K表示固有渗透率;计算时,k表示相对渗透率,m表示源项,液态水和气体源项的表达式分别为:
其中Sv-l表示气态水和液态水之间的相变,Ja和Jc分别表示阳极和阴极催化层内单位体 积的电化学反应速率。Ja和Jc由Butler-Volmer方程计算得到:
其中分别表示阳极和阴极标准状态下的参考反应速率;T为温度;分别表示氢气和氧气摩尔浓度,上标ref表示参考值。αa和αc分别表示阳极和阴极电化学反应的转移系数;R为通用气体常数;ηact表示活化过电势;该电势由催化层中电子电势和离子电势计算得出。电子电势和离子电势的守恒方程分别为:
其中分别代表有效电子电导率和有效离子电导率,Se和Sion代表源项。
获得气压和液压分布后,根据Leverett方程,可求得气体扩散层、微孔层和催化层中液态水分布:
Pc=Pg-Pl 2-10
其中σ为表面张力系数;θ为多孔介质接触角;求解过程中,根据液压流量在两个平面内的不同平均固有渗透率,从而获得液态水体积分数分布:
其中Kin表示气体扩散层纤维束平面内固有渗透率;Kthrough则表示垂直纤维束平面内的固有渗透率。
(2.2)通过构建气体组分守恒方程,获得气体扩散层、微孔层和催化层中水蒸气分布:
其中Ci表示气体摩尔浓度;表示气体有效扩散系数;Si表示源项,下标i表示氢气、氧气和水蒸气。
(3)求解质子交换膜燃料电池流道内液态水分布,同时,气相内部求解组分守恒方程获得水蒸气分布。
(3.1)采用VOF模型中质量守恒方程和动量守恒方程分别为:
其中:ρ表示密度;v表示速度矢量;g表示重力加速度常数;S表示质量源项,Fs表示由表面张力引起的源项。
流道内的壁面吸附作用通过改变边界处的曲率实现:
其中分别表示壁面垂直方向和切线方向的单位向量。
然后构建液相体积分数守恒方程获得流道内液态水分布:
αgl=1 3-5
其中αg和αl分别表示气相和液相体积分数,Sl表示液相源项。
(3.2)构建气相中气体组分守恒方程获得流道内水蒸气分布:
求解过程中,气相入口边界定义为质量入口边界,出口边界定义为压力出口且出口压力恒定为一个大气压。阴阳极入口气体质量流量ma和mc分别为:
其中
其中
式中ξa和ξc分别表示阳极和阴极进气化学计量比;Iref表示参考电流密度;分别表示阳极和阴极入口面积;分别表示阳极和阴极活化面积;分别表示阳极和阴极入口压力;RHa和RHc分别表示阳极和阴极进气相对湿度。
(4)定义流道与气体扩散层交界面两侧质量守恒
在气体扩散层与流道接触面两侧定义气压、液压及气体浓度的数据交换,实现流道与气体扩散层气液两相的质量守恒。具体为将流道中获得的气压、根据气压和液相体积分数获得的液压和气体摩尔浓度赋值到气体扩散层一侧,而将气体扩散层中得到的通量转换为源项赋值到流道一侧进行计算。计算方程为:
通量F需等效为流道与气体扩散层接触面附近厚度为δ的区域内的源项:
其中S为式3-1中气液两相质量源项,Sl为式3-4中液相源项。
(5)在整体计算域内构建能量守恒方程,获得温度分布:
其中Cp,g和Cp,l分别表示气相和液相的比热容;u表示速度;keff表示有效导热系数;ST为源项。计算过程中,设置计算域边界温度为恒定值。
根据步骤(1)~(5)构建的守恒方程,建立三维数值模型进行求解计算,可最终获得设定工况下质子交换膜燃料电池流道、气体扩散层、微孔层和催化层中气态水和液态水分布以及质子交换膜和催化层中膜态水分布。
本发明所提出的对单个质子交换膜燃料电池内部水分布进行整体求解的方法,流道内采用气液两相流模型(VOF)进行求解,将表面张力和壁面吸附的影响考虑在内。在气体扩散层(GDL)、微孔层(MPL)和催化层(CL)中,采用均匀多孔介质假设,以液压、气压为求解变量,基于Leverett方程求解液态水体积分数。在此过程中,将GDL纤维束平面内和垂直纤维束平面内的各向异性考虑在内。在质子交换膜及催化层内部,求解膜态水守恒方程。在催化层内部定义膜态水与液态水或气态水的相变过程保证燃料电池内部整体水质量守恒,同时在GDL和流道接触面两侧定义一定形式的数据交换来保证气液两相的质量守恒。
本发明的特点以及带来的有益效果是:
(1)利用现有两相流模型对质子交换膜燃料电池流道内的气液两相流动进行模拟计算,考虑表面张力和壁面吸附作用,并通过特定方法将其与膜电极(气体扩散层、微孔层、催化层及质子交换膜)内的气液两相流动进行结合,提供了一种通过数值计算预测给定工况下单个质子交换膜燃料电池内部整体水分布的方法。
(2)将质子交换膜和催化层中存在的膜态水电渗拖拽作用采用对流项形式计算,更加符合实际情况。
(3)对于气体扩散层、微孔层和催化层中液态水体积分数的求解先液压和气压方程然后基于Leverett方程计算的方式获得,充分反映了液态水在不同多孔介质界面处存在的突变现象;并且提出了在气体扩散层中考虑固有渗透率各向异性时,用于计算液态水体积分数的平均固有渗透率的计算公式。
(4)流道入口质量边界条件下,计算氢气和氧气浓度时,入口压力在计算过程中根据流道内流动情况实时更新,而不是假设为一个大气压(传统做法),更能准确反映质子交换膜燃料电池流道内流动情况。
附图说明
图1质子交换膜燃料电池流道与MEA气压、液压及气体浓度数据交换示意图。
图2质子交换膜燃料电池计算域。
图3表示质子交换膜及催化层中膜态水含量分布。
图4表示气体扩散层、微孔层及催化层中液态水体积分数分布情况。
图5表示流道中液态水体积分数分布情况。
具体实施方式
以下通过具体实施例对本发明的方法步骤过程作进一步的说明,需要说明的是本实施例是叙述性的,而不是限定性的,不以此限定本发明的保护范围。
质子交换膜燃料电池内部水分布预测方法,用于构建三维数值模型进行求解计算,具体方法通过5个步骤予以实现:
(1)构建计算质子交换膜及催化层中,膜态水λ分布情况的守恒方程
质子交换膜(主要成分为Nafion电解质)在燃料电池工作过程中会吸收一定量的水来保持较高的离子传导率,以这种形态存在的水称为膜态水,同时催化层中也存在一定含量的Nafion电解质,所以膜态水λ的求解区域包括质子交换膜(Membrane)和催化层两部分。其守恒方程为:
不同于之前的膜态水守恒方程,该方程中对电渗拖拽作用采用与实际更加相符的对流项进行计算。式中:ρmem(kg m-3)为质子交换膜干态密度,EW(kg mol-1)为当量质量;ω为Nafion电解质体积分数(质子交换膜中为1);F为法拉第常数96487C mol-1;Jm(A m-2)为矢量离子电流密度;为有效扩散系数;Smw(mol m-3s-1)为源项,由膜态水和液态水之间相变源项Sd-l和质子交换膜两侧阴阳极催化层之间由压力差引起的源项Sp两部分组成。Smw的具体表达式为:
(2)构建计算气体扩散层、微孔层及催化层中水分布情况的守恒方程
(2.1)通过构建液压Pl和气压Pg守恒方程来获得毛细压力Pc。在本发明中,将质子交换膜燃料电池中的气体扩散层、微孔层及催化层均视为均匀多孔介质。由于液态水体积分数在 不同多孔介质交界面处存在突变现象,选择先求解液压方程和气压方程,然后根据Leverett方程提供的多孔介质中毛细压力与液态水体积分数s的关系式获得其内部液态水分布。同时考虑气体扩散层纤维束平面方向和垂直纤维束平面方向的各向异性,微孔层和催化层则认为是各向同性。
液压和气压守恒方程分别为:
其中下标l和g分别代表液态水和气体;ε表示孔隙率;ρ(kg m-3)表示密度,K(m2)表示固有渗透率;k表示相对渗透率,μ(kg m-1s-1)表示动力粘度,P(Pa)表示压力,m(kg m-3s-1)表示源项,液态水和气体的表达式分别为:
其中Sv-l表示气态水和液态水之间的相变,MH2、MO2和MH2O(kg mol-1)分别表示氢气、氧气和水的摩尔质量,Ja和Jc分别表示阳极和阴极催化层内单位体积的电化学反应速率,Ja和Jc由Butler-Volmer方程计算得到:
其中上标分别表示阳极和阴极标准状态下的参考反应速率;T(K)为 温度;分别表示氢气和氧气摩尔浓度,上标ref表示参考值;αa和αc分别表示阳极和阴极电化学反应的转移系数;R为通用气体常数8.314J mol-1K-1;ηact(V)表示活化过电势;该电势由催化层中电子电势和离子电势计算得出,电子电势和离子电势的守恒方程分别为:
其中电子电势守恒方程计算域为极板、气体扩散层、微孔层和催化层,离子守恒电势计算域为催化层和质子交换膜。分别代表有效电子电导率和有效离子电导率。在气体扩散层内考虑有效电导率在纤维束平面和垂直纤维束平面有效电子电导率的各向异性,Se和Sion(A m-3)代表源项,分别为。
求解电子电势时,在阳极极板表面定义总体过电势ηtotal(V),即可逆电势Er与输出电压Vout之间差值,在阴极极板表面定义参考电势0V。
ηtotal=Er-Vout
获得气压和液压分布后,根据Leverett方程,可求得气体扩散层、微孔层和催化层中液态水分布:
Pc=Pg-Pl 2-10
其中σ(N m-1)为表面张力系数;θ为多孔介质接触角。由于气体扩散层纤维束平面内 和垂直纤维束平面内固有渗透率的不同,在求解s时根据液压流量在两个平面内的不同计算平均固有渗透率。
在气体扩散层内计算气压和液压分布时,考虑固有渗透率在纤维束平面(in-plane)和垂直纤维束平面(through-plane)内的各向异性,表达式为:
其中Kin表示气体扩散层纤维束平面内固有渗透率;Kthrough则表示垂直纤维束平面内的固有渗透率。
(2.2)通过构建气体组分守恒方程,建立气体扩撒层、微孔层和催化层中氢气、氧气和水蒸气守恒方程:
其中Ci(mol m-3)表示气体摩尔浓度;表示气体有效扩散系数;计算过程中考虑其在气体扩散层中纤维束平面和垂直纤维束平面各向异性,Si(mol m-3s-1)表示源项,下标i表示氢气、氧气和水蒸气。
(3)采用VOF两相流模型求解获得质子交换膜燃料电池流道内液态水分布。求解过程 中,将表面张力和流道壁面吸附作用考虑在内,同时求解气相内组分守恒方程获得流道内反应气体及水蒸气分布。
(3.1)VOF模型中质量守恒方程和动量守恒方程分别为:
其中:ρ表示密度;v(m s-1)表示速度矢量;g表示重力加速度常数9.8m s-2;S(kg m-3s-1)表示质量源项,Fs表示由表面张力引起的源项。
其中:I表示单位矩阵,流道内的壁面吸附作用通过改变边界处的曲率实现;分别表示壁面垂直方向和切线方向的单位向量。
流道内的壁面吸附作用通过改变边界处的曲率实现:
然后构建液相体积分数守恒方程获得流道内液态水分布:
αgl=1 3-5
其中αg和αl分别表示气相和液相体积分数,Sl(kg m-3s-1)表示液相源项,
(3.2)构建气相中气体组分守恒方程获得流道内水蒸气分布:
气相气体组分守恒方程为:
求解过程中,气相入口边界定义为质量入口边界,出口边界定义为压力出口且出口压力恒定为一个大气压,阴阳极入口气体质量流量ma和mc分别为:
其中
其中
式中ξa和ξc分别表示阳极和阴极进气化学计量比;Iref(A m-2)表示参考电流密度;分别表示阳极和阴极入口面积;分别表示阳极和阴极活化面积; 分别表示阳极和阴极入口压力(在计算过程中根据流道内流动情况实时更新);RHa和RHc分别表示阳极和阴极进气相对湿度,
(4)定义流道与气体扩散层交界面两侧质量守恒
由于电极和流道内气液两相分布求解方式的不同,为同时获得质子交换膜燃料电池一定工作条件下的水分布情况,需要保证流道与气体扩散层气液两相的质量守恒。本实施例通过在气体扩散层与流道接触面两侧定义气压、液压及气体浓度的数据交换实现:具体做法为将流道一侧的气体压力、液体压力(通过气体压力与流道内液相体积分数经Leverett方程计算得到)和各气体组分浓度赋值到气体扩散层一侧作为边界条件,然后将气体扩散层侧的通量F(包括对流动量和扩散通量两部分)返回流道一侧,具体如附图1所示。计算方程为:
数据交换过程中,气体扩散层一侧求解得到的通量F需等效为气体扩散层与流道交界面附近(流道一侧)厚度为δ的区域内的源项进行计算:
(5)在整体计算域内构建能量守恒方程,获得温度分布:
在质子交换膜燃料电池工作过程中,会不断产生热量,导致各部分温度分布并不相同,而上述计算过程中很多参数与温度有关。因此,本发明在整体计算域内求解能量守恒方程以获得温度分布:
其中Cp,g和Cp,l(J kg-1K-1)分别表示气相和液相的比热容;u(m s-1)表示速度;keff(W m-1K-1)表示有效导热系数。计算时考虑其在气体扩散层中纤维束平面内和垂直纤维束平面内的各向异性。在此过程中,计算域边界温度均设置为恒定值。ST(W m-3)为源项,表达式如下所示:
其中ΔSa和ΔSc(J mol-1K-1)分别为阳极和阴极电化学反应熵变,h(J mol-1)为水相变过程中潜热值。
根据步骤(1)~(5)构建的守恒方程,建立三维数值模型进行求解计算,可最终获得设定工况下质子交换膜燃料电池流道、气体扩散层、微孔层和催化层中气态水和液态水分布以及质子交换膜和催化层中膜态水分布。
具体计算实施例
计算给定输出电压、进气加湿程度等工况下质子交换膜燃料电池内部水分布时,首先需要构建计算域,包括阴阳极极板(BP)、流道(Channel)、气体扩散层、微孔层、催化层以及质子交换膜,如附图2所示。选取流道出入口截面宽0.8×10-3m,高1.0×10-3m,极板总高度和总宽度均为1.5×10-3m,GDL和MPL厚度分别为1.9×10-4和2.0×10-5m,阳极和阴极催化层厚度分别为5.0×10-6和1.0×10-5m,计算域总长0.01m,活化面积1.5×10-5m2
质子交换膜及催化层中膜态水计算:
本发明选用Nafion 212质子交换膜,厚度为5.08×10-5m,干态密度ρmem=1980kgm-3,当量质量EW=1.1kg mol-1,膜渗透率Kmem=2.0×10-20m2,催化层中电解质含量ω=0.21。
利用Bruggemann公式对膜态水扩散系数进行修正,修正系数为1.5:
膜态水与液态水相变源项Sd-l(kg m-3s-1)表达式如下所示:
式中γd-l(s-1)表示相变率,λeq表示平衡状态下的膜态水含量,表达式如下:
其中a表示水活度:
Psat(Pa)为水蒸气饱和蒸气压,是温度的单值函数:
阴阳极催化层中因压差引起的源项Sp(kg m-3s-1)表达式为:
其中:分别为阳极和阴极催化层平均液压,μl(kg m-1s-1)为液态水动力粘度,δmem和δcl(m)分别为膜和催化层厚度。
气体扩散层、微孔层和催化层:
气体扩散层、微孔层和催化层孔隙率ε依次为0.7、0.5和0.3,接触角θ依次为120°、120°、100°,固有渗透率K依次为1.0×10-12、1.0×10-12、1.0×10-13m2,表面张力系数σ=0.625N m-1
相对渗透率表达式为:
kl=s3、kg=(1-s)3
动力粘度系数表达式为:
μl=2.414×10-5×10247.8/(T-140)
水蒸气和液态水之间的相变源项Sv-l(kg m-3s-1)表达式为:
其中Kv-l(s-1)表示水蒸气液化速率,Kl-v(s-1)表示液态水蒸发速率,在本发明中,均取100s-1,Csat(mol m-3)表示饱和态水蒸气浓度,由饱和蒸气压计算得到:
计算离子电势和电子电势时,电化学反应参数如下:阴阳极参考电化学反应速率阴阳极反应气体参考浓度 转移系数αa=αc=0.5,输出电压Vout=0.6V。
有效电导率计算表达式如下所示,其中微孔层和催化层经Bruggemann修正,在气体扩散层内考虑纤维束平面和垂直纤维束平面电子电导率的各向异性。
气体有效扩散系数计算表达式如下所示(渗流阈值εp=0.11):
(3)流道内部:
考虑如下工况,流道入口阴阳极进气均完全加湿,即RH=100%,阴阳极出口压力均为1atm,阳极进气化学计量比为ξa=5.0,阴极化学计量比ζc=2.0,参考电流密度Iref=15000A m-2,温度T=353.15K。
(4)流道与气体扩散层交界面
计算过程中,质子交换膜燃料电池膜电极和流道的计算独立进行,将流道中计算得到的气压、气体摩尔浓度和根据气压以及液相体积分数获得的液压赋值到气体扩散层一侧,并作为膜电极计算过程中的边界条件,膜电极中计算完成后将气体扩散层中得到的通量以源项的形式赋值到流道一侧作为流道中计算过程中的源项,通过数据交换,可以保证质子交换膜燃料电池流道和膜电极质量守恒。
(5)整个计算域:
求解整体计算域能量方程时,定义极板、气体扩散层、微孔层、催化层和质子交换膜比热容分别为:1580、568、3300、3300、833J kg-1K-1,极板、微孔层、催化层和质子交换膜有效导热系数分别为20、1.0、1.0和0.95W m-1K-1,在气体扩散层内,纤维束平面内导热系数为21W m-1K-1,垂直纤维束平面内为1.7W m-1K-1
同时,阳极和阴极熵值变化分别为ΔSa=130.68J mol-1K-1,ΔSc=32.55J mol-1K-1,膜态水与液态水之间相变潜热及水蒸气与液态水之间相变潜热均为h=40650J mol-1
根据上述方法,通过构建三维数值模型可计算得到一定工况下质子交换膜燃料电池内部水分布情况,部分结果如附图3、4、5所示。

Claims (1)

1.质子交换膜燃料电池内部水分布预测方法,用于构建三维数值模型进行求解计算,其特征是:所述方法步骤包括:
(1)构建计算质子交换膜及催化层中,膜态水λ分布情况的守恒方程
其中:ρmem为质子交换膜干态密度,EW为当量质量;ω为Nafion电解质体积分数;F为法拉第常数;Jm为矢量离子电流密度;为有效扩散系数;Smw为源项,由膜态水和液态水之间相变Sd-l和质子交换膜两侧压力不同导致的渗透Sp两部分组成,Smw的具体表达式为:
(2)构建计算气体扩散层、微孔层及催化层中水分布情况的守恒方程
(2.1)通过构建液压Pl和气压Pg守恒方程来获得毛细压力Pc,然后根据Leverett方程提供的多孔介质中毛细压力与液态水体积分数s的关系式获得气体扩散层、微孔层及催化层中液态水分布情况,
液压和气压守恒方程分别为:
其中下标l和g分别代表液态水和气体;ε表示孔隙率;K表示固有渗透率;计算时,k表示相对渗透率,m表示源项,液态水和气体源项的表达式分别为:
其中Sv-l表示气态水和液态水之间的相变,Ja和Jc分别表示阳极和阴极催化层内单位体积的电化学反应速率,Ja和Jc由Butler-Volmer方程计算得到:
其中分别表示阳极和阴极标准状态下的参考反应速率;T为温度;分别表示氢气和氧气摩尔浓度,上标ref表示参考值;αa和αc分别表示阳极和阴极电化学反应的转移系数;R为通用气体常数;ηact表示活化过电势;该电势由催化层中电子电势和离子电势计算得出,电子电势和离子电势的守恒方程分别为:
其中分别代表有效电子电导率和有效离子电导率,Se和Sion代表源项,
获得气压和液压分布后,根据Leverett方程,可求得气体扩散层、微孔层和催化层中液态水分布:
Pc=Pg-Pl 2-10
其中σ为表面张力系数;θ为多孔介质接触角;求解过程中,根据液压流量在两个平面内的不同平均固有渗透率,从而获得液态水体积分数分布:
其中Kin表示气体扩散层纤维束平面内固有渗透率;Kthrough则表示垂直纤维束平面内的固有渗透率,
(2.2)通过构建气体组分守恒方程,获得气体扩散层、微孔层和催化层中水蒸气分布:
其中Ci表示气体摩尔浓度;表示气体有效扩散系数;Si表示源项,下标i表示氢气、氧气和水蒸气,
(3)求解质子交换膜燃料电池流道内液态水分布,同时,气相内部求解组分守恒方程获得水蒸气分布,
(3.1)采用VOF模型中质量守恒方程和动量守恒方程分别为:
其中:ρ表示密度;v表示速度矢量;g表示重力加速度常数;S表示质量源项,Fs表示由表面张力引起的源项,
流道内的壁面吸附作用通过改变边界处的曲率实现:
其中分别表示壁面垂直方向和切线方向的单位向量,
然后构建液相体积分数守恒方程获得流道内液态水分布:
αgl=1 3-5
其中αg和αl分别表示气相和液相体积分数,Sl表示液相源项,
(3.2)构建气相中气体组分守恒方程获得流道内水蒸气分布:
求解过程中,气相入口边界定义为质量入口边界,出口边界定义为压力出口且出口压力恒定为一个大气压,阴阳极入口气体质量流量ma和mc分别为:
其中
其中
式中ξa和ξc分别表示阳极和阴极进气化学计量比;Iref表示参考电流密度;分别表示阳极和阴极入口面积;分别表示阳极和阴极活化面积;分别表示阳极和阴极入口压力;RHa和RHc分别表示阳极和阴极进气相对湿度,
(4)定义流道与气体扩散层交界面两侧质量守恒
在气体扩散层与流道接触面两侧定义气压、液压及气体浓度的数据交换,实现流道与气体扩散层气液两相的质量守恒,具体为将流道中获得的气压、根据气压和液相体积分数获得的液压和气体摩尔浓度赋值到气体扩散层一侧,而将气体扩散层中得到的通量转换为源项赋值到流道一侧进行计算,计算方程为:
通量F需等效为流道与气体扩散层接触面附近厚度为δ的区域内的源项:
其中S为式3-1中气液两相质量源项,Sl为式3-4中液相源项。
(5)在整体计算域内构建能量守恒方程,获得温度分布:
其中Cp,g和Cp,l分别表示气相和液相的比热容;u表示速度;keff表示有效导热系数;ST为源项,计算过程中,设置计算域边界温度为恒定值,
根据步骤(1)~(5)构建的守恒方程,建立三维数值模型进行求解计算,可最终获得设定工况下质子交换膜燃料电池流道、气体扩散层、微孔层和催化层中气态水和液态水分布以及质子交换膜和催化层中膜态水分布。
CN201710098474.8A 2017-02-23 2017-02-23 质子交换膜燃料电池内部水分布预测方法 Active CN106960970B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710098474.8A CN106960970B (zh) 2017-02-23 2017-02-23 质子交换膜燃料电池内部水分布预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710098474.8A CN106960970B (zh) 2017-02-23 2017-02-23 质子交换膜燃料电池内部水分布预测方法

Publications (2)

Publication Number Publication Date
CN106960970A true CN106960970A (zh) 2017-07-18
CN106960970B CN106960970B (zh) 2019-09-03

Family

ID=59481545

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710098474.8A Active CN106960970B (zh) 2017-02-23 2017-02-23 质子交换膜燃料电池内部水分布预测方法

Country Status (1)

Country Link
CN (1) CN106960970B (zh)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107480372A (zh) * 2017-08-11 2017-12-15 天津大学 质子交换膜燃料电池阳极尾气循环模型建立方法
CN107490769A (zh) * 2017-07-28 2017-12-19 电子科技大学 一种燃料电池水分布可视化在线检测系统
CN108763737A (zh) * 2018-05-25 2018-11-06 天津大学 质子交换膜燃料电池全电池三维两相模型建立方法
CN108875183A (zh) * 2018-06-06 2018-11-23 天津大学 利用膜加湿器实现燃料电池进气加湿的建模方法
CN109145483A (zh) * 2018-09-04 2019-01-04 宁夏冰核科技有限公司 基于三对角矩阵算法的燃料电池高维模型实时仿真技术
CN109143087A (zh) * 2018-07-27 2019-01-04 武汉理工大学 一种基于pemfc阴极流道传质性能的评价方法
CN109827866A (zh) * 2019-02-28 2019-05-31 东南大学 一种质子交换膜燃料电池阴极侧含水量的软测量方法
CN109830713A (zh) * 2019-01-16 2019-05-31 清华大学 燃料电池水含量检测方法、装置及计算机可读存储介质
CN109841879A (zh) * 2019-01-16 2019-06-04 清华大学 燃料电池水含量估计系统、方法、计算机设备及存储介质
CN109921069A (zh) * 2017-12-12 2019-06-21 中国科学院大连化学物理研究所 一种直接液体燃料电池阴极水含量的测定方法
CN110413941A (zh) * 2019-07-26 2019-11-05 西安交通大学 燃料电池输入输出特性的相似原理分析方法
CN111180767A (zh) * 2020-01-07 2020-05-19 天津大学 具有流体不均匀分配效应的燃料电池堆建模方法
CN111199122A (zh) * 2020-01-02 2020-05-26 西南交通大学 基于多物理场的质子交换膜燃料电池水淹故障诊断方法
CN112909303A (zh) * 2020-12-21 2021-06-04 天津大学 质子交换膜燃料电池瞬态实时模型的建立方法
CN113889646A (zh) * 2021-09-23 2022-01-04 中汽创智科技有限公司 燃料电池催化层的性能分析方法、装置、设备及存储介质
CN113991147A (zh) * 2021-10-27 2022-01-28 广东省武理工氢能产业技术研究院 一种质子交换膜燃料电池的快速活化系统
CN114050294A (zh) * 2021-11-08 2022-02-15 西安交通大学 基于OpenFOAM平台质子交换膜氢燃料电池瞬态工况数值模拟方法
CN114186438A (zh) * 2021-12-27 2022-03-15 天津大学 用于制氢的质子交换膜电解池性能预测模型的建立方法
CN114447378A (zh) * 2022-02-10 2022-05-06 温州大学 一种质子交换膜燃料电池的参数优化方法
CN114580322A (zh) * 2022-05-06 2022-06-03 中汽研新能源汽车检验中心(天津)有限公司 一种燃料电池水传输电拖拽效应的离散化建模方法
CN114707365A (zh) * 2022-06-02 2022-07-05 中汽研新能源汽车检验中心(天津)有限公司 一种燃料电池气体扩散层水-热-质传输仿真方法
CN114843555A (zh) * 2022-05-09 2022-08-02 中汽创智科技有限公司 一种燃料电池气体扩散层的仿真方法及系统
CN116207304A (zh) * 2023-04-27 2023-06-02 中汽数据有限公司 燃料电池用膜加湿器建模方法、膜加湿器模块及模型

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007323866A (ja) * 2006-05-30 2007-12-13 Toyota Motor Corp 燃料電池システム及びその制御方法
CN101536230A (zh) * 2006-09-07 2009-09-16 住友化学株式会社 单电池的耐久性评价方法、耐久性评价装置、耐久性评价程序以及燃料电池的单电池
CN104009247A (zh) * 2014-05-04 2014-08-27 华中科技大学 一种固体氧化物燃料电池电堆温度分布估计方法
CN104040772A (zh) * 2011-12-13 2014-09-10 利兹大学 质子交换膜燃料电池

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007323866A (ja) * 2006-05-30 2007-12-13 Toyota Motor Corp 燃料電池システム及びその制御方法
CN101536230A (zh) * 2006-09-07 2009-09-16 住友化学株式会社 单电池的耐久性评价方法、耐久性评价装置、耐久性评价程序以及燃料电池的单电池
CN104040772A (zh) * 2011-12-13 2014-09-10 利兹大学 质子交换膜燃料电池
CN104009247A (zh) * 2014-05-04 2014-08-27 华中科技大学 一种固体氧化物燃料电池电堆温度分布估计方法

Cited By (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107490769A (zh) * 2017-07-28 2017-12-19 电子科技大学 一种燃料电池水分布可视化在线检测系统
CN107480372A (zh) * 2017-08-11 2017-12-15 天津大学 质子交换膜燃料电池阳极尾气循环模型建立方法
CN107480372B (zh) * 2017-08-11 2019-09-03 天津大学 质子交换膜燃料电池阳极尾气循环模型建立方法
CN109921069A (zh) * 2017-12-12 2019-06-21 中国科学院大连化学物理研究所 一种直接液体燃料电池阴极水含量的测定方法
CN109921069B (zh) * 2017-12-12 2021-03-30 中国科学院大连化学物理研究所 一种直接液体燃料电池阴极水含量的测定方法
CN108763737A (zh) * 2018-05-25 2018-11-06 天津大学 质子交换膜燃料电池全电池三维两相模型建立方法
CN108875183A (zh) * 2018-06-06 2018-11-23 天津大学 利用膜加湿器实现燃料电池进气加湿的建模方法
CN108875183B (zh) * 2018-06-06 2023-04-18 天津大学 利用膜加湿器实现燃料电池进气加湿的建模方法
CN109143087A (zh) * 2018-07-27 2019-01-04 武汉理工大学 一种基于pemfc阴极流道传质性能的评价方法
CN109145483A (zh) * 2018-09-04 2019-01-04 宁夏冰核科技有限公司 基于三对角矩阵算法的燃料电池高维模型实时仿真技术
CN109841879A (zh) * 2019-01-16 2019-06-04 清华大学 燃料电池水含量估计系统、方法、计算机设备及存储介质
CN109830713A (zh) * 2019-01-16 2019-05-31 清华大学 燃料电池水含量检测方法、装置及计算机可读存储介质
CN109827866A (zh) * 2019-02-28 2019-05-31 东南大学 一种质子交换膜燃料电池阴极侧含水量的软测量方法
CN109827866B (zh) * 2019-02-28 2022-05-10 东南大学 一种质子交换膜燃料电池阴极侧含水量的软测量及闭环控制方法
CN110413941A (zh) * 2019-07-26 2019-11-05 西安交通大学 燃料电池输入输出特性的相似原理分析方法
CN111199122A (zh) * 2020-01-02 2020-05-26 西南交通大学 基于多物理场的质子交换膜燃料电池水淹故障诊断方法
CN111199122B (zh) * 2020-01-02 2023-05-05 西南交通大学 基于多物理场的质子交换膜燃料电池水淹故障诊断方法
CN111180767A (zh) * 2020-01-07 2020-05-19 天津大学 具有流体不均匀分配效应的燃料电池堆建模方法
CN111180767B (zh) * 2020-01-07 2022-03-08 天津大学 具有流体不均匀分配效应的燃料电池堆建模方法
CN112909303A (zh) * 2020-12-21 2021-06-04 天津大学 质子交换膜燃料电池瞬态实时模型的建立方法
CN113889646A (zh) * 2021-09-23 2022-01-04 中汽创智科技有限公司 燃料电池催化层的性能分析方法、装置、设备及存储介质
CN113889646B (zh) * 2021-09-23 2023-09-22 中汽创智科技有限公司 燃料电池催化层的性能分析方法、装置、设备及存储介质
CN113991147A (zh) * 2021-10-27 2022-01-28 广东省武理工氢能产业技术研究院 一种质子交换膜燃料电池的快速活化系统
CN113991147B (zh) * 2021-10-27 2023-12-05 广东省武理工氢能产业技术研究院 一种质子交换膜燃料电池的快速活化方法
CN114050294B (zh) * 2021-11-08 2023-04-14 西安交通大学 基于OpenFOAM平台质子交换膜氢燃料电池瞬态工况数值模拟方法
CN114050294A (zh) * 2021-11-08 2022-02-15 西安交通大学 基于OpenFOAM平台质子交换膜氢燃料电池瞬态工况数值模拟方法
CN114186438A (zh) * 2021-12-27 2022-03-15 天津大学 用于制氢的质子交换膜电解池性能预测模型的建立方法
CN114447378B (zh) * 2022-02-10 2023-11-07 温州大学 一种质子交换膜燃料电池的参数优化方法
CN114447378A (zh) * 2022-02-10 2022-05-06 温州大学 一种质子交换膜燃料电池的参数优化方法
CN114580322A (zh) * 2022-05-06 2022-06-03 中汽研新能源汽车检验中心(天津)有限公司 一种燃料电池水传输电拖拽效应的离散化建模方法
CN114843555A (zh) * 2022-05-09 2022-08-02 中汽创智科技有限公司 一种燃料电池气体扩散层的仿真方法及系统
CN114843555B (zh) * 2022-05-09 2024-04-09 中汽创智科技有限公司 一种燃料电池气体扩散层的仿真方法及系统
CN114707365A (zh) * 2022-06-02 2022-07-05 中汽研新能源汽车检验中心(天津)有限公司 一种燃料电池气体扩散层水-热-质传输仿真方法
CN116207304B (zh) * 2023-04-27 2023-08-01 中汽数据有限公司 燃料电池用膜加湿器建模方法、膜加湿器模块及模型
CN116207304A (zh) * 2023-04-27 2023-06-02 中汽数据有限公司 燃料电池用膜加湿器建模方法、膜加湿器模块及模型

Also Published As

Publication number Publication date
CN106960970B (zh) 2019-09-03

Similar Documents

Publication Publication Date Title
CN106960970B (zh) 质子交换膜燃料电池内部水分布预测方法
CN108763737B (zh) 质子交换膜燃料电池全电池三维两相模型建立方法
CN106848351B (zh) 质子交换膜燃料电池性能预测模型建立的方法
CN112909303B (zh) 质子交换膜燃料电池瞬态实时模型的建立方法
CN106951677B (zh) 质子交换膜燃料电池瞬态过程模型建立的方法
Shao et al. Comparison of self-humidification effect on polymer electrolyte membrane fuel cell with anodic and cathodic exhaust gas recirculation
Fontana et al. Study of the effects of flow channel with non-uniform cross-sectional area on PEMFC species and heat transfer
Chen et al. Heat and mass transfer of a planar membrane humidifier for proton exchange membrane fuel cell
CN107480372B (zh) 质子交换膜燃料电池阳极尾气循环模型建立方法
Liu et al. Asymptotic analysis for the inlet relative humidity effects on the performance of proton exchange membrane fuel cell
Chen et al. Numerical study on water transfer characteristics under joint effect of placement orientation and flow channel size for PEMFC with dead-ended anode
CN107681181A (zh) 一种燃料电池的性能诊断方法
CN110429306A (zh) 一种燃料电池水平衡测试设备和方法
Wang et al. Performance enhancement of polymer electrolyte membrane fuel cells with a hybrid wettability gas diffusion layer
CN111914414A (zh) 质子交换膜电解池三维全电池模型的建立方法
CN108875166A (zh) 采用电化学氢气泵回收燃料电池阳极氢气的建模方法
Kitahara et al. Water vapor exchange system using a hydrophilic microporous layer coated gas diffusion layer to enhance performance of polymer electrolyte fuel cells without cathode humidification
Wang et al. Comfort index evaluating the water and thermal characteristics of proton exchange membrane fuel cell
Su et al. CFD investigating the effects of different operating conditions on the performance and the characteristics of a high-temperature PEMFC
CN108875183A (zh) 利用膜加湿器实现燃料电池进气加湿的建模方法
Ouyang et al. Energy, exergy and anti-vibration assessment of microfluidic fuel cell with a novel two-phase flow model
Masaeli et al. Performance studies of a membrane-based water and heat exchanger using serpentine flow channels for polymer electrolyte membrane fuel cell application
Pan et al. Control-oriented modeling of gas purging process on the cathode of polymer electrolyte membrane fuel cell during shutting down
Song et al. AI-based proton exchange membrane fuel cell inlet relative humidity control
Zhang et al. Three-Dimensional Multi-Scale Simulation for Large-Scale Proton Exchange Membrane Fuel Cell

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant