High order fluid model for streamer discharges. II. Numerical solution and investigation of planar fronts

The high order fluid model developed in the preceding paper is employed here to study the propagation of negative planar streamer fronts in pure nitrogen. The model consists of the balance equations for electron density, average electron velocity, average electron energy and average electron energy flux. These balance equations have been obtained as velocity moments of Boltzmann's equation and are here coupled to the Poisson equation for the space charge electric field. Here the results of simulations with the high order model, with a PIC/MC (Particle in cell/Monte Carlo) model and with the first order fluid model based on the hydrodynamic drift-diffusion approximation are presented and compared. The comparison with the MC model clearly validates our high order fluid model, thus supporting its correct theoretical derivation and numerical implementation. The results of the first order fluid model with local field approximation, as usually used for streamer discharges, show considerable deviations. Furthermore, we study the inaccuracies of simulation results caused by an inconsistent implementation of transport data into our high order fluid model. We also demonstrate the importance of the energy flux term in the high order model by comparing with results where this term is neglected. Finally, results with an approximation for the high order tensor in the energy flux equation is found to agree well with the PIC/MC results for reduced electric fields up to 1000 Townsend, as considered in this work.Comment: 26 pages, 11 figure

Third-order transport coefficients for localised and delocalised charged-particle transport

We derive third order transport coefficients of skewness for a phase-space kinetic model that considers the processes of scattering collisions, trapping, detrapping and recombination losses. The resulting expression for the skewness tensor provides an extension to Fick's law which is in turn applied to yield a corresponding generalised advection-diffusion-skewness equation. A physical interpretation of trap-induced skewness is presented and used to describe an observed negative skewness due to traps. A relationship between skewness, diffusion, mobility and temperature is formed by analogy with Einstein's relation. Fractional transport is explored and its effects on the flux transport coefficients are also outlined.Comment: 11 pages, 3 figure

Ab-initio electron scattering cross-sections and transport in liquid xenon

Ab-initio electron - liquid phase xenon fully differential cross-sections for electrons scattering in liquid xenon are developed from a solution of the Dirac-Fock scattering equations, using a recently developed framework [1] which considers multipole polarizabilities, a non-local treatment of exchange, and screening and coherent scattering effects. A multi-term solution of Boltzmann's equation accounting for the full anisotropic nature of the differential cross-section is used to calculate transport properties of excess electrons in liquid xenon. The results were found to agree to within 25% of the measured mobilities and characteristic energies over the reduced field range of 10^{-4} to 1 Td. The accuracies are comparable to those achieved in the gas phase. A simple model, informed by highly accurate gas-phase cross-sections, is presented to transform highly accurate gas-phase cross-sections to improve the liquid cross-sections, which was found to enhance the accuracy of the transport coefficient calculations.Comment: 26 pages, 9 figures. arXiv admin note: text overlap with arXiv:1503.0037

Positron transport: the plasma-gas interface

Motivated by an increasing number of applications, new techniques in the analysis of electron transport have been developed over the past 30 years or so, but similar methods had yet to be applied to positrons. Recently, an in-depth look at positrontransport in pure argon gas has been performed using a recently established comprehensive set of cross sections and well-established Monte Carlo simulations. The key novelty as compared to electron transport is the effect of positronium formation which changes the number of particles and has a strong energy dependence. This coupled with spatial separation by energy of the positron swarm leads to counterintuitive behavior of some of the transport coefficients. Finally new results in how the presence of an applied magnetic field affects the transport coefficients are presented.This work was performed under MNTRS Project No. 141025

A comparison of 3D particle, fluid and hybrid simulations for negative streamers

In the high field region at the head of a discharge streamer, the electron energy distribution develops a long tail. In negative streamers, these electrons can run away and contribute to energetic processes such as terrestrial gamma-ray and electron flashes. Moreover, electron density fluctuations can accelerate streamer branching. To track energies and locations of single electrons in relevant regions, we have developed a 3D hybrid model that couples a particle model in the region of high fields and low electron densities with a fluid model in the rest of the domain. Here we validate our 3D hybrid model on a 3D (super-)particle model for negative streamers in overvolted gaps, and we show that it almost reaches the computational efficiency of a 3D fluid model. We also show that the extended fluid model approximates the particle and the hybrid model well until stochastic fluctuations become important, while the classical fluid model underestimates velocities and ionization densities. We compare density fluctuations and the onset of branching between the models, and we compare the front velocities with an analytical approximation

Spatially hybrid computations for streamer discharges: II. Fully 3D simulations

We recently have presented first physical predictions of a spatially hybrid model that follows the evolution of a negative streamer discharge in full three spatial dimensions; our spatially hybrid model couples a particle model in the high field region ahead of the streamer with a fluid model in the streamer interior where electron densities are high and fields are low. Therefore the model is computationally efficient, while it also follows the dynamics of single electrons including their possible run-away. Here we describe the technical details of our computations, and present the next step in a systematic development of the simulation code. First, new sets of transport coefficients and reaction rates are obtained from particle swarm simulations in air, nitrogen, oxygen and argon. These coefficients are implemented in an extended fluid model to make the fluid approximation as consistent as possible with the particle model, and to avoid discontinuities at the interface between fluid and particle regions. Then two splitting methods are introduced and compared for the location and motion of the fluid-particle-interface in three spatial dimensions. Finally, we present first results of the 3D spatially hybrid model for a negative streamer in air

Streamer propagation in the atmosphere of Titan and other N2:CH4 mixtures compared to N2:O2 mixtures

Streamers, thin, ionized plasma channels, form the early stages of lightning discharges. Here we approach the study of extraterrestrial lightning by studying the formation and propagation of streamer discharges in various nitrogen-methane and nitrogen-oxygen mixtures with levels of nitrogen from 20% to 98.4%. We present the friction force and breakdown fields Ek in various N2:O2 (Earth-like) and N2:CH4 (Titan-like) mixtures. The strength of the friction force is larger in N2:CH4 mixtures whereas the breakdown field in mixtures with methane is half as large as in mixtures with oxygen. We use a 2.5 dimensional Monte Carlo particle-in-cell code with cylindrical symmetry to simulate the development of electron avalanches from an initial electron-ion patch in ambient electric fields between 1.5Ek and 3Ek. We compare the electron density, the electric field, the front velocities as well as the occurrence of avalanche-to-streamer transition between mixtures with methane and with oxygen. Whereas we observe the formation of streamers in oxygen in all considered cases, we observe streamer inceptions in methane for small percentages of nitrogen or for large electric fields only. For large percentages of nitrogen or for small fields, ionization is not efficient enough to form a streamer channel within the length of the simulation domain. In oxygen, positive and negative streamers move faster for small percentages of nitrogen. In mixtures with methane, electron or streamer fronts move 10-100 times slower than in mixtures with oxygen; the higher the percentage of methane, the faster the fronts move.Comment: 34 pages, 11 figures, 1 tabl