Aiaa 98-0328: January 12-15, 1998 / Reno, NV
Aiaa 98-0328: January 12-15, 1998 / Reno, NV
Aiaa 98-0328: January 12-15, 1998 / Reno, NV
R=19980022683 2018-07-25T15:48:05+00:00Z
AIAA 98-0328
and
Stephen P. Wilkinson
NASA Langley Research Center
Hampton, VA
For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics
1801 Alexander Bell Drive, Suite 500, Reston, Virginia 20191-4344
AIAA 98-0328
and
Stephen P. Wilkinson t
NASA Langley Research Center, Hampton, VA 23681
Abstract Introduction
Low speed wind tunnel data have been The use of magnetohydrodynamics
acquired for planar panels covered by a (MHD) to control the turbulent viscous
uniform, glow-discharge surface plasma drag due to aerodynamic boundary layer
in atmospheric pressure air known as the flow has received considerable attention
One Atmosphere Uniform Glow over the years. Most concepts have been
Discharge Plasma (OAUGDP). based on ionized flow around a
Streamwise and spanwise arrays of magnetized hypersonic vehicle, or on
flush, plasma-generating surface achieving such a plasma with ion seeding
electrodes have been studied in laminar, techniques. Numerous examples ate
transitional, and fully turbulent boundary found in patents awarded in the 1960's
layer flow. Plasma between symmetric (see references 1 and 2). Emphasis has
streamwise electrode strips caused large been placed on the magnetohydrodynainic
increases in panel drag, whereas approach in hydrodynamics due to the
asyrnmetric spanwise electrode electrically conducting nature of seawater
configurations produced a significant and perceived high economic or
thrust. Slnol/e wire flow visualization performance payoffs. However, in terms
and mean velocity diagnostics show the of a net energy balance, performance
prilnary cause of the phenomena to be a enhancement has proven elusive. An
combination of mass transport and extensive review of pure MHD methods
vortical structures induced by strong for drag reduction (i.e., those based
paraelectric electrohydrodynamic (EHD) exclusively on the cross product of the
body forces on the flow. local current density and magnetic
...................... induction, j x B), through 1989 was
* Professor, Dept. of Electrical compiled by Tsinober (ref. 3). Several
Engineering, Associate Fellow AIAA current investigations are discussed in
references 4 and 5.
' Graduate Research Assistant, Dept. of
Physics and UTK Plasma Sciences An alternative to MHD flow control
LaboIatou which has received far less attention in
the field of boundary layer research is
-i- Aerospace Engineer, Fluid Modeling based on the electric field alone, or
and Control Branch, Senior Member electrohydrodynamic (EHD) control. In
AIAA p_mially ionized gases the electric field
itself, or the pm-aelectric effects associated
1
Copyright ©1998 by the American Institute of Aeronautics and Astronautics. Inc. No copyright is asset'led in the United States
under Title 17, U. S Code. The U. S. Government has a royalty-frec license to exercise all rights under the copyright claimed
herein for Governmental purposes All other rights are reserved by tile copyright owner
with an electric field gradient, can be used charge for practicable engineering values
to accelerate ions and, via particle of magnetic and electric field strengths.
collisions (mobility drift), the neutral gas This is an important point in view of
(references 6, 7). In the past, a difficulty potential aerospace flight applications.
with the EHD approach, especially in The maximum practical magnetic field
non-hypersonic flight applications, is from permanent magnets which can be
generating an energy-efficient ionized expected in flush-mounted, non-
flow near the surface. This report obstructive surface application is
presents experimental data on the first estimated to be no more that about B =
aerodynamic application of a new EHD 0.1 Tesla. While higher values are
method based upon the One Atmosphere obtainable with electromagnets, their
Uniform Glow Discharge Plasma or Joulean dissipation (or superconducting
OAUGDP (refs. 6 to 13), that may refrigeration energy requirements) would
change this situation. seriously compromise any net energy
saving in, for instance, a drag reduction
The OAUGDP is a novel, surface- application. The minimum electric field
generated, atmospheric, RF (radio required to generate an OAUGD plasma
frequency) plasma. The concept of the in air is about E = 1.0 MV/meter (10
device is discussed in detail later. The kV/cm). Assuming a typical commercial
primary feature that distinguishes it from transport flight velocity of U = 300 m/sec,
prior RF plasmas is its efficient ability to the force ratio on a singly charged particle
create a uniform glow discharge at is given by the quotient -- E/UB -- 3.3
atmospheric pressure on an extended fiat xl04. In other words, the electric force on
surface. The electrodes required to do this such a charged particle is more than four
have characteristics which lend orders of magnitude greater than the
themselves to practical engineering maximum practicable magnetic force.
applications, such as simplicity and
robustness. They should also be To examine the ratio of boUv forces, the
inexpensive and reliable. Given this magnitude of the electrical current and
capability, the goals of the present EHD charged particle number densities rnust be
study employing the OAUGDP for considered as well. For the plasma
laminar or turbulent boundary layer considered in this report, a charged
control are twofold: 1)-demonstrate the particle number density, N = 1.0x 10_T/m 3
generation of EHD forces with is characteristic. A maximum current
magnitudes sufficient to alter boundary density corresponding to the glow-to-arc
layer flow dynamics, and 2) demonstrate transition, J = 10 4 A/n-l 2, is assumed as a
that such forces constitute a useful flow value not likely to be exceeded in any
control mechanism. glow discharge plasma application. The
body force ratio is then given by the
Before introducing the OAUGDP and the quotient rB = qN<E/JB where q is the
current flow control study, however. electronic charge. This yields rn = 16, or
some additional discussion of pure EHD an EHD body force more than one order
controls will help to show why this of magnitude greater than that of the
approach has been chosen. An interesting MHD body force.
feature of EHD controls, which perhaps
is not generally appreciated, is that the Another fundamental advantage of EHD
electrostatic force on a charged p,'uticle forces is that the electric field can do work
can be significantly larger than the on the charged particles and. through
magnetic force on the same moving strong collisional coupling at one
2
American Institute of Aeronautics and Astronautics
atmosphere,on the aerodynamicflow pertinent independent variables, which
itself. A static magnetic field of force include the electric field strength (E),
always operates orthogonally to the electrode separation distance (d), type of
chargedparticle velocities,and therefore gas, pressure (p), and RF electric field
cando no work on the particlesor the frequency (Vo). A relation among these
flow. For aerodynamic flow control variables is given by Roth (section 12.5.2
applications,it is evidentthatEHD is the in reference 6)
preferred approach. The obvious
questionsare how to effectivelyproduce Vo _ E/(pd) (1)
therequisiteelectricallycharged medium
atoneatmosphere,andhow to configure for the case of a parallel plate geometry.
and drive the electric fields to produce A planar strip geometry will have a
effectsthat may be useful in such areas similar but more complicated relation due
as drag reduction,heat transfer,lift, or to the arched field lines, but the same
flow separation. qualitative functional dependencies would
be expected to prevail. The electric field E
An adequatenumber densityof charged in Equation 1 may be approximated by
particles can be produced in an the electrode potential, V, with E=V/d.
atmospheric glow discharge using a Provided the operating parameters are in
recentlydevelopedtechnologyproprietary accordance with Equation 1, the
to the University of TennesseeResearch OAUGDP will function at one
Corporation.While undergrantfrom the atmosphere and produce a stable, steady-
Air Force Office of Scientific Research, state glow discharge. A plasma thickness
the first author was successful in of one or two millimeters at power
producing the OAUGDP (Ref. I1), densities well below one watt per cubic
which is an extremely uniform, low centimeter was typical for the current
frequencyRF glow dischargethat does experiments.
notrequireeithera vacuum environment
or the mega- or gigahertz supply Equation 1 does not represent a finely
frequencies typical of industrial RF tuned phenomenon and the parameters
plasmas.The OAUGDP operateson the can vary over a useful range while
principle of the charge-trapphzg maintaining the existence and uniformity
mechanism. Charge trapping refers to a of the plasma2 If any of the parameters
specific, constrained, periodic oscillation deviate significantly from Equation 1,
of ions and/or electrons along electric field however, either the OAUGDP will cease
lines between a pair of fiat electrodes to function, or its uniformity will degrade
which are characteristically side-by-side in into a filamentary discharge. The
flat panel aerodynamic applications. This sensitivity of the OAUGDP to variations
electrostatic trapping may reduce plasma in the independent input parameters is a
polarization, keep ions from knocking current subject of investigation at the
secondary electrons off the instantaneous University of Tennessee's Plasma
cathode (which may initiate avalanches or Sciences Laboratory.
breakdown), and prevent ions from
healing the cathode surface and initiating a The magnitudes of the parameters in
glow-to-arc transition. Equation 1 for bench top demonstration
of the OAUGDP are easily attainable.
Based on straightforward Lorentzian For instance, a frequency of several
electrodynamic analysis of the plasma, the kilohertz, an rms voltage of several
charge trapping mechanism identifies the kilovolts, and a planar strip separation
3
American Institute of Aeronautics and Astronautics
distanceof 5 or 10 mm are adequate to principle, the lowest possible energy cost
initiate the plasma at atmospheric of producing an ion-electron pair in a
pressure. The OAUGDP is not hard- plasma source, and compares very
starting, and does not require external favorably with the energy cost of other
initiation with a Tesla coil or spark gap. atmospheric plasma sources, such as
While the dissipative (or plasma) current plasma torches or arcjets, for which the
in the OAUGDP is small (about 0.030 energy cost is about 10,000 eV/ion-
amp rms in these experiments), without electron pair.
special impedance matching the reactive,
non-dissipative current can be large Regarding applications, the OAUGDP is
(approximately 0.4 amp rms) and the quenched by liquid water, although it
power source must be sized accordingly. recovers rapidly from a water spray. The
usual ranges of atmospheric, climatic
The absence of any large dissipative humidity conditions are acceptable
currents due to filamentary breakdown or although high, near condensing levels
arcing in the OAUGDP plasma allows it have not been investigated. Only dry,
to operate at low power levels, consistent high altitude applications are currently
with the possibility of net energy savings envisioned.
in flight boundary layer flow control or
drag reduction applications. For example, The OAUGDP is fundamentally different
a characteristic boundary layer viscous from ion wind concepts that rely on a
dissipation value for a long range corona discharge as an ion source. Malik,
commercial transport has been estimated et al. (ref. 14)used the ion wind technique
to be roughly 5000 watts per square in a fiat plate DC "brush" discharge
meter (737-class airplane at cruise fashion and were able to secure a small
conditions). By comparison, in bench top reduction in measured drag of about 5%
tests, the OAUGDP can operate with a for a turbulent boundary layer flow at a
power of 320 W/m 2 or less based on the length Reynolds number of
measured, non-reactive power and the approximately one million. Research was
surface area covered by the plasma. later abandoned, however, due to inability
While there is no evidence or claim at this to scale the operation of the hardware to
time that such a low power level can flight conditions. More recently, EI-
effectively control, say, a turbulent Khabiry and Colver (ref. 15) were able to
boundary flow at high Reynolds number produce up to 50% or more viscous drag
flight conditions, the energy cost of reduction in very low Reynolds number
sustaining a uniform layer of glow flows (on the order of 105) using a corona
discharge plasma over a large area is discharge between spanwise wires on a
nonetheless very low. flat surface for both DC and low
frequency (60 Hz) AC excitation. Each
This low energy cost occurs for a of these techniques is probably limited to
fundamental reason: the OAUGDP has
low Reynolds number applications due to
been shown to be a glow discharge, limitations on scaling the corona
created twice during each RF cycle (see discharge effect to higher flow velocities.
refs. 12, 13). As a glow discharge, the The OAUGDP, however, is more readily
ionization process in the instantaneous scaleable and has the potential to function
cathode region occurs at the Stoletow at much higher Reynolds numbers.
point, which is about 81 electron-
volts(eV) per ion-electron pair for air With an efficient source of surface
(Ref. 6, Section 8.3.4). This is, in plasma, the challenge becomes how to
4
American Institute of Aeronautics and Astronautics
effect a useful EHD flow control out. Any small remaining residual force
mechanism in a boundary layer, is well within the linear range of the
particularly a turbulent boundary layer. instrument and is accounted for in the no-
Initial investigations were aimed at flow drag tare readings.
understanding the basic response of a
boundary layer to several simple, planar The smoke wire was 0.1 mm diameter
electrode configurations that can be used type 304 stainless steel and was stretched
to produce the OAUGD plasma. These across the width of the test section at a
consist of streamwise and spanwise variable height above the wall. A weight
arrays of flush-mounted strip electrodes and pulley arrangement kept the wire taut
on a fiat panel, all at the same RF potential during heating. It was powered by a
and phase with respect to a ground plane variable DC power supply with a 100 vdc
or electrode on the opposite side of the maximum output (typical range at 4 m/s
panel. was 40-50 vdc). The "smoke" was the
vapor of common mineral oil. Smoke
Experimental Apparatus wire photographs were obtained by firing
an electronic flash during the vertical
Low speed wind tunnel tests of panels blanking period of a full frame,
with the OAUGD plasma were conducted monochromatic digital video camera (8-
in the NASA Langley 7xl 1 Inch Low bit resolution, 768 by 484 pixels), at a
Speed Wind Tunnel (7xll) to determine variable delay time after energizing the
the basic response of boundary layer smoke wire. The delay time was
flow to the plasma for a few simple panel determined by trial and error. Video
configurations. The 7xll is a closed pixel data were downloaded from the
return, unpressurized air facility with a digital camera to a computer for
test section 178H x 279W x 914L processing.
millimeters. A 305 x 279 millimeter
central portion of the lower test section For velocity profiles, a slender, tapered
wall was used for testing. Tests included total pressure pitot tube was traversed
the directly measured viscous drag of fiat across the boundary layer height
plate panels with the OAUGD plasma downstream of the energized plasma
generated on the surface, vertical (wall- panels. The tip was fabricated from
normal) boundary layer pitot pressure flattened, stainless steel hypodermic
profiles measured a short distance tubing. The tip height was 0.28 mm and
downstream of the panels, and smoke the width was 0.65 mm. The probe was
flow visualization tests. The air-bearing far enough downstream of the energized
drag balance used and a general view of panel to prevent any electrical arcing to
the tunnel test section is shown in Figures the instrument. The initial height of the
la and lb. The semi-catenaries shown in probe above the wall was set by
Figure l a are high voltage power leads to monitoring electrical contact between the
the plasma panels. They consist of brass- probe and metallic wall. The probe was
ball utility chains (commonly used for raised through the boundary layer with an
light switch pull chains, etc.) and were automated stepping motor-driven slide
chosen for their extreme flexibility, mechanism in 0.5 mm increments. A
electrical conductivity, and lack of any typical profile was acquired quickly (in
sharp, corona-producing features. By about 30 seconds) to prevent heating the
exerting equal and opposite horizontal panels, which could cause their adhesive
forces on the drag balance, the forces due backing to weaken and release. Pitot
to the power leads approximately cancel differential pressure was measured
5
American Institute of Aeronautics and Astronautics
between the probe and a static pressure underneath the panel to the other terminal.
port of the tunnel wall with a high The parallel electrode strips on top of the
accuracy capacitive or piezoelectric gauge. panel were generally at high voltage,
while the lower electrode was grounded,
Figure 2 shows a plan view sketch of a although configurations with the opposite
typical panel. The panels were polarity would also produce plasma and
constructed from conventional dielectric the effects reported below. A high
printed circuit board material (woven- voltage (up to 5.4 kV), low frequency RF
glass/epoxy construction, 0.75 mm thick, (up to 20 kHz) power supply was used
double-sided, I ounce copper coating). with its transformer output connected
The plasma-generating electric field lines directly to the panel without a special
arch over the upper surface of the board impedance matching network.
(where the plasma is generated) and
traverse the board thickness. In the more Figure 4 is a plan view photograph of a
recent designs, an array of electrode strips panel energized (plasma activated but with
was etched on the top (flow) side of the no flow) and is representative of the
board and the bottom surface left as a technique. The 0.5 mm solid, horizontal,
uniform copper plane as illustrated in dark strips are the parallel copper
Figure 2. Alternately, an asymmetric electrodes. The gray-scale regions to
array of electrode strips was etched on the either side of the electrodes are the
bottom of the panel when it was desired OAUGD plasma. The plasma was
to accelerate or decelerate the boundary visually extremely uniform.
layer flow. The geometrical
configurations of the various panels used For the drag tests, the panel was mounted
in this study are shown in Figures 3a on an air bearing drag balance located
through 3c. below the tunnel test section, with the
panel forming the central section of the
For all tests, the flow passed over the lower wall. The boundary layer flow was
copper electrodes with no additional tripped near the outlet of the tunnel's
dielectric coating. Since the OAUGDP contraction with a 1.07 mm circular rod
charge trapping mechanism operates on on the test wall 575 mm upstream of the
displacement rather than real electrical leading edge of the panel. Small (0.25
currents, this surface can, if desired, be ram) gaps around the test panels allowed
covered with a thin insulating and/or them to float freely on the drag balance.
protective layer without qualitatively A pressure control box around the test
affecting the results reported herein. The section allowed the static pressure in the
circuit board was attached to a 12.7 mm test section to be matched to the control
thick fiberglass backing board (type G- box pressure. This minimized errors in
10) with double-sided adhesive tape to drag measurements by reducing flow in
make the panel structurally rigid but still the gaps surrounding the panels.
capable of being disassembled. The
designation code and electrode Procedures and Results
dimensions of the various panels reported
on in this paper are listed in Table 1. Data for streamwise and spanwise
electrode orientations were acquired, as
The parallel electrode strips on the top of well as paired comparison drag data for
the panel were bussed together and both the plasma-energized and
connected to one power supply terminal unenergized (approximate smooth flat
and the lower plane or electrodes plate drag) conditions. Data were also
6
American Institute of Aeronautics and Astronautics
taken on panels with asymmetric arrays unsteady oscillations) as evidenced by
of electrodes such as those shown in smoke wire pathline visualization and the
Figure 3b to study the acceleration and absence of any turbulent breakdown in
deceleration of the flow in the boundary diagnostic hot wire signals.
layer, and the consequent drag decrease or
increase (respectively) compared to the Representative results from a panel with
unenergized fiat plate. Data on drag symmetric electrodes, each electrode a
increase or decrease were measured as copper strip 0.5 mm wide with centers
parametric functions of the flow velocity spaced 10.5 mm apart, are shown in
(up to 26 m/s), electrode voltage (up to Figure 5a for the streamwise electrode
5.4 kV rms), and RF frequency (from orientation (panel C7-A), and in Figure 5b
500 to 8000 Hz). for the spanwise electrode orientation
(panel C7-C). Each of these displays the
The direction and magnitude of the expected power-law Reynolds number
paraelectric plasma-induced acceleration dependence for the "plasma off"
of the flow is determined by the direction condition. Note the change in slope of the
of the electric field gradients, and these are "plasma off" curve in Figure 5a or 5b in
in turn strongly influenced by the the range of 7-8m/s, corresponding to
orientation and details of the electrode transition from laminar to turbulent flow.
geometry. The preliminary data reported For the "plasma on", streamwise
here are for unoptimized electrode electrode case, a substantial increase in
geometries. It is anticipated that with drag is observed. This is due to several
additional modeling studies, geometrical factors. As will be shown, the plasma
optimization will increase the magnitude excitation for velocities below about 7 m/s
of the effects reported at a given set of (laminar region) trips the flow to full
plasma operating parameters. In addition, turbulence, partially explaining the drag
the electrodes in this study were energized increase in that region. The drag increase
with a single phase of RF excitation. This persists, however, to the highest attainable
produces EHD body forces which are the velocity of the wind tunnel indicating that
result of averaging attractive and repulsive more than just flow tripping is involved.
forces over the RF cycle, a second order For the "plasma on", spanwise electrode
effect. Much stronger effects should be case, a smaller drag increase is produced
possible when adjacent electrodes are and only in the laminar/transitional region.
excited with polyphase RF power, The difference in behavior between the
providing a DC electric field parallel to the two cases along with evidence presented
surface, afirst order EHD effect. later in the paper suggests the formation
of strong, EHD-driven, streamwise
In this paper, data are presented for three vortical structures in the boundary layer
principal cases: laminar data, for which for the streamwise-oriented electrode
the wind tunnel flow was laminar before case.
encountering the panel; transitional data
corresponding to about 75% intermittency The very small differences in surface
at the upstream edge of the model, and configuration among different panels did
fully turbulent data. Since the boundary not measurably affect (beyond the
layer flow was tripped upstream of the intrinsic precision of the data) the drag for
panel, there was actually no case of the unenergized panels reported in this
completely undisturbed laminar flow. At paper. Despite the small roughness
low tunnel velocities, however, the flow introduced by the copper electrodes on the
was laminar (but with occasional panel surfaces, relative to the energized
7
American Institute of Aeronautics and Astronautics
. - . . _,?/- = _fj +._ ",,
8
American Institute of Aeronautics and Astronautics
Vertical boundarylayer velocity profiles between the pressure sensor and A/D
werealsomeasuredwith a total pressure converter ranges. The trend of the data is
probe on severalpanelswith symmetric valid.)
as well as asymmetric electrode
configurations.Figure8 presentsvelocity The profiles corroborate the drag and
profilesfor Panel C7-A (with symmetric, smoke wire data. For the streamwise
streamwise electrodes) one-half way electrode case, there is a substantial
betweentwo adjacentelectrodes,for the retardation of the profile affecting the
caseof laminar (a), transitional(b) and entire boundary layer. This increases the
fully turbulent (c) flow at the panel boundary layer momentum deficit and
leadingedge. The probe tip was located qualitatively corresponds to the large
approximately one boundary layer increase observed in the drag in Figure 5a.
thicknessdownstreamof the model over For the spanwise electrode configuration
the smooth aft filler plate of the lower shown in Figure 10, a significant effect is
wall. (A metallicaft platewasusedfor the evident only in the laminar regime, with a
profile measurements to aid probe initial similar effect on the drag (Figure 5b). For
height determination; for drag the smoke wire flow visualization, the
measurements, a non-metallic plate was eruption of vortical structures observed in
used to minimize electrostatic drag error.) Figures 6a and 6b appears to be consistent
Figure 9 presents similar data, also from with the flow retardation observed in the
the streamwise electrode configuration, velocity profiles of Figures 8 and 9.
with the pitot probe directly behind one of
the streamwise electrodes. Figure 10 Figure 11 shows the instantaneous RF
shows the profiles downstream of voltage and current for panel C1-C
spanwise oriented electrodes on panel C7- operated at an rms voltage of 1.4
C. kilovolts, and a frequency of 2.5 kilohertz.
The voltage was measured at the power
The profiles for the streamwise case supply output with a high voltage probe
(Figures 8 and 9) show a dramatic having the requisite frequency response.
alteration of the flow due to interaction The current through the high voltage
with the plasma that diminishes with power cable was measured with a high
increasing velocity. There is a large bandwidth, toroidal current transformer
acceleration of the flow near the wall and with a sensitivity of 1 volt/amp. The
a retardation farther out. The cases of the noisy region at the positive peaks of the
probe between and behind the electrodes current waveform represents the plasma
are qualitatively similar, but differ in initiation, during which a classical, "DC",
magnitude. Smoke wire (e.g., Figures normal glow discharge briefly exists
6(a,b)) and hot wire diagnostics show that between the electrodes (Refs. 12, 13). The
the energized, streamwise electrode plasma ignition appears only once per
patterns effectively trip the flow, and that cycle for the model and conditions
any between-electrode/behind-electrode portrayed in Figure 11. For most models
differences are largely mixed out at the studied during these tests, however,
end of the panel. For the spanwise case in plasma ignition occurred twice per cycle
Figure 10, the effect is largely limited to (see Introduction). There was a noticeable
the laminar flow condition, with little variability in the current waveforms for
effect in the transitional case and virtually the various panels and excitation voltages,
no discernible effect in the turbulent case. which are the subject of ongoing study.
(The step-wise appearance of the data in
Figure 10(a) is an error due to a mismatch
9
American Institute of Aeronautics and Astronautics
: •7:%:) ,5 :!,_j! _il_ : ,,,=_s_.
10
American Institute of Aeronautics and Astronautics
propulsive. For the caseof the co-flow commonly used white flow marker
orientation, however, a trend exists below chemical) injected manually in a slow,.
10 m/s indicating a clear Reynolds steady stream from a plastic squeeze
number dependence. The plasma has been bottle. The plasma is not visible in Figure
noted in all cases to trip the boundary 16 due to the strong illumination required
layer so the Reynolds number for the smoke. The paraelectric forcing in
dependency shown in Figure 14 could be Figure 16b causing the jet to deflect
more boundary layer trip related than towards the electrode is evident.
turbulence modification related.
Nonetheless, this finding along with other In terms of a phenomenology, the flow of
data presented in this report point to the the smoke and the air which it marks
possibility of using the newly discovered responds to paraelectric EHD effects in
EHD forcing to target and control the following way. In Figure 16b, the
boundary layer turbulence. flow is drawn downward by a low
pressure above the low electric field
Model E6-C was not optimized for the gradient region of the plasma, entrained in
EHD force. While the predominant the ion-driven plasma flow toward the
plasma forms on the upper surface over region of high electric field gradient, and
the lower surface electrode, flow forced outward by the region of high
visualization has shown that a small (plasma stagnation) pressure along the
amount of plasma forms on the opposite surface of the panel. The flow is rapidly
edge of the upper surface electrode due to accelerated away from the region of high
field lines wrapping around to the lower gas pressure and high electric field
electrode. The net effect is to have a large gradient (primarily to the left of the
EHD force in one direction (downstream electrode due to the asymmetry but also to
in the co-flow case) and a smaller force in a lesser degree to the right as well). This
the opposite direction. effect is responsible for the blowing
velocity profiles illustrated in Figure 15.
The asymmetric panel E6-C was
mounted in the wind tunnel without flow, The behavior shown in Figure 16b is
but with the pitot tube positioned at the consistent with a pure paraelectric effect
same location used in Figures 8-10. The on the plasma and on the flow which it
resulting blowing velocity profiles are entrains. It is not a classical case of
shown in Figure 15 for electrode voltages dielectrophoresis, although similarities
of 3, 4, and 5 KV rms. Maximum exist. Dielectrophoresis refers to the
plasma-induced velocities up to 4.0 forces on neutral, polarizable, dielectric
meters/sec were observed. Particularly material when subjected to a spatially
interesting were the induced velocities of non-uniform or a time-varying electric
up to 0.5 meters/sec at distances at least 3 field (ref 16). In the current case, no
cm from the wall, which occurred for "all smoke or air movement is observed until
driving voltages. sufficient voltage is reached for the
plasma to initiate. This indicates a
Figure 16a and 16b are photographs of different phenomenon than
the influence of the OAUGDP on a dielectrophoretic behavior alone. It is
laminar jet of smoke injected above a clear that the underlying mechanism for
single, asymmetric electrode arrangement. neutral gas movement in the presence of
The test was conducted in a still air the OAUGD plasma warrants further
chamber. The "smoke" in this case was theoretical and experimental study.
actually titanium tetrachloride (a
11
American Institute of Aeronautics and Astronautics
"4 ;_<. . _ ,,,_'.
12
American Institute of Aeronautics and Astronautics
layer temperature downstream of an the ability to move a neutral gas with
energized model showed only a small EHD forcing to reduce or enhance drag,
temperature rise of several degrees or significantly alter the velocity profile
Celsius. A more pertinent question of the boundary layer. Beyond this,
would be the magnitude of the localized future efforts need to be directed at
electron temperature within the plasma passing through Bushnell's second filter,
and its impact on the observed that of technological feasibility. This
phenomena. However, this is beyond the entails demonstrating such factors (where
scope of current investigations. not already demonstrated) as simplicity,
economy, retrofittability, mechanical
Figures 8 and 9 show that the effect on passivity, and robustness and reliability.
the plasma is spread across the entire It also means demonstrating the ability to
boundary layer for the streamwise simulate the processes and mechanisms
symmetric electrode case. It seems clear in such a way as to make possible
that a major vortex-dominated developmental work on small inexpensive
mechanism is in play. This is evidenced models in ground facilities. In addition,
by direct manipulation of the streamwise much future work needs to be done to
flow by EHD forces in the (initially) characterize, parameterize, and understand
laminar smoke wire data shown in Figure the physical processes both in the
6. OAUGDP and with respect to the
paraelectric EHD effects responsible for
Future Plans the plasma-flow interaction.
13
American Institute of Aeronautics and Astronautics
observed flow phenomena attributed to Engineering Conf., Forum on Control of
electrohydrodynamic (EHD) forcing of Transitional and Turbulent Flows, July 7-
the flow by a paraelectric RF body force. 11, 1996.
14
American Institute of Aeronautics and Astronautics
(i ;_(_;__>: i__!i_:"_ _,_,_, _::_-:
!> <_._t :_<_v_:.__% __ _:_¢_,:.
13. Massines, F.;Ben Gadri, R.;Rabehi, 16. Pohl, H. A.: "Dielectrophoresis, The
A.; Decomps, Ph.; Segur, P.; and Behavior of Neutral Matter in
Mayoux, Ch.: "Mechanisms of a Glow Nonuniform Electric Fields", Cambridge
Discharge at Atmospheric Pressure University Press, 1978, ISBN 0-521-
Controlled by Dielectric Barrier", Journal 21657-5, pp. 6-18.
of Applied Physics. (accepted and in
press, 1997). 17. Jung, W., Mangiavacchi, N. and
Akhavan, R.: "Suppression of Turbulence
14. Malik, M. R., Weinstein, L.M., and in Wall-Bounded Flows by High
Hussani, M. Y., "Ion Wind Drag Frequency Spanwise Oscillations", Phys.
Reduction", AIAA Paper 83-0231, 1983. Fluids A. vol. 4(8), pp 1605-1607, 1992.
15
American Institute of Aeronautics and Astronautics
Power leads
Figure la. NASA Langley 7xl 1 Inch Low Speed Wind Tunnel. Plasma panel
(unenergized) on drag balance in wind tunnel test section (front and top walls of test
section removed for clarity)
Tunnel Contraction
Tunnel
' [ 915 [ diffuser
""--__ N\\\\\\\\\\\\\\\\\\\\\\\\"L .......
I
Airflow _ll_* 177.8
I_ 305
I--.- 305
16
All dimensions mm
46 270
14.3
'_2_-" ..................................
i
! II
_!I
bottom planar
electrode
(dashed outline)
i
[---) ;i
273
J !i
1 mm bus bar !
__
/
_ ................................
I
i
14.3
x" /
2 mm bus b." 26 electrodes 0.5 mm width, 10.5 mm center spacing
17
dielectric
circuit boar_e__,_._._n.__
(b)__
NOT TO SCALE
0.5 mm
_k___
Electrode
Figure 4. Portion of Model C7-C with plasma. E=3 kVrms, F=3 kHz.
(Original photograph taken with 35mm 400 ISO color print film,
approximately 10 second exposure.)
18
30 I l I I I I I I I I
I0
m
m
{0
E
1
_3b
a
0.1
Plasma Off
m
0.01 I I I I I I I II I
I 0 20 30
Velocity (m/s)
30
I I I I I I I I I I
10
m
O3
E
Plasma On
0")
v
tm
0.1
Plasma Off
m
0.01 I I I I I I I I I I
10 20 30
Velocity (m/s)
19
Fforward bus bar
10.5 mm
Figure 6(a,b,c). Smoke wire flow visualization of panel C7-A (streamwise upper
electrodes with planar lower electrode) at two excitation voltages. U_ = 4 m/s. Smoke
wire at Y=5 mm (u/U_- 0.65). Images digitally enhanced with unsharp mask filter.
20
leading tip of electrode
I I I I
E 0.8
03
0.6
o
o_
t_
c_ 0.4
o
O
09
m 0.2
0 I I I
1000 2000 3000 4000 5OOO
Voltage(V rms)
Figure 7. Electrostatic drag for panel C 1-B (streamwise upper electrodes over planar
lower electrode); F= 1.5 kHz, U_ = 0.
21
35
U=4m/s
E25 .............
::.:..v.::...P.!.a...s._.
"',,, 20 .......................
_-.......................
!.......................
_......................
. _....
(a)
15 .......................
_........................
i.......................
_.....................
a
>- 10 .......................
_........................
_........... '......
35 , ,
i U=7m/s
30
.........
___-0_ ................
_..................
_...................
8
........... ..............
:...i
.......................
i.....................
i i i i; w _
(b)
15
>- 20 ......................
10 _........................
:: i .......................
_........................
:: .,.
..............
:: _ ::
5 :: f
35
i ! i
U=26rrv'sl
3O
25
E
E
20
(c) ,-
15
121
10
0
0.0 0.2 0.4 0.6 0.8 1.0
Velocity ratio (u/U)
--- 25
E
E
v
_) 20
(b
(a) t"
._ 15
a
.............
--.-,-.-i...._!.a...,.....=...._
.............
I'"
>" 10
.......................
i............
_-i ...........
_-i.-.-
0
0.0 0.2 0.4 0.6 0.8 1.0
Velocity ratio (u/U)
35 i i
_ U=7m/s _
3O
(b)
E
E
L,)
t-
121
25
2O
15
....i.
".'.i.'i.
i.I"i..i......,
...............
,,-.-.,
............
......................
......................................
'...!....
.........
10
0
0.0 0.2 0.4 0.6 0.8 1.0
35 ! i
U=26m/s
3O
.............
_ i_ias_iSn
..................
......................
i...................
_'"
25 ..........
':_.I:..P'as_O"
...............
...........
E
E
v
20
(c) t-
.o_ 15
>- 10
0
0.0 0.2 0.4 0.6 0.8 1.0
23
35 !
u=4m/s
30
................
_"Pias'_" On
---v-'_ Plasma _ ,
!
......................................
•.-.. 25
__.....1_.
8 2o
(a) ..................
!........................
i.......................
_.......................
i,-.-
._._ r-
C
_ 15
i5
_- 10 ........................
i........................
i.......................
!.......................
iL:_/_
i i i _IF _"
5 i _....___---_ i
0 i#r _i-- i
0.0 0.2 0.4 0.6 0.8 1.0
Velocity ratio (u/U)
35
U=7rn/s
30
25
.................
_Pia_s__Oii
---v--i Plasma Off
......................
i............................
i................
i i ii]
E
g
@ 20
(b) e'-
.__ 15
o
_" 10
0
0.0 0.2 0.4 0.6 0.8 1.0
Velocity ratio (u/U)
35
U=26m/si
3O ..............
_ i_ias_On .................................................................
---v.i- Plasma Off
25 ....................... i ....................... .......................................................
E
E
20
,,J
(c) t-"
15
0
>- 10
0
0.0 0.2 0.4 0.6 0.8 1.0
Figure lOa.,b.,c.) Pitot tube velocity profiles for Panel C7-C (with symmetric,
span-wise electrodes; 3.0 kHz, 5.0 kV). Pitot tube located 28 mm downstream
of last electrode.
24
i i i i
2 - V/lO00,"l .l_ -J
/ _ , I I/lO / /\ IMh \ ]
><
O
., ,,("J"\)l
- / f_\ ,' 'Y.D_t"I...
ID
r,.)
O
ID
-1
_
la ma
Formation
ex0
O
> -2
t I I I
0 2re 4_
Phase
0.2
I I I I I
Figure 12. Electrostatic drag for panel E6-C in still air. Electrodes are spanwise,
asymmetric with lower electrode staggered downstream relative to the upper electrode.
Frequency: 3 kHz. Negative values correspond to a reduction in drag.
25
15
I ! I I I I
10 m
_o Plasma On:
E Counter-Flow
v
13t) N
m
m
n n
0.1 I I I t t J I I t
2 4 6 8 10 20 3(
Velocity (m/s)
1.1 I I I I I I
-_ _ A_?,
_'0 0.9
,,,j Lf -', -
,_ 0.8
i_
c5
0.7
f _C °-_,°w
_- 0.6
-- B-- - Counter-Flow
0.5 I I I I I I
0 5 10 15 20 25 30
Velocity (m/s)
Figure 14. Drag change due to OAUGDP for panel E6-C in Co-flow and Counter-flow
configurations (See Figure 13). Voltage and frequency were 4kV rms and 3kHz
respectively.
26
35 I ' I
30
25
20
E
v
E 15
>-
10
0 - 3 kV 4kV 5 kV _
, I , I , I I
0 1 2 3 4
Velocity (m/s)
Figure 15. OAUGDP-induced velocity (blowing) along panel E6-C surface, normal
to the electrodes. Frequency was 3 kHz. Measurements conducted in still air.
27
laminarjet of titanium
tetrachloride vapor
~15mm
/_
(a) Plasma Off
28