Gas-liquid turbulent two phase flow along a vertical flat plate with gas evolution was numerically simulated with a two-fluid model and Plandtl’s mixing length theory using the SIMPLER method, The calculation domaint was assumed to be two-dimensional and gas bubbles to he evolved uniformly from the plate surface, The computed velocity profiles were affected by void fraction on the wall and the bubble diameters, but;not by the bubble depature velocity from the plate surface. Both the gas and liquid velocity profiles and the computed wall shear stress with the present model were qualitatively and quantitatively coincident with those measured using LDV for the vertical flat plate-electrode with H-2 or O-2 gas evolution.