Fang, Li and Tao (2021) built a two-dimensional pore-scale lattice Boltzmann framework for vapor generation, condensation, and liquid-water transport across the entire cathode of a proton exchange membrane fuel cell (PEMFC), but with a prescribed temperature field, one-way coupling from reaction to water, and idealized microstructures. The open problem is the fully coupled three-dimensional version: a pore-scale cathode model that solves the energy equation, two-way couples multicomponent transport, phase change, and electrochemical reaction in a realistic heterogeneous microstructure, applied to investigate liquid-water evolution (vapor generation, condensation, transport, and flooding) in the cathode under that coupling, including the effects of temperature distribution, GDL wettability, and microstructure. Each coupling element exists in a separate published model, but none combines them, so the problem remains open.

Background

A proton exchange membrane fuel cell (PEMFC) generates electricity by oxidizing hydrogen at an anode and reducing oxygen at a cathode, separated by a proton-conducting membrane. The cathode is where performance is usually lost. It comprises a catalyst layer (CL, a few micrometres of carbon-supported platinum and ionomer), a gas diffusion layer (GDL, a 100-300 micrometre carbon-fiber paper or cloth, usually rendered hydrophobic by polytetrafluoroethylene, PTFE), sometimes a microporous layer (MPL) between the two, and gas flow channels bounded by solid ribs (lands) that press on the GDL. The oxygen reduction reaction (ORR) at the cathode produces water vapor; if that vapor condenses and accumulates as liquid it blocks oxygen pathways (flooding) and covers catalyst sites, whereas too little water dries the membrane. Managing this balance is a central PEMFC design problem, and it is intrinsically thermal: near the normal operating temperature of 80 °C, small temperature shifts strongly change where vapor condenses, and the liquid-to-vapor density ratio is of order 2000.

Continuum two-phase fuel-cell models treat the porous layers with Darcy-type averaged equations that require semi-empirical closures, such as capillary pressure-saturation curves, which are poorly known for GDL materials. Pore-scale models instead resolve the actual solid microstructure and the interfaces inside it. The lattice Boltzmann method (LBM), which evolves particle distribution functions on a lattice rather than discretizing the Navier-Stokes equations directly, dominates this field because it handles complex geometry and moving liquid-gas interfaces naturally; in the pseudopotential (Shan-Chen) variant, mesoscopic interaction forces produce liquid-vapor coexistence and phase change, and wall forces encode wettability (contact angle). Realistic liquid/vapor density ratios are demanding for such methods, and the coexistence behavior must be verified against an equation of state (for example Peng-Robinson with a Maxwell construction).

Against this background, Fang, Li and Tao (2021) developed a two-dimensional pore-scale LBM framework for the entire cathode (CL, GDL, and gas channel) coupling: (i) a single-component multiphase water model with condensation at a density ratio of about 2000 at operating temperature; (ii) oxygen transport with reactive boundary conditions on the CL pore surfaces, combining diffusive consumption and reaction-induced advective flux; and (iii) a prescribed temperature field, either uniform or an imposed through-plane gradient. Over idealized regular and reconstructed GDL microstructures with assigned wettability patterns, they concluded that a cooler gas-channel side moves condensation away from the CL and relieves CL flooding, that a uniform temperature concentrates condensation inside the CL, that condensed liquid throttles oxygen transport and deactivates covered catalyst sites (oxygen being sparingly soluble in liquid water), and that GDL wettability decides whether condensate drains to the channel or is trapped.

The authors themselves delimited what the study does not do: the domain is two-dimensional; the temperature field is prescribed instead of emerging from a solved energy equation with conjugate heat transfer; the water dynamics is single-component, so the influence of oxygen on the dynamic transport of water is neglected and the reaction-to-water coupling is one-way; and the microstructures are idealized. The open question is the fully coupled three-dimensional problem that removes these restrictions simultaneously.

Problem Statement

Develop, validate, and apply a pore-scale model of the complete PEMFC cathode (catalyst layer, GDL, with microporous layer if present, and gas channel, including the rib where relevant) that simultaneously includes, in three spatial dimensions:

  1. a solved energy equation, so that the temperature field emerges from conjugate heat transfer between solid and fluid phases and from reaction heat, instead of being prescribed;
  2. multicomponent gas transport (at least oxygen and water vapor in a carrier gas) that is two-way coupled to the liquid-water dynamics, meaning the influence of oxygen transport and of the electrochemical reaction on water generation, condensation, and motion, as well as the blocking of oxygen by liquid, are all represented, removing the one-way simplification;
  3. dynamic liquid-vapor phase change (condensation and evaporation) at the density ratio of PEMFC operating conditions (of order 1000-2000 near 80 °C);
  4. the oxygen reduction reaction resolved on the catalyst-layer pore surfaces, not collapsed into a boundary condition of a homogenized medium;
  5. a realistic static heterogeneous microstructure, imaged or stochastically reconstructed, with spatially varying wettability, rather than idealized regular structures.

Using such a model, investigate liquid-water evolution (vapor generation, condensation, transport, and flooding) in the cathode under this fully coupled three-dimensional setting and under more general operating conditions and microstructural heterogeneity, including the effects of temperature distribution, GDL wettability, and microstructure. The conclusions of the two-dimensional prescribed-temperature study of Fang, Li and Tao (2021) (relocation of condensation and flooding by through-plane temperature gradients, the trapping-versus-drainage role of GDL wettability, and the deactivation of liquid-covered catalyst sites), summarized in the Background, motivate this application.

The problem is complete when a model containing all the couplings above has been developed, validated in the ordinary sense of pore-scale multiphase modeling, and applied to investigate liquid-water evolution in the cathode under the fully coupled three-dimensional setting, covering the effects of temperature distribution, GDL wettability, and microstructure. Partial models that omit one or more of the couplings, for example a 3D isothermal two-phase reactive model or a non-isothermal model without phase change, repeat known component work and do not answer the question.

References

  1. Fang, Wen-Zhen; Li, Jin; Tao, Wen-Quan (2021). Two-dimensional pore-scale investigation of liquid water evolution in the cathode of proton exchange membrane fuel cells. Numerical Heat Transfer, Part A: Applications 79(4):261-277. DOI: 10.1080/10407782.2020.1845558. https://doi.org/10.1080/10407782.2020.1845558
  2. Kopanidis, A.; Theodorakakos, A.; Gavaises, M.; Bouris, D. (2011). Pore scale 3D modelling of heat and mass transfer in the gas diffusion layer and cathode channel of a PEM fuel cell. International Journal of Thermal Sciences 50:456-467. DOI: 10.1016/j.ijthermalsci.2010.11.014. https://doi.org/10.1016/j.ijthermalsci.2010.11.014
  3. Molaeimanesh, G.R.; Akbari, M.H. (2014). A three-dimensional pore-scale model of the cathode electrode in polymer-electrolyte membrane fuel cell by lattice Boltzmann method. Journal of Power Sources 258:89-97. DOI: 10.1016/j.jpowsour.2014.02.027. https://doi.org/10.1016/j.jpowsour.2014.02.027
  4. Molaeimanesh, Gholam Reza; Akbari, Mohammad Hadi (2015, online 2014). A pore-scale model for the cathode electrode of a proton exchange membrane fuel cell by lattice Boltzmann method. Korean Journal of Chemical Engineering 32:397-405. DOI: 10.1007/s11814-014-0229-6. https://doi.org/10.1007/s11814-014-0229-6
  5. Molaeimanesh, Gholam Reza; Dahmardeh, Masoud (2021). Pore-scale analysis of a PEM fuel cell cathode including carbon cloth gas diffusion layer by lattice Boltzmann method. Fuel Cells 21:208-220. DOI: 10.1002/fuce.202000191. https://doi.org/10.1002/fuce.202000191
  6. Jithin, M.; Siddharth, Saurabh; Das, Malay K.; De, Ashoke (2017). Simulation of coupled heat and mass transport with reaction in PEM fuel cell cathode using lattice Boltzmann method. Thermal Science and Engineering Progress 4:85-96. DOI: 10.1016/j.tsep.2017.09.004. https://doi.org/10.1016/j.tsep.2017.09.004
  7. Chen, Li; Kang, Qinjun; He, Yaling; Tao, Wenquan; Zelenay, Piotr (2013). Pore-Scale Simulation of Water Vapor Generation, Condensation, and Transport in the Cathode of a Proton Exchange Membrane Fuel Cell. ECS Meeting Abstracts MA2013-02:1254. DOI: 10.1149/ma2013-02/15/1254. https://doi.org/10.1149/ma2013-02/15/1254
  8. Ye, Shengjie; Hou, Yuze; Li, Xing; Jiao, Kui; Du, Qing (2023, online 2022). Pore-Scale Investigation of Coupled Two-Phase and Reactive Transport in the Cathode Electrode of Proton Exchange Membrane Fuel Cells. Transactions of Tianjin University 29:1-13. DOI: 10.1007/s12209-021-00309-4. https://doi.org/10.1007/s12209-021-00309-4
  9. Guo, Lingyi; Chen, Li; Zhang, Ruiyuan; Peng, Ming; Tao, Wen-Quan (2022). Pore-scale simulation of two-phase flow and oxygen reactive transport in gas diffusion layer of proton exchange membrane fuel cells: Effects of nonuniform wettability and porosity. Energy 253:124101. DOI: 10.1016/j.energy.2022.124101. https://doi.org/10.1016/j.energy.2022.124101
  10. Zhou, Chongbo; Guo, Lingyi; Chen, Li; Tian, Xin; He, Tiefeng; Yang, Qinghua (2021). Pore-Scale Modeling of Air-Water Two Phase Flow and Oxygen Transport in Gas Diffusion Layer of Proton Exchange Membrane Fuel Cell. Energies 14(13):3812. DOI: 10.3390/en14133812. https://doi.org/10.3390/en14133812
  11. Wang, Hao; Yang, Xiaoxing; Yang, Guogang; Zhang, Guoling; Li, Zheng; Li, Lingquan; Huang, Naibao (2024). Impact of Gas Diffusion Layer Compression on Electrochemical Performance in Proton Exchange Membrane Fuel Cells: A Three-Dimensional Lattice Boltzmann Pore-Scale Analysis. Nanomaterials 14(24):2012. DOI: 10.3390/nano14242012. https://doi.org/10.3390/nano14242012
  12. Jiang, Ziheng; Yang, Guogang; Li, Shian; Shen, Qiuwan; Liao, Jiadong; Wang, Hao; Zhang, Guoling; Li, Zheng; Liu, Zheng (2023). Pore scale study of heat and mass transfer in ice-containing gas diffusion layer of polymer electrolyte membrane fuel cells. Journal of Power Sources 580:233347. DOI: 10.1016/j.jpowsour.2023.233347. https://doi.org/10.1016/j.jpowsour.2023.233347
  13. Jiang, Ziheng; Yang, Hua; Huang, Xueli; Qi, Cong; Zheng, Haitao; Liao, Jiadong (2024). Study on heat and mass transfer behavior of coupled electrochemical reaction in porous structure of proton exchange membrane fuel cell. Materials Today Communications 41:110536. DOI: 10.1016/j.mtcomm.2024.110536. https://doi.org/10.1016/j.mtcomm.2024.110536
  14. Grunewald, Jonathan B.; Goswami, Navneet; Mukherjee, Partha P.; Fuller, Thomas F. (2021). Two-Phase Dynamics and Hysteresis in the PEM Fuel Cell Catalyst Layer with the Lattice-Boltzmann Method. Journal of The Electrochemical Society 168:024521. DOI: 10.1149/1945-7111/abe5e8. https://doi.org/10.1149/1945-7111/abe5e8
  15. Goswami, Navneet; Grunewald, Jonathan B.; Fuller, Thomas F.; Mukherjee, Partha P. (2024). Probing multiphase reactive transport interactions in the polymer electrolyte fuel cell catalyst layer degradation. Electrochimica Acta 486:144143. DOI: 10.1016/j.electacta.2024.144143. https://doi.org/10.1016/j.electacta.2024.144143
  16. Belgacem, Najib; Prat, Marc; Pauchet, Joel (2017). Coupled continuum and condensation-evaporation pore network model of the cathode in polymer-electrolyte fuel cell. International Journal of Hydrogen Energy 42:8150-8165. DOI: 10.1016/j.ijhydene.2017.01.184. https://doi.org/10.1016/j.ijhydene.2017.01.184
  17. Carrère, Pierre; Prat, Marc (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. Lai, Tao; Qu, Zhiguo; Zhang, Jianfei (2025). Pore-scale study of liquid water transport in gas diffusion layers with in-plane non-uniform distributed pore size of polymer electrolyte membrane fuel cell. Applied Energy 379:124933. DOI: 10.1016/j.apenergy.2024.124933. https://doi.org/10.1016/j.apenergy.2024.124933