“…This procedure can be used to optimize the GDL porous structure, flow field patterns, and operation conditions of PEFCs. Liquid-gas two phase modeling [17,18], modeling, modeling of microporous layers [19][20][21], and expansion to 3D modelings are future studies.…”