Illustrative Case Study on the Performance and Optimization of Proton Exchange Membrane Fuel Cell

Modeling is a powerful tool for the design and development of proton exchange membrane fuel cells (PEMFCs). This study presents a one-dimensional, two-phase mathematical model of PEMFC to investigate the two-phase transport process, gas species transport flow and water crossover fluxes. The model reduces the computational time for PEMFC design with guaranteed accuracy. Analysis results show that the concentration and activation overpotentials of the cell decrease with the increase of operation pressure, which result in enhanced cell performance. Proper oxygen stoichiometry ratio in the cathode decreases the cell activation overpotential and is favorable for performance improvement. The cell ohmic resistance correspondingly increases with the increase of catalyst layer thickness, which leads to a deteriorated cell performance. The improvement on cell performance could be facilitated by decreasing the membrane thickness. Predicted results show that the present model is a useful tool for the design optimization of practical PEMFCs.


Introduction
A fuel cell is an electrochemical energy converter that directly converts the chemical energy of fuel into direct current electricity [1,2].Proton exchange membrane fuel cells (PEMFCs) are considered as ideal clean energy and power sources due to their high-power density and fast start-up at low temperatures [3,4].Remarkable advancements have been achieved in the cost reduction, stack performance and operation life of PEMFCs in the past decades.
On the basis of the studies on heat and mass transfer mechanisms, which ultimately affect cell performance, modeling plays an important and vital part in fuel cell design and cell performance promotion.Weber and Newman [5] established a 1D model to investigate thermal and water management.They obtained the optimal operating temperature and maximum power density in accordance with an external heat transfer coefficient.Caisheng Wang and M.H. Nehrir developed 1D dynamic models of solid oxide fuel cell [6] and PEMFC [7] in MATLAB/Simulink.Their simulation results demonstrated the accuracy of the models under steady-state and transient conditions.Fuller et al. [8] and Jung and Nguyen [9] investigated the thermal management and water hydrobalance by using a 2D model.Dannenberg et al. [10] developed an along-the-channel mass and heat transfer model based on simplified Bulter-Volumer and Stefan-Maxwell equations.Wen et al. [11] built 3D thermal models to investigate the irregular distributions of thermal sources and temperature in a single fuel cell.Berning et al. [12][13][14] established 3D models by using CFX user-defined process and customized iteration routine to obtain the detailed thermal and water management.Liu et al. [15] presented a simplified non-isothermal, single-phase model by using CFD-Micromesh, which contains 200,000 mesh elements.Askarzadeh and Coelho [16] used a backtracking search algorithm combined with Burger's chaotic map to determine the effect of each parameter on the polarization curve.Their method was effective in obtaining the nonlinearity influences caused by operational states in electrochemical processes.For the parameter selection problem of PEMFC modeling, Chakraborty et al. [17] developed a differential-evolution-based solution for PEMFC stack modeling and their results were statistically significant regardless of the algorithmic parameter settings used.
Water is produced on the cathode side as the product of electrochemical reaction.Water flooding occurs, which results in low current densities and poor cell performance when water is inappropriately removed from the cell, especially from the cathode.Many researchers have investigated water management in PEMFCs.Nguyen and Knobbe [18] developed a method of sequentially exhausting each individual cell.Excess water was effectively removed by providing adequate gas.Qin et al. [19] reported their studies on the transport and removal of liquid water droplet at the membrane-electrode assembly (MEA) by using a volume of fluid model.Water removal could be effectively facilitated with the hydrophilic plate based on its in-between surface contact angle of the bottom channel and MEA surfaces.Liu et al. [20] found that the local current density could be promoted due to the local transport enhancement of the reactant gas.The promotion was induced through the blockage effects caused by stack baffles.Jung et al. [21] added hydrophilic SiO 2 particles to the anode catalyst layer and found that the SiO 2 facilitates the hydration and water removal on the anode and cathode, respectively.
The main design philosophy of the channels is by ensuring uniform gas distributions and gas diffusions and facilitating the efficiency of water removal.The common channel designs are straight, crisscross and serpentine.Diffusion coefficient is crucial in the pathway from the flow field channels to the catalyst layer.A series of serpentine channels was compared and evaluated in reference [22] to obtain an appropriate pressure drop.Nguyen [23] presented a new crisscross flow field that avoided water flooding.The reactant gases were forced to flow through the gas diffusion layer (GDL) and entered the electrodes due to the dead-ended inlet and outlet throughout where processed water could be removed.However, this condition increased the parasitic power and reduced the physical strength of the MEA.
Existing studies have focused on 2D and 3D modeling.However, these models are complicated and time consuming in most cases, which are unsuitable to the rapid design and high-efficiency optimization of PEMFCs.In this study, a 1D analytical model is built to obtain the diverse PEMFC performance influenced by structural parameters and operational limits.Thus, an integrated design program is developed, in which material selection, parameter adjustment, structure optimization and cell performance prediction are efficiently and accurately conducted.

Physical Problem Description
Figure 1 shows the schematic of a proton exchange membrane fuel cell.A single PEMFC is divided into an anode electrode and a cathode electrode by a proton-conducting membrane [24].Hydrogen flows through porous electrodes to the anode catalyst layer to react and dissociate into protons and electrons at the anode side.The protons flow through the membrane to the cathode catalyst layer and the electrons released in the anode are transferred through an external circuit to the cathode.The chemical reaction expression at the anode is expressed as follows: At the cathode side, the air flows through the similar porous electrodes to the cathode catalyst layer, in which oxygen reacts with the electrons from the external circuit and the protons passing through the membrane producing water and heat [25].The cathodic reaction expression is expressed as follows: The total reaction formula is given by: A single hydrogen and air PEMFC is considered as the computational domain, which includes bipolar plate, straight flow channels, gas diffusion layers (GDLs), catalyst layers (CLs) and membrane.The major direction of mass transfer is normal to the membrane surface.The work includes major transport that affects the PEMFC behavior and ignores the minor transport, which has negligible effects on cell performance.
A single hydrogen and air PEMFC is considered as the computational domain, which includes bipolar plate, straight flow channels, gas diffusion layers (GDLs), catalyst layers (CLs) and membrane.The major direction of mass transfer is normal to the membrane surface.The work includes major transport that affects the PEMFC behavior and ignores the minor transport, which has negligible effects on cell performance.A fuel cell model includes macro-scale and micro-scale transport processes.Construction of a precise fuel cell model is a remarkable challenge due to the lack of experimental data.Hence, several assumptions are made in establishing a convenient model.The mass and heat transport in an isothermal model are considered to be 1D processes.It is assumed that the flow is steady state and laminar due to the small flow channel size and low flow velocities.All gases are considered to be ideal gas.The contact resistance between different layers is ignored.The porous media, such as GDLs, CLs and membranes, are considered to be isotropic and homogeneous.

Gas Transport
1D gas transport along the proton transfer direction is investigated.Similar to the assumptions in many literature [26,27], only diffusive transport occurs in the proton-transfer direction within the GDLs, CLs and membrane.Various gas reactants and water vapor concentrations at the CL could be obtained by considering 1D mass transport normal to the membrane surface.1D mass transport includes major transport processes, which improves the calculation efficiency.
Humidified hydrogen is supplied to the anode and air is fed to the cathode.The gas species concentrations at the CL are obtained by considering mass transfer along the proton transfer direction.The species conservation equations are listed in Table 1.A fuel cell model includes macro-scale and micro-scale transport processes.Construction of a precise fuel cell model is a remarkable challenge due to the lack of experimental data.Hence, several assumptions are made in establishing a convenient model.The mass and heat transport in an isothermal model are considered to be 1D processes.It is assumed that the flow is steady state and laminar due to the small flow channel size and low flow velocities.All gases are considered to be ideal gas.The contact resistance between different layers is ignored.The porous media, such as GDLs, CLs and membranes, are considered to be isotropic and homogeneous.

Gas Transport
1D gas transport along the proton transfer direction is investigated.Similar to the assumptions in many literature [26,27], only diffusive transport occurs in the proton-transfer direction within the GDLs, CLs and membrane.Various gas reactants and water vapor concentrations at the CL could be obtained by considering 1D mass transport normal to the membrane surface.1D mass transport includes major transport processes, which improves the calculation efficiency.
Humidified hydrogen is supplied to the anode and air is fed to the cathode.The gas species concentrations at the CL are obtained by considering mass transfer along the proton transfer direction.The species conservation equations are listed in Table 1.

Species Species Conservation Equation Region
Hydrogen (c Flow channel

CLs
As shown in Table 1, c i k denotes gas molar concentration at the interface and A ch−GDL is the contact area between the flow channel and GDL.The species concentration equation is corrected by Sherwood number, Sh due to the negligible convection effect.Sh is the ratio of convective mass transfer to mass diffusion.It is equal to 2.3 in the laminar flow channel in much of the literature because the diffusion transport is dominant in the proton transfer direction [28][29][30].D eff i,k is the effective diffusivity of species i in electrode k (GDL, CL).On the basis of Bruggeman correlation, the expression of gas effective diffusivity can be expressed as [31]: where D i is the intrinsic diffusivity [32]: 15 ) 1.5 15 ) 1.5 D vap = 2.982 × 10 −5 ( T 333.15 ) 1.5 Assuming that the gas concentration between the GDL and CL is linear.

Liquid Water Transport
Water is produced at the cathode side due to electrochemical reaction.The water used for anode humidification is transported through the electrolyte to the cathode due to electroosmotic drag phenomenon.The model considers the transport of liquid water in the porous media.The liquid water conservation equations are presented in Table 2.

Species Liquid Water Conservation Equation Region
Liquid water Anode side of the membrane As shown in Table 2, M H 2 O is the molar mass of water molecule and is equal to 0.018 kg/mol.c lq,j i refers to the liquid water concentrations, in which i refers to the interface between adjacent domains and j is the anode or cathode.I cat,lq is the current density by producing liquid water.D lq refers to the diffusivity of liquid water [33]: where K 0 is the intrinsic permeability.K lq is the permeability of liquid transport in the electrode considering the blockage of liquid water.Capillary pressure p c affects liquid water transport and is related to liquid water volume fraction (s lq ), which can be expressed as [33]: (44.9s lq −14.41) + 2.09e (−22.2slq +7.13) in CL −2395 − 2431e (92.36s lq −52.37) + 2431e (−0.0088s lq +0.005) in GDL (10) n d denotes the electro-osmotic drag coefficient and is related with water content λ in the membrane [34]: S vap−lq,CL and S vap−lq,GDL represent the source terms, which are due to the evaporation or condensation between liquid water and vapor and they are calculated by:

Cell Performance
Effective work from the fuel cell could be obtained only when current is generated [35].In fact, the actual output voltage is remarkably lower than that of the theoretical voltage due to irreversible losses.These irreversible losses (also called overpotentials) include activation, ohmic and concentration losses [36].The output voltage of single cell due to irreversible losses is determined by: With the reference pressure and temperature values, the Nernst equation can be modified by [37]: where p H 2 refers to the partial pressure of hydrogen.The gas partial pressures changes with current density [17,38] and p H 2 could be evaluated by: ) − RH ano p sat (15) Considering that the model is assumed to be isobaric, p ano and p cat are equal to operating pressure P opt .The formula of the partial pressure of oxygen p O 2 could be calculated as: Activation overpotential is due to the sluggish kinetics of electrochemical reaction [39].This loss dominates at low current density.Among the two activation overpotentials, anodic activation loss could be ignored compared with cathodic activation loss because the oxygen reduction reaction requires a high overpotential [28].The Butler-Volmer equation is used to describe the relationship between current density and activation overpotential.Considering that the solution of Butler-Volmer equation is difficult, it could be simplified to Tafel equation.Thus, the activation overpotential could be given by [40]: Ohmic overpotential is caused by cell electrical resistance.Resistance losses are due to the ionic resistance in the electrolyte and electronic resistance flowing through the electrodes.On the basis of Ohm's law, the voltage loss arising from the resistance of charge transport could be expressed as: The contact resistance between different layers is ignored in which the total cell resistance is equivalent to the sum of the internal resistances of various regions.R MEM , R CL , R GDL , R ch and R BP are the potential drops.The potential drop is attributed to the ohmic resistance of different regions.The internal resistance of each region could be calculated by: where κ i refers to the conductivity of each region i (membrane, CL, GDL, flow channel and bipolar plate).The ion electrical conductivity of the membrane could be calculated by using an empirical correlation [29].
where λ is the water content in the membrane and is evaluated according to [41]: where water activity a could be defined as [27]: where p sat is the saturation pressure of water at the operating temperature and is expressed as follows [42]: Concentration overpotential is attributed to mass transport limitation.Limiting current density (I D ) is produced [31] when the concentrations of reactants at the catalyst reaction sites are reduced to zero in the limiting case.The causes of reactant depletion increase the reaction rates and reduce the reactant transport rates.The concentration overpotential is expressed as the following equation [34]: Table 3 lists the value for the parameters used in the present model.These parameters are related to the materials and the electrochemical properties of the PEMFC.

Boundary Conditions
A constant stoichiometry ratio mode is applied for the two reactants (hydrogen and oxygen).The concentration of gases i (H 2 , O 2 , vapor) at the inlet is calculated as: where c i in refers to the species molar concentration at the inlet.P i represents the inlet partial pressure of gas i.Considering that the reactants are consumed in the electrochemical reaction, the species molar concentration at the outlet can be calculated as follows [29]: where A ch-GDL is the interface area between the flow channel and GDL.ε GDL is the porosity of GDL.
A in is the inlet area of the flow channel.V i in is the gas flow velocity: where ST represents the stoichiometry ratio, which is defined as the ratio of the amount of reactant supplied to the amount of reaction to generate the current density.A rea refers to the reaction area of CL.Species concentration distribution is assumed to be linear in the flow channel and the average value of species molar concentration can be expressed by: The concentration of liquid water is considered to be constant in the flow channel because the straight flow channel of PEMFC is short.

Results and Discussion
The numerical calculation of the model is implemented in MATLAB.The comparison among the predicted polarization curve, model predictions in Ref. [43] and published experimental results in Ref. [44] is presented.In this section, the effects of various operating (operating pressure, oxygen stoichiometry ratio) and design parameters (CL thickness, membrane thickness) on the transport processes and cell performance are analyzed.The input physical parameters and operating conditions are shown in Table 4.The parameters in Table 4 are considered as the default values of a validation model to verify the accuracy of the proposed model.For parametric studies, one of the parameters is changed and the others are kept as their default values.Several numerical parameters are introduced in the proposed model for calculation.These parameters, such as the effective diffusivity of gas species, capillary diffusivity of liquid water and ion electrical conductivity of the membrane, can be calculated by using empirical correlations (Equations ( 5)-( 7), ( 8)-( 10) and ( 22)), which are taken from references [25,31,32].The cell performance can be obtained by solving the modeling equations, including the species conservation, liquid water conservation and electrochemical kinetic equations.The polarization curves of the proposed model are plotted together with other experimental data [44] and numerical model results [43] at a cell temperature of 80 • C and gas pressure of 1.0 atm.The fuel cell experiments completed by Kim K.H. et al. in Ref. [40] is a single cell with fabricated MEA of 0.4 mg/cm 2 platinum loading for the anode and cathode layers and 25 cm 2 active electrode area.The cathodic and anodic gas flow rates are controlled to obtain stoichiometric ratios of 1.5 and 2.0, respectively.The relative humidity of cathode gases and hydrogen are fixed at 100%.The single cell is operated at 80 • C under 1.0 atm pressure conditions.At the same conditions, the polarization characteristics of the model and experiment are compared.As shown in Figure 2, the present model predictions, other numerical results and experimental data are in good agreement with a maximum relative deviation of less than 5%.These results verify the accuracy of the 1D model.
Several numerical parameters are introduced in the proposed model for calculation.These parameters, such as the effective diffusivity of gas species, capillary diffusivity of liquid water and ion electrical conductivity of the membrane, can be calculated by using empirical correlations (Equations ( 5)-( 7), ( 8)-( 10) and ( 22)), which are taken from references [25,31,32].The cell performance can be obtained by solving the modeling equations, including the species conservation, liquid water conservation and electrochemical kinetic equations.The polarization curves of the proposed model are plotted together with other experimental data [44] and numerical model results [43] at a cell temperature of 80 °C and gas pressure of 1.0 atm.The fuel cell experiments completed by Kim K.H. et al. in Ref. [40] is a single cell with fabricated MEA of 0.4 mg/cm 2 platinum loading for the anode and cathode layers and 25 cm 2 active electrode area.The cathodic and anodic gas flow rates are controlled to obtain stoichiometric ratios of 1.5 and 2.0, respectively.The relative humidity of cathode gases and hydrogen are fixed at 100%.The single cell is operated at 80 °C under 1.0 atm pressure conditions.At the same conditions, the polarization characteristics of the model and experiment are compared.As shown in Figure 2, the present model predictions, other numerical results and experimental data are in good agreement with a maximum relative deviation of less than 5%.These results verify the accuracy of the 1D model.Polarization curves of present prediction compared with model prediction [43] and experimental data [44].

Operating Pressure
The influence of operating pressure on the performance of PEMFC is evaluated.Figure 3 shows the cell performance with respect to various operating pressures.The gas inlet pressures at the anode and cathode are kept the same.As shown in Figure 3a, a high output voltage can be obtained at high pressure.This condition is because the partial pressure of the reactants increases with the increment of operating pressure, as indicated in the Nernst equation.
The enhancement of cell performance could be explained by investigating the response of voltage loss.Figure 3b presents three voltage losses changing with current density under three different operating pressures.The results suggest that the better performance at a high pressure is mainly due to the drop in activation overpotential.The reduction percent in activation overpotential approximately attains 25% when the pressure increases from 1 atm to 3 atm and current density is 1.0 Although the formula of the activation overpotential does not show any pressure term, the variation of reactant concentration with pressure affects the change of activation overpotential.
As shown in Figure 3c,d, a high operating pressure indicates high reactant diffusion rates.Thus, the reactants could be rapidly replenished and the concentration of reactants is increased at the reactive site to facilitate the electrochemical processes.Furthermore, the concentration and activation overpotentials of the cell are decreased at a high pressure, which result in an enhanced cell performance.Compared with activation overpotential, the variation of concentration overpotential is negligible.[43] and experimental data [44].

Operating Pressure
The influence of operating pressure on the performance of PEMFC is evaluated.Figure 3 shows the cell performance with respect to various operating pressures.The gas inlet pressures at the anode and cathode are kept the same.As shown in Figure 3a, a high output voltage can be obtained at high pressure.This condition is because the partial pressure of the reactants increases with the increment of operating pressure, as indicated in the Nernst equation.
The enhancement of cell performance could be explained by investigating the response of voltage loss.Figure 3b presents three voltage losses changing with current density under three different operating pressures.The results suggest that the better performance at a high pressure is mainly due to the drop in activation overpotential.The reduction percent in activation overpotential approximately attains 25% when the pressure increases from 1 atm to 3 atm and current density is 1.0 A/cm −2 .Although the formula of the activation overpotential does not show any pressure term, the variation of reactant concentration with pressure affects the change of activation overpotential.
As shown in Figure 3c,d, a high operating pressure indicates high reactant diffusion rates.Thus, the reactants could be rapidly replenished and the concentration of reactants is increased at the reactive site to facilitate the electrochemical processes.Furthermore, the concentration and activation overpotentials of the cell are decreased at a high pressure, which result in an enhanced cell performance.Compared with activation overpotential, the variation of concentration overpotential is negligible.Considering that the expression derived for the ohmic overpotential has no relation with operating pressure, no change is observed in the ohmic overpotential in terms of operating pressure.
In general, the cell performance could be improved by increasing the operating pressure.However, a high pressure costs additional energy from the air compressor, which ultimately offsets the voltage increments.
Considering that the expression derived for the ohmic overpotential has no relation with operating pressure, no change is observed in the ohmic overpotential in terms of operating pressure.
In general, the cell performance could be improved by increasing the operating pressure.However, a high pressure costs additional energy from the air compressor, which ultimately offsets the voltage increments.

Oxygen Stoichiometry Ratio
Polarization curves and ohmic resistance at different current densities under various oxygen stoichiometries are shown in Figure 4. Slight oxygen stoichiometry ratio has a slight influence on cell performance, as shown in Figure 4a.This finding is due to the process simplification adopted in the present model.In general, the cell potential improves with the increase of oxygen stoichiometry due to the augmented oxygen concentration in the CL reaction sites.
A high oxygen stoichiometry ratio represents a high air flow rate, which rapidly replenishes the oxygen consumed by electrochemical reactions.Thus, the activation overpotential of the cell could be suppressed.Furthermore, the high flow rate of reactants leads to an effective water removal from the gas transfer channel, which prevents flooding by blocking the gas transfer channel.Then, a low water flooding in the GDL leads to a small gas transfer resistance to the CL.Therefore, the activation and concentration overpotentials could be reduced by increasing oxygen stoichiometry for a specific current density, as presented in Figure 4b.A large amount of energy should be consumed to compress the air to achieve a high oxygen stoichiometry.

Oxygen Stoichiometry Ratio
Polarization curves and ohmic resistance at different current densities under various oxygen stoichiometries are shown in Figure 4. Slight oxygen stoichiometry ratio has a slight influence on cell performance, as shown in Figure 4a.This finding is due to the process simplification adopted in the present model.In general, the cell potential improves with the increase of oxygen stoichiometry due to the augmented oxygen concentration in the CL reaction sites.

Thickness of CL and Membrane
The influence of various CL thickness on the PEMFC polarization curve from the present model is investigated and displayed in Figure 5. Figure 5a shows that a thick CL causes the improvement on the cell performance at a low current density.The activation overpotential of the cell remarkably drops with the increase of the effective reaction area and CL thickness.The increment of ohmic overpotential is less than that of the decrease of activation overpotential in which the cell performance could be enhanced.A high oxygen stoichiometry ratio represents a high air flow rate, which rapidly replenishes the oxygen consumed by electrochemical reactions.Thus, the activation overpotential of the cell could be suppressed.Furthermore, the high flow rate of reactants leads to an effective water removal from the gas transfer channel, which prevents flooding by blocking the gas transfer channel.Then, a low water flooding in the GDL leads to a small gas transfer resistance to the CL.Therefore, the activation and concentration overpotentials could be reduced by increasing oxygen stoichiometry for a specific current density, as presented in Figure 4b.A large amount of energy should be consumed to compress the air to achieve a high oxygen stoichiometry.

Thickness of CL and Membrane
The influence of various CL thickness on the PEMFC polarization curve from the present model is investigated and displayed in Figure 5. Figure 5a shows that a thick CL causes the improvement on the cell performance at a low current density.The activation overpotential of the cell remarkably drops with the increase of the effective reaction area and CL thickness.The increment of ohmic overpotential is less than that of the decrease of activation overpotential in which the cell performance could be enhanced.

Thickness of CL and Membrane
The influence of various CL thickness on the PEMFC polarization curve from the present model is investigated and displayed in Figure 5. Figure 5a shows that a thick CL causes the improvement on the cell performance at a low current density.The activation overpotential of the cell remarkably drops with the increase of the effective reaction area and CL thickness.The increment of ohmic overpotential is less than that of the decrease of activation overpotential in which the cell performance could be enhanced.
As shown in Figure 5b, the cell output potential decreases with an increment in the CL thickness due to the limited mass diffusion rates in the CL.The results are consistent with those published in Ref. [28].The variations of cell potential are mainly attributed to the variations of ohmic overpotential.The increment of CL thickness causes a long distance for the electrons and ions to arrive at the CL and the transport distance of reactants becomes long.The prolonged transport distance for reactants, electrons and ions leads to a large ohmic resistance and a low cell performance.Although the effective reaction area increases with the CL thickness, the PEMFC activation overpotential is actually not affected by the variation of CL thickness.The performance of the cell could be obviously improved with a thin CL when it operates at a high current density.Figure 6 presents the variation of the cell output and voltage losses at three different membrane thickness.The improvement on cell performance is obtained by decreasing the membrane thickness, as shown in Figure 6a. Figure 6b indicates that the ohmic overpotential has a remarkable influence on cell potential.A thin membrane leads to a small membrane resistance for a specific current density and results in a low ohmic overpotential because the resistance for proton transport to the reaction sites is reduced.Furthermore, the increasing membrane thickness has less remarkable impact on the ohmic As shown in Figure 5b, the cell output potential decreases with an increment in the CL thickness due to the limited mass diffusion rates in the CL.The results are consistent with those published in Ref. [28].The variations of cell potential are mainly attributed to the variations of ohmic overpotential.The increment of CL thickness causes a long distance for the electrons and ions to arrive at the CL and the transport distance of reactants becomes long.The prolonged transport distance for reactants, electrons and ions leads to a large ohmic resistance and a low cell performance.Although the effective reaction area increases with the CL thickness, the PEMFC activation overpotential is actually not affected by the variation of CL thickness.The performance of the cell could be obviously improved with a thin CL when it operates at a high current density.
Figure 6 presents the variation of the cell output and voltage losses at three different membrane thickness.The improvement on cell performance is obtained by decreasing the membrane thickness, as shown in Figure 6a. Figure 6b indicates that the ohmic overpotential has a remarkable influence on cell potential.A thin membrane leads to a small membrane resistance for a specific current density and results in a low ohmic overpotential because the resistance for proton transport to the reaction sites is reduced.Furthermore, the increasing membrane thickness has less remarkable impact on the ohmic overpotential compared with the increasing CL thickness.This condition is because the CL thickness has dominant effects on oxygen transport from the flow channel to the cathode CL and the membrane thickness has no relation with oxygen transport.Hence, fuel cells with optimum CL and membrane thickness are designed to reduce the cost and improve the cell performance.overpotential compared with the increasing CL thickness.This condition is because the CL thickness has dominant effects on oxygen transport from the flow channel to the cathode CL and the membrane thickness has no relation with oxygen transport.Hence, fuel cells with optimum CL and membrane thickness are designed to reduce the cost and improve the cell performance.

Conclusion
A 1D two-phase mathematical model of PEMFC is presented in this study.The model reproduces the effects of the operation conditions on PEMFC performance.The results of the proposed model have a good agreement with model prediction and experimental data.The model could reduce computational time for the design of PEMFCs with guaranteed accuracy.The results indicate that a high operating pressure results in a high diffusion rates of oxygen through the GDL to the CL and the ohmic and activation losses are reduced to enhance the cell performance.The cell performance is continuously improved with the increase of oxygen stoichiometry ratio from 1 to 3 because considerable oxygen are diffused in the CL for electrochemical reaction, which result in a decreased activation loss.Decreasing the catalyst and membrane thickness generally improves the cell performance.The transport distance of ion, electron and reactants becomes long with the increase of catalyst and membrane thickness, which leads to a large ohmic resistance associated with deteriorated cell performance.

Figure 1 .
Figure 1.Schematic of a proton exchange membrane fuel cell.

Figure 1 .
Figure 1.Schematic of a proton exchange membrane fuel cell.

Table 2 .
Liquid water conservation equations.