Modeling Streamer Discharge in Air Using Implicit and Explicit Finite Difference Methods with Flux Correction

Implementing a computationally efficient numerical model for a single streamer discharge is essential to understand the complex processes such as lightning initiation and electrical discharges in high voltage systems. In this paper, we present a streamer discharge simulation in air, by solving one-dimensional (1D) drift diffusion reaction (DDR) equations for charged species with the disc approximation for electric field. A recently developed fourth-order space and time-centered implicit finite difference method (FDM) with a flux-corrected transport (FCT) method is applied to solve the DDR equations, followed by a comparative simulation using the well-established explicit FDM with FCT. The results demonstrate good agreement between implicit and explicit FDMs, verifying their reliability for streamer modeling. The total electrons, total charge, streamer position, and hence the streamer bridging time obtained using the FDMs with FCT agree with the same streamer computed in the literature using different numerical methods and dimensions. The electric field is obtained with good accuracy due to the inclusion of image charges representing the electrodes in the disc method. This accuracy can be further improved by introducing more image charges. Both implicit and explicit FDMs effectively capture the key streamer behavior, including the variations in charged particle densities and electric field. However, the implicit FDM is computationally more efficient.