**Chemical Kinetics in Air Plasmas at Atmospheric Pressure**

Claudia Riccardi and Ruggero Barni

*Department of Physics G. Occhialini, University of Milano-Bicocca, Milan, Italy*

#### **1. Introduction**

184 Chemical Kinetics

The authors are deeply grateful to all their collaborators and students that over the years have partecipated in their research studies, contributing to the success of the enterprise, or more poetically adventure, of the PlasmaPrometeo Center of Excellence. We could enlist M.Piselli, S.Zanini, P.Tosi, A.Quintini, F.Picciolo, C.Broggi, R.Benocci, M.Orlandi,

Barni, R.; Broggi, C. ; Benocci, R. & Riccardi, C. (2006). Chemical kinetics of an

Barni, R.; Quintini, A. ; Piselli, M. & Riccardi, C. (2009). Experimental study of hydrogen

Barni, R. & Riccardi, C. (2009). Hydrocarbon plasma reforming through intermittent spark

discharges, *High Temperature Material Processes*, Vol. 13, pp. 119-127

Davies, D.K. ; Kline, L.E. & Bies, W.E. (1989) *Journal of Applied. Physics*, Vol. 65, pp. 3311 Kline, L.E. ; Partlow, W.D. & Bies, W.E. (1989) *Journal of Applied. Physics*, Vol. 65, pp. 70-77 Gruen, D.M. ; Zuiker, C.D. ; Krauss, A.R. & Pan, X. (1995). *Journal of Vacuum Science &*

Kline, L.E. ; Partlow, W.D. & Bies, W.E. (1988) *Journal of Applied. Physics*, Vol. 65, pp. 70-77

Lieberman, M.A. & Lichtenberg, A.J. (1994). *Principles of Plasma Discharges and Materials*

Smith, J.M. (1987). *Introduction to Chemical Engineering Thermodynamics*, McGraw-Hill, New

Tosi, P. ; Bassi D. ; Brunetti, B. & Vecchiocattivi, F. (1995). *International Journal of Mass*

argon/methane plasma in a hydrogen reforming reactor, *European Physical Journal*

plasma reforming by intermittent spark discharges, *Journal of Applied. Physics*, Vol.

A.Brandirali just to mention only a few involved in the topics of this chapter.

Angus, J.C. & Hayman, C.C. (1988). *Science,* Vol. 241, pp. 913

Benson, S.W. (1982). *Thermochemical Kinetics*, Wiley, New York US Chantry, P.J. (1987). *Journal of Applied. Physics*, Vol. 62, pp. 1141

Li, H.Q. ; Zou, J.J. & Lin, C.J. (2005) *Progress in Chemistry*, Vol. 17, pp. 69

Mutsukura, N. & Yoshida, K. (1996). *Diamond Related Materials,* Vol. 5, pp. 919

Paulmier, T. & Fulcheri, L. (2005). *Chemical Engineering Journal*, Vol. 106, pp. 59 Press, W.H. (1985). *Numerical Recipes*, Cambridge University Press, Cambridge UK Raizer, Y.P. (1991). *Gas Discharge Physics,* Springer-Verlag, Heidelberg Germany Rhallabi, A. & Catherine, Y. (1991) *IEEE Transactions on Plasma Science*, Vol. 19, pp. 270 Riccardi, C. ; Barni, R. & Fontanesi, M. (2001). *Journal of Applied. Physics*, Vol. 90, pp. 3735 Riccardi, C. ; Barni, R. ; Fontanesi, M. & Tosi, P . (2000). Gaseous precursors of diamond-like Carbon films in CH4/Ar plasmas, *Chemical Physics Letters*, Vol. 392, pp. 66-70 Riccardi, C. ; Barni, R. ; Sindoni, E.; Fontanesi, M. & Tosi, P. (2001). *Vacuum*, Vol. 61, pp. 211 Semaniak, J.; Andfersen, L.H.; Heber, O.; Kella, D., Pedersen, H.B., Schmidt, H.T. & Zajfman,

*- Applied. Physics*, Vol. 35, pp. 135-143

*Technology A*, Vol. 13, pp. 1628

*Processing*, J.Wiley & Sons, New York US

Mc Daniel,E.W. (1993). *Atomic Collisions,* Wiley, New York US

D. (1998). *Astrophysical Journal,* Vol 498, pp. 886-895

*Spectrometry and Ion Processes*, Vol. 149/150, pp. 345-353

Shih, H.C. (1993). *Diamond Related Materials,* Vol. 2, pp. 531

York US

103, pp. 063302.1-9

**5. Acknowledgment**

**6. References**

The chapter, after a general introduction concerning plasmas produced in air at atmospheric pressure, describes the chemical kinetics happening in the plasma gas-phase, focusing on the nitrogen oxides mechanism of formation and elimination. After a general review of the chemical kinetic pattern, we discuss dielectric barrier discharges and their use for NOx remediation. In fact, research concerning the use of plasmas for the promotion of chemical reactions has emerged as one of the most prominent application field. Gas treatment aimed to toxic waste disposal or to pollutant emission reduction is an important field where plasma processing can prove to be a very effective technology. We have studied the characteristics of plasma reactors based on dielectric barrier discharges (DBD) operating in a streamer regime, as well as those based on high voltage electrical discharges operated in an intermittent spark regime. Both kind of discharges have been considered for the promotion of chemical reactions aimed to nitrogen oxides removal form a diluted gas flow. Detailed modeling of the chemical kinetics happening in an air plasma gas-phase has been investigated in our previous research. The evolution of the gas-phase composition during and after the interaction in the discharge region was calculated under different operating conditions. The results have been used to predict the NOx reduction achievable and to study the role played by the different discharge parameters in order to determine the optimal interaction between plasma and gas mixture. Comparison between streamer and spark discharges could be performed too.

#### **2. Discharges at atmospheric pressure**

Let's consider a system where at least two electrodes (conductive, in general metal surfaces) are separated by a gap filled with air at atmospheric pressure. Then one of the electrode electrical potential is increased respect to the second, which we could suppose, without any lack of generality to be grounded. An electric field arises in the gas gap, whose magnitude and space distribution depend on the potential applied and the geometry of the system. At large enough values of the field a discharge could spark in the gas. When an electrical breakdown happens in atmospheric pressure air the discharge quickly evolves towards to a condition allowing large electric current flows and strong heating of the gas. Self-sustaining discharges at such high pressure could not develop as the usual glow regime, since heating

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 187

occurs and a steady plasma state is formed, independent on time and, in first approximation, determined only by the effective, i.e. rms, electric field value and not from the wave frequency (Raizer, 1991). Then modeling as well as the pattern of the chemical kinetics associated with these plasmas could be grasped from the trends observed in d.c. corona discharges. The oscillating electric field however prevent most of the charges to be lost at the walls or the electrodes by drifts alone and the plasma state is less influenced by badly known parameters linked to electrode surface and material properties as in general it happens in d.c. configurations. This make the realization of suitable plasma sources in form

The main characteristic of a Dielectric Barrier Discharge (DBD) device lays in the setup of the electrodes. In particular the setup includes the presence of a dielectric layer within the discharge gap distance d (in general in the 0.1-5 mm range for an atmospheric plasma) insulating at least one of the electrodes. The dielectric presence forbids the flow of a direct current in the gap, preventing the formation of arcs or sparks, as well as of a steady plasma state. To sustain the discharge an a.c. voltage or a pulsed one is applied to the electrodes at a frequency ranging from several hundreds of hertz to a few hundreds of kHz. While many different kinds of DBD exist, in general, for a such high value of the product pd of gas pressure p and the electrode gap width d, the discharge operates in a streamer regime (Raizer, 1991). Just to make a little example, let's consider a discharge gap d of 1 mm supplied by an HV voltage of 5 kV, which makes an electric field of 5 MV/m whether it would be uniform. Using the data reported in Fig. 2.2.1 it easy to see that the ionization length is well below 0.1 mm. Actually the product α·d turns out 21.3 exceeding the threshold value of 18-20 , known as the Meek criterion, for streamer discharge onset (Raizer, 1991).

But what characterises most the steamer regime? These discharges are built up of many repeated microdischarges, each leading to the formation of narrow discharge filament (with an electron density in the range 1012-1017 cm-3) of some nanoseconds duration (Kogelschatz, 2003). In particular a plane electrode DBD at near atmospheric pressure which operates in streamer regime, under an a.c.voltage supply displays an uniform glow filling the gap. At a closer inspection the glow reveals itself as built up with a large number of current filaments appearing through the gap at seemingly random locations and not exactly at the same time during an active phase which happens each voltage half-cycle (Kogelschatz, 2003). The uniformity is just brought in by the high rate of repetition and superposition of such microdischarges. Depending on the operating conditions, at each half-cycle a new streamer head could propagate in the opposite direction starting from the same location when the voltage is reversed. On the other hand, during the active phase, the streamers cannot propagate on the same path, being inhibited by the space charge accumulated on the dielectric surfaces. This has important consequences on the chemical kinetics and its modeling since streamers can be considered as isolated and repeating at frequency which is the double of the a.c. voltage frequency. Another key feature is that, at least in small gaps like the one we are considering, they appear to be very thin (radius R ~ 100 μm) while their average distance is much larger (a few mm, of the order of the gap width) as discussed in literature (Raizer, 1991). So, as a first approximation one could chose to simulate the chemical kinetics induced in a single isolated streamer and its repetition at a fixed

frequency. The results of this approach are examined in the following.

of plasma jets, pencils or needles much easier and user-friendly (Mariotti, 2010).

**2.2 Dielectric barrier discharges**

and ionization instability leads to the contraction of the discharge electrical current channel and, if the system could sustain it, the discharge evolves to a so called arc. Under such conditions the air in the gap becomes very hot, reaching a few thousands of K and electrical current could exceed hundreds Ampers (Raizer, 1991).

These kinds of discharges have interesting chemical kinetics in their own, being almost thermal equilibrium systems, but their behaviour as chemical kinetics reactors bears more resemblances with that of a hot gas than of a plasma. Moreover such devices require very large power and the high temperatures imply a considerable stress of the materials used, which somewhat limits their use as processing tools in material chemistry.

#### **2.1 Corona and radiofrequency discharges**

More interesting from the point of view of plasma chemical kinetics is another kind of discharge, the so called corona. However corona discharges occur only whether the electric field is sharply non-uniform in space, which means that the electrode setup should reflect such condition. The easiest way to accomplish this is when one of the electrode size is much smaller than the inter-electrode or gap distance. Point or wire electrodes are the best instances of that. Corona are self-sustaining discharges, producing plasmas in steady conditions and burning in a stable manner for quite long periods (Raizer, 1991) . Then chemical kinetics have time to evolve and an interesting pattern of radical formation and oxidation could develop in the air of the gap. On the other hand, corona discharges are inherently disuniform in space, requiring geometry to be taken into account and complicating the task of modeling. The result is, in a sense, to obscure partially the true chemical kinetics effects.  Then their use in material processing is somewhat limited, even if, mainly in non specialized literature, the term corona is often employed and abused for other kind of discharges and setups. In the following sections we will discuss briefly the chemical kinetics associated with their plasmas. Although the use of an oscillating voltage as supply could make several differences in the gas breakdown and in the discharge formation, there are striking similarities in the plasma state which is produced. Indeed at the somewhat large density of an atmospheric pressure gas the oscillation frequency is order of magnitude smaller than the elastic collision frequency (in air at STP this is about 3 THz or 3x1012 Hz). Under such conditions averaging

Fig. 2.1.1. A radiofrequency atmospheric pressure air plasma produced by a plasma needle.

occurs and a steady plasma state is formed, independent on time and, in first approximation, determined only by the effective, i.e. rms, electric field value and not from the wave frequency (Raizer, 1991). Then modeling as well as the pattern of the chemical kinetics associated with these plasmas could be grasped from the trends observed in d.c. corona discharges. The oscillating electric field however prevent most of the charges to be lost at the walls or the electrodes by drifts alone and the plasma state is less influenced by badly known parameters linked to electrode surface and material properties as in general it happens in d.c. configurations. This make the realization of suitable plasma sources in form of plasma jets, pencils or needles much easier and user-friendly (Mariotti, 2010).

#### **2.2 Dielectric barrier discharges**

186 Chemical Kinetics

and ionization instability leads to the contraction of the discharge electrical current channel and, if the system could sustain it, the discharge evolves to a so called arc. Under such conditions the air in the gap becomes very hot, reaching a few thousands of K and electrical

These kinds of discharges have interesting chemical kinetics in their own, being almost thermal equilibrium systems, but their behaviour as chemical kinetics reactors bears more resemblances with that of a hot gas than of a plasma. Moreover such devices require very large power and the high temperatures imply a considerable stress of the materials used,

More interesting from the point of view of plasma chemical kinetics is another kind of discharge, the so called corona. However corona discharges occur only whether the electric field is sharply non-uniform in space, which means that the electrode setup should reflect such condition. The easiest way to accomplish this is when one of the electrode size is much smaller than the inter-electrode or gap distance. Point or wire electrodes are the best instances of that. Corona are self-sustaining discharges, producing plasmas in steady conditions and burning in a stable manner for quite long periods (Raizer, 1991) . Then chemical kinetics have time to evolve and an interesting pattern of radical formation and oxidation could develop in the air of the gap. On the other hand, corona discharges are inherently disuniform in space, requiring geometry to be taken into account and complicating the task of modeling. The result is, in a sense, to obscure partially the true chemical kinetics effects.  Then their use in material processing is somewhat limited, even if, mainly in non specialized literature, the term corona is often employed and abused for other kind of discharges and setups. In the following sections we will discuss briefly the chemical kinetics associated with their plasmas. Although the use of an oscillating voltage as supply could make several differences in the gas breakdown and in the discharge formation, there are striking similarities in the plasma state which is produced. Indeed at the somewhat large density of an atmospheric pressure gas the oscillation frequency is order of magnitude smaller than the elastic collision frequency (in air at STP this is about 3 THz or 3x1012 Hz). Under such conditions averaging

Fig. 2.1.1. A radiofrequency atmospheric pressure air plasma produced by a plasma needle.

which somewhat limits their use as processing tools in material chemistry.

current could exceed hundreds Ampers (Raizer, 1991).

**2.1 Corona and radiofrequency discharges**

The main characteristic of a Dielectric Barrier Discharge (DBD) device lays in the setup of the electrodes. In particular the setup includes the presence of a dielectric layer within the discharge gap distance d (in general in the 0.1-5 mm range for an atmospheric plasma) insulating at least one of the electrodes. The dielectric presence forbids the flow of a direct current in the gap, preventing the formation of arcs or sparks, as well as of a steady plasma state. To sustain the discharge an a.c. voltage or a pulsed one is applied to the electrodes at a frequency ranging from several hundreds of hertz to a few hundreds of kHz. While many different kinds of DBD exist, in general, for a such high value of the product pd of gas pressure p and the electrode gap width d, the discharge operates in a streamer regime (Raizer, 1991). Just to make a little example, let's consider a discharge gap d of 1 mm supplied by an HV voltage of 5 kV, which makes an electric field of 5 MV/m whether it would be uniform. Using the data reported in Fig. 2.2.1 it easy to see that the ionization length is well below 0.1 mm. Actually the product α·d turns out 21.3 exceeding the threshold value of 18-20 , known as the Meek criterion, for streamer discharge onset (Raizer, 1991).

But what characterises most the steamer regime? These discharges are built up of many repeated microdischarges, each leading to the formation of narrow discharge filament (with an electron density in the range 1012-1017 cm-3) of some nanoseconds duration (Kogelschatz, 2003). In particular a plane electrode DBD at near atmospheric pressure which operates in streamer regime, under an a.c.voltage supply displays an uniform glow filling the gap. At a closer inspection the glow reveals itself as built up with a large number of current filaments appearing through the gap at seemingly random locations and not exactly at the same time during an active phase which happens each voltage half-cycle (Kogelschatz, 2003). The uniformity is just brought in by the high rate of repetition and superposition of such microdischarges. Depending on the operating conditions, at each half-cycle a new streamer head could propagate in the opposite direction starting from the same location when the voltage is reversed. On the other hand, during the active phase, the streamers cannot propagate on the same path, being inhibited by the space charge accumulated on the dielectric surfaces. This has important consequences on the chemical kinetics and its modeling since streamers can be considered as isolated and repeating at frequency which is the double of the a.c. voltage frequency. Another key feature is that, at least in small gaps like the one we are considering, they appear to be very thin (radius R ~ 100 μm) while their average distance is much larger (a few mm, of the order of the gap width) as discussed in literature (Raizer, 1991). So, as a first approximation one could chose to simulate the chemical kinetics induced in a single isolated streamer and its repetition at a fixed frequency. The results of this approach are examined in the following.

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 189

In numerous problems of plasma chemistry, electro-physics, gas lasers physics and atmospheric physics a detailed consideration of the kinetics of non-equilibrium plasmachemical processes in nitrogen-oxygen mixtures is essential (Raffaele-Addamo, 2003; Esena, 2005). Another relevant topic concerns the need to examine the ecological consequences of projects involving a prolonged use of powerful electromagnetic radiation in the atmosphere. Channels determining the balance of neutral and charged particles in a non-equilibrium plasma are highly diverse, making plasma chemical kinetics vastly different from that of an ordinary hot gas. In plasmas, the main chemical processes proceed on the basis of reactions with nitrogen and oxygen atoms being formed by electron impact dissociation of the initial molecules. In some cases, electronically excited molecules play an important role, as for instance, the influence of oxygen molecules in the state a' on the process of ozone formation. This makes relevant any work undertaken in an attempt to form a kinetic scheme of non-

As discussed above, a sensible starting point to simulate the chemical kinetics induced in a DBD is to study what happens in a single isolated streamer and its repetition at a fixed frequency. For the sake of simplicity we choose to discuss dry air. However the validity of the model can be extended to general electrodes configuration and more extensively to any plasma devices operating in a streamer discharge mode. The inherent difficulty encountered in developing a model for air discharges at high pressure is that one must solve a highly nonlinear coupled set of differential equations which may be both spatially and temporally dependent. In DBDs, moreover, the time scales involved span many order of magnitude, from the very fast connected with streamer formation and propagation to the very slow connected with ion recombination, chemical kinetics of neutral species and diffusion. So several approximation are worthy in order to make tractable the problem. From the chemical kinetics point of view, the process can be seen as consisting of an almost instantaneous phase during which energetic electrons produce mostly ions and atomic radicals, followed by a phase in which electrons cool down almost instantaneously and reach the set temperature or a temperature so low, that further dissociations are stopped. In the latter phase, slow chemical reactions between neutral and charged species happen in the gas phase, while they diffuse outside the streamer channel. Although existing models do predict that the shape as well as the peak value of the electric field pulse depends on the position along the streamer path, the differences are not large, apart from the region near the dielectric surfaces (Kulikovsky, 1997). As a first approximation, one could neglect such differences and thus let's suppose uniformity along the streamer path. From the mathematical point of view, then the model formulation used treats the streamer channel as a cylinder, with a circular cross-section. The radial profile is determined by diffusion and the only parameter describing evolution is time. The actual value of the streamer radius is treated as a free parameter of the simulation. We assume also a square pulse shape for the time evolution of the electric field ahead of the streamer. The pulse is therefore characterized by the time length (here 2 ns, based on the model in (Kulikovsky, 1998)) and by the peak electric field Emax. Although it appears as a very rough approximation, it embodies all the

**3. Chemical kinetics in atmospheric pressure air plasmas**

equilibrium discharge in nitrogen-oxygen mixtures (Kossyi, 1992).

**3.1 Modeling of the air kinetics in dielectric barrier discharges**

relevant dynamics for our aims.

Fig. 2.2.1. The ionization inverse length coefficient in atmosperic pressure air

Several models (see for instance Kulikovsky, 1997-1998) and a few experimental observations (Malik, 2008) of streamer development and propagation exist by now. Despite many differences, the overall picture appears to be the same in all kind of DBDs. The microdischarges develops in three steps. After a phase of charge accumulation in the head, coming from previously deposited charges or from electron multiplication (Raizer 1991), the stremer head starts moving towards cathode. Streamer propagation is very fast (v ~ 105 - 106 m/s) and is sustained by the electron avalanche multiplication in the strong field ahead of the streamer (Emax ~ 5-15 MV/m). After the head has reached the cathode, spreading its charge on the dielectric surface, in the gap a discharge filament is left. The streamer body consists of a thin channel of weakly ionized gas, almost quasi-neutral (n+ ~ n- ~ 1013-1016 cm-3) and driven by a small or just vanishing electric field. After the streamer head reaches the dielectric barrier insulating the cathode, the channel allows the flow of a limited discharge current, which is rapidly quenched by the space charge accumulating at the dielectric surface. At high pressures electrons moves along the electric field at the drift velocity and almost instantaneously (the corresponding relaxation times for velocity and energy, dictated by collisions with neutrals, are τ ~ 0.3 ps, τE ~ 50 ps) reach a mean energy determined by the local amplitude of the field. On the other hand ions inertia keeps them almost frozen and substantially cold, that is at the set temperature. At a fixed position along the streamer path, the electric field rises quickly as the streamer head is approaching and then decay almost as fast to the much smaller value it assumes in the streamer channel. This has important consequences on the chemical kinetics and its modeling too. Indeed the high energy bond-dissociation and ionization processes happens only through electron impact interactions, since neutral and ions are mostly cold. Then a simple cold plasma approximation is valuable. Moreover tha active phase happens only in a very short time when the streamer head is propagating. Since propagation is so fast, chemical kinetics could be considered freezed and happening independently at each position along the streamer channel (Siliprandi, 2008).

**0 5 10 15 20 25 30 35 40**

**Electric field strength (MV/m)**

Several models (see for instance Kulikovsky, 1997-1998) and a few experimental observations (Malik, 2008) of streamer development and propagation exist by now. Despite many differences, the overall picture appears to be the same in all kind of DBDs. The microdischarges develops in three steps. After a phase of charge accumulation in the head, coming from previously deposited charges or from electron multiplication (Raizer 1991), the stremer head starts moving towards cathode. Streamer propagation is very fast (v ~ 105 - 106 m/s) and is sustained by the electron avalanche multiplication in the strong field ahead of the streamer (Emax ~ 5-15 MV/m). After the head has reached the cathode, spreading its charge on the dielectric surface, in the gap a discharge filament is left. The streamer body consists of a thin channel of weakly ionized gas, almost quasi-neutral (n+ ~ n- ~ 1013-1016 cm-3) and driven by a small or just vanishing electric field. After the streamer head reaches the dielectric barrier insulating the cathode, the channel allows the flow of a limited discharge current, which is rapidly quenched by the space charge accumulating at the dielectric surface. At high pressures electrons moves along the electric field at the drift velocity and almost instantaneously (the corresponding relaxation times for velocity and energy, dictated by collisions with neutrals, are τ ~ 0.3 ps, τE ~ 50 ps) reach a mean energy determined by the local amplitude of the field. On the other hand ions inertia keeps them almost frozen and substantially cold, that is at the set temperature. At a fixed position along the streamer path, the electric field rises quickly as the streamer head is approaching and then decay almost as fast to the much smaller value it assumes in the streamer channel. This has important consequences on the chemical kinetics and its modeling too. Indeed the high energy bond-dissociation and ionization processes happens only through electron impact interactions, since neutral and ions are mostly cold. Then a simple cold plasma approximation is valuable. Moreover tha active phase happens only in a very short time when the streamer head is propagating. Since propagation is so fast, chemical kinetics could be considered freezed and happening independently at each position along the streamer

Fig. 2.2.1. The ionization inverse length coefficient in atmosperic pressure air

**10<sup>2</sup>**

channel (Siliprandi, 2008).

**10<sup>3</sup>**

**10<sup>4</sup>**

**Ionization Inverse length**

 **(m-1**

**)**

**10<sup>5</sup>**

**10<sup>6</sup>**

#### **3. Chemical kinetics in atmospheric pressure air plasmas**

In numerous problems of plasma chemistry, electro-physics, gas lasers physics and atmospheric physics a detailed consideration of the kinetics of non-equilibrium plasmachemical processes in nitrogen-oxygen mixtures is essential (Raffaele-Addamo, 2003; Esena, 2005). Another relevant topic concerns the need to examine the ecological consequences of projects involving a prolonged use of powerful electromagnetic radiation in the atmosphere. Channels determining the balance of neutral and charged particles in a non-equilibrium plasma are highly diverse, making plasma chemical kinetics vastly different from that of an ordinary hot gas. In plasmas, the main chemical processes proceed on the basis of reactions with nitrogen and oxygen atoms being formed by electron impact dissociation of the initial molecules. In some cases, electronically excited molecules play an important role, as for instance, the influence of oxygen molecules in the state a' on the process of ozone formation. This makes relevant any work undertaken in an attempt to form a kinetic scheme of nonequilibrium discharge in nitrogen-oxygen mixtures (Kossyi, 1992).

#### **3.1 Modeling of the air kinetics in dielectric barrier discharges**

As discussed above, a sensible starting point to simulate the chemical kinetics induced in a DBD is to study what happens in a single isolated streamer and its repetition at a fixed frequency. For the sake of simplicity we choose to discuss dry air. However the validity of the model can be extended to general electrodes configuration and more extensively to any plasma devices operating in a streamer discharge mode. The inherent difficulty encountered in developing a model for air discharges at high pressure is that one must solve a highly nonlinear coupled set of differential equations which may be both spatially and temporally dependent. In DBDs, moreover, the time scales involved span many order of magnitude, from the very fast connected with streamer formation and propagation to the very slow connected with ion recombination, chemical kinetics of neutral species and diffusion. So several approximation are worthy in order to make tractable the problem. From the chemical kinetics point of view, the process can be seen as consisting of an almost instantaneous phase during which energetic electrons produce mostly ions and atomic radicals, followed by a phase in which electrons cool down almost instantaneously and reach the set temperature or a temperature so low, that further dissociations are stopped. In the latter phase, slow chemical reactions between neutral and charged species happen in the gas phase, while they diffuse outside the streamer channel. Although existing models do predict that the shape as well as the peak value of the electric field pulse depends on the position along the streamer path, the differences are not large, apart from the region near the dielectric surfaces (Kulikovsky, 1997). As a first approximation, one could neglect such differences and thus let's suppose uniformity along the streamer path. From the mathematical point of view, then the model formulation used treats the streamer channel as a cylinder, with a circular cross-section. The radial profile is determined by diffusion and the only parameter describing evolution is time. The actual value of the streamer radius is treated as a free parameter of the simulation. We assume also a square pulse shape for the time evolution of the electric field ahead of the streamer. The pulse is therefore characterized by the time length (here 2 ns, based on the model in (Kulikovsky, 1998)) and by the peak electric field Emax. Although it appears as a very rough approximation, it embodies all the relevant dynamics for our aims.

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 191

where K are the reaction rates for the gas phase reactions and D is the diffusion coefficient. Due to the supposed uniformity, diffusion along the streamer axis is neglected. Adsorption on the dielectric barrier is included but, because of the much smaller extension of the streamer radius respect to the discharge gap, it turns out to be negligible. The choice of species included in the model was based on existing experimental information based on emission spectroscopy or mass spectroscopy. A sub-set of 8 neutral and 9 charged species was included in a preliminary study. Then the code was updated to include 20 neutral and

As already discussed, to simulate DBD plasmas a cold plasma approximation is suitable. So reaction rates for ion as well as neutral species are evaluated at the set temperature, that is 300 K. The effect of larger temperatures is addressed below. Since thermal energy is much smaller than the molecular bond energy, the initial step in the decomposition is presumed to involve mainly electron-impact dissociation. Electron-impact cross-sections used in this study were obtained from literature collected in (Riccardi, 2001; Barni, 2005), apart from attachment rates which are taken from recent experiments. Rate constants have been evaluated assuming a Maxwellian energy distribution function for electron, described by a single parameter, their temperature Te. Rate constants for neutral gas-phase, charge exchange and ion recombination reactions have been taken from literature referenced in (Riccardi, 2001; Guerra 1997) or estimated based on approximating techniques (Benson, 1982). In this report we assume a substantial vibrational excitation of nitrogen molecules, which allows direct dissociation from atomic O radicals (Guerra, 1997). A total of 121 reactions have been taken into account in the simulation. The diffusion loss rates are estimated by using a formulation due to (Chantry, 1987) in cylindrical coordinates. The radial density profile of the streamer is approximated in term of the lowest order normal mode and is integrated away. This formulation allows substituting the Laplacian operator with an effective diffusion length and makes the system of balance equations only time dependent. Diffusion of air from the outside the streamer is also included as a source term. Ion diffusion loss rates was corrected by enforcing the plasma quasi-neutrality, following an approach proposed for low pressure electronegative plasmas (Lichtenberg, 2000). Balance equations containing all the processes described above form a system of coupled differential equations which has been integrated to obtain the time evolution by using an adaptive Runge-Kutta routine. Results from the calculations are discussed below and have been

As a first step we start discussing the temporal evolution of the discharge gas phase using a sort of reference values for the free parameters. That is electron temperature was 4 eV, electron density reached in the microdischarge channel is 2x1015 cm-3 in a time length of 2 ns and the streamer radius was 100 µm. At first, we neglect diffusion losses, in order to look directly the evolution pattern of the chemical kinetics alone. Neutral as well charged species

As expected, the charged phase shows a very fast increase, almost exponential, in almost all species, as the streamer gets formed. Afterwards, as diffusion is not included, the system evolves under recombination alone, giving the overall decreasing trend as the inverse of time (that is linear in the log-log plot). The composition however displays some important

densities are displayed as a function of time respectively in figures on a log-log scale.

excited and 12 ionic states.

partially presented in literature (Barni, 2005; Barni, 2010).

**3.2 Chemical kinetics pattern in air**

In building our model and performing numerical simulations we have found more sensible to employ two different parameters, which are related to those describing the electric field ahead of the streamer. The first one, arises naturally when we have observed that since the electron density rises almost exponentially as the streamer develops, the total amount of radicals produced by electron impact does not depend strongly on the time length of the pulse or on the initial electron density but only on the actual value reached by the electron density at the end of the pulse. We have treated the value reached by the electron density ne at the end of the streamer passing as the relevant parameter of the simulation. As reference we consider a density of 1015 cm-3 (Braun, 1991; Kulikovsky, 1998). However we have varied ne in a broad range in order to investigate how it can affect the gas-phase chemical kinetics.

After imposing that the electron density reached in the streamer channel after its formation were the same, we also checked that the length of the electric field pulse does not make any significant difference in a broad range from 0.2 to 5 ns. This range is enough to cover the different streamer head velocities reported in DBDs (Raizer, 1991). As it could be guessed, most of the chemical kinetics is simply frozen during such fast times and the total production of species through electron impact reactions during the pulse depends on the electron reaction rates (influenced by the electric field amplitude) and on the electron density which is reached at the end of the pulse rather than on the pulse length.

As for the second parameters, the electric field strength Emax, we have found more sensible to employ the electron temperature instead. In fact it was simpler to express the electron rate constants, calculated from their cross-sections and evaluated assuming a Maxwellian energy distribution function for electrons, described by a single parameter, their temperature Te.

Then it was possible to collect reaction rates in form of simple functions of a single parameters with Arrhenius formulas, which are quite common in literature and databases of chemical kinetics and reaction rates. The use of electron temperature as a parameter does not prevent the capability to make comparison with other existing models, since it could be related, in a one to one correspondence, with experimental informations on the streamer electric field. Indeed electron temperature is trivially connected with the mean electron energy, which is determined by the local electric field in the Boltzmann equation (Raizer, 1991). This is sufficient to make straightforward a direct comparison between this simulation and other existing ones or experimental data. In the following we considered as reference an electron temperature value of 4 eV (Kulikovsky, 1998).

Under such simplified description of the streamer development, we could model the subsequent evolution of the gas phase in a standard way, using the continuity equations for each chemical species and solving a system of mono-dimensional first-order differential equations easily and quickly tackled by numerical integration (Riccardi, 2000). From a chemical engineering point of view, indeed it means that the model can be formulated as a well-mixed reactor (Benson, 1982). The gas-phase composition in the reactor is determined by the chemical reactions among the reactive species and the transport processes. The time evolution of the concentration of the different N species in the gas phase is determined by integrating each balance equation for the density nk of the kth species:

$$\frac{dn\_k}{dt} = \sum\_{\substack{i < j = 1 \ j = 1}}^N K(i + j \to k) \ n\_{\hat{i}} n\_{\hat{j}} - \sum\_{\substack{i < j = 1 \ i \neq j}}^N K(k + i \to j) \ n\_{\hat{k}} n\_{\hat{i}} + \ |D\_{\hat{k}} \nabla^2 n\_{\hat{k}}| \tag{1}$$

In building our model and performing numerical simulations we have found more sensible to employ two different parameters, which are related to those describing the electric field ahead of the streamer. The first one, arises naturally when we have observed that since the electron density rises almost exponentially as the streamer develops, the total amount of radicals produced by electron impact does not depend strongly on the time length of the pulse or on the initial electron density but only on the actual value reached by the electron density at the end of the pulse. We have treated the value reached by the electron density ne at the end of the streamer passing as the relevant parameter of the simulation. As reference we consider a density of 1015 cm-3 (Braun, 1991; Kulikovsky, 1998). However we have varied ne in a broad range in order to investigate how it can affect the gas-phase chemical kinetics. After imposing that the electron density reached in the streamer channel after its formation were the same, we also checked that the length of the electric field pulse does not make any significant difference in a broad range from 0.2 to 5 ns. This range is enough to cover the different streamer head velocities reported in DBDs (Raizer, 1991). As it could be guessed, most of the chemical kinetics is simply frozen during such fast times and the total production of species through electron impact reactions during the pulse depends on the electron reaction rates (influenced by the electric field amplitude) and on the electron

density which is reached at the end of the pulse rather than on the pulse length.

reference an electron temperature value of 4 eV (Kulikovsky, 1998).

integrating each balance equation for the density nk of the kth species:

1 , 1

As for the second parameters, the electric field strength Emax, we have found more sensible to employ the electron temperature instead. In fact it was simpler to express the electron rate constants, calculated from their cross-sections and evaluated assuming a Maxwellian energy distribution function for electrons, described by a single parameter, their temperature Te.

Then it was possible to collect reaction rates in form of simple functions of a single parameters with Arrhenius formulas, which are quite common in literature and databases of chemical kinetics and reaction rates. The use of electron temperature as a parameter does not prevent the capability to make comparison with other existing models, since it could be related, in a one to one correspondence, with experimental informations on the streamer electric field. Indeed electron temperature is trivially connected with the mean electron energy, which is determined by the local electric field in the Boltzmann equation (Raizer, 1991). This is sufficient to make straightforward a direct comparison between this simulation and other existing ones or experimental data. In the following we considered as

Under such simplified description of the streamer development, we could model the subsequent evolution of the gas phase in a standard way, using the continuity equations for each chemical species and solving a system of mono-dimensional first-order differential equations easily and quickly tackled by numerical integration (Riccardi, 2000). From a chemical engineering point of view, indeed it means that the model can be formulated as a well-mixed reactor (Benson, 1982). The gas-phase composition in the reactor is determined by the chemical reactions among the reactive species and the transport processes. The time evolution of the concentration of the different N species in the gas phase is determined by

*dn N N <sup>k</sup> Ki j k nn K k i j n n D n i j ki k k dt i j i j*

<sup>2</sup> () ()

(1)

where K are the reaction rates for the gas phase reactions and D is the diffusion coefficient. Due to the supposed uniformity, diffusion along the streamer axis is neglected. Adsorption on the dielectric barrier is included but, because of the much smaller extension of the streamer radius respect to the discharge gap, it turns out to be negligible. The choice of species included in the model was based on existing experimental information based on emission spectroscopy or mass spectroscopy. A sub-set of 8 neutral and 9 charged species was included in a preliminary study. Then the code was updated to include 20 neutral and excited and 12 ionic states.

As already discussed, to simulate DBD plasmas a cold plasma approximation is suitable. So reaction rates for ion as well as neutral species are evaluated at the set temperature, that is 300 K. The effect of larger temperatures is addressed below. Since thermal energy is much smaller than the molecular bond energy, the initial step in the decomposition is presumed to involve mainly electron-impact dissociation. Electron-impact cross-sections used in this study were obtained from literature collected in (Riccardi, 2001; Barni, 2005), apart from attachment rates which are taken from recent experiments. Rate constants have been evaluated assuming a Maxwellian energy distribution function for electron, described by a single parameter, their temperature Te. Rate constants for neutral gas-phase, charge exchange and ion recombination reactions have been taken from literature referenced in (Riccardi, 2001; Guerra 1997) or estimated based on approximating techniques (Benson, 1982). In this report we assume a substantial vibrational excitation of nitrogen molecules, which allows direct dissociation from atomic O radicals (Guerra, 1997). A total of 121 reactions have been taken into account in the simulation. The diffusion loss rates are estimated by using a formulation due to (Chantry, 1987) in cylindrical coordinates. The radial density profile of the streamer is approximated in term of the lowest order normal mode and is integrated away. This formulation allows substituting the Laplacian operator with an effective diffusion length and makes the system of balance equations only time dependent. Diffusion of air from the outside the streamer is also included as a source term. Ion diffusion loss rates was corrected by enforcing the plasma quasi-neutrality, following an approach proposed for low pressure electronegative plasmas (Lichtenberg, 2000). Balance equations containing all the processes described above form a system of coupled differential equations which has been integrated to obtain the time evolution by using an adaptive Runge-Kutta routine. Results from the calculations are discussed below and have been partially presented in literature (Barni, 2005; Barni, 2010).

#### **3.2 Chemical kinetics pattern in air**

As a first step we start discussing the temporal evolution of the discharge gas phase using a sort of reference values for the free parameters. That is electron temperature was 4 eV, electron density reached in the microdischarge channel is 2x1015 cm-3 in a time length of 2 ns and the streamer radius was 100 µm. At first, we neglect diffusion losses, in order to look directly the evolution pattern of the chemical kinetics alone. Neutral as well charged species densities are displayed as a function of time respectively in figures on a log-log scale.

As expected, the charged phase shows a very fast increase, almost exponential, in almost all species, as the streamer gets formed. Afterwards, as diffusion is not included, the system evolves under recombination alone, giving the overall decreasing trend as the inverse of time (that is linear in the log-log plot). The composition however displays some important

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 193

temperature). NO production is quite fast, the fastest among other molecules, but even so its density stays under 0.1% by the end of the streamer formation time. Approach to equilibrium is much slower. Between 10 and 100 µs oxygen atoms decline and molecules accumulate, with ozone already dominating. With the disappearance of O, all nitrogen oxides decline, leaving only ozone (besides air) in the gas phase. Nitrogen atoms are much less reactive and are removed only late, in tens ms, mainly through recombination to N2, as their linear decay in log-log plot clearly reveals. In real world, however such behaviour is difficult to observe. The relative narrow sizes of streamer channel and electrode gap make diffusion enter the game relatively early after only some ms. A more realistic modelization is therefore provided in the following pictures, with a more extended reaction dataset and where again neutral as well charged species densities are displayed as a function of time on a log-log scale. As for the neutral species, during the time length needed for streamer development mainly atomic N and O species are created, from electron impact-dissociation (the former more than the latter at such an electron temperature). NO production is quite fast and its density approaches 0.02% of N atoms one by the end of the pulse. Ozone formation is slower but almost catches up NO. On the other hand NO2, NO3 and N2O production is clearly much slower and become appreciable only on a larger scale time, even more so for N2O5. After the pulse end, in the temporal range from 2 to 20 ns, O and NO slowly increase till saturation is reached. N density reaches a value of 9x1015 cm-3, whereas O grows to 3x1015 cm-3 and NO to a value of 4x1013 cm-3. The dissociation degree is therefore rather limited, up to 5x10-4. Ozone, on the contrary, keeps increasing eventually reducing O and NO and peaking afterwards. Starting from 0.1 ms, the diffusion term becomes relevant and all species density starts to decrease. The gas phase is practically cleared out after a few ms. Ozone density exceeds atomic oxygen after about 10 µs, at a level of 1.3x1015 cm-3, and peaks after 30 µs, at a level of 2x1015 cm-3 about 100 ppm. The increase is triggered by the partial balancing between production from O + O2 and loss due to dissociative attachment

<sup>10</sup>-11 10-10 <sup>10</sup>-9 <sup>10</sup>-8 10-7 10-6 10-5 10-4 10-3 10-2 <sup>10</sup><sup>4</sup>

Fig. 3.2.3. Time evolution of the neutral species density in the gas-phase of a streamer in an air DBD. After a brief transient due to streamer formation, the neutral gas-phase reaches a definite composition in about 0.1 ms and it is then dissipated by diffusion in a few ms.

Time (s)

Chemical Composition of the Reactor Gas in an Air DBD

 N N2 O O2 O3 NO N2O NO2 NO3 N2O5

Density (cm^-3)

Fig. 3.2.1. Evolution of the charged density in the gas-phase of a streamer in an air DBD.

Fig. 3.2.2. Time evolution of the neutral species density in the gas-phase of a streamer in an air DBD. After a brief transient due to streamer formation, the neutral gas-phase reaches an equilibrium with the formation of ozone in a few ms.

changes. First atomic ions are removed immediately due to charge exchange reactions on nitrogen and oxygen molecules. Then electrons are lost, due to prevalence of attachment for cold electrons in the streamer channel. Their evolution reflects the shift between atomic and molecular oxygen negative ions. Later on, a sort of equilibrium in the positive ion composition is reached. It reflects the strength of interactions with main neutral molecules. The major ions are then negative atomic oxygen O and positive molecular oxygen ion O2 +, with traces of nitrogen and nitrogen oxides ions. As for the neutral species, during the time length needed for streamer development mainly atomic N and O species are created, from electron impact-dissociation (the former more than the latter at such an electron

 N+ N2+ N4+ O+ O2+ NO+ O- O2- Electron

 N N2 O O3 NO N2O NO2 O2

Chemical Composition of the Ion Species in an Air DBD

10-11 10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2 10-1 <sup>100</sup> <sup>104</sup>

10-11 10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2 10-1 <sup>104</sup>

Fig. 3.2.2. Time evolution of the neutral species density in the gas-phase of a streamer in an air DBD. After a brief transient due to streamer formation, the neutral gas-phase reaches an

changes. First atomic ions are removed immediately due to charge exchange reactions on nitrogen and oxygen molecules. Then electrons are lost, due to prevalence of attachment for cold electrons in the streamer channel. Their evolution reflects the shift between atomic and molecular oxygen negative ions. Later on, a sort of equilibrium in the positive ion composition is reached. It reflects the strength of interactions with main neutral molecules. The major ions are then negative atomic oxygen O- and positive molecular oxygen ion O2+, with traces of nitrogen and nitrogen oxides ions. As for the neutral species, during the time length needed for streamer development mainly atomic N and O species are created, from electron impact-dissociation (the former more than the latter at such an electron

Time (s)

Fig. 3.2.1. Evolution of the charged density in the gas-phase of a streamer in an air DBD.

Chemical Composition of the Reactor Gas in an Air DBD

Time (s)

equilibrium with the formation of ozone in a few ms.

Density (cm^-3)

Density (cm^-3)

temperature). NO production is quite fast, the fastest among other molecules, but even so its density stays under 0.1% by the end of the streamer formation time. Approach to equilibrium is much slower. Between 10 and 100 µs oxygen atoms decline and molecules accumulate, with ozone already dominating. With the disappearance of O, all nitrogen oxides decline, leaving only ozone (besides air) in the gas phase. Nitrogen atoms are much less reactive and are removed only late, in tens ms, mainly through recombination to N2, as their linear decay in log-log plot clearly reveals. In real world, however such behaviour is difficult to observe. The relative narrow sizes of streamer channel and electrode gap make diffusion enter the game relatively early after only some ms. A more realistic modelization is therefore provided in the following pictures, with a more extended reaction dataset and where again neutral as well charged species densities are displayed as a function of time on a log-log scale. As for the neutral species, during the time length needed for streamer development mainly atomic N and O species are created, from electron impact-dissociation (the former more than the latter at such an electron temperature). NO production is quite fast and its density approaches 0.02% of N atoms one by the end of the pulse. Ozone formation is slower but almost catches up NO. On the other hand NO2, NO3 and N2O production is clearly much slower and become appreciable only on a larger scale time, even more so for N2O5. After the pulse end, in the temporal range from 2 to 20 ns, O and NO slowly increase till saturation is reached. N density reaches a value of 9x1015 cm-3, whereas O grows to 3x1015 cm-3 and NO to a value of 4x1013 cm-3. The dissociation degree is therefore rather limited, up to 5x10-4. Ozone, on the contrary, keeps increasing eventually reducing O and NO and peaking afterwards. Starting from 0.1 ms, the diffusion term becomes relevant and all species density starts to decrease. The gas phase is practically cleared out after a few ms. Ozone density exceeds atomic oxygen after about 10 µs, at a level of 1.3x1015 cm-3, and peaks after 30 µs, at a level of 2x1015 cm-3 about 100 ppm. The increase is triggered by the partial balancing between production from O + O2 and loss due to dissociative attachment

Fig. 3.2.3. Time evolution of the neutral species density in the gas-phase of a streamer in an air DBD. After a brief transient due to streamer formation, the neutral gas-phase reaches a definite composition in about 0.1 ms and it is then dissipated by diffusion in a few ms.

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 195

the influence of the different streamer parameters we have performed several simulations changing their values in a broad range. The overall temporal evolution appears to be similar to the one reported in figures, even if some particulars are affected. In next figures the maximum density of a few neutral species is displayed as a function of the peak electron temperature. We made simulations varying Te from 2 to 4.5 eV with a step of 0.25 eV. We found that a minimum electron temperature of 2.5 eV is needed in order to allow electron density to rise during the streamer formation time. This can be understood since attachment prevails on ionization at small electron temperatures. The general behaviour shows that the maximum density decreases by increasing the electron temperature. This is due to the fact that ionization rate rises faster than the dissociation rate as the electron temperature is increased. Indeed, those simulations were performed keeping the peak electron density constant and since it means that the ionization grade stays constant then the overall dissociation grade should be diminished. The effects are more pronounced for species produced by secondary processes after the streamer formation, such as O3, NO2 and N2O than for atomic radicals such as N and O. We also checked that the length of the electric field pulse does not make any significant difference in a broad range from 0.2 to 5 ns. Indeed most of the chemical kinetics is simply frozen during such fast times and the total production of species through electron impact reactions during the pulse depends on the electron temperature (through the reaction rates) and on the electron density which is reached at the end of the pulse rather than on the pulse length. We have also performed several simulations changing the peak electron density ne from 1011 cm-3 to 1016 cm-3. Again, the temporal dynamics appears to be similar, with the same behavior of nitrogen oxides, ozone and atomic radicals and ions. In figure the maximum density of a few neutral species is displayed as a function of the peak electron density. All the maximum densities appear to increase increasing the electron density. The increase is approximately linear on the log-log scale of the figure. However the effect is much more pronounced for the heavier nitrogen oxides, thus increasing their concentration in the late gas phase, where ozone percentage in excess of 1000 ppm is reported. The strong sensitivity to such a badly known parameter of the discharge is clearly

2,5 3,0 3,5 4,0 4,5

Electron temperature (eV)

NO2

NO

NO3

(a)

N2 O5

O

O3

N2 O

Fig. 3.2.5. Dependence of the peak density reached by different neutral species on the electron temperature and maximum density in the gas-phase of a streamer in an air DBD.

<sup>1011</sup> <sup>1012</sup> <sup>1013</sup> <sup>1014</sup> <sup>1015</sup> <sup>1016</sup> 1017104

Peak electron density (cm-3)

NO3

NO

N

O

O3

N2 O5

Peak density (cm-3

)

N

(b)

N2 O

NO2

Peak density (cm-3

)

Fig. 3.2.4. Time evolution of the charged species density in the gas-phase of a streamer in an air DBD. After the transient producing an almost exponential increase of charged particles, electrons are quickly removed and an electron-less cold plasma is slowly dissipated due to recombination and then diffusion in about 10 ms.

and charge exchange reactions. Since charged particles are destroyed faster than neutrals due to the high rate of ion-ion recombination, ozone keeps on growing, until atomic O is finished and then declines due to diffusion. NO2 and NO3 appear to follow a similar trend, with density three order of magnitude lesser and a peak value after tens of µs, at a level of 2x1012 cm-3. N2O5 also has a similar trend, with much lower density and a slight delay. three order of magnitude lesser and a peak value after tens of µs, at a level of 2x1012 cm-3. On the other hand, N2O shows a partial equilibrium within a hundred of ns, then it grows to levels similar to the other nitrogen oxides. Its dynamics is closely controlled by fast dissociative attachment and charge exchange reactions in particular on atomic oxygen ions. So its concentration cannot rise until most of the electrons are removed. N2O density appears to take over NO2 after 100 µs and reaches a maximum of 7x1011 cm-3 soon later.

On the other hand, the temporal dynamics of charged species seems to be faster than that of neutrals. This is due to the higher rates of ion-ion recombination and of electron attachment. All the charged species increase almost exponentially during the streamer development time (2 ns) and then the total charge density begins to decrease because of recombination. In presence of oxygen, after the pulse ends, when the electron temperature decreases to less than 2 eV, electrons are quickly removed due to attachment processes. O- is the majority negative ion, except for a small time window between 10 and 300 ns. This happens because, while at Te=4 eV the dominant attachment process on oxygen is dissociative to O- , at set temperature prevails direct attachment to O2-. Then electrons disappear and soon after charge exchange reactions on oxygen atoms removes O2- and produces again the more stable O-. Positive charged species composition is mainly determined by charge exchange reactions on N2 and O2. They quickly remove atomic ions and leave O2 + as the majority ion, with a 10% of N2+ and a few % of N4 + and NO+. This situation is reached in about 100 ns. Later on positive charged species density shows the same course: a power-law decrease in time (that is linear on the log-log plot) dictated by the ion-ion recombination process, followed by an exponential decay due to ambipolar diffusion after 0.5 ms. In order to study

 N+ N2+ N4+ O+ O2+ NO+ O- O2- Electron

, at set

Chemical Composition of the Ion Species in an Air DBD

10-11 10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2 <sup>104</sup>

Fig. 3.2.4. Time evolution of the charged species density in the gas-phase of a streamer in an air DBD. After the transient producing an almost exponential increase of charged particles, electrons are quickly removed and an electron-less cold plasma is slowly dissipated due to

and charge exchange reactions. Since charged particles are destroyed faster than neutrals due to the high rate of ion-ion recombination, ozone keeps on growing, until atomic O is finished and then declines due to diffusion. NO2 and NO3 appear to follow a similar trend, with density three order of magnitude lesser and a peak value after tens of µs, at a level of 2x1012 cm-3. N2O5 also has a similar trend, with much lower density and a slight delay. three order of magnitude lesser and a peak value after tens of µs, at a level of 2x1012 cm-3. On the other hand, N2O shows a partial equilibrium within a hundred of ns, then it grows to levels similar to the other nitrogen oxides. Its dynamics is closely controlled by fast dissociative attachment and charge exchange reactions in particular on atomic oxygen ions. So its concentration cannot rise until most of the electrons are removed. N2O density appears to

On the other hand, the temporal dynamics of charged species seems to be faster than that of neutrals. This is due to the higher rates of ion-ion recombination and of electron attachment. All the charged species increase almost exponentially during the streamer development time (2 ns) and then the total charge density begins to decrease because of recombination. In presence of oxygen, after the pulse ends, when the electron temperature decreases to less than 2 eV, electrons are quickly removed due to attachment processes. O- is the majority negative ion, except for a small time window between 10 and 300 ns. This happens because, while at Te=4 eV the dominant attachment process on oxygen is dissociative to O-

temperature prevails direct attachment to O2-. Then electrons disappear and soon after charge exchange reactions on oxygen atoms removes O2- and produces again the more stable O-. Positive charged species composition is mainly determined by charge exchange reactions on N2 and O2. They quickly remove atomic ions and leave O2+ as the majority ion,

Later on positive charged species density shows the same course: a power-law decrease in time (that is linear on the log-log plot) dictated by the ion-ion recombination process, followed by an exponential decay due to ambipolar diffusion after 0.5 ms. In order to study

+ and NO+. This situation is reached in about 100 ns.

take over NO2 after 100 µs and reaches a maximum of 7x1011 cm-3 soon later.

Time (s)

recombination and then diffusion in about 10 ms.

with a 10% of N2+ and a few % of N4

Density (cm^-3)

the influence of the different streamer parameters we have performed several simulations changing their values in a broad range. The overall temporal evolution appears to be similar to the one reported in figures, even if some particulars are affected. In next figures the maximum density of a few neutral species is displayed as a function of the peak electron temperature. We made simulations varying Te from 2 to 4.5 eV with a step of 0.25 eV. We found that a minimum electron temperature of 2.5 eV is needed in order to allow electron density to rise during the streamer formation time. This can be understood since attachment prevails on ionization at small electron temperatures. The general behaviour shows that the maximum density decreases by increasing the electron temperature. This is due to the fact that ionization rate rises faster than the dissociation rate as the electron temperature is increased. Indeed, those simulations were performed keeping the peak electron density constant and since it means that the ionization grade stays constant then the overall dissociation grade should be diminished. The effects are more pronounced for species produced by secondary processes after the streamer formation, such as O3, NO2 and N2O than for atomic radicals such as N and O. We also checked that the length of the electric field pulse does not make any significant difference in a broad range from 0.2 to 5 ns. Indeed most of the chemical kinetics is simply frozen during such fast times and the total production of species through electron impact reactions during the pulse depends on the electron temperature (through the reaction rates) and on the electron density which is reached at the end of the pulse rather than on the pulse length. We have also performed several simulations changing the peak electron density ne from 1011 cm-3 to 1016 cm-3. Again, the temporal dynamics appears to be similar, with the same behavior of nitrogen oxides, ozone and atomic radicals and ions. In figure the maximum density of a few neutral species is displayed as a function of the peak electron density. All the maximum densities appear to increase increasing the electron density. The increase is approximately linear on the log-log scale of the figure. However the effect is much more pronounced for the heavier nitrogen oxides, thus increasing their concentration in the late gas phase, where ozone percentage in excess of 1000 ppm is reported. The strong sensitivity to such a badly known parameter of the discharge is clearly

Fig. 3.2.5. Dependence of the peak density reached by different neutral species on the electron temperature and maximum density in the gas-phase of a streamer in an air DBD.

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 197

N2 O2 NO O3

NO2 N2 O

N2 O2

NO3

NO

O

)

)

10-11 10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2 10-1 100

N2 O5

O

NO3

N2 O

N

N2 O5

A similar trend is reported for N2O which is however better oxidised.

oxidation of biological matter for sterilization purposes (Moreau, 2008).

N

time (s)

10-11 10-10 10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2 10-1 100

time (s)

Fig. 3.3.2. Time evolution of the neutral species density in the gas-phase of a streamer. NO2 is removed to mainly ozone, the transition happening within a few tens of microseconds.

too, but the initial nitrogen oxide is not completely removed and about 2% is not disposed.

The effect of the repetition of streamer discharges during subsequent cycles of a.c. voltage was studied and displayed in Figure 3.3.4 at 10 kHz. A small number of discharges is needed to remove NO under a few %. To reach an equilibrium in the gas-phase, here producing a constant flux of ozone, consuming mainly atmospheric oxygen, requires a larger number of cycles, in the order of thousands, yet in times less than seconds. In this case substantial amounts of radicals in the form of atomic nitrogen and oxygen are present and could be delivered to materials exposed to the plasma phase, for instance to achieve

Fig. 3.3.1. Time evolution of the neutral species density in the gas-phase of a streamer. NO is removed to mainly ozone, the transition happening within a few tens of microseconds.

<sup>O</sup> NO <sup>3</sup> <sup>2</sup>

disappointing from the point of view of driving the optimization process with chemical kinetics modeling, but it offers also a clue for an indirect insight of the streamer properties from this strong link between ne and the composition of radical fluxes outside the streamer. We have also investigated the effects due to the streamer radius Rs. Such a parameter affects mainly the diffusion process. We have varied Rs between 0.05 to 1 mm. The smaller Rs is, the faster diffusion transports radicals and ions outside the streamer. If it happens fast enough it interferes with secondary chemical processes leading to transformation of atomic radicals in O3, NO2 and N2O. A slower diffusion gives time to the system to approach the chemical equilibrium in composition. However equilibrium is approached in times exceeding a few hundreds milliseconds, which is definitely too long respect to any reasonable diffusion time.

#### **3.3 Reactivity of different nitrogen oxides**

Plasma remediation is an efficient and promising technology to destroy toxic and greenhouse gases (Manheimer, 1997). These techniques are potentially lower in cost and more efficient than conventional thermal incineration. In these plasmas, the majority of electrical energy is expanded in heating electrons, rather than heating gas, as discussed previously. The understanding of the process and the optimization of the efficiency benefit largely from computational studies of the dynamics of the plasma gas-phase. In particular the air chemical kinetics of high-pressure dielectric barrier discharges. Such systems are considered in the context of volatile organic compound (VOC) and NOx remediation. The air plasma kinetic processes are discussed with the goal of providing insight for optimizing efficiencies. Indeed the most extensively investigated applications of this type are used to treat dilute concentrations of toxic molecules in airstreams. Potential advantages of the highpressure, nonthermal plasma approach include the highly energy efficient selectivity provided by the plasma chemistry and its capability for minimizing secondary-waste production. Atmospheric pressure operation is preferred for high-throughput waste processing (Shoyama, 2007). Indeed the removal of NO from environment is a relevant problem concerning pollution control (Amoroso, 2008). NO stays in the low atmosphere for several days and enters cycles producing ozone, causing smog. Nitrogen oxides are mainly produced as by-product in combustion processes in free air. For instance in diesel engine exhaust NO can reach 200-500 ppm and NO2 50 ppm (Srinivasan, 2007). So we applied our simulations to air contaminated by similar level of nitrogen oxides. In Figures 3.3.1 and 3.3.2 we display the results of a simulation, with the same standard parameters but starting with 500 ppm of NO, that is a density of 1.3x1016 cm-3. No diffusion of neutrals is considered here, to let the system reach the chemical equilibrium, which however happens only at late times, both considering discharge repetition times (several kHz generally) and diffusion times for the considered streamer radius (a few ms). After the end of the discharge at 2 ns, NO is completely removed, oxidated initially by O atoms and then, when atomic radicals disappear, consumed in reactions with ozone and NO2. In the final state it is substituted mainly by ozone (about 1000 ppm).

A substantial amount of NO2 (192 ppm) and N2O (11 ppm) is produced too, while other oxides are not formed with an appreciable concentration. The transition from oxygen atoms to ozone happens after about 10 μs, while NO disappears only after 10 ms. More details are discussed in literature (Barni, 2010). Starting from NO2 as initial contaminant produces a different kinetic pattern. Oxidation and substitution with ozone in the final state is achieved

disappointing from the point of view of driving the optimization process with chemical kinetics modeling, but it offers also a clue for an indirect insight of the streamer properties from this strong link between ne and the composition of radical fluxes outside the streamer. We have also investigated the effects due to the streamer radius Rs. Such a parameter affects mainly the diffusion process. We have varied Rs between 0.05 to 1 mm. The smaller Rs is, the faster diffusion transports radicals and ions outside the streamer. If it happens fast enough it interferes with secondary chemical processes leading to transformation of atomic radicals in O3, NO2 and N2O. A slower diffusion gives time to the system to approach the chemical equilibrium in composition. However equilibrium is approached in times exceeding a few hundreds milliseconds, which is definitely too long respect to any

Plasma remediation is an efficient and promising technology to destroy toxic and greenhouse gases (Manheimer, 1997). These techniques are potentially lower in cost and more efficient than conventional thermal incineration. In these plasmas, the majority of electrical energy is expanded in heating electrons, rather than heating gas, as discussed previously. The understanding of the process and the optimization of the efficiency benefit largely from computational studies of the dynamics of the plasma gas-phase. In particular the air chemical kinetics of high-pressure dielectric barrier discharges. Such systems are considered in the context of volatile organic compound (VOC) and NOx remediation. The air plasma kinetic processes are discussed with the goal of providing insight for optimizing efficiencies. Indeed the most extensively investigated applications of this type are used to treat dilute concentrations of toxic molecules in airstreams. Potential advantages of the highpressure, nonthermal plasma approach include the highly energy efficient selectivity provided by the plasma chemistry and its capability for minimizing secondary-waste production. Atmospheric pressure operation is preferred for high-throughput waste processing (Shoyama, 2007). Indeed the removal of NO from environment is a relevant problem concerning pollution control (Amoroso, 2008). NO stays in the low atmosphere for several days and enters cycles producing ozone, causing smog. Nitrogen oxides are mainly produced as by-product in combustion processes in free air. For instance in diesel engine exhaust NO can reach 200-500 ppm and NO2 50 ppm (Srinivasan, 2007). So we applied our simulations to air contaminated by similar level of nitrogen oxides. In Figures 3.3.1 and 3.3.2 we display the results of a simulation, with the same standard parameters but starting with 500 ppm of NO, that is a density of 1.3x1016 cm-3. No diffusion of neutrals is considered here, to let the system reach the chemical equilibrium, which however happens only at late times, both considering discharge repetition times (several kHz generally) and diffusion times for the considered streamer radius (a few ms). After the end of the discharge at 2 ns, NO is completely removed, oxidated initially by O atoms and then, when atomic radicals disappear, consumed in reactions with ozone and NO2. In the final state it is substituted

A substantial amount of NO2 (192 ppm) and N2O (11 ppm) is produced too, while other oxides are not formed with an appreciable concentration. The transition from oxygen atoms to ozone happens after about 10 μs, while NO disappears only after 10 ms. More details are discussed in literature (Barni, 2010). Starting from NO2 as initial contaminant produces a different kinetic pattern. Oxidation and substitution with ozone in the final state is achieved

reasonable diffusion time.

**3.3 Reactivity of different nitrogen oxides**

mainly by ozone (about 1000 ppm).

Fig. 3.3.1. Time evolution of the neutral species density in the gas-phase of a streamer. NO is removed to mainly ozone, the transition happening within a few tens of microseconds.

Fig. 3.3.2. Time evolution of the neutral species density in the gas-phase of a streamer. NO2 is removed to mainly ozone, the transition happening within a few tens of microseconds.

too, but the initial nitrogen oxide is not completely removed and about 2% is not disposed. A similar trend is reported for N2O which is however better oxidised.

The effect of the repetition of streamer discharges during subsequent cycles of a.c. voltage was studied and displayed in Figure 3.3.4 at 10 kHz. A small number of discharges is needed to remove NO under a few %. To reach an equilibrium in the gas-phase, here producing a constant flux of ozone, consuming mainly atmospheric oxygen, requires a larger number of cycles, in the order of thousands, yet in times less than seconds. In this case substantial amounts of radicals in the form of atomic nitrogen and oxygen are present and could be delivered to materials exposed to the plasma phase, for instance to achieve oxidation of biological matter for sterilization purposes (Moreau, 2008).

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 199

temperature within the device to be simulated, the software has been changed adding the temperature dependence in several reactions involving neutral species, rates have been collected through different public databases as coefficients for the Arrhenius equation. The main results was that at low temperature it has been observed a cut down in nitrogen monoxide concentrations, which has been replaced by ozone as principal product. At high temperature nitrogen monoxide remains unperturbed through the discharge while ozone disappears. Thus gas phase temperature, and then the control of a strictly cold plasma

Fig. 3.4.1. Final composition of the neutral species density in the gas-phase of a streamer. NO is removed to ozone only at low temperatures, the transition happening at about 600 K.

Then we consider application of the modeling to a steady plasma conditions, like the one in a corona discharge. In order to study direct chemical kinetics effect we have chosen an ideal setup, under a well-mixed reactor approximation, and with an uniform plasma background. Apart from geometry (here 3 mm corona electrode is discussed), we have parametrised our results in terms of two physical values, the effective electron density and temperature. As discussed previously, the former affects mostly the velocity of the chemical kinetics evolution, whereas the latter affects the relative weight of the different electron impact. In order to study the influence of the different parameters we have performed several simulations changing their values in a broad range. The overall temporal evolution appears to be similar to the one reported in figure, even if some particulars are affected. The corona discharge produce mainly NO reaching a steady state after a few hundreds of microseconds. There is a transition from the initial phases of the discharge, where N and NO are the major neutrals. All neutral molecules are accumulating in the gas-phase until 10 microseconds. Then nitrogen atoms start to decline, followed a little after by ozone. The final composition is then reached, with about 90% of NO as final product. In next figure the maximum density of a few neutral species is displayed as a function of the peak electron temperature and

condition is crucial for efficient NOx remediation.

**4. Chemical kinetics in air steady plasmas**

Fig. 3.3.3. Time evolution of the neutral species density in the gas-phase of a streamer. N2O is removed to mainly ozone, the transition happening within a few tens of microseconds.

Fig. 3.3.4. Evolution of the neutral species density in the gas-phase of a streamer under subsequent repetitions at a frequency of . N2O is removed to mainly ozone, the transition happening within a few tens of microseconds.

#### **3.4 Effect of gas temperature**

Reaction rates for ion as well as neutral species depends on their temperature, which could depart from the set, that is 300 K. Many neutral gas reactions are quite affected by the gas temperature, so it is interesting to use the theoretical framework we have developed to study what is the impact of a hotter gas phase. In any case, since thermal energy is much smaller than the molecular bond energy, the initial step in the decomposition is presumed to involve again mostly electron-impact dissociation. An idea of the impact the gas temperature has on the chemical kinetics is reported in figure.  Time evolution of the neutral species density in the gas-phase of a streamer. N2O is removed to mainly ozone, the transition happening within a few tens of microseconds.

100 101 102 103

repetition

Reaction rates for ion as well as neutral species depends on their temperature, which could depart from the set, that is 300 K. Many neutral gas reactions are quite affected by the gas temperature, so it is interesting to use the theoretical framework we have developed to study what is the impact of a hotter gas phase. In any case, since thermal energy is much smaller than the molecular bond energy, the initial step in the decomposition is presumed to involve again mostly electron-impact dissociation. An idea of the impact the gas temperature has on the chemical kinetics is reported in figure. Taking into account the air

Fig. 3.3.4. Evolution of the neutral species density in the gas-phase of a streamer under subsequent repetitions at a frequency of . N2O is removed to mainly ozone, the transition temperature within the device to be simulated, the software has been changed adding the temperature dependence in several reactions involving neutral species, rates have been collected through different public databases as coefficients for the Arrhenius equation. The main results was that at low temperature it has been observed a cut down in nitrogen monoxide concentrations, which has been replaced by ozone as principal product. At high temperature nitrogen monoxide remains unperturbed through the discharge while ozone disappears. Thus gas phase temperature, and then the control of a strictly cold plasma condition is crucial for efficient NOx remediation.

Fig. 3.4.1. Final composition of the neutral species density in the gas-phase of a streamer. NO is removed to ozone only at low temperatures, the transition happening at about 600 K.

#### **4. Chemical kinetics in air steady plasmas**

Then we consider application of the modeling to a steady plasma conditions, like the one in a corona discharge. In order to study direct chemical kinetics effect we have chosen an ideal setup, under a well-mixed reactor approximation, and with an uniform plasma background.

Apart from geometry (here 3 mm corona electrode is discussed), we have parametrised our results in terms of two physical values, the effective electron density and temperature. As discussed previously, the former affects mostly the velocity of the chemical kinetics evolution, whereas the latter affects the relative weight of the different electron impact. In order to study the influence of the different parameters we have performed several simulations changing their values in a broad range. The overall temporal evolution appears to be similar to the one reported in figure, even if some particulars are affected. The corona discharge produce mainly NO reaching a steady state after a few hundreds of microseconds. There is a transition from the initial phases of the discharge, where N and NO are the major neutrals. All neutral molecules are accumulating in the gas-phase until 10 microseconds. Then nitrogen atoms start to decline, followed a little after by ozone. The final composition is then reached, with about 90% of NO as final product. In next figure the maximum density of a few neutral species is displayed as a function of the peak electron temperature and

Chemical Kinetics in Air Plasmas at Atmospheric Pressure 201

saturation at the higher end of the range. The effect is much more pronounced the heavier nitrogen oxides, thus increasing their concentration in the late gas phase. Ozone and

A chemical kinetics modeling of atmospheric pressure air discharges in the streamer regime of DBDs as well as in corona steady state plasmas has been presented. The model was used to evaluate the prospect of nitrogen oxides removal from the gas phase. The results obtained show that the chemical system could remove NO and to some extent also other nitrogen oxides up to a certain maximum initial concentration. Under realistic conditions for NO this could be as high as 800 ppm. The system produces mainly ozone through secondary reactions in the gas-phase. This theoretical study shows that important and significant aspects determined by the chemical kinetics effects could be investigated through simulations. Then they could prove worthy in order to study the applications of plasmas at high pressure in streamer regime (for example DBDs) and their capabilities as gas

The authors are deeply grateful to their collaborators and students that over several years have partecipated in their research studies, in the laboratory activities and in discussion about plasma science during working hours and not, contributing to the success of the PlasmaPrometeo Center of Excellence. We could enlist M.Piselli, S.Zanini, P.Esena, I.Nardi, E.Grimoldi, F.Fumagalli, I.Biganzoli just to mention a few involved in the topics of this

Amoroso, A.; Beine, H.J.; Esposito, G.; Perrino, C.; Catrambone, M. & Allegrini, I. (2008).

Barni, R. & Riccardi, C. (2010). Perspective of NOx removal from numerical simulation of

Barni, R.; Esena, P. & Riccardi, C. (2005). Chemical kinetics simulation for atmospheric

Barni, R.; Esena, P. & Riccardi, C. (2005). Chemical kinetics simulations of an atmospheric pressure plasma device, *Surface & Coatings Technology,* Vol. 200/1-4, pp. 924-927

Esena, P.; Riccardi, C.; Zanini, S.; Poletti, G. & Orsini, F. (2005). *Surf. & Coatings Techn*, Vol.

Braun, D.; Gibalov, V.I. & Pietsch, G.J. (1991). *Journal of Physics D*, Vol. 24, pp. 564

Guerra, V. & Louriero, J. (1997). *Plasma Sources Science & Technology*, Vol. 6, pp. 373

Kogelschatz, U. (2003). *Plasma Chemistry and Plasma Processing*, Vol.23, pp.1

non-thermal atmospheric pressure plasma chemical kinetics, *High Temperature*

pressure air plasmas in a streamer regime, *Journal of Applied. Physics*, Vol. 97, pp.

*Water Air and Soil Pollution*, Vol. 188, pp. 81

Benson, S.W. (1982). *Thermochemical Kinetics*, Wiley, New York US

Chantry, P.J. (1987). *Journal of Applied. Physics*, Vol. 62, pp. 1141

*Material Processes*, Vol. 14, pp. 205-210

nitrogen atoms decrease, becoming comparable to atomic oxygen.

**5. Conclusions**

processing units.

chapter.

**7. References**

073301.1-7

1/4, pp. 664

**6. Acknowledgment**

density. We made simulations varying Te from 0.5 to 2 eV with a step of 0.25 eV. We found that above 1 eV the density of nitrogen oxides produced in the gas-phase becomes independent from the temperature. Ozone and atomic radicals are minority showing a broad minimum around 1 eV. We have also performed several simulations changing the peak electron density ne from 108 cm-3 to 1012 cm-3. Again, the temporal dynamics appear to be similar, with the same behavior of nitrogen oxides, ozone and atomic radicals and ions.

Fig. 4.1. Time evolution of the neutral species density in the gas-phase of an air corona. The neutral gas-phase reaches an equilibrium in about one second, the transition happening between 1 and 100 ms.

Fig. 4.2. Dependence of the peak density reached by different neutral species on the electron temperature and density in the gas-phase of an air corona.

All the nitrogen oxides densities appear to increase increasing the electron density. The increase is approximately linear on the log-log scale of the figure, showing a trend to saturation at the higher end of the range. The effect is much more pronounced the heavier nitrogen oxides, thus increasing their concentration in the late gas phase. Ozone and nitrogen atoms decrease, becoming comparable to atomic oxygen.

#### **5. Conclusions**

200 Chemical Kinetics

density. We made simulations varying Te from 0.5 to 2 eV with a step of 0.25 eV. We found that above 1 eV the density of nitrogen oxides produced in the gas-phase becomes independent from the temperature. Ozone and atomic radicals are minority showing a broad minimum around 1 eV. We have also performed several simulations changing the peak electron density ne from 108 cm-3 to 1012 cm-3. Again, the temporal dynamics appear to be similar, with the same behavior of nitrogen oxides, ozone and atomic radicals and ions.

Chemical Composition of the Reactor Gas in an Air corona

10-6 1x10-5 1x10-4 10-3 10-2 10-1 10<sup>0</sup> 10<sup>1</sup>

Time (ms)

(b)

NO

N2

O3

<sup>10</sup><sup>9</sup> <sup>10</sup><sup>10</sup> <sup>10</sup><sup>11</sup> <sup>10</sup><sup>12</sup> <sup>10</sup><sup>11</sup>

O

Peak electron density (cm-3

<sup>O</sup> NO2

N

Peak density (cm-3

)

)

Fig. 4.1. Time evolution of the neutral species density in the gas-phase of an air corona. The neutral gas-phase reaches an equilibrium in about one second, the transition happening

Fig. 4.2. Dependence of the peak density reached by different neutral species on the electron

All the nitrogen oxides densities appear to increase increasing the electron density. The increase is approximately linear on the log-log scale of the figure, showing a trend to

between 1 and 100 ms.

(a) NO

N2 O

Peak density (cm-3

)

 N N2 O O3 NO N2O NO2 O2

0,5 1,0 1,5 2,0

Electron temperature (eV)

temperature and density in the gas-phase of an air corona.

O3

N

NO2

O

Density (cm^-3)

A chemical kinetics modeling of atmospheric pressure air discharges in the streamer regime of DBDs as well as in corona steady state plasmas has been presented. The model was used to evaluate the prospect of nitrogen oxides removal from the gas phase. The results obtained show that the chemical system could remove NO and to some extent also other nitrogen oxides up to a certain maximum initial concentration. Under realistic conditions for NO this could be as high as 800 ppm. The system produces mainly ozone through secondary reactions in the gas-phase. This theoretical study shows that important and significant aspects determined by the chemical kinetics effects could be investigated through simulations. Then they could prove worthy in order to study the applications of plasmas at high pressure in streamer regime (for example DBDs) and their capabilities as gas processing units.

#### **6. Acknowledgment**

The authors are deeply grateful to their collaborators and students that over several years have partecipated in their research studies, in the laboratory activities and in discussion about plasma science during working hours and not, contributing to the success of the PlasmaPrometeo Center of Excellence. We could enlist M.Piselli, S.Zanini, P.Esena, I.Nardi, E.Grimoldi, F.Fumagalli, I.Biganzoli just to mention a few involved in the topics of this chapter.

#### **7. References**

**10**

*Canada*

**The Chemical Kinetics of**

David M. Holloway

 **Shape Determination in Plants**

*Mathematics Department, British Columbia Institute of Technology,*

Plants are integral to our lives, providing food, shelter and the air we breathe. The shapes that plants take are central to their functionality, tailoring each for its particular place in the ecosystem. Given the relatively large and static forms of plants, it may not be immediately apparent that chemical kinetics is involved in, for example, distinguishing the form of a spruce tree from that of a fern. But plants share the common feature that their shapes are continuously being generated, and this largely occurs in localized regions of cell division and expansion, such as the shoot and root apical meristems at either end of a plant's main axis; these regions remain essentially embryonic throughout the life cycle. The final regular structure of a plant, such as the arrangement of leaves along the main stalk, may seem to follow an overall spatial template; but in reality the spatial patterning is occurring at relatively short range, and it is the temporal unfolding of this small scale patterning which generates the plant's form. A key part of understanding plant morphogenesis, or shape generation, therefore, is to understand how the molecular determinants of cell type, cell division and cell expansion are localized to and patterned within the actively growing regions. At this scale, transport processes such as diffusion and convection are obvious components of localization, for moving molecules to the correct places; but the reaction kinetics for molecular creation, destruction and interaction are also critical to maintaining

It was Turing (1952) who first combined the ideas of reaction kinetics with diffusion into a theory of spatial pattern formation – how to establish, and maintain, molecules at specific concentrations in specific locations. The Turing mechanism involves two mutually interacting (activating and inhibiting) intermediates (which he termed 'morphogens'), which diffuse at different rates. While his paper chiefly introduced the kinetics required for pattern formation from homogeneous initial conditions (in particular, finding critical values of rate constants for this to occur), Turing was clearly thinking of the implications in embryonic tissues, as evidenced by his title, "The chemical theory of morphogenesis". In the 1960's, Prigogine and collaborators presented the first specific nonlinear chemical mechanism which could form stable spatial concentration waves (the 'Brusselator'; Prigogine & Lefever, 1968), as well as providing a thermodynamic basis for spontaneous pattern formation (see Nicolis & Prigogine, 1977). In the 1970's, reaction-diffusion (RD) mechanisms were introduced into biology by Meinhardt (see Meinhardt, 1982), Harrison (see Harrison, 1993),

the molecular identity and the size regulation of the active regions.

**1. Introduction**


5. Conclusion

Within the framework of the theory of row-column determinants previously introduced by the author, we get determinantal representations (analogs of Cramer's rule) of partial solutions to the system of two-sided quaternion matrix equations A1XB1=C1, A2XB2=C2, and its special cases with 1 and 2 one-sided matrix equations. We use previously obtained by the author determinantal representations of the Moore-Penrose inverse. Note to give determinantal representations for all above matrix systems over the complex field, it is obviously needed to substitute all row and column determinants by usual determinants.

[6] Navarra A, Odell PL, Young DM. A representation of the general common solution to the matrix equations A1XB1=C1 and A2XB2=C2 with applications. Computers & Mathematics

Cramer's Rules for the System of Two-Sided Matrix Equations and of Its Special Cases

http://dx.doi.org/10.5772/intechopen.74105

19

[7] Wang QW. The general solution to a system of real quaternion matrix equations. Computers & Mathematics with Applications. 2005;49:665-675. DOI: 10.1016/j.camwa.2004.12.00 [8] Kyrchei I. Cramer's rule for generalized inverse solutions. In: Kyrchei I, editor. Advances

[9] Aslaksen H. Quaternionic determinants. Mathematical Intelligence. 1996;18(3):57-65. DOI:

[10] Cohen N, De Leo S. The quaternionic determinant. Electronic Journal of Linear Algebra.

[11] Dieudonne J. Les determinants sur un corps non-commutatif. Bulletin de la Société

[13] Cayley A. On certain results relating to quaternions. Philosophical Magazine. 1845;26:141- 145. Reprinted in The collected mathematical papers. Cambridge Univ. Press. 1889;1:123-126

[14] Moore EH. On the determinant of an hermitian matrix of quaternionic elements. Bulletin

[15] Dyson FJ. Quaternion determinants. Helvetica Physica Acta. 1972;45:289-302. DOI: 10.5169/

[16] Chen L. Definition of determinant and Cramer solution over the quaternion field. Acta

[17] Kyrchei I. Cramer's rule for quaternion systems of linear equations. Fundamentalnaya i

[18] Kyrchei I. The theory of the column and row determinants in a quaternion linear algebra. In: Baswell AR, editor. Advances in Mathematics Research 15. New York: Nova Sci. Publ;

[19] Kyrchei I. Determinantal representations of the Moore-Penrose inverse over the quaternion skew field and corresponding Cramer's rules. Linear Multilinear Algebra. 2011;59(4):

[20] Kyrchei I. Explicit representation formulas for the minimum norm least squares solutions of some quaternion matrix equations. Linear Algebra and its Applications. 2013;438(1):

[21] Kyrchei I. Determinantal representations of the Drazin inverse over the quaternion skew field with applications to some matrix equations. Applied Mathematics and Computation.

with Applications. 2001;41:929-935. DOI: 10.1016/S0898-1221(00)00330-8

in Linear Algebra Research. New York: Nova Sci. Publ; 2015. pp. 79-132

[12] Study E. Zur theorie der linearen gleichungen. Acta Mathematica. 1920;42:1-61

10.1007/BF03024312

seals-114,385

2012. pp. 301-359

2000;7:100-111. DOI: 10.13001/1081-3810.1050

of the American Mathematical Society. 1922;28:161-162

Mathematica Sinica. 1991;7:171-180. DOI: 10.1007/BF02633946

Mathématique de France. 1943;71:27-45

Prikladnaya Matematika. 2007;13(4):67-94

413-431. DOI: 10.1080/03081081003586860

136-152. DOI: 10.1016/j.laa.2012.07.049

2014;238:193-207. DOI: 10.1016/j.amc.2014.03.125


**Rijin Cheng 1,2,3, Hua Zhang 1,3,\* and Hongwei Ni 1,\***

Received: 12 September 2019; Accepted: 13 October 2019; Published: 16 October 2019

**Abstract:** In most cases, arsenic is an unfavorable element in metallurgical processes. The mechanism of arsenic removal was investigated through roasting experiments performed on arsenopyrite-bearing iron ore. Thermodynamic calculation of arsenic recovery was carried out by FactSage 7.0 software (Thermfact/CRCT, Montreal, Canada; GTT-Technologies, Ahern, Germany). Moreover, the arsenic residues in dust ash were recovered by roasting dust ash in a reducing atmosphere. Furthermore, the corresponding chemical properties of the roasted ore and dust ash were determined by X-ray diffraction, inductively coupled plasma atomic emission spectrometry, and scanning electron microscopy, coupled with energy-dispersive X-ray spectroscopy. The experimental results revealed that the arsenic in arsenopyrite-bearing iron ore can be removed in the form of As2O3(g) in an air or nitrogen atmosphere by a roasting method. The efficiency of arsenic removal through roasting in air was found to be less than that in nitrogen atmosphere. The method of roasting in a reducing atmosphere is feasible for arsenic recovery from dust ash. When the carbon mass ratio in dust ash is 1.83%, the arsenic removal products is almost volatilized and recovered in the form of As2O3(g).

**Keywords:** arsenopyrite; arsenic removal; mechanism; roasting; arsenate; dust ash; arsenic recovery