“…The computational domain is meshed using a structured hexahedral mesh and the relative residuals of all variables are set to 1:0 Â 10 À5 . The numerical simulation was carried out on Cathode reference exchange current density, A m À2 i ref,c 10,000 [9] Anode and cathode transfer coefficient α a , α c 0.5 [9] Anode and cathode activation energy, kJ mol À1 E act,a , E act,a 62.836, 24.359 [22] Contact angle, θ 80 [35] Surface tension, N m À1 σ 0.0625 [35] Porosity (GDL, CL) in anode ε GDL , ε CL 0.5, 0.25 [34] Porosity (GDL, CL) in cathode ε GDL , ε CL 0.5, 0.25 [34] Electrolyte volume fraction in CLs ε m 0.26 [27] Absolute permeability of GDLs and CLs, m 2…”