What electrode hole size to use in zero gap alkaline water electrolysers?

The electrodes in modern electrolysers for the production of green hydrogen have holes to let bubbles escape, but what size should these holes be? Despite more than half a century since its first usage, there is no clear guidance on this question available in the open literature.

Therefore, we performed an extensive experimental campaign under industrially relevant conditions with 17 electrodes, with varying hole sizes and shapes. We found that by far the most important feature determining the electrode performance was the hole size, not its shape, or the electrode thickness or other features like the presence of pillars.

We found that for most relevant current densities and atmospheric pressure a hole size of several millimeter is optimal. Much larger holes lead to larger cell voltages as the current becomes more inhomogeneously distributed, increasing the resistance. But particularly much smaller holes lead to excessive overpotentials.

To study why this is, we looked at the behavior of gas bubbles, as they leave the holes from the back, but also from the front, through a transparent membrane. It turned out that if the holes become so small that many times per second bubbles are generated that can fill these holes, the chance of them clogging these holes becomes large. Holes substantially smaller than a millimetre tend to clog very frequently. This clogging leads to large gas films forming between the electrode and the diaphragm, inactivating large parts of the electrode surface, forcing the reaction to the backside, increasing the resistance.

Millimeter-sized holes are largely immune to this clogging and perform much better, which explains why commercial expanded metal electrodes have holes in this size range.

J.W. Haverkort, A.S. Aghdam, E.J.B. Craye. The optimal electrode hole size in zero gap alkaline water electrolysis: A combined electrochemical, theoretical, and bubble imaging approach. Int. J. Hydrogen Energy (2025), 171, 150919.

What keeps the gas fraction from approaching 100%?

Of course, the volume fraction of hydrogen gas bubbles cannot exceed 100%. Consider spherical gas bubbles of equal size, the maximum theoretical packing fraction is much lower at about 74 %, or even roughly 64 % for a random packing. The space between bubbles can be taken up by smaller bubbles and the space between those smaller bubbles again by smaller bubbles, so that theoretically there is no strict limit below 100%. Also, gas bubbles could coalesce to form a big bubble filling the entire electrolyser space so that the gas fraction could approach 100%.

However, it turns out that bubbles in strong electrolytes do not coalesce so easily and prefer to stay at a small distance from each other. This makes that these hydrogen and oxygen bubbles behave somewhat like solid particles. Therefore, we introduced into our models what in granular matter modelling is called solid pressure. As the gas fraction increases, a repulsive pressure avoids a further increase.

Figure: Simulations of the gas fraction near a zero-gap electrode. Red is about 60%.

We find that a maximum gas fraction around 65 % corresponds well with experimental results on the bubble-induced resistance. In this way, the extra ingredient of a solid pressure helps simulations correspond more to reality and as a result converge more easily as extremely high gas fractions are avoided.

W.L. van der Does, N. Valle, J.W. Haverkort, Multiphase alkaline water electrolysis simulations: The need for a solid pressure model to explain experimental bubble overpotentials, Int. J. Hydrogen Energy (2025), 102, p 295-303

Can less catalyst be more effective?

CO2 and water can be electrochemically converted to the more valuable and industrially useful mixture of CO and H2 (syngas). However, as more CO2 is converted relatively more H2 is produced. Therofore, a tempting strategy is to use less catalyst downstream and produce relatively more of the desireable CO at lower catalyst costs.

We first made a simple model that we optimised for our self-defined “effectiveness” giving

On the left is the average CO2 reaction rate coefficient <k> times the residence time L/U. Here q is the ratio between the inlet CO2 and H2 reaction rate constants.

By using less catalyst downstream, <k> decreases and a smaller residence time is better. However...the same effectiveness is obtained!

A similar conclusion is obtained for a more elaborate model including more physical and economical complexity. For a given residence time the catalyst distribution can be optimised using our freely available Matlab script. But in this case the improvements in effectiveness can be larger than those obtained by simply increasing the residence time, while at the same time using less catalyst!

Blake, J. W., Haverkort, J. W., & Padding, J. T. (2024). Less is more: Optimisation of variable catalyst loading in CO2 electroreduction. Electrochimica Acta, 145177.

What determines mass transfer at gas evolving electrodes?

Mass transfer at electrodes with gaseous products is excellent, because bubbles that coalesce or leave the electrode mix the fluid very well. Additionally, rising bubbles set the electrolyte into motion causing additional flow transport. Traditionally, these two contributions, kμ and kf, have been simply added. To improve upon this oversimplification, a new simple ‘addition rule’ was devised

which takes into account the surface coverage θ of bubbles and a single empirical micromixing parameter a.

Next, this formula was compared with all available literature data for water electrolysis, revealing large differences. For hydrogen in alkaline electrolytes, bubbles hardly mix and the mass transfer is determined by natural convection. For hydrogen bubbles in sulphuric acid, on the other hand, natural convection can be usually neglected and for other cases, both mechanisms can be important.

The reason why there are these large differences in the parameter a is hypothesized to be due to solutal Marangoni convection, which repels hydrogen bubbles away from each other and the electrode in alkaline electrolytes while it attracts hydrogen bubbles towards each other and the electrodes in sulphuric acid.

Haverkort, J. W. (2024). A general mass transfer equation for gas-evolving electrodes Int. J. Hydrogen Energy74, 283-296.

Electrolysers, Fuel Cells and Batteries: Analytical Modelling

Electrochemical engineering deals with electrochemical devices like electrolysers, fuel cells, and batteries. While several excellent books exist in this long-standing and still growing field, their focus is usually on chemistry or  phenomenology. In this textbook, we focus on mathematical modelling of the physical phenomena involved. Instead of resorting to numerical modelling, the aim is to derive simplified analytical models that maximise understanding.

Porous electrodes, ion mass transport, and multiphase flow are central themes in this book. Examples include modelling the water saturation in a fuel cell diffusion layer, the gas fraction and current distribution in an alkaline water electrolyser, the potential distribution in a binary electrolyte inside porous battery electrode, and the concentration distribution in the flow channel of a redox flow battery. This makes for a diverse, challenging, and stimulating journey, for both students and researchers.

Order a hard-copy
Read online
Download a free pdf

Two-dimensional, transient, yet analytical model for capacitive deionization

A promising way to remove low salt concentrations from brackish water is capacitive deionization (CDI). Instead of removing the water from the salt, as for example in reverse osmosis, in CDI the salt ions are removed from the water by applying an electric field. The ions are removed from the main flow and stored in the electric double layers of the porous electrodes flanking the main flow channel. Despite the fact that this is an inherently two-dimensional and transient process, usually modelled with several partial differential equations, we managed to simplify the problem to two coupled ordinary differential equations and obtained an explicit analytical solution. Here is a comparison of the salt concentration from our analytical model (left) with that from a comprehensive computational COMSOL model (right):

The power of such an analytical solution lies in making optimization much easier. Instead of seven individual physical and geometrical parameters, time, and a spatial coordinate, our solution primarily depends on a single dimensionless number that is a combination of these parameters.

We find that the optimal porous electrodes are roughly six times thinner than the spacer. By minimizing the energy losses and maximizing the amount of water processed, we find that an optimal design can increase the latter metric by an order of magnitude compared to typical values in the literature.

J.W. Haverkort, B. Sanderse, J.T. Padding, and J.W. Blake (2024). An analytical
ow-by capacitive deionization mode. Desalination 582, 117408.

What gap and pore size maximize the energy efficiency of electrochemical flow cells?

Various electrochemical cells including microfluidic fuel cells, membraneless redox flow batteries, and microfluidic (CO2) electrolysers include a channel between their two electrodes. The wider this channel, the larger the ohmic drop. The thinner this channel, the larger the pressure drop. We obtained a simple analytical formula for the optimal channel width that minimizes the combined associated energy dissipation, which we verified using numerical simulations.

To increase surface area while avoiding diffusion limitations, sometimes porous flow-through electrodes are used. Here, a similar optimization can be performed. Smaller pore sizes give more reactive surface area, but also a larger pressure drop. Also here an analytical formula was obtained, which works well even for the popular interdigitated flow configuration.

We compared our formulas with values from various papers in the literature and found that typically an order of magnitude too large channels and pores are used. Therefore, significant improvements in energy efficiency can be obtained by further miniaturization.

The flow velocity (arrows) and reactant concentration (color) inside a repeating unit of an interdigitated flow field (left). The combined activation and pumping dissipation from simulations (diamonds) and the obtained analytical expression (solid line) along with the optimal (dashed vertical line) volumetric surface area a (roughly the inverse of pore size).

Bhadra, A., & Haverkort, J. W. (2023). The optimal electrode pore size and channel width in electrochemical flow cells. Journal of Power Sources579, 233240.

How thick is the bubble layer near an electrode ?

Hydrogen bubbles generated at an electrode of a water electrolyzer rise due to their buoyancy and set the surrounding liquid into motion. Alternatively, the bubbles make the mixture lighter causing it to rise – similar to the hot fluid near a heated plate. An important difference with this well-studied ‘thermal natural convection’ case is that the gas fraction has a certain maximum, below one, which influences the plume in a previously unknown way. This understanding is important, for example, to describe heat and mass transport. Knowing the bubble layer thickness is also useful to properly choose the dimensions of an electrolyzer.

Therefore, we made an analytical model to describe the plume thickness and associated velocity profile. To validate this model, we performed computational simulations. A comparison between the analytical and computational models is shown in the below figure. The resulting formulas will be very useful for scaling up electrolyzers and understanding the effect of height and current density on gas hold-up and heat transport.

Reasonable agreement is obtained between the velocity profiles of our new analytical model (solid lines) and those of a computational mixture model (arrows).

Rajora, A., & Haverkort, J. W. (2023). An analytical model for the velocity and gas fraction profiles near gas-evolving electrodes. International Journal of Hydrogen Energy.

What would be the performance of a 1 m tall CO2 electrolyzer?

CO2-electrolysis is a promising way of making carbon-based fuels and chemicals from CO2, instead of from fossil fuels. Affordable electrolyzers that can efficiently convert CO2 and water into for example syngas (CO and H2) hold great promise for transforming the way carbon-based products can be industrially made in the future. There has been enormous research progress over the last few years, typically in lab set-ups with dimensions of a few centimeter. In a recent publication we tried to answer the question what would happen in a much taller electrolyzer. We find that much more of the CO2 entering the electrolyser is converted and also the fraction of the desirable CO in the outlet stream is increased. Unfortunately, our simulations show that this seemingly positive effect arises because more CO2 undergoes undesirable side-reactions with the electrolyte. As a partial remediation, we show that by varying the amount of catalyst along the electrolyzer this problem can be substantially reduced.

Joseph W. Blake, Vojtecȟ Konderla, Lorenz M. Baumgartner, David A. Vermaas, Johan T. Padding, and J. W. Haverkort, Inhomogeneities in the Catholyte Channel Limit the Upscaling of CO2 Flow Electrolysers, ACS Sustainable Chem. Eng. 2023, 11, 7, 2840–2852

What is the maximum height of a membraneless electrolyser?

Water electrolyzers without a membrane have the potential to make green hydrogen more energy efficiently and cheaper. Perhaps the simplest type of membraneless electrolyzer consists of two vertical parallel plate electrodes with upwards electrolyte flow separating hydrogen and oxygen bubbles, avoiding the formation of an explosive mixture. The faster the flow, the thinner the bubble plumes and the closer the electrodes can be placed together, resulting in a higher energy efficiency. Natural convection can only provide modest velocities. Forced flow can provide higher velocities, but to avoid turbulent mixing of the bubble plumes, these are limited to similarly modest laminar flow velocities.

A comparison of the electrolyte velocity obtained numerically (arrows) with our analytical model (solid line)

To determine how tall a membraneless electrolyzer can be made while still avoiding overlap between the oxygen and hydrogen gas plumes, we developed an analytical model that we verified with a more complete computational model and validated with experimental data from the literature. Based on our model we show that natural convection can allow safe and efficient atmospheric membraneless electrolysers up to about 5-10 cm height, while forced flow adds another 10 cm. At higher pressure, or by inducing smaller bubbles, taller or more energy efficient electrolysers of this type can be made.

Rajora, A., & Haverkort, J. W. (2022). An Analytical Multiphase Flow Model for Parallel Plate Electrolyzers. Chemical Engineering Science, 117823.

See this previous post for an alternative type of membraneless electrolyzer.