Pulse discharge in the dielectric cell: simulation via PIC method
2D electrostatic PIC code for simulation of the pulse discharge in the dielectric cell is described. The first simulation results (discharge current temporal dependence, electric potential spatial distribution, electrons' energy distribution) for the discharge in Ne - Xe mixture are presented.
Saved in:
| Date: | 2007 |
|---|---|
| Main Authors: | , , |
| Format: | Article |
| Language: | English |
| Published: |
Національний науковий центр «Харківський фізико-технічний інститут» НАН України
2007
|
| Series: | Вопросы атомной науки и техники |
| Subjects: | |
| Online Access: | http://dspace.nbuv.gov.ua/handle/123456789/110510 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Pulse discharge in the dielectric cell: simulation via PIC method / O.I. Kelnyk, O.V. Samchuk, I.O. Anisimov // Вопросы атомной науки и техники. — 2007. — № 1. — С. 148-150. — Бібліогр.: 5 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
irk-123456789-110510 |
|---|---|
| record_format |
dspace |
| spelling |
irk-123456789-1105102017-01-05T03:04:20Z Pulse discharge in the dielectric cell: simulation via PIC method Kelnyk, O.I. Samchuk, O.V. Anisimov, I.O. Low temperature plasma and plasma technologies 2D electrostatic PIC code for simulation of the pulse discharge in the dielectric cell is described. The first simulation results (discharge current temporal dependence, electric potential spatial distribution, electrons' energy distribution) for the discharge in Ne - Xe mixture are presented. Описано двовимірний код для електростатичного моделювання імпульсного розряду методом частинок у комірках. Наводяться перші результати моделювання розряду у суміші неону та ксенону(часова залежність розрядного струму, просторовий розподіл потенціалу, розподіл електронів по енергіях). Описывается двумерный код для электростатического моделирования импульсного разряда методом частиц в ячейках. Приводятся первые результаты моделирования разряда в смеси неона и ксенона (временная зависимость разрядного тока, пространственное распределение потенциала, распределение электронов по энергиям). 2007 Article Pulse discharge in the dielectric cell: simulation via PIC method / O.I. Kelnyk, O.V. Samchuk, I.O. Anisimov // Вопросы атомной науки и техники. — 2007. — № 1. — С. 148-150. — Бібліогр.: 5 назв. — англ. 1562-6016 PACS: 52.80.Tn, 52.90.+z http://dspace.nbuv.gov.ua/handle/123456789/110510 en Вопросы атомной науки и техники Національний науковий центр «Харківський фізико-технічний інститут» НАН України |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| language |
English |
| topic |
Low temperature plasma and plasma technologies Low temperature plasma and plasma technologies |
| spellingShingle |
Low temperature plasma and plasma technologies Low temperature plasma and plasma technologies Kelnyk, O.I. Samchuk, O.V. Anisimov, I.O. Pulse discharge in the dielectric cell: simulation via PIC method Вопросы атомной науки и техники |
| description |
2D electrostatic PIC code for simulation of the pulse discharge in the dielectric cell is described. The first simulation results (discharge current temporal dependence, electric potential spatial distribution, electrons' energy distribution) for the discharge in Ne - Xe mixture are presented. |
| format |
Article |
| author |
Kelnyk, O.I. Samchuk, O.V. Anisimov, I.O. |
| author_facet |
Kelnyk, O.I. Samchuk, O.V. Anisimov, I.O. |
| author_sort |
Kelnyk, O.I. |
| title |
Pulse discharge in the dielectric cell: simulation via PIC method |
| title_short |
Pulse discharge in the dielectric cell: simulation via PIC method |
| title_full |
Pulse discharge in the dielectric cell: simulation via PIC method |
| title_fullStr |
Pulse discharge in the dielectric cell: simulation via PIC method |
| title_full_unstemmed |
Pulse discharge in the dielectric cell: simulation via PIC method |
| title_sort |
pulse discharge in the dielectric cell: simulation via pic method |
| publisher |
Національний науковий центр «Харківський фізико-технічний інститут» НАН України |
| publishDate |
2007 |
| topic_facet |
Low temperature plasma and plasma technologies |
| url |
http://dspace.nbuv.gov.ua/handle/123456789/110510 |
| citation_txt |
Pulse discharge in the dielectric cell: simulation via PIC method / O.I. Kelnyk, O.V. Samchuk, I.O. Anisimov // Вопросы атомной науки и техники. — 2007. — № 1. — С. 148-150. — Бібліогр.: 5 назв. — англ. |
| series |
Вопросы атомной науки и техники |
| work_keys_str_mv |
AT kelnykoi pulsedischargeinthedielectriccellsimulationviapicmethod AT samchukov pulsedischargeinthedielectriccellsimulationviapicmethod AT anisimovio pulsedischargeinthedielectriccellsimulationviapicmethod |
| first_indexed |
2025-07-08T00:41:24Z |
| last_indexed |
2025-07-08T00:41:24Z |
| _version_ |
1837037293395771392 |
| fulltext |
148 Problems of Atomic Science and Technology. 2007, 1. Series: Plasma Physics (13), p. 148-150
PULSE DISCHARGE IN THE DIELECTRIC CELL:
SIMULATION VIA PIC METHOD
O.I. Kelnyk, O.V. Samchuk, I.O. Anisimov
1Taras Shevchenko National University of Kyiv, Radio Physics Faculty,
Volodymyrs’ka Str. 64, 01033, Kyiv, Ukraine, e-mail: oles@univ.kiev.ua
2D electrostatic PIC code for simulation of the pulse discharge in the dielectric cell is described. The first simulation
results (discharge current temporal dependence, electric potential spatial distribution, electrons' energy distribution) for
the discharge in Ne - Xe mixture are presented.
PACS: 52.80.Tn, 52.90.+z
1. INTRODUCTION
Gas discharges in the dielectric cells are common for
many laboratory and industrial applications. Many of
these applications, especially related to the discharges of
microscopic sizes, present difficulties for experimental
investigations (see, e.g., [1]). For that case, computer
simulation can be very useful for the investigation of the
processes in gas discharges. For now, most common ap-
proach to the computer simulation of the gas discharges is
based on the numerical solution of the kinetic equations
for the elementary processes [2-4]. This method is natural
for the stationary homogeneous plasma systems, but it
does not fit well for the non-stationary gas discharges
inside the small dielectric cells. Another methods, such as
solving of the hydrodynamic equations (see, e.g. [5]), take
into account the plasma inhomogeneities, but do not con-
sider the kinetic effects.
In this work the pulse gas discharge in the dielectric
cell is studied via computer simulations using Large Par-
ticles in Cells method.
2. COMPUTATION ALGORITHM
2D electrostatic PIC model has been applied in the
simulation code. That code was developed for the PC
platform (MS Windows) with user friendly interface (see
Fig.1).
At each simulation step, equations for the electric po-
tential and field were solved on the mesh with variable
step using the matrix sweep method. This method is based
on the solving of the finite differences’ equation set in the
shape of matrix three-diagonal equations. The method has
a good accuracy for non-uniform spatial meshes so the
simulation gives the reasonable results even for large
(about 106) amount of time steps. The time performance
of matrix sweep method is sufficiently decreased for PIC
simulations because main volume of calculations must be
performed only once at the first simulation step.
Based on the values of electric field, the new values of
the particles’ coordinates and velocities are found out
from the motion equations.
3. PROCESSING OF THE PARTICLES’
COLLISIONS
Key part of the simulation code is the processing of
the particles’ collisions. According to the code purpose, it
is devoted to the simulation of weakly ionized plasma that
is typical for gas discharge devices like PDP. So all sorts
of elementary processes taken into account can be divided
Fig.1. Program simulator window for the elementary processes taken into account
mailto:oles@univ.kiev.ua
149
in two different classes – collisions with neutrals and col-
lisions with other particles. Neutrals in ground (non-
excited) state are not treated as sorts of large particles but
form a background. Based on the free path of the particles
and probabilities of all possible elementary processes with
ground state neutrals, these processes are simulating using
the Monte Carlo method. The elementary processes taken
into account are non-elastic collisions, excitations on dif-
ferent levels, ionization, recombination, photon emission
and absorption (the radiation transport is also considered).
The collisions between the particles and cell walls are
also considered as well as the secondary emission from
these walls.
4. SIMULATION PARAMETERS
Simulation was carried out for the cell of plasma dis-
play panel. Simulation parameters are given below. Dis-
charge cell with dielectric walls has dimensions
500×200 µm, partial pressures of neutral gases are
450 Torr for neon and 50 Torr for xenon. Direct driven
voltage of 200 V is applied upon 200 µm cell side and is
turned on at t=0. The initial portion of large particles’
contains 50 electrons, 25 ions Ne+ and 25 ions Xe+. Simu-
lation time step was 10-12 s. Large particle sorts taken into
account included electrons, Ne and Xe ions and Ne and
Xe excited particles with excitation energy enough to
excite the phosphor – simulation was carried out for the
cell of plasma display panel. We made 1 million steps of
simulation (that correspond to real driven voltage pulse
length) and controlled the amount of radiated photons.
5. SIMULATION RESULTS
Using the code mentioned above, computer simulation
of the pulse discharge in the dielectric cell was carried
out. The results for the plasma density and electric field
spatial distribution and for the electron energy distribution
are in good accord with the respective experimental re-
sults.
Current temporal dependency on Fig.2 corresponds to the
known facts about the gas discharge in the dielectric cell.
Such a discharge is initiated by the electron avalanche that
appears in the electric field of driven voltage (front of the
current pulse on Fig.2). Moving in this field, charged parti-
cles can reach the isolated electrodes and adsorb on the di-
electric surfaces, so those planes are charging and counter
voltage is appearing. When that voltage compensates the
driven voltage, the discharge initiating field disappears and
discharge starts to extinguish (current pulse back front). The
duration of discharge current pulse is determined by the time
of the charging of cell electric capacitance (about 2⋅10-13 F
for that case) by this current. Maximum charge for the 200 V
voltage must be about 4⋅10-11 Cl. So the duration of almost
triangular current pulse with 300 µA magnitude must be
about 250 ns that corresponds to duration of the pulse pre-
sented on Fig.2. The shape of the discharge pulse is also in
good accordance with the experimental dependence [1].
Evolution of the electric potential spatial distribution
during the duration of discharge current pulse is shown on
Fig.3. One can see that discharge positive column appears
near the positive electrode at the beginning of this pulse
(Fig.3a).
Fig.2. Temporal dependencies of driven voltage and
discharge current during the simulation
a
b
c
Fig.3. Electric potential spatial distribution during the
simulation: a) t=0.51ns, b) t=0.52ns, c) t=0.54ns
150
Then positive column expands quickly (Fig. 3b), so the
duration of the current pulse front is rather small – about
10 times smaller then the entire pulse duration. At the
moment directly after the discharge current reaches its
maximum value (Fig. 3c), one can see that positive col-
umn is quite close to the negative electrode, so the dielec-
tric surface of electrodes’ isolation is charged quickly,
that tends to discharge extinguishing.
Fig.4. shows the evolution of the electron energy dis-
tribution during the discharge current pulse duration. The
small initial portion (about 100 large particles) assigns at
the beginning of simulation (t = 0) with random velocities
distributed uniformly. As more of new particles appear
inside the cell due to the elementary processes, such as
ionization and excitation, electron energy distribution
tends to be closer to Maxwellian shape.
a
b
c
Fig.4. Electron energy distribution during the simulation:
a) t=0.51 ns, b) t=0.52 ns, c) t=0.54 ns
On Fig.4 one can see the electron energy distributions for
the same moments as the potential dependencies on Fig.3.
At the beginning of current pulse (Fig. 4a), about 58000
electron large particles are distributed in a Maxwellian-
like shape, but with diffused maximum. At the middle of
pulse front (Fig. 4b – about 190000 electron large parti-
cles) and, especially, for the moment near the discharge
current maximum (Fig. 4c – about 280000 electron large
particles) the electron energy distribution practically cor-
responds to Maxwellian law.
6. CONCLUSIONS
1. Two-dimensional code for simulation of weakly
ionized plasma systems (such as discharge devices,
plasma display panels etc.) is developed and tested for the
case of pulse discharge in dielectric cell.
2. Spatial distribution of the potential inside dielectric
cell during the discharge current pulse has a positive col-
umn region that is quickly expanding and, finally, fills
almost the entire cell. Then the discharge extinguishes.
3. Electron energy distribution during the current
pulse changes towards the Maxwellian shape.
ACKNOWLEDGEMENTS
Authors thank S.M.Levitsky (Taras Shevchenko Na-
tional University of Kyiv) and Sung Chun Choi (LG Elec-
tronics Inc.) for fruitful discussion.
This work was financially supported by LG Electron-
ics Inc.
REFERENCES
1. J.P. Boeuf. Plasma display panels: physics, recent
development and key issues // J. Phys. D: Appl. Phys.
2003, v. 36, p.R53-R79.
2. M. Surendra. Radiofrequency discharge benchmark
model comparison // Plasma Sources Sci. Technol. 1995,
v. 4, N 1, p.56-73.
3. M.G. Zubrilin, G.G. Kalyuzhna, I.A. Popov,
A.I. Tschedrin. Comparative Characteristics of Excimer
XeCl Laser Based on He/Xe/HCl and He/Xe/CF2Cl2 //
Ukr. Phys. Journ. 2005, v. 50, N5, p.442-447 (In Ukrainian).
4. K. Hassouni, G. Lombardi, X. Duten, G. Hangelaar,
F. Silva, A. Gicquel, T.A. Grotjohn, M. Capitelli,
J. Ropcke. Overview of the different aspects in modelling
moderate pressure H2 and H2/CH4 microwave discharges
// Plasma Sources Sci. Technol. 2006, v.15, N1 p.117.
5. V.V. Osipov, V.V. Lisenkov. Formation of the self-
maintained volumetric gas discharge// ZhTF. 2000, v.70,
N10, p.27-33 (In Russian).
:
. , . , .
. (
, ,
).
:
. , . , .
. (
, , ).
|