Electrochemical and Physical Principles

Electrochemical and Physical Principles#

Thermodynamics of Water Electrolysis#

Water electrolysis splits water into hydrogen and oxygen using electrical energy:

\[ \mathrm{H_2O(l)} \rightarrow \mathrm{H_2(g)} + \tfrac{1}{2}\,\mathrm{O_2(g)} \]

The standard Gibbs free energy change (\(\Delta G^{0}\)) and enthalpy (\(\Delta H^{0}\)) for this reaction at \(25\,^{\circ}\mathrm{C}\) are approximately:

\[ \Delta G^{0} = 237.2\,\mathrm{kJ/mol} \]
\[ \Delta H^{0} = 285.8\,\mathrm{kJ/mol} \]

The minimum required cell voltage (reversible voltage) is:

\[ E_{\text{rev}}^{0} = \frac{\Delta G^{0}}{nF} \approx 1.23\,\mathrm{V} \]

The thermoneutral voltage is:

\[ E_{\text{tn}}^{0} = \frac{\Delta H^{0}}{nF} \approx 1.48\,\mathrm{V} \]

The thermoneutral voltage depends on temperature (see the original schematic for its variation).

Nernst Equation#

The reversible cell voltage under non-standard conditions is given by the Nernst equation:

\[ E = E^{0} + \frac{RT}{nF} \ln\!\left( \frac{a_{\mathrm{H_2}}\, a_{\mathrm{O_2}}^{1/2}}{a_{\mathrm{H_2O}}} \right) \]

Assuming liquid water activity \(a_{\mathrm{H_2O}} = 1\) and using gas partial pressures, this simplifies to:

\[ E = E^{0} + \frac{RT}{nF} \ln\!\left( p_{\mathrm{H_2}}\, p_{\mathrm{O_2}}^{1/2} \right) \]

Faradaic and Energy Efficiency#

Faradaic efficiency (\(\eta_F\)) relates delivered charge to the amount of product:

\[ \eta_F = \frac{nF Q_{\text{actual}}}{I \cdot t} \]

Energy efficiency compares the energy content of hydrogen to electrical input:

\[ \eta_E = \frac{\Delta G^{0}}{U_{\text{cell}} \cdot nF} \]

Electrochemical Kinetics and Overpotentials#

Real electrolyzer operation requires a cell voltage above the reversible value because of kinetic, ohmic, and transport losses:

\[ U_{\text{cell}} = E_{\text{rev}} + \eta_{\text{act,a}} + \eta_{\text{act,c}} + \eta_{\Omega} + \eta_{\text{mt}} \]

where:

  • \(\eta_{\text{act,a}}\) and \(\eta_{\text{act,c}}\) are activation overpotentials (anode and cathode)

  • \(\eta_{\Omega}\) is the ohmic overpotential (membrane, contacts, plates)

  • \(\eta_{\text{mt}}\) is the mass-transport or concentration overpotential

Butler-Volmer Kinetics#

At each electrode, current density and activation overpotential are linked by the Butler-Volmer equation:

\[ j = j_0 \left[\exp\!\left(\frac{\alpha_a F\eta_{\text{act}}}{RT}\right) - \exp\!\left(-\frac{\alpha_c F\eta_{\text{act}}}{RT}\right)\right] \]

with exchange current density \(j_0\) and transfer coefficients \(\alpha_a,\alpha_c\).

For PEM water electrolysis, the oxygen evolution reaction (OER, anode) is usually the dominant kinetic limitation, while the hydrogen evolution reaction (HER, cathode) is faster.

Tafel Approximation#

At moderate-to-high overpotential, Butler-Volmer simplifies to a Tafel form:

\[ \eta_{\text{act}} \approx \frac{2.303RT}{\alpha F}\log_{10}\!\left(\frac{j}{j_0}\right) \;=\; b\,\log_{10}\!\left(\frac{j}{j_0}\right) \]

where \(b\) is the Tafel slope (V/dec). This form is commonly used for parameter identification from polarization curves.

Ohmic Overpotential#

Ohmic loss is often modeled as:

\[ \eta_{\Omega} = j\,R_{\text{ASR}} \]

where \(R_{\text{ASR}}\) is the area-specific resistance (membrane plus electronic/contact contributions). Membrane hydration and temperature strongly affect this term.

Mass-Transport Overpotential#

At high current density, reactant/product transport limitations can be represented by:

\[ \eta_{\text{mt}} = -\frac{RT}{nF}\ln\!\left(1-\frac{j}{j_{\text{lim}}}\right) \]

where \(j_{\text{lim}}\) is an effective limiting current density. In PEM electrolysis, this term becomes more significant at high utilization, inadequate water feed, or poor gas removal.

Typical PEM Parameter Ranges#

Representative values used in system-level PEM models are:

  • Current density: \(1\) to \(3\,\mathrm{A\,cm^{-2}}\) (up to higher values in advanced stacks)

  • Cell voltage under load: about \(1.8\) to \(2.2\,\mathrm{V}\)

  • OER Tafel slope: typically \(40\) to \(80\,\mathrm{mV\,dec^{-1}}\)

  • HER Tafel slope: typically \(20\) to \(40\,\mathrm{mV\,dec^{-1}}\)

  • Area-specific resistance (ASR): often around \(0.1\) to \(0.3\,\Omega\,\mathrm{cm^2}\) (technology- and condition-dependent)

These ranges vary with catalyst loading, membrane thickness, temperature, pressure, and water management.

Transport Phenomena in PEM Cells#

Transport processes strongly couple electrochemistry, ohmic losses, and durability in PEM electrolyzers. The main coupled phenomena are water transport, proton transport, gas transport, and heat transport.

Water and Proton Transport in the Membrane#

Proton conduction in the membrane depends on hydration state and temperature. A common local relation is:

\[ i_{\mathrm{H^+}} = -\kappa_{\mathrm{m}}(\lambda, T)\,\nabla \phi_{\mathrm{m}} \]

where \(\kappa_{\mathrm{m}}\) is membrane proton conductivity, \(\lambda\) is membrane water content, and \(\phi_{\mathrm{m}}\) is ionic potential.

Membrane hydration evolves through competing mechanisms:

  • Electro-osmotic drag (water carried from anode to cathode with protons)

  • Back-diffusion (water diffusion from wetter to drier side)

  • Hydraulic permeation under pressure gradients

Insufficient hydration increases membrane resistance and accelerates degradation risk; excessive water accumulation can exacerbate gas transport limitations.

Reactant/Product Gas Transport#

At the anode, liquid water must reach catalyst sites while generated oxygen must be removed. At the cathode, hydrogen removal and water crossover management are critical.

Mass-transport limitations can be represented by an effective concentration overpotential, and they become important at high current density, high stoichiometric utilization, or inadequate channel/porous-layer design.

For modeling, gas-side transport is often represented with:

  • Channel-scale pressure-drop and species-balance equations

  • Effective diffusivity models in porous transport layers (PTLs)

  • Interfacial source terms linking electrochemical rates to local composition

Gas Crossover and Purity#

Hydrogen and oxygen can permeate through the membrane by diffusion and pressure-driven transport. Crossover increases with pressure differential and can affect:

  • Product gas purity

  • Faradaic efficiency at low current density

  • Safety margins due to combustible mixtures

Consequently, pressure control and membrane condition monitoring are central in high-pressure PEM operation.

Thermal Transport and Water Management#

Cell temperature is determined by Joule heating, activation losses, and heat removal through flow fields and plates. Temperature rise can improve kinetics and conductivity up to a point, but excessive thermal gradients may accelerate aging.

A practical control objective is to maintain a balanced operating window in which:

  • Membrane stays sufficiently hydrated

  • Reactant supply avoids starvation

  • Gas evacuation prevents local accumulation

  • Temperature remains within stack design limits

Faradaic Efficiency and Energy Efficiency#

Efficiency metrics quantify how effectively electrical input is converted into hydrogen.

Faradaic Efficiency#

Faradaic efficiency reflects the fraction of current that produces hydrogen:

\[ \eta_F = \frac{\dot n_{\mathrm{H_2,actual}}}{\dot n_{\mathrm{H_2,ideal}}} = \frac{\dot n_{\mathrm{H_2,actual}}}{I/(2F)} \]

In ideal operation \(\eta_F \approx 1\), but it can decrease due to gas crossover, parasitic side reactions, and measurement uncertainties. In PEM systems, faradaic efficiency is often high at nominal load and can drop at very low current density where crossover is relatively more important.

Energy Efficiency#

Voltage-based energy efficiency can be written as:

\[ \eta_E = \eta_F\,\frac{E_{\mathrm{rev}}}{U_{\mathrm{cell}}} \]

or, using thermodynamic quantities,

\[ \eta_E = \frac{\dot n_{\mathrm{H_2}}\,\Delta G}{P_{\mathrm{el}}} \]

where \(P_{\mathrm{el}} = U_{\mathrm{cell}} I\).

Because \(U_{\mathrm{cell}}\) includes all overpotentials, reducing activation, ohmic, and transport losses directly improves efficiency.

Influence of Operating Conditions#

  • Current density: increasing current usually increases production rate but also raises overpotentials, which lowers voltage efficiency.

  • Temperature: moderate temperature increase generally improves kinetics and membrane conductivity, reducing voltage losses; too high temperatures can affect durability.

  • Pressure: higher pressure reduces downstream compression demand and changes reversible voltage, but can increase crossover and impose stronger mechanical constraints.

  • Water feed and flow management: poor distribution can induce local dry-out or flooding, increasing both ohmic and transport losses.

In practice, optimal operation is a multi-objective trade-off between hydrogen production rate, electrical efficiency, gas purity, and long-term stack health.

Contents