A robust numerical framework for positive streamer modelling based on electro-hydrodynamic equations coupled with Poisson and Helmholtz differential equations for the photoionization process is presented. The proposed multi-layer meshing scheme in a 2D non-axisymmetric finite-element model along with a hybrid meshing technique presented in part I of this series paper for negative streamers provide high accuracy, spatial resolution, and capability to present the major features of both positive and negative streamers. In addition, the presented model is utilized to simulate multi positive and negative streamers propagation in a non-uniform electric field in the air. The main characteristics of the positive and negative streamers including the morphology, distribution pattern of space charges, local electric field, diameter, length, and velocity are presented, discussed, and compared with the experiment. Moreover, the impacts of initial seed density and voltage on the propagation of streamers are presented and explained. The branching mechanism arising from Laplacian instability and its impact on the streamer parameters such as tip electric field and dominant charge density is explained.