An adaptive grid refinement strategy for the simulation of negative streamers

The evolution of negative streamers during electric breakdown of a non-attaching gas can be described by a two-fluid model for electrons and positive ions. It consists of continuity equations for the charged particles including drift, diffusion and reaction in the local electric field, coupled to the Poisson equation for the electric potential. The model generates field enhancement and steep propagating ionization fronts at the tip of growing ionized filaments. An adaptive grid refinement method for the simulation of these structures is presented. It uses finite volume spatial discretizations and explicit time stepping, which allows the decoupling of the grids for the continuity equations from those for the Poisson equation. Standard refinement methods in which the refinement criterion is based on local error monitors fail due to the pulled character of the streamer front that propagates into a linearly unstable state. We present a refinement method which deals with all these features. Tests on one-dimensional streamer fronts as well as on three-dimensional streamers with cylindrical symmetry (hence effectively 2D for numerical purposes) are carried out successfully. Results on fine grids are presented, they show that such an adaptive grid method is needed to capture the streamer characteristics well. This refinement strategy enables us to adequately compute negative streamers in pure gases in the parameter regime where a physical instability appears: branching streamers.

[1]  Willem Hundsdorfer,et al.  Ebert and Hundsdorfer Reply , 2002 .

[2]  R. Nemzek,et al.  Television Image of a Large Upward Electrical Discharge Above a Thunderstorm System , 1990, Science.

[3]  Jens Lang,et al.  Adaptive Multilevel Solution of Nonlinear Parabolic PDE Systems - Theory, Algorithm, and Applications , 2001, Lecture Notes in Computational Science and Engineering.

[4]  Cornelis Vuik,et al.  How fast the Laplace equation was solved in 1995 , 1997 .

[5]  U. Inan,et al.  Electrical discharge from a thundercloud top to the lower ionosphere , 2002, Nature.

[6]  P. Williams,et al.  Two‐dimensional studies of streamers in gases , 1987 .

[7]  L. C. Pitchford,et al.  Electrohydrodynamic force and aerodynamic flow acceleration in surface dielectric barrier discharge , 2005 .

[8]  W. Hundsdorfer,et al.  Branching of negative streamers in free flight. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  A. Kulikovsky,et al.  RAPID COMMUNICATION: The role of the absorption length of photoionizing radiation in streamer dynamics in weak fields: a characteristic scale of ionization domain , 2000 .

[10]  Willem Hundsdorfer,et al.  Spontaneous branching of anode-directed streamers between planar electrodes. , 2001, Physical review letters.

[11]  Shengtai Li,et al.  Stability of Moving Mesh Systems of Partial Differential Equations , 1998, SIAM J. Sci. Comput..

[12]  D. Hampton,et al.  Preliminary results from the Sprites94 aircraft campaign: 1 , 1995 .

[13]  Kj Krzysztof Ptasinski,et al.  A high-temperature pulsed corona plasma system for fuel gas cleaning , 2004 .

[14]  B. H. Crichton,et al.  Gas discharge physics , 1996 .

[15]  Bardsley,et al.  Simulation of negative-streamer dynamics in nitrogen. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[16]  Roland A. Sweet,et al.  A direct method for the solution of poisson's equation with neumann boundary conditions on a staggered grid of arbitrary size , 1976 .

[17]  E. M. Veldhuizen,et al.  Electrical discharges for environmental purposes : fundamentals and applications , 2000 .

[18]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[19]  Lothar Schäfer,et al.  Regularization of moving boundaries in a laplacian field by a mixed Dirichlet-Neumann boundary condition: exact results. , 2005, Physical Review Letters.

[20]  W. Hundsdorfer,et al.  The multiscale nature of streamers , 2006, physics/0604023.

[21]  K. Schoenbach,et al.  Bacterial decontamination of water by means of pulsed corona discharges , 2001, PPPS-2001 Pulsed Power Plasma Science 2001. 28th IEEE International Conference on Plasma Science and 13th IEEE International Pulsed Power Conference. Digest of Papers (Cat. No.01CH37251).

[22]  M. Elta,et al.  Large‐signal time‐domain modeling of low‐pressure rf glow discharges , 1987 .

[23]  W. Saarloos,et al.  PROPAGATION AND STRUCTURE OF PLANAR STREAMER FRONTS , 1997, patt-sol/9702006.

[24]  Martin A. Gundersen,et al.  Premixed flame ignition by transient plasma discharges , 2002 .

[25]  S. Pancheshnyi,et al.  Two-dimensional numerical modelling of the cathode-directed streamer development in a long gap at high voltage , 2003 .

[26]  Golub Gene H. Et.Al.  Matrix Computations, 3rd Edition , 2007 .

[27]  Manuel Arrayás,et al.  Stability of negative ionization fronts: Regularization by electric screening? , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[28]  E. M. Veldhuizen,et al.  Streamer branching in a short gap: the influence of the power supply , 2002 .

[29]  W. C. Finney,et al.  Combined removal of SO/sub 2/, NO/sub x/, and fly ash from simulated flue gas using pulsed streamer corona , 1989 .

[30]  P. Colella,et al.  A Conservative Finite Difference Method for the Numerical Solution of Plasma Fluid Equations , 1999 .

[31]  W. Saarloos,et al.  Front propagation into unstable states : universal algebraic convergence towards uniformly translating pulled fronts , 2000, cond-mat/0003181.

[32]  Phillip Colella,et al.  Numerical Solution of Plasma Fluid Equations Using Locally Refined Grids , 1997 .

[33]  P. Wesseling.  Principles of Computational Fluid Dynamics , 2000 .

[34]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[35]  J. G. Verwer,et al.  A static-regridding method for two-dimensional parabolic partial differential equations , 1991 .

[36]  U. Ebert,et al.  Diffusion correction to the Raether–Meek criterion for the avalanche-to-streamer transition , 2006 .

[37]  Ningyu Liu,et al.  Effects of photoionization on similarity properties of streamers at various pressures in air , 2006 .

[38]  Gene H. Golub,et al.  Matrix computations , 1983 .

[39]  Yu. K. Bobrov,et al.  Application of high-resolution schemes in the modeling of ionization waves in gas discharges , 1998 .

[40]  J. Dutton,et al.  A survey of electron swarm data , 1975 .

[41]  Ningyu Liu,et al.  Effects of photoionization on propagation and branching of positive and negative streamers in sprites , 2004 .

[42]  A. Kulikovsky.  Comment on "Spontaneous branching of anode-directed streamers between planar electrodes". , 2002, Physical review letters.

[43]  U. Ebert,et al.  Streamer branching rationalized by conformal mapping techniques. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  S. Starikovskaia,et al.  Role of photoionization processes in propagation of cathode-directed streamer , 2001 .

[45]  Eduard M. Bazelyan,et al.  Lightning Physics and Lightning Protection , 2000 .

[46]  Jeroen Wackers,et al.  A nested-grid direct Poisson solver for concentrated source terms , 2005 .

[47]  Jan G. Verwer,et al.  Algorithm 758: VLUGR2: a vectorizable adaptive-grid solver for PDEs in 2D , 1996, TOMS.

[48]  U. Inan,et al.  Telescopic imaging of sprites , 2000 .

[49]  D. M. Sloan,et al.  Numerical Solution of Fisher's Equation Using a Moving Mesh Method , 1998 .

[50]  Jeroen Wackers.  A nested-grid finite-difference Poisson solver for concentrated source terms , 2003 .

[51]  J. Sid Clements,et al.  Formation of chemical species and their effects on microorganisms using a pulsed high voltage discharge in water , 1994, Proceedings of 1994 IEEE Industry Applications Society Annual Meeting.