Numerical Simulations of Linear Single Dielectric Barrier Discharge Plasma Actuators

1
Numerical Simulations of Linear Single Dielectric
Barrier Discharge Plasma Actuators

• Daniel A. Reasor, Dr. Raymond P. LeBeau
• University of Kentucky
• Dr. George P. Huang
• Wright State University
• Dr. Ylidirim B. Suzen
• North Dakota State University

This work is supported by Kentucky NASA EPSCoR
through the Kentucky Space Grant Consortium.
32nd Annual Dayton-Cincinnati Aerospace Sciences
Symposium
2
Plasma Actuators
M. Post, T. Corke, Flow Control with Single
Dielectric Barrier Plasma Actuators AIAA Paper
2005-4630, 35th AIAA Fluid Dynamics Conference
and Exhibit, 2005.
3
The Linear DBD Plasma Actuator
plasma region
flow
body force
exposed electrode
dielectric
embedded electrode
V
ac
4
Model of DBD Plasma Actuator
From Maxwells Equations
Assume
Gausss Law states
Note that
permittivity of free space
where
5
Model of DBD Plasma Actuator
Combining yields two equations

is the Debye length (length scale of plasma),
Where
is the net charge density.

Enloe, McLaughlin, et al. Mechanisms and
Responses of a Single Dielectric Barrier Plasma
Actuator Geometric Effects, AIAA Journal, Vol.
42, No. 3, 2004.
6
Potential Due to Electric Field BCs
Fluid Region
Outer Boundaries
Dielectric Region
7
Charge Density BCs
Fluid Region
Outer Boundaries
Dielectric Region
8
Effects of the Plasma Actuator
Body Force Vector
body force
9
Model Parameters
where
Air
Material Properties
Kapton
Empirical Value
Experimental Input
Arbitrary Parameter
10
Kentucky Fluid Clusters
KFC6A 23 AMD Athlon 64x2 4600 KFC6I 24
Intel Core 2 Duo e6400 KFC5 47 AMD Athlon 64
3200 KFC4 48 AMD XP 2500 KFC3 32 Intel
Pentium 4 2.8GHz
11
UNCLE

• Originally written by Dr. George P. Huang while
  at UK
• 2D/3D unstructured, finite volume based,
  incompressible Navier-Stokes solver
• 2nd order upwind scheme for advection terms
• 2nd order centered difference for diffusion terms
• 2nd order fully implicit time discretization
• Momentum - collocated grid with Rhie and Chow
  interpolation method
• Pressure - cell-centered-based method similar to
  SIMPLE
• Turbulence - one-equation Spalart-Allmaras and
  two-equation Menter SST
• 2nd order centered difference for additional
  potential equations

12
Dimensions of Actuator
13
Grid Generation
Grid Generation using Gambit
80,000 grid points
14
Potential Due to Electric Field
UNCLE Results
GHOST Results
x(m)
Y.B. Suzen, P.G. Huang, J.D. Jacob, D.E. Ashpis
Numerical Simulations of Flow Separation Control
in Low-Pressure Turbines using Plasma Actuators.
45th AIAA Aerospace Sciences Meeting and Exhibit,
Reno NV, 2007.
15
Normalized Net Charge Density
UNCLE Results
GHOST Results
Y.B. Suzen, P.G. Huang, J.D. Jacob, D.E. Ashpis
Numerical Simulations of Flow Separation Control
in Low-Pressure Turbines using Plasma Actuators.
45th AIAA Aerospace Sciences Meeting and Exhibit,
Reno NV, 2007.
16
Preliminary Cross Flow Results
Crossflow results with Plasma Actuator
Crossflow results without Plasma Actuator
17
Preliminary Crossflow Results
Crossflow results with Plasma Actuator
Crossflow results without Plasma Actuator
18
Future Ongoing Work

• Verify results through comparisons with
  experimental and numerical data
• quantitative comparisons with experiments in
  quiescent flow

• Simulate separation control on airfoils/turbine
  blades

• Expand to three dimensions
• Implement and compare with experimental data from
  the annular plasma actuator synthetic jet (PASJ)
  by Jacob, Santhanakrishnan.

19
Questions?
20
Continued
Results of the use of plasma actuators at various
positions on the Flying-Wing 1303 UAV Model
increase in lift at high angles of attack
M. Patel, T. Ng, S. Vasudevan Plasma Actuators
for Hingeless Aerodynamic Control of an Unmanned
Air Vehicle' AIAA Paper 2006-3495, 3rd AIAA Flow
Control Conference, San Francisco, California,
June 5-8, 2006.
21
Continued
M. Post, T. Corke, Flow Control with Single
Dielectric Barrier Plasma Actuators'' AIAA Paper
2005-4630, 35th AIAA Fluid Dynamics Conference
and Exhibit, Toronto, Ontario, June 6-9, 2005.
22
Comments on the model

• considers the potential due to electric field
• considers the changes in net charge density
• does not model the change in temperature due to
  the plasma
• model of charge density above embedded electrode
  is preliminary but
• does include the dependence on input voltage

W. Shyy, B. Jayaraman, A. Anderson Modeling of
glow discharge-induced fluid dynamics. Journal
of Applied Physics Volume 92, Number 11, December
1, 2002.
23
Why Plasma Actuators?
initially still air
Can the resultant cause a buoyancy force to be
present?
location of the body force
G.I. Font, S. Jung, C.L. Enloe, T.E. McLaughlin,
W.L. Morgan, J.W. Baughn Simulation of the
Effects of Force and Heat Produced by a Plasma
Actuator on Neutral Flow Evolution.''
24
Continued
F. Thomas, A. Kozlov, T. Corke Plasma Actuators
for Bluff Body Flow Control.'' AIAA Paper
2006-2845, 3rd AIAA Flow Control Conference, San
Francisco, California, June 5-8, 2006.