Plasma Assisted Combustion:

Yiguang Ju

Princeton University

2017. 6. 13

Princeton Combustion Summer School

Copyright ©2017 by Yiguang Ju. This material is not to be sold, reproduced or distributed without prior written
permission of the owner, Yijuang Ju

1
 Sang Hee Won
Associate professor,
Acknowledgement Prof. Richard B Miles
Princeton
Univ. South Carolina

Timothy Ombrello, Prof. Walter R. Lempert

Senior research
engineer Dr. Andrey Starikovskiy
Prof. Haibao Mu
AFRL Princeton
XiAn Jiaotong Univ.

Prof. Min Suk Cha

Wenting Sun
KAUST
Assistant professor,
Georgia Tech Prof. Christophe Laux
Ecole Centrale Paris
Prof. Anne Bourdon
Ecole Centrale Paris
Joseph Lefkowitz Prof.Qi Chen
Research fellow Beijing Jiaotong Univ.
AFRL
Prof. Igor V. Adamovich
Ohio State University

Aric Rousso Prof. Rechard Yetter

Graduate student Penn-state Univ.

Xingqian Mao Prof. Svetlana Starikovskaya

Prof. Haixin Wang Ecole Polytechnique
Visiting student
Timothy Chen Beihang Univ.
Graduate student
 Lecture contents and review articles
1. Introduction and plasma discharge
2. Plasma Assisted Combustion and Applications in Engines
3. Effects of plasma on ignition, flame propagation, and minimum ignition energy
4. Effects of electric field on combustion: Joule heating and ionic wind
5. Physics and kinetic studies of plasma assisted combustion
6. Modeling of plasma assisted combustion
7. Perspectives of future research in plasma assisted combustion

Review papers of plasma assisted combustion

1. Ju, Y. and Sun, W., 2015. Plasma assisted combustion: dynamics and chemistry. Progress in Energy and Combustion Science, 48,
pp.21-83.
2. Starik AM, Loukhovitski BI, Sharipov AS, Titova NS. 2015 Physics and chemistry of the influence of excited molecules on combustion
enhancement. Phil. Trans. R. Soc. A 373: 20140341.
3. Igor V Adamovich and Walter R Lempert, 2015, Challenges in understanding and predictive model development of plasma-assisted
combustion, Plasma Physics and Controlled Fusion, Volume 57, Number 1.
4. Starikovskiy A, Aleksandrov N. Plasma assisted ignition and combustion. Prog. Energy Combust. Sci. 2013;39:61–110.
5. Starikovskaia SM. Plasma assisted ignition and combustion. 2006; J. Phys. D: Appl. Phys. 39:R265–R299.
6. CO Laux, TG Spence, CH Kruger, RN Zare, Optical diagnostics of atmospheric pressure air plasmas, Plasma Sources Science and
Technology 12 (2), 125
7. A Fridman, S Nester, LA Kennedy, A Saveliev, O Mutaf-Yardimci, Gliding arc gas discharge, Prog. Energy Combust. Sci. 25 (2), 211-231
3
 1. Introduction: Plasma assisted combustion
1814 – W.T. Brande. Phil.Trans.Roy.Soc, 104, 51. (Electric field-flame interaction)

1860 Étienne Lenoir used an electric spark plug in his gas

engine, the first internal combustion piston engine.
(Spark ignitor for engines)

1948, Calcote, 3rd Symposium on Combustion and Flame, and Explosion

Phenomena (Vol. 3, No. 1, pp. 245-253) (Ionic wind)

(a) XH2,PJ = 0.3, PIN = 3.0 kW, PIN total = 4.3 kW

1981: Kimura, L, et al., Combustion and Flame, Vol. 42, No. 3,

pp. 297- 305 (Plasma jet in supersonic combustion)

(b) XH2,PJ = 0.5, PIN = 6.0 kW, PIN total = 8.2 kW

1998: Starikovskaia, S.M., Starikovskii, A.Y. and Zatsepin, D.V., Journal of
Physics D: Applied Physics, 31(9), p.1118. ]Anikin N B and Marchenko N 2005
(Nanosecond discharge).

(a) Hot diffusion flame

2015: Won, S.H., Jiang, B.,
2013: Leonov, S.B., Firsov, A.A., Shurupov, M.A., Diévart, P., Sohn, C.H. and
Michael, J.B., Shneider, M.N., Miles, R.B. and Popov, Ju, Y., Proc Combust Inst,
N.A., 2012. Physics of Plasmas, 19(12), p.123502. 35(1), pp.881-888. (Plasma
(laser guiding plasma discharge) (b) Cool diffusion flame assisted cool flames)
 Plasma Assisted Combustion: a multi-disciplinary and multi-physics problem
Plasma discharge

Ionic wind
O2+, N2+ Ions/electrons Instability

Plasma Physics Temperature NO, O3 Radicals Fuel

O, H, OH Int. species fragments
Electric field increase
Joule heating
Electron collision reactions
Excited H2 , CO
N2*, N2(v) species CH4
O2 (a1Δg) CH2O
 Charged species
 Excited species  Ionic wind

Plasma Thermal Kinetic Transport

combustion
Chemical Kinetics studies
Flame Dynamics
Reaction pathways
Reaction rates Traditional
 Extinction Combustion Enhancement
 Ignition
Heat release rate combustion
 Flame speed
studies

5
Ju and Sun: Plasma assisted combustion, Progress of Energy & Combustion Science, 2015
Applications of plasma assisted combustion
Scramjet
engine

Mild Plasma assisted New engine

Combustion combustion technology

Low Fuel/CO2 Cool

Emissions Reforming Flames

6
 2. Plasma Discharges
Plasma: A partially ionized, quasi-neutral charged mixture
in which electrons and ions are separately free.
Frequency: DC, AC, RF, MW, Pulsed…
• Non-thermal (Non-Equilibrium plasma) Tgas ~ 300K-2000K
Cold plasma: Ttrans = Trot < Tvib < Te
Temperature: Low pressure – DC, RF glow discharges
Atmospheric pressure – DBD, Microwave, Corona & Micro-plasmas
• Thermal plasma (Equilibrium) Tgas ~ 2000K-20000K
Hot plasma: Ttrans = Trot = Tvib = Te

Discharge processes: Electric field, Corona, glow, arc

Discharge types: Electron-beam, Corona, Dielectric barrier discharge (DBD),

gliding arc, arc, micro discharge, surface discharge… 7
 Plasma frequency
Perturbation of a neutral plasma
+ - + - + - + + - + - -
+ - + - + - + + - + - - Total charge number: 𝑄 = 𝑒𝑉𝑛𝑒
+ - + - + - + + - + - - 𝑄 𝑒𝑉𝑛 𝑒𝑛
Electric field (between two slabs): E=𝜀𝐴= 𝜀𝐴 𝑒 = 𝜀 𝑒 x
+ - + - + - + + - + - -
Coulomb force on an electron: F=-eE
+ - + - + - + + - + - -
+ - + - + - + + - + - - ε：permittivity
x
𝑑2𝑥 𝑒2𝑛𝑒 If the electron density is 109 cm-3 , the
Equation of electron motion: F=mea 𝑚𝑒 2 = − 𝑥
𝑑𝑡 𝜀 frequency is about 300 MHz. Therefore,
plasma is very fast to restore charge
𝑒2𝑛𝑒
The frequency of electron plasma oscillation is: 𝜔𝑝 = 𝑚𝑒𝜀
=9000 𝑛𝑒 (Hz) neutral properties.

The electron plasma frequency is critical to the propagation of electromagnetic wave in plasma.
If the electromagnetic wave frequency (ω) is less than ωp, electrons in the plasma will response ω
and extracts energy from the electric field and reflect the incident wave. If ω>ωp, electrons in ωp
plasma can not response and the electric field will transmit through the plasma without
reflection. Therefore, for a given ω, there is a critical plasma electron number density (cm-3):
 p 2 me
At microwave frequency of 2.45 GHz, if the electron density is 7.5×1010 cm-3 , ωp > 2.45Gz, ne,c 
microwave will not penetrate to the interior of the plasma, but the plasma surface with reflection. e82
 Mean free path and collisional frequency
Mean free path:
B
u: mean velocity 
8k BT λ
v: relative mean velocity u   vf (u )du  , v  2u
d: neutral particle diameter 0
m d u
A
Cross section area : d 2
The collision number per unit time : d 2 vn Mean free path
traveling distance per unit time u 1
the mean free path    
number of collisions per unit time d vn
2
2d 2n
T=300K, p=1 atm, Molecule diameter: d=3.5A0, the molecule number density:

1.013  10 5
=2.45 10
19
n=p/(kBT)= 1.38  10 23  300  2.5  10 / m
25 3
/ cm 3
1

2 (3.5 10 )  2.5 10
10 2 =0.075 m
25  c  u /   500 / 0.075  10 6  6.7  10 9 / sec
ue 4
For electron neutral molecule collisions in weakly ionized gas, e,n  
n (d 2 / 4)ue d 2 n
d 2 n 8k BTe
Collision frequency  e,n  ue / e,n 
4 me 9
 Debye Shielding and plasma sheath
• Shielding effect: the free charges move towards a perturbing
objective and neutralize the perturbing electric field effect in a
characteristic distance of D. D is the Debye length.

+ +
+ D
+ +
+ Cathode
+ E~0 Sheath

+Q
+ D Potential distribution
+ + +
+ + + +
+ E>0 + - - - - -
+ + +
+ + Neutral plasma

+ Ion bombardment
+ + +

How large is Debye length?

10
Maxwell equation e
E  ( ni  n e )
 E: electric field,
B  0 B: magnetic field
B j: current
E  
t ε: permittivity
1 E μB: permeability
  B  B j 
c 2 t
j  e(n i v i - n e v e )  e(n i - n e ) v  e(n iVi - n eVe ) v i  v  Vi
Mean + relative velocity
Diffusion with an external electric field: diffusion velocity and drifting velocity

niVi   Di   ni  i ni E Di Charge diffusivities,

μi electron and ion mobilities
neVe   De  ne  e ne E
For weakly ionized plasma: ni << n; ne << n
i e
Einstein relation 
Di k BT
e e

De k BTe 11
 Debye length
𝑒
In the Maxwell’s equation: 𝛻 ∙ 𝑬 = (ni-ne) ε: the permittivity of the plasma
𝜀

𝑒
For a steady state problem: 𝑬 = −𝛻𝜑 𝛻𝟐𝝋 = − (ni-ne)
𝜀

𝒆𝝋 1 ni ne
Debye length 𝛻𝟐 𝒌𝑩𝑻𝒆 =− 𝟐(
- ) + +
λ𝑫 n0 n0 + D
+ +
𝒌𝑩𝑻𝒆𝜀 +
λ𝑫 = + E~0
n0𝑒2 +
+ +Q
+ E>0 +
+
In air, if Te = 1000K and n0 = 1013cm-3, we have λ D = 6.9 × 10−5cm. +
+
+ +
+ +
Therefore, plasma is almost quasi-neutral everywhere.

The equation means that the net charge potential will decrease exponentially in a length scale of λ D

12
 Ambipolar diffusion (steady state neutral plasma)
In one-dimensional plasma (zero flow velocity): The ion and electrons fluxes are,

i  niVi   Di   ni  i ni E
e  neVe   De  ne   e ne E
In steady state and quasi-neutral plasma:

( De  Di )ne
i  e E- : ambipolar electric field
ne ( i   e )

Therefore: i  e  - Di ni  i ni E
( De  Di )ni
 - Di ni  i ni In non-equilibrium plasma, Te is much
ni ( i   e )
greater than Ti, the ambipolar diffusivity is
  Dambi ni
much higher than ion diffusivity.

( i De   e Di )  T
Dambi   Di  i De  Di (1  e ) ( i  e )
( i   e ) e Ti 13
Energy transfer in Plasmas

Second
Photoionization electrons

Heating

David Staack, 2016 14

6
Energy transfer in Plasmas: fast heating and vibrational energy relaxation
time-resolved and spatially-resolved measurements of N2 vibrational temperature

Two-step thermalization

Figure 1. Experimental and predicted temperature and N2 Figure 2. Experimental and predicted temperature during and
vibrational temperature during and after a ns pulse after a ns discharge pulse in an H2-air mixture (ϕ=0.14) between
discharge in air between two spherical electrodes 1 cm two spherical electrodes 0.9 cm apart at 40 Torr, plotted together
apart at 100 Torr. with predicted number density of electronically excited N2
molecules and Tv(N2).
Igor V Adamovich and Walter R Lempert, 2015, Challenges in understanding and predictive model development of
15
plasma-assisted combustion, Plasma Physics and Controlled Fusion, Volume 57, Number 1.
 Plasma and Plasma Properties
Ionization processes: Thermal ionization, electron impact ionization, photo-ionization for second electrons
hν + O2= O2+ + e
e + O2= 2e + O2+ +
M+O2= M+O2++ e
Electron quenching processes: recombination and attachment
e + O2= O2- Second
e + O2+=O(3P) + O(1D) 0 electron
x d
Plasma temperature: Electron temperature, vibrational and rotational temperature
Electron temperature: 1 eV = 11600 K = 1.6 ×10−19 Joules.

Equilibrium and non-equilibrium plasma:

Equilibrium: Distribution function: Boltzmann
Q(1, E ,V )   g i exp( i / k BT ) Temperature: Te ≈ Tv ≈ Tn
i.
Non-equilibrium: Distribution function: non-Boltzmann
𝜕𝑓 𝑒
+ 𝒗 ∙ 𝛻𝑓 − 𝑬 ∙ 𝛻𝑣 𝑓 = 𝐶 𝑓 Temperature: Te >> Tv>>Tn 16
𝜕𝑡 𝑚
 A few examples
Te ≈ Tv ≈ Tn Near Equilibrium plasma

arc
Te > Tv>Tn
Tn~n*10,000 K

Non-Equilibrium
plasma
Gliding arc
Te >> Tv>>Tn
Tn~n*1000 K

Corona
Tn~n*100 K
17
 𝑑𝑛𝑒 𝐵𝑝 +
= α𝑛𝑒 α = 𝐴𝑝𝑒𝑥𝑝 − 𝑛𝑒
Electron impact avalanche 𝑑𝑥 𝐸
A,B: constants
E: electric field
𝑛𝑒=eα𝑥 p: pressure
ne: electron number density 0
x d
α: The 1st Townsend coefficient, inverse of net ionization length scale. It is determined
exponentially by E/p or E/N
If α > 0, electron avalanche phenomenon.

Breakdown voltage:
The minimum voltage between two electrodes that causes an arc. At the
breakdown voltage, the rates of ionization and dissociative attachment

Few collision
becomes equal.

Paschen's law: The nonlinear dependence of breakdown voltage is to due to

electron impact avalanche via collisioinal energy transfer
More collisions
The breakdown field for atmospheric air~28.7 kV cm−1

Mean free path (1 atm, air): Molecules: 0.1 µm, Electron-molecule: 5.5 µm
Lieberman, Michael A.; Lichtenberg, Allan J. (2005)
pdmin=1 torr cm at 760 Torr, dmin=13.2 µm, twice of the mean free path.
18
How to produce uniform plasma at high pressure?
 Plasma Discharges
Streamer discharge: a non-thermal narrow filamentary discharge channels formed at the initial stage of a spark
breakdown by a high voltage pulse (1-100 ns). A streamer has a streamer head (space charge) with a high reduced electric
field (~100 Td, 5-10 kV cm−1 for air at atmospheric pressure) followed by a streamer channel with lower electric field and
higher conductivity (charge number density). Formation of a streamer discharge occurs when the electric field in the
streamer head is at the same magnitude or greater than a critical external electric field (4.4 kV cm−1 in air at atmospheric
pressure for a positive streamer, 8–12.5 kV cm−1 for negative streamer). Streamers are fundamental components in many
kinds of discharges such as the dielectric-barrier discharges, corona discharge, and spark. It is widely used in industrial
ozone production, biomedical treatment, plasma assisted combustion, pollution control. Note: a microwave streamer is a
hot plasma not a streamer.

Meek and Loeb criterion: Streamer is formed once the total number of electrons in
Streamer Positive streamer the electron avalanche is so large that their space charge field becomes comparable
cathode anode to E0, the avalanche-to-streamer transition occurs. α𝑑~18 − 20 and ne= 1013 cm-1
E0
Stream propagation: A space charge wave, which can penetrate into neutral gas
ℎν with a velocity much higher than the electron drift velocity, up to a fraction of the
+-
speed of light.
E + -
+ - neutral +
+ -
Space charge Energy balance: Power input by external field: N 𝑞 𝐸 𝑉
0 𝑒 0 𝑠
Streamer head 𝑑𝑛0
Power consumed in ionization: 𝑑𝑡 𝑄𝑒
α𝑑~18 − 20 19
Positive and negative streamers:
Propagation of negative streamer requires
a much stronger space charge field.

Fig. (a) Geometry of the simulation domain. (b) Propagation of growing and decaying positive streamers in an external field of
10 kV cm−1. Both positive streamers are initiated from a Gaussian distributed plasma cloud with a peak density of 1020 m−3 and a
characteristic size σ0 of 0.05 mm. The radius of the spherical electrode Rsph is 0.5 mm. The only difference is that in the left panel
the spherical electrode has a potential Usph = 3.5 kV, whereas in the right panel Usph = 3.2 kV. (c) Propagation of negative
streamers in an external field of 20 kV cm−1. For both negative streamers, the initial plasma cloud has a peak density of 1018 m−3
and a characteristic size of 0.10 mm. The electrode radius Rsph = 1.0 mm, and in the left panel Usph = 4.0 kV, whereas in the right
panel Usph = 3.4 kV
Qin et al., J. Phys. D: Appl. Phys. 47 (2014) 435202 (9pp)

Streamer propagation vs. Flame propagation

Flame front: auto-ignition and diffusive heat transfer, self-supported propagation
Streamer: ionization and space charge transfer, external field supported propagation

20
 Corona Discharge
An discharge around a highly curved conducting electrode induced by a high electric field, but the
external electric field is not high enough to cause a breakdown or arc. Widely used in ozone generation.

Positive corona: electrons are attracted to curved positive electrode and have enough energy to cause
electron avalanche. Electron energy is high, density is low.

Negative corona: ions are attracted to curved negative electrode. The photon emissions via ion-
bombardment on electrode surface cause electron avalanche. Electrons have lower energy but higher
density.

(a) (b) (c)

Fig. Pulsed corona discharge and positive streamer development: CCD photos of the point-wire discharge in air using 5μs optical gate. Applied voltages: (a) at
7.5 kV, (b) and (c) at 12.5 kV. For (a) and (b) the semiconductor switch is used, for (c) the spark gap. The electron temperature is about 5-10 eV. By E M van
Veldhuizen and W R Rutgers, J. Phys. D: Appl. Phys. 35 (2002) 2169–2179 PII: S0022 21
Dielectric barrier discharge (DBD) and NS DBD
40 Torr/AR, NS BDB
A discharge that occurs between electrodes with at least one electrode is covered by dielectric materials. It is a corona
discharge with a dielectric electrode. The existence of dielectric barrier limits the current and restricts transition of DBD
discharge to arcing. DBD discharge often has filamentary micro discharge structures and is physically behaving like an
incomplete streamer breakdown. DBD discharge has low electron number density and high electron energy and been widely
used in ozone generators.

Rectangular quartz channel 22 mm x 10

mm in cross-section and 280 mm in
Nanosecond DBD discharge in air: 20 kV, 10 kHz, pulse N10. Left:
length. Rectangular copper electrodes, 15 Front view; Right: side view
x 60 mm. High-voltage pulses 20 kV on
the high voltage electrode, 25 ns duration Conclusion: NS discharge in DBD geometry in air is non-uniform. Initial electrical
at the half-amplitude, up to 20 KHz. field’s distribution and thermal ionization instability development form the non-
uniform energy distribution in the discharge. This non-uniformity can play a key
Andrey Starikovskiy et al., 2014, AIAA-paper
role in kinetic experiments in this type of the discharge.
22
 DC Glow Discharge (high special uniformity and volumetric)
A self sustained weakly ionized volumetric (non-
filamentary) discharge supported by the
secondary electron emission from the cathode. It
has three distinctive structures: Negative glow,
Faraday dark space, and positive column.

Princeton Plasma Physics Laboratory

http://en.wikipedia.org/wiki/File:Electric_glow_discharge_schematic.png

The glow discharge is stable in a low pressure, but it is possible

to stabilize such a plasma at atmospheric pressure if three
requirements are met: (i) use of a source frequency of over 1
kHz, (U) insertion of a dielectric plate (or plates) between the
two metal electrodes, (iii) use of a helium dilution gas.
23
 Positive column
Faraday dark space
Negative glow
5
FIG. 2. 10 ns exposure time photograph of the gap taken when
the discharge current is maximum. The gap length is 5 mm and
the cathode is located at the bottom.

Cathode is at the bottom

FIG. 1 100 ns exposure time photographs of the gap taken during the discharge initiation, the discharge current being
periodic. The number on the current wave form a) corresponds to the number on the left side of the picture b) and indicates
the time when the picture was taken. The gap length is 5 mm. In each picture, the cathode is located at the bottom.24Francoise
Massines et al., J. Appl. Phys., Vol. 83, No. 6, 15 March 1998
 Atmospheric pressure DC glow discharge

David Staack, Bakhtier Farouk, Alexander

Gutsol and Alexander Fridman, Plasma
Sources Sci. Technol. 14 (2005) 700–711

Figure 2. Images of glow discharge in atmospheric Figure 3. Image of the glow discharge in atmospheric
pressure air at (a) 0.1 mm, (b) 0.5 mm, (c) 1mm pressure hydrogen. Positive column and negative glow
and (d) 3mm electrode spacing are visible. In addition standing striations are visible in the
positive column.
25
 Transition from micro glow discharge to equilibrium arc discharge
•Rotational temperature (Trot) increases with vibrational temperature (Tvib) decreases with increase in pressure.
• Above 100 psi, they are measured to be within 500K of each other which is equal to the uncertainty in Tvib fitting.

David Staack,
26 UTAM
 Spark Discharge
A small volume, high temperature, and high current equilibrium arc initiated by a
high voltage breakdown discharge (~10 ns). It has high current (1-1000 A), low
voltage (10-100 V), and low electron temperature (~1 eV). Spark discharge is
widely used in gasoline engines. The role of spark discharge is to create high
temperature environment for ignition. Laser ignition is also to create a spark.

Plasma torch
Plasma torch is also a continuous electric arc. It is high temperature near
equilibrium plasma. It is widely used in ignition and materials processing.
The temperature, power, and electron number density is very high. It
mostly places a thermal effect in dissociating reactants and accelerating
chemical reactions.

Meghnad Saha derived an equation for the relative number of atoms in

each ionization state in an equilibrium plasma:
Ei 1  Ei
ni 1 2 g i 1  2me k BT 
3/ 2

   e k BT

ni ne g i  h 2

It depends on the number density of electrons, ne. This is because as the number density of
27
electrons increase, the electric field decreases and thus lower the ionization state.
 Gliding arc
a gliding non-equilibrium electric discharges invented by Lesueur et al. [1]. The main distinctive aspect of the gliding
arc is a high level of non-equilibrium with both high electron temperature (1-2 eV) and high electron density as well as
high gas temperature (~2000 K). It can be inexpensively generated under near-to-atmospheric pressures.

Fig. 2 Pictures of the gliding arc plasma system with the (a) side view of central
electrode, (b) top view of system, and (c) time integrated top view photograph of
the magnetic gliding arc creating a plasma disk to quasi-uniformly activate the
flow. The numbers in (a) and (b) indicate the path of the gliding arc from
initiation, point 1, to arc rotation/elongation, points 2 and 3, and final arc
stabilization, point 4.

Fig. 1 Left: Schematic of a traditional gliding arc plasma discharge with

Fig. 3 Short exposure grayscale photograph of the
the numbers corresponding to the sequence in time evolution of the arc as
magnetic gliding arc discharge once stabilized at the
it moves along the electrodes (Ombrello and Ju). Right: direct image of a
largest gap, with the cathode spot (CS) and positive
gliding arc time trajectory [Courtesy from Dr. Z.S. Li at Lund University].
column (PC) shown.
[1]. H. Lesueur, A. Czemichowski and J. Chapelle, Frenchpatent 2 639 172.
[2] A Fridman, S Nester, LA Kennedy, A Saveliev, O Mutaf-Yardimci, Gliding
[3] Ombrello et al., AIAA Journal 2006. 28
arc gas discharge, Prog. Energy Combust. Sci. 25 (2), 211-231
 1   T 
 r (T )    (T ) E
2
Energy conservation equation (1)
r r  r 

Electrical conductivity.  (T ) /  0  exp(  E0 / k BT )

2 2

2
Temperature: T
Effective electric field strength: E0
Conductive arc heat loss per unit length from the solution
T
W  2r (T )  16 (T0 )k BT0 / E0
2

r
From Ohm’s law: V0  RI  Wl / I
, We have: I  (V0  V02  4WlR ) / 2 R

Corresponding to steady and unsteady gliding arc.

W
Electric field: E  2WR /(V0  V02  4WlR )
I

Critical condition: V02  4WlR  0

lcrit  V02 /(4WR), I crit  V0 /(2 R),

Vcrit  V0 / 2, Wcrit  V0 /(4 R),
2

29
 Gliding arc voltage

(a)

Fig. 2 Left: Plot of the increase in electric field in plasma after the transition point in a gliding arc discharge [8]. Right: three sequential
frames of gliding arc images recorded by a high-frame-rate camera, showing the conversion from a glow-type discharge to a much
brighter spark-type discharge [7]. (Courtesy from Dr. Z.S. Li at Lund University)

30
 Gliding arc dynamics and radial production

Short-cut

Fig. 3 Left: A short-cut event recorded at 20 kHz framing rate using an exposure time of 13.9 μs. The short-cut current path is indicated by the arrow in the frame of t
= 50 μs. Right: Three typical single-shot OH PLIF images of a gliding arc using an exposure time of 2 µs, at two flow rates (a) 17.5 SLM, (b) 42 SLM. The typical
thickness of the OH distribution is labelled in the images with unit of centimeters [6, 7] (Courtesy from Dr. Z.S. Li at Lund University)

Magnetic gliding arcs

31
 RF and Microwave discharges
In DC and AC discharges, electrical power is delivered to plasma by moving electrons/ions to the electrodes across the
cathode and anode sheaths. When the electrical frequency is very high like RF and MW, the time required for charge particles
to move across the sheath becomes comparable or longer than the wave period of electrical field . Therefore, the interaction
between electrical field and plasma is exclusively by charge displacement current, not by a directed current to electrodes.
Therefore, it can be delivered without requirement of an electrode in contact with plasma by a sheath.

RF & MW plasma coupling

• Inductive coupling: via oscillating magnetic field
• Capacitive coupling: via oscillating electric field

RF discharge (10k-100M Hz) Microwave discharge (1G-300G Hz)

  p Particle interaction  (ne  7 1010 cm )   p

1
Collective interaction
Low pressure-1 atm Low pressure-high pressure
Field wavelength: meters Field wave length: 12.24 cm at 2.45 GHz
Lower electron energy (1-2 eV) Higher electron energy (5-15eV)
Some sheath No high voltage sheath

32
 Breakdown condition of microwave discharge
dne
Electron production, attachment, and diffusion:  ne ( i  a )  D 2 ne
dt
νi: ionization rate, νa: attachment rate

1 v2 8eTe
Diffusivity of electrons (no Ambipolar diffusion): D  lv  
3 3 c 3 c me
Introducing a diffusion length scale:  : characteristic diffusion length of electrons

dne D
Electron production, attachment, and diffusion:  ne ( i  a )  ne
dt 2
D
( i  a  )t
2
ne  ne 0 e

D
Breakdown condition:  i (E / N )   a (E / N ) 
2
33
 Microwave discharge for ignition and flames

Miles et al., Princeton microwave resonator Ikeda et al., Imagineering Inc.

Qiang Wang et al, APPLIED PHYSICS LETTERS 104, 074107 (2014) 34
 Micro-discharge 1

Microscale Discharge micro tips Microscale Discharges in Liquids

Courtesy by David Staack

Applications
• Largescale surface ignition
• Crude Oil Fuel Reforming
Micro-discharge in CH4/He • Medical Treatments
Microdischarge between ceramic
spheres (Tomohiro Nozaki, 2015)
250 torr (Princeton, 2017) • Plasma catalysts
Electrode
• Aerodynamic Control
• High pressure materials processing

Electrode
Preliminary tests of single (top) and four channels (bottom) micro-discharge using single a RF power supply. 35
The channel is 76mm Χ 26 mm with a gap distance about 0.5 mm, (Princeton, 2017)
Nanostructured discharge
Parallel Plate electrodes
 Dielectrics(Al2O3 of 0.6 mm with AAO
 AAO film of 50 μm, holes of 200nm
 p= 60Torr
 d= 1 inches
 rf power =1~100W

Findings
 α to ϒ discharge mode transition is observed
 Under α discharge mode, dielectric materials has small
effect on the discharge
 Under ϒ discharge mode, the discharge on AAO area is
difficult to transit from α to ϒ mode than Al2O3
Possibility of plasma control using nanostructures

Haibo Mu and Yiguang Ju, 2017, unpublished work 36

 Equilibrium & Non-equilibrium Various Plasmas

Micro dis. Nanosec

Electron temperature, eV
10 Corona
Corona RF DBDDBD
Glow
DC, MW
Gliding Arc
1
Arc/Spark
Arc

0.1 Flame
Flame MHD

1010 1015
Electron number density, 1/m3
37
Ju and Sun: Plasma assisted combustion, Progress of Energy & Combustion Science, 2015
 2. Plasma Assisted Combustion and Applications in Engines
Yiguang Ju
Princeton University

1. PAC in scramjet engines and pulse detonation engines

2. PAC in spark ignition engines
3. PAC in flame stabilizations for gas turbine engines
4. PAC in ignition, combustion, and emission control

Combustion lab.
 Motivation
Development of High Speed/Low Emission Power Systems & Synfuels

•Enhanced combustion efficiency and flame stability

•Predictive new engine design with alternative fuels for low emissions
Princeton University

Advanced Gas Turbines

(Low NOx, after burner, renewable fuels)
 Ignition and Flame Stabilization

Pulse Detonation Engines

 Ignition and DDT

Internal combustion engines:

 Ignition timing control
Conventional discharge?

Combustion lab.
 1. PAC for Ramjets and Scramjets

X43: Hydrogen
Silane: ignition enhancer
Princeton University

Mach 10

X-51 (AFOSR)
• Mach 4-8
• Hydrocarbon fuels
(JP-8, JP-7)
•Flow time scale <1 ms

Flow Residence Time

<< 1
Chemical Reaction Time

Difficulty in ignition, flame stabilization, and

combustion completion Combustion lab.
 Solutions:
Increase residence time or reduce combustion time

•Cavity
Princeton University

•Oblique shocks
Niioka et al.1998

•Plasma
Takita et al.
Kimura et al. 1981
Masuya et al. 1993
Combustion lab.
 Plasma-Assisted Combustion

Ignition & Flame speed ~(a, Dfuel, , Φ, Rc, e -Ea Tf , …)

Princeton University

Plasma discharge
Mixing, flow

Ionic wind
Ions/electrons
O2 +

Temperature Fuel
O, NO Radicals fragments
increase
Excited H2, CH4
species C 2 H4
N2 *(A)
O2 (a1Δg)
Thermal Enhancement Kinetic enhancement Transport enhancement

Combustion lab.
 O Radical Production by Plasma on ignition and flame extinction

Takita and Ju (2006)

Princeton University

He/O2 = 0.45:0.55 He/O2 = 0.38:0.62

Effect of PAC on methane flame ignition and extinction, Sun W. et al. 2010, 2011 Combustion lab.
 Plasma jet ignition enhancement
Electron-impact reactions (CH4, O2, N2, H2)
e + CH4 = CH4++2e
e + CH4 = CH3 + H + e
e + O2 = 2O + e
e + N2 = 2N + e
Princeton University

N+O2 = NO+N
NO + HO2 = NO2 + OH
:

30
z
Nozzle 170
M=1.8
320
Pyrex Glass Window
(a) XH2,PJ = 0.3, PIN = 3.0 kW, PIN total = 4.3 kW
y
Plasma

30
Main Stream x
Jet

Plasma
Fuel Injector
Torch
Feedstock
(b) XH2,PJ = 0.5, PIN = 6.0 kW, PIN total = 8.2 kW
Fig. 2.1 H2 ignition by plasma torch in a M=2.3 flow and the
effect of total heat addition on pre-combustion shock wave, XH2
is hydrogen mole fraction in H2/N2 plasma torch; Pin: plasma
torch electric power, Pin total: total heat addition (Pin+H2 enthalpy
flux) [192]

Effect of mixing ratio of N2/O2 feedstock on wall pressure

increase due to combustion of fuel injected at Xi = 24 mm in
experiment, H2/air mixture, Takita, Abe, Masuya, and Ju,
2006, Combustion lab.
 DBD/plasma jet ignition
Princeton University

Fig. 2.2a Schematic of the test section with

DBD and plasma torch [181]

Wall pressure distributions when H2 fuel was injected

from xi = −40 mm under simultaneous operations of
Fig. 2.2b Direct photographs of DBD plasma in M=2.0 DBD device and PJ torch. (a) Pin = 1.75 kW.
supersonic flow [181] Combustion lab.
(b) Pin = 2.4 kW. (c) Pin = 3.3 kW. (d) Pin = 4.05 kW.
 A nanosecond pulsed plasma discharge, arc jet vitiated
ignition
Princeton University

Fig. 2.4 Schematic of the cavity model [182] Fig. 2.5 OH PLIF images of a cavity flame in supersonic flows
of two different enthalpies: (a) without the plasma and (b) with
plasma at M=2.9, Jn~4 of H2 jet, and (c) without plasma and
Equivalence
(d) with plasma at M=2.6 Jn~3.5 of H2 jet [182]
ratio
0.15 1)
Partially-
0.48 Premixed
Flame
0.93
2) Cavity
1.34 Cavity Flame
Holding
1.95

2.26 3) Non-
Premixed
2.72
Flame
3.51 ISO 200, Exp. 3
ms
N2 arc jet heated air free stream: T0 = 2,500 K, P0 = 1 Combustion
Arc jet vitiated ignition, Do lab.
et al. 2013
Princeton University Gliding arc flame stabilization

Fig. 2.6 Schematic of experimental setup and electrodes arrangement

[184]

Fig. 2.7 Left: discharge without fuel injection.

Right: discharge interaction with H2 injection [184]

large volume and forced ignition Combustion lab.

 Subcritical streamer microwave discharge

The discharge in the undercritical field can be initiated,

for example, at the location of a cylindrical MW
Princeton University

vibrator in the EM beam. Conditions of local electrical

breakdown E ≥ Ecr in this case are realized at the
ends of the vibrator.

vibrator

I. Esakov et al., IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 37,

NO. 12, DECEMBER 2009
Esakov et al. AIAA-paper-2005 Combustion lab.
 Microwave flame stabilization in a high speed flow (200 m/s) as a
preburner

cylindrical MW vibrator
Princeton University

Van Wie et al., AIAA 2006-1212

Fig. 2.9 Schematic of experimental setup [186]

Precombustor Combustion lab.

 Nanosecond plasma ignition in PDE
Princeton University

Fig. 2.22 Left: A valve less PDE setup at the Naval Postgraduate School. This type of architecture requires a
booster and its anticipated applications are missiles or rockets. Right: Comparison of ignition delays for C 2H4/air
mixture using spark plug and transient plasma igniter [178]

7
MSD

Flame-Development Time [ms]

NRP discharge, f = 1-5 kHz
6
NRP discharge, f = 10-40 kHz

0
1 10 100 1000
Total Energy [mJ]

Fig. 2.23 (a) PDE engine facility at the Air Force Research Lab at Wright-Patterson Air Force Base, (b) Schlieren imaging of nanosecond pulsed
discharge igniter in CH4/air mixture, Φ=1, (c) Schlieren imaging of nanosecond pulsed discharge igniter in CH4/air mixture, Φ=0.8 [62]
Combustion lab.
 1.E+2 1400
O atom, 5 kHz O atom, 10 kHz
O atom, 40 kHz T, 5 kHz
0.01 ms 2 kHz 40 kHz T, 10 kHz T, 40 kHz 1200
1.E+0

1000

Mole Fraction

Temperature
0.5 ms
1.E-2
800

1.0 ms 1.E-4
600

400

2.0 ms 1.E-6
200

a.
4.0 ms 1.E-8
0 0.05 0.1 0.15 0.2 0.25
0

Time (ms)

6.0 ms

Figure 6. a) Computed atomic oxygen concentration and temperature

Figure 3: Ignition kernel development for 5 as a function of time with 1% oxygen dissociation repeated at 5 kHz,
pulses of 3.2 mJ per pulse. Left images: 10 kHz, and 40 kHz frequencies for stoichiometric methane-air
mixtures at 850 K. b) . Computed atomic oxygen concentration and
pulse repetition frequency of 2 kHz. Right temperature as a function of time with 0.1%, 0.5%, and 1% oxygen
images: pulse repetition frequency of 40 dissociation repeated at 40 kHz frequency for stoichiometric
methane-air mixtures at 850 K.
kHz. • Increase of ignition kernel volume
• Reactivate chemical radicals before quenching
Lefkowitz, J.K., Guo, P., Ombrello, T., Won, S.H., Stevens, C.A., Hoke, J.L., Schauer, F. and Ju, Y.,
2015. Combustion and Flame, 162(6), pp.2496-2507.
 2. PAC for IC engines

Pressure Sensor Thermocouple

Princeton University

Carburetor

Air Flow Sensor

Exhaust
Fig. 1 2HP FUJI-IMVAC 34 is used
current for both US and Foreign UAVs.

Throttle Rod Choke Rod

PU-Imagineering Inc. Combustion lab.

Princeton University Transient corona discharge

Disk electrode & streamers

Corona enhanced ignition

Gundersen et al., 2003 Combustion lab.

Princeton University

Fig. 2.19 Left: streamers generated by a single 370 mJ, 56 kV, 54 ns pulse
(maximum E/N~400 Td) in air (10 s gate time); Right: flame propagation from
multiple ignition sites at the base of the streamers after a single pulse in F=1.1
C2H4/air mixture (1 ms gate time) [81]

Fig. 2.20 Images of flame development in F=1.1 C2H4/air mixture, 6 ms after ignition. A 300 ms gate
time was used with equal sensitivity for both images and 996×990 resolution. Left: spark ignition
using a standard 105 mJ, 10 ms, 15 kV spark ignition system and a spark plug with a 1 mm gap. Right:
transient plasma ignition using a 70 mJ, 12 ns, 54 kV pulse with a 6 mm gap [81]

1. Shiraishi T, Urushihara T, Gundersen MA. A trial of ignition innovation of gasoline engine by nanosecond
pulsed low temperature plasma ignition. J. Phys. D: Appl. Phys. 2009;42:135208. Combustion lab.
 Microwave and nanosecond plasma assisted ignition
t1 t2

Large
Princeton University

Spark heat
Loss,
Small
volume
Fig.1 Current spark ignition plug: large heat loss, small volume
t1 t2

Microwave Less
Heat Loss

Nanosec pulses Larger

volume
O, OH, NO, C2H4… production

microwave repetitive nanosecond ignition with radical production,

Increased volume, less heat loss
Combustion lab.
 Princeton University

Spark, Microwave, Gliding arc

Microwave antenna
Gliding arc power (Imaging Eng. Inc)
generator
Spark
Nanosecond MGA
pulsed plasma NSD
generator electrodes

Microwave
pulsed power Synchronization
generator pulse
generator

OH* Comparison between Spark and MW ignition.

Imagineering Inc.

Combustion lab.
Princeton University

Ikeda et al., Imagineering Inc.

Combustion lab.
Princeton University

Fig. 2.16(a) direct photograph of plasma assisted 34 cc Fuji engine test setup and (b) the comparison of
limits of stable engine operating conditions with and without microwave (MW) discharge at 2000 rpm [69] .
Lefkowitz, J.K., Ju, Y., Tsuruoka, R. and Ikeda, Y., 2012. A study of plasma-assisted ignition in a small
internal combustion engine. AIAA paper-2012-1133. Combustion lab.
 Microwave/spark ignition
Spark plug Microwave
ignition (Φ) ignition
1 bar >2 1.6
2 bar 1.8 1.0
4 bar 0.9 0.7
6 bar 0.9 0.7
Princeton University

8 bar 0.9 0.7

Table 1. The lean burn limits at different initial pressures

The pressure curve of MW ignition at 8 bar

Q Wang et al., Applied Physics Letters 103, 204104 (2013); doi:

10.1063/1.4830272 Combustion lab.
 Microwave/spark ignition in engine
Princeton University

Fig. 2.12 Comparison of C3H8 flame images

in a compression-expansion engine using
conventional spark plug and microwave
enhanced spark plug, F=1, initial pressure 600
kPa, initial engine speed 600 rpm [172]

Fig. 2.13 The effect of spark ignition and microwave enhanced spark ignition on COVIMEP, fuel consumption and exhaust
emission [172]

Fig. 2.14 the SI and SI+MW modes as a function of equivalence ratio at an initial pressure of 1.08 bar and
300 K (a) for FDT, (b) for FRT [70]

Combustion lab.
Why do we see a significant extension of lean burn in engines with microwave?

Fig. 2.11 Comparison of ignition using spark plug (left), microwave

(middle), gliding arc (right) (Photos were provided by Knite Inc. and
Imagineering Inc.) [197]

Ignition to flame transition

(critical radius, Rc)
10 1
(a), Le=1.2, h=0.0
i Q
?
Flame Propagating Speed, U

g
10 0 a
b
Role of plasma:
j
e

-1
• Mainly increase the initial ignition volume,
10
Q=0.0 Rc; not increase flame speed!
Q=0.1
Q=0.15

h
Q=0.6
• The thermal effect is not very large.
f d c
-2
10
10
-1
Rc 10
0
10
Flame Radius, R
1
10
2

Critical ignition radius

Ignition by nanosecond surface dielectric barrier discharge (SDBD)
S M Starikovskaia

16 (CH4:O2, =1) + 76% Ar

Pressure Pd, atm 14 (CH4:O2, =0.5) + 75% Ar

(CH4:O2, =0.3) + 75% Ar

12

(n-C4H10:O2, =1) + 76% Ar

10

(n-C4H10:O2, =1) + 76% N2

8

No autoignition
6
600 650 700 750 800 850 900 950 1000
Temperature Tc, K
S M Starikovskaia, J. Phys. D: Appl. Phys. 47 (2014) 353001 (34pp)
S.A. Stepanyan, M.A. Boumehdi, G. Vanhove, P. Desgroux, N.A. Popov, S.M. Starikovskaia, Comb. & Flame, 162 (2015) 1336-1349
 Dissociation and fast gas heating via
electronic excitation of molecular nitrogen
N. Popov, 2011, J. Phys. D: Appl. Phys.
44, 285201 (16pp)

N2 + e -> N2 (C3Pu) + e
N2 (C3Pu) + O2 -> N2 + O + O(1D)
2500

2000
Temperature, K

1500

1000 Capillary ns discharge in air,

P=25 Torr, T0=300K
1.5 mm ID, 80 mm length
500 U=10 kV, T=20 ns
O. Dutuit, N. Carrasco, R. Thissen et al. 2012
The Astrophysical J. Suppl. Series, 204/2
0 500 1000 1500 2000
Time, ns 26
 Autoignition vs plasma ignition in RCM
at PTDC=15 bar and TC=970 K, (CH4:O2)+76%Ar

E=0.1-5 mJ, 100 “kernels”

Autoignition Plasma assisted ignition, U=-24 kV
PTDC=14.7 atm PTDC=14.7 atm
40 TC=972 K 40 TC=972 K

30 30
Pressure, atm

Pressure, atm
20 20

discharge
10 10 initiation
(blue step)

0 0
0 50 100 150 200 0 50 100 150 200
Time, ms Time, ms

S.A. Stepanyan, M.A. Boumehdi, G. Vanhove, P. Desgroux, N.A. Popov, S.M. Starikovskaia,
Comb. & Flame, 162 (2015) 1336-1349
27
 Pressure trace and corresponding fast
imaging of flame propagation
40

35

30

Pression / bar
25

20

15

10
discharge
5 initiation

0
200 205 210 215 220
Time / ms

Pressure detector
0.8 ms 1 ms 1.2 ms 1.4 ms

CH4:O2, ER=1 + 70% Ar,

TC=947 K, PTDC=15.4 bar

1.6 ms 1.8 ms 2 ms

S.A. Stepanyan, M.A. Boumehdi, G. Vanhove, P. Desgroux, N.A. Popov, S.M. Starikovskaia,
28
Comb. & Flame, 162 (2015) 1336-1349
 Experiments in n–C7H16:O2:N2
Autoignition (black) and plasma ignition (red)

No discharge No discharge
With discharge With discharge
PTDC=1,6 bar PTDC=2,3 bar
1.6 5
TC=648 K TC=646 K
V = 46 kV 4 V = 46 kV
1.2
Pressure, bar

Pressure, bar
3
0.8
Discharge initiation 2

0.4 Discharge initiation

1

0.0 0
100 200 300 400 100 200 300 400
Time, ms Time, ms

The discharge is able to modify gradually a cool flame

(U increase or P increase) and to initiate a 2-stage flame
29
 Flame Initiation in H2/Air ER=0.5, P=6 bar

Second regime of ignition: Polarity: U>0

Ignition along the
perimeter of HV electrode Energy deposition
W= 4.8 mJ

Quasiuniform ignition around HV electrode. Streamer discharge.

Pressure 6 bar, Temperature 300 K.

S.A. Shcherbanev, N.A. Popov, S.M. Starikovskaia, Comb. & Flame, 176 (2017) 272-284 30
 Flame Initiation in H2/Air ER=0.5, P=6 bar

Third regime of ignition: Polarity: U>0

Ignition along the Energy deposition
discharge channels W= 12 mJ

Ignition along the channels. Filamentary discharge.

Pressure 6 bar, Temperature 300 K.

S.A. Shcherbanev, N.A. Popov, S.M. Starikovskaia, Comb. & Flame, 176 (2017) 272-284 31
Princeton University Controlled plasma discharge for volumetric ignition

(a)

1
cm

Fig. 3.16 Arc produced flow instability and jets [238]

(b)

(c)

Laser ignition and laser guided

discharge control, Miles et al. 2013

Fig. 2.21 Direct photograph of a prototype laser igniter

showing breakdown in air at multiple points [200]
Combustion lab.
 Lean flame stabilization demonstrations
• MINI-PAC Bluff-body stabilized flame (propane or methane, 1 bar, 11 kW)

• LEL: reduced by 10%

• Plasma power = 75 W

• TWO-STAGE SWIRLED INJECTOR (Propane air, 1 bar, 52 kW)

• LEL: reduced
100mm

from 0.4 to 0.11

• Plasma power = 300 W

• AERODYNAMIC INJECTOR (MERCATO, Kerosene/air, 3 bar, 200 kW)

• LEL: reduced from 0.44 to 0.21

• Plasma power = 1 kW
33
Princeton University Nanosecond discharge on fuel lean flame stabilization

Laux et al., 2007 Combustion lab.

 Plasma assisted fuel reforming
Princeton University

400 0C

Ozaki, 22nd ISPC, 2015

Combustion lab.
 Technical questions:
1. Plasma can do a lot of “magics” in combustion
Princeton University

enhancement. Does it really have any “kinetic

merits” on combustion enhancement?

2. How does plasma kinetically enhance ignition, flame

speed, and minimum ignition energy?

3. What are the reaction pathways of plasma assisted

combustion?

Combustion lab.
Lecture 3 Effects of plasma on ignition, flame propagation, burning limits, and
the minimum ignition energy

Yiguang Ju

• The impact of plasma on Ignition and ignition limits

• Flame propagation and the effects of heat loss and stretch

• Extinction, quenching distance, and flammability limits

• The Minimum ignition energy and the critical flame initiation radius

• The effect of electric field on flame propagation speed

 3.1.1 Ignition and ignition limits
Auto-ignition Considering an auto-ignition problem at constant pressure, p, at initial temperature of T0, and fuel
mass fraction of YF0.
ρ:density
dT
Conservation equation C p  BQYF e  E / RT T: temperature
dt Y: mass fraction P, T0, YF0
dY
 F   BYF e  E / RT Q: heat release per unit mass of fuel
dt E: activation energy
T (0)  T0 , YF (0)  YF 0
Adding the mass and energy equations: d
 (C pT / Q Y F)  0
dt

Normalization: C pT / Q Y F C pT0 / Q Y F 0 C pTad / Q  0

  T / T0
q  QY F 0/ C pT0
Tad  T0  QY F 0/ C p  T0 (1  q)
  E / RT 0 C pT0
  t / 0 Y FY F 0C p (T  T0 ) / Q  (1  q   )
Q
 0  ( Be  E / RT ) 1
0

 E / RT  E / RT0 E / RT0  E / RT  E / RT0  ( 1) / 

e e e e e
 d 1
Normalized equation:  (1  q   )e  ( 1) / 
d 
 (0)  1
Asymptotic theory : in the limit of   ,

a small change in temperature will lead to dramatic change in the reaction rate, therefore, in this limiting case, we have
  1  O (1 /  )
T
Define:  1   /    1/ q
d
 qe T0
d
 (0)  0
Solution: tig
   ln(1  q )
Ignition time :  ig  1 / q; exponential growth
2
RC p T0
t ig   ig 0  q 1 ( Be  E / RT0 ) 1  e E / RT0
BQY F 0E
Plasma effects on homogenous ignition (B, E, T):
1. Increase reaction rate B; 2. Reduce activation energy E, 3. change temperature (heat loss or addition)
Auto-ignition with heat loss or heat addition Assume: heat addition or loss is a small perturbation O(1/β):

dT 4R 2 h(T  T0 )
C p  BQYF e d
 E / RT

4R 3 / 3
dt
 qe  H H
3ht0
dY d
 F   BYF e  E / RT RC p
dt
T (0)  T0 , YF (0)  YF 0
 (0)  0

T h>0
Plasma heat addition/loss will shorten/extend the ignition delay time h<0
T0

tig
3.1.2 Plasma chemistry for radical production and heating
Electron impact ionization/dissociation/excitation
>10 eV e +O2 =O++O+2e (R1a) Radical production
~10 eV e +O2 =O+O(1D) (R1b) Excitation
~1 eV e +O2 =O2(1Δg)+e (R1c)
0.2-2 eV e +O2 =O2(v)+e (R1d)

Electron ion recombination, attachment, charge transfer

e+O2+ =O+O(1D) (R2a) Recombination/fast heating
O2+ +O2- =2O2 (R2b) Recombination/fast heating
e+O2 +M = O2- +M (R2c) Attachment
H2O+N2+ =H2O ++N2 (R2d) Charge Transfer

Dissociation and energy transfer by ions and excited species

N2(A,B,C)+O2 =O+O+N2 (R3a)
O(1D)+H2 = OH+H (R3b)
H+ O2(1Δg)= O+OH (R3c)
N++O2= O++NO (R3d)
CH3+HO2(v)=CH2O+OH (R3e)
N2(v=5) +N2 = N2(v=3) + N2 (R3f) Slow heating
N2(v) + HO2 → N2 + HO2(v) (R3g)
 Kinetic ignition enhancement by radiation production

5 Chain initiation
Princeton University

H2 +O2 → H+HO2
H2:O2=2:1 Slow
T=1000 K H2 +O2 → OH+OH
Chain branching and propagation
Ignition time (s)

10
-4
H+O2 → OH + O
O+H2 → OH + H
5
H
OH + H2 → H2O+H
OH

O
Chain-termination
H+O2+(M) → HO2 +(M)
H+OH+(M) → H2O+(M)
-5
10 -7 -6 -5 -4
10 10 10 10
Mole fraction of radicals added into mixture

Combustion lab.
 Kinetic effect by NO production on counterflow ignition
N2

H2 &
N2
Temperature &
16
15
Species Measurements
14

12
13 • FTIR, PLIF, Rayleigh
N2 11 N2
10 8 9
6 7
5
4
Fuel Fuel

3 H2/N2
Diffusion Flame
Air 2 Air

1. Silicon Controlled Rectifier, 2. Silicon carbide

heater, 3. R-type thermocouple, 4. Fuel injection
spacer 5. MGA plasma power supply, 5. MGA device,
6. MGA power supply, 7. Cathode, 8. Anode, 9.
Air/H2/CH4
Magnets, 10. Gliding arc initiation wire, 11. MGA, 12.
Insulator, 13. Nozzle with N2 co-flow, 14. K-type 7
thermocouple & FT-IR probe, 15. Diffusion flame, 16.
Water-cooled nozzle with N2 co-flow.
 Plasma assisted ignition: H2 Ignition by gliding arc
1025
NP + NF NP + 1% H2
1000 NP + 2% H2 P + NF
P + 1% H2 P + 2% H2
Ignition Temperature, K Comp. NP + NF Comp. P + NF
975 Comp. P + 2% H2 Comp. NP + 2% H2

950

925
H+O2+H2O=HO2+H2O
900

875 NO+HO2=NO2+OH
NO2+H=NO+OH
850

825
175 200 225 250 275 300 325
-1
Strain Rate, s

Plasma catalytic effects reduce H2 ignition temperature 8

(Ombrello, T., Ju, Y. and Fridman, A., 2008. AIAA journal, 46(10), pp.2424-2433.)
Plasma can break the conventional explosion limit
HO2+H=OH+OH d[H]/dt→infinity

2k1
1 Not explosive
k2 [M ]

H+O2 → OH + O
1 explosive
H+O2+(M) → HO2 +(M)

OH+O (R1)
H+O2
+(M) HO2 (R2)

• Radical and heat production by plasma can extend the explosion limit.
 Ignition Chemistry: Elementary chain reactions of CH4-O2 system

Chain initiation:
CH4 +O2 → CH3 +HO2
Princeton University

CH4+(M) → CH3 +H+(M) Slow

Chain-branching and propagation

H+O2 → OH + O Opportunity of plasma
CH3 +O2 → CH3O+O
CH3 +O2 → CH2O+OH Slow
CH3 +HO2 → CH3O+OH
CH3O + O2 →CH2O+ HO2
CH3O + M →CH2O+ H+M
CH2O +X →HCO+XH (X=H, OH, O, HO2)
HCO+M→CO+H+M
HCO+O2→HO2+CO
CO+HO2→CO2+OH

Termination reaction
H+O2+M → HO2 +M
CO+OH→CO2+H Combustion lab.
 Kinetics of the ignition: CH4:O2:Ar mixture
(T5=1530 K, n5=5x1018 cm-3)
Autoignition Plasma Assisted Ignition
-1 O2 10
-1 O2 H2O
10
2500 2500
CO2
CH4 CH4

Temperature, K
-2 -2
10 10 H2O

Mole fraction
Mole fraction

Temperature, K
OH,H,CO2,O O CO
10
-3
10
-3 H CH2O
CH3
2000 2000
-4
1x10 CH3,H2O,H2,CO
-4
1x10 OH
CO2
-5 -5
1x10 1x10
H2
-6
10 1500 10
-6
1500
-1 0 1 2 3 4 -1 0 1 2 3 4
10 10 10 10 10 10 10 10 10 10 10 10
Time, s Time, s

Plasma assisted ignition is characterized by:

– slow increase of gas temperature
– developed kinetics of intermediates
– partial fuel conversion during induction time

I N Kosarev, N L Aleksandrov, S V Kindysheva, S M Starikovskaia, A Yu Starikovskii,

Combustion and Flame, 154 (2008) 569-586
11
 Plasma assisted ignition: experiments and
numerical modeling: (CH4-C5H12):O2 + 90% Ar
Auto Exp, C2H6
Auto Calc, C2H6
5 PAI Exp, C2H6
10 PAI Calc, C2H6
, , , C3H8
, , , C4H10 CH4, auto, 0.4-0.7 atm

Ignition delay time, s

4 , , , C5H12
10

3
10
CH4, auto, 2 atm

2
10
C2H6-C5H12,
1
10 CH4-C5H12,
auto, 0.2-0.7 atm
PAI, 0.2-0.7 atm
0
Shock tube/nanosecond dsicharge experiments 10
0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85
-1
1000/T, K
I N Kosarev, N L Aleksandrov, S V Kindysheva, S M Starikovskaia, A Yu Starikovskii,
Combustion and Flame, 156 (2009) 221-233 12
Plasma kinetic effect on CH4 ignition (gliding arc)
1500
Heated Air (Fotache, Kreutz and Law, 1997)
1450 Heated Air (experiment)
MGA (experiment)
1400 Heated Air (model)
MGA (model)

Ignition Temperature, K
1350

1300

1250

1200

1150

1100

1050

1000
150 200 250 300 350 400

Strain Rate, s-1

Plasma catalytic effects reduce CH4 ignition temperature 13
(Ombrello, T., Ju, Y. and Fridman, A., 2008. AIAA journal, 46(10), pp.2424-2433.)
AIAA paper-2007-1025
 3.1.3 Plasma Assisted Combustion: The impact of plasma on the ignition and extinction S-curve
The effect of kinetic enhancement (μs ~ ms, 800-1200 K)

New “S-curve” by Plasma assisted combustion for

small molecule fuel such as H2, CH4

Extinction 15
7x10
Temperature

Plasma generated O2=34%

15
CH4
species: O2=62%

OH number density (cm )

6x10

-3
O, H, O2(a∆g) …
15
5x10
Plasma Smooth
the classical S-curve 15
Ignition 4x10 Transition
15
Extinction
3x10
Residence time plasma S-curve
15
2x10
Scramjet, afterburner
15
1x10
•Strong kinetic enhancement at intermediate temperature Ignition
•Less effect at high temperature 0.05 0.10 0.15 0.20 0.25 0.30 0.35
Fuel mole fraction

Non-thermal plasma dramatically enhances ignition Sun et al. Proc. Comb. Inst. 34, 2010, Combust. Flame 2011, 2012
Ombrello et al. 2008
chemistry, but less impact on flame speed/extinction limit!
Plasma Activated Low Temperature Combustion for large hydrocarbon fuels
Two-stage ignition: n-heptane
1500
Low temperature ignition Hot ignition
Kinetic effect Thermal effect

1200

Temperature (K)
>1100 K
High
H+O2=O+OH
O+H2=H+OH
900 800-1100 K
H2O2=2OH
2HO2=H2O2+O2 Intermediate
HCO+O2=CO+HO2
CH2O+X=HCO+XH 500-800 K
R+O2=RO2
600 RO2→QOOH →R’+OH Low
O2QOOH →R’’+2OH
1 2
300 Large molecules Fuel fragments Small molecules
0.0 0.1 0.2
Time (sec)
More kinetics effect of PAC at low temperature combustion?
 Nanosecond plasma assisted low temperature ignition of
dimethyl ether ignition in a diffusion counterflow flame
P = 72 Torr, a= 250 1/s, f = 24 kHz P = 72 Torr, a= 250 1/s, f = 34 kHz,
XO2=40%, varying Xf XO2=60%, varying Xf

5
5 6x10 increase
6x10 LTC decrease
5
5x10
5
Extinction
5x10 LTC

CH2O PLIF (a.u.)

CH2O PLIF (a.u.)

5
HTC
4x10 5
4x10
5
3x10 5
HTC 3x10
5
2x10 Hot Ignition
5
increase 2x10
5
1x10 decrease
5
0.00 0.02 0.04 0.06 0.08 0.10 0.12 1x10
Fuel mole fraction
0.00 0.02 0.04 0.06 0.08 0.10 0.12
S-Curve Fuel mole fraction
New ignition/extinction curve without
extinction limit
Radical production by plasma can activate LTC at much shorter timescale.
Sun, W., Won, S.H. and Ju, Y., 2014. Combustion and Flame, 161(8), pp.2054-2063. 16
Plasma activated low temperature combustion pathway

H+O2=OH+O
Plasma activated high temperature
combustion pathway

LTC

Plasma activated
low temperature
combustion pathway

O+RH → R+OH
Radical production by plasma
R → R’’+2OH
O+RH → R’’+ 3OH
Flow reactor studies of non-equilibrium plasma-assisted oxidation of n-alkanes

Tsolas, N., Lee, J.G. and Yetter, R.A., 2015. Phil. Trans.
R. Soc. A, 373(2048), p.20140344.
 Ignition enhancement by transient corona discharge

2-10KV, 20-200ns
Princeton University

Disk electrode & streamers

Gundersen et al.
•Increased volume
•Transient discharge

Combustion lab.
 Ignition delay time: corona discharge vs. spark
Princeton University

Radical production
CnHm+e=CnHm-1+H*+e
O2+e=O(1D)+O(3P)+e

Large ignition volume?

Liu J, Wang F, Lee L, Ronney P, Gundersen M. In42nd AIAA Aerospace Sciences Meeting and Exhibit 2004 (p. 837).
Combustion lab.
3.2.1 Adiabatic flame propagation
A flame is a self-propagating auto-ignition and thermal diffusion front. The propagation speed of a one-dimensional
flame front relative to the far field unburned mixture is the flame speed.
Governing equations
Fuel/air ignition
SL T
dT d 2T YF Heat conduction
uC p   2  q , (1)
dx dx Reaction zone
dY F d 2YF
u  D  W (2)
dx dx 2 YF  2 P  Q xf
1st order reaction
T ()  T0 , YF ()  YF 0   BYF e  E / RT
T ()  Tad , YF ()  0


Enthalpy conservation outside diffusion zone:  [ Eq.(1) / q  Eq.(2) / W ]dx



qY F , 0
C pT0 Y F ,0 C pTad Tad  T0 
  C pW
In reaction-diffusion zone: q W q
xf 

 [ Eq.(1) / q  Eq.(2) / W ]dx and neglect convection terms T  Tad


1 Y
x Tad  T0 Le Y F , 0

Tad  T0 Tad  T0 E YF 
Define:  ;  ;  T  Tad (1   /  )
Tad Tad RTad Y F ,0 Le
 d 2YF W
In reaction-diffusion zone (neglect convection):   
dx 2 D

dYF 2 W
Rewriting: d[ ]  2 dYF
dx D

Integrating from flame front to a location x in the reaction zone:

x Y  
dYF 2 W F W YF , 0 Le W YF , 0 Le

xf
d[
dx
]  2 
D YF 
dYF  2
D   f   d  2
D   f   BYF e  E / RT d F

WB f YF , 0 Le 2  E / RTad 0  WB f YF , 0 Le 2  E / RTad

 2
D
(

) e  e d  2
D
(

) e 1

dYF WB f YF , 0 Le 2  E / RTad

 2
2
[ ]x f  ( ) e (1)
dx D 
Integrating from unburned region to flame front to find the fuel concentration gradient at flame front:
xf xf
dY F d 2YF dY F

 u
dx
  [ D

dx 2
 W ]dx uY F ,0 D[
dx
]x f  (2)

WB
Mass burning rate: m 2  ( u) 2  2 Le e  E / RTad U ad   n / 21 , n : reaction order
 2C p

Flame speed is affected by Le, B, E, and T. How does plasma affect flame speed?
3.2.2 Flame propagation speed with heat loss and addition:
dT d 2T   (T  T0 ) /(Tad  T0 ) d d 2 H
uC p   2  q  4 Kp (T 4  T04 ) (1) m   W 
dx dx y Y F /Y F 0 dX dX 2 
dY F d 2 YF X  x / x ref dy 1 d2y
u  D  W (2) m  W
dx dx 2 dX Le dX 2
m  u / U ad
2   (  1) 
 / C p T0 W y exp  
T ()  T0 , YF ()  YF 0 x ref  2 Le 1   (  1) 
dT () dYF () U ad 4KpTad4 
 0 H
dx dx x ref C pU ad C pU ad
t ref 
U ad
Kp : Planck mean absorption coefficient

Outer solution (convection diffusion zone):

Fuel/air
1  exp(mLeX ) X0 YF T
y0  
0 X0
 exp(mX ) X0 Convection-diffusion Reaction-diffusion
0  
1 X0 X=0
 Adding the mass and energy equation and integrate from upstream boundary to flame front:

0 0 0
dy d 1 d 2 y d 2 H dy d 0 1 dy d 0 H


(m
dX
m
dX
)dX   (

Le dX 2

dX 2


)dX m(
dX
m
dX
)   ( 
Le dX dX
)   


)dX

Perturbation: assume heat loss or addition only perturb the temperature and mole fraction in O(1/β)

y  y 0  y1 ,    0  1
0
 1 dy d  H
Rewrite the equation above: m[1 (0 )  1 ()]  

 
 Le dX dX  m
0
 1 dy d  d 0  (0) 
 Le dX  dX   0,
0
 exp  1  or 1 (0)  ln m 2
Using the jump condition across the reaction zone: 0 dX  2 

 1 dy d 
0
H  d 
0
H  d1 
0 Fuel/air
m ln m  
2
     
H 0- 0+
 Le dX dX  m  dX  m  dX  m YF T
Find the perturbation in the burned gas zone:
 d1  H
 dX    Convection-diffusion Reaction-diffusion
  0 m

Flame speed with heat loss/addition: m 2 ln m 2  2 H How does heat loss/addition affect flame?
Extinction limit and flammability limit:

m  u / U ad

Normalized burning velocity, m

1.0
m ln m  2 H
2 2
4KpTad4  0.8
H
-1/2
e Extinction
C pU ad C pU ad 0.6 Limit

0.4
• For a given mixture with a constant adiabatic flame speed, the increase
of heat loss will reduce the flame speed and lead to flame extinction at 0.2
-1
2H=1/e and the normalized flame speed at extinction limit is e-1/2 e
0.0
0.0 0.1 0.2 0.3 0.4 0.5
• For a given heat loss intensity (e.g. Kp), as the mixture fuel concentration
Normalizedheat
heatloss H
loss,2H
decreases, the normalized heat loss H will increases. Therefore, at a Normalized
critical fuel concentration, 2H becomes 1/e, and no flame is available
Fig. The dependence of the normalized burning
below this fuel concentration. This defines the lean flammability limit. velocity on the normalized radiative heat loss
of a one-dimensional planar flame.
• How does plasma can change the flammability limit?
Quenching diameter: Heat recirculation Convection
Heat losses
For a flame propagating into a tube, the heat loss from the flame to the wall

wall is governed by the convective heat transfer to the wall Fuel Uf

Air
Heat loss to the wall  4 f2
2H   Nu wall
Total chemical heat release d2

Flame speed: 4 f2

1.0
m ln m   
2 2
Nu
d2

0.9
Quenching diameter:  4 f2 1

0.8
2H  Nu 
d2 e

0.7
d 0  2 eNu f

0.6
0.1 0.2 0.3 0.4 0.5
2 f / d
Fig. Burning rate (solid line) and normalized flame
d0: the minimum diameter in which a laminar flame can propagation speed U (U=m in this figure) plotted against
Propagate. the ratio of flame thickness to channel width (d) for
selected values of reduced heat transfer coefficient (k)
with in a quiescent, two dimensional channel flow.
(Matalon et al. 2003)
How does plasma discharge affect the quenching diameter?
Effect of chemistry and transport on diffusion flame extinction:
Radical Index and Transport-Weighted Enthalpy
• Extinction of diffusion flames are governed by 600 n-decane
n-heptane
n-nonane
iso-octane

• Fuel loading: fuel concentration (mole/cm3)

n-propyl benzene toluene

Extinction strain rate aE [1/s]

500 1,2,4-trimethyl benzene 1,3,5-trimethyl benzene
n-decane model [24] n-heptane model [24]
• Energy content: heat of combustion (cal/mole) iso-octane model [12, 25]
toluene model [12, 26]
n-propyl benzene model [9]

• Transport (diffusivity): molecular weight of fuel 400

and diluent
• Chemical kinetics: Radical index 300

• Transport-weighted Enthalpy (TWE) 200

• Normalize the energy content and transport

effects 100

• TWE = [fuel]  Hc  (MWf/MWn2)-0.5 Tf = 500 K and To = 300 K

• Universal correlation of extinction limits has been

0
0.02 0.06 0.1 0.14 0.18 0.22
Plasma can
derived with Radical index. [Won et al, CNF 2012] Fuel mole fraction Xf change the
500
n-decane
iso-octane
n-nonane
n-propyl benzene
n-heptane
toluene
500
n-decane chemistry and
n-nonane

Extinction strain rate aE [1/s]

transport
Extinction strain rate aE [1/s]

1,2,4-trimethyl benzene 1,3,5-trimethyl benzene R² = 0.97

n-heptane
400
400 n-alkanes iso-octane
aromatics
n-propyl benzene
toluene
300 1,2,4-trimethly benzene
300
1,3,5-trimethly benzene

200
200

100
100

iso-alkane
Tf = 500 K and To = 300 K
Tf = 500 K and To = 300 K
0
0
0.5 1 1.5 2
0.5 1 1.5 2 2.5 3

[Fuel]Hc(MWfuel/MWnitrogen)-1/2 [cal/cm3] Ri[Fuel]Hc(MWfuel/MWnitrogen )-1/2 [cal/cm3]

Does plasma really kinetically enhance flame speed? Gliding arc on flame extinction experiment
350

300 78W

250 60W
Strain Rate, 1/s

200
44W
33W
150

0W
100
Bundy et al.
Puri & Seshadri
No Plasma
50 33 Watts
44 Watts
60 Watts
78 Watts
0
19 20 21 22 23 24 25 26
Percent Methane Diluted in Nitrogen

1.80E+16
0 Watts, a=83.3 1/s
H2/N2
48 Watts, a=183 1/s
Number Denisty of OH

78 Watts, a=127.7 1/s

Diffusion Flame
Computation
1.20E+16

6.00E+15 0W 48 W 78 W

0.00E+00
-0.4 -0.2 0 0.2 0.4 Air/H2/CH4
Distance Between Nozzles, cm
Role of plasma: mainly thermal effect
Ombrello, T., Qin, X., Ju, Y., Gutsol, A., Fridman, A. and Carter, C., 2006. AIAA journal, 44(1), pp.142-150.
 Plasma Generated Active Species
Lifetime vs. Pressure
ramjets & ICE’s PDE’s gas turbines

Species Lifetime
scramjets

long lifetime
O2(a1Δg)
10-1000 times
more reactive
than O2

1 atm Pressure

What are the effects of O3, O2(a1Δg), O, … on flame propagation?

 Enhancement of flame Speed by plasma generated O3
0.5 14
0 ppm O3 592 ppm O3
1110 ppm O3 1299 ppm O3 12

10
0.4

Enhancement [%]
Slifted [m/s]
8
1299 ppm O3
6
1110 ppm O3
10 mm 0.3
4
Nozzle Tip
592 ppm O3
2
Lifted C3H8/O2/N2 flames
0.2 0
0 0.005 0.01 0.015
Mixture fraction gradient dY F /dR
(~ 1/axial distance)

b
Flame speed extraction S L  Slifted
u
Ombrello, T., Won, S.H., Ju, Y. and Williams, S., 2010. Part I:
Effects of O 3. Combustion and flame, 157(10), pp.1906-1915.
 Kinetic Thermal Enhancement Mechanism by O3
14
SL (O3 decomposing in pre-heat zone)
Flame Speed Enhancement [%]

SL (O3 to O2 far upstream of pre-heat zone)

12
Extrapolated Enhancement (experiment)

Flame Pre-Heat Zone

10 O2
O3 O3
or other stable
O species
8 C3H7+OH
C3H8 Products+H2O
Reactants +HEAT
6
Kinetic, Curvature and Stretch Effects
4
Radical production by plasma
2 O3+O3  O2+O2+O2 mainly accelerates heat release
(not change branching).
0
Kinetic-thermal effect!
0 1000 2000
Concentration of O3 [ppm]
Kinetic Effect by O2(a1Δg) on flame propagation

O2 (a1Δg) + H = OH+O fast

O2 + H = OH +O slow

O2 (a1Δg) at 0.98 eV
O2 (b1Σg+) at 1.6 eV
Lifted flame experimental system ignition Lifted Flame
P2 system

Intensity
w/o O3
camera
φ=1

T3
w/ O3
oxidizer
T2 T1
Wavelength vacuum vacuum fuel
254 nm pump pump

C3H8 or C2H4

O3

O3 FTIR
254nm lower light O3 Hg light O2 P1
O3 O3 microwave
intensity O3 emission NOx power supply
NO
OO33 Ar
3-way
Detector valve
O3
Notch Hg Light OO3
3
Filter Absorption
Cell
Flow

O2(a1Δg) ICOS Cavity

Q(12) Experimental measurement

Flow Flow Q(12) Curve fit

Cavity Enhanced
Absorption (GA)

Cross-Section
(x10-23) [cm2]
PD Off-Axis ICOS Cavity

Mirror Mirror
Lens Lens

Computer Diode Laser

6636.16 6636.20 6636.24
Frequency [cm-1]
 O2(a1Δg) Enhancement of C2H4 Flame Speed
[O2(a1Δg)], ppm ΔHL, mm O2 (a1Δg) + H = OH+O fast
3137 4.76
Far less than O2 + H = OH +O slow
4470 6.82
4627 6.83
Hydrocarbon quenching?
5098 7.31
≈ 5000 ppm O2(a1Δg)  2-3 % Lifted Flame Speed Enhancement

6000
SDO (w/ NO)
Microwave Power = 80 Watts
SDO (w/o NO)
5000 O3 (w/o NO)
Concentration [ppm]

4000 Energy Coupling Into Flow

≈ 1 eV to produce O2(a1Δg)
3000

2000
The kinetic effect of O3 and O2(a1Δg)
1000 on flame speeds is small.
Nozzle Tip
0 Ombrello, T., Won, S.H., Ju, Y. and Williams, S., 2010.
4 5 6 7 8 Part II: Effects of O 2 (a 1 Δ g). Combustion and
Change of Flame Liftoff Height, ΔHL [cm] Flame, 157(10), pp.1916-1928.
Effect of O production in nanosec plasma on
flame extinction
15.24 mm × 22 mm

10 mm

10 mm away from exit

8000
20 & 28 mm ID
6000 f=40 kHz
FWHM= 6 ns
Voltage (V)

4000

2000 f = 5~50 kHz

0
E/N~10-15 Vcm2
-2000

-20 -15 -10 -5 0 5 10 15 20

Power~0.7 mJ
Time (ns) 35
 Atomic O measurement (TALIF) and effect of extinction limit
Atomic oxygen concentration (1015 cm-3)

525

Extinction strain rate (1/s)

14
12 450
10
375
8
6 300
4 no plasma
225 with plasma (f=5 kHz)
computation (T=348 K)
2 computation (T=398 K)
with plasma (f=20 kHz) computation (T=528 K)
0 150 heated flow (T=398 K) computation (T=528 K)
heated flow (T=528 K) 2000 ppm O addition
0 5 10 15 20 25 30 35 40 45
Pulse repetition frequency (kHz) 0.30 0.31 0.32 0.33 0.34 0.35 0.36
Fuel mole fraction Xf
“O production has minor kinetic effect on flame extinction!”
Ar diluted CH4/O2 diffusion flame:
Sun, W., Uddi, M., Ombrello, T., Won, S.H., Carter, C. and Ju, Y.,
XO2=0.28, Peak voltage= 7 kV, P= 60 Torr 2011. Proceedings of the Combustion Institute, 33(2), pp.3211-3218.
36
 Plasma enhances flame propagation of cool flames Flame
1500
500-950 K Tmax~700 K
(a) Cool flame
Hot flame ignition
1200 Ignition

Temperature (K)
(a) Cool diffusion flame
CH2O Tmax~2000 K
900
(b) Hot flame
1300-2200 K
600
N-heptane cool flame ignition (b) Hot diffusion flame

Diffusion Flames Premixed Flames

300
The existence of cool 0.0 0.1 0.2
flame dramatically Time (sec)
change the burning
limits on both lean Won & Ju et al., Proc. Combust. Inst. 35, 2015
and rich sides. Reuter, Won, & Ju, Submitted to 36 symposium, 2016
Ju et al., Combustion and Flame, 2015.
How will plasma
affect cool flames?
When the chemistry
is slow, plasma has a
chance!
 Plasma Assisted Combustion: a multi-disciplinary and multi-physics problem
Plasma discharge

Ionic wind
O2+, N2+ Ions/electrons Instability

Plasma Physics Temperature NO, O3 Radicals Fuel

O, H, OH Int. species fragments
Electric field increase
Joule heating
Electron collision reactions
Excited H2 , CO
N2*, N2(v) species CH4
O2 (a1Δg) CH2O
 Charged species
 Excited species  Ionic wind

Plasma Thermal Kinetic Transport

combustion
Chemical Kinetics studies
Flame Dynamics
Reaction pathways
Reaction rates Traditional
 Extinction Combustion Enhancement
 Ignition
Heat release rate combustion
 Flame speed
studies

38
Ju and Sun: Plasma assisted combustion, Progress of Energy & Combustion Science, 2015
3.2.2 Flame speed, extinction and flammability limits with by flame stretch (Le)
dT d 2T 1 dA f
 C p ax
dx
  2  QYF , 0 Be  E / 2 RTad  ( x  x f )  Qr (T  T0 )
dx Flame stretch: a
 ax
dY F
 D
d 2Y F
 YF , 0 Be  E / 2 RTad  ( x  x f )
A f dt
u  ax
2
dx dx
dT dYF
a  du / dx
(0)  0 (0)  0
dx dx
T ()  T0 , YF (  )  YF , 0
T  T0
  ,
Tad  T0  4K p Tad3 
Qr  4 K p (Tad3  T03 ) H  u=-ax
y  YF / YF , 0 C pU ad C pU ad
a  du / dx
Tad  T0  QYF , 0 / C p X  x / x ref , 
Le 
Tad  T0  QYF , 0 / C p C p D
U ad  Be  E / 2 RT
ad

 : Dirac delta function U ad  Be  E / 2 RT ad at ref

 X,
E Tad  T0  / C p 2
 x ref 
RTad Tad U ad
t ref  x ref / U ad Potential flow (outside)

In the limit of large β

d d 2 2 H 2  ( f 1) / 2
2    e  (   f )  0
d d 2
a  a
Here the stretch rate a is non-dimensional
dy F 1 d 2 yF 2  ( f 1) / 2
2   e  (   f )  0
d Le d 2
a
Perturbation
T  T 0  T 1 /   ..., YF  YF0  YF1 /   ..., Le  1  l /   ... p  T 1  YF ,
1

Outer solution in convection-diffusion zone

 f  f

   e dt /  e dt ,
t 2 t 2
y  1   e dt /  e t dt
t 2
at   ( f , )
2
0 0

   

 0  1, y0  0 at   (0, f )

Jump conditions
d 2 2  ( f 1) / 2
(1) Integrating  e  (   f )  0 from flame front (-) to end of reaction zone (+)
d 2
a


d   2  1f /2
   e 0
 d  a

(2) Integrating the summation of mass and energy equation in reaction-diffusion zone,

 dp dy F 
0

 l  0
 d d  
p  0,   0, and 
Governing equation for perturbed variables:
2 H 0
0
dp d 2 p d 2 yF
2  l  0
d d 2 d 2 a 
  a
Solution: e m   f 2a
pf / 2
,
2 f g1
1 1 2H  2 H
Enthalpy change p f  l (   f  )   f g1 g 2  ( g 3  I 1 / g1 )
2

2 g1 a a
1 (1 n 2 ) 2f
g1   e dn

1 ( n 2 1) 2f
g2   e dn
0
(1 n 2 ) 2f
1 e 1
g3   dn
 n 1
2

(1 k 2 ) n 2 2f
1 (1 n 2 ) 2f 1 e
1
I1   e  k 2  1 dkdn


If  f  ,  2 f g1  1,
2
m 2 ln m 2  2 H
Flame propagation speed with stretch and heat loss: sublimit combustion
2.0 20
• At low fuel concentration (YF=0.029) below the flammability limit,
YF=0.029 YF=0.0295
1.5
15 flame can exist in a narrow range of stretch rate bounded by a
Le=0.9
radiation extinction limit and a stretched extinction limit.
Xf

Xf
10
1.0
• As the fuel concentration increase close to the flammability limit,
5
0.5
there exist two flame islands, respectively, close and away from the
stagnation plane.
0
0.0
0.1 1
a
10 0.1 1 a 10 • As the fuel concentration further increases to slightly above the
25 20 flammability limit, there exist both planar flame at zero stretch rate
Y =0.0296
F
YF=0.03 and a near stagnation flame island at Lewis number below unity.
20
15
Xf
Xf

15 • As the equivalence ratio becomes above flammability limit, the two

10 flame islands merge together and a stretched flame can becomes a
10 planar flame as the stretch rate decreases.
5
5

0 0
0.1 1 a 10 0.1 1 10 Ju, Y. and Minaev, S., 2002. Proceedings of the Combustion Institute, 29(1),
a
pp.949-956.
Dependence of flame location on stretch at Le=0.9
Microgravity experiments

Maruta, K., Yoshida, M., Ju, Y. and Niioka, T., 1996, Symposium
(International) on Combustion (Vol. 26, No. 1, pp. 1283-1289). Elsevier.
Numerical simulation: detailed chemistry
CH4-O2-N2-He

1450
0.469 =0.48
.
h Le=1.4
1400 0.46 .
i
l
T0=1358
. . . .e .f

Flame temperature (K)

1350
k. d
.bc .g 0.469 m
1300
0.46
0.45
1250 .a
.
j
T
DF NSF NF WF
1200
x
0 10 20 30 40
-1
Stretch rate (s )
The G-Curve G-curve(Le < Lecr)
Temperature curve of
1D planar propagating flame Stretched
flame

Standard limit WF

G
A
0 F
D
C

'
0
E(FSWSF limit)
Limit of NSF B
Stretch rate, a
 10
3
The G-curve (Le < Lecr )
CH4/AIR A
Le=0.967

2
Stretch limit of normal flame
10
Stretch rate at extinction (1/s)

, experiment

1
10

B
D Jump limit of weak flame

G
10
0 Φ0
C
Radiation limit
F When a mixture has a low Lewis number, the flammability (Φ0)
of NSF
region can be extended significantly by stretch!
Radiation limit of weak flame
E
-1
10
0.4 0.488 0.6 0.8 1
Ju, Y., Guo, H., Maruta, K. and Liu, F., 1997..JFM, 342, pp.315-334.
Equivalence ratio Guo, H., Ju, Y. and Niioka, T., 2000. CTM, 4(4), pp.459-475.
Φ0
 The K-curve (Le > Lecr )
Now we can understand the experimental
3
10 data on the figure below
A
C3H8/Air
Stretch rate at extinction limit (1/s)

10
2 B

, Experiment

1
10 C

0
10 G

Standard limit F
E
-1
10
0.4 0.6 0.8 1 How does plasma change the flammable region of stretched
Equivalence ratio
flames?
3.3 Plasma effect on the minimum ignition energy and the critical flame initiation radius

Flammability limit?
Internal combustion engine, microwave
Lefkowitz et al. 2012, Ikeda et al. 2009

Spark Microwave gliding arc

Why does a flammable mixture can not be ignited by a spark for a small engine or at lower pressure?
 Puzzle of high altitude relight:
an unresolved ignition problem or a flame problem?

Altitude Engine
instability Flow
[1/p]
speed
Flow speed

Flame speed ~ pn/2-1

Flight speed

Is the flame speed really a problem for relight?

 Ignition spark • What governs the ignition & Eig?
to a flame • What are the chemistry and
transport effects?
Q ? δ

Thermal diffusivity oxygen

Le  
Mass diffusivity Jet fuel

•Eig,min: Defined by flame thickness, δ (make a guess)?

B. Lewis and Von Elbe (1961), Ronney, 2004, Glassman (2008)
4 3 1 1
Eig   C p (Tad  T )  3  3 / 2
3 Su Le
volume heat capacity Larger fuel molecules  smaller Eig

•Eig,min: Defined by stable “flame ball” size?

Zeldovich et al. (1985), Champion et al. (1986)
4
Eig  RZ C p (Tad  T ) ~ Le
3

3 Larger fuel molecules  larger Eig

 Assumptions and simplification:
C ~ 1-1/r • 1D quasi-steady state, Constant properties
Temperature
T* • One-step chemistry
Fuel concentration
• Center energy deposition
T ~ 1/r T T 1  2 T
T• Q t
U 
r r 2 r
(r
r
)  H  r  0, r 2 T / r  Q, Y 0
r  R, T  Tf , Y 0
1
Interior filled
Y Y Le  2 Y
with combustion
Reaction zone U  2 (r )  r  , T  0, Y 1
t r r r r
products

~~
H f0
H  ~~ ~0 ~ ~
Heat & C P S u (Tad  T )
Fuel & oxygen products
diffuse inward diffuse outward Z Tf 1 
  exp    ( r  R)
 2   (1   )T f 
~ ~
rf
r Le: Lewis number
r  ~0 , R  ~0
f f Z: activation energy
Flame speed: effect of flame radius, heat addition and Lewis number σ: density ratio
1  2 ULeR

Z Tf  1  U: flame speed
  T f  Q  Q  R e /  e
 2  ULe
d  exp   Q: ignition energy
Le R  2   (1   )T f  Ω: analytic functions

Chen, Z. and Ju, Y., 2007. Combustion Theory and Modelling,11(3), pp.427-453.
 1
10
The Critical Ignition Radius f i (a), Le=1.0, h=0.0

Flame Propagating Speed, U

0
10
adiabatic a
1.6 (h=0.0)
Flame propagating speed, U

-1
10
1.2
Le=0.5
Q ? 10 -2
h Q=0.00
Q=0.05
Q=0.092
0.8 Q=0.10
0.8 Q=0.20
1.0 g
1.2 e c b
0.4 U=0: 1.4
2.0
10
-3
-2 -1 0 1 2
10 10 10 10 10
Flame Extinction Flame Radius, R
O
ball O O
limit
0.0 O O O
-1 0
10 10 10 1 10 2
Flame radius, R 10 1
(a), Le=1.2, h=0.0
i
1. The critical ignition size and energy is governed by

Flame Propagating Speed, U

two different length scales: g
10 0 a
b
•Flame ball size (small Le) j

•Extinction diameter (large Le) e

-1
10
2. With ignition energy, there is a critical flame Q=0.0
Q=0.1
initiation radius, below which, ignition will fail even Q=0.15
Q=0.6
the mixture is above the flammability limit. f h d c
-2
10 -1 0 1 2
10 10 10 10
Chen & Ju, Comb. Theo. Modeling, 2007 Flame Radius, R
Ignition by heat and radical deposition (qt=0.05)
0
10
10
0
1:
1: q
q cc =
= 0.0
0.0
0 LeZ =2:1.0 Radical Only
10 2: q cc = 0.5
q = 0.5
LeF =3:1.0
qc = = 0.675
0.675 1: q c = 0.0
3: q (b)
c
q t =4: 0.0 (b)
2: q c = 0.4
4: q
q cc =
= 0.7
0.7
5: q = 0.73 3: q c = 0.8
-1 6 5: q cc = 0.73
6: q c = 1.0 4: q c = 1.0
10
10-1 5
10-1 5: q cLe
= 1.2
Le = 2.2
F = 2.2
F
4 Le
LeZZ == 1.0
1.0

UU
3 q
q tt = 0.05
=

U
st 1 0.05
1nd flame
2 flame 0.1
2 1
2
bifurcation
2
1 3
3
-2 5
10 -2

U
10-2 2
10 2 0.05
3
3 4
4
4 4 4
5 4 4 3
5 0
3
3 50
5
0.05
6 0.1
R
0.15 0.2
-3
10 -3
10-3 -1-1
10 -1 0 1 2
1010
10 10
100
10
0 10
101
10
1 10
102
10
2

R
R
R

LeF = 2.2
Critical flame initiation radius

53
Chen et al. 2011
 Minimum Ignition Energy vs. Critical ignition radius:
impacts of flame chemistry and transport
2.5
2
2.4
2.0

Minimun ignition power, Q min 2.3

1.5 Fuel Mean Radical

2.2 1.9 molecular Index
Le = 2.1 weight
1.8 = Le
1 2.0 JP8 POSF 6169 153.9 0.80
1.9
1.8
1.7
Activation energy SHELL SPK 136.7 0.85
1.7 1.6 POSF 5729
0.5 Z = 10
1.6
1.5 1.5 Z = 13
1.4
3
1.4
@ 1 atm
0 2.5
Unburned Temperature = 450 K
0 500 1000 1500 2000 2500 Fuel/Air (21% O2) mixture
3
Cube of critical flame radius, R C

Critical radius [cm]

2

Chen, Burke, Ju, Proc. Comb. Inst. Vol.33, 2010 1.5

0.5 JP8 POSF 6169

SHELL SPK POSF 5729

0
Won, Santer, Dryer, Ju, 2012 0.6 0.7 0.8 0.9 1
Equivalence ratio 
Flame Initiation/Propagation: experimental confirmation
• Experimental setup
• Chamber with high speed Schlieren imaging
• 10 cm radius chamber (2.5 cm flame radius is used)
• Pressure rise less than 3% when Rf = 2.5cm
• Unburned gas temperature = 450 K
• Critical flame initiation radius (stretch rate) and steady state flame
speed can be measured.
• Schlieren visualization / High speed camera (15000 fps)

P TC
Heated Oven
N2/O2/
tube Fan
Vacuum
Liquid fuel pump Pressure
injection release
TC electrodes tank

Vaporization
chamber
TC: Thermocouple Heater Heater
P: Pressure gauge
55
Outwardly Propagation Flames
• Outwardly propagating flames
• n-Decane/Air at ϕ 0.7, 1 atm, 400 K 10
0
1: q c = 0.0
2: q c = 0.5

• Schlieren imaging 15000 fps 1 cm 6

3: q c = 0.675
4: q c = 0.7
5: q c = 0.73
(b)

10-1 6: q c = 1.0
LeF = 2.2
LeZ = 1.0

U
q t = 0.05
1
2
3
-2
10
2
3
4 4 4
5
3 5 6
-3
10 -1 0 1 2
10 10 10 10
R

2.5 180

160

Flame speeds dRf / dt [cm/s]

2.0
Flame radius Rf [cm]

Flame speed 140

1.5

120

1.0 Flame radius

100

0.5
80

0.0 60
0 2 4 6 8 10 12 14 16 18
56
Time [ms]
Flame trajectory / Flame regimes
• Outwardly propagating flame trajectory Rapid rise
5.7 ms
• Flame speed Sb = dRf /dt 220
2 ms

• Stretch rate K = (2/Rf)  (dRf /dt)1 200 (a) Ignition

180 Regime II
• Three distinct flame regimes 160

Flame speed, Sb [cm/s]

140 Critical radius
• Regime I 120

• Spark assisted ignition kernel 100 Regime I Regime III

• Regime II 80
0.0 0.5 1.0 1.5 2.0 2.5
Flame radius, Rf [cm]
• Transition from ignition kernel 220
to normal flame 200
(b)
• Weak flame regime 180 Ignition
• Regime III 160 Linear
extrapolation Regime II
140
• Self-sustained stable 120
propagating flame 100
• Consistent with previous study2 80 0 100 200 300 400 500 600
-1
Stretch rate K [s ]
< n-Decane/Air at ϕ 0.7, 1 atm, 400 K >
• Laminar flame speed / Critical radius
1R.A. Strehlow, L.D. Savage, Combust. Flame 31 (1978) 209–211. 57
2D. Bradley, C.G.W. Sheppard, I.M. Suardjaja, R. Woolley, Combust. Flame 138 (2004) 55–77.
 Critical Flame Initiation Radius
• Experiments have been done
50

45
for JP8, SPK, IPK, and HRJ
tallow at lean conditions. 40

Sl, cm/s
35
JP-8 POSF 6169
30
SHELL SPK POSF 5729

• Reactivity orders from critical 25 HRJ Tallow POSF 6308

radius measurements 20
SASOL IPK POSF 7629

• HRJ tallow ~ SPK ~ JP8 > IPK 2.4

0.7 0.8 0.9 1.0
Equivalence Ratio
1.1 1.2

2.2 JP-8 POSF 6169

• Consistent results to diffusion 2.0 SHELL SPK POSF 5729

HRJ Tallow POSF 6308
flame extinction in TWE 1.8

Rc, cm
SASOL IPK POSF 7629
1.6
1.4
1.2
1.0
0.7 0.8 0.9 1.0 1.1 1.2
Equivalence Ratio
Opportunity: volume ignition by plasma can enhance ignition to
Kim, H.H., Won, S.H., Santner, J., Chen, Z. and Ju, Y., 2013. Proceedings
flame transition and reduce ignition failure. of the Combustion Institute, 34(1), pp.929-936. 58
 Subsonic Ignition Tunnel Utilized
to Elucidate Fundamental Interactions
• Subsonic Wind Tunnel
◦ Premixed methane/air at room temperature and pressure
◦ U = 1 - 10 m/s
◦ Re = 6,000 - 24,000
◦ Optical access through windows on three sides
• Transient Plasma Systems Pulsed Power Supply
◦ 10 ns FWHM
◦ Pulse repetition frequency (PRF) up to 330 kHz
◦ Peak voltage of 10 kV into 50 Ω resistor
◦ Maximum Energy Per Pulse ≈3 mJ

• Electrodes
◦ Lanthanated tungsten
◦ Pin-to-pin configuration
◦ Micrometer controlled inter-electrode gap distance
◦ Tip angle of 20° Courtesy of Timothy Ombrello
59
 Effect of Time Scale of Energy Deposition
Fixed Total Energy and Varying Pulse Repetition Frequency (PRF)

CH4-Air, φ = 0.6, U = 10 m/s, D = 2 mm, and N = 20

Fully-Coupled Partially-Coupled Decoupled

300 kHz 100 kHz 20 kHz 10 kHz 5 kHz 3.3 kHz 2.5 kHz 2 kHz 1 kHz
3.3 µs 10 µs 50 µs 100 µs 200 µs 300 µs 400 µs 500 µs 1000 µs

Three Distinct Regimes Identified 60

Effect of ignition kernel size on ignition probability

Larger ignition size, leaner mixture ignition

J.K. Lefkowitz, T. Ombrello / Combustion and Flame 180 (2017) 136–147

 Effect of Inter-Pulse Time
and Number of Pulses
CH4-Air, φ = 0.6, U = 10 m/s, D = 2 mm, and N = 20

40

30

MIP (W)
20

10

0
0 5 10 15 20
Number of Pulses
MIP = Minimum Ignition Power
(determined for 50% ignition probability)

De-
Fully-Coupled Partially-Coupled coupled

Ignition probability is dependent on PRF (inter-pulse time), not total energy deposition!
• Increasing power deposition rate (high PRF) is a superior method to ensure ignition
• In partially-coupled regime, more pulses increases ignition probability, but not to 100%
• In decoupled regime, ignition probability is a linear function of number of pulses
62
 How Does This Translate to a More Realistic Flow
Implications in a Recirculating Turbulent
Reactive Flow: Mach 2 Cavity
6.6 cm

1.65 cm
M=2
2.54 cm
1.9 cm
Time to Ignition Time to Ignition for
for NPHFD Capacitive Discharge
22.5°
Capacitive Discharge

Steady-State
Chemiluminescence
NPHFD (300 kHz)

63
 Time to Ignition for Lean Cavity (~Φ=0.8)
Energy Deposition of 50-800 mJ

Factor of 7 Difference in Energy

Deposition, But Same Ignition Time
Directly Ties to The Subsonic
Benchtop Experiments to Highlight
Synergy Between Pulses and the
Effect on Flame Growth Rates

Approximately 1
Cavity Cycle Time

Drastic Change in Ignition

Time Below ~ 100 mJ

64
 Summary of plasma ignition in a reactive flow
1. Synergy Between Pulses at High Frequency Changes the Way Ignition is
Approached in a Reactive Flow
2. Time Scale of Energy Deposition Process is Critical
• Can couple into thermal, kinetic, and flow effects

Power Drives Ignition Process, Not Energy

3. Ignition Probability and Flame Growth Rate Can be Enhanced, But There is
an Optimization
• Depends upon the flow velocity, inter-pulse time, reactivity of mixture

4. Connections to Realistic Environment Demonstrated

• Scramjet cavity ignition

More Details of What is Contained in These Slides Can Be Found in the Following References

o J.K. Lefkowitz, T. Ombrello, “An Exploration of Inter-Pulse Coupling in Nanosecond Pulsed High Frequency Discharge Ignition,”
Combustion and Flame, 180 (2017) 136-147.
o J.K. Lefkowitz, T. Ombrello, “Reduction of Flame Development Time Using Nanosecond-Pulsed High-Frequency Discharges in Flowing
Mixtures,” 10th U. S. National Combustion Meeting, April 23-26, 2017 College Park, Maryland.
o T. Ombrello, J.K. Lefkowitz, S.D. Hammack, C. Carter, K. Busby, “Scramjet Cavity Ignition Using Nanosecond-Pulsed High-Frequency
Discharges,” 10th U. S. National Combustion Meeting, April 23-26, 2017 College Park, Maryland.

65
 Summary: The impact of plasma on fundamental combustion properties:
How does plasma assist combustion? Ignition, Flame speed/limit, Emin
Flame speed and propagation Ignition to flame transition
Ignition/extinction S-curve
(Flammability limit) (critical radius, Rc)
10 1
Extinction (a), Le=1.2, h=0.0
Temperature, K

Flame temperature, K
i Q
High ?

Flame Propagating Speed, U

1500 temperature 10 0
g
a

flame b
1200 j
Plasma e

Ignition 10
-1

Plasma Q=0.0
Q=0.1
Q=0.15
Flow residence time Q=0.6

Φ0 Equivalence ratio Φ0,r 10

-2
f h d c
10
-1
Rc 10
0
10
Flame Radius, R
1
10
2

Critical ignition radius

• Increase flame speed • Make ignition kernel > Rc

• Shorten ignition time
• Extend flammability limit • Accelerate ignition to flame transition
• Extend extinction limit

WB
2
RC p T0
t ig  e E / RT0
S L  2 Le e  E / RTad E min  C p (Tad  T0 ) Rc3
BQY F 0E  2 C p
Minimum ignition energy
Ignition delay Mass burning rate
 Summary

1. Plasma has both kinetic and thermal effects on ignition enhancement.

2. Plasma has only minor kinetic effect on flame propagation speed at high
temperature. The main effect for the extension of extinction limit is thermal.
3. Plasma may have strong kinetic effect on cool flame propagation speed and
limits.
4. Plasma can cause fuel fragmentation and reduce the fuel Lewis number, thus
enhance flame speed via the Lewis number effect (Transport).
5. The minimum ignition energy is governed by a Critical Radius. Plasma can
create a large volumetric discharge greater than the critical radius to reduce
the minimum ignition energy, especially at low pressure and fuel lean
conditions.
 Lecture 4 Electric Field Effect on Flames: Ionic wind and Joule heating
Yiguang Ju

Flame is a weakly ionized plasma. It produces ionic wind, electron heating under an electric field

CH + O = CHO+ + e Flame Calcote (1963)

+ - + - + -
weak plasma
CHO+ + H2O = H3O+ + CO ne,ni~1012/cm3
-
DC field Brande (1814)
+ + +
- - - Electron heating:σE2
Ionic wind~10 m/s
+
AC field
Electron heating
+ - + - + -

t combustion  t mw  t e  n
t combustion ~ t AC  t e  n  Ea 
Timescales: S L  exp   1
t combustion  t DC  t e  n 
 2 RT f 

 Ionic wind and Joule heating by an electric field
Flow induced by ion collision with neutral molecules in a flame and
corona with an electric field. The ionic wind velocity can be 10 m/s
which significantly modifies the near electrode flow field.
 
Momentum transfer between: mi vin (Vi  U )
vin : ion - neutral molecule collision frequency

Methane, Φ = 1.0, Air Flow = 30 slm, Fuel Flow = 3.2 slm Flow velocity ≈ 1.0
m/s, Voltage: 0 V vs. t 2000 V, The anode-cathode gap was kept constant
at 40mm, Ganguly et al., 2008
3
Electron-molecule collision energy transfer: E   k i e i ne k B (Te  Tg )
2

2 i
Joule heating: E 2
[1] Robinson, M., 1962, “A History of the Ionic Wind,” American Journal of Physics, Vol. 30, pp. 366-372. 2
 Ionic wind
Calcote, 3rd Symposium on Carleton and Weinberg, Lawton, Mayo, Weinberg,
Combust. Flame, Explosion Phenomena(1948) Nature 330 (1987) Proc. Roy. Soc. A 303 (1968)

F = qE

Min Suk Cha

 Ionic wind: Mechanism
E E : Electric field
F=qE q : Charge • In flames, most of ions are positive ions.
q • Electron mobility is high (smaller mass than ions), its motion
U q : Characteristic ionic wind velocity (m/s)
is reduced by the motion of positive ions.
  qV 
1/ 2
V : Voltage change (V)
• Therefore, ionic wind by negative ions and electrons is
U q     : gas density (kg/m 3 )
  
smaller than that by the positive ions.
 q : Charge density (C/m 3 )

Comparable positive/negative charge carrier Heavier positive charge carrier

O2- H3O+ electron H3O+

- + + + + + +
-
- - - + + + + + +
-
- - - + + + + + +
+ +
- - + + +
- - - - + + + + +
+
+ + +
- - - - + + + + cathode + + + +

cathode
-
anode

+ e

anode
- - - + + + + +
- - + + + + +
-
- e- - + + + + e + + + +
+
- - + +
- - -
- + + + +
- - + + + + + + + +
- - - + + +
- - +
+ + + + + +
- - - - - + + + +
+ + + + +
+
- - - + + + + + +
- -

U U
 Axisymmetric jet flames: transverse DC fields

L = 50 mm
Premixed Non-Premixed
Laser

E
Electrode

Electrode

80 mm
30 mm

6 mm -16kV 16kV -16kV 16kV

L / 2 mm
High Most ions are positive.
Voltage
CH4 or
CH4 – air mixture
Park and Cha, Combust. Flame, submitted (2017)
 Axisymmetric jet flames: transverse DC fields
Premixed Nonpremixed
Vertical Horizontal Vertical Horizontal
H = 2 mm

0 kV
Ground

−8 kV
6 mm

Fig. 2a Fig. 2b

−16 kV
10 mm

E
25 (a) Vertical 1.4 12.5 (b) Horizontal 0.6

14 mm

0 Nozzle

0 0 −12.5 0
−20 Nozzle 20 −20 0 20
[mm] [mm]
Min Suk Cha
 Counterflow nonpremixed flames: DC fields
+ Positive ion
(a) Anode − Negative ion
Electron
Neutral molecule
Movement of ion (+)
E = V/d O02kV
+ N2 Cathode Movement of ion (−)
(b) Bulk flow

− 0.5 kV
(c)

− 1.6 kV
o Drastic change in flow field
o Formation of the dark zone (d)

Fuel + N2
o Flame acts as a source of flow
− 2.4 kV
o Double stagnation planes

Park and Cha et al., Combust. Flame, 168 (2016)

 Counterflow nonpremixed flames: AC fields
Park and Cha et al.
Under preparation
No field 2 kV, 10 Hz

2 kV, 100 Hz 2 kV, 1000 Hz

Propagating edge flames in counterflow: DC

Propagation of nonpremixed edge flames through DC fields

GND O2/inert
(a) (c) (e) Propagating direction
Field line

Stagnation Plane propagating LPF

Propagating direction DF Field directio
Field line direction
RPF YF
–2 kV Ud [cm/s] = 126 kV 162 +2 kV 131
HV fuel/inert

Stagnation Plane

(b) (d) (f)

Tran and Cha , Combust. Flame, 173 (2016)

Propagating edge flames in counterflow: DC

E = V/d o Ionic wind and

secondary flow
modification is the
most important
factor
o Reduced
displacement
speeds
o Rather unaffected
propagation speeds

Combust. Flame 173:114(2016)

 Propagating edge flames in counterflow: AC
Tran and Cha
Proc. Combust. Inst., 36 (2017)
o Wavy motion is
kV Ud [cm/s] = 162 closely related with
fAC.
Stagnation Plane o Ionic wind and
secondary flow
2 kV, 1 Hz Propagating direction 131 modification is the
most important
factor
Field line
o Reduced
2 kV, 50 Hz 135 250 displacement
CH4/O2/N2 = 1/2/5.5 0 kV 1 kV 2 kV
speeds
Ud
o Rather unaffected
Speeds [cm/s]
200 Uedge
propagation speeds
2 kV, 100 Hz 135
150

100
2 kV, 2000 Hz 129 (a)
250
50
C1H /O2/N2 =10
3 8
1/5/14.3 100 1000 Proc. Combust. Inst. 36: (2017)
200 Applied fAC [Hz]
 DC electric field on flame stability

• Ionic wind
• Corona effect
• Instability
• Electron heating
Fig. 2.29 Image sequence of a propane/air flame with an equivalence ratio of 1.2. The applied dc
voltage was slowly increased (left to right), leading to the flame blowing off the burner [65].

Why?
Fuel: Methane
Instability Growth rate: Φ = 1.0
Air Flow = 30 slm
Fuel Flow = 3.2 slm
[n i q e E  (  u   b ) g ] Wisman, D., Ryan, M., Carter, C. and Ganguly, B.,
 k Flow velocity ≈ 1.0 m/s
u  b 0 V - 2000 V (DC)
2008. In 46th AIAA Aerospace Sciences Meeting and
Exhibit (p. 1400).
 Reduction of emission via DBD electric field and discharge
on a diffusion flame
– E-field makes a flame shorter through ionic wind effect
– As soon as a discharge lights up
Princeton University

• No yellow luminosity
– Reduction of soot particles PAH
PLIF
• Onset of PAHs is suppressed by DBD

Soot suppression / Enhanced reaction rate

PAH
PLIF
0 kV 4 kV 6 kV 8 kV 9 kV 11 kV 14 kV
brush corona weak streamer strong streamer

Cha et al., Combust. Flame 141:438, 2005) Combustion lab.

 Effect on microwave electric field on flame speed enhancement
Princeton University

• Microwave frequency is 2.45 GHz

• Three stub tuner to tune the cavity
• Actual Q (5-1000)

Mass flow rate = 5744 st. cm3/min

Exit velocity = 54 cm/s, Equivalence ratio = 0.70
Zaidi, S., Stockman, E., Qin, X., Zhao, Z., Macheret, S., Ju, Y., Miles, R., Sullivan, D. and Kline, J., 2006, In 44th AIAA
Aerospace Sciences Meeting and Exhibit (p. 1217). Combustion lab.
 Estimated vs Experimental Results for Laminar Flame Speed Enhancement
Princeton University

1.1
mm
MW Off MW On
Sref = 29.6 cm/s Sref = 35.7 cm/s

Fraction of eletron heating and normalized heat relase

1.0 50

Electron number density (10 1/cm )

3
In flames, microwave field qr/qr,max
0.8 40

10
mainly heats the electrons
CH4-air
and raises flame temperature =1.0
0.6 30

ne
0.4 20
qe0-x/qrt

0.2 10

0.0 0
0.0 0.2 0.4 0.6 0.8
X (cm) Combustion lab.
 Planar FRS Measurements
30 W Pulsed Microwave Enhancement
30 kW-peak, 1 ms, 1000 Hz Pulses

CH4/air, f = 0.76
Princeton University

• Observed Flame Speed Effect

Flame Shifted
Coordinates – Increase of 0.25 L/D units  ~6% flame
speed enhancement
• 80 K in post flame temperature (CW
saw ~150 K)
• Energy deposited within a few mm of
reaction zone
Temperature [K]

No MW Pulsed MW
80 K  12 W energy deposition  50% magnetron coupling efficiency

Combustion lab.
Stockman, E.S., Zaidi, S.H., Miles, R.B., Carter, C.D. and Ryan, M.D., 2009. Combustion and Flame,
16 156(7), pp.1453-1461.
 Combustion enhancement by a gliding arc (Joule heating and kinetic effect)

Non-Equilibrium (a)
4
Princeton University

Critical
Point
3

Near
Equil.
2
(Yardimici et al. 1999).

(-) 1 (+) R
Breakdown/
Arc Initiation

Combustion lab.
 Short cut event and OH measurements in a gliding arc (air)
Princeton University

Short-cut

Left: A short-cut event recorded at 20 kHz framing rate using an exposure time of 13.9 μs. The short-cut current path is
indicated by the arrow in the frame of t = 50 μs. Right: Three typical single-shot OH PLIF images of a gliding arc using an
exposure time of 2 µs, at two flow rates (a) 17.5 SLM, (b) 42 SLM. The typical thickness of the OH distribution is labelled
in the images with unit of centimeters (Courtesy from Dr. Z.S. Li at Lund University)

The combination of thermal heating and radical production (high electron density and high electron energy) of a
non-thermal gliding arc can enhance ignition and flame stabilization in both thermal and kinetic ways.
1. Sun, Z.W., Zhu, J.J., Li, Z.S., Aldén, M., Leipold, F., Salewski, M., Kusano, Y., Optics Express. 2013, 21 (5) 6028-6044.
2. Zhu, J., Sun, Z., Li, Z., Ehn, A., Aldén, M., Salewski, M., Leipold, F., Kusano, Y., Dynamics, Journal of Physics D: Applied Physics. 2014, 47 (29)
295203. Combustion lab.
 Summary: Electric field effect on flames
1. Ionic wind
• Low frequency electric field generates ionic wind flowing to both electrodes from a flame due to positive and negative charge
carriers.
• Ionic wind can reduce soot/NOx formation due to the change of mixing and flame temperature.
• Ionic wind may induce flame instability due to the force field.
• Ionic wind also modifies flame speed and reduces flame temperature due to increased heat losses from the flame zone.

2. Joule heating
• Electric field generates Joule heating in the flame zone and at the downstream of the flame.
• The electron Joule heating can enhance flame speed via the increase of flame temperature.
• Microwave Joule heating in flames is not energy efficient because much of the energy absorbed by in the burned gas.

3. Radical production by strong electric field

• When the electric field is above the breakdown threshold, a gliding arc or corona can produce radicals to enhance ignition via kinetic
pathway.
• A gliding arc has high temperature and high electronic energy and density, which lead to both thermal and non-thermal enhancement
effects on flames.

It is necessary to understand the kinetic effect of non-thermal plasma at high E/N on combustion
Lecture 5 Chemistry and Kinetic Studies of Plasma-Assisted Combustion

Yiguang Ju

• Important chain-initiation and branching reactions in combustion

• Plasma chemistry and timescales
• Impact of plasma chemistry on combustion
• Diagnostics of plasma properties and chemistry in PAC
 1. Important combustion reactions
1. Chain initiation and propagation reactions
RH+ O2 → R+HO2 High Temperature (>1100 K) slow
RH+HO2 → R+H2O2 High pressure/low temperature (>550 K) slow
R+HO2 → RO+OH High pressure/low temperature (>550 K) slow

2. important branching reactions at different temperatures

H+ O2 → O+OH High Temperature (>1100 K) Fast
H2O2 → 2OH Intermediate temperature (800-1100 K) Slow
R → RO2→QOOH → O2QOOH →R’’+2OH Low temperature (300-800K) Slower

Plasma assisted combustion:

O+RH → R+OH
e+O2 → e+ 2O
R → R’’+2OH
O+RH → R’’+ 3OH Faster

e+O2 → e+ O2(a1Δg) H+O2(a1Δg) → OH+O Faster

Plasma provides new reaction pathways to accelerate chain reaction processes

 Interaction of plasma chemistry with reaction kinetics of
large alkanes w/wo in plasma assisted combustion
Fuel(RH)
+OH Reaction rate
RO2
+O2 Small C2H3/CH2O Transition state theory
R alkene
+O2 k1 k2
RO2* QOOH HO2 H/HCO A  B  AB *  C  D
+O2+(M)
R*, O(1D)

O2QOOH H2O2 CO/CO2 k BT q AB * * E *  E A B

+O2 k (T )  exp( )
h q A * qB * k BT
2OH
Plasma
How does plasma affect elementary rate
I e, R*, N2*, O2*,O*
II O2(v), O2(a1Δg) constant?
RH(v), R(v), N2(v), O2(v), HO2(v)
e.g. at 800 K
A schematic of the key reaction pathways for high
pressure fuel oxidation of at different temperatures O2 (a1Δg) + H = OH+O Fast
(blue arrow: Below 700K; yellow arrow: 700-1050 O2 + H = OH+O Slow
K; red: above 1050K). Green: plasma activated
pathway CH3 +O2(v) → CH2O+OH Fast
CH3 +O2 → CH2O+OH Slow
 2. Plasma chemistry and timescales of
kinetic processes

ttr telec
Ion/Molecule
Kinetics
tfp trot
Ion-Ion, Ion-Molecular

EEDF Excitation /
Quenching
Electron Kinetics
tfp Ionization Recombination

Combustion
ttr tvib Molecules

tfp
Fig. 1.5 Schematic of timescales and key kinetic
Processes
trot Radicals pathways at different stages of plasma assisted
ignition and combustion.

10-14 10-12 10-10 10-8 10-6 10-4 10-2 s Ju and Sun, PECS, 2015

Courtesy of Andrey Starikovskiy

 Potential Energy Curves of O2
O2(B3Su-), 8.4 eV
smax = 1.0 A2 (9.4 eV)
DE ~ 1 eV
O2(3Pg), 5.6 eV
smax = 0.16 A2 (12 eV)
DE ~ 1.5 eV
O2(A3Su+), 4.5 eV
smax = 0.18 A2 (6.6 eV)
E, eV

O2 (a1Δg) O2 (b1Σg+) at 1.6 eV

O2 (a1Δg) at 0.98 eV

r, nm
Electron impact reaction is a function of electron energy distribution (E/N)
 Electron impact reaction cross sections-O2

1. effective 7. O2(a1) 8. O2(b1) 9. O2(A3Su+), 4.5 eV

2. rotational excitation 10. O+O 11. O+O(1D) 12. O+O(1S)
3-6. O2(v1) - O2 (v4) 13. O2+
 Potential Energy Curves of N2

N2(A3Su+), 6.2 eV
smax = 0.08 A2 (10 eV)

N2(B3Pg), 7.35 eV
smax = 0.20 A2 (12 eV)
E, eV

N2(C3Pu), 11.03 eV
smax = 0.98 A2 (14 eV)

Threshold energy diagram

r, nm Electron impact reaction is a function of electron energy distribution
Energy Transfer of non-equilibrium excitation in
Plasma Discharge
N2:O2:H2 = 4:1:2
1
N2(el)
N2(v)

H2(v) H2(el) ion

Energy loss fraction
0,1 H2(rot)
Rot+tr
O2(dis)
O2(v)

0,01
O2(4.5 eV)
O2(a+b)

O2 (a1Δg)
1E-3
1 10 100 1000
E/N, Td

Physics of Nonequilibrium Systems Laboratory

 Influence of Electronic and Vibrational Excitation on
Combustion Kinetics

N2:O2:H2 = 4:1:2
1
N2(v)
N2(el) N2 + e = N2(C3) + e
H2(v) H2(el) ion N2(C3) + O2 = N2 + O + O
Energy loss fraction

0,1
Rot+tr
H2(rot)
O2 + e = O + O + e
O2(dis)
O2(v)

0,01
O2(a+b)
O2(4.5 eV) N2 + e = N2(v) + e
N2(v) + HO2 = N2 + HO2(v)
1E-3
1 10 100 1000 HO2(v) = O2 + H
E/N, Td

Physics of Nonequilibrium Systems Laboratory

 Influence of Vibrational Excitation on Low-
Temperature Kinetics: H2O2 Decomposition

Measured and calculated OH decay time. P = 1 atm.

a) 3%H2 + air; b) 0.3%C4H10 + air.
Physics of Nonequilibrium Systems Laboratory
 PRINCETON
University
Effect of “Hot” Atoms on Active Species Production
in High-Voltage Pulsed Discharges
Nonequilibrium distributions of neutral species are formed in different physical situations.
In laboratory experiments and in the terrestrial atmosphere, there are numerous collisional
processes in which translationally energetic (superthermal) atoms with energies much
above thermal energies are produced.

Potential energy curves Momentum transfer cross Cross sections for scattering of H
and hot atoms formation section for the H-H2 scattering atoms with H2, O2, CH4 and N2
[15]
[33]
[23] (a)
This work

2
cm

2
cm
-16
Cross section, 10
1
10

-16
Cross section, 10
10

H-H2 (el)
H-O2 (el)
H-CH4 (el)
H-N2 (el)
0
10 H+O2=O+OH (new)
H+CH4=CH3+H2 (new-1)
H+CH4=CH3+H2 (new-2)

1
0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5 6,0 6,5 7,0
0,1 1
Energy, eV
E, eV

Direct electron-impact dissociation

e + O2 → e + O2* O2* → 2O(3P,1D) + 1.3 eV
e + H2 → e + H 2* H2* → 2H(1S) + 4.5 eV
e + CH4 → e + CH4* CH4* → CH3 + H + 3.5 eV
 PRINCETON
University
Effect of “Hot” Atoms on Active Species Production
in High-Voltage Pulsed Discharges
Average energy of H atoms Species produced during energy Plasma-assisted oxidation
in various gaseous mixtures degradation of one “hot” H atom in H2-O2 mixture
H CH3 1.8 2H2-O2 mixture
(a) 1

Species amount per one hot atom

T=300 K; P=1 atm
1.6
H2 [H2O]
OH

Production, ppm/10 eV/molec

0 1.4 [H2O2]
10

CH4:O2:N2=1:2:8 -1
10 O 1.2
Energy, eV

-5
1.0
CH4:O2=1:2
0.8
HO2
-2 H2O 0.6
10
-1
H2:O2=1:9 10
H2:O2=2:1 0.4

0.2

-3
10 0.0
0 10 20 30 40 50 0.4 0.5 0.6 0.7 0.80.9 1 2 3 4 5 6 7 8

s
)
D

(h

al
at
Number of collisions Initial H atom energy, eV

(1

(1
H

ic
St
O

ad
d
)+

R
un

o
(h

ro

N
H

G
Plasma-assisted oxidation
 Analysis of the effect of formation of "hot" atoms with excessive
in CH4-O2 mixture translational energy shows the important role of these processes
0,9

0,8
CH4-2O2 mixture
T=300 K; P=1 atm
in formation of active radicals.
[H2]
 The density of radicals produced in discharge plasma can be
Production, ppm/10 eV/molec

0,7 [H2O]
[H2O2]
0,6
[CH2O] several times higher than that produced in the absence of high-
[CH3OH]
-5

0,5
[CH3O2H] energy atoms.
0,4
 The effect plays a fundamental role in the formation of the initial
0,3

0,2
distribution of active species in combustible mixtures and can
0,1 greatly influence the kinetics of ignition and oxidation at low gas
0,0 temperatures.
)

s
)
D

(h

al
at
(1

(1
H

ic
St
O

ad
d
+

R
un
)

o
(h

ro

N
H

G
Gas heating at high E/N E/N = 103 Td

50

Fast Gas Heating

Fractional power, %
40
Electron-ion recombination e + O2+ → O + O* + ΔE
15 -3
Ion-ion recombination O2 + O2 + M→ 2O2 + M + ΔE
- + ne0=10
14
cm 1 atm
-3
ne0=10 cm 1 atm
Collisional energy transfer 30 15
ne0=10
-3
cm 300 Tor
Electronically-excited species
N2(A,B,C,a) + O2  2O + DE 14
ne0=10
-3
cm 300 Tor

Hot atom and molecule O2* → 2O(3P,1D) + DE 20

0,0 0,1 0,2 0,3 0,4 0,5
Mole fraction of O2

Slow Gas Heating Fractional Electron Power Transferred Into

Heat in N2:O2 Mixtures
Vibrational energy relaxation N2(v) + M = N2(v-1) + M +DE
High oxygen, faster gas heating!
 Princeton Plasma Combustion Kinetics
Major Pathways
Ar O2 N2 H2 CxHyOz

Ionization by electron impact. k = f(E/N)

Ar+ O2+ N2 + H2+ CxHyOz+,…, CxH1Oz+

Ar, N2, O2 O2, CxHyOz N2, O2, CxHyOz O2, CxHyOz CxHyOz
H2, CxHyOz H- transfer
Charge transfer, negative and complex ions formation

Ar2+, N4+, O4+, N2O2+, NH2+, H3+, HO2+, H3O+; O-, O2-, O3-, O4-; CxHyOz+,…, CxH1Oz+
Electron-ion recombination Ion-ion recombination Molecule-ion reactions
O2+, O4+, CxHyOz+ O2- + N2+; O2- + CxHyOz+ O2- + H; O- + H2
electron detachment
Electronically- “Hot” atoms and
excited molecules formation Ionic chains
particles formation
Oh, Hh, Nh, O2h, H2h
Low-Temperature
Fast Gas Heating O(1D), O(1S), N(2D), H(n=2) Reactions
Andrey Starikovskiy
 Princeton Plasma Combustion Kinetics
Major Pathways
Ar O2 N2 H2 CxHyOz

Vibrational levels excitation by electron impact. k = f(E/N)

N2(vib) H2(vib)

VT relaxation Reactions of vibrationaly excited molecules

N2(v) + O; N2(v) + H2 H2(v) + O → H + OH(v)

N2(v) + H2O; N2(v) + CxHy H2(v) + OH → H2O + H

Slow Gas Heating Formation of vibrationaly-excited products

Energy transfer to reagents Energy transfer to buffer

N2(v) + HO2 → N2 + HO2(v) OH(v) + N2 → OH + N2(v)

Reactions of vibrationaly excited molecules

HO2(v) → H + O2 OH(v) + H2 → H2O + H

Typical plasma reactions for radical production and heating
Electron impact ionization/dissociation/excitation
>10 eV e +O2 =O++O+2e (R1a) Radical production
~10 eV e +O2 =O+O(1D) (R1b) Non-equilibrium excitation
~1 eV e +O2 =O2(1Δg)+e (R1c)
0.2-2 eV e +O2 =O2(v)+e (R1d)

Electron ion recombination, attachment, charge transfer

e+O2+ =O+O(1D) (R2a) Recombination/fast heating
O2+ +O2- =2O2 (R2b) Recombination/fast heating
e+O2 +M = O2- +M (R2c) Attachment
H2O+N2+ =H2O ++N2 (R2d) Charge Transfer

Dissociation and energy transfer by ions and excited species

N2(A,B,C)+O2 =O+O(1D)+N2 (R3a)
O(1D)+H2 = OH+H (R3b)
H+ O2(1Δg)= O+OH (R3c)
N++O2= O++NO (R3d)
CH3+HO2(v)=CH2O+OH (R3e) Slow heating
N2(v=5) +N2 = N2(v=3) + N2 (R3f)
N2(v) + HO2 → N2 + HO2(v) (R3g)
What are the major species produced by plasma?
Time

Pressure

•Long lifetime species? NO, O3, O2(a1Δg)

•Short lifetime plasma generated species? O, N2 (A,B,C)*
17
 3. Impact of plasma chemistry on combustion
Question: When will electron impact dissociation process become important in combustion?

Fig. 3.5: Rate constants (a) and reaction flux (b) for reactions for dissociation by electron
impact at electric field values equal to 200 Td and 500 Td and chain branching reactions.

Ju and Sun, PECS, 2015.

Comparison of the reaction rates of electron impact and excited species for radical production

(Ground) (Ground)

Important radical
production channels

(Ground)

S M Starikovskaia, J. Phys. D: Appl. Phys. 47 (2014) 353001

A. M. Starik, B. I. Loukhovitski, A. S. Sharipov and N. S. Titova, 2016, Phil. Trans. R. Soc. A 373: 20140341
 PAC: how does plasma change the branching reactions in combustion?

1500
Low temperature ignition Hot ignition
Kinetic effect Thermal effect

1200

Temperature (K)
>1100 K
High Temp.
H+O2=O+OH
O+H2=H+OH
900 850-1100 K
H2O2=2OH
2HO2=H2O2+O2 Intermediate
HCO+O2=CO+HO2 Temp.
CH2O+X=HCO+XH
600 R+O2=RO2 500-850 K
RO2→QOOH →R’+OH Low Temp.
O2QOOH →R’’+2OH
t1 t2
Fuel fragments Small molecules
300 Large molecules
0.0 0.1 0.2
Time (sec)

Schematic of kinetic and thermal enhancement pathways of plasma assisted combustion for liquid fuels at high, intermediate, and low temperature, respectively

Y. Ju and W. Sun, Prog. Energy Combust. Sci., 2015

Plasma activated Cool Flames :A new way to burn with plasma
1 n-heptane


Ignition delay time (s)

1 atm
5 atm
10 atm
20 atm
0.1
Plasma

0.01
Extinction
Temperature

0.8 1.0 1.2 1.4 1.6

Plasma 1000K/T
HTC generated

LTC
LTC Ignition

t2 t1 Residence time
t2<< t1

Plasma activated LTC at much shorter

time, lower pressure….

We can create cool flames even at

1 atm or below? 21
 Observation of plasma activated self-sustaining Cool Flames
Fuel/N2 @ 550 K Tf~1900 K Tf~650 K
Heated N2 @ 550 K

Stagnation
plane (a) Hot diffusion flame (b) Cool diffusion flame
N2 @ 300 K Fig. 2 Hot and cool n-heptane
diffusion flames at the same condition
Oxidizer @ 300 K
with plasma discharge 2400
nC7H16/N2 vs O2 or O2/O3
HF branch

Maximum temperature Tmax [K]

in counterflow burner
Fig. 1 Schematic of experimental setup Xf = 0.05,Tf = 550 K, and To = 300 K
2000
Extinction limit of
conventional hot diffusion flame
(HFE)
1600
without O3

Won, S.H., Jiang, B., Diévart, P., Sohn, C.H.

Transition to
and Ju, Y., 2015. Proceedings of the
1200
hot flame Extinction limit of
Combustion Institute, 35(1), pp.881-888. with O3 cool diffusion flame
HTI (CFE)
800
CF branch

LTI

400
Extinction/instability
0.1 1 10 100 1000 10000
Strain rate a [s-1]
 Plasma assisted Self-Sustaining
Premixed/partially premixed Cool Flames

• Cool diffusion flames

– n-Heptane/O2/O3
– Won et al., Proc. Combust. Inst. 2015

• Cool premixed flames

– DME/O2/O3
– Reuter et al., Combust. Flame 2016

• Cool partially premixed flames

– DME/O2/O3
– Reuter et al., Proc. Combust. Inst. 2017
(hopefully)
S. H. Won et al., Proc. Combust. Inst. 35 (2015) 881-888
23 C. B. Reuter et al., Combust. Flame (2016), in press
 4. Diagnostics of plasma physics and chemistry in PAC

1. Measurements of plasma properties and kinetic processes

2. Plasma assisted ignition and combustion with active species production

3. Kinetic studies of plasma assisted combustion

1. Measurements of Plasma Properties: electron density and temperature Thomson scattering
𝑦 Power of scattering: 𝑃𝑠 ∝ 1 − sin2 𝜃 cos 2 𝜙0
𝐸𝑖0 1
𝐸𝐿 𝑑𝜎𝑒 2𝜆0 𝜃 2𝑘𝐵 𝑇𝑒 2
Number of photo-electrons: 𝑁𝑠 = Δ𝐿 𝑛𝑒 𝜂 Δ𝜆1 = sin
ℎ𝜈0 𝑑Ω 𝑒
𝑐 2 𝑚𝑒
𝐴𝑒 𝑑𝜎𝑁2 2
𝐴𝑁2 𝑑Ω 𝑐 2 𝑚𝑒 Δ𝜆1 𝑒
𝑛 = 532𝑛𝑚
𝑛𝑁2 𝑓𝐽=6 𝑇𝑒 =
𝜙0 𝑘0 𝑒 𝑑𝜎𝑒 8𝑘𝐵 sin2 𝜃 2 𝜆0
𝜃 𝑑Ω
𝑘𝑠
𝑧 k0: Laser beam direction, ks: Scattering signal wave vector
𝑥 𝐸𝑖0 : Polarized electric field, scattering is rotationally symmetric about 𝐸𝑖0 .
x-z plane: the plane of observation
θ: the scattering angle relative to the laser beam.
𝜙0 : angle between observation plane and the polarization angle.
𝐴𝑒 and 𝐴𝑁2 : integrated intensities of the Thomson and Raman spectra
𝑑𝜎𝑒 𝑑𝜎𝑁
and 𝑑Ω2 the Thomson and N2 Raman scattering cross sections
𝑑Ω
𝑓𝐽=6 : the fraction of N2 molecules in the J = 6 rotational state
Δ𝜆1 𝑒 : the half 1/e width of the Gaussian broadening profile
EL: laser energy, η: optical efficiency, ΔL: length of observed scattering segment.
• H. Van der Meiden, "Thomson scattering on low and high temperature plasmas", Ph.D, Technische Universiteit Eindhoven, 2011.
• A. Roettgen, "Vibrational Energy Distribution, Electron Density and Electron Temperature Behavior in Nanosecond Pulse Discharge Plasmas by Raman and
Thomson Scattering", Ph.D, The Ohio State University, 2015.
Thomson Scattering Experimental Setup and Calibration

Roettgen (2015): use Rotational Raman Scattering for

calibration using the J = 6 → 8 transition of N2 at P = 100 Torr

Timothy Chen, Princeton, 2017

Decoupling Raman and Thomson Signals

A. van Gessel, E. Carbone, P. Bruggeman and J. van der Mullen, Plasma Sources Science and Technology, vol. 21, no. 1, p. 015003, 2012.
 Filtered Thomson Scattering:
ne , Te, and EEDF inference
Rayleigh
8.0E+04
7.0E+04 scattering blocked
He, 200 Torr 6.0E+04
5.0E+04 Gaussian
4.0E+04 Fit
3.0E+04
10 mm 2.0E+04
1.0E+04
0.0E+00
525 530 535 540
Wavelength (nm)
N2 Raman scattering
60000

50000

Intensity [a.u.]
40000

• Electron density: area under Thomson scattering spectrum 30000

20000
• Electron temperature: spectral linewidth
10000
• Gaussian scattering lineshape: Maxwellian EEDF 0
528 532 536
• Raman scattering rotational transitions in N2 used for absolute calibration Wavelength [nm]

Courtesy of Prof. Igor V. Adamovich

 Thomson Scattering Spectra
Ns pulse discharge in H2-He and O2-He, P=100 Torr

Thomson signal
16000 Gaussian fit 30000 Synthetic spectrum
Experiment
25000
12000
Intensity [Counts]

Intensity [a.u.]
20000
8000
15000

10000
4000
5000
0
0
524 528 532 536 540 526 528 530 532 534 536 538

Wavelength [nm] Wavelength [nm]

5% H2-He, 10% O2-He

ne = 1.5∙1014 cm-3, Te = 2.0 eV ne= 1.7·1013 cm-3, Te= 1.6 eV, T=350
 Electron Density and Electron Temperature
Ns pulse discharge in O2-He

Experimental ne
4 Predicted ne 4
10% O2-He
Experimental Te
Predicted Te
3 3

ne [1014 cm-3]

Te [eV]
2 2

1 1

0 0
0 200 400 600 800
Time [ns]

• “Double maxima” in ne, Te : two discharge pulses ≈ 400 ns apart

• Electron temperature in the afterglow Te ≈ 0.3 eV (controlled by superelastic collisions)

• Modeling predictions in good agreement with data

• Measurements in air are more challenging (strong interference from N2, O2 Raman scattering)
Measurements of electron number density

Helium

K. Sasaki et al., Contrib. Plasma Phys. 55, No. 8, 563 – 569 (2015)

2v0 2 ln(2)k BTe  

Thomson scattering profile (Broadening): v D'  sin  
c me 2 θ: the scattering angle relative to the laser beam.
 d 
Ae : integrated Thomson scattering signal intensity
  AHe,J : integrated He J level Raman transition intensity
Ae  d  532 nm
ne  n He f J f,J : relative J level population fraction in distribution function
AHe, J  d 2 
 
 d 
 Measured electron energy distribution, temperature, and number density

Fig.1Typical electron energy distribution Fig.2 Pressure dependences of (a) the electron density and (b) the electron
function measured by laser Thomson temperature. Values observed at three delay times after the initiation of the
scattering in a microwave helium plasma at microwave power are plotted.
a pressure
of 0.3 MPa.

• Thomson scattering is its weak scattering signal intensity owing to the low number density of free electrons in the plasma.
• Strong interference from Rayleigh scattering as well as plasma emission

K. Sasaki et al., Contrib. Plasma Phys. 55, No. 8, 563 – 569 (2015)
Femtosecond Localized E-Field
Measurement (FLEM)
• In a centrosymmetric medium, second • Described as a third order nonlinear
harmonic generation is impossible process:
• I(2 ω)∝ N2(EExt)2(IPump)2
• Applying an electric field destroys that • I(2 ω) : Second Harmonic Intensity
symmetry allowing for E-Field • N: Number Density
measurements • Eext: Applied Field to be Measured
• IPump: Pump Beam Intensity
• Benefits of FLEM Method:
• Signal scales as E2
• Works well at higher pressure
• Time resolution determined by pump beam duration
• Non-resonant method works in any species and gas mixtures
• Spatial resolution determined by beam focusing parameters

Courtesy of Prof. Richard Miles, Princeton Supported by the Army Research Office grant W911NF-15-1-0236 under Dr. Matthew Munson.
Femtosecond Localized E-
Field Measurement (FLEM)
• Sub-breakdown electric
field applied
• SHG response, pump
intensity, current and
voltage monitored
• Quadratic dependence
verified
• We have measured down to  Voltage rise time of 20 ns
100V/cm in room air
 Collected and analyzed ~30,000
 Fs laser pulse acts as a δ function compared with individual waveforms
ns high voltage pulses
 Determine when laser pulse arrives
 Temporal resolution determined by oscilloscope with respect to high voltage pulse
rather than a physical limit
 Bin and averaged into discreet time
values
Courtesy of Prof. Richard Miles, Princeton
 Electric Field Measurements in 2-D Ns Pulse
Discharge in Atmospheric Air
Laser beam
locations

Voltage [kV]
2.5 Current [A] 2.5 3
Coupled energy [mJ]

0.0 0.0
2

-2.5 -2.5
1

-5.0 -5.0
0

-7.5 -7.5
-100 -50 0 50 100 150 200 250
Time [ns]

• Ns pulse discharge between a high-voltage electrode and a thin quartz plate

• Discharge gap 0.6 - 1.0 mm, two-dimensional geometry, diffuse plasma
• Time-resolved electric field measured at multiple locations in the discharge gap
 “Curtain Plasma” Images, Negative Polarity Pulse

Front view, 100 ns gate

Laser beam
locations

Side view, 2 ns gate Top view, 2 ns gate

• Surface ionization wave plasma ~ 200 μm thick, wave speed ~ 0.03 mm/ns

• Electric field measured by picosecond four-wave mixing (calibration by electrostatic field)

• Time resolution 2 ns, spatial resolution across laser beam ~ 100 μm

• Objective: electric field mapping in ns pulse discharges in high-pressure fuel-air mixtures

 Electric Field Vector Components
in a Surface Ionization Wave Discharge
Forward Voltage
8 Current 40 30
breakdown - Ex
Absolute field
25

Electric field [kV/cm]

6 Actual field 30 Ey

Electric field [kV/cm]

4 20
(Ex2 + Ey2)1/2
20
-U [kV], I [A]

2 10
15
0 0
10
-2 -10
HV electrode 5
-4 -20
Reverse
breakdown Laser beam 0
-6 -30
-100 0 100 200 300 400 locations -100 -50 0 50 100 150 200
Time [ns] 150 μm from surface Time [ns]

• Initial field offset (at t < 0): charge accumulation on dielectric from previous pulse
• Field follows applied voltage rise, increases until “forward breakdown”
• After breakdown, field reduced due to charge accumulation on dielectric
• Field is reversed after applied voltage starts decreasing
• Away from HV electrode, field peaks later (Ey before Ex): surface ionization wave
• Measurements in a hydrogen-air diffusion flame underway
Plasma property measurements using H2/Ar emission lines
3
D1 ,3 D 2 ,3 D3 , 3
P0 ,3 P1 ,3 P2 , 3
S1 , 1
D 2 ,1 P1 ,1 S0 ,
3s23p54p1
L-S coupling
3s23p6↔ 3s23p54S1
15.75eV ionization
L  l1  l2  1  0  1

14.7eV J  L  S ,... L  S  2,1,0

4d Term : 3
P0 ,3 P1 3 P2

S  S1  S 2  1 / 2  1 / 2  0
13.3eV 3d J  L  S ,... L  S  1
4p Term : 1P1

4s
3s23p6↔ 3s23p54p1

L  l1  l2  1  1  0,1,2 (any quantum number)

1 1
S  s1  s2    0,1
2 2
J  L  S ,... L  S  2,1,0, 1, 2,1,0, 3,2,1
S=0 S=1
Term : 1S0 ,3 S1 , 1
P1 ,3 P0 ,3 P1 ,3 P2 , 1D 2 ,3 D1 ,3 D 2 ,3 D3
 Stark broadening of hydrogen lines and Ar optical emission line-ratio method

Schematic diagram of the experimental setup. The spatially resolved optical

measurement system is shown on the left bottom. On the right bottom is a zoom-in
figure showing the stainless steel needle tips and the discharge gap.
Xi-Ming Zhu, James L Walsh, Wen-Cong Chen1 and Yi-Kang Pu, 15
J.
Phys. D: Appl. Phys. 45 (2012) 295201 (11pp)
Experimentally measured electron densities in a high-pressure
nanosecond pulsed microplasma (Ar/Ne = 700/30 Torr, discharge Uncertainties in the ne measurement (%) using the
current lasts for about 100 ns, pulse period 1 ms). In the legend on Stark roadening method with Ar 696.5 nm line, Hα
the right top, ‘line ratio’ refers to the line-ratio method and ‘Stark line and Hβ line (for ne > 1016 cm−3) and that using
broadening’ refers to the Stark broadening method using Ar 696.5 the line-ratio method.
nm, Hα and Hβ lines with a single-Voigt fitting procedure. In the
legend on the left bottom, ‘centre’ and ‘edge’ denote ne,centre and
ne,edge obtained with double-Voigt fitting. The solid line shows a Xi-Ming Zhu, James L Walsh, Wen-Cong Chen1 and Yi-Kang Pu, J.
function, ne = 6 × 1018 × exp(−(t/0.15)0.22), where ne and tare in Phys. D: Appl. Phys. 45 (2012) 295201 (11pp)
units of cm−3 and ns, respectively.
16
 Measurements of temperature, vibrational level populations
and fast heating using picosecond CARS

Air, P=100 Torr

2 mm
10
mm

t= 1-10 μs (frames are 1 μs apart)

ns pulse discharge and afterglow: Air vs. nitrogen, P=100 Torr

• Compression waves formed by “rapid” heating, on sub-acoustic time scale, τacoustic ~ r / a ~ 2 μs

• What processes control other features of temperature rise (e.g. “slow” heating”)?
 Comparison with modeling predictions in air:
vibrational kinetics and temperature rise

• Strong vibrational excitation in the discharge, N2(v=0-8)

• Tv(N2) rise in early afterglow: V-V exchange, N2(v) + N2(v=0) → N2(v-1) + N2(v=1)

• Tv(N2) decay in late afterglow: V-T relaxation, N2(v) + O → N2(v-1) + O , radial diffusion

• “Rapid” heating: quenching of N2 electronic states, N2(C,B,A,a) + O2 → N2(X) + O + O

• “Slow” heating: V-T relaxation, N2(X,v) + O → N2(X,v-1) + O

• “Rapid” heating: pressure overshoot , compression wave formation

• NO formation: dominated by reactions of N2 electronic states, N2* + O → NO + N

A. Montello, Z. Yin, D. Burnette, I.V. Adamovich, and W.R Lempert, Journal of Physics D: Applied Physics 46 (2013) 464002
 Single-shot measurement rotational and rovibrational energy distributions
by Hybrid fs/ps coherent anti-Stokes Raman scattering (CARS) spectroscopy
fs ps

(Four wave mixing and fs broadband dual pumping)

ωp1: Rovibrational Raman transition (Q-branch, Δv=+1, ΔJ =0)
ωp2: pure rotational Raman transition (S-branch, Δv=0, ΔJ =+2)
ωprobe: frequency-narrowed ps probe pulse

The He∕N2 dielectric barrier discharge

Dedic, C.E., Meyer, T.R. and Michael, J.B., 2017. Single-shot ultrafast coherent anti-Stokes Raman scattering of vibrational/rotational
nonequilibrium. Optica, 4(5), pp.563-570.
Experimental measurements of Plasma chemistry and Kinetic Processes

NRP discharge in
4.5 mm
air at 1000 K, 1 atm:
NRP spark • Measured quantities:
• 10-ns pulse grounded
discharge
• O atoms: TALIF with absolute calibration (Xe)
• 5.7 kV electrode • N2 (A): CRDS
• N2 (B) and N2 (C): OES
• 10 kHz • Temperature: OES
• Gap: 4 mm • Electron density: Hb Stark broadening

• 670 mJ/pulse Preheated air

at 1000 K

Courtesy of Prof. Christophe Laux

 Measurements of V, I, temperature, densities
40 6 300
5 250

Voltage (V)

E/N [Td]
Current [A]
30
4 V Iconduction 200
20 3 150
10 2 100
1 50
0 0 0

Temperature [K]
2500 Temperature
Ultrafast heating: from N2(C-B)
hheating =21±5%
900 K in 20 ns 2000 from N2(B-A)
1500
18
1.2x10
18 3
Absolute densities [cm-3]

1.0x10 O ( P) density
Ultrafast 8.0x10
17

dissociation of O 17
26.0x1017 hdiss. = 35±5%
4.0x10
17
2.0x10
0
17
10
10
16
N2(B)
15
10
14
10 N2(A)
10
13
N2(C)
12
10
-10 0 10 20 30 40 50
Time (ns) Rusterholtz et al, J. Phys.D, 46, 464010, Dec 2013
 Summary of processes involved in
flame stabilization by NRP discharges
Chemical effects:
RH + O  R + OH
O2 Oxidation
2O

N2(X) N2(A) O2
e- N2(B) N2(X) + 2 O + E
N2(C)

T

Thermal effects
2-step mechanism (Popov, 2001):
N2 + e → N2* + e (N2* = N2 A, B, C, …) 5 μs after pulse
Thresholds: 6.2, 7.4, 11.0 eV
(Xu et al., APL. 99, 121502, 2011)

N2* + O2 → N2 + O + O + T
T = 1.0, 2.2, 5.9 eV
 2. Measurements of Chemical Processes in Plasma Assisted Ignition and Combustion

O atom measurements by using TALIF

O atom mole fraction

Air
5.0E-5
Air-ethylene, =0.5 Atomic O production
600
4.0E-5 525

Extinction strain rate (1/s)

3.0E-5 O (3P) 450
375
2.0E-5
300

1.0E-5 225
no plasma
150 with plasma (f=5 kHz)
with plasma (f=20 kHz)
0.0E+0
1.0E-7 1.0E-6 1.0E-5 1.0E-4 1.0E-3 1.0E-2 0.30 0.31 0.32 0.33 0.34 0.35 0.36
Time, seconds Fuel mole fraction Xf
O atom formation in a plasma discharge of O atom formation in a ns plasma Extension of extinction limit by
air and air-C2H4 mixture in a flow reactor discharge of methane/air plasma discharge
counterflow flames

J. Uddi et al. 2009 W. Sun et. Al. 2010

 OH measurements in a flow reactor:
Plasma chemical reactions result in ignition
End View
Pulse #10 Pulse #100
H2 – air, ϕ=0.3
T0=500 K, P=100 torr
C2H4 – air, ϕ=0.3
T0=500 K, P=100 torr

50 pulses

H2-air, ϕ=0.4

Short burst: OH transient rise and decay

Long burst: plasma assisted ignition, Tignition ≈ 700 K < Tauto-ignition ≈ 900 K
 Plasma Assisted Combustion: Change of ignition and extinction S-curve
The effect of kinetic enhancement (μs ~ ms, 800-1200 K)
New “S-curve” by Plasma assisted combustion for
small molecule fuel such as H2, CH4

15
7x10
O2=34%
Extinction 15
O2=62% CH4

OH number density (cm )

6x10

-3
Temperature

Plasma generated 5x10

15

species: 15
Smooth
4x10 Transition
O, H, O2(a∆g) … Extinction
15
3x10
plasma S-curve
Plasma the classical S-curve 2x10
15

Ignition 1x10
15

Ignition

Residence time 0.05 0.10 0.15 0.20 0.25 0.30 0.35

Fuel mole fraction

Scramjet, afterburner P = 72 Torr, f = 24 kHz, a = 240 1/s

•Strong kinetic enhancement at intermediate temperature
•Less effect at high temperature

Sun et al. Proc. Comb. Inst. 34, 2010, Combust. Flame 2011, 2012
Ombrello et al. 2008
 Plasma assisted low temperature combustion
Methane vs. Dimethyl ether (DME)

P = 72 Torr
f = 24 kHz

OH* emission ~310 nm

30 ms gate
Laser beam
OH, CH2O PLIF

25.4 mm

E = 7500 V/cm, E/N ~ 900 Td

Peak Voltage
= 7.8 KV Power ~ 17 W (repetitive pulses)

26
 OH PLIF measurements in Dimethyl ether (DME) Ignition

S-shaped ignition and extinction curves

DME vs. CH4
Top burner Direct image
OH
fluorescence
at Q1(6)

Flame
Bottom burner (fuel)

OH density vs. fuel mole fraction XO = 0.55, P = 72 Torr, f = 24 kHz, for DME (a = 250 1/s) and CH4 (a = 400 1/s,) as the fuel,
respectively (solid square symbols: increasing XF, open square symbols: decreasing XF)

What is the role of plasma before ignition of DME?

 Plasma activated LTC: change of S-Curve

P = 72 Torr, a= 250 1/s, f = 24 kHz P = 72 Torr, a= 250 1/s, f = 34 kHz,

XO2=40%, varying Xf XO2=60%, varying Xf

5
5 6x10 increase
6x10 LTC decrease
5
5x10
5
Extinction
5x10 LTC

CH2O PLIF (a.u.)

CH2O PLIF (a.u.)

5
HTC
4x10 5
4x10
5
3x10 5
HTC 3x10
5
2x10 Hot Ignition
5
increase 2x10
5
1x10 decrease
5
0.00 0.02 0.04 0.06 0.08 0.10 0.12 1x10
Fuel mole fraction
0.00 0.02 0.04 0.06 0.08 0.10 0.12
S-Curve Fuel mole fraction
New ignition/extinction curve without
extinction limit

Radical production by plasma can activate LTC at much shorter timescale, lower
pressure and temperature; and enable new flame regimes 28
 Flow reactor studies of plasma assisted
low temperature reaction

Comparison of thermal (1, 10 atm) and plasma (1 atm)

propane fuel consumption as a function of temperature.

Nicholas Tsolas, Jong Guen Lee and Richard A. Yetter, 2015, Phil. Trans. R. Soc. A 373: 20140344.
Plasma activated Cool Flames :A new way to burn with plasma
1 n-heptane


Ignition delay time (s)

1 atm
5 atm
10 atm
20 atm
0.1
Plasma

0.01
Extinction
Temperature

0.8 1.0 1.2 1.4 1.6

Plasma 1000K/T
HTC generated

LTC
LTC Ignition

t2 t1 Residence time
t2<< t1

Plasma activated LTC at much shorter

time, lower pressure….

We can create cool flames even at

1 atm or below? 30
2. Plasma activated self-sustaining Cool Diffusion Flames:
A new way to burn

Fuel/N2 @ 550 K
Tf~1900 K
Heated N2 @ 550 K

Stagnation
plane (a) Hot diffusion flame
N2 @ 300 K

Oxidizer @ 300 K
with plasma discharge Tf~650 K
Fig. 1 Schematic of experimental setup
(b) Cool diffusion flame

Fig. 2 Hot and cool n-heptane

diffusion flames at the same condition

Won et al., 2014

 3. Multispecies diagnostics and kinetic modeling
Fuel(RH)
+OH
+O2 Small C2H3/CH2O
RO2 R
e +O2=O+O(1D) +e alkene
+O2
H+O2(1Δg) =O+OH QOOH HO2 H/HCO
O(1D)+RH =OH+R +O2+(M)
N2(A,B,C)+O2=O+O+N2 O2QOOH H2O2
+O2 CO/CO2
N2(v)+HO2 =OH+O+N2
R(v,*)+O2=RO+OH 2OH
Plasma
=??? Slow e, R*, N2*, O2*
R(*), R(v), N2(v), O2(v)
O3+M =O+O2+M

A schematic of the key reaction pathways for high pressure fuel

oxidation of at different temperatures
(blue arrow: Below 700K; yellow arrow: 700-1050 K;
red arrow: above 1050K).
 Question to aks:
How does the key plasma reactions affects n-heptane dissociation and oxidation
in the first 10 ms with efifferent excitation processes involving Ar and O2?

e +Ar =Ar*+e
e +O2 =O+O(1D)
e +O2 =e + O2(v)
e +O2 =O2-
e +C7H16 =H+R
......
Ar* +C7H16 =?
O(1D) +C7H16 =?
O2(v) +C7H16 =?
O2(v) +C7H15 =RO2(v)
......
 Plasma chemistry reactor
• Reactor
• Pressure: 60 Torr
• Initial Temperature: 300 K
• Flow speed: 40 cm/s Mini-Herriott cell
• Reactor size: 45 x 14 x 152 mm3 showing 24 pass
configuration
• Nanosecond repetitively
pulsed discharge: FID GmbH Vacuum
Chamber
FPG 30-50MC4
• Peak Voltage: ~7 kV
• Pulse Duration: 12 ns
FWHM
Reactor
• Continuous mode: 0 – 5 kHz
• Pulse burst mode: 150
pulses, 30 kHz
Ge Etalon
• Quartz double dielectric Macor
Quartz
barrier: 1.6 mm thickness Wall
Wall
Flip Mirror
• Gap distance: 14 mm
Collimating Mirror
Lenses

34
 Experimental Apparatus
Laser inlet
purge tube N2 Purge
Box
Electrode
Connection QCL Laser

Vacuum
Chamber

Observation
Window

Alignment
Detector Laser

To Vacuum
Direct and ICCD Images of Plasma Discharge in a Reactor

Stoichiometric mixtures: C2H4/O2 with 75% AR, 60 Torr, Vmax= 6 kV

•Direct Image: 1 kHz, 3.6 mJ/pulse, 2 s exposure time.

•ICCD images: Gate time = 100 ns, Gain = 250

Cathode
Direct 1000 Hz
Anode
ICCD
1000 Hz

2000 Hz

3000 Hz
 Experiment/Model Comparison (C2H4)
Oxidation Pyrolysis
500 600
Oxidation C2H2 Pyrolysis C2H2
Oxidation CH4 Pyrolysis CH4
Oxidation H2O C2H2 Model
C2H2 Model CH4 Model
CH4 Model 500
400 H2O Model
Mole Fraction (ppm)

Mole Fraction (ppm)

400
300

300

200
200

100
100

0 0
0 0.005 0.01 0.015 0 0.005 0.01 0.015
Time from first pulse (s) Time from first pulse (s)

• Conditions for Oxidiation case: • Conditions for pyrolysis case:

– Stoichiometric C2H4/O2 – 0.9375 Ar/0.0625 C2H4
– 25% Reactants, 75% Argon
– V=40 cm/s
– V=40 cm/s
– P=60 Torr – P=60 Torr
– Ti=300 K – Ti=300 K
– 150 Pulses, 30 kHz, 10 kV – 150 Pulses, 30 kHz, 10 kV
– E/n = 37 Td – E/n = 37 Td
• Acetylene measurements in the pyrolysis • Methane is greatly over-predicted by the
experiment were used to match the E/n ratio model
of model calculations
 Continuous Plasma - Oxidation

• GC and in situ IR-LAS

measurements 1000000
– Continuous plasma
– Stoichiometric, 25% 100000

Mole Fraction (ppm)

Reactant
• Excellent agreement 10000

between GC and IR-

LAS 1000

• Major Species 100

– CH2O
O2 C2H4
– CO CO CH2O
10
– C2H2 C2H2 C2H2 LAS
CH4 CH4 LAS
– CH3OH, C2H4O CH3OH/C2H4O
1
– CH4 0 1000 2000 3000 4000 5000
Plasma Frequency (Hz)

Plasma activated low temperature fuel oxidation is an important process

38
Lefkowitz et al. 2014
 Validation of plasma combustion chemistry of C2H4 oxidation
C2H4/O2/Ar: 6.25/18.75/75

HP-Mech (a) USC-Mech (b)

Fig. 16 Measurements and predictions of C2H2, CH4, and H2O concentrations after 150 pulses at 30
kHz repetition rate for a mixture of C2H4/O2/Ar: 6.25/18.75/75 by using (a) HP-Mech and (b) USC-
Mech II.
 Ethylene Oxidation Pathways
+ Ar(+) 13%
+ OH 15% C2H4 C2H3+

CH2CH2OH + Ar* 5%
+ e- 30% + e- 65%
+ O2 100% LTC
O2C2H4OH
HTC C2H2 C2H
100%
2 CH2O + OH + O 11% + O 13% + O2 100%

+ H +M 31%
+ H 21% HCO + CO
H + CH2CHO CH3+ HCO
C2H5
+ O 21% + O2 46% + O2 + M 85%

+ O2 + M 97% CH20 + HCO CO + CH2O + OH CH3O2

+ X 95%
C2H5O2 M = Third body collider Blue = Plasma CH3O
X = Radical Red = High temperature, + X 96%
+ HO2 98% Green= Low temperature

C2H5O2H CH2O
PAC activates C2H4 low temperature chemistry
Challenges: Plasma activated CH4/O2/He oxidation at 400 K

250

CH2O Experiment
CH2O Model
CH4
Mole Fraction (ppm) 200

+ O(1D) 2%
150

100

CH3 + OH CH2OH + H
50

+ O2 100%
0
0
Pulse
5
Burst10 15 20
Time (ms)

CH2O + H2/HO2
 300 pulse burst, 8.75 kV peaks voltage, 30
kHz pulse repetition frequency
 Stoichiometric mixture, 75% diluent, 60
Torr, 300 K initial temperature
Lefkowitz, J.K., Guo, P., Rousso, A. and Ju, Y., 2015. Phil. Trans. R. Soc. A, 373(2048), p.20140333.
 Challenges: Time dependent Formaldehyde
Measurement
• CH2O production slightly
greater for stoichiometric 70

case
60
– Similar linear trend as in

Mole Fraction (ppm)

n-heptane consumption 50
• Lack of further CH2O production
40
after pulse burst: indicates that
only plasma produced species are
responsible for fuel oxidation 30

• Mode significantly under- 20

predicts measurement by ϕ=1
10 ϕ = 1, Modeling
a factor of 40. ϕ = 0.5
• Missing pathways for low 0 ϕ = 0.5, Modeling
0 0.002 0.004 0.006 0.008 0.01
temperature CH2O
Time (s)
formation

Rousso, Aric, Suo Yang, Joseph Lefkowitz, Wenting Sun, and Yiguang Ju.
" Proceedings of the Combustion Institute, Vo.36, 2017, Pages 4105–4112
42
Cross-sections database available for electron-molecule collisions

Andrey Starikovskiy and Nikolay Aleksandrov, AIAA paper-2017-1977

 Non-equilibrium plasma kinetics
A new reaction pathway of plasma assisted low temperature combustion
via excited RO2*(v)
Vibrational and electronically excited O2(v) collides with fuel radical (R) forming highly
energized RO2* in which the vibrational energies are quickly redistributed due to strong
coupling between different vibration states. These RO2*, comparing to those formed by
ground state oxygen with R, carry much higher internal energies that enable them
overcome the barrier TS1, TS2, and TS3 much easier to produce the bimolecular product
HO2+Alkene and OH+Ether. Therefore, the ignition processes/species profiles of the
mixture can be significantly different from the ground state system.

• Increase the rate

R+O2*(v,e) • Change branching ratios (TSi)
R(v)*+O2 Energized • Modify pressure dependence
RO2*(v)
TS3
R+O2
TS2
TS1 CH2O + OH + R”

HO2+alkene
QOOH
RO2 OH + RO (cyc-Ether)
Rousso, Aric, Suo Yang, Joseph Lefkowitz, Wenting Sun, and Yiguang Ju.
" Proceedings of the Combustion Institute, Vo.36, 2017, Pages 4105–4112
 O(1D) reaction kinetics: Photolysis reactor for elementary kinetic study
O(1D)/O3/O2/CH3OH/Ar mixture
time dependent measurements of OH, HO2, CH2O,…

H2O:1338.55 cm-1
CH2O:1726.8 cm-1 913 mm long, 40 mm diameter, multi-pass (21)
OH: 3568 cm-1 Herriott cell, a19.17 meter optical path length
HO2: 1397 cm-1
Ju, Y., Lefkowitz, J.K., Reuter, C.B., Won, S.H., Yang, X., Yang, S., Sun, W., Jiang, Z. and Chen, Q., 2016. Plasma Chemistry and Plasma
Processing, 36(1), pp.85-105.
 Experimental results and model validation: updated model
O(1D)+CH3OH→CH2OH+OH knew=1.5x1014 (R27)
O(1D)+CH3OH→HOCHOH+H knew=0.5x1014 (R28) (in mole-cm3-s)
O(1D)+CH3OH→CH2O+H2O knew=1.0x1014 (R29)

Time-resolved mole fraction of H2O in the 266 nm laser photolysis of 0.224% CH3OH 1.91%
O2 and 596 ppm O3 in Ar mixture with the variation of CH3OH flow rate compared to model
simulations for 1.0 ml/hour flow of CH3OH. ○: Experimental measurement; ― : simulation
using the original model. ― : simulation with updated reaction rates [38].
Ju, Y., Lefkowitz, J.K., Reuter, C.B., Won, S.H., Yang, X., Yang, S., Sun, W., Jiang, Z. and Chen, Q., 2016. Plasma Chemistry and Plasma
Processing, 36(1), pp.85-105.
 Potential energy surface (PES) of O(1D)+CH3OH
using M062X/cc-pvtz level

New channel?

Ju, Y., Lefkowitz, J.K., Reuter, C.B., Won, S.H., Yang, X., Yang, S., Sun, W., Jiang, Z. and Chen, Q., 2016. Plasma Chemistry and Plasma
Processing, 36(1), pp.85-105.
HO2/OH using mid-IR Faraday Rotational Spectroscopy

VRMS ( )   GP0 sin 2   RMS

+Bfield

Laser Lock-In Amplifier

Paramagnetic (radical) species

HO2 energy levels

Zeeman splitting
Absorption

Dispersion


 Experimental results: HO2/OH measurements

Signal
HO2
Sensitivity OH

 detection limit  1 ppmv / 𝐻𝑧 3 detection limit = 20 ppbv / 𝐻𝑧

DME flow reactor

model validation

Implication:
RO2→QOOH→O2QOOH uncertainty
HCO+O2=HO2+CO reaction uncertainty and HCO formation pathway?
Bremfield et al., 2013, JPC letters, 2013; Kurimoto et al. 2014
 Summary
• Time-resolved, spatially-resolved, in-situ laser diagnostics of electric field, electron
density, and electron temperature, excited and radical species greatly enhanced the
understanding of plasma kinetic and chemical process in PAC.

• Production of O and O(1D), O2(singlet), N2(*), and N2(v) by the plasma is the major
processes in kinetically enhancement of combustion.

• Fast and slow heating in PAC is important, but the energy transfer processes are very
complicated.

• Plasma activated low temperature combustion pathways and enable cool flame
formation, but existing mechanisms have large uncertainties, especially for large
hydrocarbons.

• Electron impact reaction cross-sections of large alkanes and reaction rates involving
O(1D) and non-equilibrium excitations are poorly known.

• The effect of vibrational species excitation on PAC is still poorly known.

Lecture 5 Chemistry and Kinetic Studies of Plasma-Assisted Combustion

Yiguang Ju

• Important chain-initiation and branching reactions in combustion

• Plasma chemistry and timescales
• Impact of plasma chemistry on combustion
• Diagnostics of plasma properties and chemistry in PAC
 1. Important combustion reactions
1. Chain initiation and propagation reactions
RH+ O2 → R+HO2 High Temperature (>1100 K) slow
RH+HO2 → R+H2O2 High pressure/low temperature (>550 K) slow
R+HO2 → RO+OH High pressure/low temperature (>550 K) slow

2. important branching reactions at different temperatures

H+ O2 → O+OH High Temperature (>1100 K) Fast
H2O2 → 2OH Intermediate temperature (800-1100 K) Slow
R → RO2→QOOH → O2QOOH →R’’+2OH Low temperature (300-800K) Slower

Plasma assisted combustion:

O+RH → R+OH
e+O2 → e+ 2O
R → R’’+2OH
O+RH → R’’+ 3OH Faster

e+O2 → e+ O2(a1Δg) H+O2(a1Δg) → OH+O Faster

Plasma provides new reaction pathways to accelerate chain reaction processes

 Interaction of plasma chemistry with reaction kinetics of
large alkanes w/wo in plasma assisted combustion
Fuel(RH)
+OH Reaction rate
RO2
+O2 Small C2H3/CH2O Transition state theory
R alkene
+O2 k1 k2
RO2* QOOH HO2 H/HCO A  B  AB *  C  D
+O2+(M)
R*, O(1D)

O2QOOH H2O2 CO/CO2 k BT q AB * * E *  E A B

+O2 k (T )  exp( )
h q A * qB * k BT
2OH
Plasma
How does plasma affect elementary rate
I e, R*, N2*, O2*,O*
II O2(v), O2(a1Δg) constant?
RH(v), R(v), N2(v), O2(v), HO2(v)
e.g. at 800 K
A schematic of the key reaction pathways for high
pressure fuel oxidation of at different temperatures O2 (a1Δg) + H = OH+O Fast
(blue arrow: Below 700K; yellow arrow: 700-1050 O2 + H = OH+O Slow
K; red: above 1050K). Green: plasma activated
pathway CH3 +O2(v) → CH2O+OH Fast
CH3 +O2 → CH2O+OH Slow
 2. Plasma chemistry and timescales of
kinetic processes

ttr telec
Ion/Molecule
Kinetics
tfp trot
Ion-Ion, Ion-Molecular

EEDF Excitation /
Quenching
Electron Kinetics
tfp Ionization Recombination

Combustion
ttr tvib Molecules

tfp
Fig. 1.5 Schematic of timescales and key kinetic
Processes
trot Radicals pathways at different stages of plasma assisted
ignition and combustion.

10-14 10-12 10-10 10-8 10-6 10-4 10-2 s Ju and Sun, PECS, 2015

Courtesy of Andrey Starikovskiy

 Potential Energy Curves of O2
O2(B3Su-), 8.4 eV
smax = 1.0 A2 (9.4 eV)
DE ~ 1 eV
O2(3Pg), 5.6 eV
smax = 0.16 A2 (12 eV)
DE ~ 1.5 eV
O2(A3Su+), 4.5 eV
smax = 0.18 A2 (6.6 eV)
E, eV

O2 (a1Δg) O2 (b1Σg+) at 1.6 eV

O2 (a1Δg) at 0.98 eV

r, nm
Electron impact reaction is a function of electron energy distribution (E/N)
 Electron impact reaction cross sections-O2

1. effective 7. O2(a1) 8. O2(b1) 9. O2(A3Su+), 4.5 eV

2. rotational excitation 10. O+O 11. O+O(1D) 12. O+O(1S)
3-6. O2(v1) - O2 (v4) 13. O2+
 Potential Energy Curves of N2

N2(A3Su+), 6.2 eV
smax = 0.08 A2 (10 eV)

N2(B3Pg), 7.35 eV
smax = 0.20 A2 (12 eV)
E, eV

N2(C3Pu), 11.03 eV
smax = 0.98 A2 (14 eV)

Threshold energy diagram

r, nm Electron impact reaction is a function of electron energy distribution
Energy Transfer of non-equilibrium excitation in
Plasma Discharge
N2:O2:H2 = 4:1:2
1
N2(el)
N2(v)

H2(v) H2(el) ion

Energy loss fraction
0,1 H2(rot)
Rot+tr
O2(dis)
O2(v)

0,01
O2(4.5 eV)
O2(a+b)

O2 (a1Δg)
1E-3
1 10 100 1000
E/N, Td

Physics of Nonequilibrium Systems Laboratory

 Influence of Electronic and Vibrational Excitation on
Combustion Kinetics

N2:O2:H2 = 4:1:2
1
N2(v)
N2(el) N2 + e = N2(C3) + e
H2(v) H2(el) ion N2(C3) + O2 = N2 + O + O
Energy loss fraction

0,1
Rot+tr
H2(rot)
O2 + e = O + O + e
O2(dis)
O2(v)

0,01
O2(a+b)
O2(4.5 eV) N2 + e = N2(v) + e
N2(v) + HO2 = N2 + HO2(v)
1E-3
1 10 100 1000 HO2(v) = O2 + H
E/N, Td

Physics of Nonequilibrium Systems Laboratory

 Influence of Vibrational Excitation on Low-
Temperature Kinetics: H2O2 Decomposition

Measured and calculated OH decay time. P = 1 atm.

a) 3%H2 + air; b) 0.3%C4H10 + air.
Physics of Nonequilibrium Systems Laboratory
 PRINCETON
University
Effect of “Hot” Atoms on Active Species Production
in High-Voltage Pulsed Discharges
Nonequilibrium distributions of neutral species are formed in different physical situations.
In laboratory experiments and in the terrestrial atmosphere, there are numerous collisional
processes in which translationally energetic (superthermal) atoms with energies much
above thermal energies are produced.

Potential energy curves Momentum transfer cross Cross sections for scattering of H
and hot atoms formation section for the H-H2 scattering atoms with H2, O2, CH4 and N2
[15]
[33]
[23] (a)
This work

2
cm

2
cm
-16
Cross section, 10
1
10

-16
Cross section, 10
10

H-H2 (el)
H-O2 (el)
H-CH4 (el)
H-N2 (el)
0
10 H+O2=O+OH (new)
H+CH4=CH3+H2 (new-1)
H+CH4=CH3+H2 (new-2)

1
0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5 6,0 6,5 7,0
0,1 1
Energy, eV
E, eV

Direct electron-impact dissociation

e + O2 → e + O2* O2* → 2O(3P,1D) + 1.3 eV
e + H2 → e + H 2* H2* → 2H(1S) + 4.5 eV
e + CH4 → e + CH4* CH4* → CH3 + H + 3.5 eV
 PRINCETON
University
Effect of “Hot” Atoms on Active Species Production
in High-Voltage Pulsed Discharges
Average energy of H atoms Species produced during energy Plasma-assisted oxidation
in various gaseous mixtures degradation of one “hot” H atom in H2-O2 mixture
H CH3 1.8 2H2-O2 mixture
(a) 1

Species amount per one hot atom

T=300 K; P=1 atm
1.6
H2 [H2O]
OH

Production, ppm/10 eV/molec

0 1.4 [H2O2]
10

CH4:O2:N2=1:2:8 -1
10 O 1.2
Energy, eV

-5
1.0
CH4:O2=1:2
0.8
HO2
-2 H2O 0.6
10
-1
H2:O2=1:9 10
H2:O2=2:1 0.4

0.2

-3
10 0.0
0 10 20 30 40 50 0.4 0.5 0.6 0.7 0.80.9 1 2 3 4 5 6 7 8

s
)
D

(h

al
at
Number of collisions Initial H atom energy, eV

(1

(1
H

ic
St
O

ad
d
)+

R
un

o
(h

ro

N
H

G
Plasma-assisted oxidation
 Analysis of the effect of formation of "hot" atoms with excessive
in CH4-O2 mixture translational energy shows the important role of these processes
0,9

0,8
CH4-2O2 mixture
T=300 K; P=1 atm
in formation of active radicals.
[H2]
 The density of radicals produced in discharge plasma can be
Production, ppm/10 eV/molec

0,7 [H2O]
[H2O2]
0,6
[CH2O] several times higher than that produced in the absence of high-
[CH3OH]
-5

0,5
[CH3O2H] energy atoms.
0,4
 The effect plays a fundamental role in the formation of the initial
0,3

0,2
distribution of active species in combustible mixtures and can
0,1 greatly influence the kinetics of ignition and oxidation at low gas
0,0 temperatures.
)

s
)
D

(h

al
at
(1

(1
H

ic
St
O

ad
d
+

R
un
)

o
(h

ro

N
H

G
Gas heating at high E/N E/N = 103 Td

50

Fast Gas Heating

Fractional power, %
40
Electron-ion recombination e + O2+ → O + O* + ΔE
15 -3
Ion-ion recombination O2 + O2 + M→ 2O2 + M + ΔE
- + ne0=10
14
cm 1 atm
-3
ne0=10 cm 1 atm
Collisional energy transfer 30 15
ne0=10
-3
cm 300 Tor
Electronically-excited species
N2(A,B,C,a) + O2  2O + DE 14
ne0=10
-3
cm 300 Tor

Hot atom and molecule O2* → 2O(3P,1D) + DE 20

0,0 0,1 0,2 0,3 0,4 0,5
Mole fraction of O2

Slow Gas Heating Fractional Electron Power Transferred Into

Heat in N2:O2 Mixtures
Vibrational energy relaxation N2(v) + M = N2(v-1) + M +DE
High oxygen, faster gas heating!
 Princeton Plasma Combustion Kinetics
Major Pathways
Ar O2 N2 H2 CxHyOz

Ionization by electron impact. k = f(E/N)

Ar+ O2+ N2 + H2+ CxHyOz+,…, CxH1Oz+

Ar, N2, O2 O2, CxHyOz N2, O2, CxHyOz O2, CxHyOz CxHyOz
H2, CxHyOz H- transfer
Charge transfer, negative and complex ions formation

Ar2+, N4+, O4+, N2O2+, NH2+, H3+, HO2+, H3O+; O-, O2-, O3-, O4-; CxHyOz+,…, CxH1Oz+
Electron-ion recombination Ion-ion recombination Molecule-ion reactions
O2+, O4+, CxHyOz+ O2- + N2+; O2- + CxHyOz+ O2- + H; O- + H2
electron detachment
Electronically- “Hot” atoms and
excited molecules formation Ionic chains
particles formation
Oh, Hh, Nh, O2h, H2h
Low-Temperature
Fast Gas Heating O(1D), O(1S), N(2D), H(n=2) Reactions
Andrey Starikovskiy
 Princeton Plasma Combustion Kinetics
Major Pathways
Ar O2 N2 H2 CxHyOz

Vibrational levels excitation by electron impact. k = f(E/N)

N2(vib) H2(vib)

VT relaxation Reactions of vibrationaly excited molecules

N2(v) + O; N2(v) + H2 H2(v) + O → H + OH(v)

N2(v) + H2O; N2(v) + CxHy H2(v) + OH → H2O + H

Slow Gas Heating Formation of vibrationaly-excited products

Energy transfer to reagents Energy transfer to buffer

N2(v) + HO2 → N2 + HO2(v) OH(v) + N2 → OH + N2(v)

Reactions of vibrationaly excited molecules

HO2(v) → H + O2 OH(v) + H2 → H2O + H

Typical plasma reactions for radical production and heating
Electron impact ionization/dissociation/excitation
>10 eV e +O2 =O++O+2e (R1a) Radical production
~10 eV e +O2 =O+O(1D) (R1b) Non-equilibrium excitation
~1 eV e +O2 =O2(1Δg)+e (R1c)
0.2-2 eV e +O2 =O2(v)+e (R1d)

Electron ion recombination, attachment, charge transfer

e+O2+ =O+O(1D) (R2a) Recombination/fast heating
O2+ +O2- =2O2 (R2b) Recombination/fast heating
e+O2 +M = O2- +M (R2c) Attachment
H2O+N2+ =H2O ++N2 (R2d) Charge Transfer

Dissociation and energy transfer by ions and excited species

N2(A,B,C)+O2 =O+O(1D)+N2 (R3a)
O(1D)+H2 = OH+H (R3b)
H+ O2(1Δg)= O+OH (R3c)
N++O2= O++NO (R3d)
CH3+HO2(v)=CH2O+OH (R3e) Slow heating
N2(v=5) +N2 = N2(v=3) + N2 (R3f)
N2(v) + HO2 → N2 + HO2(v) (R3g)
What are the major species produced by plasma?
Time

Pressure

•Long lifetime species? NO, O3, O2(a1Δg)

•Short lifetime plasma generated species? O, N2 (A,B,C)*
17
 3. Impact of plasma chemistry on combustion
Question: When will electron impact dissociation process become important in combustion?

Fig. 3.5: Rate constants (a) and reaction flux (b) for reactions for dissociation by electron
impact at electric field values equal to 200 Td and 500 Td and chain branching reactions.

Ju and Sun, PECS, 2015.

Comparison of the reaction rates of electron impact and excited species for radical production

(Ground) (Ground)

Important radical
production channels

(Ground)

S M Starikovskaia, J. Phys. D: Appl. Phys. 47 (2014) 353001

A. M. Starik, B. I. Loukhovitski, A. S. Sharipov and N. S. Titova, 2016, Phil. Trans. R. Soc. A 373: 20140341
 PAC: how does plasma change the branching reactions in combustion?

1500
Low temperature ignition Hot ignition
Kinetic effect Thermal effect

1200

Temperature (K)
>1100 K
High Temp.
H+O2=O+OH
O+H2=H+OH
900 850-1100 K
H2O2=2OH
2HO2=H2O2+O2 Intermediate
HCO+O2=CO+HO2 Temp.
CH2O+X=HCO+XH
600 R+O2=RO2 500-850 K
RO2→QOOH →R’+OH Low Temp.
O2QOOH →R’’+2OH
t1 t2
Fuel fragments Small molecules
300 Large molecules
0.0 0.1 0.2
Time (sec)

Schematic of kinetic and thermal enhancement pathways of plasma assisted combustion for liquid fuels at high, intermediate, and low temperature, respectively

Y. Ju and W. Sun, Prog. Energy Combust. Sci., 2015

Plasma activated Cool Flames :A new way to burn with plasma
1 n-heptane


Ignition delay time (s)

1 atm
5 atm
10 atm
20 atm
0.1
Plasma

0.01
Extinction
Temperature

0.8 1.0 1.2 1.4 1.6

Plasma 1000K/T
HTC generated

LTC
LTC Ignition

t2 t1 Residence time
t2<< t1

Plasma activated LTC at much shorter

time, lower pressure….

We can create cool flames even at

1 atm or below? 21
 Observation of plasma activated self-sustaining Cool Flames
Fuel/N2 @ 550 K Tf~1900 K Tf~650 K
Heated N2 @ 550 K

Stagnation
plane (a) Hot diffusion flame (b) Cool diffusion flame
N2 @ 300 K Fig. 2 Hot and cool n-heptane
diffusion flames at the same condition
Oxidizer @ 300 K
with plasma discharge 2400
nC7H16/N2 vs O2 or O2/O3
HF branch

Maximum temperature Tmax [K]

in counterflow burner
Fig. 1 Schematic of experimental setup Xf = 0.05,Tf = 550 K, and To = 300 K
2000
Extinction limit of
conventional hot diffusion flame
(HFE)
1600
without O3

Won, S.H., Jiang, B., Diévart, P., Sohn, C.H.

Transition to
and Ju, Y., 2015. Proceedings of the
1200
hot flame Extinction limit of
Combustion Institute, 35(1), pp.881-888. with O3 cool diffusion flame
HTI (CFE)
800
CF branch

LTI

400
Extinction/instability
0.1 1 10 100 1000 10000
Strain rate a [s-1]
 Plasma assisted Self-Sustaining
Premixed/partially premixed Cool Flames

• Cool diffusion flames

– n-Heptane/O2/O3
– Won et al., Proc. Combust. Inst. 2015

• Cool premixed flames

– DME/O2/O3
– Reuter et al., Combust. Flame 2016

• Cool partially premixed flames

– DME/O2/O3
– Reuter et al., Proc. Combust. Inst. 2017
(hopefully)
S. H. Won et al., Proc. Combust. Inst. 35 (2015) 881-888
23 C. B. Reuter et al., Combust. Flame (2016), in press
 4. Diagnostics of plasma physics and chemistry in PAC

1. Measurements of plasma properties and kinetic processes

2. Plasma assisted ignition and combustion with active species production

3. Kinetic studies of plasma assisted combustion

1. Measurements of Plasma Properties: electron density and temperature Thomson scattering
𝑦 Power of scattering: 𝑃𝑠 ∝ 1 − sin2 𝜃 cos 2 𝜙0
𝐸𝑖0 1
𝐸𝐿 𝑑𝜎𝑒 2𝜆0 𝜃 2𝑘𝐵 𝑇𝑒 2
Number of photo-electrons: 𝑁𝑠 = Δ𝐿 𝑛𝑒 𝜂 Δ𝜆1 = sin
ℎ𝜈0 𝑑Ω 𝑒
𝑐 2 𝑚𝑒
𝐴𝑒 𝑑𝜎𝑁2 2
𝐴𝑁2 𝑑Ω 𝑐 2 𝑚𝑒 Δ𝜆1 𝑒
𝑛 = 532𝑛𝑚
𝑛𝑁2 𝑓𝐽=6 𝑇𝑒 =
𝜙0 𝑘0 𝑒 𝑑𝜎𝑒 8𝑘𝐵 sin2 𝜃 2 𝜆0
𝜃 𝑑Ω
𝑘𝑠
𝑧 k0: Laser beam direction, ks: Scattering signal wave vector
𝑥 𝐸𝑖0 : Polarized electric field, scattering is rotationally symmetric about 𝐸𝑖0 .
x-z plane: the plane of observation
θ: the scattering angle relative to the laser beam.
𝜙0 : angle between observation plane and the polarization angle.
𝐴𝑒 and 𝐴𝑁2 : integrated intensities of the Thomson and Raman spectra
𝑑𝜎𝑒 𝑑𝜎𝑁
and 𝑑Ω2 the Thomson and N2 Raman scattering cross sections
𝑑Ω
𝑓𝐽=6 : the fraction of N2 molecules in the J = 6 rotational state
Δ𝜆1 𝑒 : the half 1/e width of the Gaussian broadening profile
EL: laser energy, η: optical efficiency, ΔL: length of observed scattering segment.
• H. Van der Meiden, "Thomson scattering on low and high temperature plasmas", Ph.D, Technische Universiteit Eindhoven, 2011.
• A. Roettgen, "Vibrational Energy Distribution, Electron Density and Electron Temperature Behavior in Nanosecond Pulse Discharge Plasmas by Raman and
Thomson Scattering", Ph.D, The Ohio State University, 2015.
Thomson Scattering Experimental Setup and Calibration

Roettgen (2015): use Rotational Raman Scattering for

calibration using the J = 6 → 8 transition of N2 at P = 100 Torr

Timothy Chen, Princeton, 2017

Decoupling Raman and Thomson Signals

A. van Gessel, E. Carbone, P. Bruggeman and J. van der Mullen, Plasma Sources Science and Technology, vol. 21, no. 1, p. 015003, 2012.
 Filtered Thomson Scattering:
ne , Te, and EEDF inference
Rayleigh
8.0E+04
7.0E+04 scattering blocked
He, 200 Torr 6.0E+04
5.0E+04 Gaussian
4.0E+04 Fit
3.0E+04
10 mm 2.0E+04
1.0E+04
0.0E+00
525 530 535 540
Wavelength (nm)
N2 Raman scattering
60000

50000

Intensity [a.u.]
40000

• Electron density: area under Thomson scattering spectrum 30000

20000
• Electron temperature: spectral linewidth
10000
• Gaussian scattering lineshape: Maxwellian EEDF 0
528 532 536
• Raman scattering rotational transitions in N2 used for absolute calibration Wavelength [nm]

Courtesy of Prof. Igor V. Adamovich

 Thomson Scattering Spectra
Ns pulse discharge in H2-He and O2-He, P=100 Torr

Thomson signal
16000 Gaussian fit 30000 Synthetic spectrum
Experiment
25000
12000
Intensity [Counts]

Intensity [a.u.]
20000
8000
15000

10000
4000
5000
0
0
524 528 532 536 540 526 528 530 532 534 536 538

Wavelength [nm] Wavelength [nm]

5% H2-He, 10% O2-He

ne = 1.5∙1014 cm-3, Te = 2.0 eV ne= 1.7·1013 cm-3, Te= 1.6 eV, T=350
 Electron Density and Electron Temperature
Ns pulse discharge in O2-He

Experimental ne
4 Predicted ne 4
10% O2-He
Experimental Te
Predicted Te
3 3

ne [1014 cm-3]

Te [eV]
2 2

1 1

0 0
0 200 400 600 800
Time [ns]

• “Double maxima” in ne, Te : two discharge pulses ≈ 400 ns apart

• Electron temperature in the afterglow Te ≈ 0.3 eV (controlled by superelastic collisions)

• Modeling predictions in good agreement with data

• Measurements in air are more challenging (strong interference from N2, O2 Raman scattering)
Measurements of electron number density

Helium

K. Sasaki et al., Contrib. Plasma Phys. 55, No. 8, 563 – 569 (2015)

2v0 2 ln(2)k BTe  

Thomson scattering profile (Broadening): v D'  sin  
c me 2 θ: the scattering angle relative to the laser beam.
 d 
Ae : integrated Thomson scattering signal intensity
  AHe,J : integrated He J level Raman transition intensity
Ae  d  532 nm
ne  n He f J f,J : relative J level population fraction in distribution function
AHe, J  d 2 
 
 d 
 Measured electron energy distribution, temperature, and number density

Fig.1Typical electron energy distribution Fig.2 Pressure dependences of (a) the electron density and (b) the electron
function measured by laser Thomson temperature. Values observed at three delay times after the initiation of the
scattering in a microwave helium plasma at microwave power are plotted.
a pressure
of 0.3 MPa.

• Thomson scattering is its weak scattering signal intensity owing to the low number density of free electrons in the plasma.
• Strong interference from Rayleigh scattering as well as plasma emission

K. Sasaki et al., Contrib. Plasma Phys. 55, No. 8, 563 – 569 (2015)
Femtosecond Localized E-Field
Measurement (FLEM)
• In a centrosymmetric medium, second • Described as a third order nonlinear
harmonic generation is impossible process:
• I(2 ω)∝ N2(EExt)2(IPump)2
• Applying an electric field destroys that • I(2 ω) : Second Harmonic Intensity
symmetry allowing for E-Field • N: Number Density
measurements • Eext: Applied Field to be Measured
• IPump: Pump Beam Intensity
• Benefits of FLEM Method:
• Signal scales as E2
• Works well at higher pressure
• Time resolution determined by pump beam duration
• Non-resonant method works in any species and gas mixtures
• Spatial resolution determined by beam focusing parameters

Courtesy of Prof. Richard Miles, Princeton Supported by the Army Research Office grant W911NF-15-1-0236 under Dr. Matthew Munson.
Femtosecond Localized E-
Field Measurement (FLEM)
• Sub-breakdown electric
field applied
• SHG response, pump
intensity, current and
voltage monitored
• Quadratic dependence
verified
• We have measured down to  Voltage rise time of 20 ns
100V/cm in room air
 Collected and analyzed ~30,000
 Fs laser pulse acts as a δ function compared with individual waveforms
ns high voltage pulses
 Determine when laser pulse arrives
 Temporal resolution determined by oscilloscope with respect to high voltage pulse
rather than a physical limit
 Bin and averaged into discreet time
values
Courtesy of Prof. Richard Miles, Princeton
 Electric Field Measurements in 2-D Ns Pulse
Discharge in Atmospheric Air
Laser beam
locations

Voltage [kV]
2.5 Current [A] 2.5 3
Coupled energy [mJ]

0.0 0.0
2

-2.5 -2.5
1

-5.0 -5.0
0

-7.5 -7.5
-100 -50 0 50 100 150 200 250
Time [ns]

• Ns pulse discharge between a high-voltage electrode and a thin quartz plate

• Discharge gap 0.6 - 1.0 mm, two-dimensional geometry, diffuse plasma
• Time-resolved electric field measured at multiple locations in the discharge gap
 “Curtain Plasma” Images, Negative Polarity Pulse

Front view, 100 ns gate

Laser beam
locations

Side view, 2 ns gate Top view, 2 ns gate

• Surface ionization wave plasma ~ 200 μm thick, wave speed ~ 0.03 mm/ns

• Electric field measured by picosecond four-wave mixing (calibration by electrostatic field)

• Time resolution 2 ns, spatial resolution across laser beam ~ 100 μm

• Objective: electric field mapping in ns pulse discharges in high-pressure fuel-air mixtures

 Electric Field Vector Components
in a Surface Ionization Wave Discharge
Forward Voltage
8 Current 40 30
breakdown - Ex
Absolute field
25

Electric field [kV/cm]

6 Actual field 30 Ey

Electric field [kV/cm]

4 20
(Ex2 + Ey2)1/2
20
-U [kV], I [A]

2 10
15
0 0
10
-2 -10
HV electrode 5
-4 -20
Reverse
breakdown Laser beam 0
-6 -30
-100 0 100 200 300 400 locations -100 -50 0 50 100 150 200
Time [ns] 150 μm from surface Time [ns]

• Initial field offset (at t < 0): charge accumulation on dielectric from previous pulse
• Field follows applied voltage rise, increases until “forward breakdown”
• After breakdown, field reduced due to charge accumulation on dielectric
• Field is reversed after applied voltage starts decreasing
• Away from HV electrode, field peaks later (Ey before Ex): surface ionization wave
• Measurements in a hydrogen-air diffusion flame underway
Plasma property measurements using H2/Ar emission lines
3
D1 ,3 D 2 ,3 D3 , 3
P0 ,3 P1 ,3 P2 , 3
S1 , 1
D 2 ,1 P1 ,1 S0 ,
3s23p54p1
L-S coupling
3s23p6↔ 3s23p54S1
15.75eV ionization
L  l1  l2  1  0  1

14.7eV J  L  S ,... L  S  2,1,0

4d Term : 3
P0 ,3 P1 3 P2

S  S1  S 2  1 / 2  1 / 2  0
13.3eV 3d J  L  S ,... L  S  1
4p Term : 1P1

4s
3s23p6↔ 3s23p54p1

L  l1  l2  1  1  0,1,2 (any quantum number)

1 1
S  s1  s2    0,1
2 2
J  L  S ,... L  S  2,1,0, 1, 2,1,0, 3,2,1
S=0 S=1
Term : 1S0 ,3 S1 , 1
P1 ,3 P0 ,3 P1 ,3 P2 , 1D 2 ,3 D1 ,3 D 2 ,3 D3
 Stark broadening of hydrogen lines and Ar optical emission line-ratio method

Schematic diagram of the experimental setup. The spatially resolved optical

measurement system is shown on the left bottom. On the right bottom is a zoom-in
figure showing the stainless steel needle tips and the discharge gap.
Xi-Ming Zhu, James L Walsh, Wen-Cong Chen1 and Yi-Kang Pu, 15
J.
Phys. D: Appl. Phys. 45 (2012) 295201 (11pp)
Experimentally measured electron densities in a high-pressure
nanosecond pulsed microplasma (Ar/Ne = 700/30 Torr, discharge Uncertainties in the ne measurement (%) using the
current lasts for about 100 ns, pulse period 1 ms). In the legend on Stark roadening method with Ar 696.5 nm line, Hα
the right top, ‘line ratio’ refers to the line-ratio method and ‘Stark line and Hβ line (for ne > 1016 cm−3) and that using
broadening’ refers to the Stark broadening method using Ar 696.5 the line-ratio method.
nm, Hα and Hβ lines with a single-Voigt fitting procedure. In the
legend on the left bottom, ‘centre’ and ‘edge’ denote ne,centre and
ne,edge obtained with double-Voigt fitting. The solid line shows a Xi-Ming Zhu, James L Walsh, Wen-Cong Chen1 and Yi-Kang Pu, J.
function, ne = 6 × 1018 × exp(−(t/0.15)0.22), where ne and tare in Phys. D: Appl. Phys. 45 (2012) 295201 (11pp)
units of cm−3 and ns, respectively.
16
 Measurements of temperature, vibrational level populations
and fast heating using picosecond CARS

Air, P=100 Torr

2 mm
10
mm

t= 1-10 μs (frames are 1 μs apart)

ns pulse discharge and afterglow: Air vs. nitrogen, P=100 Torr

• Compression waves formed by “rapid” heating, on sub-acoustic time scale, τacoustic ~ r / a ~ 2 μs

• What processes control other features of temperature rise (e.g. “slow” heating”)?
 Comparison with modeling predictions in air:
vibrational kinetics and temperature rise

• Strong vibrational excitation in the discharge, N2(v=0-8)

• Tv(N2) rise in early afterglow: V-V exchange, N2(v) + N2(v=0) → N2(v-1) + N2(v=1)

• Tv(N2) decay in late afterglow: V-T relaxation, N2(v) + O → N2(v-1) + O , radial diffusion

• “Rapid” heating: quenching of N2 electronic states, N2(C,B,A,a) + O2 → N2(X) + O + O

• “Slow” heating: V-T relaxation, N2(X,v) + O → N2(X,v-1) + O

• “Rapid” heating: pressure overshoot , compression wave formation

• NO formation: dominated by reactions of N2 electronic states, N2* + O → NO + N

A. Montello, Z. Yin, D. Burnette, I.V. Adamovich, and W.R Lempert, Journal of Physics D: Applied Physics 46 (2013) 464002
 Single-shot measurement rotational and rovibrational energy distributions
by Hybrid fs/ps coherent anti-Stokes Raman scattering (CARS) spectroscopy
fs ps

(Four wave mixing and fs broadband dual pumping)

ωp1: Rovibrational Raman transition (Q-branch, Δv=+1, ΔJ =0)
ωp2: pure rotational Raman transition (S-branch, Δv=0, ΔJ =+2)
ωprobe: frequency-narrowed ps probe pulse

The He∕N2 dielectric barrier discharge

Dedic, C.E., Meyer, T.R. and Michael, J.B., 2017. Single-shot ultrafast coherent anti-Stokes Raman scattering of vibrational/rotational
nonequilibrium. Optica, 4(5), pp.563-570.
Experimental measurements of Plasma chemistry and Kinetic Processes

NRP discharge in
4.5 mm
air at 1000 K, 1 atm:
NRP spark • Measured quantities:
• 10-ns pulse grounded
discharge
• O atoms: TALIF with absolute calibration (Xe)
• 5.7 kV electrode • N2 (A): CRDS
• N2 (B) and N2 (C): OES
• 10 kHz • Temperature: OES
• Gap: 4 mm • Electron density: Hb Stark broadening

• 670 mJ/pulse Preheated air

at 1000 K

Courtesy of Prof. Christophe Laux

 Measurements of V, I, temperature, densities
40 6 300
5 250

Voltage (V)

E/N [Td]
Current [A]
30
4 V Iconduction 200
20 3 150
10 2 100
1 50
0 0 0

Temperature [K]
2500 Temperature
Ultrafast heating: from N2(C-B)
hheating =21±5%
900 K in 20 ns 2000 from N2(B-A)
1500
18
1.2x10
18 3
Absolute densities [cm-3]

1.0x10 O ( P) density
Ultrafast 8.0x10
17

dissociation of O 17
26.0x1017 hdiss. = 35±5%
4.0x10
17
2.0x10
0
17
10
10
16
N2(B)
15
10
14
10 N2(A)
10
13
N2(C)
12
10
-10 0 10 20 30 40 50
Time (ns) Rusterholtz et al, J. Phys.D, 46, 464010, Dec 2013
 Summary of processes involved in
flame stabilization by NRP discharges
Chemical effects:
RH + O  R + OH
O2 Oxidation
2O

N2(X) N2(A) O2
e- N2(B) N2(X) + 2 O + E
N2(C)

T

Thermal effects
2-step mechanism (Popov, 2001):
N2 + e → N2* + e (N2* = N2 A, B, C, …) 5 μs after pulse
Thresholds: 6.2, 7.4, 11.0 eV
(Xu et al., APL. 99, 121502, 2011)

N2* + O2 → N2 + O + O + T
T = 1.0, 2.2, 5.9 eV
 2. Measurements of Chemical Processes in Plasma Assisted Ignition and Combustion

O atom measurements by using TALIF

O atom mole fraction

Air
5.0E-5
Air-ethylene, =0.5 Atomic O production
600
4.0E-5 525

Extinction strain rate (1/s)

3.0E-5 O (3P) 450
375
2.0E-5
300

1.0E-5 225
no plasma
150 with plasma (f=5 kHz)
with plasma (f=20 kHz)
0.0E+0
1.0E-7 1.0E-6 1.0E-5 1.0E-4 1.0E-3 1.0E-2 0.30 0.31 0.32 0.33 0.34 0.35 0.36
Time, seconds Fuel mole fraction Xf
O atom formation in a plasma discharge of O atom formation in a ns plasma Extension of extinction limit by
air and air-C2H4 mixture in a flow reactor discharge of methane/air plasma discharge
counterflow flames

J. Uddi et al. 2009 W. Sun et. Al. 2010

 OH measurements in a flow reactor:
Plasma chemical reactions result in ignition
End View
Pulse #10 Pulse #100
H2 – air, ϕ=0.3
T0=500 K, P=100 torr
C2H4 – air, ϕ=0.3
T0=500 K, P=100 torr

50 pulses

H2-air, ϕ=0.4

Short burst: OH transient rise and decay

Long burst: plasma assisted ignition, Tignition ≈ 700 K < Tauto-ignition ≈ 900 K
 Plasma Assisted Combustion: Change of ignition and extinction S-curve
The effect of kinetic enhancement (μs ~ ms, 800-1200 K)
New “S-curve” by Plasma assisted combustion for
small molecule fuel such as H2, CH4

15
7x10
O2=34%
Extinction 15
O2=62% CH4

OH number density (cm )

6x10

-3
Temperature

Plasma generated 5x10

15

species: 15
Smooth
4x10 Transition
O, H, O2(a∆g) … Extinction
15
3x10
plasma S-curve
Plasma the classical S-curve 2x10
15

Ignition 1x10
15

Ignition

Residence time 0.05 0.10 0.15 0.20 0.25 0.30 0.35

Fuel mole fraction

Scramjet, afterburner P = 72 Torr, f = 24 kHz, a = 240 1/s

•Strong kinetic enhancement at intermediate temperature
•Less effect at high temperature

Sun et al. Proc. Comb. Inst. 34, 2010, Combust. Flame 2011, 2012
Ombrello et al. 2008
 Plasma assisted low temperature combustion
Methane vs. Dimethyl ether (DME)

P = 72 Torr
f = 24 kHz

OH* emission ~310 nm

30 ms gate
Laser beam
OH, CH2O PLIF

25.4 mm

E = 7500 V/cm, E/N ~ 900 Td

Peak Voltage
= 7.8 KV Power ~ 17 W (repetitive pulses)

26
 OH PLIF measurements in Dimethyl ether (DME) Ignition

S-shaped ignition and extinction curves

DME vs. CH4
Top burner Direct image
OH
fluorescence
at Q1(6)

Flame
Bottom burner (fuel)

OH density vs. fuel mole fraction XO = 0.55, P = 72 Torr, f = 24 kHz, for DME (a = 250 1/s) and CH4 (a = 400 1/s,) as the fuel,
respectively (solid square symbols: increasing XF, open square symbols: decreasing XF)

What is the role of plasma before ignition of DME?

 Plasma activated LTC: change of S-Curve

P = 72 Torr, a= 250 1/s, f = 24 kHz P = 72 Torr, a= 250 1/s, f = 34 kHz,

XO2=40%, varying Xf XO2=60%, varying Xf

5
5 6x10 increase
6x10 LTC decrease
5
5x10
5
Extinction
5x10 LTC

CH2O PLIF (a.u.)

CH2O PLIF (a.u.)

5
HTC
4x10 5
4x10
5
3x10 5
HTC 3x10
5
2x10 Hot Ignition
5
increase 2x10
5
1x10 decrease
5
0.00 0.02 0.04 0.06 0.08 0.10 0.12 1x10
Fuel mole fraction
0.00 0.02 0.04 0.06 0.08 0.10 0.12
S-Curve Fuel mole fraction
New ignition/extinction curve without
extinction limit

Radical production by plasma can activate LTC at much shorter timescale, lower
pressure and temperature; and enable new flame regimes 28
 Flow reactor studies of plasma assisted
low temperature reaction

Comparison of thermal (1, 10 atm) and plasma (1 atm)

propane fuel consumption as a function of temperature.

Nicholas Tsolas, Jong Guen Lee and Richard A. Yetter, 2015, Phil. Trans. R. Soc. A 373: 20140344.
Plasma activated Cool Flames :A new way to burn with plasma
1 n-heptane


Ignition delay time (s)

1 atm
5 atm
10 atm
20 atm
0.1
Plasma

0.01
Extinction
Temperature

0.8 1.0 1.2 1.4 1.6

Plasma 1000K/T
HTC generated

LTC
LTC Ignition

t2 t1 Residence time
t2<< t1

Plasma activated LTC at much shorter

time, lower pressure….

We can create cool flames even at

1 atm or below? 30
2. Plasma activated self-sustaining Cool Diffusion Flames:
A new way to burn

Fuel/N2 @ 550 K
Tf~1900 K
Heated N2 @ 550 K

Stagnation
plane (a) Hot diffusion flame
N2 @ 300 K

Oxidizer @ 300 K
with plasma discharge Tf~650 K
Fig. 1 Schematic of experimental setup
(b) Cool diffusion flame

Fig. 2 Hot and cool n-heptane

diffusion flames at the same condition

Won et al., 2014

 3. Multispecies diagnostics and kinetic modeling
Fuel(RH)
+OH
+O2 Small C2H3/CH2O
RO2 R
e +O2=O+O(1D) +e alkene
+O2
H+O2(1Δg) =O+OH QOOH HO2 H/HCO
O(1D)+RH =OH+R +O2+(M)
N2(A,B,C)+O2=O+O+N2 O2QOOH H2O2
+O2 CO/CO2
N2(v)+HO2 =OH+O+N2
R(v,*)+O2=RO+OH 2OH
Plasma
=??? Slow e, R*, N2*, O2*
R(*), R(v), N2(v), O2(v)
O3+M =O+O2+M

A schematic of the key reaction pathways for high pressure fuel

oxidation of at different temperatures
(blue arrow: Below 700K; yellow arrow: 700-1050 K;
red arrow: above 1050K).
 Question to aks:
How does the key plasma reactions affects n-heptane dissociation and oxidation
in the first 10 ms with efifferent excitation processes involving Ar and O2?

e +Ar =Ar*+e
e +O2 =O+O(1D)
e +O2 =e + O2(v)
e +O2 =O2-
e +C7H16 =H+R
......
Ar* +C7H16 =?
O(1D) +C7H16 =?
O2(v) +C7H16 =?
O2(v) +C7H15 =RO2(v)
......
 Plasma chemistry reactor
• Reactor
• Pressure: 60 Torr
• Initial Temperature: 300 K
• Flow speed: 40 cm/s Mini-Herriott cell
• Reactor size: 45 x 14 x 152 mm3 showing 24 pass
configuration
• Nanosecond repetitively
pulsed discharge: FID GmbH Vacuum
Chamber
FPG 30-50MC4
• Peak Voltage: ~7 kV
• Pulse Duration: 12 ns
FWHM
Reactor
• Continuous mode: 0 – 5 kHz
• Pulse burst mode: 150
pulses, 30 kHz
Ge Etalon
• Quartz double dielectric Macor
Quartz
barrier: 1.6 mm thickness Wall
Wall
Flip Mirror
• Gap distance: 14 mm
Collimating Mirror
Lenses

34
 Experimental Apparatus
Laser inlet
purge tube N2 Purge
Box
Electrode
Connection QCL Laser

Vacuum
Chamber

Observation
Window

Alignment
Detector Laser

To Vacuum
Direct and ICCD Images of Plasma Discharge in a Reactor

Stoichiometric mixtures: C2H4/O2 with 75% AR, 60 Torr, Vmax= 6 kV

•Direct Image: 1 kHz, 3.6 mJ/pulse, 2 s exposure time.

•ICCD images: Gate time = 100 ns, Gain = 250

Cathode
Direct 1000 Hz
Anode
ICCD
1000 Hz

2000 Hz

3000 Hz
 Experiment/Model Comparison (C2H4)
Oxidation Pyrolysis
500 600
Oxidation C2H2 Pyrolysis C2H2
Oxidation CH4 Pyrolysis CH4
Oxidation H2O C2H2 Model
C2H2 Model CH4 Model
CH4 Model 500
400 H2O Model
Mole Fraction (ppm)

Mole Fraction (ppm)

400
300

300

200
200

100
100

0 0
0 0.005 0.01 0.015 0 0.005 0.01 0.015
Time from first pulse (s) Time from first pulse (s)

• Conditions for Oxidiation case: • Conditions for pyrolysis case:

– Stoichiometric C2H4/O2 – 0.9375 Ar/0.0625 C2H4
– 25% Reactants, 75% Argon
– V=40 cm/s
– V=40 cm/s
– P=60 Torr – P=60 Torr
– Ti=300 K – Ti=300 K
– 150 Pulses, 30 kHz, 10 kV – 150 Pulses, 30 kHz, 10 kV
– E/n = 37 Td – E/n = 37 Td
• Acetylene measurements in the pyrolysis • Methane is greatly over-predicted by the
experiment were used to match the E/n ratio model
of model calculations
 Continuous Plasma - Oxidation

• GC and in situ IR-LAS

measurements 1000000
– Continuous plasma
– Stoichiometric, 25% 100000

Mole Fraction (ppm)

Reactant
• Excellent agreement 10000

between GC and IR-

LAS 1000

• Major Species 100

– CH2O
O2 C2H4
– CO CO CH2O
10
– C2H2 C2H2 C2H2 LAS
CH4 CH4 LAS
– CH3OH, C2H4O CH3OH/C2H4O
1
– CH4 0 1000 2000 3000 4000 5000
Plasma Frequency (Hz)

Plasma activated low temperature fuel oxidation is an important process

38
Lefkowitz et al. 2014
 Validation of plasma combustion chemistry of C2H4 oxidation
C2H4/O2/Ar: 6.25/18.75/75

HP-Mech (a) USC-Mech (b)

Fig. 16 Measurements and predictions of C2H2, CH4, and H2O concentrations after 150 pulses at 30
kHz repetition rate for a mixture of C2H4/O2/Ar: 6.25/18.75/75 by using (a) HP-Mech and (b) USC-
Mech II.
 Ethylene Oxidation Pathways
+ Ar(+) 13%
+ OH 15% C2H4 C2H3+

CH2CH2OH + Ar* 5%
+ e- 30% + e- 65%
+ O2 100% LTC
O2C2H4OH
HTC C2H2 C2H
100%
2 CH2O + OH + O 11% + O 13% + O2 100%

+ H +M 31%
+ H 21% HCO + CO
H + CH2CHO CH3+ HCO
C2H5
+ O 21% + O2 46% + O2 + M 85%

+ O2 + M 97% CH20 + HCO CO + CH2O + OH CH3O2

+ X 95%
C2H5O2 M = Third body collider Blue = Plasma CH3O
X = Radical Red = High temperature, + X 96%
+ HO2 98% Green= Low temperature

C2H5O2H CH2O
PAC activates C2H4 low temperature chemistry
Challenges: Plasma activated CH4/O2/He oxidation at 400 K

250

CH2O Experiment
CH2O Model
CH4
Mole Fraction (ppm) 200

+ O(1D) 2%
150

100

CH3 + OH CH2OH + H
50

+ O2 100%
0
0
Pulse
5
Burst10 15 20
Time (ms)

CH2O + H2/HO2
 300 pulse burst, 8.75 kV peaks voltage, 30
kHz pulse repetition frequency
 Stoichiometric mixture, 75% diluent, 60
Torr, 300 K initial temperature
Lefkowitz, J.K., Guo, P., Rousso, A. and Ju, Y., 2015. Phil. Trans. R. Soc. A, 373(2048), p.20140333.
 Challenges: Time dependent Formaldehyde
Measurement
• CH2O production slightly
greater for stoichiometric 70

case
60
– Similar linear trend as in

Mole Fraction (ppm)

n-heptane consumption 50
• Lack of further CH2O production
40
after pulse burst: indicates that
only plasma produced species are
responsible for fuel oxidation 30

• Mode significantly under- 20

predicts measurement by ϕ=1
10 ϕ = 1, Modeling
a factor of 40. ϕ = 0.5
• Missing pathways for low 0 ϕ = 0.5, Modeling
0 0.002 0.004 0.006 0.008 0.01
temperature CH2O
Time (s)
formation

Rousso, Aric, Suo Yang, Joseph Lefkowitz, Wenting Sun, and Yiguang Ju.
" Proceedings of the Combustion Institute, Vo.36, 2017, Pages 4105–4112
42
Cross-sections database available for electron-molecule collisions

Andrey Starikovskiy and Nikolay Aleksandrov, AIAA paper-2017-1977

 Non-equilibrium plasma kinetics
A new reaction pathway of plasma assisted low temperature combustion
via excited RO2*(v)
Vibrational and electronically excited O2(v) collides with fuel radical (R) forming highly
energized RO2* in which the vibrational energies are quickly redistributed due to strong
coupling between different vibration states. These RO2*, comparing to those formed by
ground state oxygen with R, carry much higher internal energies that enable them
overcome the barrier TS1, TS2, and TS3 much easier to produce the bimolecular product
HO2+Alkene and OH+Ether. Therefore, the ignition processes/species profiles of the
mixture can be significantly different from the ground state system.

• Increase the rate

R+O2*(v,e) • Change branching ratios (TSi)
R(v)*+O2 Energized • Modify pressure dependence
RO2*(v)
TS3
R+O2
TS2
TS1 CH2O + OH + R”

HO2+alkene
QOOH
RO2 OH + RO (cyc-Ether)
Rousso, Aric, Suo Yang, Joseph Lefkowitz, Wenting Sun, and Yiguang Ju.
" Proceedings of the Combustion Institute, Vo.36, 2017, Pages 4105–4112
 O(1D) reaction kinetics: Photolysis reactor for elementary kinetic study
O(1D)/O3/O2/CH3OH/Ar mixture
time dependent measurements of OH, HO2, CH2O,…

H2O:1338.55 cm-1
CH2O:1726.8 cm-1 913 mm long, 40 mm diameter, multi-pass (21)
OH: 3568 cm-1 Herriott cell, a19.17 meter optical path length
HO2: 1397 cm-1
Ju, Y., Lefkowitz, J.K., Reuter, C.B., Won, S.H., Yang, X., Yang, S., Sun, W., Jiang, Z. and Chen, Q., 2016. Plasma Chemistry and Plasma
Processing, 36(1), pp.85-105.
 Experimental results and model validation: updated model
O(1D)+CH3OH→CH2OH+OH knew=1.5x1014 (R27)
O(1D)+CH3OH→HOCHOH+H knew=0.5x1014 (R28) (in mole-cm3-s)
O(1D)+CH3OH→CH2O+H2O knew=1.0x1014 (R29)

Time-resolved mole fraction of H2O in the 266 nm laser photolysis of 0.224% CH3OH 1.91%
O2 and 596 ppm O3 in Ar mixture with the variation of CH3OH flow rate compared to model
simulations for 1.0 ml/hour flow of CH3OH. ○: Experimental measurement; ― : simulation
using the original model. ― : simulation with updated reaction rates [38].
Ju, Y., Lefkowitz, J.K., Reuter, C.B., Won, S.H., Yang, X., Yang, S., Sun, W., Jiang, Z. and Chen, Q., 2016. Plasma Chemistry and Plasma
Processing, 36(1), pp.85-105.
 Potential energy surface (PES) of O(1D)+CH3OH
using M062X/cc-pvtz level

New channel?

Ju, Y., Lefkowitz, J.K., Reuter, C.B., Won, S.H., Yang, X., Yang, S., Sun, W., Jiang, Z. and Chen, Q., 2016. Plasma Chemistry and Plasma
Processing, 36(1), pp.85-105.
HO2/OH using mid-IR Faraday Rotational Spectroscopy

VRMS ( )   GP0 sin 2   RMS

+Bfield

Laser Lock-In Amplifier

Paramagnetic (radical) species

HO2 energy levels

Zeeman splitting
Absorption

Dispersion


 Experimental results: HO2/OH measurements

Signal
HO2
Sensitivity OH

 detection limit  1 ppmv / 𝐻𝑧 3 detection limit = 20 ppbv / 𝐻𝑧

DME flow reactor

model validation

Implication:
RO2→QOOH→O2QOOH uncertainty
HCO+O2=HO2+CO reaction uncertainty and HCO formation pathway?
Bremfield et al., 2013, JPC letters, 2013; Kurimoto et al. 2014
 Summary
• Time-resolved, spatially-resolved, in-situ laser diagnostics of electric field, electron
density, and electron temperature, excited and radical species greatly enhanced the
understanding of plasma kinetic and chemical process in PAC.

• Production of O and O(1D), O2(singlet), N2(*), and N2(v) by the plasma is the major
processes in kinetically enhancement of combustion.

• Fast and slow heating in PAC is important, but the energy transfer processes are very
complicated.

• Plasma activated low temperature combustion pathways and enable cool flame
formation, but existing mechanisms have large uncertainties, especially for large
hydrocarbons.

• Electron impact reaction cross-sections of large alkanes and reaction rates involving
O(1D) and non-equilibrium excitations are poorly known.

• The effect of vibrational species excitation on PAC is still poorly known.

Lecture 6 Plasma Kinetic Modeling

Yiguang Ju

Contribution of Xingqian Mao

 6.1 Electron impact reactions
• E/N (electric field / molecule number density):
Electron impact reactions depend on the electron energy
distribution function (or electron temperature) and the electron
collisional cross-section areas of a molecule. The electron energy
is controlled by the electric field strength, E/N.
electronic e
excitation
AB* +

CH4:O2:He=1:2:9 vibrational AB(v) + e

excitation

dissociation
e + AB A + B* + e

ionization
AB+ + 2e

CH4:O2:N2=1:2:3.76 ionization
dissociation
A+ + B + 2e
Energy loss fraction: Fraction for net energy loss per
unit time in each individual collision process k.
Energy loss coefficients(eV m3/s) attachment
 AB-
2me k T f 0
For elastic collisions: K k    [ k ( 2 f 0 + B )]d  ……
Mk 0 e 
For inelastic collisions: K k  U k kk   (2e / me )1/ 2 Mk: particle mass of target particles of collision process k.
𝑈𝑘 : threshold energy of inelastic collision process k.
Cross-sections database available for electron-molecule collisions

Andrey Starikovskiy and Nikolay Aleksandrov, AIAA paper-2017-1977

 Electron impact electronic excitation reactions

electronic
e + AB
excitation
AB* + e

e + N2 → e + N2(A, B, a’, C)

e + O2 → e + O2(𝑎1Δg),
O2(𝑏1Σg+),
O2 *

Reaction rate coefficients for each

individual collision process k (m3/s)

kk     k f 0 d
0

𝜀: electron energy.
f0: isotropic part of electron distribution
function, corresponding to zeroth-order term of
Threshold energy diagram
spherical harmonics expansion in velocity Impact on combustion
space.
𝜎 k: cross section areas of electron-neutral N2(A,B,a’,C) + O2 → N2 + 2O
collision process k. H + O2(𝑎1Δg) → OH + O
• Cross section area data comes from Lxcat database. (www.lxcat.net) H + O2 = OH + O
 Electron impact dissociation reactions
dissociation
e + AB A + B* + e

e + O2 → e + O + O Impact on combustion
e + O + O(1D)
e + O + O(1S)

e + CH4 → e + CH3 + H
e + CH2 + H2
e + CH + H2 + H
e + C + 2H2
New channel?

Threshold energy diagram

 Electron impact ionization reactions
ionization
e + AB AB+ + 2e

ionization
e + AB A+ + B + 2e
dissociation

e + O2 → 2e + O2+

e + N2 → 2e + N2+

e + CH4 → 2e + CH4+
2e + CH3+ + H
Electron and ion production

Impact on combustion: radical production and gas heating

Threshold energy diagram
dissociation recombination
N2++ CH4 → N2 + H + CH3+ e + O2 + → O + O
Electron attachment reactions
attachment
e + AB AB-

e + O2 → O 2 -

e + O2 → O + O-

Impact on combustion ~5 orders of

O2- + H → OH- +O magnitude
O2 + H = OH + O faster at 1000K
 Electron impact vibrational excitation reactions
vibrational
e + AB
excitation
AB(v) + e

e + N2 → e + N2(v) e + O2 → e + O2(v)
e + CH4 → e + CH4(v) e + CO2 → e + CO2(v)

Vibrational modes of CH4

Mode Energy(eV)
v1 Symmetric stretching 0.362
v2 Twisting 0.190
v3 Asymmetric stretching 0.374
v4 Scissoring 0.162

Vibrational modes of O2 Vibrational modes of N2 Vibrational modes of CO2

Energy(eV) Energy(eV)
v1 0.19 v1 0.29
𝜈 1:symmetric valence
v2 0.38 v2 0.59 vibrations;
v3 0.57 … … 𝜈 2: double degenerate
symmetric deformation
v4 0.75 v8 2.35 vibrations;
𝜈 3:asymmetric valence
vibrations.
1. Mi-Young Song, et al. J. Phys. Chem. Ref. Data, 44 (2015).
2. Tomas Kozak and Annemie Bogaerts. Plasma Sources Sci. Technol., 23 (2014).
3. Alexander Fridman. Plasma Chemistry, (2008).
 Reactions involving vibrationally excited species

V-T relaxation
AB(v=n) + M → AB(v=n-1) + M
e + AB
electron vibrational
excitation V-V relaxation.
AB(n) + C(m-1) → AB(n-1) + C(m)
M(v), M* + AB N2(v), O2(v),
excited CO2(v), RH(v),…
Chemical
species AB(v) + C → AC + B
reaction
excitation consumption

• VT relaxation: Energy transfer between vibrational and translational degrees of freedoms.

• VV relaxation: Energy transfer between vibrational states of different molecules.

• Vibrationally excited reaction: A reaction between a vibrational excited molecule with a

neutral molecule. It will increase the reaction rate, add new reaction pathways, and modify
heat release rate.
 Reaction rate of vibrationally excited molecules:
Fridman-Macheret α-Model
An effective decrease of activation energy
FA BC
Ea  Ev
FA BC  FAB C
FA BC  1 Ea (1)

FAB C  2 Ea (2)

 1Ea (1) Ea (1) 1 /  2  1

 
 1Ea (1)   2 Ea (2) Ea (1)  ( 1 /  2 ) Ea (2)
Ea (1)
  (1)
Ea  Ea (2)

𝐹𝐴+𝐵𝐶 𝐹𝐴𝐵+𝐶 : are characteristic slopes of the terms A+BC and AB+C.
𝛾: reverse radii of corresponding exchange forces.
A + BC → AB + C 1: forward reaction direction.
2: reverse reaction direction.
A + BC*(Ev) → AB + C
Ea   Ev • The efficiency α of vibrational energy is the highest for
k ( Ev )  AT exp( n
) strongly endothermic reactions with activation energies
T
closest to the reaction enthalpy (close to 100%).
If Ea>0, the overall
activation energy • The efficiency of vibrational energy is the lowest for
decreases exothermic reactions without activation energies (close to 0).
Alexander Fridman. Plasma Chemistry, (2008).
Vibrational-translational (V-T) relaxation for diatomic molecules

AB(v=n) + M → AB(v=n-1) + M
AB(v=1) + M → AB + M
Values of parameters for rate coefficients of the
processes of N2(v1) + M → N2 + M
 V-T relaxation rate constants of diatomic
molecules M n m A B C D
N2 1 1 7.8x10-12 218 690 1
k10 (cm3 / s )
H2 1 2/3 4.9x10-12 167.1 394 1
1
B C  E10 
 AT n exp(  )  1  D  exp(  ) H2O 1 0 2.5x10-15 21.18 0 0
 T 
1/3 m
T T CO2 1 1 1.1x10-12 218 690 1

The probability of the vibrational transition n in one collision

Em  En
Pnm  Pmn exp( )
k BT Values of parameters for rate coefficients of the
Em  En processes of O2(v1) + M → O2 + M
knm  kmn exp( )
k BT M n m A B C D
Single-quantum transitions with probabilities O2 1 0 1.35x10-12 137.9 0 1

Pn 1, n  (n  1) P10 H2 1 0 2.69x10-12 91.5 0 1

He 1 0 4.54x10-15 60.85 0 1
kn 1,n  (n  1)k10 Ar 1 2 3.14x10-12 173.1 6.2x105 1

M. Capitelli, et al. Plasms Kinetics in Atmospheric Gases, (2000).

V-T relaxation and V-V relaxation for non-diatomic molecules
 V-T relaxation rate constants calculation  V-V exchange rate constants
by Schwartz, Slawsky, and Herzfeld (SSH
theory) AB(n) + CD(m-1) → AB(n-1) + CD(m)
F ( n ) F ( nm )
kn ,n 1  k1,0 Z n
F ( 1 ) knm,n1,1m  k1,0
0,1
Zn Zm
F ( 11 )
1  xe
Zn  n Vibrational energy transfer reactions of CO2
1  nxe
Reaction xe(10-3)
1 2  2
F ( )   3  exp(   )  exp(  ) CO2(va) + M → CO2 + M
2 3  3 0.0 a

0.32E  1/2 CO2(v1) + M → CO2(va) + M 3.7 b

n  ( )
 Tg CO2(v1) + M → CO2(vb) + M 1.0 b
E  En  En 1 CO2(v1) + M → CO2(vc) + M -15.6 b
  17.5 / r0 CO2(v1) + CO2 → CO2(vb) + CO2(va) 2.8
𝑍𝑛 : scaling factor.
CO2(v1) + CO2 → CO2(va) + CO2(vb) 17.6
𝑥𝑒 : the anharmonicity of the energy levels.
𝛾𝑛 : a parameter which is a measure of adiabaticity of the
CO2(v1) + CO → CO2 + CO(v1) 5.25;6.13
reaction.
𝛼 : a parameter of the exponential repulsive potential
between the colliding species. a. Multiply by 1.0, 0.7 and 0.7 for M=CO2, CO and O2
𝜇 : reduced mass of the collision species. b. Multiply by 1.0, 0.3 and 0.4 for M=CO2, CO and O2
𝑟0 : the
。 radius parameter of the Lennard-Jones Potential. 1. M. Capitelli, et al. Plasms Kinetics in Atmospheric Gases, (2000).
3.94 A for CO2. 2. Tomas Kozak and Annemie Bogaerts. Plasma Sources Sci. Technol., 23 (2014).
3. Alexander Fridman. Plasma Chemistry, (2008).
4. M. Capitelli. Nonequilibrium Vibrational Kinetics,(1986).
 6.2 Simulation of plasma assisted combustion
Challenges:
• Multi-physics problems: photons, electrons, electronic and
vibrational excitations
• Electromagnetic field, acoustic waves, shockwaves, ignition and
combustion waves
• Multi-length scale: Sheath (Debye length), diffusion & reaction
zones, mixing layer, far field effects
• Multi-timescale: from plasma, reactive flow, to combustion (ps-ms)
• Multi-species and stiffness of reactions
• Multi-dimension and far non-equilibrium (not Maxwellian
distribution)
• …
Kinetic description of plasma
– distribution function
– kinetic equation
– collisional energy transfer

Multi-fluid description of plasma

– Fluid conservation equations
– Input of transport coefficients and rate coefficients, e.g. fitted as functions of E/N
– Different energy modes (electrons, vibration, neutral temperature)
– Coupling between neutral and charged particles
 Plasma Modeling: A Hybrid ZDPlasKin-CHEMKIN Model (0D)
 BOLSIG+: incorporated in ZDPlasKin[1], a
computer program for the numerical solution of
the Boltzmann equation for electrons in weakly
ionized gases in uniform electric field.
Boltzmann equation for electron energy
distribution function (EEDF) in a plasma is:
f e
   f  E   v f  C[ f ]
t m
𝑓 : electron distribution in six-dimensional phase space.
𝜐 : velocity coordinates.
𝛻v : the velocity-gradient operator.
C : the rate of change in 𝑓 due to collisons.
Assume: a steady state problem and the electric field and
the collision probabilities are spatially uniform, on the
scale of the collisional mean free path (elastic collision),
using the two-term approximation expansion in spherical
coordinates
f (v,cos , z , t )  f 0 (v, z , t )  f1 (v, z , t )cos
 ZDPlasKin[1]: A zero-dimensional Plasma Kinetics
solver which is a Fortran 90 module designed to follow
isotropic part anisotropic part
the time evolution of the species densities and gas
(temporal dependence) (spatial dependence)
temperature in a non-thermal plasma with an arbitrarily
complex chemistry. 𝑣 : the magnitude of the velocity.
𝜃 : the angle between the velocity and the field direction.
[1].S. Pancheshnyi, B. Eismann, G.J.M. Hagelaar, L.C. Pitchford, Computer code ZDPlasKin,
http://www.zdplaskin.laplace.univ-tlse.fr (University of Toulouse, LAPLACE, CNRS-UPS-INP, z : the position along this direction.
Toulouse, France, 2008).
[2]. G.J.M. Hagelaar, et al. Plasma Sources Sci. Technol., 14 (2005). BOLSIG+ solver: https://www.bolsig.laplace.univ-tlse.fr
 Plasma Modeling: A Hybrid ZDPlasKin-CHEMKIN Model
 Time evolutions of the species and temperature
of plasma reactions using ZDPlasKin
dN i jmax
  Qij (t )
dt j 1

Deposited energy: Pext  Pgas  Pelec  Pchem

3 d( N eTe )
Translational degree of Pelec  kB
freedom of electrons: 2 dt
1 d ( NT )
Translational degree of Pgas  kB
freedom of the gas:  1 dt
dNi
Internal degree of Pchem    i
freedom of the gas: i dt
 Time evolutions of the species and temperature
in combustion using CHEMKIN
• A combined ns-plasma and DC discharge to 
dYi
 iWi
dt
control electronic and vibrational excitation and
i max
dT
selective reactivity in plasma. C p    eiiWi
dt i 1

 Coupling of ZDPlasKin and CHEMKIN

in 1 -in Tn 1  Tn
 icombustion  iplasma  Tcombustion  Tplasma
t t

i* -in T *  Tn
 Tcombustion
 icombustion t
t
Tn 1  T *
in 1 -i*  Tplasma
 plasma
t
t
i

1. A.E. Lutz, et al. SEKIN, (1998).

2. P.N. Brown, et al. J Sci. Stat. Comput., 10 (1989). 𝜌𝑖∗ and 𝑇 ∗ are the intermediated density and temperature between
3. J.K. Lefkowitz, Y. Ju, et al. Phil. Trans. Soc. A, (2015).
the calculation of CHEMKIN and ZDPlasKin.
0-Dimensional model approximation
• Plane to plane
• Homogenous
• Square wave of the E/N

Kinetics model
• Species 100

• Reactions
Plasma kinetics: 541(ZDPlasKin)
HP-mech: 595 (C1-C2) (CHEMKIN) Mean electron energy

Modelling conditions
• Pressure 60torr
• Temperature 373K
• Frequency 30kHz
• Pulse number 300
• NSD 180Td
• DC 1Td, 5Td, 10Td, 20Td
• Mixture 8.333% CH4, 16.667% O2, 75% He

http://engine.princeton.edu/mechanism/HP-Mech.html Electron number density

Time evolution of the vibrationally excited species
Reaction rates in the first 20 pulses at 20Td
 1D Plasma Modeling
 Kinetic model of Nagaraja et al. (2015)
 𝛁 ⋅ 𝝐𝛁𝜙 = −𝑒 𝑛+ − 𝑛− − 𝑛𝑒
𝜕𝑛𝜖
 +𝛁 ⋅ 𝑱𝜖 = 𝑄𝜖
𝜕𝑡
𝜕𝑛𝑘
 +𝛁 ⋅ 𝑱𝑘 = 𝜔𝑘
𝜕𝑡
𝜕𝜌 𝜕𝜌𝑢𝑖
 + = 0
𝜕𝑡 𝜕𝑥𝑖
𝜕𝜌𝑢𝑖 𝜕(𝜌𝑢𝑖 𝑢𝑗 ) 𝜕𝑝 𝜕𝜏𝑖𝑗
 + =− + + 𝐹𝑖𝐸𝐻𝐷
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗
𝜕𝜌𝐸 𝜕[ 𝜌𝐸+𝑝 𝑢𝑖 ] 𝜕𝑞𝑖 𝜕 𝑢𝑖 𝜏𝑖𝑗
 + =− + + 𝑄 𝐽𝐻
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑖 𝜕𝑥𝑗

 Transport coefficients of electron and rate

coefficients of electron impact reactions: BOLSIG
𝜙: potential (voltage) 𝑢𝑖 : 𝑉𝑒𝑙𝑜𝑐𝑖𝑡𝑦
𝑛𝑘 : 𝑛𝑢𝑚𝑏𝑒𝑟 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑜𝑓 𝑘𝑡ℎ 𝑠𝑝𝑒𝑐𝑖𝑒𝑠 𝐸: 𝑇𝑜𝑡𝑎𝑙 𝑒𝑛𝑒𝑟𝑔𝑦
𝑛𝜖 : electron energy density, 𝑛𝑒 *E𝑒

Multi-scale modeling
 Adaptive time-stepping: small (10-13 - 10-12 s)
during each discharge pulse; larger (10-10 s) in the Fig. Time evolution of short lived
gap between 2 consecutive pulses. electronically excited species after a
Nagaraja, S.,Yang, V. and Adamovich, I., 2013.. Journal of single nanosecond pulse in air at 60 Torr
Physics D: Applied Physics, 46(15), p.155205. and 300 K
 Multi-timescale modeling method

Multi-timescale (MTS) Method

t
1 
   dYi  k
 i    
 Yi  dt 
Yk  K k e
40
t = 0.1 ms
35 t = 0.2 ms
30 t = 0.3 ms

Number of species
ΔtF 25
Fastest
Group
Fast species 20
ΔtM 15
Medial
Groups 10
5
0
-1 -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 -12 -13
ΔtS
Slowest Log10(characteristic time / s)
Group

0 1 2 Δt
Time
Diagram of multi time scale scheme
ΔtF is the time step of the fastest group, ΔtM is the time step of the medial
group, and ΔtS is the time step of the slowest group

Gou, Sun, Chen and Ju., Combust. Flame, 2010, 2013

 Hybrid Multi-Time Scale (HMTS) Method

• HMTS method: Initial calculation

Calculate characteristic
times

Calculate group number

and time steps

Group calculation

Parameters update

No Yes Output
End ? results

X.L. Gou, W.T. Sun, Z. Chen, Y.G. Ju, Combust. Flame 157 (6) (2010) 1111–1121.
21
 HTMS Validation: Homogeneous ignition
Ignition in
homogeneous
n-decane/Air 121 species (M. Chaos, IJCK,2007) mixture

2.8 2.0
temperature VODE
5 2.6

Ignition delay time (ms)

10 CO2 MTS
1.5 1atm

Temperature (1000K)
0
2.4 HMTS
10 OH
Mass fraction

2.2 1.0
-5
10 ODE 2.0
0.5
10
-10 MTS 1.8 20atm

-15
HMTS 1.6 0.0
10
1.4
C10H22 -0.5
-20
10 1.2 0.5 0.6 0.7 0.8 0.9
0 1 2 3 4 5
1000/Initial temperature (1/K)
Time (0.1 ms)

Temperature and species profiles Ignition delay time for n-decane-air

Comparison of computation efficiency: same mechanism
with HMTS, HMTS/DAC vs. VODE, VODE/DAC
120
Real Fuel 2_Reduced_425
P = 1.0 atm; Φ = 1.0; T0 = 400 K
100
PFA Time w. DAC
CPU Time

80

Solver's Time (VODE

60
or HMTS)

40 Other Time

20
?
0

1D unsteady Flame propagation

 Methodology
Correlated DAC (CO-DAC) method

• Correlated groups Time correlation

Space correlation

tn+1

tn
How to choose
criteria?

Sun, W., Gou, X., El-Asrag, H.A., Chen, Z. and Ju, Y., 2015. Multi-timescale and correlated dynamic adaptive chemistry modeling of
ignition and flame propagation using a real jet fuel surrogate model. Combustion and Flame, 162(4), pp.1530-1539.

24
 Spherically Propagating Flame
• Computation time

Transport

~50%

Sun, W., Gou, X., El-Asrag, H.A., Chen, Z. and Ju, Y., 2015. Multi-timescale and correlated dynamic adaptive chemistry modeling of
ignition and flame propagation using a real jet fuel surrogate model. Combustion and Flame, 162(4), pp.1530-1539.

25
 Spherically Propagating Flame
• Temperature and Radicals Profile
Zoom in

Shifted 0.3%

9th U.S. National Combustion Meeting, Cincinnati

May 17 – 20, 2015

26
 Correlated Dynamic Adaptive Chemistry and Transport (CO-DACT) Method

• CO-DACT Time correlation

Space correlation

tn+1

tn
How to choose
criteria?
Sun, W., and Ju, Y., 2016. A Multi-timescale and Correlated Dynamic Adaptive Chemistry and Transport (CO-DACT) Method for
27
Computationally Efficient Modeling of Jet Fuel Combustion with Detailed Chemistry and Transport. Combustion and Flame, submitted.
Phase Parameters & Similarity Criteria

• How to choose criteria?

• ηk, λk and Djk T
(T, N2, O2, Fuel, H2O, H2, CO2, CO, CH2O)T
• η, λ and Dmk Xi
• Major species: of N2, O2, Fuel, H2O, H2, CO2, CO, CH2O, C2H4
• In most combustion systems, the summation of these major reactants and
productions account for at least 95% of the molar fraction in total.
• CH2O is dominating in low temperature region.
• Criteria: T T 0
X N 2 X N0 2
X O2 X O0 2
X Fuel X Fuel
0

d  X H 2O X H0 2O 
Error in transport
X H 2 X H0 2
~ O(ε)
X CO2 X CO
0
2
If extrapolated,
X CO X CO
0

X CH 2O X CH
0 Error ~ O(ε2)
2O 

28
 Spherically Propagating Flame
• Computation efficiency:

HMTS/CO-DACT
HMTS/CO-DAC

More than 200 time faster

Sun, W., and Ju, Y., 2017. A Multi-timescale and Correlated Dynamic Adaptive Chemistry and Transport (CO-DACT) Method for
Computationally Efficient Modeling of Jet Fuel Combustion with Detailed Chemistry and Transport. Combustion and Flame, 2017.
Comparison: CPU Time Dependence on species number

t=0.028*Ns2.84

t=0.0002*Ns3.29

t=0.62*Ns1.16

CH Dimethyl ether N-heptane

4
H2

, Ns

Opportunity: adaptive HMTS/G-scheme

 2D fluid model for plasma assisted combustion

Input of transport coefficients

and rate coefficients, e.g.
fitted as functions of E/N

Anne Bourdon, Sumire Kobayashi, Zdenek Bonaventura,Fabien Tholin and Nikolay Popov, Kaust
Research Conference: New Combustion Concepts, March 6-8, 2017, KAUST
Anne Bourdon, Sumire Kobayashi, Zdenek Bonaventura,Fabien Tholin and Nikolay Popov,
Kaust Research Conference: New Combustion Concepts, March 6-8, 2017, KAUST
Time histories of species and temperature

Key reactions for gas heating

• Rapid decrease of ne at the end of the voltage pulse, then a

much slower decrease up to 100 ns
• O(P), H and OH are mainly produced in the post-discharge by
dissociative quenching reactions

Anne Bourdon, Sumire Kobayashi, Zdenek Bonaventura,Fabien Tholin and Nikolay Popov,
Kaust Research Conference: New Combustion Concepts, March 6-8, 2017, KAUST
 7. Future research of plasma assisted reactive flow
Yiguang Ju, 2017
1. Game changers in PAC applications
• Engines (ICEs: lean burn, Turbine engine: relight, ignition)
• CO2 capture and chemicals : CO2 utilization and methane reforming (plasma catalyst)
• Bio-medicine
• Materials synthesis
• …
2. Game changers in plasma control
• High pressure volumetric discharge
• Selective excitation (electronic and vibrational)
• Selective species/radicals production
• Low cost and low electronic noise
• …
3. Fundamental Research of PAC
• Plasma properties: Electric field, electron number density, excitation states, non-equilibrium temperatures
• Plasma physics: Energy transfer processes between different excited states
• Plasma chemistry: Low temperature kinetic pathways, non-equilibrium kinetics
• Kinetic process: Key species, reaction rates, and cross-section areas for large fuel molecules
• Multi-dimensional modeling tools