We study a microfluidic flow model where the movement of several charged species is coupled with electric field and the motion of ambient fluid. The main numerical difficulty in this model is the net charge neutrality assumption which makes the system essentially overdetermined. Hence we propose to use the involutive and the associated augmented form of the system in numerical computations. Numerical experiments on electrophoresis and stacking show that the completed system significantly improves electroneutrality constraint conservation and recovers analytical results while a direct implementation of the initial model fails.
