In protonic ceramic fuel cells (PCFCs), the H2 and H2O partial pressures in the anode side vary along the gas flow direction, particularly during high-fuel-utilization operation. This results in different species conductivities and overpotentials, which renders the quantitative evaluation of local performance based on H2 and H2O mole fractions important. A numerical model is developed that can reflect changes in species conductivities and overpotentials resulting from varying gas mole fractions. The numerical current density–voltage characteristic results obtained numerically agree well with those obtained via measurement in a wide range of H2–H2O gas conditions. Numerical modeling confirms a tradeoff, i.e., higher supply humidity levels result in lower electrolyte resistances but higher concentration overpotentials. The maximum output is achieved at a modest humidification level of 20%. The developed numerical model is particularly useful for optimizing the operating conditions of PCFCs with high fuel utilization.