Fuel Cell Modeling

Category: Coupled Analysis – Electrochemistry | Revised 2026-10-01

Theory: theoretical voltage minus three losses

Overview

🙋

Why does a fuel cell deliver only about 0.6 V rather than the theoretical 1.2 V?

🎓

Because current flow incurs several losses. First, part of the voltage drives electrode reactions (activation loss); in PEM cells the cathode oxygen reaction is slow, costing about 0.4 V. Next, resistance in the membrane and between parts drops voltage (ohmic loss). And at high current oxygen and hydrogen can't reach reaction sites fast enough, so voltage falls sharply (concentration loss). Summing these against current density gives the polarization curve, the basis of cell design and operating choices. All lost voltage becomes heat, which also follows from this curve.

Polarization-curve model

$$ V = E_0 - b\ln\frac{i}{i_0} - i\,R_{ASR} + c\,\ln\!\left(1 - \frac{i}{i_L}\right),\qquad b = \frac{RT}{\alpha F} $$

E₀ is theoretical voltage (1.19 V at 80 °C), b the Tafel slope (0.0304 V with α = 1, 70 mV per decade), i₀ exchange current density (10⁻⁶ A/cm²), R_ASR area-specific resistance (0.10 Ω·cm²), i_L limiting current density (2.0 A/cm²) and c the concentration coefficient (0.05 V).

Power, efficiency and heat

$$ P = V i,\qquad \eta_{LHV} = \frac{V}{1.25\ \mathrm{V}},\qquad Q_{heat} = (1.25 - V)\,i\ \ (\text{per area, water as vapour}) $$
Coffee Break Trivia

From spacecraft to cars

Grove demonstrated the fuel-cell principle in 1839, but practical use came with 1960s space programmes: Gemini carried PEM fuel cells and Apollo alkaline ones, supplying electricity and drinking water. Membrane and catalyst advances later raised power density, and fuel-cell cars went on sale in 2014. Cutting platinum catalyst while limiting cathode activation loss remains a major research challenge.

Worked examples

Example 1: voltage and losses by current density (80 °C)

Current densityCell voltageActivationOhmicConcentrationPowerEfficiency (LHV)
0.1 A/cm²0.827 V0.350 V0.010 V0.003 V0.083 W/cm²66%
0.5 A/cm²0.726 V0.399 V0.050 V0.014 V0.363 W/cm²58%
1.0 A/cm²0.635 V0.420 V0.100 V0.035 V0.635 W/cm²51%
1.5 A/cm²0.538 V0.433 V0.150 V0.069 V0.807 W/cm²43%
1.8 A/cm²0.457 V0.438 V0.180 V0.115 V0.822 W/cm²37%

Peak power is 0.831 W/cm² at 1.70 A/cm² (0.488 V).

Example 2: a drying membrane raises resistance

Current densityASR 0.10 Ω·cm²ASR 0.15 Ω·cm²
1.0 A/cm²0.635 V0.585 V
1.5 A/cm²0.538 V0.463 V

Delivering 80 kW at 1.0 A/cm² (0.635 V) needs about 12.6 m² of total cell area and releases about 78 kW of heat in operation (water leaving as vapour).

🙋

If it releases about as much heat as power, cooling must be as hard as for an engine car.

🎓

Less heat than an engine, but in some ways harder to handle. Fuel cells run near 80 °C, so the temperature difference to ambient is small and rejecting the same heat needs a large radiator. Engines dump much heat in exhaust, but fuel-cell exhaust is cool, so most heat must go through coolant. Temperature is also tightly tied to membrane water: too hot and the membrane dries, resistance rises and, as in Example 2, voltage drops and heating grows; too cold and product water condenses, blocking channels and raising concentration loss. Fuel cell models therefore couple electrochemical reactions, heat transfer and water transport (vapour and liquid) to resolve distributions within the cell.

Modelling workflow

  1. Fit exchange current density, resistance and limiting current to measured polarization curves.
  2. Include channels, diffusion layers, catalyst layers and membrane to solve gas and water flow.
  3. Derive heat generation from current distribution and solve temperature including cooling.
  4. Update membrane resistance with temperature and water content, coupling electrochemistry with heat and water.
  5. Compare predicted current, temperature and water distributions with measurements.
Coffee Break Trivia

“A cell whose outlet side underperformed”

Current-distribution measurements on one fuel cell showed current falling sharply only near the air outlet. Coupled analysis revealed humidity rising along the channel from product water until it condensed near the outlet, flooding the diffusion layer and starving oxygen, while the inlet ran dry with high membrane resistance. Revising air flow, humidification and coolant direction to even out temperature and humidity made current uniform and raised cell voltage.

Common mistakes

Mistakes and fixes

MistakeEffectFix
Constant temperatureDrying and flooding missedCouple heat and water
Constant membrane resistanceWrong voltage at high currentFunction of water content
Average current onlyLocal degradation missedResolve distributions
Not deriving heat from lossesInsufficient coolingHeat from voltage loss
Fitting a single pointWrong at other currentsFit the whole curve
🙋

I'd like to learn related topics.

Related Simulators

Experience the theory with interactive simulators in this field

All Simulators

Related fields

Structural AnalysisElectromagnetic Field AnalysisThermal Analysis
Rate this article
Thank you for your feedback!
Helpful
More details
Report error
Helpful
0
More details
0
Report error
0