Numerical and analytical studies of critical radius in spherical and cylindrical geometries for corona discharge in air and CO2-rich environments

In order to determine the most effective geometry of a lightning rod, one must first understand the physical difference between their current designs. Benjamin Franklin's original theory of sharp tipped rods suggests an increase of local electric field, while Moore et al.'s (2000) studies of rounded tips evince an increased probability of strike (Moore et al., 2000; Gibson et al., 2009).In this analysis, the plasma discharge is produced between two electrodes with a high potential difference, resulting in ionization of the neutral gas particle. This process, when done at low current and low temperature can create a corona discharges, which can be observed as a luminescent emission. The Cartesian geometry known as Paschen, or Townsend, theory is particularly well suited to model experimental laboratory scenario, however, it is limited in its applicability to lightning rods. Franklin's sharp tip and Moore et al.'s (2000) rounded tip fundamentally differ in the radius of curvature of the upper end of the rod. As a first approximation, the rod can be modelled as an equipotential conducting sphere above the ground. Hence, we expand the classic Cartesian geometry into spherical and cylindrical geometries. In this work we explore the effects of shifting from the classical parallel plate analysis to spherical and cylindrical geometries more adapted for studies of lightning rods or power lines. Utilizing Townsend's equation for corona discharge, we estimate a critical radius and minimum breakdown voltage that allows ionization of the air around an electrode. Additionally, we explore the influence of the gas in which the discharge develops. We use BOLSIG+, a numerical solver for the Boltzmann equation, to calculate Townsend coefficients for CO2-rich atmospheric conditions. This allows us to expand the scope of this study to other planetary bodies such as Mars (Hagelaar, 2005). We solve the problem both numerically and analytically to present simplified formulas per each geometry and gas mixture. The development of a numerical framework will ultimately let us test the influence of additional parameters such as background ionization, initiation criterion, and charge conservation on the values of the critical radius and minimum breakdown voltage.