“…Firstly, to reduce the tedious computation associated with the 2D continuity equations, a fully implicit discretization technique is adopted to discard the restriction of the CFL condition; furthermore, variable time step technique is taken from previous work. [10,11] Secondly, the efficient photoionization (important for positive corona) model [15,16] is first adopted to restore the physical process of corona current pulse production, leading to just a 2.9% increase of computation-cost for the simulation of a 10-cm rod-plane gap; whereas, in previous 2D simulations, photoionization was usually ignored [5,7,10,11] or replaced by a uniform background photoionization [8,17] for its high computation complexity, which was evidenced as "the consideration of photoionization costs nearly 90 percent of the computing time in the whole simulation." [17] The rest of the present paper is organized as follows: firstly, a detailed description of the methodology is presented, which includes the governing equations, boundary conditions, initial conditions, discretization, etc; secondly, the method is applied to a rod-plane example, and a comparison between the simulated results and the measured results is made to verify the validity of the proposed method; thirdly, the convergence of the present method and the physical implications of the results are discussed; finally, some conclusions are derived.…”