a b s t r a c tThis work is devoted to the modeling of two phase flows arising in typical electrolysis devices. A numerical mixture model is used in order to resolve the two dimensional bubble plumes evolving along the electrodes. Plumes thickness sensitivity is studied for various parameters, such as bubble diameter, electrolyte viscosity, electrochemical cell geometry and current density. Using thermal buoyancy driven flow analogy, a dimensionless Rayleigh-like number Ra f ;e is defined to predict the behavior of the wallbounded gas convection between two vertical facing electrodes. Different bubbles dispersion mechanisms are observed depending on two-phase flow dynamics and physical properties of the mixture. The effect of forced convection in the channel is also investigated. A scaling law for plume thickness evolution for a large range of Prandtl-equivalent number values is proposed. These results show that the bubble plume can be efficiently controlled by an imposed electrolyte velocity.

International audienceDirect Numerical Simulations of a laminar two-phase flow into a vertical channel are investigated. An Eulerian–Lagrangian approach allows tracking each bubble position with a four-way coupling strategy, i.e. taking into account bubble-fluid and bubble-bubble interactions. The flow configuration has been chosen to highlight the buoyancy effects due to significant values of void fraction (high numbers of bubbles); hence the bubbles collisions and wall effects are the critical parameters to ensure the dispersion of the bubble plume. The DNS approach is self-consistent and does not rely on closure relations or empirical correlations for describing the collective bubble dynamics. It is found that the DNS predicts well the behavior of the bubble plume and its back effect on the liquid phase when compared with a mixture model and experimental data. The elastic nature of collisions, the sensitivity of the mean and RMS values of velocities and void fraction to the mesh quality are explored