The 2013 two-dimensional model of Yin, Hsuen and Hsu computes liquid-water saturation in the cathode gas diffusion layer (GDL) of a proton exchange membrane fuel cell (PEMFC) from a pseudo-phase-equilibrium closure, but it is cathode-only, isothermal, treats membrane water as a fixed boundary flux, and rests on an empirical smoothing parameter delta fitted as 0.03. The authors themselves list the unresolved directions toward a comprehensive simulation tool: coupling membrane hydration with variable electro-osmotic drag and back-diffusion, a non-isothermal formulation, inclusion of ohmic loss, GDL compression and finite catalyst-layer thickness, a physical justification of delta in terms of GDL material properties, and experimental validation including the anode side.

Background

A proton exchange membrane fuel cell (PEMFC) converts hydrogen and oxygen to electricity through a polymer electrolyte membrane that conducts protons but not electrons. Water management is its central design conflict: the membrane must stay hydrated to conduct protons, but water produced at the cathode and dragged through the membrane by proton flux can flood the porous carbon-paper gas diffusion layer (GDL) that carries reactant gases to the catalyst layer (CL), starving the cell of oxygen. Predicting this balance requires two-phase (vapor plus liquid) porous-media models in which liquid saturation - the local volume fraction of pores filled with liquid water - follows from coupled vapor diffusion, capillary transport, and evaporation/condensation.

Mainstream two-phase models close the phase-change term either by finite-rate evaporation/condensation kinetics (non-equilibrium models, e.g. Khajeh-Hosseini-Dalasm et al. 2010; Ramos-Alvarado et al. 2013) or by assuming local thermodynamic equilibrium between vapor and liquid. The finite-rate route has a well-known weakness: the interfacial mass-transfer coefficients for hydrophobic GDLs are unknown and treated as adjustable constants. The pseudo-phase-equilibrium approach of Hsuen and Yin (2011), extended to two dimensions by Yin, Hsuen and Hsu (2013), takes the equilibrium route in smoothed form: instead of resolving interfacial kinetics, liquid saturation is obtained from a phase-equilibrium function (Eq. (5) of the 2013 paper) that spreads the sharp vapor-liquid transition over a window controlled by an empirical parameter delta (a fitted value of 0.03 was used), so that saturation emerges as a continuous function of the local degree of oversaturation of the water vapor. Their 2013 model solved the cathode GDL in two dimensions under varying feed humidification, but with acknowledged simplifications: the GDL is isothermal and isobaric, the membrane is not modeled (water crossing it enters only as a boundary flux with a constant electro-osmotic drag coefficient n_d = 1.59; electro-osmotic drag is the transport of water molecules bound to protons moving through the membrane), the catalyst layer is infinitesimally thin, GDL compression and electronic (carbon-phase) ohmic losses are neglected, and no direct comparison with measured polarization curves (voltage-current characteristics) or in-plane water distributions was made.

The terminology for the missing pieces: back-diffusion is the diffusive transport of water through the membrane from cathode to anode that partially counteracts electro-osmotic drag; a non-isothermal formulation would let local temperature vary, which matters because the saturated vapor mole fraction depends exponentially on temperature (Clausius-Clapeyron relation) so that small thermal gradients strongly redistribute condensation. Neighboring approaches that already treat parts of this territory include non-equilibrium two-phase full-cell models (e.g. the pseudo-three-dimensional, two-phase, non-isothermal model of Rizvandi and Yesilyurt 2019, and the many models reviewed by Jiao and Li 2011) and pore-network models, which discretize the actual GDL pore space and resolve condensation and capillary transport at the pore scale (Hinebaugh and Bazylak 2010; Straubhaar et al. 2015; Carrère and Prat 2019), including validation against X-ray tomographic water distributions (Agaesse et al. 2016) and in-operando synchrotron imaging studies of water in operating GDLs (Arlt et al. 2017). These provide both the competing closures and the experimental data that any comprehensive extension of the pseudo-phase-equilibrium framework would have to engage.

Problem Statement

As stated by the source authors, the unresolved problem is to extend the two-dimensional pseudo-phase-equilibrium model of the cathode GDL (Yin, Hsuen and Hsu, 2013) toward a comprehensive, predictive PEMFC water-management simulation tool. Concretely, the model must be developed so that: (1) membrane hydration/dehydration is coupled to the cathode GDL solution - water transport across the membrane treated beyond a fixed boundary flux, with back-diffusion and hydration-dependent electro-osmotic drag replacing the constant coefficient nd=1.59n_d = 1.59 - and the anode side is included so the complete cell water balance is predicted; (2) the formulation is non-isothermal, with the saturated vapor mole fraction entering the pseudo-phase-equilibrium function expressed as a function of local temperature, so that thermally driven phase change and transport at high current densities are captured; (3) the neglected real-geometry effects - ohmic loss in the carbon (GDL solid) phase, deformation of the GDL under compression, and the finite thickness of the catalyst layer - are represented; (4) the adjustable smoothing parameter δ\delta of the pseudo-phase-equilibrium function (fitted as δ=0.03\delta = 0.03) is given a physical justification or a predictive correlation in terms of GDL material properties such as pore-size distribution, wettability, and porosity; and (5) the resulting model is directly validated against experimental polarization and water-distribution data under controlled humidification. These components serve one common objective: a physically grounded, experimentally validated water-management model built on the pseudo-phase-equilibrium closure, and all five are jointly required for it - completion is judged on the integrated model, so delivering a strict subset of the components (for example, the grounded δ\delta closure alone or an uncoupled cathode treatment) does not by itself resolve the task.

References

  1. Ken-Ming Yin, Hsiao-Kuo Hsuen, Chih-Hsuan Hsu (2013). Two-dimensional pseudo-phase-equilibrium approach on the cathode gas diffusion layer of proton exchange membrane fuel cell under varying degrees of humidification. Journal of Power Sources 234, 223-233. DOI: 10.1016/j.jpowsour.2013.01.168. https://doi.org/10.1016/j.jpowsour.2013.01.168
  2. Hsiao-Kuo Hsuen, Ken-Ming Yin (2011). A pseudo-phase-equilibrium approach for the calculation of liquid water saturation in the cathode gas diffuser of proton-exchange-membrane fuel cells. International Journal of Hydrogen Energy 36, 5487-5499. DOI: 10.1016/j.ijhydene.2011.02.047. https://doi.org/10.1016/j.ijhydene.2011.02.047
  3. Hsiao-Kuo Hsuen, Ken-Ming Yin (2012). Performance equations of proton exchange membrane fuel cells with feeds of varying degrees of humidification. Electrochimica Acta 62, 447-460. DOI: 10.1016/j.electacta.2011.12.079. https://doi.org/10.1016/j.electacta.2011.12.079
  4. Hsiao-Kuo Hsuen, Ken-Ming Yin (2011). Performance equations for a polymer electrolyte membrane fuel cell with unsaturated cathode feed. Journal of Power Sources 196, 218-227. DOI: 10.1016/j.jpowsour.2010.06.058. https://doi.org/10.1016/j.jpowsour.2010.06.058
  5. K.-M. Yin, H.-K. Hsuen (2013). Mathematical Model of Proton Exchange Membrane Fuel Cell with Consideration of Water Management. Fuel Cells 13, 1213-1225. DOI: 10.1002/fuce.201300006. https://doi.org/10.1002/fuce.201300006
  6. Ken-Ming Yin, Hsiao-Kuo Hsuen, Yu-An Kung (2014). A comparison of numerical methods on the two-phase flow in PEMFCs. Journal of the Taiwan Institute of Chemical Engineers 45, 2358-2369. DOI: 10.1016/j.jtice.2014.07.001. https://doi.org/10.1016/j.jtice.2014.07.001
  7. Ken-Ming Yin, Hsiao-Kuo Hsuen, Yu-An Kung, Bo-Syun Cheng (2016). One-dimensional comparison of numerical approaches on two-phase flow in the membrane electrode assembly of PEMFC. Journal of the Taiwan Institute of Chemical Engineers 63, 133-142. DOI: 10.1016/j.jtice.2016.02.021. https://doi.org/10.1016/j.jtice.2016.02.021
  8. Ken-Ming Yin, Bo-Syun Cheng, Kuo-Wei Chiang (2016). Non-uniform agglomerate cathode catalyst layer model on the performance of PEMFC with consideration of water effect. Renewable Energy 95, 191-201. DOI: 10.1016/j.renene.2016.04.015. https://doi.org/10.1016/j.renene.2016.04.015
  9. Sanggyu Kang (2015). Quasi-three dimensional dynamic modeling of a proton exchange membrane fuel cell with consideration of two-phase water transport through a gas diffusion layer. Energy 90, 1388-1400. DOI: 10.1016/j.energy.2015.06.076. https://doi.org/10.1016/j.energy.2015.06.076
  10. Omid Babaie Rizvandi, Serhat Yesilyurt (2019). A pseudo three-dimensional, two-phase, non-isothermal model of proton exchange membrane fuel cell. Electrochimica Acta 302, 180-197. DOI: 10.1016/j.electacta.2019.02.018. https://doi.org/10.1016/j.electacta.2019.02.018
  11. Kui Jiao, Xianguo Li (2011). Water transport in polymer electrolyte membrane fuel cells. Progress in Energy and Combustion Science 37, 221-291. DOI: 10.1016/j.pecs.2010.06.002. https://doi.org/10.1016/j.pecs.2010.06.002
  12. N. Khajeh-Hosseini-Dalasm, Kazuyoshi Fushinobu, Ken Okazaki (2010). Phase change in the cathode side of a proton exchange membrane fuel cell. Journal of Power Sources 195, 7003-7010. DOI: 10.1016/j.jpowsour.2010.04.089. https://doi.org/10.1016/j.jpowsour.2010.04.089
  13. Bladimir Ramos-Alvarado, Abel Hernandez-Guerrero, Michael W. Ellis (2013). Non-equilibrium two-phase model of the air-cathode of a PEM fuel cell based on GDL experimental water transport characteristics. Journal of Power Sources 232, 376-388. DOI: 10.1016/j.jpowsour.2012.10.091. https://doi.org/10.1016/j.jpowsour.2012.10.091
  14. J. Hinebaugh, A. Bazylak (2010). Condensation in PEM Fuel Cell Gas Diffusion Layers: A Pore Network Modeling Approach. Journal of The Electrochemical Society 157, B1382. DOI: 10.1149/1.3467837. https://doi.org/10.1149/1.3467837
  15. Benjamin Straubhaar, Joel Pauchet, Marc Prat (2015). Water transport in gas diffusion layer of a polymer electrolyte fuel cell in the presence of a temperature gradient. Phase change effect. International Journal of Hydrogen Energy 40, 11668-11675. DOI: 10.1016/j.ijhydene.2015.04.027. https://doi.org/10.1016/j.ijhydene.2015.04.027
  16. Tristan Agaesse, Adrien Lamibrac, Felix N. Büchi, Joel Pauchet, Marc Prat (2016). Validation of pore network simulations of ex-situ water distributions in a gas diffusion layer of proton exchange membrane fuel cells with X-ray tomographic images. Journal of Power Sources 331, 462-474. DOI: 10.1016/j.jpowsour.2016.09.076. https://doi.org/10.1016/j.jpowsour.2016.09.076
  17. Pierre Carrère, Marc Prat (2019). Impact of non-uniform wettability in the condensation and condensation-liquid water intrusion regimes in the cathode gas diffusion layer of proton exchange membrane fuel cell. International Journal of Thermal Sciences 145, 106045. DOI: 10.1016/j.ijthermalsci.2019.106045. https://doi.org/10.1016/j.ijthermalsci.2019.106045
  18. Tobias Arlt, Merle Klages, Matthias Messerschmidt, Joachim Scholta, Ingo Manke (2017). Influence of artificially aged gas diffusion layers on the water management of polymer electrolyte membrane fuel cells analyzed with in-operando synchrotron imaging. Energy 118, 502-511. DOI: 10.1016/j.energy.2016.10.061. https://doi.org/10.1016/j.energy.2016.10.061