well as a Fast Equilibrium to Non-Equilibrium TItansition (which we call FENE-. TRe) between them [4,5]. During the FENETRe an ionization instability creates.
V o l . 89 (1996)
ACTA
PHYSIC A POLONICA
A
No. 5-6
S P E C T R A L A N D E L E C T R I C A L DIAGNOSTICS OF GLIDING A R C A . CZERNICHOWSKI, H . NASSAR, A . RANAIVOSOLOARIMANANA
G R E M I , University of Orleans, 45067 Orleans cedex 2, France A . A . FRIDMAN
Kurchatov Institute, 123182 Moscow, Russia M . SlMEK
Institute of Plasma Physics, 182 11 Prague, Czech Republic K . MUSIOL, E . PAWELEC
Jagiellonian University, Reymonta 4, 30-059 Cracow, Poland AND L . DlTTRICHOVA
Technical University, 616 69 Brno, Czech Republic (Received January 9, 1996)
Using spectroscopic and electric measurements, vibrational and rotational molecular gas temperatures as well as free electron temperature and concentration were determined in different regions of a time-periodical type, atmospheric pressure non-equilibrium low current gliding arc. It was shown that this discharge includes an initial quasi-equilibrium zone, with the quasi-equilibrium temperature of 4 to 6 kK, and the non-equilibrium zone with the electron temperature about 10 kK, the vibrational temperature about 3 kK, rotational and translational temperatures from 1 to 1.5 kK. The transition between two mentioned zones coincides with the phenomenon of the arc "length explosion" already observed in moderate-current gliding arc. PACS numbers: 33.10.Jz, 52.70.Kz, 52.25.Kn 1. I n t r o d u c t i o n In the previous papers [1,2] we described some properties of the so-called gliding arc (or GlidArc), a quasi periodic electrical discharge, which has been used for several applications in chemistry and environment protection [3]. This powerful (up to 40 k W ) discharge can be easily obtained between at least two divergent electrodes, with a gas flowing between them at pressures ranging from (595)
A. Czernichowski et al.
596
soft vacuum to several atmospheres. It was also explained how this discharge can produce a relatively cold non-equilibrium plasma with a complicated time-space structure including quasi-equilibrium (hot) and non-equilibrium (cold) periods, as well as a Fast Equilibrium to Non-Equilibrium TItansition (which we call F E N E TRe) between them [4,5]. During the F E N E T R e an ionization instability creates conditions for existence of a relatively cold plasma and at this moment, the phenomenon of the arc "length explosion" takes place. More detailed theory about the F E N E T R e is in preparation. Previous spectral measurements of such gliding discharges already indicated striking differences between rotational, vibrational and electron temperatures. A low DC-current (0.15 to 0.25 A ) GlidArc, burning in slightly hydrogenated argon, was studied by Sember and Svoboda [6]. In their experiment the mean excitation temperature was measured from the Boltzmann plot of A r l lines and the electron density N was determined from the hydrogen 11^ line broadening. In the discharge initiation point, these quantities were equal to 4.3 ± 0.4 k K and 3 x 1 0 c m , respectively. In the developed discharge zone the temperature and N values decreased to 3.6 ± 0.2 k K and 2 X 1 0 c m , respectively. The electron concentrations were therefore much higher than expected for the plasma in the local thermodynamic equilibrium (LTE) state at the measured temperatures. The N values, derived form the Up broadening, corresponded rather to the temperature of 7 k K (in the L T E state). When more hydrogen was added to the discharge (2%),the electron density in the initiation and main part of the GlidArc increased to 1 x 1 0 and 5 x 1 0 c m respectively, while the excitation temperatures changed only a little. e
14
- 3
14
- 3
e
e
15
14
- 3
Another diagnostics of the gliding AC-discharge burning in air, was recently presented by Tesar et al. [7]. For the average arc current of 0.1 A , authors described the following discharge areas between horn-shaped electrodes: — a disruptive spark discharge (SD) zone in the narrowest part of electrodes separation and closely above it (w 30 m m in the gas flow ^direction), with a high electric field strength and current density in the discharge channel, — a disintegrating corona and spark discharges ( D D ) , expanded by the blowing stream of air, having the length of about 70 mm. A l l temperatures in the SD area, determined from the N 2 and O H molecular emission spectra, were about 500 K higher than temperatures in the D D zone. For example, in the SD zone the vibrational temperature T determined from the population of the upper C / 7 state of N 2 , was equal to 3740 K , while the rotational temperature T determined from the same N 2 band was equal to 2600 K . The corresponding values for the D D zone were equal to 3200 and 2060 K , respectively. In the present paper, the non-equilibrium properties of the low-current arc will be proved by spectroscopic measurements of the N 2 second positive system and electrophysical observations. v
3
U
r
2. E x p e r i m e n t a l m e a s u r e m e n t s a n d results The gliding arc device is schematically shown in F i g . 1. Two diverging blade-shaped steel electrodes are placed in a relatively fast air flow (> 10 m/s) as
Spectral and Electrical Diagnostics of Gliding Arc
597
NON-EQUILIBRIUM Z O N E
FENETRe
QUASI-EQUIUBRIUM ZONE
Electrode"
Fig. 1.
Gliding arc electrical scheme.
indicated in the figure. The electrodes are connected to a commercial high-inductance neon-transformer (50 Hz, 10 k V at open circuit, 0.13 A at the short circuit). The spark discharge starts at the place of the shortest distance (1-2 mm) between two electrodes. During the time of about 1 |is a low resistance plasma channel is formed which causes a fall of the voltage between electrodes. This small plasma filament, in the state close to the L T E , is then dragged up by the gas flow and the arc column length increases together with the voltage. The arc current is almost constant during this quasi-equilibrium period but the arc voltage is increasing. The total arc power increases also almost linearly with the discharge length. This evolution is terminated, when the arc length approaches to a critical value, for which the growing plasma resistance becomes equal to the external one, and the maximum value of discharge power corresponds to the half of generator maximal power. A t the same time, the voltage, the electric field and power approach their maximum values. When the moving arc length exceeds the critical value, the arc heat losses continue to increase. But, since the electrical power cannot increase any more, it is no longer possible to sustain the arc in the quasi-thermodynamic equilibrium state. The gas temperature falls rapidly and the plasma conductivity is maintained by the high value of electron temperature and stepwise ionization. After this fast transition ( F E N E T R e , see F i g . 1), the gliding arc continues its evolution but in non-equilibrium state. The specific heat losses in this regime become much smaller than in equilibrium one. Therefore, the discharge length can increase up to other critical value, when the arc breaks and a new discharge starts again from the initial breakdown. The final breaking point is defined by the geometry of electrodes, by the flow conditions and characteristics of the power supply. In the case of an A C high-inductance power supply, one can observe a phenomenon of the multiple zero-current crossings, without the discharge interruption, during several periods of the discharge evolution. About 80% of the total gliding arc power is dissipated after the F E N E T R e , during the non-equilibrium period. It means that gliding arc may really be useful as an efficient and very simple generator of the non-equilibrium plasma [4,5].
598
A. Czernichowski et al.
To show directly the non-equilibrium nature of the gliding arc burning in air at atmospheric pressure, the vibrational ( T ) and rotational (T ) temperatures were measured in the initial (quasi-equilibrium) and final (post-FENETRe) discharge regions. The rotational temperature is supposed to be close to the kinetic gas temperature To. The electron temperature T and free electron concentration N were evaluated from electric parameters and other discharge characteristics for both, the initial and final discharge regions (see below). The "length explosion" phenomenon during transition stage ( F E N E T R e ) has been already demonstrated [5] by high speed camera pictures of another medium-current gliding arc burning in air. v
r
e
e
2.1. Vibrational and rotational temperatures Because of quasi-periodic nature of the gliding arc (with 5 to 30 ms as a quasi-period) the time integrated emission spectroscopy was used to obtain a sufficiently high signal. Light was collected by two lenses, imaged onto the input slit (0.15 mm) of the 1 m McPherson monochromator. The monochromator was equipped with the 1200 G / m m ruled grating, blazed at 300 nm and an U V sensitive photomultiplier as the light detector. The photomultiplier signal was averaged by the P A R boxcar averager. The boxcar output signal was stored in a 4-channel digital oscilloscope. Two observation areas were selected: the quasi-equilibrium "before the F E N E T R e " zone and non-equilibrium "after the F E N E T R e " . A n analysis of time integrated air emission spectra shows intense molecular bands in the region from 200 nm to 400 nm (with maximum around 350 nm). Band shapes and their relative intensities depend strongly on the observation position. The N 2 , N " , N O and O i l bands were detected, the last one was due to an ambient air humidity which was used as plasma-forming gas. Some of bands had quite low intensities in the "after F E N E T R e " region, and were merged in the signal noise or even disappeared (e.g. the No band). The spectra were registered during a slow scanning (1 nm/min) of the monochromator. 4
In order to trigger the boxcar, while the arc filament image was passing the monochromator input slit, the light reflected from the input slit blades was collected by a quartz fiber and detected by a second photomultiplier. Output from the second photomultiplier was used to trigger the boxcar at the moment of maximum light intensity. The boxcar gate width ranged from 10 to 100 /JS depending on the length of the light pulses. 3
3
Bands of the 2nd positive system of nitrogen molecule C V7„-Z? 77 in At; = ± 2 , ± 1 , and 0 sequences were studied in different positions along the axis of the discharge. The monochromator spectral band pass (0.1 to 0.2 nm) was low enough to recover partially resolved structure of the emission band for determination of vibrational and rotational temperatures. Spectra of the 2nd positive system were analyzed by the numerical procedure searching for the best T and T temperature values to fit theoretically calculated spectrum to the experimental band shape. The monochromator instrumental function was measured by means of a low-pressure mercury lamp spectrum. g
v
r
Spectral and Electrical Diagnostics of Gliding Arc
599
Two different fitting procedures were used: first one (marked as A) was elaborated in the Institute of Plasma Physics, Czech Academy of Sciences, and second one (B) which originates from the University of Orleans. No correction for the spectral sensitivity of the optical system has been done. It was not necessary because of the proper combination of the photomultiplier type and grating efficiency, this spectral sensitivity may be considered as constant for analysis of the applied rotational spectrum. Using the A procedure [8,9] the rotational temperatures were found by fitting the calculated spectrum to the beginning of At> = ± 2 , ± 1 or 0 sequences, after the numerical smoothing of data and background subtraction. The calculation model was based on high-resolution analysis of the rotational system by Iioux et al. [10,11] and Honl-London factors [12]. The (0,2), (0,1) and (0,0) N bands were used for T determination in the "before F E N E T R e " region (almost the same T were obtained from all bands), whereas only (0,1) band could be analyzed for the " p o s t - F E N E T R e " region. For vibrational temperatures, the data evaluation technique consisted of the band-profiles generation for given T and instrumental function, with subsequent subtraction from experimental spectra in order to find the best fit. Resulting relative populations of the C 7 7 state in two vibronic levels of v = 0 and i; = 1 determined the T value. Using the B procedure [13-16], the rotational and vibrational temperatures were found by unique fit over a large complex spectrum of the overlapping N 2 bands covering Av = ± 2 , ± 1 and 0 sequences, without: preliminary numerical smoothing, baseline subtraction and assumption that apparatus function is known. The band simulation was based on the same method as in procedure A. Generally, our program permits to simulate a molecular band without limits in wave numbers v' and v", but in the case of the C / 7 - i ? 7 7 g transitions and in order to gain a time, we have simulated only the overlapped (0,0), (1,1), (2,2), (3,3), (4,4), (0,1), (1,2), (1,0), (2,1), (0,2), (1,3), (2,4) and (3,5) bands between the selected signal wave numbers. To determine the T and T from a given measured noisy spectrum, the simulated spectrum was compared point-by-point using a standard best-fit criterion. Four parameters were fitted: T , T , the apparatus width (supposed to be Gaussian) and background level. 2
r
r
r
3
u
v
3
3
u
t
v
r
v
To test the noise sensitivity of the determined temperature, we have simulated noisy "experimental" spectra: for the apparatus profile corresponding to the Gaussian width of AX = 7 c m , the noise-to-signal ratio of 0.025, 0.05 and 0.1 and T = T = 3 k K . The applied best-fit criterion was found to be the good one, since the relative error of the derived temperatures had the same order of magnitude as the noise-to-signal ratio in the applied simulated "experimental" spectra: ATjT = 0.015, 0.025 or 0.04, respectively. - 1
t
v
Figure 2 presents four measured and four corresponding calculated spectra (procedure B) for two characteristic regions of the gliding arc. The vibrational and rotational temperatures, found in initial equilibrium and final non-equilibrium zones of the gliding arc are given in Table. The experimental results obtained from the first fitting procedure are marked as A , those derived from the second procedure are marked as B. From Table one can see that in the initial quasi-equilibrium zone, the vibrational and rotational temperatures are almost the same and quite high,
600
A. Czernichowski et al. 3000 K I V . 3300 K
OX - 6 8 cm'
1
Simulated Spectrum .01
i
i
«
«
5. i
30000
30500 TR -
1100 K TV = 2200 K 0 1 = 6.8 c m "
Simulated Spectrum
27700
27900
28100
2B300
Wavenumber (cm-') Fig.
2.
A n example
of the simulated
a n d observed
spectra
both
for before a n d
p o s t - F E N E T R e z o n e s i n t w o different r a n g e s o f w a v e n u m b e r s .
T A B L E Experimental
characteristics of the gliding
non-equilibrium
Zone
T
arc equilibrium a n d
zones.
e
T
v
T =T r
A
0
e
Quasi-equilibrium 10
13
cm-
3
3900 ± 100 K 1 0
13
cm"
3
Model A
6 kK
3300 ± 200 K 2600 ± 300 K
Model B
6 kK
3900 ± 100 K
Non-equilibrium Model A
10 k K
2600 ± 300 K
1700 ± 4 0 0 K
10
12
cm"
3
Model B
10 k K
2200 ± 200 K
1100 ± 300 K
10
12
cm"
3
Spectral and Electrical Diagnostics of Gliding Arc
601
about 4 k K . In the final non-equilibrium zone, the vibrational temperature remains on the high level of about 3 k K , but the rotational temperature falls to a relatively low level of 1-1.5 k K . 2.2. Electron temperature and concentration Free electron parameters were determined for the same regions of the gliding arc where spectroscopic measurements were done. These parameters were estimated from measurements of instantaneous values of the discharge voltage and current. Electrical signals corresponding to arc voltage and current were registered and stored by a four-channel digital oscilloscope. Signals from the photomultiplier, used for described above plasma spectroscopy, were collected simultaneously by the same oscilloscope, in order to have temperatures, voltage and current measured at the same moment. Values of the near electrodes voltage drop were measured separately. It is interesting that the voltage drop in the initial, "before F E N E T R e " region, is usually lower than that in the "after F E N E T R e " region. It corresponds to the non-equilibrium nature of the second region. Spatial characteristics of the discharge, lengths and cross-sections, were estimated also at the same moments and places of the gliding arc as the spectroscopic and electrophysical measurements. The electric field E was determined from the measured values of discharge voltages and lengths. The reduced electric field E/n — EkTo/p was then calculated from known pressure p and translational temperature To = T (obtained from spectroscopic experiments). The electron temperature was then determined from its known dependence on the E/n ratio in the air [17,18]. Corresponding results are given in Table. The electron concentration N was derived from discharge current density and electron mobility (which is estimated also as known function of the specific electric field E/n in air) [17,18]. Corresponding results are also presented in Table. One can see from Table that in the quasi-equilibrium zone of the gliding arc the electron temperature is not much different from other temperatures, and the electron concentration is relatively high. Our results are therefore in complete accordance with the theoretical model considered in Refs. [4] and [5]. In the non-equilibrium zone, the electron temperature T w 10 k K is much higher than the vibrational one T « 3 k K , which in turn is higher than rotational and translational ones: T ss To ~ 1 to 1.5 k K . Such a hierarchy of temperatures, as well as relatively low electron concentration, with respect to the equilibrium plasma, are in good agreement with the theory of the non-equilibrium plasma [17, 18]. On the other hand, such plasma parameters are very typical of powerful non-equilibrium discharges like a microwave plasma [18]. r
e
e
v
r
3. C o n c l u s i o n s The complicated time-space structure of the gliding arc was experimentally investigated. It was shown by spectroscopic and electrophysical measurements that the initial quasi-equilibrium ("before F E N E T R e " ) plasma of the gliding arc has a temperature of the order 4-6 k K . The share of dissipated electrical power in this
602
A. Czernichowski el al.
zone can be about 20-25% of the total power, which is a relatively small part of the total discharge energy. It was proved experimentally that "post-FENETRe" plasma has strong non-equilibrium parameters. According to the spectroscopic and electrophysical diagnostics, the electron temperature was there about 10 k K , the molecular gas vibrational temperature was about 3 k K and was consequently much higher than the rotational and translational ones (1 to 1.5 k K ) . Similar results were shown independently in Refs. [6] and [7]. 12
- 3
Electron concentration after the F E N E T R e phase is about 1 0 c m , smaller than in the initial zone, but the main part of the gliding arc power (up to 80%) could be dissipated just in this non-equilibrium zone. The transition between two mentioned discharge zones includes the phenomenon called the gliding arc "length explosion". This spectacular phenomenon illustrates the F E N E T R e phase itself and was experimentally observed with a high speed camera [5,19,20]. The gliding arc, the continuous powerful atmospheric pressure non-equilibrium discharge, can be successfully applied for different chemical reaction [3] with the high energy efficiency and selectivity as well as for example non-equilibrium microwave discharge. It can be illustrated already by high efficiency of C O 2 dissociation, I I 2 S decomposition, and sulfur reduction from S O 2 in gliding arc. Moreover, the moderate gas temperature (7b = 1-1.5 k K ) permits to perform in the gliding arc chemical reactions with complicated kinetics (for example such as N 0 oxidation to N O + N O 2 , or dehalogenation of concentrated organic), which could not be effectively stimulated in the room temperature non-equilibrium plasma. 2
References [l] H. Lesueur, A . Czernichowski, J. Chapelle, Coll. Phys., Suppl. au J. Phys. 51, C5-49 (1990). [2] II. Lesueur, A . Czernichowski, J. Chapelle, Coll. Phys., Suppl. au J. Phys. 51, C5-57 (1990). [3] A . Czernichowski, Pure Appl. Chem. CC, 1301 (1994). [4] V . D . Rusanov, A.S. Petrusev, B . V . Potapkin, A . A . Fridman, A . Czernichowski, J. Chapelle, Sov. Phys.-Doklady 332, 306 (1993). [5] A . A . Fridman, A . Petrousov, J. Chapelle, J . M . Cormier, A . Czernichowski, H. Lesueur, J . Stevefelt, J. Phys. Ill (France) 4, 1449 (1994). [6] V . Sember, P. Svoboda, in: Proc. 16th Symp. on Plasma Physics and Technology, Prague 1993, p. 230. [7] C. Tesar, J . Janca, V. Kapicka, in: Proc. Xlth Int. Symp. on Physics of Switching Arc, Brno 1994, p. 74. [8] M . Simek, Acta Phys. Slov. 42, 83 (1992). [9] S. DeBenedictis, G . Dilecce, M . Simek, Acta Phys. Univ. Comeniae 33, 159 (1992). [10] F. Roux, F. Michaud, M . Vervloet, J. Mot. Spectrosc. 07, 143 (1989). [11] F. Roux, F. Michaud, M . Vervloet, J. Mol. Spectrosc. 158, 270 (1993). [12] I. Kovacs, Astrophys. J. 145, 643 (1966).
Spectral and Electrical Diagnostics of Gliding Arc
603
[13] A . Czernichowski, J. Phys. D 20, 559 (1987). [14] H . Nassar, A . Czernichowski, J. Phys. (France) 51, C5-289 (1990). [15] H . Nassar, A . Czernichowski, Acta Phys. Pol. A 84, 215 (1993). [16] H . Nassar, A . Czernichowski, in: Proc 11th Int. Symp. on Plasma Chemistry, ISPC-11, Loughborough (England) 1993, Vol. 1, p. 487. [17] Y.P. Raizer, The Physical Fundamentals of Gas Discharge Processes, Nauka, Moskva 1980 (in Russian). [18] V . D . Rusanov, A . A . Fridman, Physics of Chemically Active Plasma, Nauka, Moskva 1984 (in Russian). [19] J . M . Cormier, F. Richard, J. Chapelle, M . Dudemaine, in: Proc. 2nd Int. Con}, on Electr. Contacts, Arcs, Apparatus and their Applications, Xi'an (China) 1993, p. 40. [20] J . M . Cormier, F. Richard, J . Chapelle, in Proc. 7th Int. Conf. on Switching Arcs Phenomena, Lodz (Poland) 1993, p. 313.