质子交换膜燃料电池(proton exchange membrane fuel cell, PEMFC)是氢能利用的一个重要方式，作为一种能量转换器，具有运行温度低、功率密度高、启动速度快、系统鲁棒性好和无污染等优点。目前，PEMFC已成功应用于汽车、便携式电源、固定式充电站等，甚至一些公司已经实现量产并达到了一定程度的商业化应用，例如日本丰田的Mirai。然而，除基础设施建设以外，PEMFC的发展还受到成本、寿命、稳定性等的制约。

为提高PEMFC的性能和解决其面临的各种难题，国内外大量学者对PEMFC的内部流场进行了研究^[1]。X F Wang等^[2]利用多孔电极中的惯性效应，通过在流道中加入阻块，并根据阻块的平行放置原则和交错放置原则设计了两种流场，分别采用平行梯形挡板和交错梯形挡板。结果表明，采用交错梯形挡板的流场可以在相邻通道之间形成稳定的压降，从而驱动反应物实现肋下对流, 流场肋下对流的增强能提高扩散层内液态水的去除率。吴明格^[3]通过在交指流场中添加排水流道，增强了排水能力。这种主动排水流场可以很好地解决交指流场的水管理问题，但是流场内不同区域的电流密度出现了较大的不均匀性，这可能是由于压降较大或者是扩散层的厚度及有效孔隙率降低引起的。C T Wang等^[4]在流场中加入二级流道，并根据二级流道的密度提出了两种生物特征流场板。Chen等^[5]研究了在叶脉型仿生流场中加入二级流道后的电池性能，仿真结果显示，改进后的仿生流场与传统的平行流场相比，具有更均匀的气体分布、更高的气体消耗比，产生更少的Ohm热。李子君等^[6]提出了一种有利于去除气体扩散层内部水分的波形流道。Hassan^[7]引入了一种波浪流场，并与平行流场、蛇形流场、交指流场进行了实验对比，利用电化学阻抗谱分析了电池运行过程中的各项损失。高攀^[8]以羽状叶脉为原型，以Murray定律为基础，设计了一种仿生流场。相较于传统流场，该仿生流场更有利于水气在电池内的扩散传输，使水气在电池内的分布更加均匀。肖勇^[9]从发明问题解决理论(theory of inventive problem solving, TRIZ)出发，采用增加非对称性原理设计了一种仿生流场，该仿生流场的非对称性二级分支结构很好地改善了组分在催化层界面的分布。Jang等^[10]提出了一种螺旋通道，这种螺旋通道产生的二次涡增强了弯曲通道内的传热传质，显著提高了燃料电池的性能。Lian等^[11]综合分析了质子交换膜燃料电池的极化和功率密度曲线、氧分布、水分布和压降，提出了基于叶脉结构的质子交换膜燃料电池流场设计的最佳角度为15°~30°。虽然已有一些有关仿生流场的研究，但事实上目前的仿生流场尚未完全发挥出其潜力^[12]。

本文通过观察银杏叶的叶脉形状(如图 1所示)，提出了一种双银杏叶形仿生流场。由于实验研究只能从宏观的角度对流场性能进行研究，而不能对流场内部参数进行定点定性研究，更难以达到参数可视化的效果，为尽可能地对流场内部参数进行研究，本文采用三维数值模拟的方法，从内部参数的角度出发，将所提出的仿生流场与五蛇流场、平行流场进行性能对比。

1 计算模型 1.1 几何模型

基于银杏叶脉结构设计的仿生流场(以下简称为仿生流场)以及传统的五蛇流场、平行流场示意图如图 2所示。可以看出，双银杏叶脉结构仿生流场主要由主流道和众多子流道两部分组成。主流道连接流场入口和出口。子流道用于将反应气体输送到远离主流道的区域。主流道与子流道之间的夹角为45°。所有流场的活化面积均为50 mm×50 mm，每个流道的宽度和深度都为1 mm，肋宽也均为1 mm。基于计算流体动力学商业软件Fluent的单元模型的计算域均由集流板、流道、气体扩散层(gas diffusion layer, GDL)、催化剂层(catalyst layer, CL)、质子交换膜(proton exchange membrane, MEM)组成。

1.2 模型假设

1) PEMFC在等温、多相和稳态条件下工作；

2) 气体在通道中处于理想、不可压缩、层流状态；

3) 气体扩散层、催化剂层为均质多孔介质，孔隙度均匀且各向同性。

1.3 控制方程

连续性方程：

$ \frac{\partial(\varepsilon \rho)}{\partial t}+\nabla \cdot(\varepsilon \rho v)=0. $

(1)

动量守恒方程：

$ \frac{\partial(\varepsilon \rho \boldsymbol{v})}{\partial t}+\nabla \cdot(\varepsilon \rho \boldsymbol{v} \boldsymbol{v})=-\varepsilon \nabla p+\nabla \cdot(\varepsilon \mu \nabla \boldsymbol{v})+S_{\mathrm{u}} . $

(2)

能量守恒方程：

$ \frac{\partial}{\partial t}\left(\varepsilon \rho c_p T\right)+\nabla \cdot\left(\varepsilon \rho c_p \boldsymbol{v} T\right)=\nabla \cdot\left(k^{\mathrm{eff}} \nabla T\right)+S_{\mathrm{h}}. $

(3)

物种守恒方程：

$ \frac{\partial\left(\varepsilon c_i\right)}{\partial t}+\nabla \cdot\left(\varepsilon \boldsymbol{v} c_i\right)=\nabla \cdot\left(D_i^{\mathrm{eff}} \nabla c_i\right)+S_i $

(4)

上述各式中：ε为多孔介质孔隙率；ρ为密度，kg/m³；v为速度矢量，m/s；t为时间，s; p为压强，Pa；μ为黏度，Pa·s；c_p为定压比热容，J/(kg·K)；T为温度，K；k^eff为有效导热系数；c_i为反应物i的浓度；D_i^eff为气体组分的扩散率，通过稀释近似方法进行计算，

$ D_i^{\mathrm{eff}}=\varepsilon^{1.5}(1-s)^{r_{\mathrm{s}}} D_i^0\left(\frac{p_0}{p}\right)^{\gamma_{\mathrm{p}}}\left(\frac{T}{T_0}\right)^{\gamma_{\mathrm{t}}} . $

(5)

其中：D_i⁰为参考温度和压强(T₀、p₀)下物种i的质量扩散率, s为饱和度，r_s为孔隙堵塞指数。方程中的参数取值为：p₀=101 325 Pa，T₀=300 K，γ_p=1.0，γ_t=1.5，r_s=2.5。

式(2)中，S_u为动量源项，在气体扩散层内可简化为Darcy定律，

$ \varepsilon \boldsymbol{v}=\frac{K_{\mathrm{gdl}}}{\mu} \nabla p. $

(6)

式中K_gdl为GDL的渗透率，m²。

式(3)中，S_h为能量源项，具体如表 1所示。

表 1中：j_sol为固体的面电流密度，j_mem为膜的面电流密度，A/m²；L为水凝结产生的潜热，J/mol；σ_sol和σ_mem分别为固体和膜的电导率，1/(Ω·m)；η_an和η_cat分别为阳极、阴极表面过电位，V；F为Faraday常数，9.65×10⁷ C/kmol；S_ld为液相与固相之间的质量交换速率，S_gl为气液两相之间的质量交换速率，kg/(m³·s)。

式(4)中，S_i为物种源项：

$ S_{\mathrm{H}_2}=-\frac{M_{\mathrm{H}_2}}{2 F} R_{\text {an }}<0, $

(7)

$ S_{\mathrm{O}_2}=-\frac{M_{\mathrm{O}_2}}{4 F} R_{\mathrm{cat}}<0, $

(8)

$ S_{\mathrm{H}_2 \mathrm{O}}=\frac{M_{\mathrm{H}_2 \mathrm{O}}}{2 F} R_{\mathrm{cat}}>0 . $

(9)

式(7)—(9)中：R_an和R_cat分别为阳极、阴极的交换电流密度，A/cm³；M为摩尔质量。

1.4 电化学模型

对于一块质子交换膜燃料电池，其工作原理为氢气与氧气(空气)之间发生氧化还原反应，其中氢气在阳极供应，空气在阴极供应。在阳极和阴极的催化剂层中分别发生以下电化学反应:

$ 阳极: \mathrm{H}_2 \longrightarrow 2 \mathrm{H}^{+}+2 \mathrm{e}^{-},\\ 阴极: \frac{1}{2} \mathrm{O}_2+2 \mathrm{e}^{-}+2 \mathrm{H}^{+} \longrightarrow \mathrm{H}_2 \mathrm{O}. $

这些反应背后的驱动力是表面过电位，即固体的相电位与电解质/膜的相电位之间的差值。因此，在燃料电池和电解模型中需求解两个势能方程：一个方程求解电子在固体间的传递；另一个方程求解H⁺或O^2-的质子输运。两个势能方程分别为:

$ \nabla \cdot\left(\sigma_{\mathrm{sol}} \nabla \phi_{\mathrm{sol}}\right)+R_{\mathrm{sol}}=0, $

(10)

$ \nabla \cdot\left(\sigma_{\text {mem }} \nabla \phi_{\text {mem }}\right)+R_{\text {mem }}=0. $

(11)

式(10)在催化剂层、多孔介质的固体网格和集流器中求解，式(11)在催化剂层和膜中求解。其中：ϕ为电势；膜的电导率为

$ {\sigma _{{\rm{mem}}}} = \beta {(0.514\lambda - 0.326)^\omega }{{\rm{e}}^{1268(1/303 - 1/T)}}. $

(12)

方程(12)中：β和ω为模型常量，用于保证模型的通用性。λ为水含量：

$ \lambda=0.043+17.18 a-39.85 a^2+36 a^3 \quad(a \leqslant 1), $

(13)

$ \lambda=14+1.4(a-1) \quad(a>1). $

(14)

式(13)、(14)中，a为水活度，

$ a = \frac{{{p_{{\rm{wv}}}}}}{{{p_{{\rm{sat}}}}}} + 2s, $

(15)

式(15)中：p_wv和p_sat分别为水蒸气分压和饱和蒸气压，Pa。水蒸气分压可以根据蒸气摩尔分数x_H2O和局部压强计算，而饱和蒸气压只与温度有关。

$ p_{\mathrm{wv}}=x_{\mathrm{H}_2 \mathrm{O}} p, $

(16)

$ \begin{aligned} \lg p_{\mathrm{sat}}=&-2.1794+0.02953(T-273.17)-\\ & 9.1837 \times 10^{-5}(T-273.17)^2+\\ & 1.4454 \times 10^{-7}(T-273.17)^3. \end{aligned} $

(17)

式(10)、(11)中：R_sol和R_mem分别为固体和膜的体电流密度，A/m³。对于固体而言，阳极侧R_sol=-R_an(R_sol < 0)，阴极侧R_sol=R_cat(R_sol>0)。对于膜而言，阳极侧R_mem=R_an(R_mem>0)，阴极侧R_mem= -R_cat((R_mem < 0)。

阳极、阴极的交换电流密度R_an和R_cat的定义为：

$ \begin{gathered} R_{\mathrm{an}}=\left(\zeta_{\mathrm{an}} j_{\mathrm{an}}(T)\right)\left(\frac{\left[\mathrm{H}_2\right]}{\left[\mathrm{H}_2\right]_{\mathrm{ref}}}\right)^{\gamma_{\mathrm{an}}} \cdot \\ \left(\mathrm{e}^{\alpha_{\mathrm{an}}^{\mathrm{cat}} F \eta_{\mathrm{an}} /(R T)}-\mathrm{e}^{-\alpha_{\mathrm{cat}}^{\mathrm{an}} F_{\mathrm{an}} /(R T)}\right), \end{gathered} $

(18)

$ \begin{aligned} &R_{\mathrm{cat}}=\left(\zeta_{\mathrm{cat}} j_{\mathrm{cat}}(T)\right)\left(\frac{\left[\mathrm{O}_2\right]}{\left[\mathrm{O}_2\right]_{\mathrm{ref}}}\right)^{\gamma_{\mathrm{cat}}}\\ &\left(-\mathrm{e}^{\alpha_{\mathrm{an}}^{\mathrm{cat}} F \eta_{\mathrm{cat}} /(R T)}+\mathrm{e}^{-\alpha_{\mathrm{cat}}^{\mathrm{cat}} F \eta_{\mathrm{cat}} /(R T)}\right) . \end{aligned} $

(19)

其中：j(T)为单位活性表面面积的参考交换电流密度，A/m²；ζ为电化学活性比表面积，1/m；[]、[]_ref分别为浓度和参考条件下的浓度，kmol/m³；γ_an为阳极浓度系数，γ_cat为阴极浓度系数；α_an^an和α_an^cat分别为阳极电极的阳极转移系数和阴极转移系数；α_cat^an和α_cat^cat分别为阴极电极的阳极转移系数和阴极转移系数；R为理想气体常数，R=8.314 J/(mol·K)。参考交换电流密度j_an(T)和j_cat(T)都与温度有关：

$ j_{\mathrm{an}}(T)=j_{\mathrm{an}}^{\mathrm{ref}} \mathrm{e}^{-\left(1-T / T_{\mathrm{an}}^{\mathrm{ref}}\right) E_{\mathrm{an}} /(R T)} , $

(20)

$ j_{\mathrm{cat}}(T)=j_{\mathrm{cat}}^{\mathrm{ref}} \mathrm{e}^{-\left(1-T / T_{\mathrm{cat}}^{\mathrm{ref}} / E_{\mathrm{cat}} /(R T)\right.}. $

(21)

其中：j^ref为参考温度下的参考交换电流密度，E为活化能，T^ref为参考温度。

式(18)、(19)中，η_an和η_cat分别为：

$ \eta_{\mathrm{an}} =\phi_{\mathrm{sol}}-\phi_{\mathrm{mem}}-U_{\mathrm{an}}^0 , $

(22)

$ \eta_{\mathrm{cat}} =\phi_{\mathrm{sol}}-\phi_{\mathrm{mem}}-U_{\mathrm{cat}}^0. $

(23)

阴阳极半电池电位U_an⁰和U_cat⁰可由Nernst方程计算：

$ U_{\mathrm{an}}^0=E_{\mathrm{an}}^0-\frac{\Delta S_{\mathrm{an}}}{2 F}\left(T-T_0\right)-\frac{R T}{2 F} \ln \left(\frac{p_{\mathrm{H}_2}}{p_0}\right), $

(24)

$ U_{\mathrm{cat}}^0=E_{\mathrm{cat}}^0+\frac{\Delta S_{\mathrm{cat}}}{2 F}\left(T-T_0\right)-\frac{R T}{2 F} \ln \left(\frac{p_{\mathrm{H}_2 \mathrm{O}}}{p_{\mathrm{sat}} \sqrt{p_{\mathrm{O}_2} / p_0}}\right) . $

(25)

其中：p_H2、p_O2、p_H2O分别代表氢气、氧气、水蒸气的分压；E⁰为可逆电位；ΔS为反应熵，J/(kg·mol·K)。

1.5 液态水输运方程

在PEMFC中，液态水呈现气态、溶解相等状态，存在于催化剂层以及电解质膜中的液态水在阴阳极的产生及输运可描述为

$ \begin{gathered} \frac{\partial}{\partial t}\left(\varepsilon_i M_{\mathrm{w}, \mathrm{H}_2 \mathrm{O}} \frac{\rho_{\mathrm{mem}}}{\mathrm{EW}} \lambda\right)+\nabla \cdot\left(i_{\mathrm{mem}} \frac{n_{\mathrm{d}}}{F} M_{\mathrm{w}}\right)= \\ \nabla \cdot\left(M_{\mathrm{w}} D_{\mathrm{w}} \nabla \lambda\right)+S_\lambda+S_{\mathrm{gd}}+S_{\mathrm{ld}}. \end{gathered} $

(26)

式(26)中: ρ_mem为干膜密度；EW为膜的等效质量, kg/kmol；膜的离子流密度i_mem=－σ_mem∇ϕ_mem；n_d为渗透阻力系数，n_d=2.5 λ/22；D_w为水分扩散系数；催化剂层中阴极副反应导致的产水速率$S_\lambda =\frac{M_{\mathrm{w}, \mathrm{H}_2 \mathrm{O}}}{2 F} R_{\mathrm{cat}} $；S_gd为气相和溶解相之间的质量变化率，kg/(m³·s)。S_gd和S_ld的计算公式分别为：

$ S_{\mathrm{gd}}=\left(1-s^\theta\right) \gamma_{\mathrm{gd}} M_{\mathrm{w}, \mathrm{H}_2 \mathrm{O}} \frac{\rho_{\mathrm{mem}}}{\mathrm{EW}}\left(\lambda_{\mathrm{eq}}-\lambda\right), $

(27)

$ S_{\mathrm{ld}}=s^\theta \gamma_{\mathrm{ld}} M_{\mathrm{w}, \mathrm{H}_2 \mathrm{O}} \frac{\rho_{\mathrm{mem}}}{\mathrm{EW}}\left(\lambda_{\mathrm{eq}}-\lambda\right). $

(28)

其中：θ为液体覆盖指数；γ_gd和γ_ld分别为气相和液相的质量交换速率常数；λ_eq为平衡含水率，

$ \begin{gathered} \lambda_{\mathrm{eq}}=0.3+6 a(1-\tanh (a-0.5))+ \\ 0.69\left(\lambda_{a=1}-3.52\right) a^{0.5}\left[1+\tanh \left(\frac{a-0.89}{0.23}\right)\right]+ \\ s \cdot\left(\lambda_{s=1}-\lambda_{a=1}\right) . \end{gathered} $

(29)

其中：λ_s=1和λ_a=1分别为水饱和度s=1和水活度ɑ=1时的平衡水含量。

多孔电极和膜中液态水输运的驱动力主要是液态水的压力梯度，

$ \frac{\partial}{\partial t}\left(\varepsilon \rho_1 s\right)=\nabla \cdot\left[\frac{\rho_{\mathrm{l}} K K_{\mathrm{r}}}{\mu_{\mathrm{l}}} \nabla\left(p_{\mathrm{c}}+p\right)\right]+S_{\mathrm{gl}}-S_{\mathrm{ld}}. $

(30)

式中：ρ_l为液态水密度，μ_l为液态水动力黏度，K为绝对渗透率，K_r为相对渗透率, p_c为毛细管压强。

在多孔气体扩散层中，相对渗透率K_r=s^b，b为相对渗透率指数；而在膜中，相对渗透率为

$ K_{\mathrm{r}}=\left(\frac{\frac{M_{\mathrm{w}, \mathrm{H}_2 \mathrm{O}}}{\rho_1} \lambda_{s=1}+\frac{\mathrm{EW}}{\rho_i}}{\frac{M_{\mathrm{w}, \mathrm{H}_2 \mathrm{O}}}{\rho_{\mathrm{l}}} \lambda+\frac{\mathrm{EW}}{\rho_i}} \cdot \frac{\lambda}{\lambda_{s=1}}\right)^2 . $

(31)

其中ρ_i为物种i的密度。

气液两相之间的传质速率根据单向扩散理论进行计算，

$ {S_{{\rm{gl}}}} = \left\{ {\begin{array}{*{20}{l}} {{\gamma _{\rm{e}}}\varepsilon s{D_{{\rm{gl}}}}\frac{{{M_{{\rm{w}}, {{\rm{H}}_2}{\rm{O}}}}}}{{RT}}p\;\ln \frac{{p - {p_{{\rm{sat}}}}}}{{p - {p_{{\rm{wv}}}}}}, \;\;\;}&{{\kern 1pt} {p_{{\rm{wv}}}} \le {p_{{\rm{sat}}}}};\\ {{\gamma _{\rm{c}}}\varepsilon (1 - s){D_{{\rm{gl}}}}\frac{{{M_{{\rm{w}}, }}{{\rm{H}}_2}{\rm{O}}}}{{RT}}p\;\ln \frac{{p - {p_{{\rm{sat}}}}}}{{p - {p_{{\rm{wv}}}}}}, }&{{p_{{\rm{wv}}}} > {p_{{\rm{sat}}}}}. \end{array}} \right. $

(32)

其中：γ_e和γ_c分别为蒸发速率系数和冷凝速率系数。D_gl的计算公式为

$ D_{\mathrm{gl}}= \begin{cases}0.365 \times 10^{-4}\left(\frac{T}{343}\right)^{2.334} \times \frac{10^5}{p}, & \text { 阴极; } \\ 1.79 \times 10^{-4}\left(\frac{T}{343}\right)^{2.334} \times \frac{10^5}{p}, & \text { 阳极. }\end{cases} $

(33)

在阴阳极气体通道与气体扩散层之间的接触面处，假设液态水只从气体扩散层流入气体流道，而不允许回流，且其通量由毛细管压力和气体通道中的动压力f_liq驱动，

$ f_{\mathrm{liq}}=\varepsilon s \cdot \max \left[\left(p_{\mathrm{c}}+\frac{1}{2} \rho V^2\right), 0\right]. $

(34)

式(34)中毛细管压强的计算公式如下：

$ p_{\mathrm{c}}=\sigma\left|\cos \theta_{\mathrm{c}}\right| \sqrt{\frac{\varepsilon}{K}} J(s), $

(35)

$ J(s)=a s-b s^2+c s^3. $

(36)

式(35)、(36)中：σ为表面张力，N/m；θ_c为接触角。a、b、c为Leverett函数系数，其默认值分别为：

$ a=1.417, b=2.12, c=1.263. $

气体流道中的液态水通过式(37)进行跟踪，

$ \frac{\partial}{\partial t}\left(\rho_1 s\right)+\nabla \cdot\left(\rho_1 \boldsymbol{v}_1 s\right)=\nabla \cdot\left(D_{\mathrm{liq}} \nabla s\right) $

(37)

式(37)中：D_liq为液态水在气体流道中的扩散系数；v_l为液态水的速度(m/s)，假设其为气体速度v_g的一部分，v_l=χv_g，χ为液气速度比。

1.6 边界条件

各区域边界条件如图 3所示。

1) 表面各边界的边界条件为：电流通量与离子流通量均为0，而固体的电导率σ_sol与膜的离子电导率σ_mem均为常数，因此有$\frac{\partial \phi_{\mathrm{mem}}}{\partial n}=0、\frac{\partial \phi_{\mathrm{sol}}}{\partial n}=0$。其中n表示x、y或z方向。

2) 阳极侧电位为0 V，即ϕ_sol=0 V；阴极侧电位由电池工作电压决定，即ϕ_sol=V_cell。

3) 为尽可能贴近真实实验，在阴阳极入口采用质量流入口边界：

$ m_{\mathrm{an}}=\frac{j \zeta_{\mathrm{an}}}{2 F w_{\mathrm{H}_2}}, \quad m_{\mathrm{cat}}=\frac{j \zeta_{\mathrm{cat}}}{4 F w_{\mathrm{O}_2}}. $

其中：ζ_an为阳极过量系数，ζ_cat为阴极过量系数，w_H2为氢气质量分数，w_O2为氧气质量分数。

4) 流场出口：假设阳极和阴极通道出口边界条件为压力出口，并赋值温度(343 K)和回流压强(1 atm，即101.325 kPa)。

5) 壁面：所有进、出口之间的壁面调整为静止壁面，流速采用无滑移条件。

2 模型验证

由于PEMFC内部涉及多种物理化学现象，如电化学反应、传热传质等，而其结构极为紧凑，且难以实现参数可视化，于是可以实现参数可视化且成本低的数值仿真研究被广泛应用于PEMFC研究中。

2.1 实验方案

实验采用群翌HS 330燃料电池测试台对五蛇流场的PEMFC进行性能测试。电池的活化区域尺寸为50 mm×50 mm，流道宽度和深度均为1 mm，所用膜电极为Nafion^®115。实验条件为：工作温度70 ℃，阳极过量系数1.5，阴极过量系数2，阳极入口气体相对湿度为80%，阴极入口气体相对湿度为80%。

2.2 仿真模型

基于SIMPLE算法，根据本文建立的三维模型使用Fluent进行仿真计算，计算过程中根据质量和能量的残差值判断解是否收敛，收敛的判别标准分别为：质量(1×10^-6)、能量(1×10^-4)。其余电化学参数及结构参数如表 2所示。

2.3 模型有效性验证

为了验证仿真模型的有效性，利用燃料电池测试台对相同工作条件下的五蛇流场PEMFC进行了实验测试。结果如图 4所示。极化曲线和功率密度曲线的计算结果与实验结果在低电流密度下较为吻合，而在高电流密度下误差逐渐增大，且在峰值功率点处误差为2.6%，而平均误差为0.5%，因此可以认为本文建立的仿真模型是准确、有效的。

2.4 网格无关性验证

为尽量节省计算资源，选取功率密度最大的点对网格尺寸与计算精度之间的关系进行研究。通过对面网格进行加密，建立了面网格尺寸分别为0.2、0.25、0.3、0.35、0.4 mm的网格。将网格文件导入Fluent中进行数值模拟，结果如图 5所示。当网格尺寸大于0.35 mm时，计算结果与实验结果之间的误差较大，可达到5.7%；而在网格尺寸减小到0.25 mm时，误差也降至2.606%。考虑到仿生流场结构特殊，本研究采用0.25 mm的网格进行计算。

3 结果与讨论 3.1 极化曲线与功率密度曲线

由于不具备平行流场和仿生流场的实验条件，因此假设各流场仿真结果的误差相同，将各流场放在同样的仿真条件下进行对比。

图 6是采用各流场的电池的极化曲线和功率密度曲线。可以看到，仿生流场相比于传统的平行流场具有更好的性能，但是其性能略低于五蛇流场。在低电流密度下，极化曲线之间的差异较小，但在中高电流密度下，流场之间的性能差异逐渐显现。

从图 6中可以看出，3种流场均在0.52 V时达到性能峰值，且五蛇流场的电流密度最高(达到了1.15 A/cm²)，而平行流场的电流密度最低(仅为0.853 A/cm²)，仿生流场(1.1 A/cm²)仅次于五蛇流场。根据计算，采用仿生流场的燃料电池在峰值功率密度上比平行流场高28.85%，而比五蛇流场低4.36%。

3.2 阴极气体流速分布及压强分布

图 7为采用不同流道的电池的内部流道与气体扩散层交界面处气体流速的分布情况。可以看出, 对于平行流场(图 7a)，气体流速只在气体出入口处最高，这可能是因为内部歧管结构引起的局部损失较高，进而导致在各歧管中气体流速都不高。

对于五蛇流场(图 7b)，其内部气体流速较高。尽管其内部同样有着歧管结构，但气体在入口处便已分配至各“蛇”中，在单个“蛇”中相当于沿着单蛇流道进行运输，沿程损失占主要部分而局部损失很小，因此动量损失较小而流速较高。由图 7b还可看出，每个“蛇”中的气体流速都随流道长度的增加而降低，这是由于沿程损失越来越高以及气体不断被消耗的缘故。

对于仿生流场(图 7c)，其内部气体流速分布情况与平行流场类似：仅在与进出口直接相通的流道中流速较高，而在其余流道中均处于较低水平。这同样是因为内部各种歧管增大了局部损失。

根据Bernoulli方程，气体流速与气体压强通常具有此消彼长的关系，即气体速度低时则气体压强高，而气体速度高时则反之。但是，燃料电池中的反应气体不断被消耗，因此随着气体逐渐远离气体入口，同样会出现压强递减的趋势，这一点在图 8中得到了验证。

如图 8a所示，由于平行流场内部气体流速处于较低水平，故其内部压强处于相对较高水平，但是由于其内部气体过于分散，故而内部压强并没有其余两种流场高。如图 8b所示，五蛇流场内部气体压强沿着每个“蛇”逐渐递减，且在出口附近压强较低，产生了回流，不利于排水。如图 8c所示，仿生流场气体压强整体处于中高水平，且在出口附近仍具有较高压强，这有利于排水。

3.3 阴极毛细管压强分布及流道内液态水分布

燃料电池内部结构紧凑。液态水在燃料电池中受到毛细管压力以及气体压力等力的作用，因此毛细管压力同样影响着燃料电池的排水能力。图 9为3种流场内部的毛细管压强分布。

由图 9a可以看出，平行流场内部毛细管压强处于较低水平，平均值仅有2.389×10^-11 Pa，而五蛇流场平均值为1.131×10^-10 Pa，仿生流场平均值为5.295×10^-11 Pa。可见，3种流场中，五蛇流场的毛细管压强最高，仿生流场次之，平行流场最低。但是，毛细管压力影响的是液态水从气体扩散层中析出的能力，而液态水在流道中的运输主要受气体压力的影响。

图 10为采用各流场的电池阴极流道内液态水的分布。平行流场和五蛇流场中液态水的量相对较多，而仿生流场中仅存在少量液态水。

3.4 氧气质量分数分布及电流密度分布

流道内液态水最终会影响流场中反应物分布及电流密度分布。图 11为各流场内部氧气质量分布云图，图 12为各流场内部电流密度分布云图。

对于平行流场而言，结合图 10a，由图 11a可以看出，平行流场在液态水分布区域的氧气质量分数极低，这可能是因为平行流场内部“水淹”现象较为严重，导致气体无法通过该区域进行扩散。由图 12a可以看出，对应区域的电流密度也处于较低水平，且电流密度整体均匀性较差。

结合图 10b、图 11b以及图 12b来看，对于五蛇流场，虽然其内部存在液态水，但是由于并未影响到气体运输，因此氧气质量分数的分布较为合理，且电流密度分布较为均匀。

结合图 10c、图 11c以及图 12c来看，对于本文设计的仿生流场，尽管其内部仅存在少量液态水，但是由于其内部歧管较多，这些液态水同样影响到了反应物在该区域的扩散。同时，“水淹”区域周边电流密度高于平均水平，消耗了更多的氧气，进一步降低了该区域的氧气质量分数。

综上所述，仿生流场的电流密度均匀性以及反应物分布均匀性均低于五蛇流场，但高于平行流场。

4 结论

为了提高PEMFC的性能及其内部反应物和电流密度的均匀性，本研究根据银杏叶脉形状设计了具有双银杏叶形仿生流场的燃料电池，并利用Fluent对其性能及内部情况与传统的五蛇流场、平行流场进行了仿真对比，主要的研究发现如下：

1) 建立三维稳态多相等温层流模型，利用五蛇流场的实验结果验证了该模型的准确性。设计了银杏叶脉结构的仿生流场，并在相同条件下对该仿生流场、平行流场、五蛇流场进行数值仿真对比，结果表明采用仿生流场的燃料电池峰值功率密度比采用平行流场的电池高28.85%，比采用五蛇流场的电池低4.36%。

2) 分析了采用不同流场的电池内部的气体流速分布、气体压强分布、毛细管压强分布以及流道内液态水分布，结果表明仿生流场比五蛇流场更有利于排出流道内液态水。

3) 分析了采用不同流场的电池内部的氧气质量分数分布，结果表明五蛇流场的内部氧气分布比仿生流场更均匀。

4) 分析了采用不同流场的电池内部电流密度分布，结果表明且五蛇流场内部电流密度较为均匀，不容易产生“热点”，尽管存在小部分区域电流密度较高，但是峰值与谷值之间的差距不大。仿生流场存在小部分区域电流密度较高，且峰值与谷值之间的差距较大，容易产生“热点”，但仍优于平行流场。