Assessment of Sensitivity to Evaluate the Impact of Operating Parameters on Stability and Performance in Proton Exchange Membrane Fuel Cells

As a highly nonlinear system, the performance of proton exchange membrane fuel cell (PEMFC) is controlled by various parameters. If the effects of all parameters are considered during the performance optimization, low working efficiency and waste of resources will be caused. The development of sensitivity analysis for parameters can not only exclude the parameters which have slight effects on the system, but also provide the reasonable setting ranges of boundary values for simulation of performance optimization. Therefore, sensitivity analysis of parameters is considered as one of the methods to optimize the fuel cell performance. According to the actual operating conditions of PEMFC, the fluctuation ranges of seven sets of parameters affecting the output performance of PEMFC are determined, namely cell operating temperature, anode/cathode temperature, anode/cathode pressure, and anode/cathode mass flow rate. Then, the control variable method is used to qualitatively analyze the sensitivity of main parameters and combines with the Monte Carlo method to obtain the sensitivity indexes of the insensitive parameters under the specified current density. The results indicate that among these parameters, the working temperature of the fuel cell is the most sensitive to the output performance under all working conditions, whereas the inlet temperature is the least sensitive within the range of deviation. Moreover, the cloud maps of water content distribution under the fluctuation of three more sensitive parameters are compared; the results verify the simulated data and further reveal the reasons for performance changes. The workload of PEMFC performance optimization will be reduced based on the obtained results.


Introduction
Proton exchange membrane fuel cell (PEMFC) is a kind of device that converts chemical energy in hydrogen directly into electrical energy. Its working principle can be concluded as the oxidation-reduction reaction of hydrogen and oxygen to generate electrons, and the transmission of electrons to form current [1,2]. The rate and degree of this process determine the performance of PEMFC, while the process is controlled by the coupling of multiple parameters [3][4][5].
The research on the influencing factors of PEMFC performance has always been a hot spot [6][7][8]. Abdallah et al. [9] obtained the optimal design parameters of PEMFC, which enables the best performance by comparing the power density curves of 30 kinds of cross-sectional flow channels under the same conditions. Zhao et al. [10,11] focused on the structural parameters of pore structures in PEMFCs and how to impact the mass transport and performance of PEMFCs. In addition to the inherent design parameters of PEMFC, the

•
To make the obtained sensitivity indexes comparable, this paper combined with the actual operating conditions to uniformly delineate the calibration parameters and their fluctuation ranges before starting the sensitivity analysis. Moreover, a reliable parameter sensitivity assessment method was designed. • The Monte Carlo method was adopted to quantitatively analyze the parameters, which have little effect on the output performance (failed to judge by the deviation percentages) under calibration conditions. The parameter sensitivity rankings under specific current density were obtained. • In addition to determining the impact of single parameter fluctuations on the output voltage, the water content distribution cloud maps of the contact surface between anode side catalyst layer (CL) and membrane are used to prove the improvement/deterioration of PEMFC performance, which can intuitively explain the phenomena.
Finally, according to the sensitivity indexes, the parameters were classified as insensitive, sensitive, and highly sensitive. The results provide a reference for the optimization of PEMFC and the controls of operating parameters.

Governing Equation
In [34], the equations of conservation of mass, momentum, energy, species, and charges were summarized into a general convection-diffusion equation, as follows: However, this equation is applicable in the steady state. In this study, an unsteady state model was constructed. Thus, a new general governing equation may be formulated, as follows: where ρ and p refer to nominal density and pressure, respectively, φ is a general variable and the velocity vector is represented by ν. The symbols p, φ, ρ, and Γ φ for various equations are listed in Table 1. Table 1. ρ, φ, Γ φ of different governing equations.

Electrochemical Reaction Model
After describing the phenomena of mass transfer, heat transfer, and electricity transfer in PEMFC using the governing conservation equation, the electrochemical reaction model of the cathode catalyst layer was constructed. According to Dickinson et al. [35], most PEMFC models use the Butler-Volmer equation to describe the electrode kinetics: where α Rd and α Ox represent the transfer coefficients of the reduction and oxidation reactions, respectively. F, T and R represent the Faraday, temperature, and gas constants, respectively. E denotes the electrode potential and E r refers to the equilibrium potential; the difference between them is referred to as the overpotential, which is usually represented by η, see Equations (4) and (5).
The overpotential of the anode is positive, which refers to electrons moving away from the electrode. Similarly, the value obtained from Equation (5) is negative. For the anode and cathode reactions to occur in the fuel cell, the Butler-Volmer equation is also effective: where j(T) is the reference exchange current density per effective surface area (A/m 2 ), ξ is the activation area ratio (1/m), [C] ref represent the ratio of the local component concentration to the reference value in the anode and cathode, respectively, γ refers to the concentration index, α an Rd and α cat Rd are the anode and cathode charge transfer coefficients of the reduction reaction, respectively, α an Ox and α cat Ox express the anode and cathode charge transfer coefficients of the oxidation reaction.

Membrane Electrode Water Transport Model
The following transport equation describes the membrane electrode liquid water saturation: where S is the liquid water saturation, which characterizes the ratio of the volume occupied by the liquid water to the volume of the void in the porous medium of the membrane electrode. λ l , M l , i e , ρ l , and ρ g , refer to the relative mobility of liquid water, the molar mass of liquid water, electrolyte phase current density, liquid water density, and gaseous water density, respectively. D c is the capillary diffusion coefficient of liquid water and is defined by the following Equation: In this equation, k l , p c , µ l , and S l represent the liquid water permeability, capillary pressure, dynamic viscosity of liquid water, and source term of water phase transition, respectively; S l is defined as where X w represents the mass fraction of water vapor on the proton membrane surface and k cond and k evap are the condensation of water and evaporation rate coefficients, respectively. The latent heat of phase transition in water is expressed as where h c is the latent heat of phase transition in water. The water content in the membrane is maintained in a dynamic equilibrium state, which is essential for managing the water content in the fuel cell under working conditions. When the cell is working, the water content in the membrane is determined using water transport by the membrane. Water can be transported in four ways in a proton membrane: through ionic conduction, electroosmotic drag, back-diffusion, and hydraulic permeability.
As one of the most significant parameters within the PEM, the proton conductivity is always represented by the symbol σ (S/cm), which is used to evaluate the transport capacity of the protons. The proton conductivity of the PEM materials with resistance can be expressed as where R S refers to the overall impedance of the membrane samples, I is the distance between the reference electrodes, and S is the cross-sectional area of the membrane samples. Electroosmotic drag is caused by the anode side of the polar water molecules and proton hydration, and (H + (nH 2 O)) is formed. The water molecules migrate from the anode to the cathode side with protons, and the corresponding water migration flux caused by the electroosmotic dragging of the proton membrane can be expressed as where i and λ are the current density and water content of the membrane, respectively. n d (λ) is the electroosmotic drag coefficient, which represents the number of water molecules carried by each proton (dragged) transported. n d (λ) is a function of the proton membrane water content (λ) and can be expressed by the following widely used empirical expression: The reverse diffusion water transfer flux is related to the water concentration gradient on both sides of the proton membrane [36,37]: where D m represents the diffusion coefficient of water molecules in the membrane and ∂c ∂y is the water concentration gradient, which is distributed on both sides of the proton membrane.
As one of the water transport methods within the membrane, hydraulic permeability is caused by the pressure difference between the anode and cathode. Compared to other transport mechanisms, the amount of water transferred through hydraulic permeability is at least one order of magnitude lower [38]. Therefore, the impact of this parameter is ignored in this paper. The source terms of various governing equations are listed in Table 2.

Boundary Conditions and Material Properties
The Dirichlet boundary condition was applied in this simulation.
(1) Anode/cathode adopted mass flow inlet (AMFR/CMFR): C a = C a,in was used to define the mass fraction of each component of the reaction gas at the inlet. (2) Channel outlet employed pressure outlet: P = P an,out , P = P ca,out was used to define the anode and cathode pressures, respectively. Moreover, anode operating pressure (AOP) and cathode operating pressure (COP) are usually controlled by back pressure of PEMFC. (3) The voltage or current density of the fuel cell model was set at the bipolar plate surface; the potential on the surface of the end plate of the anode and cathode bipolar plates can be expressed as φ an The no-slip wall boundary was applied to the wall interface between the fluid and solid computing domains. (5) The operating temperature (COT) is expressed by the wall temperature of PEMFC [9,39], which can be independent of the anode inlet temperature (AIT) and cathode inlet temperature (CIT).

Properties of reactants:
(1) Stoichiometric ratio of hydrogen to air (excess coefficient): In a fuel cell that uses high-purity hydrogen (99.999%) as fuel and air as oxidant, the excess coefficient of the anode can be expressed as m H 2 (16) where m in H 2 is the actual flow supplied by the anode hydrogen and . m H 2 represents the theoretical demand flow when the target refers to the current density.
Similarly, the cathode excess coefficient can be expressed as (2) Relative humidity: To decrease the ohmic resistance, humidifying the reactant gas is necessary. The relative humidity (RH) is usually used to characterize the degree of inlet humidity and can be defined as where P in sat refers to the saturated vapor pressure corresponding to the humidification temperature of the gas inlet and P sat (T) is the saturated vapor pressure corresponding to the operating temperature of cell.

Geometric Model and Meshing
At least nine computational fields exist in a three-dimensional model of a complete PEMFC-bipolar plates sites in the cathode and anode, respectively, cathode/anode flow channels, cathode/anode gas diffusive layers (GDL), cathode/anode catalyst layers (CL), and the proton exchange membrane (PEM). Owing to the feature of parallel channel symmetry, to reduce the amount of calculation, half of a complete PEMFC model with a single straight flow channel was truncated along the symmetric plane. A 3D model design software known as UG (Unigraphics NX, developed by Siemens) was used in this study, and the computational domain model covered all directions. In this study, the grid division strategy of a regular hexahedral structure was adopted; compared to unstructured grids such as tetrahedral prisms, it calculates quickly, converges easily, and controls the grid nodes conveniently.
Meshing in ICEM is depicted in Figure 1a. Four sets of regular hexahedral structural grids were prepared with the number of grid elements 100,468, 128,600, 140,400, 158,400 to validate the grid independence. The current density under the given working voltage can comprehensively reflect the degree of electrochemical reaction inside the cell, the multicomponent transportation of the reaction gas, and the performance of charge travel. Thus, this study validated the grid independence with the current density corresponding to the operating voltage of 0.7 V under certain working conditions, which set the intake and exhaust pressures of the cathode/anode as 2 atm, operating temperature as 80 • C, and humidity as 100%, as depicted in Figure 1b. The current density corresponding to the number of grid elements, 100,468, was approximately 0.386 A/cm 2 ; the number of grid elements should be avoided because of its larger deviation compared to other grids. Comprehensively considering the computational stability and efficiency, the number of grid elements of 140,400 was eventually elected as the calculated number of grid elements. The details of the configuration are presented in Table 3.

Model Validation
A variety of phenomena that can be expressed by the polarization curve, such as convection and diffusion transport performance of reactive gas in PEMFC, charge transfer resistance in the membrane and solid phases, as well as the state of electrochemical reaction, have a comprehensive impact on cell performance. Therefore, it is usual to contrast the tendency of the polarization curve in simulation with that in the experiment.
In this paper, the working parameter was set the same as that in [9]; thereafter, the numerical simulation was started. The simulated results of this paper are compared with the simulated results obtained by Mohammedi et al., as can be seen in Figure 2. The numerical results are basically consistent with the polarization curve and power density curve in literature; compared with the middle-current-density region, the numerical agreement is better in the high-and low-current-density region. There are the larger differences in the ohmic region, which means more serious ohmic polarization loss in the simulated model of this paper. The main reason is that some physical properties of PEMFC are not provided in the reference, especially the electronic conductivity, ionic conductivity, and anode/cathode reference concentration. According to the literature [40,41], the ionic conductivity has remarkable effects on the resistance of PEMFC. Moreover, the reference concentration of anode and cathode and the permeability coefficient and viscosity of the mixed gas all have effects on the output performance of PEMFC. However, the error is within an acceptable range, which can be used to predict the performance.

Parameter Setting
The Monte Carlo method with multiparameter sensitivity analysis was used in this part of the experiment, under the common calibration condition of vehicle-used PEMFC, which is aimed at quantitatively analyzing the impact of fluctuating parameters (±5%) on the sensitivity of output performance.
The common calibration conditions used in cars include the idle condition and rated condition. The idle condition means starting the fuel cell engine but not driving the car. At this time, the energy provided by the fuel cell system is required only to maintain the normal working of the subsidiary parts.
In this part, two types of working conditions are used as inlet conditions in the single cell. Under a state of low power, the output current density of the cell is approximately 0.1 A/cm 2 . The rated condition means that the cell is running at the rated output power and the output current density is approximately 0.6 A/cm 2 . A current density of 0.3 A/cm 2 is chosen as the intermediate working condition. Table 4 lists the calibration parameters and their range of fluctuation. Table 5 provides the inlet boundary conditions.  The impact of fluctuation of the operating parameters (the corresponding current densities are 0.1 A/cm 2 , 0.3 A/cm 2 , and 0.6 A/cm 2 ) within a certain range on the cell output voltage under common operating conditions of automobiles is studied in this paper and the results can be seen in Section 3.3. The calibration parameters and their fluctuation ranges are displayed in Table 4. Under these conditions, the behavior of ∆V, which represents the deviation of the PEMFC output voltage as well as the percentage change in voltage deviation, was investigated. The percentage change in the voltage deviation is defined as |(∆V)/V| × 100%, where V is the calibration condition voltage and ∆V = V test − V.

Simulation Results
The pictures of polarization curves and voltage deviations are presented in different parts of Section 3.3 to analyze the effect of different parameters and their deviations on performance. The relationship between the cell operating temperature and voltage is negative correlation, which means that at the same current density, the deviation of voltage is negative and larger as the current density increases (see Figure 3a). The main reason is that when other parameters are kept constant, the lifting COT increases the saturated vapor pressure of water. The relative humidity (RH) of the reactant gas decreases, which weakens the hydration of the proton membrane, reduces the water content in the membrane, increases the ohmic resistance, and lowers the output voltage. Moreover, under the same deviation of temperature, as the current increases, increasingly more water is produced, which increases the ohmic losses. Therefore, the value of the amplitude of negative deviation at 0.6 A/cm 2 is the largest.

Deviation of Cell Operating Temperature
To quantify the effect of OCT fluctuation on the output voltage under three values of current density (0.1 A/cm 2 , 0.3 A/cm 2 , 0.6 A/cm 2 ) in Figure 3a, Figure 3b shows the voltage deviations (∆V) and their deviation percentages corresponding to the temperature deviations. It can be seen that the overall amplitude of fluctuation is smooth when the current density is 0.1 A/cm 2 . This is because when the current density is low, there is less water produced by electrochemical reaction and less water accumulated in the porous media layer, which will cause a decrease in the amplitude of performance changes. The lifting and declining amplitudes of the voltage are low compare with high current density. Figure 4 shows the impact of deviation of the cathode operating pressure (COP) on the deviation of the output voltage and the deviation percentage under various current densities. This is in agreement with the upper segment results. COP is positively correlated with the fuel cell output voltage at the same current density, that is, the voltage deviation is positive and increases the same as the increase of the COP at the same current density (see as Figure 4a). According to the data in Figure 4b, the change trend of output voltage with cathode pressure under three current densities is similar under three kinds of current density and the amplitudes of deviation under different cathode pressure increase with the increase of current density. The magnitude of deviation/percentage in Figure 4b indicates that the impact of cathode pressure fluctuation on PEMFC performance is second only to COT.

Deviation of Cathode Operating Pressure
The main reason for this phenomenon may be due to the supplied gas for the cathode inlet being air (oxygen only accounts for 21%). The fluctuated operating cathode pressure has a great influence on the oxygen concentration in cathode CL, which in turn affects the performance of PEMFC. As for why there are obvious differences of polarization curves among five kinds of COP at high current density, it may be explained that changing the pressure can change the limiting current density, which means that only the high current density (the region close to the limiting current density) can be significantly influenced.

Deviation of Anode Pressure
The impact of deviation of the anode operating pressure (AOP) on the polarization curves and the output deviation/deviation percentage depicted in Figure 5a is lower compared with Figure 5b. However, the fluctuating trends of the output voltage and the corresponding percentages shown in Figure 5b are basically the same as that in Figure 5a. Same as the theory of air transfer in the cathode, the increase in anode pressure promotes the convection diffusion of hydrogen molecules in the porous layer of the anode, which enhances the mass transfer performance of hydrogen, increases the exchange current density, and reduces the loss of concentration, finally improving the performance of the fuel cell. Meanwhile, the oxygen required in the cathode will increase to participate in the reaction with hydrogen in time. However, the increase of anode pressure will prevent the oxygen mass transfer in the cathode side, that is, the oxidation-reduction reaction will be hindered. Therefore, the influence of AOP fluctuation on the output voltage is inferior to the fluctuation of COP.  Figure 6 depicts when the calibration intake parameters are maintained constant, how the deflection in the anode mass flow rate (AMFR) impacts the deviation of the output voltage, and the percentage of deviation under various current densities. Figure 6a,c indicates that the relationship between the AMFR and the cell performance is a negative correlation, that is, at the same current density, the cell output voltage decreases as a function of the deviation of the AMFR. Meanwhile, the amplitude of the voltage deviation increases with an increase in the deviation of the AMFR. To explore the reason for the large changes in output voltage in Figure 7 when the current density is high, the cloud maps of water content distribution in the contact region between the anode CL and GDL when the fluctuations of AMFR are ±2.5%, ±5%, respectively, and current density is 0.6 A/cm 2 are shown in Figure 7. Although there is a minimal current density that occurs when the AMFR is 6 × 10 −7 kg/s, the overall distribution of water content within the interface seems to more evenly compared with the values of AMFR, which are 5.85 × 10 −7 kg/s and 5.7 × 10 −7 kg/s, respectively. Especially when the mass flow rate is 5.7 × 10 −7 kg/s, large areas of extreme water content are distributed in different areas of the interface. With the AMFR increasing, high water content areas occur accompanied with the increase of water content in the central areas of the interface. The anode is drier than the cathode due to redox reaction; the excessive low water content in anode will decrease the performance of PEMFC, so the performance of PEMFC can be improved by enhancing the anode water content.

Deviation of Cathode Mass Flow Rate
The effects of the fluctuations of cathode mass flow rate (CMFR) on the polarization curves are too small to be observed clearly. Comparing both sides of the y-axis in Figures 6b and 8, it can be found that the influence of the fluctuated CMFR on the output voltage is in the same order of magnitude as that of fluctuated AMFR. Moreover, the deviations caused by the fluctuated CMFR are not affected by current density, especially when the load current is relatively large. On the one hand, the improvement of air flow rate bonds to increase the oxygen concentration in CL, which will promote the redox reactions in three-phase areas. On the other hand, with the flow rate increasing, the water removal rate of the cathode side increases. Finally, the performance of PEMFC is reduced because of the dehydration of membrane. Therefore, the trend of the output voltage deviation depends on the dominant party.

Deviation of Anode/Cathode Inlet Temperature
The effect of deviations of AIT and CIT on the output voltage deviation and deviation percentage are shown in Figures 9 and 10, respectively. Both the AIT and CIT are negatively correlated with the output performance of the fuel cell, which means that at the same current density, the output voltage of the cell increases (decreases) with a decrease (increase) in the anode/cathode intake temperature.  The temperature of the reactants can affect the performance of the fuel cell by affecting the relative humidity of the inlet gas. According to the literature [42], the higher inlet temperature corresponds to higher water saturation vapor pressure, the change in absolute humidity is ignored, and the relative humidity of the inlet gas decreases, which ultimately leads to a decrease in PEMFC performance. It can be seen that the larger effect of the AIT deviation on the PEMFC performance, when the current density is 0.1 A/cm 2 from Figures 9 and 10, is because, compared with the cathode, the anode is drier and requires humidification.

Process of Sensitivity Analysis
Above contents discuss how the operating parameters of PEMFC within a certain range of positive and negative deviations impact the deviation of the output voltage and the percentage deviation. According to the pictures of polarization curves and parameter deviations present in Section 3.3, a comparison of the sensitivity of the operating parameters has been obtained. However, there is a nuance of the parameters' sensitivity rankings among the insensitivity indexes.
To further quantify the impact of fluctuation on the insensitive working parameters on the output voltage of the fuel cell under various working current densities, based on the Monte Carlo method, multiparameter sensitivity analysis (MSA) is used to obtain the insensitivity indexes of the output voltage caused by the operating conditions under various current densities. The process of MSA can be expressed as: (1) Determining the insensitive parameters: T an , T ca , S an , S ca , as seen in Table 4, and defining the deviation range of the standard value as ±5%. (2) Randomly generating 12 sets of uniform operating parameter values within the range of positive and negative deviation of each parameter, then bringing in the calculated model and obtaining the voltage, which corresponds to the current density as 0.1 A/cm 2 , 0.3 A/cm 2 , 0.6 A/cm 2 , respectively. (3) Substituting the output voltage from step 2) into the function. Finally, the sensitivity of static evaluation parameters was obtained from Equation (19).
(4) where f (i), K(i), X 0 (i), and X c (i) represent the sum of the output voltage deviation, the parameter sensitivities, the reference voltage, and calculated voltage, respectively, when the current density is i. k is the zoom factor and n refers to the number of simulation times when the current density is i.

Results of Analysis
According to the calculation of Equation (19), under different current densities, the sensitivity indexes obtained by four operating parameters within a certain range of deviation (±5%) are depicted in Figure 11. Under the current density that represents the three types of vehicle operating conditions, the sensitivity index of the output voltage was measured. The insensitivity indexes of the operating parameters were observed to be highly nonlinear under various working conditions. Among all the parameters, the working temperature of the fuel cell T cell is the most sensitive under three types of conditions (e.g., idle condition, intermediate condition, rated condition). Quantity differences of y-axis between Figures 4b and 5b show that the sensitivity index of cathode pressure P ca is far more than the sensitivity index of anode pressure P an . Moreover, with the increase in current density, both the sensitivity indexes of P ca and Pan improve. Even for these insensitive parameters (S an , S ca , T an , and T ca ), according to the z-axes of Figure 11a,b, the differences among these sensitivity indexes can be far more than 100 times. By comprehensive comparisons with Figures 3-10, it can be observed that the working parameters tend to have higher sensitivity when the PEMFC is in rated condition.
According to the sensitivity indexes, the parameter sensitivities under various conditions can be ranked as Table 6. Table 6. Parameter sensitivity rankings.

Operating Conditions Rankings
Idle condition T cell > P ca > P an > S an > S ca > T an > T ca Intermediate condition T cell > P ca > P an > S ca > S an > T ca > T an Rated condition T cell > P ca > P an > S an > S ca > T ca > T an The parameters are sorted into three groups, which are: very sensitive T cell ; normal sensitive P ca , P an ; insensitive T ca , T an , S ca , S an .

Cloud Maps Analysis
The sensitivities of operating parameters of PEMFC are evaluated by digitized form in previous contents. Some sensitive parameters: T cell , P ca , P an , which are evaluated in Section 3.4.2, are used to obtained cloud maps in this section. The effect of fluctuations of these three parameters on performance of PEMFC can reflect on the cloud map, which represents the water content distribution on the interface between membrane and anode CL. Analyzing these cloud maps can more directly explain the fluctuation of operating parameters about how to influence the PEMFC performance.
Due to the too long fuel cell body, only the same part of different cloud maps is intercepted to compare. It can be seen from Figure 12, as the cell temperature decreases, the water content in the interface between the membrane and anode-side CL continues to increase and the water distribution tends to be even. The water content of interfaces increases with the increase of inlet pressure, both in Figures 13 and 14. During the reaction process of PEMFC, water flows from the anode side to cathode side under the impact of electroosmotic migration, which causes the lack of water on the anode side of PEMFC. There is a sharp drop of PEMFC performance when the water content on the anode side is too low; only when the drying effect of the anode side is alleviated can the performance be improved. Moreover, under the fluctuation of the very sensitive parameter T cell , there is the largest difference in the water content distribution on the anode side. The cloud map of water content distribution becomes more and more similar under the fluctuation of parameters with lower sensitivity index. The cloud maps of water content distribution in

Conclusions
A set of unified parameter fluctuation ranges and evaluation methods, which were combined with the actual operating conditions of PEMFC, has been constructed in this study to avoid the evaluation errors caused by inconsistent standards and to make the results meaningful. The simulation was conducted in a 3D, non-isothemal PEMFC model. Furthermore, a multiparameter sensitivity analysis method combined with the Monte Carlo method was used to evaluate the sensitivity of the working parameter fluctuations (±5%) to the output performance. Among these parameters, the temperature of the cell is the most sensitive under all working conditions and the sensitivity far exceeds other parameters. However, within the range of deviation, the inlet temperature in both the anode and the cathode has a minimal impact on the output voltage. Moreover, with the increase in current density, the sensitivity indexes of the working pressures in both the anode and the cathode increase, whereas those of the mass flow rate of the anode/cathode change disorderly. Generally, the cathode parameters are more sensitive than the anode parameters, especially under high current density. These results provide a reference for the performance optimization of the PEMFC and a control strategy for the operational parameters.
In future studies, the interaction between the parameters will be considered in order to analyze a more complete working process of the PEMFC.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature and Abbreviation
List