Corona discharge is a significant problem in the operation of high-voltage transmission and distribution systems, particularly for polymer insulators. Numerical simulation has become an effective tool for investigating the underlying physical mechanisms and optimizing the design of insulators. In this paper, we present a two-dimensional numerical simulation study on corona discharge plasma affecting the surface behavior of polymer insulators. The simulation was performed with the Comsol Multiphysic software and is based on the finite element method and the fluid plasma model, which considers ionization, recombination, and the transport of plasma species. The numerical results are analyzed to study the spatial and temporal characteristics of the corona discharge and its effect on the surface behavior of polymer insulators. The results show that the electric field is affected not only by the volume charge density but also by the surface charge density, which in turn depends on the densities of the charge carriers migrating on the insulator surface. However, the electric field drops drastically when one or two grading rings are installed. But one grading ring is not enough to limit the discharge.