质子交换膜燃料电池参数敏感性的确定方法
专利汇可以提供质子交换膜燃料电池参数敏感性的确定方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种 质子交换膜 燃料 电池 参数敏感性的确定方法，包括计算模型的建立和对 燃料电池 参数敏感性的分析两个部分：第一部分分为5个步骤，包括：电池性能以 输出 电压 为指标；确定欧姆损失；确定活化损失；反应气体的浓度计算；以及电池 水 管理。第二部分为7个步骤，包括：选取待分析的参数；确定各个参数敏感性指标为；确定各个参数的取值范围；将样本空间标准化；将标准化的样本空间离散；随机取样策略；以及计算敏感性指标。能够定量确定各个工况下参数的敏感性顺序，找出电池运行时哪些参数是重要参数，哪些是次要参数，能够指导电池的设计优化和实际开发，最终为提高电池性能和 稳定性 提供依据。，下面是质子交换膜燃料电池参数敏感性的确定方法专利的具体信息内容。
1.质子交换膜燃料电池参数敏感性的确定方法，其特征是：所述确定方法包括计算模型的建立和对燃料电池参数敏感性的分析两个部分：
(1)计算模型建立的方法步骤为：
1.1电池性能以输出电压为指标
Eout＝Erev-ηohm-ηact (1-1)
其中：Eout为输出电压，Erev为可逆电压，ηohm和ηact分别为欧姆损失和活化损失，Erev可逆电压由能斯特方程求得：
其中：ΔG为吉布斯自由能变；F为法拉第常数；ΔS为熵变；R为理想气体常数；T为工况温度；Tref为参考温度； 和 分别为阳极催化层氢气压力和阴极催化层氧气压力；
1.2确定欧姆损失
欧姆损失包括极板、多孔介质层和质子交换膜造成的欧姆损失之和，即：
其中ηohm,P、ηohm,por和ηohm,m依次分别为极板、多孔介质层、和质子交换膜造成的欧姆损失；I为电流密度； 分别为流道极板和多孔介质各层传输电子的面电阻；
分别为催化层和质子交换膜内传输质子的面电阻，面电阻的求解通式如下：
Ω＝δ/σeff (1-4)
其中δ为厚度；σeff为有效电导率。
1.2.1对于扩散层或微孔层或催化层，电子的有效电导率经修正为：
式中： 表示电子的有效电导率；σs为电子固有电导率；ε为孔隙率，
1.2.2催化层和离子交换膜的质子电导率：
式中： 为催化层中质子有效电导率；Xm为催化层内电解质Nafion体积分数；σm为质子交换膜Nafion的质子电导率，
σm取决于Nafion中含水量：
其中：λ为水含量，水含量由膜内水的活度a决定：
对于催化层，水的活度为：
acl＝RH+2s (1-9)
其中RH为催化层内气体的相对湿度，s为催化层孔隙内液态水体积分数，
对于质子交换膜，水活度aaver近似等于阳极催化层和阴极催化层内水活度的平均值：
1.3确定活化损失
活化损失的解析解：
其中ηact,ano，ηact,cat分别代表阳极和阴极活化过电势；αaαc分别为阳极和阴极的电荷传输系数；n为单位反应中传输的电子数；j0,ref为参考电流密度； 分别为参考氢气浓度和参考氧气浓度； 分别为参考氢气浓度和参考氧气浓度，
1.4反应气体的浓度计算
电池内多孔介质结构中氢气、氧气扩散传输方式遵循菲克定律：
阳极和阴极各包含三个求解域，分别是扩散层、微孔层、催化层，
稳态时阳极催化层氢气扩散方程为：
其中 是微孔层、催化层交界面处氢气浓度； 是催化层、质子交换膜交界面
处氢气浓度； 为阳极催化层内平均氢气有效扩散系数，由Bruggemann修正得
δCL为催化层厚度，
阳极催化层氢气平均浓度：
稳态时阴极催化层氧气扩散方程为：
其中 是微孔层、催化层交界面处氧气浓度； 是催化层、质子交换膜交界面
处氧气浓度； 为阴极催化层内氧气有效扩散系数，
阴极催化层氧气平均浓度：
扩散层、微孔层区域的反应气体控制方程可类似列出，然后结合阳极流道内氢气浓度和阴极流道内氧气浓度的边界条件，可求得催化层内反应气体真实浓度，
1.5电池水管理
水跨膜运输方式包含电渗拖拽、膜态水扩散和压差扩散三种形式，膜内压差扩散项相对较小，模型构建中忽略了膜内压差扩散项的影响，
电渗拖拽效应表现为质子跨膜运输，同时会拖拽一定量的水从阳极到阴极，电渗拖拽系数nd为伴随每个质子由阳极到阴极跨膜的水分子数目：
膜态水扩散系数Dm的计算方法如下：
阳极催化层为水蒸气扩散，水守恒方程如下：
其中Jvap为水蒸气传输通量，方向为从流道流向膜；cvap,MPL-CL是阳极微孔层、催化层交界面水蒸气浓度；cvap,CL-PEM是催化层、质子交换膜交界面水蒸气浓度； 为催化层内水蒸气有效扩散率；ρdry为干态膜密度；EW为质子交换膜的当量质量；λacl，λccl分别为阳极和阴极催化层膜态水含量；
阴极产水，催化层内为液态水传输，水守恒方程如下：
其中ρl为液态水密度； 为水摩尔质量；sccl为阴极催化层液态水体积分数；εccl为阴极催化层孔隙率；Kl,cl为催化层水的渗透率；μl为水的动力粘度； 是阳极催化层、质子交换膜交界面的液压； 是微孔层、催化层交界面的液压；Jl为液态水流通量，扩散层，微孔层区域的水控制方程可类似列出，假设流道内无液态水，结合阳极流道内水蒸气浓度和阴极流道与扩散层交界面处液压等于一个大气压的边界条件，求得阳极各层水蒸气浓度和阴极各层交界面液压，
由Leverett方程得出多孔介质内毛细压力pc和液态水体积分数s的关系：
Pc＝Pg-Pl (1-24)
其中σlq表面张力系数；θ为多孔介质接触角，由求得的毛细压力Pc可以反求出电池内各部分液态水体积分数s，
由上述1.1-1.5五个步骤可以最终求得电池的输出电压Eout
(2)燃料电池参数敏感性分析的步骤为：
2.1选取待分析的参数：
包括电化学参数、物性参数、几何参数等，记为p1,p2,…,pi,…,pk，下标i为参数的标号，k为参数总数，则有关系如下：
Eout＝f(p1,p2,...,pk) (2-1)
2.2确定各个参数敏感性指标为：
EEi用于衡量参数i的敏感度，Δ为参数pi在取值区间的增量，为了得到整个参数范围内平均意义上的敏感度，需要在整个参数范围内随机取点，依次计算不同点处的敏感性指标，最后做平均处理作为参数的整体敏感性指标，
2.3确定各个参数的取值范围：
每个参数在各自的取值区间服从均匀分布：pi～U[ai,bi]，i取1～k，ai和bi分别是参数i的取值下限和上限，这k个参数共同确定出整个样本空间，该样本空间为k维空间，
2.4将样本空间标准化
每个参数pi经过变量替换为归一化的无量纲参数Xi，使得Xi～U[0,1]:
其中， 为参数pi的概率分布函数，它由概率密度函数 求得，无量纲处理Xi～
U[0,1]，计算实际参数和无量纲参数的对应关系如下：
pi＝Xi×(bi-ai)+ai (2-4)
处理后敏感度指标计算式为：
2.5将标准化的样本空间离散
将每个标准化的参数Xi离散为L层，则各个参数的取值集合为{0,1/(L-1),2/(L-1),…,
1}，参数的增量Δ取值集合为{1/(L-1),2/(L-1),…,1-1/(L-1)}，为使得等概率取样本空间上的样本点，分层数L一般取偶数，Δ取(L/2)/(L-1)，
2.6随机取样策略
将样本空间离散后需要在空间上随机取点，首先在样本空间上随机确定一个基准点，记为x*，接下来由该基准点生成第一个样本点x(1)，x(1)仍在样本空间内，它由x*随机朝着某个参数或者多个参数方向增加或减小一个增量Δ得到，
随后，由第一个样本点x(1)随机朝着某个参数方向i增加或减小一个增量Δ得到x(2)：x(2)＝x(1)+eiΔ或者x(2)＝x(1)-eiΔ，其中ei为单位方向向量，方向为参数i增大的方向，按照上述方法依次遍历各个参数方向得到一个样本序列x(1)，x(2)，…，x(k+1)，该序列中每相邻两个点仅在某个参数方向上增加了一个增量Δ，因此每相邻两个点可以得到某个参数的一个敏感性指标EEi，而每生成一个样本序列可以计算一次每个参数的EEi，按照上述方式随机生成r个样本序列，则可以将每个参数求解r次敏感性指标EEi，对r个EEi求平均值，从而可以得到各个参数在整个样本空间平均意义上的敏感性指标，样本矩阵B*可以由下式表示：
B*＝(Jk+1,1x*+ΔB)P* (2-6)
*
ΔB＝(Δ/2)[(2B-Jk+1,k)D+Jk+1,k] (2-7)
其中：x*为随机生成的基准点；Jk+1,1为(k+1)×1阶的全1矩阵；P*为k×k阶随机生成的置换矩阵，矩阵的每一行和每一列有且仅有一个元素为1，其它元素为0，按照置换矩阵的性质，从矩阵第一行读到最后一行，P*给定了这些参数移动的顺序；ΔB为增量矩阵，包含每个*
参数方向的增量信息；Jk+1,k为(k+1)×k阶的全1矩阵；D为随机生成的方向矩阵，是一个k维对角矩阵，且对角元素仅为1或者-1，1表示沿着对应的参数方向增加Δ，-1表示沿着对应参数方向减小Δ；B是严格下三角矩阵，且下三角部分全为1：
2.7计算敏感性指标
由2.6生成的样本矩阵B*的每一行确定了一个样本点，从矩阵第一行到k+1行依次对应样本点 分
别将其代入模型：
Eout＝f(p1,p2,...,pk)＝f(X1,X2,...,Xk) (2-9)
(1) (2) (k+1)
计算出k+1个函数值{f ,f ,…,f }，则第j个样本矩阵计算的第i个参数的敏感性指标为：
其中x(l)和x(l+1)为两个相邻样本点，它们仅在参数i方向上有增量±Δ，若x(l+1)和x(l)相比仅在参数i方向上增加了Δ，则分母取正号；若减小Δ，则分母取负号，
综合r个样本矩阵得到的各个参数的敏感性指标，可以定义敏感性系数如下：
根据上述步骤可以完成对燃料电池多个参数同时进行敏感性分析，确定这些参数的敏感性顺序，最终确定出各个电流密度下哪些参数是重要参数。
质子交换膜燃料电池参数敏感性的确定方法
技术领域
背景技术
发明内容
具体实施方式
p2,...,pk)，计算出k+1个函数值{f(1),f(2),…,f(k+1)}＝{0.990167，0.989931，0.991195，
0.991475，0.99192，0.991564，0.986391，0.986384，0.986387，0.986834，0.984845，
0.976077，0.975801}，则该样本矩阵计算的第i个参数的敏感性指标计算式为：
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 红外气体传感器浓度信号处理方法 | 2020-05-19 | 737 |
| 用于机载确定燃料挥发性的方法 | 2020-05-15 | 735 |
| 一种同时测定多温度下建材甲醛散发特性参数的阶跃温升密闭散发法 | 2020-05-22 | 880 |
| 利用电容模型的改进的阴极压力控制 | 2020-05-24 | 1006 |
| 一种获取往复式活塞压缩机阀腔压力脉动的方法 | 2020-05-14 | 762 |
| 用流体分析仪确定流体特性 | 2020-05-26 | 49 |
| 储气洞室的气密性计算方法 | 2020-05-19 | 381 |
| 以气体测量不规则物体体积的装置及方法 | 2020-05-11 | 901 |
| 一种质子交换膜燃料电池系统状态估计方法 | 2020-05-26 | 27 |
| 气体分子运动论模拟实验仪 | 2020-05-16 | 412 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
