 A 2D axisymmetric electrochemical model of a lithium experimental cell is solved using a finite element method. The model is considering the coupled effect of momentum, electric, kinetic and mass transfer phenomena. The dense diaphragm used in the setup, which is not hydraulically permeable, separates the cell into two regions: 1 (a) turbulent region, between the anode and the diaphragm, and 2 (a) laminar region between the diaphragm and the cathode. The k-epsilon model is used to solve the turbulent flow resulting from bubbles generation at the anode. A two-phase flow model is also developed to simulate the volume fractions of bubbles and electrolyte in the cell. The non-uniform bubbles distribution over the anode surface, derived from the non-uniform current distribution, has been added to the two-phase model. The effects of the diaphragm length, position and porosity on the electric and two-phase flow fields are simulated. In fact, the diaphragm position and length both influence the current distribution at the surface of electrodes and the velocity distribution in the cell, all of which influence ohmic and kinetic overpotentials. The results show that up to 40% of energy can be saved when running the lithium electrolysis cell with a shorter porous diaphragm located as far as possible from the anode. The maximum current density, found at the bottom corner of anode, is higher when the diaphragm is longer and when it is closer to the anode.