Analysis of Discharge Parameters and Spectroscopic Diagnostic of DBDs Pooja Gulati Plasma Device Technology, Microwave Tubes Division CSIR-Central Electronics Engineering Research Institute (CSIR-CEERI) Pilani Rajasthan- 333031 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data, Trieste - Italy
CSIR-Central Electronics Engineering Research Institute (CSIR-CEERI), Pilani Foundation was laid on 1953. Around 450 Employees Leading Research Institute in India in the field of Electronics Devices
Major Research Areas at CSIR-CEERI Pilani 1.Microwave Tubes Gyrotron Klystron Magnetron TWT Plasma Devices Technology 2.Electronic Systems Agri-Electronics Embedded System Digital System Power Electronics 3.Semiconductor Hybrid Microcircuits IC Design MEMS and Microsensors anotechnology & Devices Photonics & optoelectronics Semiconductor Material & Tech.
Activities in Plasma Devices Group High Power Plasma Switches: Thyratrons & Pseudospark VUV/UV Excimer Sources based on DBD: Biomedical Applications Surface Treatment Water Purification (jointly with EERI) Plasma Cathode Electron Gun: Electron and Ion Source Plasma Assisted Microwave Sources: Plasma TWT, Pasotron Penning Discharge Devices Ion Sources and VUV Spectroscopy
Organization of Presentation Motivation and objectives of the work. Introduction What is Dielectric Barrier Discharges (DBDs)? Advantages and open areas of research and applications of DBDs. Experiments Experimental Setup and testing Results and discussions. Conclusion
Motivation In recent time it has been observed that Dielectric Barrier Discharge (DBD) based micro-discharges and micro-arraydischarge plasmas can produce ultraviolet radiation in germicidal wavelength range UV-C(200-280nm) and VUV(100-200nm), UV-B(280-315nm), UV-A (315-400nm) that can effectively treat impure water, and also can be used for medical and other industrial applications. Our group is working in this area and I am motivated in the spectroscopic studies related to this technology which can in future transit me for large scale plasma related spectroscopic analysis.
Objectives To investigate and analyze the characteristics of discharge patterns occurring in the volume discharge (VD) configuration of DBDs filled with inert gases. The traditional metallic diagnostic technique is not useful in the very small geometry of the proposed DBD configuration. Hence to derive the internal electrical and plasma parameters with the help of electrical analysis and spectroscopic diagnostics is the key component of the objectives.
Introduction What is Dielectric Barrier Discharge(DBDs)? Dielectric barrier discharges (DBDs), also known as silent discharges or barrier discharges, are generated in discharge configurations with at least one dielectric barrier between the electrodes.
Possible DBD Geometries (i). Volume discharge, (ii). Surface discharge (iii). Coplanar discharge Typical dielectric barrier discharge configurations
Fabricated Geometries
Geometrical Design and Parameters C d1 C g C d2 Dielectric Thickness= 1 mm Electrode diameter =36 mm Electrode Thickness = 1mm Gas gap = 2 mm Dielectric Material Used= Quartz Pressure of Gas=100mbar Gas Used = Helium Dielectric barrier capacitance C d1 = C d2 = 20.48 pf
Experimental Setup Base pressure=1x10-4 mbar Gas filling assembly is used to fill gas at different pressures. Working Pressure~100 mbar Sinusoidal voltage supply up to 2kV peak with frequencies from 30 to 90 khz has been used. Applied voltage and the total current are measured using high voltage probe and Rogowski-type current transformer. Schematic View of experimental setup Oscilloscope and visible spectrometer are interfaced with computer.
Testing & Characterization of He DBD 600 400 200 Va It 0.03 0.02 0.01 Experimental Setup for the DBD Source At breakdown voltage, the discharge begins with some filaments distributed on the dielectric wall. Increasing the applied voltage little bit, number of filaments increases and for further increase in voltage, the discharge finally get diffused. Va(V) 0-200 -400-600 10 15 20 25 30 35 40 45 Time(u sec) 0.00 It(mA) -0.01-0.02-0.03
Using Kirchoff s theorem for the model, we obtain the following equations V = V V (1) I = I I I = I I a d + Total current through DBD and displacement current through gap dvd dvg I dbd = Cd (4) I dg = C g (5) dt dt Diff. (1) with respect to time and replacing (4) and (4) in (1), dva 1 = dt C g ( I dbd I dis I ) + C Rearranging (6), we will get I dis Equivalent electrical circuit of DBDs C g dva ( t ) = (1 + ) I dbd C g C dt d 1 V = + V g dbd d d I dbd dt Vm0 Cd 1 = g Va I dbd dt Vm0 Cd m 1 T / 2 dbd 0 = I dt 2C d 0 tc dbd + The values of dielectric and gas gap voltages are, Where V m0 is memory voltages, ( In case of sinusoidal excitation) V (6) (7) (8) (9) sc (2) dbd dis + dg (3) In case of sinusoidal excitation (10) U Pal et al, J. Phys. D: Appl. Phys. vol. 42, 045213 (8pp), 2009. U Pal et al, J. Phys.: Confe. Ser.208, 012142, 2010.
Results and discussions Volatage (V) 400 200 0-200 Va Vd Vg Vm Input Parameters: Gap capacitance C g (11.30 pf) Dielectric barrier capacitance C d (20.48 pf) Curren nt (ma) -400 14 7 0-7 It Idbd Idis. Equations are used from references: U Pal et al, J. Phys. D: Appl. Phys. vol. 42, 045213 (8pp), 2009. U Pal et al, J. Phys.: Confe. Ser.208, 012142, 2010. Power (W) -14 4 2 0 Psup. Pdis. Experimental waveforms of dynamic processes occurring in gap (gas: Helium at f = 34.5 khz) for the Parallel plate DBD Geometry. at 100 mbar. -2 15 20 25 30 35 40 Time (µs)
Spectroscopic Results 4.5x10 3 4.0x10 3 HeI 7065.1Å 3.5x10 3 Intensity(a.u) 3.0x10 3 2.5x10 3 2.0x10 3 1.5x10 3 HeI 4921.9 Å HeI 5015.6Å HeI 3888.6Å HeII 6559.7Å HeI 5875.6 Å HeI 6678.1Å HeI 7281.3Å 1.0x10 3 5.0x10 2 He I 3888.6 Å (2 3 S-3 3 P) He I 4921.9 Å (2 1 P-4 1 D) 3000 4000 5000 6000 7000 8000 W avelength (Å) eutral Helium Lines: He I 5015.6 Å (2 1 S-3 1 P) He I 5875.6 Å (2 3 P-3 3 D) He I 6678.1 Å (2 1 P-3 1 D) He I 7065.1 Å (2 3 P-3 3 S) He I 7281.3 Å (2 1 P-3 1 S)
Collisional-Radiative (CR) Model Intensities of the He I lines are calculated using collisionalradiative (CR) model based ADAS code [H. P. Summers, ADAS users manual, JET IR 06 (Abingdon: JET Joint undertaking) (1994)]. With an assumption that the average electron density and temperature in an emission length x, the photon intensity I (λ ul ) of a spectral line can be written from the CR-model as, ~ ~ I ~ ( λ ) = PEC ( x ) + PEC ul recombining e i + excitation From CR-model the ground state populations of atoms and ions is given by, d g d i = = α CR Under steady-state approximation, dt dt i e S CR g e ~ e ~ ( g x ) S α CR CR = i g
Under equilibrium condition α S = 0 and So, the condition Purely Ionizing Condition CR i e CR g e i g S = α CR CR S CR For ionizing plasma is >> 1 In true sense the ionizing plasma condition holds well when i / g << S CR / α CR [Fujimoto T. and Sawada K., IFS-DATA-39 (1997)]. α CR S CR and for recombining plasma is << 1 HeI α CR
Intensity line ratios Under ionizing condition the term PEC be negligibly small and the line intensity from level u to level l is expressed as, I ~ ( λ ul ) = PEC excitation ~ ~ recombining I( λ ul e ~ ( g x) ) e ~ ( i x) is taken to for a transition The significance of the line ratio technique is that the experimentally observable intensity ratio of two lines (which is not directly dependent on, and ) can be easily i g e obtained from the code as the ratio of corresponding photon emission rate coefficients is given by, I I 1 2 = PEC PEC 1 2 ( ( e e, T e, T e ) )
Temperature & Density Calculation Temp. Sensitive Intensity Ratios: 7281.3/ 7065.7, 5049/ 4713.1 Density Sensitive Intensity Ratios: 6678.1/ 7281.3, 4921.9/ 5047.7 Calculated Values: Electron Temp. = (6.5±0.5) ev, Density = (3.5±1.5) х10 11 cm -3 Ref. Summers H P 1994 ADAS Users Manual JET IR 06 (Abingdon: JET Joint Undertaking). R. Prakash, et al J. Appl. Phys. vol. 97, no.4, p.043301, 2005.
Simulated results of the electron Density using OOPIC-Pro 6.00x10 11 120nsec 105ns 95ns e (cm -3 ) 4.50x10 11 3.00x10 11 1.50x10 11 0.00 0 5 10 15 20 25 30 35 40 Diameter of the electrode (mm) The statistical mechanics have shown that the many small perturbation (errors) that affect a physical system almost always force the measurement to follow the Gaussian distribution. It is usually referred to as simply the normal distribution. Based on this, the average distribution of electron plasma density is derived for entire system geometry using OOPIC-Pro simulation code and if we take line average of the saturated density it would give nearly similar results to the spectroscopic diagnostic measurements, which are in agreement to each other.
Conclusion The homogeneous type of discharge has been observed at 100mbar operating pressure for a fixed frequency 34.5 khz in parallel plate DBD cell filled with helium gas. The dynamic evolution of the process in the gap provides the useful information about the electrical characterization of the DBD source. The electron plasma temperatures and electron plasma density obtained for present configuration at 100mbar gas pressure are typically(6.5±0.5)evand(3.5±1.5)х10 11 cm -3 respectively. The existence of such density and temperature in this source is useful for existence of higher metastable states which needs to be further investigated.
Acknowledgement Dr. Ram Prakash, CSIR-CEERI, Pilani Mr. U. Pal, CSIR-CEERI, Pilani and all other group members
Thank You