Direct numerical simulations of electroconvection instability near an ion-selective surface are conducted using a mesoscopic lattice Boltzmann method (LBM). An electrohydrodynamic model of ion transport and fluid flow is presented. We numerically solve the Poisson-Nernst-Planck equations for the electric field and the Navier-Stokes equations for the flow field. The results cover Ohmic, limiting, and overlimiting current regimes, and they are in good agreement with the asymptotic analytical solution for the relationship between current and voltage. The influences of different ion transport mechanisms on the voltage-current relationship are discussed. The results reveal that the electroconvection mechanism is as important as other ion transport mechanisms in electrohydrodynamic flow. By comparing the contribution of different regions in the numerical domain, we find that the flow in the extended space charge layer is dominated by electroconvection. We also study the influences of multiple driving parameters, and the electrohydrodynamic coupling constant plays a dominant role in triggering convective instability. The flow pattern and ion concentration distribution are described in detail. Moreover, the route of flow from steady state to periodic oscillation and then to chaos is discussed.