Numerical Simulation of the Two-Phase Flow Dynamics in an Alkaline Water Electrolyzer Considering Difference in Bubble Diameter

Abstract

For more efficient hydrogen production, predicting the characteristics of bubble flow inside a water electrolyzer is critical since bubbles affect ion transport, and electrode coverage by bubbles influence overpotentials and the entire electrolysis performance. Transport of bubble significantly depends on its size. In this study, a CFD model of an alkaline water electrolyzer which considers difference in sizes of generated bubbles was developed. The model treats bubbles of different sizes as different gas phases, which are independent of each other. The total bubble mass source for a given current density is divided into phases of each bubble size according to experimental visualizations. Drag, virtual mass, lift, wall lubrication, and turbulent dispersion forces are considered as interaction between gas and liquid phase.The following three sizes of bubbles were considered: diameters d = 36 μm (small), 75 μm (medium), and 208 μm (large). The entire calculation domain represents the anode side of a zero-gap monopolar alkaline water electrolyzer. Figure 1 shows the volume fraction fields of the three gas phases (i.e., bubble sizes) at a constant current density of 0.6 A/cm2. The bubbles are generated on the electrode surface at the z=0 plane. As the liquid phase forms a circulating flow inside the electrolyzer, small bubbles are transported along the circulating downstream flow due to smaller buoyancy force compared to their viscous drag. In contrast, large bubbles flow out directly to the upper outlet because buoyancy force dominates. This result points out the need to consider the difference in bubble diameter for designing flow field inside alkaline water electrolyzer.[Acknowledgements]This study was based on results obtained from the Development of Fundamental Technology for Advancement of Water Electrolysis Hydrogen Production in Advancement of Hydrogen Technologies and Utilization Project (P14021) commissioned by the New Energy and Industrial Technology Development Organization (NEDO).Fig. 1. Total gas phase volume fraction field obtained from the sum of the volume fraction field for each size bubble, and its experimental visualization. Figure 1