Numerical and analytical studies of critical radius in new geometries for corona discharge in air and CO2-rich environments Jacob A. Engle, Jeremy A. Riousset Department of Physical Sciences, Center for Space and Atmospheric Research (CSAR), Embry-Riddle Aeronautical University, Daytona Beach, FL CEDAR 2017 ([email protected]) Abstract II. Model Formulation In this work, we focus on plasma discharge produced between two electrodes with a high potential Objectives Geometry Cartesian Spherical Cylindrical difference, resulting in ionization of the neutral gas particles and creating a current in the gas • Apply Paschen theory to Cartesian, 푥2 푅2 푅2 medium. This process, when done at low current and low temperature can create corona and “glow” Analytical 훼eff푑푥 = ln(푄) 훼eff푑푟 = ln(푄) 훼eff푑푟 = ln(푄) 푥1 푅1 푅1 discharges, which can be observed as a luminescent, or “glow,” emission. The parallel plate geometry spherical, and cylindrical geometries; 푥1 = 0 R2 →c; V(R2) = 0 R2 →c; V(R2) = 0 −퐵푝 used in Paschen theory is particularly well suited to model experimental laboratory scenario. V(c) = 0; V(c) = 0 • Obtain analytical expressions for critical 훼eff(퐸) = 퐴푝푒 퐸 However, it is limited in its applicability to lightning rods and power lines (Moore et al., 2000). −퐵푝 −퐵푝 푑 = 푥 − 푥 훼 (퐸) = 퐴푝푒 퐸 훼 (퐸) = 퐴푝푒 퐸 Franklin’s sharp tip and Moore et al.’s rounded tip fundamentally differ in the radius of curvature of 2 1 eff eff radius and Stoletov’s point; 휕푉 = 0: Stoletov′s point 휕푉 ′ 휕푉 ′ the upper end of the rod. Hence, we propose to expand the classic Cartesian geometry into spherical 휕푑 = 0: Stoletov s point = 0: Stoletov s point • Develop numerical models for 휕푅1 휕푅1 geometries. In a spherical case, a small radius effectively represents a sharp tip rod, while larger, 푐2 푐 퐸 푅1 = 퐸푐 2 퐸 푅1 = 퐸푐 centimeter-scale radius represents a rounded, or blunted tip. Experimental investigations of lightning- Cartesian, spherical, and cylindrical 푅1 푅1 like discharge are limited in size. They are typically either a few meters in height, or span along the geometries; ground to allow the discharge to develop over a large distance. Yet, neither scenarios account for the 푥2 푅2 푅2 Numerical 훼eff푑푥 = ln(푄) 훼eff푑푟 = ln(푄) 훼eff푑푟 = ln(푄) change in pressure, which conditions the reduced electric field, and therefore hardly reproduce the • Verify numerical models and analytical 푥1 푅1 푅1 푥 = 0 R →c; V(R ) = 0 R →c; V(R ) = 0 condition of discharge as it would occur under normal atmospheric conditions (Gibson et al, 2009). In 1 2 2 2 2 solutions with experimental data; νi 퐸 −νa(퐸) νi 퐸 −νa(퐸) νi 퐸 −νa(퐸) 훼eff(퐸) = 훼eff(퐸) = 훼eff(퐸) = this work we explore the effects of shifting from the classical parallel plate analysis to spherical and 휇e 퐸 퐸 휇e 퐸 퐸 휇e 퐸 퐸 • Establish the differences between sharp 휕푉 휕푉 cylindrical geometries more adapted for studies of lightning rods and power transmission lines, 푑 = 푥2 − 푥1 = 0: Stoletov′s point = 0: Stoletov′s point respectively. Utilizing Townsend’s equation for corona discharge, we estimate a critical radius and 휕푉 휕푅1 휕푅1 and blunt tipped rods for corona = 0: Stoletov′s point minimum breakdown voltage that allows ionization of neutral gas and formation of a glow corona 휕푑 Boltzmann solver Boltzmann solver 2 푐 Boltzmann solver 푐 퐸 푅1 = 퐸푐 around an electrode in air. Additionally, we explore the influence of the gas in which the discharge discharges 퐸 푅 = 퐸 푅1 1 푐 푅 2 develops. We use Bolsig, a numerical solver for the Boltzmann equation, to calculate Townsend • Generalize to any atmosphere using a 1 • 푝 = 푁푘 푇 coefficients for CO2-rich atmospheric conditions (Hagelaar and Pitchford, 2005). This allows us to 퐵 Assumptions: 푁 explore the feasibility of a glow corona on other planetary bodies such as Mars. We calculate the Boltzmann solver (Hagelaar and • 퐸 푅 = 퐸 푐 = 퐸 ≈ 30 0 kV/cm (Earth) 1 푐 푁 breakdown criterion both numerically and analytically to present simplified formulae per each Figure 5: Fit of the exponential approximation for 훼eff(E) for coefficients obtained Pitchford, 2005). • ∇. 푬 = 휌0 = 0 geometry and gas mixture. from: Morrow and Lowke (1997), Hagelaar and Pitchford (2005). I. Introduction III. Results Coefficients and Stoletov’s points Cartesian solutions • A and B coefficients derived from the exponential fit accurately predict the Corona Discharge −퐵푝 minimum voltages (Table 2); • Electrical discharge around a • Critical electric field: 퐸 d = ln 푄 ln • Numerical, analytical, and experimental data are all in excellent agreement conductor due to electric field; 퐴푝푑 in the recreated Cartesian solution; • Weakly ionized gas responsible for • Minimum breakdown voltage: • CO2 dominated atmospheres have a higher critical electric field than air at glow at visible wavelengths; comparable densities; • Hypothesized to promote the −퐵푝푑 • Mars minimum breakdown voltages are lower than Earth due to low formation of upward connecting 푉 푑 = ln 푄 Martian atmospheric pressure (0.6% P ). ln Earth leaders in lightning discharges. 퐴푝푑 Coefficients Raizer Morrow and Bolsig+ Bolsig+ (Mars) Figure 1: Glow Coronas form on the edges of a 푒퐵 (1991) Lowke (1997) (Earth) powerline transformer (Berkoff, 2005). • Stoletov’s point: 푉min = ln 푄 퐴 (Earth) (Earth) Electron Avalanche • CO and air solutions taken at STP A (1/cm/Torr) 15 7.7 9.29 33.44 The process of electron avalanching 2 is similar between various types of • Boltzmann equation solver (Bolsig) B (V/cm/Torr) 365 274.7 295.18 430.07 discharges: • Table 2: Exponential approximation coefficients (A and B) from Figure 5 found from fitting: • Initial step of a discharge; Comparison with experimental data −퐵푝 • 훼eff(퐸) = 퐴푝푒 퐸 Release of secondary electrons in • Convergence of solutions near Stoletov’s points electron-neutral collision; Stoletov’s Point: 푽퐦퐢퐧 (V) Analytical Numerical • Figure 6 → : Critical electric field and breakdown voltage to meet Secondary electrons with enough the initiation criteria. The critical voltage curves are plotted at STP KE to repeat the process; to be consistent with the conditions of the experimental data. Cartesian (Air) 348.2 350.9 • Avalanche criteria: (Raiser, 1991) Cartesian (CO ) 517.6 603 푅2 4 Figure 2: Visual representation of the process of an 2 훼eff 푑푟= ln(Q) ≈ 18-20; Q = 10 Spherical solutions 푅1 electron avalanche in Townsend’s breakdown model. This 4퐵(ln 푄 +퐴푝푟)2 Spherical (Earth) 1414 1709 can also be referred to as a Cartesian case (Gewartowski et • Critical electric field: 퐸 푟 = 2 2 Types of Discharges al., 1965). 휋푝퐴 푟 Spherical (Mars) 475.4 603.1 Cylindrical (Earth) 1426 1132 Parameter Glow Corona Streamer Leader • Minimum breakdown voltage: Temperature ∼300 K ∼300 K ≳5000 K Cylindrical (Mars) 584.3 469.8 4퐵(ln 푄 +퐴푝푟)2 Electron energy 1-2 eV 5-15 eV 1-2 eV 푉 푟 = Table 3: The minimum breakdown voltages for each geometry and atmosphere; also known as 휋푝퐴2푟 Stoletov’s points. Electric field 0.2-2.7 kV/cm 5-7.5 kV/cm 1-5 kV/cm 16퐵 • Stoletov’s point: 푉 = ln 푄 . Electron density 2.6×108 cm-3 5×1013-1015 cm-3 4×1014 cm-3 min 휋퐴 IV. CONCLUSIONS Table 1: Characteristics for types of discharge at sea level in Earth’s atmosphere, adapted from (Gibson et • Largest error due to Taylor expansion of The results and conclusions obtained in this work can be summarized as al, 2009). Gauss error function follows: • A new model for calculations of the critical radius and minimum breakdown voltage • Boltzmann equation solver (Bolsig) for Corona discharge in Cartesian, spherical, and cylindrical geometries is presented; • Highest minimum breakdown voltage • The model is validated using classic Paschen theory and experimental data in air from Meek and Craggs (1978) and CO2 from Stumbo (2013); Figure 7 → : Critical electric field and breakdown voltage to meet • We expand classic Paschen theory into an analytical solution for spherical and the initiation criteria. The critical voltage curves are plotted at STP cylindrical geometry; for each planetary body (Mars and Earth). • Our numerical model and the analytical solution show excellent agreement with (A) (B) (C) experimental data; • The significantly lower pressure on Mars compared to Earth lowers the minimum Figure 3: (A) A Wartenberg wheel with glow coronas forming at the tip of each spindle (Berkoff, 2005); (B) Cylindrical solutions breakdown voltage required to create corona discharge. Streamers forming a sprite phenomenon (courtesy of H. H. C. Stenbaek-Nielsen); (C) Lightning channels as an example of leader discharge (Whetmore, 2016). 퐵(ln 푄 +퐴푝푟) Acknowledgements: • Critical electric field: 퐸 푟 = −퐵푝 This work is supported by the Embry-Riddle Aeronautical Application to Martian Studies 퐴푟(1−푒 퐸0 ) University Office of Undergraduate Research (ERAU) and Motivations: Earth Analogy: the Center for Space and Atmospheric Research (CSAR). • Minimum breakdown voltage: • Potential hazard due to arcing on • Tribocharging in Martian dust storms REFERENCES landers and rovers; akin to Earth sandstorms; • E. Berkoff, The Corona Discharge: Its Properties and Uses - Colorado Wire & Cable (2015). 퐵ln 푄 • W. M. Farrell and M. D. Desch. Is there a Martian atmospheric electric circuit? J. Geophys. Res., • Interference with sensitive external • Charge separation due to ln( ) 퐴퐸0푟 106:7591–7596, 4 (2001). doi: 10.1029/2000JE001271. systems and data measurements; sedimentation & gravitation; 푉 푟 = 퐵푝푟 퐴푝푟 • J. Gewartoski, J. W. Watson, H. Alexander. Principles of Electron Tubes: Including Grid-controlled Tubes, 푊( +1) • Possibility of electrical shortage and • Integration in the Martian global ln 푄 Microwave Tubes and Gas Tubes, D. Van Nostrand Co., Inc. (1965). • A. S. Gibson, J. A. Riousset, and V. P. Pasko. Minimum breakdown voltages for corona discharge in failure.