Nanosecond Pulsed Dielectric Barrier Discharge in Air

Nanosecond Pulsed Dielectric Barrier Discharge in Air

Dielectric Barrier Discharges (DBD) are stable discharges capable of generating reactive non-equilibrium plasmas at high pressures.  Operation of these discharges at atmospheric pressure is particularly attractive for a number of applications including vacuum-chamber free materials processing, chemical processing of gas streams, and plasma actuators for aerodynamic flow control.  In non-ideal gases such as air, DBD’s produce large volume plasmas through the formation of a multitude of thin streamer channels that can fill the available space in the discharge.  DBD’s are essentially pulsed discharges that are produced by applying a high voltage oscillating waveform or a series of repeated sharp-rising pulses at an electrode in an electrode pair system.  One or both of the electrodes must be covered by a dielectric solid layer so as to prevent a direct

conducting channel between the electrodes.  The dielectric layer traps charge from the plasma which in turn limits a large current density in the plasma.  The self-limiting behavior of the dielectric barrier keeps the currents low enough that glow-to-arc (GAT) instabilities are prevented, thereby maintaining a stable non-equilibrium discharge at high pressures.  The self-limiting behavior also means that the discharge is necessarily pulsed with the formation of streamer channels that are quenched rapidly during a single pulse.  An oscillating excitation is therefore necessary to maintain a semi-continuous discharge.

This post discusses the simulation of a prototypical DBD in air.  The VizGlow Plasma Modeling Software Package is used.

The geometry for the DBD is shown in Fig. 1 and comprises a bare powered electrode(red), a surrounding ground electrode (green) that is covered by a dielectric layer (black).  A sharp-rising triangular positive pulse is applied at the electrode as shown in Fig. 2.  The pulse duration last about 8 ns and therefore the discharge is called a “nanosecond pulsed DBD”.  The DBD is generated at atmospheric pressure in air.  Since the high pressure imposes a very fine resolution requirement, over 80,000 cells are required in the mesh.  The mesh comprises pure quadrilateral cells.  Simulation domain is 2.5 mm x 0.5 mm (see Fig. 1).  The air is represented by a finite-rate chemistry with 21 gas-phase reactions among 11 species (electrons E, oxygen radical O, nitrogen molecule N2, oxygen molecule O2, nitrogen dimer ion N2+, oxygen dimer ion O2+, nitrogen cluster ion N4+, oxygen cluster ion O4+, ion complex O2+N2, oxygen dimer negative ion O2-, oxygen ion O-).  Since air species comprise radiative states that emit in the ultra-violet, photoionization becomes an important mechanism by which electrons are generated and therefore this mechanism is included in the simulation.

The pulsed excitation results in the formation of a positive streamer channel that propagates from the powered electrode to the grounded electrode. However, the presence of the dielectric barrier on the grounded electrode results in the deposition of a positive charge on the dielectric surface.  This locally trapped positive charge causes a self-induced electric field that pushes the streamer further down the dielectric surface resulting in a streamer channel that propagates significant distances from away from the powered electrode.  Essentially, the dielectric charge trapping is the cause for the large volume filling property of the DBD.  Figure 3 shows a series of snap-shots of the electron density in the streamer channel formed by the nanosecond pulsing of the DBD in air.  Note that the streamer channel is thin (~ 100 microns).  However, a much broader region of electron densities is observed surrounding the streamer.  This is essentially a cloud of electrons resulting from photoionization surrounding the streamer.  Note that the streamer propagates about 1 mm from the powered electrode in about 7 ns and then extinguishes once the power is turned off.

Figure 4 shows the evolution of electrostatic potential in the discharge.  The potential starts out looking like a vacuum potential profile during the early part of the transient (~ 1 ns).  However as the streamer develops the self-consistent nature of the electric field induced by the streamer channel becomes evident.  The streamer head (leading edge of the propagating streamer) is electronegative due to an excess of electrons that distorts the potential imposed by the powered electrode.  At the end of the transient the potential collapses with the quenching of the plasma.

Finally, a comparison is made of the streamer propagation with the photoionization and without the photoionization.  Figure 5 shows propagation of the streamer for the two cases at three different times during the transient.  Note that the case without the photoionization results in a much weaker streamer that propagates a very short distance away from the powered electrode before being quenched at the end of the pulse.  Also note that the absence of photoionization results in a strongly confined electron density profile without the cloud of electrons that surrounds that streamer in the case with photoionization.  This direct comparison emphasizes the role of photoionization in streamer propagation in air DBD’s.

VizGlow provides robust capability for the high-fidelity simulation of high-pressure non-equilibrium discharges such as DBD’s.  Simulation of DBD’s is a particularly challenging problem because of the highly disparate time scales involved in the problem and the potentially wide range of length scales imposed by the small region of discharge formation in the context of a larger application domain.  Furthermore, the high pressure results in very restrictive resolution requirements since the thin propagating streamer must be resolved accurately in order to model the DBD discharge.  This kind of modeling is very much the state-of-the-art in plasma computational modeling community.

VizGlow is part of the Overviz Simulation Suite and is provably fast, robust, easy-to-use, and is currently a leading industrial plasma simulation tool.

DC Non-Transferred Plasma Torch

DC Non-Transferred Plasma Torch

VizSpark is used to model a non-transferred direct current plasma torch in pure argon gas assuming axi-symmetry along the central axis.  A steady state thermal plasma with temperatures from 15-20,000 Kelvin forms a supersonic, under expanded jet in the ambient region with visible Mach diamonds.

Plasma torches are devices that produce a steady thermal plasma jet.  There are many applications for plasma torches including arc cutting, wielding, plasma spraying and waste treatment.

A plasma torch works by passing a high current between a working gas, typically an inert gas such as argon, helium, hydrogen or a mixture of inert gases.  The high current arc heats the gas mixture to high temperatures (15-20,000 K) which then exits the device forming a jet of high temperature plasma.  Plasma torches typically use direct-current (DC) power source applied through a central (hot) powered electrode.  A plasma torch can either user a transferred arc, where the working material located outside the device acts as the anode, or a non-transferred arc where the anode is located inside the torch housing.

The simulation demonstrated here is for a direct-current, non-transferred arc using pure argon as the working gas.

Fig. 1 shows a typically plasma torch device as well as the domain and mesh for the simulation presented here.  For this simulation, an axi-symmetric assumption was made to simplify the mesh and reduce the computational cost.  The domain consists of a 2.5 mm radius tube that expands to a 5 mm.  A 2.0 mm radius pin powered electrode is placed along the axis in the wider tube region.  The walls of the tube act as the anode.  Argon gas is fed into the tube from the bottom boundary at a total pressure of 3 bar, and the ambient outer region is kept at a pressure of 1 bar.

Fig. 2 shows the plasma temperature (top) and gas pressure (bottom) after a quasi-steady state has been reached.  Gas temperatures exceed 20, 000 Kelvin inside the device while gas temperatures in the ambient jet are on the order of 15,000 Kelvin.  The gas leaves the torch housing at supersonic flow velocities (4000 m/s) with pressures higher than the ambient.   This results in an under expanded jet with Mach diamonds clearly visible in Fig.2 (middle) and Fig. 2 (bottom).

VizSpark software fully support parallel capabilities using a domain decomposition approach.  Fig. 3 compares the mesh and resulting temperature for the same simulation conditions using 1 core (serial) and 48 cores (parallel).  Note that the results between the serial and parallel simulation are identical between demonstrating that there is no loss in correctness when using parallel vs serial.

Fig. 4 is a plot of the simulation speedup (strong scaling) for the 80,000 cell plasma torch problem for processor counts varying between 1 and 48.  8 times speedup can be obtained for 24 cores and 16 times speedup is achieved when 48 cores are used.

VizSpark Plasma Modeling Software Package is part of the Overviz framework suite which provides an intuitive interface to set-up a project to be solved using VizSpark, manipulate multiple projects for parametric studies.  VizSpark is provably fast, robust, and easy-to-use software and currently a leading industrial plasma simulation tool.