In the present work a new lattice Boltzmann (LB) framework has been developed to study the electroosmotic flows in a 2-D flat microchannel. The governing equations are presented in the continuum model, while a set of equivalent equations in LB model is introduced and solved numerically. In particular, the Poisson and the Nernst–Planck (NP) equations are solved by two new lattice evolution methods. In the analysis of electroosmotic flows, when the convective effects are not negligible or the Electric Double Layers (EDLs) have overlap, the NP equations must be employed to determine the ionic distribution throughout the microchannel. The results of these new models have been validated by available analytical and numerical results. The new framework has also been used to examine the electroosmotic flows in single and parallel heterogeneous microchannels.