In this paper, the corona discharge process of the bar-plate gap at −1 kV DC voltage is simulated using a two-dimensional axisymmetric plasma module. We analyze the variation of air negative corona discharge current, and the distribution morphology of microparticles in different discharge stages in detail. The significance of plasma chemical reactions at some typical time and the distribution characteristics of heavy particles are investigated according to reaction rates. Results show that, in the current rising stage, the collision ionization reactions (e.g., R1 and R2) and electron adsorption reaction (e.g., R3) play a major role, which lead to the increase in charged particles and the formation of an electron avalanche. In the current drop stage, all reaction rates decreased, except for collision ionization and electron attachment, partial charge transfer reactions (e.g., R8, R10, R11, and R14), and composite reactions (e.g., R16, R17, and R18), which come into play and gradually reduce the number of charged ions in the gap. In the current stabilizing stage, the main chemical reactions are composite reactions (e.g., R16 and R17), then the corona discharge ends. For the heavy particle distribution, O2+ and O4+ are the main positive ions, O2− is the most abundant negative ions, and the neutral particles are mainly O.