1. Introduction
The development of an electrical discharge and the subsequent plasma formation is strongly dependent on the operating pressure and gas composition. Non-thermal plasmas, including low-pressure DC glows, RF discharges, and atmospheric pressure discharges, have many applications and serve as enabling technology in critical manufacturing processes [
1,
2,
3,
4,
5]. The transient phase of an atmospheric pressure discharge consists of an avalanche that leads to the streamer phase. If left uninterrupted, it will eventually lead to an arc formation [
6]. In most applications, the discharge is terminated: In repetitive nano-second pulse discharges, the pulse widths are short enough to avoid arc formation, and in dielectric barrier discharges, the charging of the dielectric quenches the micro discharge. From the perspective of plasma chemistry applications in such diverse areas as plasma medicine and plasma-assisted combustion, the understanding of reactive species production during the streamer phase is important [
7,
8]. The applications of non-thermal plasmas are based on the production of excited species, photo emission, and reactive radicals at or near ambient temperature [
9,
10,
11,
12]. These types of plasma are partially ionized gases, where the free electrons and heavy ions gain energy from the electric field and undergo elastic or inelastic collision with background particles or walls, leading to energy loss. Due to the large difference in the mass ratio, most of the energy gain from the field is by electron transport across potential gradients. In certain applications with a changing gas composition, a computationally efficient model for the streamer development will help the implementation of complex plasma chemistry phenomena such as repetitively pulsed plasma-assisted combustion [
7].
Non-thermal plasmas can reach a steady state in discharges, such as low-pressure DC glows or radio-frequency discharges, and others are classed as transient discharges, such as an atmospheric pressure streamer discharge. Near atmospheric pressure, electrical discharges produce spatially and temporally varying space charges that substantially alter the applied electric field and impose constraints on numerical models due to sharp spatial gradients. Non-local electron kinetics play an important role in low-pressure capacitively coupled RF discharge and low-pressure DC glows. We are particularly interested in studying the breakdown of atmospheric pressure gases. Non-thermal or cold plasma at atmospheric pressure forms the basis for many applications, including manufacturing, plasma medicine, disinfection, etc. [
1,
2,
12].
The streamer mechanism was first proposed by Raether [
13] and Leob and Meek [
14] to explain electrical breakdown at high pressures. Since then, a large number of both experimental and numerical studies have resulted in a better understanding of its formation and propagation [
15,
16,
17,
18,
19,
20,
21,
22,
23]. Theoretical efforts are constrained by the fact that the mathematical description of space charge-dominated transport is difficult to deal with because of the sharp density and field gradients. Most approaches to the modeling of streamers can be lumped under kinetic or fluid approaches. In kinetic models, the Boltzmann equation, coupled with the Poisson equation, is solved for the phase space of electrons. Alternately, Monte Carlo simulations with a detailed particle transport using collisional cross-sections are used to determine the particle phase space. To cover the full six-dimensional phase space is computationally expensive, and this approach has found very little use in general. A computationally tractable hybrid approach is particle-in-cell Monte Carlo simulations (PIC-MCCs), where a large number of electrons are followed in the phase space using Monte Carlo techniques used to simulate collisions, and the electric field is obtained from the solution of Poisson’s equation from the charged particle densities lumped in an appropriate mesh [
21]. The 3-D PIC-MCCs have been particularly successful in modeling stochastic fluctuations, leading to the branching of streamers being observed experimentally [
18]. Several articles have been published exploring the validity of the modeling approach for plasma fluid models: Nijdam et al. have reported numerical modeling including the pros and cons of the particle-in-cell (PIC) and fluid models [
19]. Kim et al. benchmarked the PIC, fluid, and hybrids models by comparing simulation results with experimental results for plasma displays, capacitively coupled plasma, and inductively coupled plasma. They concluded that despite progress in the modeling and simulation of low-temperature plasma, these models still need improvements [
20].
In plasma fluid models, the plasma hydrodynamics is described by macroscopic quantities such as electron density, drift velocity, and mean electron energy. The fluid models can be theoretically constructed by taking the velocity/energy moments of the Boltzmann equation. The first three moments gives the particle, momentum, and energy conservation equations to describe the plasma hydrodynamics, and depending on the number of moments considered, appropriate closure approximations are required. The first-order drift–diffusion model based on local field approximation (LFA) has been used with some success in predicting and reproducing experimental results of the formation and propagation of streamer channels [
15,
16,
17,
18,
19,
20,
21,
22,
23]. It has been reported that the assumption of local equilibrium of the electron energy deviates significantly at fast-changing ionization fronts with steep density and field gradients [
19]. In the hydrodynamic description of non-thermal plasma, the electron transport and collisional rate coefficients are commonly parameterized by the local field or the local mean energy. For most commonly used gases, the transport and collision rates can be readily obtained from the two-term solution of the Boltzmann equation from electron impact cross-section data [
24,
25].
The second-order drift–diffusion model considers electron energy transport where the parameters are based on local mean energy approximation (LMEA), and several reports concluded that this model gives better results at streamer fronts. Luque and Ebert reviewed density models for streamer discharge simulation, detailing their physical foundation, their range of validity, and the most relevant algorithm employed in solving them [
21]. Markosyan et al. compared plasma fluid models for one-dimensional streamer ionization fronts and compared them to the PIC model [
22]. They found the local energy approximation and a higher-order model were in better agreement with the PIC simulations, and the local field approximation gave reasonably close results. Gruber et al. examined the local field and local energy models for the simulation of low-pressure DC glow and capacitively coupled RF discharges at low pressures (10 and 100 Pa) in argon and oxygen. They concluded that the LFA method is not recommended for gas discharge modeling in general at this pressure due to the inadequacy of the drift–diffusion approximation, and their results should be checked against experimental data or benchmark approaches [
23].
We are interested in examining the plasma fluid models that would be suitable for investigating streamers near atmospheric pressures under ambient conditions. Although several publications have investigated this question, there are few results on the impact of the fluid model parametrization on excited species production. Our approach is not so much to replicate experimental results or streamer branching but come up with a suitable model to predict the important characteristics of a streamer that can inform the modeling of applications such as plasma medicine and plasma-assisted combustion. The purpose of the current paper is to understand under what condition are local field approximation and local mean energy approximation valid parameters for the transport and rate coefficients. In this paper, we simulate streamer development and propagation in nitrogen by using three different parametrization schemes and compare the important characteristics, such as excited species generation. Most applications are under ambient conditions; therefore, the study of nitrogen gas can serve as a good model.
2. Fluid Models for Streamer Discharges
During the transient phase, the heavy particles do not gain energy in this short period, and the neutral gas and ion are at or near room temperature. Also, in the time scale of interest (a few ns), the ions can be considered to be stationary compared to the lighter electrons. Both the first-order and the second-order fluid models for a non-attaching gas include the following particle conservation equations. In the first-order model, the parameter
is the local reduced electric field, E/N, and in the second-order model, it is the local mean electron energy,
ε [
15,
22].
The quantity
S represents various ion/electron source or sink mechanisms such as photoionization, recombination, attachment, or remnant space charge in repetitive discharges. In the absence of a magnetic field and assuming the velocity of the electrons is large compared to the slow species and the plasma is isothermal, the particle flux can be obtained from the momentum conservation equation and is given by [
15,
19]
where
and
are the electron and positive ion density, respectively,
is the electron mobility,
is electron diffusion coefficient, and
is the ionization frequency.
In a slowly varying electric field where the magnetic field can be neglected, the electric field
is obtained from the solution of the Poisson equation [
15].
where
is the unsigned electron charge and
is the free space permittivity.
In the second-order fluid model, the mean electron energy,
, is determined from the energy conservation equation [
22],
where
and
are the electron energy-dependent collision rate coefficient and energy loss per electron per collision for the
jth collision process. The first term in the right-hand side of Equation (5) is the convective term, the second term is the energy gained from the electric field, and the third term is the energy loss due to inelastic collisions. This form of the energy conservation equation is derived by assuming the electron pressure tensor is isotropic [
23]. A later part of this paper examines the contribution of each of these terms to the time evolution of the electron energy.
In streamer discharges, the ionization front propagates at speeds several times higher than the local drift velocity and experiences steep spatial gradients due to a rapid growth in ionization. It has been suggested that the first-order model is not adequate to describe the formation and propagation of these discharges, and a second-order model that calculates the local mean energy should be used [
24]. The local mean energy can then be used as a parameter to estimate the transport parameters and rate coefficients.
In this article, three different parametrizations shown in
Table 1 are investigated to understand the differences between the schemes in predicting streamer characteristics. The LMEA and the hybrid methods require the solution of the energy conservation equation. We introduce a new parameterization scheme (hybrid) where the mobility and diffusion are determined from the local electric field, which is readily available, and the electron impact rates such as ionization and excitation from the local electron energy. The justification for proposing this scheme is from a previous set of studies that show that the electron drift tracks the local electric field more closely compared to the electron energy [
25].
Table 1.
Parametrization schemes.
Table 1.
Parametrization schemes.
Local Field Approximation (LFA) | | | |
---|
Local Mean Energy Approximation (LMEA) | , | | |
Hybrid | , | | |
3. Results and Discussion
The simulations were performed in two dimensions with azimuthal symmetry. The transport parameters and rate coefficients were determined using the open-source Boltzmann solver, BOLSIG+, with Lxcat nitrogen cross sections [
26,
27]. The Boltzmann solver solves for the electron energy distribution function for a given reduced electric field (E/N), which also provides the corresponding electron energy. The transport coefficients and rates can then be parameterized with either the reduced field or the electron energy.
The discharge consists of a gap of 5 mm filled with atmospheric-pressure nitrogen. The results presented here are for an applied step voltage of
(186 Td), which resulted in an anode-directed streamer. The set of Equations (1)–(5) was solved numerically using the finite difference method. For all density variables, homogeneous Neuman boundary conditions, i.e., zero derivatives at the electrode boundaries, are prescribed [
22]. A more physical boundary condition would require detailed knowledge of the electrode and secondary electrode processes. In such discharges, the cathode sheath has very little impact on the bulk properties of the streamer, which is of interest in most atmospheric pressure plasmas. The Flux-Corrected Transport (FCT) method proposed by Boris and Book was used for the convective term of both the electron density and electron energy equations [
28,
29]. This method is particularly suitable for handling the steep density and field gradients encountered in streamer propagation. The contributions from the diffusion terms and from the impact ionization terms in Equation (1) were calculated at each time step and added to the convective term. The details can be found in reference [
15]. A uniform grid spacing of 500 was used both in the z and r directions, giving a spatial resolution of 10 μm. The details of the method as applied to streamers have been extensively reported [
15,
16,
17].
The following boundary conditions were used for the solution of the Poisson equation.
The Poisson equation is solved for the electric potential by the Successive Over-Relaxation (SOR) method [
3]. This is an iterative method and converges rapidly, as there is small perturbation from the previous time step. For our simulation, a relaxation factor ω = 1.9475 gave the fastest convergence. The maximum relative error at any grid point was set to 10
−5 for the convergence criterion.
The results presented here are for anode-directed streamers. The results were very similar for cathode-directed streamers. In order to bypass the avalanche phase, an initial neutral plasma is placed at the cathode, which represents the space charge formed. Details of this method and the background seeding of electrons to represent ionization can be found in references [
15,
16].
Figure 1,
Figure 2 and
Figure 3 show the contour plots of electron density, electron energy, and electric field contours at different times for the fluid model with hybrid parametrization. These plots are typical of steamer propagation, where the streamer tip shows high gradients for electron density, electron energy, and electric field. In the streamer bulk, away from the tip, the electron energy and electric field is fairly constant.
The three different types of parametrizations result in similar contours plots. For comparison, the on-axis density and the electric field were plotted for streamer formation using the three different parametrizations at a nearly equal time and are shown in
Figure 4 and
Figure 5, respectively. Since the time steps for the simulation are determined from the current Courant–Friedrichs–Lewy criteria, the comparison plots are not exactly at the same instant [
30]. The electron density and the electric field as determined from three different parametrizations give very similar results with only a minor difference, which has very little impact on the development and propagation of the streamer.
The streamer speed as a function of time is shown in
Figure 6 for the three different parametrizations. Again, we see very little difference in the speed, irrespective of the method used for solving the fluid equations. The speed of the ionization front increases with time, as the field enhancement at the tip of the streamer also increases with time as the streamer propagates. The higher the field enhancement, the quicker is the plasma density build-up due to electron impact ionization. The magnitude of the velocities of the streamer in nitrogen is close to experimentally reported values: Wagner reported an anode-directed velocity of 0.4 × 10
6 m/s at 156 Td in atmospheric-pressure nitrogen, and Chalmers et al. reported an anode-directed velocity of 0.1 to 0.4 × 10
6 m/s in the range of 126 to 156 Td [
31,
32].
The spatial and temporal evolution of the energy obtained from the solution of the energy equation were compared to the energy predicted by the local electric field using the equilibrium relationship between the reduced electric field (E/N) in Td and electron energy,
in nitrogen as obtained from the solution of the Boltzmann equation shown below.
The axial plot of the energy for the three different parameterizations is shown in
Figure 7,
Figure 8 and
Figure 9 for the LFA, hybrid, and LMEA, respectively. Remarkably, for all three parametrizations the electron energy agrees very well with the local electric field prediction at the streamer tip where the density and field gradients are the highest. There is a slight difference in the bulk of the streamer due to slight electron cooling predicted by the energy equation, which is discussed in a later section. The very steep rise in the electron density and the electric field in front of the streamer introduces some numerical oscillations, particularly in the later stages of the simulation. These oscillations will not have a significant impact in determining the excited species concentration.
In most applications where streamer-type discharges are used to generate excited species, an estimation of the species concentration is important in determining the subsequent plasma chemical pathways. Shown in
Figure 10 is the G-factor, which is the number of radicals produced per 100 eV of electrical energy input, for three of the nitrogen excited states:
,
, and
[
33]. The three model parametrizations predict very similar G-factors, which remain fairly constant with time. The current increases with time as the streamer propagates along the axial direction. Therefore, the electrical energy input increases with time, and radical densities increase with time, although the G-factor does not change with time.
The contribution to the energy change at spatial points along the axis was studied by looking at the first term on the right-hand side of the energy conservation equation, which is the convective term, and the net energy gain, which is the sum of the second (energy gained from the electric field) and the third (energy lost due to inelastic process) terms. The energy gain from the electric field and the inelastic loss are almost equal and an order of magnitude higher than the convective contribution.
Figure 11 shows the relative contribution due to the convective term, gain–loss term, and the net energy density rate for the hybrid and LMEA parametrizations. The LFA is not shown because the energy equation is not relevant to that method. At the streamer tip, there is a net transport of electron energy due to convection from the back to the front due to electron transport. The net effect is very similar in both parametrizations, and there is a net cooling behind the streamer tip. This cooling effect is seen in the energy plots (
Figure 7,
Figure 8 and
Figure 9) where the energy in the bulk shows a slight decline further away from the cathode.
Our results suggest that the electron energy in a nitrogen streamer at atmospheric pressure quickly reaches equilibrium with the electric field, and the convective transport of energy does not have a significant impact. The various characteristics of the streamer discussed here show very little dependence on the parameterization. Markosyan et al. have performed similar studies in a 1-D streamer model. They compared different fluid models with a 1-D PIC simulation. Although the 1-D models have inherent shortcomings due to the gross simplification in estimating the spatial profile of the electric field, their findings are similar to what we observe in our simulations. A similar conclusion was reached by Wang et al., where they looked at particle and fluid models for streamer discharges in air [
34]. Li et al. looked at the deviations from the LFA in negative streamer heads in nitrogen streamer heads when compared to particle models. They concluded that the largest discrepancy is in the leading edge of the front where the electron density is very low, and the electric field is not screened [
35]. As such, this result will have a minimum impact on estimating the overall radical generation in a streamer discharge. Our work confirms the choice of parametrization has very little impact on the species generation in the discharges studied.