2. Microplasma actuator

investigated in applications for separation flow control [5–8] and noise reduction [8, 9]. The plasma actuator has an operating principle based on the electrohydrodynamic (EHD) phenomenon occurring due to the momentum transfer from ions accelerated by an electric field to neutral molecules by collision. Dielectric barrier discharge (DBD) and corona discharge are among the most common types of plasma actuators. A single DBD plasma actuator can induce a flow up to 7 m/s, and with a multiple DBD plasma actuator design, the value of induced flow can reach 11 m/s [10]. Various applications of flow control require different types of plasma actuators. For high-speed flow control, an induced flow speed of more than 7 m/s is necessary; thus instead of using a single DBD plasma actuator, the corona discharge could be used [11–15]. In the case of turbulent boundary-layer control for skin-friction drag reduction, millimeter-size discharge gap DBD plasma actuators were energized at peak-to-peak voltages of about 7 kV [16–19]. Research studies regarding the applications of plasma actuators involve turbulent boundary-layer separation control, steady airfoil leading-edge separation control, oscillating airfoils dynamic stall control, and circular cylinder wake control. According to various researchers, high values of the induced flow are desired. High values of the flow are obtained conventionally by energizing the plasma actuators at tens of kilovolts which are difficult to insulate and for which a large sized power supply is necessary. An effective actuation effect requires also a higher EHD force density that can be achieved using a micrometer order discharge gap plasma actuator, which lowers the discharge voltage and consequently requires lower power. Micro-sized plasma actuators were used for separation flow control and drag reduction [20–22]. A microplasma actuator was developed for flow control [23]. Similar electrode configurations were described in [24–28] but the required discharge voltages are more than 2 kV. Microplasma is a type of dielectric barrier discharge nonthermal plasma and its driven voltage is around only 1 kV. This technology could be used as a replacement of conventional technologies for surface treatment of polymers, indoor air treat-

Measurements of microplasma actuators are difficult due to their micrometer size discharge gaps; thus, we have developed a numerical simulation of the induced flow based on the Suzen-Huang [33, 34] and Orlov [35] models. Due to the light emission from microplasma, observation and measurement of flow is difficult. Numerical simulations of the plasma actuator are carried out using the plasma fluid model and particle in cell model [36–39]. A simplified phenomenological model which does not model the species transport equations but can replicate the effects of the actuator in the air is the Suzen-Huang model [40]. Results close to the experimental data were obtained by various researchers who developed numerical simulations based on this model. It is less computationally expensive than solving the species

The microplasma actuator developed in our laboratory is thin and flexible and can be attached to any surface. To energize the actuator, small-sized power supplies are necessary; thus, a potential use for this actuator could be on drones. The electronic switching adds a greater flexibility in order to obtain flows in various directions. In comparison with macro-plasma actuators where the power supplies are bulky, the small-sized power supply necessary for a microplasma actuator adds little weight to the drone; also, less electrical insulation is necessary.

ment, biomedical applications, or flow control [29–32].

transport equations [40, 41].

4 Actuators

A microplasma actuator was developed for flow control. Due to its small size, the experimental results were difficult to obtain near the active electrodes; thus, a numerical simulation was developed in order to simulate the flow and add additional information about the flow. The experimental and numerical simulation results showed the capability of the microplasma actuator to control the flow.

#### 2.1. Experimental study of microplasma actuator

#### 2.1.1. Characteristics of microplasma actuator

A schematic of the typical construction of a plasma actuator and the induced air flow is shown in Figure 1. A pulse or alternating high voltage is applied to the two electrodes with a dielectric layer in between; thus, a plasma is generated at the surface of the electrodes [42]. An electrohydrodynamic (EHD) phenomenon occurs because ions are accelerated by the electric field and furthermore collide with neutral molecules; thus, a momentum transfer occurs and air flow is generated [21, 43].

A microplasma actuator that can be driven by a lower voltage of less than 2 kVpp was developed and investigated. Owing to their low discharge voltage, the applied high voltage could be controlled easily using semiconductor switches. The structure of the microplasma actuator is shown in Figure 2 [44]. A dielectric layer consisted of a polymer film of 25-μm thickness sandwiched in between a high voltage electrode and a grounded electrode. Voltages less than 2 kVpp are enough to generate a microplasma due to the thickness of the 25-μm discharge gap. A pulse high voltage power supply was used to energize the actuator with the schematic shown in Figure 3. Four field-effect transistors (FETs) are used with a DC high voltage power supply. The microplasma actuator is energized by a positive pulse voltage while FETs 1 and 4 are in the ON state and by negative pulse voltage while FETs 2 and 3 are in the ON state, respectively.

Figure 4 shows the experimental setup used to visualize and investigate the air flow induced by the microplasma actuator. The microplasma actuator was set on a Z stage. A high voltage probe (Tektronix, P6105A) connected to an oscilloscope (Tektronix, TDS 3014) was used to

Figure 1. Typical construction of a plasma actuator. High voltage is applied to the electrodes with a thin dielectric layer in between.

Figure 2. Geometry of the microplasma actuator. The four high voltage energized electrodes can be driven independently. (a) Top-view. (b) Cross-sectional view.

1280 1024-pixel resolution and 4000 Hz recording frequency. The area captured by the

Figure 4. Experimental setup for visualization and measurement of the air flow induced by the microplasma actuator.

of 1 mm in the x-axis and 0.5 mm in the y-axis. For the PTV method, steady state images of two

The Stokes number was calculated in order to estimate the ability of the incense smoke to

tracing accuracy with error less than 1% [46]. Ozone concentration due to the microplasma was measured by an ozone monitor (Seki Electronics, SOZ-3300). The measurements were carried out by placing the microplasma actuator inside a glass tube which was connected to the ozone

The microplasma actuator could be energized by bipolar pulse as shown in Figure 5 with T period and also positive and negative pulses as shown in Figure 6(a) and (b). For both positive and negative pulses at a frequency of 2 kHz, the pulse width was 4 μs. The calculated energy consumption by integrating the discharge power waveform shown in Figure 6 for a single pulse is shown in Figure 7. In order to have microplasma discharge, a voltage higher than 0.5 kV was necessary. The values of energy were almost the same for positive and negative pulses respectively; thus, at 1.4 kV, the calculated value of energy was about 800 μJ for both

Energy consumption was measured for a discharge voltage of 1.4 kV at various pulse frequencies for both polarities as shown in Figure 8. At 1 Hz, the energy consumption was about 650 μJ. With the increase of frequency, the single pulse energy consumption increased. This could be explained by the charges in plasma that due to the electric field were moved and trapped at the surface of the dielectric. These surface charges [48, 49], when the next opposite polarity is applied, will increase the electric field strength; furthermore, higher energy consumption for a single pulse was obtained with the increase of frequency. Due to recombination

positive and negative pulses. This corresponds to an average power of 1.6 W.

or diffusion, the surface charge will decrease in time.

follow the flow field [47], and the value was in the order of 10<sup>5</sup>

. The results were analyzed by the PTV method having a resolution

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

http://dx.doi.org/10.5772/intechopen.75802

7

; this is an acceptable flow

camera was 34 11 mm<sup>2</sup>

continuous frames were used.

monitor at a distance of 90 cm.

Figure 3. Schematic of pulsed high voltage source with semiconductor switches.

measure the discharge voltage. For the discharge current measurement, a current probe (Tektronix, P6021) connected to the oscilloscope was used. By multiplying the values of the voltage and current waveform measured with an oscilloscope, the discharge power was obtained for the pulse voltage. Furthermore, in order to calculate the discharge energy, the obtained discharge power waveform was integrated. Particle tracking velocimetry (PTV) was used as the method to measure the induced air flow by the microplasma actuator [45]. The microplasma actuator was placed inside an acrylic box and incense smoke with particles 0.3 μm in diameter [46] was inserted inside the box. For the measurement of the incense smoke particles, a laser particle counter (Kanomax, Model 3886) was used. For the flow measurements, a laser was used (Nd:YVO4 532-nm) to irradiate the particles. The flow of tracer particles was recorded by a high-speed camera (Red Lake, Motionscope M3) having Dielectric Barrier Discharge Microplasma Actuator for Flow Control http://dx.doi.org/10.5772/intechopen.75802 7

Figure 4. Experimental setup for visualization and measurement of the air flow induced by the microplasma actuator.

1280 1024-pixel resolution and 4000 Hz recording frequency. The area captured by the camera was 34 11 mm<sup>2</sup> . The results were analyzed by the PTV method having a resolution of 1 mm in the x-axis and 0.5 mm in the y-axis. For the PTV method, steady state images of two continuous frames were used.

The Stokes number was calculated in order to estimate the ability of the incense smoke to follow the flow field [47], and the value was in the order of 10<sup>5</sup> ; this is an acceptable flow tracing accuracy with error less than 1% [46]. Ozone concentration due to the microplasma was measured by an ozone monitor (Seki Electronics, SOZ-3300). The measurements were carried out by placing the microplasma actuator inside a glass tube which was connected to the ozone monitor at a distance of 90 cm.

The microplasma actuator could be energized by bipolar pulse as shown in Figure 5 with T period and also positive and negative pulses as shown in Figure 6(a) and (b). For both positive and negative pulses at a frequency of 2 kHz, the pulse width was 4 μs. The calculated energy consumption by integrating the discharge power waveform shown in Figure 6 for a single pulse is shown in Figure 7. In order to have microplasma discharge, a voltage higher than 0.5 kV was necessary. The values of energy were almost the same for positive and negative pulses respectively; thus, at 1.4 kV, the calculated value of energy was about 800 μJ for both positive and negative pulses. This corresponds to an average power of 1.6 W.

measure the discharge voltage. For the discharge current measurement, a current probe (Tektronix, P6021) connected to the oscilloscope was used. By multiplying the values of the voltage and current waveform measured with an oscilloscope, the discharge power was obtained for the pulse voltage. Furthermore, in order to calculate the discharge energy, the obtained discharge power waveform was integrated. Particle tracking velocimetry (PTV) was used as the method to measure the induced air flow by the microplasma actuator [45]. The microplasma actuator was placed inside an acrylic box and incense smoke with particles 0.3 μm in diameter [46] was inserted inside the box. For the measurement of the incense smoke particles, a laser particle counter (Kanomax, Model 3886) was used. For the flow measurements, a laser was used (Nd:YVO4 532-nm) to irradiate the particles. The flow of tracer particles was recorded by a high-speed camera (Red Lake, Motionscope M3) having

Figure 3. Schematic of pulsed high voltage source with semiconductor switches.

Figure 2. Geometry of the microplasma actuator. The four high voltage energized electrodes can be driven independently.

(a) Top-view. (b) Cross-sectional view.

6 Actuators

Energy consumption was measured for a discharge voltage of 1.4 kV at various pulse frequencies for both polarities as shown in Figure 8. At 1 Hz, the energy consumption was about 650 μJ. With the increase of frequency, the single pulse energy consumption increased. This could be explained by the charges in plasma that due to the electric field were moved and trapped at the surface of the dielectric. These surface charges [48, 49], when the next opposite polarity is applied, will increase the electric field strength; furthermore, higher energy consumption for a single pulse was obtained with the increase of frequency. Due to recombination or diffusion, the surface charge will decrease in time.

Figure 5. Waveform of pulsed high voltage applied to the microplasma actuator.

Figure 6. Waveforms of discharge voltage and corresponding discharge current and power. (a) Positive pulse. (b) Negative pulse.

Ozone (O3) is harmful to humans due to its high oxidation ability, which is next to fluorine. The World Health Organization Air Quality Guidelines for Europe in its second edition imposed a value of 0.05 ppm as an 8-hours daily average concentration for ozone. The microplasma actuator generates ozone and the concentration of generated ozone was measured by energizing the actuator while it was rolled and inserted into an 8 mm diameter glass tube. An 1.5 L/min airflow was set inside the tube while the actuator was energized. The microplasma actuator was energized by a bipolar pulse with a frequency of 2 kHz. In Figure 9 the characteristics of ozone concentration versus power consumption is shown. The microplasma actuator has an ozone yield of about 28 g/kWh, which is lower than what other researchers reported [50, 51]. Besides the fact that ozone causes damages to tissues, a higher concentration of ozone can cause the degradation of electrodes [52]. Thus, a low ozone yield is necessary for flow control applications [53].

2.1.2. Flow control with microplasma actuator

The micro-scale plasma actuator, besides the requirement of low discharge voltage, has the advantage of easy integration due to its small size. A microplasma actuator similar to the one

Figure 8. Energy consumption of a single pulse versus pulse frequency for a discharge voltage of 1.4 kV.

Figure 7. Energy consumption of single pulse versus discharge voltage for a pulse frequency of 2000 Hz.

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

http://dx.doi.org/10.5772/intechopen.75802

9

Figure 7. Energy consumption of single pulse versus discharge voltage for a pulse frequency of 2000 Hz.

Figure 8. Energy consumption of a single pulse versus pulse frequency for a discharge voltage of 1.4 kV.

#### 2.1.2. Flow control with microplasma actuator

Ozone (O3) is harmful to humans due to its high oxidation ability, which is next to fluorine. The World Health Organization Air Quality Guidelines for Europe in its second edition imposed a value of 0.05 ppm as an 8-hours daily average concentration for ozone. The microplasma actuator generates ozone and the concentration of generated ozone was measured by energizing the actuator while it was rolled and inserted into an 8 mm diameter glass tube. An 1.5 L/min airflow was set inside the tube while the actuator was energized. The microplasma actuator was energized by a bipolar pulse with a frequency of 2 kHz. In Figure 9 the characteristics of ozone concentration versus power consumption is shown. The microplasma actuator has an ozone yield of about 28 g/kWh, which is lower than what other researchers reported [50, 51]. Besides the fact that ozone causes damages to tissues, a higher concentration of ozone can cause the degradation of electrodes [52]. Thus, a low ozone yield is

Figure 6. Waveforms of discharge voltage and corresponding discharge current and power. (a) Positive pulse. (b)

Figure 5. Waveform of pulsed high voltage applied to the microplasma actuator.

necessary for flow control applications [53].

Negative pulse.

8 Actuators

The micro-scale plasma actuator, besides the requirement of low discharge voltage, has the advantage of easy integration due to its small size. A microplasma actuator similar to the one

The flow was visualized using particle tracking velocimetry (PTV) [45]. Sub-micron size incense smoke was used as tracer particle and a Nd YVO4 532 nm laser was utilized to visualize the flow as shown in Figure 4. The phenomena induced by the microplasma actuator,

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

http://dx.doi.org/10.5772/intechopen.75802

11

In the left part of the actuator (x < 0), the strips which were energized were marked as HV1 while the other strips were at floating potential as shown in Figure 10. The energized electrodes in the right part of the actuator (x > 0) were marked as HV4 while the other strips were at floating potential. The electrodes below the dielectric layer were grounded. The duty ratio of the actuator is defined as D by subtracting from 100% the duty ratio of the voltage applied to

A duty ratio of D = 80% represents a duty ratio of 20% for the HV4 electrodes and a duty ratio

The waveform of the applied voltage with D = 20% is shown in Figure 11. The peak of the AC applied voltage was 1.4 kV at a 20 kHz frequency. The burst frequency for the applied voltage

In the initial stages of the discharge, vortexes appeared near the active electrodes as shown in Figure 12 at 2.5 and 5 ms, respectively. The duty ratio was D = 20%; thus, the vortexes near the HV1 electrodes are rotating counter-clockwise and the vortexes near the HV4 electrodes are rotating clockwise. The flow velocity is higher for the HV1 electrodes since the duty ratio is 20% (HV1 electrodes are energized for a shorter time compared with HV4 electrodes). Gradu-

We have developed a numerical simulation code based on the Suzen model [33, 34]. According

Figure 11. Waveform of applied voltage. The duty ratio is the ratio between the time in which the voltage is on and the

the HV4 electrodes. The duty ratios of HV1 and HV4 have a sum of 100%.

ally, the vortexes join together and leftward flow is established up to 50 ms.

2.2. Numerical simulation study of microplasma actuator

was measured using a high-speed camera.

of 80% for the HV1 electrodes.

to the model, the EHD force is:

sum of the time on plus time off.

was 4 kHz.

Figure 9. Ozone concentration measured while driving the microplasma actuator at 1.5 L/min. air flow.

Figure 10. Microplasma actuator. Schematic image. Top-view showing (a) 20 strip-like electrodes and (b) cross-sectional view.

shown in Figure 2 was developed. The schematic of a DBD microplasma actuator is shown in Figure 10 [32, 41]. It consists of 20 strip-like electrodes with a 200 μm width and 16 μm thickness (top-side electrode) which are placed above a plate-like electrode (bottom-side electrode) with a dielectric layer of 25 μm thicknesses in between. Electrodes were made of copper and the dielectric layer was made of resin film. The gap is in micrometer order; thus, if only 1.4 kV is applied, the value of the obtained electric field is high (~107 –108 V/m). This can be considered a low voltage and thus is easier to control and insulate and furthermore requires smaller sized power supplies. The resin film is a flexible polymer material that makes the microplasma actuator suitable for bending and attaching to various shapes, thus making the actuator suitable for flow control applications.

The flow was visualized using particle tracking velocimetry (PTV) [45]. Sub-micron size incense smoke was used as tracer particle and a Nd YVO4 532 nm laser was utilized to visualize the flow as shown in Figure 4. The phenomena induced by the microplasma actuator, was measured using a high-speed camera.

In the left part of the actuator (x < 0), the strips which were energized were marked as HV1 while the other strips were at floating potential as shown in Figure 10. The energized electrodes in the right part of the actuator (x > 0) were marked as HV4 while the other strips were at floating potential. The electrodes below the dielectric layer were grounded. The duty ratio of the actuator is defined as D by subtracting from 100% the duty ratio of the voltage applied to the HV4 electrodes. The duty ratios of HV1 and HV4 have a sum of 100%.

A duty ratio of D = 80% represents a duty ratio of 20% for the HV4 electrodes and a duty ratio of 80% for the HV1 electrodes.

The waveform of the applied voltage with D = 20% is shown in Figure 11. The peak of the AC applied voltage was 1.4 kV at a 20 kHz frequency. The burst frequency for the applied voltage was 4 kHz.

In the initial stages of the discharge, vortexes appeared near the active electrodes as shown in Figure 12 at 2.5 and 5 ms, respectively. The duty ratio was D = 20%; thus, the vortexes near the HV1 electrodes are rotating counter-clockwise and the vortexes near the HV4 electrodes are rotating clockwise. The flow velocity is higher for the HV1 electrodes since the duty ratio is 20% (HV1 electrodes are energized for a shorter time compared with HV4 electrodes). Gradually, the vortexes join together and leftward flow is established up to 50 ms.

#### 2.2. Numerical simulation study of microplasma actuator

shown in Figure 2 was developed. The schematic of a DBD microplasma actuator is shown in Figure 10 [32, 41]. It consists of 20 strip-like electrodes with a 200 μm width and 16 μm thickness (top-side electrode) which are placed above a plate-like electrode (bottom-side electrode) with a dielectric layer of 25 μm thicknesses in between. Electrodes were made of copper and the dielectric layer was made of resin film. The gap is in micrometer order; thus, if only

Figure 10. Microplasma actuator. Schematic image. Top-view showing (a) 20 strip-like electrodes and (b) cross-sectional

Figure 9. Ozone concentration measured while driving the microplasma actuator at 1.5 L/min. air flow.

considered a low voltage and thus is easier to control and insulate and furthermore requires smaller sized power supplies. The resin film is a flexible polymer material that makes the microplasma actuator suitable for bending and attaching to various shapes, thus making the

–108 V/m). This can be

1.4 kV is applied, the value of the obtained electric field is high (~107

actuator suitable for flow control applications.

view.

10 Actuators

We have developed a numerical simulation code based on the Suzen model [33, 34]. According to the model, the EHD force is:

Figure 11. Waveform of applied voltage. The duty ratio is the ratio between the time in which the voltage is on and the sum of the time on plus time off.

Figure 12. PTV results for four active electrodes (left side) and images of the flow for the entire actuator (right side). D = 20% up to 50 ms and D = 70% for more than 50 ms.

$$
\overrightarrow{f} = \rho\_{\vec{c}} \cdot \overrightarrow{E} \text{ ,} \tag{1}
$$

where V represents the potential. Gauss' law is written as:

and furthermore:

If we consider:

conditions:

We can write Eq. (8) as:

charge density is expressed as [33, 34]:

Two independent equations can be written as:

Furthermore, the body force is calculated by

∇ ε� E � �!

where ε represents the permittivity that was calculated as the product of relative permittivity ε<sup>r</sup> and the permittivity of free space ε0. Using the potential V and the Debye length λD, the

> <sup>¼</sup> �<sup>1</sup> λ2 D

Using Eqs. (1) and (5), the body force is calculated. In the Suzen model, it is assumed that the gas particles are weakly ionized; thus, the potential V can be decoupled in a potential due to

<sup>∇</sup> <sup>ε</sup><sup>r</sup> ð Þ¼ � <sup>∇</sup><sup>φ</sup> �r<sup>c</sup>

<sup>¼</sup> �<sup>1</sup> λ2 D

<sup>∇</sup> <sup>ε</sup><sup>r</sup> � <sup>∇</sup>r<sup>c</sup> ð Þ¼ <sup>r</sup><sup>c</sup>

The permittivity between the dielectric and air was considered as the harmonic mean between dielectric permittivity taken as εrd = 4 and air permittivity εrair = 1; thus, the electric field was conserved [33]. Neumann boundary conditions were used for Eq. (7) as outer boundary

!

ε0

λ2 D

!

rc ε0

the external electric field Ø and a potential due to the net charge density φ [33, 34]:

rc ε0

f ! ¼ rc� E ! ¼ r<sup>c</sup> (3)

http://dx.doi.org/10.5772/intechopen.75802

13

V: (5)

V ¼ ϕ þ φ (6)

<sup>∇</sup> <sup>ε</sup><sup>r</sup> � <sup>∇</sup><sup>ϕ</sup> � � <sup>¼</sup> <sup>0</sup> (7)

: (8)

φ (9)

¼ rcð Þ �∇φ (11)

(10)

∇ð Þ¼� ε � ∇V r<sup>c</sup> (4)

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

where f represents the body force per unit volume, r<sup>c</sup> represents the net charge density, and E represents the intensity of the electric field. In this model, magnetic forces were considered. Furthermore, the electric field:

$$
\overrightarrow{E} = -\nabla V \tag{2}
$$

where V represents the potential. Gauss' law is written as:

$$\nabla \left( \varepsilon \cdot \overrightarrow{E} \right) = \rho\_c \tag{3}$$

and furthermore:

$$
\nabla(\varepsilon \cdot \nabla V) = -\rho\_c \tag{4}
$$

where ε represents the permittivity that was calculated as the product of relative permittivity ε<sup>r</sup> and the permittivity of free space ε0. Using the potential V and the Debye length λD, the charge density is expressed as [33, 34]:

$$\frac{\rho\_c}{\varepsilon\_0} = \left(\frac{-1}{\lambda\_D^2}\right) V. \tag{5}$$

Using Eqs. (1) and (5), the body force is calculated. In the Suzen model, it is assumed that the gas particles are weakly ionized; thus, the potential V can be decoupled in a potential due to the external electric field Ø and a potential due to the net charge density φ [33, 34]:

$$V = \phi + \varphi \tag{6}$$

Two independent equations can be written as:

$$\nabla(\varepsilon\_r \cdot \nabla \phi) = 0 \tag{7}$$

$$\nabla(\varepsilon\_r \cdot \nabla \varphi) = \frac{-\rho\_c}{\varepsilon\_0}.\tag{8}$$

If we consider:

f ! ¼ rc� E !

E !

Furthermore, the electric field:

12 Actuators

D = 20% up to 50 ms and D = 70% for more than 50 ms.

where f represents the body force per unit volume, r<sup>c</sup> represents the net charge density, and E represents the intensity of the electric field. In this model, magnetic forces were considered.

Figure 12. PTV results for four active electrodes (left side) and images of the flow for the entire actuator (right side).

, (1)

¼ �∇V (2)

$$\frac{\rho\_c}{\varepsilon\_0} = \left(\frac{-1}{\lambda\_D^2}\right)\wp \tag{9}$$

We can write Eq. (8) as:

$$\nabla(\varepsilon\_r \cdot \nabla \rho\_c) = \frac{\rho\_c}{\lambda\_D^2} \tag{10}$$

Furthermore, the body force is calculated by

$$
\overrightarrow{f} = \rho\_c \cdot \overrightarrow{E} = \rho\_c(-\nabla \phi) \tag{11}
$$

The permittivity between the dielectric and air was considered as the harmonic mean between dielectric permittivity taken as εrd = 4 and air permittivity εrair = 1; thus, the electric field was conserved [33]. Neumann boundary conditions were used for Eq. (7) as outer boundary conditions:

$$\frac{\partial \varphi}{\partial n} = 0 \tag{12}$$

The microplasma actuator shown in the experimental part had 20 strip-like exposed electrodes. These conditions while keeping a high mesh density necessary for accurate results will also increase the simulation time; thus, we chose for the numerical simulation, 12 exposed electrodes and 6 covered electrodes as shown in Figure 13. The energized electrodes were labeled HV1 (3 electrodes) and HV4 (3 electrodes). In the experimental part, the electrodes were energized at an AC waveform having an amplitude of 1.4 kV and a frequency of 20 kHz. The signal was modulated at 4 kHz with duty ratios of 20, 30, 70, and 80%. Considering that the effective value of 1.4 kV peak is 1 kV, we chose for computational reasons for the simulation that the applied voltage is a positive pulse signal with a peak value of 1 kV. Duty ratios for the

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

http://dx.doi.org/10.5772/intechopen.75802

15

Eqs. (7), (10), and (11) were solved to obtain the potential of the external electric field Ø, the potential of the net charge density φ, and the body force. These were calculated before solving the Navier-Stokes, as shown in Figure 15. Near the active electrodes HV1 and HV4, high values of the body force were calculated. For developing the numerical simulation code, the Julia programming language was used [54]. The same duty ratios were used as in the experimental part: first case when electrodes HV1 were energized at 20% duty ratio and HV4 electrodes were energized at 80%; thus, the actuator duty ratio D = 20%, and second case where HV1 electrodes were energized at 70% duty ratio and HV4 electrodes were energized

positive pulse were 20, 30, 70, and 80% as shown in Figure 14.

Figure 14. Simulation voltage waveforms. Positive pulse with duty ratio of 70, 30, 80, and 20%.

Dirichlet boundary conditions were used for Eq. (10) as the outer boundary conditions:

$$
\rho\_c = 0.\tag{13}
$$

In the numerical simulation code, in the dielectric, the values of u and v were 0. Using Eq. (10), the charge distribution over the encapsulated electrode was calculated considering the covered electrodes as the source charge. A value of rc = 0.00751 C/m<sup>3</sup> was used for the source charge. The same values as Suzen [34] were used for the Debye length with λ<sup>D</sup> = 0.00017 m for air and λ<sup>D</sup> = ∞ for the dielectric. Good agreement with experimental data was obtained for these values. After solving Eq. (11), the value of the body force was inserted in the Navier-Stokes equations:

$$\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = \frac{-1}{\rho} \frac{\partial p}{\partial x} + v \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} \right) + Fx \tag{14}$$

$$\frac{\partial v}{\partial t} + \mu \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} = \frac{-1}{\rho} \frac{\partial p}{\partial y} + \nu \left( \frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} \right) + Fy \tag{15}$$

$$\frac{\partial^2 p}{\partial \mathbf{x}^2} + \frac{\partial^2 p}{\partial y^2} = -\rho \left( \frac{\partial \mu}{\partial \mathbf{x}} \frac{\partial \mu}{\partial \mathbf{x}} + 2 \frac{\partial \mu}{\partial y} \frac{\partial \mathbf{v}}{\partial \mathbf{x}} + \frac{\partial \mathbf{v}}{\partial y} \frac{\partial \mathbf{v}}{\partial y} \right) \tag{16}$$

where u and v represent the components of the flow velocity on the x-axis and y-axis, r is the fluid density, p is the pressure, and υ is the kinematic viscosity. The dynamic viscosity μ is written as:

$$
\mu = \mathfrak{g} \,\, \mathfrak{v} \tag{17}
$$

We have used for the air density, r = 1.177 kg/m<sup>3</sup> and for the kinematic viscosity, <sup>υ</sup> = 1.57 � <sup>10</sup>�<sup>5</sup> m2 /s; thus, the calculated dynamic viscosity <sup>μ</sup> = 1.8 � <sup>10</sup>�<sup>5</sup> kg/ms. The computational geometry is shown in Figure 13. The dimensions of the grid were 14 � 14 mm with 561 � 561 grid points. For solving Eqs. (7) and (10), the convergence parameter was 10�<sup>8</sup> . A projection method in primitive variables on a collocated mesh was used for solving the Navier-Stokes equations.

Figure 13. Computational geometry. Upper electrodes HV1 and HV4 were energized and the encapsulated electrodes were grounded.

The microplasma actuator shown in the experimental part had 20 strip-like exposed electrodes. These conditions while keeping a high mesh density necessary for accurate results will also increase the simulation time; thus, we chose for the numerical simulation, 12 exposed electrodes and 6 covered electrodes as shown in Figure 13. The energized electrodes were labeled HV1 (3 electrodes) and HV4 (3 electrodes). In the experimental part, the electrodes were energized at an AC waveform having an amplitude of 1.4 kV and a frequency of 20 kHz. The signal was modulated at 4 kHz with duty ratios of 20, 30, 70, and 80%. Considering that the effective value of 1.4 kV peak is 1 kV, we chose for computational reasons for the simulation that the applied voltage is a positive pulse signal with a peak value of 1 kV. Duty ratios for the positive pulse were 20, 30, 70, and 80% as shown in Figure 14.

∂φ

Dirichlet boundary conditions were used for Eq. (10) as the outer boundary conditions:

equations:

14 Actuators

written as:

<sup>υ</sup> = 1.57 � <sup>10</sup>�<sup>5</sup> m2

Navier-Stokes equations.

were grounded.

∂u ∂t þ u ∂u ∂x þ v ∂u <sup>∂</sup><sup>y</sup> <sup>¼</sup> �<sup>1</sup> r ∂p ∂x þ υ

∂v ∂t þ u ∂v ∂x þ v ∂v <sup>∂</sup><sup>y</sup> <sup>¼</sup> �<sup>1</sup> r ∂p ∂y þ υ

∂2 p ∂x<sup>2</sup> þ ∂2 p <sup>∂</sup>y<sup>2</sup> ¼ �<sup>r</sup>

In the numerical simulation code, in the dielectric, the values of u and v were 0. Using Eq. (10), the charge distribution over the encapsulated electrode was calculated considering the covered electrodes as the source charge. A value of rc = 0.00751 C/m<sup>3</sup> was used for the source charge. The same values as Suzen [34] were used for the Debye length with λ<sup>D</sup> = 0.00017 m for air and λ<sup>D</sup> = ∞ for the dielectric. Good agreement with experimental data was obtained for these values. After solving Eq. (11), the value of the body force was inserted in the Navier-Stokes

> ∂u ∂x ∂u ∂x þ 2 ∂u ∂y ∂v ∂x þ ∂v ∂y ∂v ∂y

where u and v represent the components of the flow velocity on the x-axis and y-axis, r is the fluid density, p is the pressure, and υ is the kinematic viscosity. The dynamic viscosity μ is

We have used for the air density, r = 1.177 kg/m<sup>3</sup> and for the kinematic viscosity,

computational geometry is shown in Figure 13. The dimensions of the grid were 14 � 14 mm with 561 � 561 grid points. For solving Eqs. (7) and (10), the convergence parameter was 10�<sup>8</sup>

A projection method in primitive variables on a collocated mesh was used for solving the

Figure 13. Computational geometry. Upper electrodes HV1 and HV4 were energized and the encapsulated electrodes

∂2 u ∂x<sup>2</sup> þ

∂2 v ∂x<sup>2</sup> þ

/s; thus, the calculated dynamic viscosity <sup>μ</sup> = 1.8 � <sup>10</sup>�<sup>5</sup> kg/ms. The

∂2 u ∂y<sup>2</sup> 

∂2 v ∂y<sup>2</sup> 

μ ¼ ϱ υ (17)

<sup>∂</sup><sup>n</sup> <sup>¼</sup> <sup>0</sup> (12)

r<sup>c</sup> ¼ 0: (13)

þ Fx (14)

þ Fy (15)

(16)

.

Eqs. (7), (10), and (11) were solved to obtain the potential of the external electric field Ø, the potential of the net charge density φ, and the body force. These were calculated before solving the Navier-Stokes, as shown in Figure 15. Near the active electrodes HV1 and HV4, high values of the body force were calculated. For developing the numerical simulation code, the Julia programming language was used [54]. The same duty ratios were used as in the experimental part: first case when electrodes HV1 were energized at 20% duty ratio and HV4 electrodes were energized at 80%; thus, the actuator duty ratio D = 20%, and second case where HV1 electrodes were energized at 70% duty ratio and HV4 electrodes were energized

Figure 14. Simulation voltage waveforms. Positive pulse with duty ratio of 70, 30, 80, and 20%.

Figure 15. Calculated electric potential, charge density, and body force: Electrodes HV1 and HV4 = 1000 V; covered electrode = 0 V. Higher values of body force were obtained near the active exposed electrodes (shown in HV1 and HV4 enlarged).

observed in the experimental results. The maximum calculated value of the flow of 0.83 m/s was near the active electrodes and the value of the flow in the diagonal direction was 0.45 m/s. The PTV measured value of the diagonal flow was 0.43 m/s, which is in agreement with

Figure 16. Flow: Flow for duty ratio D = 20%: Initial vortexes appeared (2 ms); vortexes increase in size (10 ms); leftward

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

http://dx.doi.org/10.5772/intechopen.75802

17

In Figure 17 are shown the simulation results for D = 70%. The electrodes HV1 are energized for a longer time compared with HV4 electrodes; thus, the flow near the HV1 electrodes is faster. The direction of the diagonal flow obtained up to 50 ms is changing gradually, from leftward to rightward as measured by the PTV method. Vortexes are observed near the HV1 electrodes up to 55 ms. A complete tilt toward the right part is observed at 80 ms. At 120 ms, it

numerical simulations. The duty ratio of the actuator after 50 ms was D = 70%.

can be considered that a steady state was achieved.

flow steady state (50 ms).

at 30%; thus, actuator duty ratio D = 70%. The results for D = 20% are shown in Figure 16. Vortexes appeared in the initial stages of the phenomena up to 5 ms. As was observed also in the experimental results, the vortexes appeared above the covered electrodes and had a counter-clockwise direction near the HV1 electrodes and clockwise direction near the HV4 electrodes. With the lapse of time at t = 10 ms, the vortexes moved up; thus, at 15 ms, it can be observed that the vortexes start to join together and create a flow directed diagonal to the left. At 50 ms, the flow is in a steady state with an extended diagonal flow toward the left as

Dielectric Barrier Discharge Microplasma Actuator for Flow Control http://dx.doi.org/10.5772/intechopen.75802 17

Figure 16. Flow: Flow for duty ratio D = 20%: Initial vortexes appeared (2 ms); vortexes increase in size (10 ms); leftward flow steady state (50 ms).

observed in the experimental results. The maximum calculated value of the flow of 0.83 m/s was near the active electrodes and the value of the flow in the diagonal direction was 0.45 m/s. The PTV measured value of the diagonal flow was 0.43 m/s, which is in agreement with numerical simulations. The duty ratio of the actuator after 50 ms was D = 70%.

at 30%; thus, actuator duty ratio D = 70%. The results for D = 20% are shown in Figure 16. Vortexes appeared in the initial stages of the phenomena up to 5 ms. As was observed also in the experimental results, the vortexes appeared above the covered electrodes and had a counter-clockwise direction near the HV1 electrodes and clockwise direction near the HV4 electrodes. With the lapse of time at t = 10 ms, the vortexes moved up; thus, at 15 ms, it can be observed that the vortexes start to join together and create a flow directed diagonal to the left. At 50 ms, the flow is in a steady state with an extended diagonal flow toward the left as

Figure 15. Calculated electric potential, charge density, and body force: Electrodes HV1 and HV4 = 1000 V; covered electrode = 0 V. Higher values of body force were obtained near the active exposed electrodes (shown in HV1 and HV4

enlarged).

16 Actuators

In Figure 17 are shown the simulation results for D = 70%. The electrodes HV1 are energized for a longer time compared with HV4 electrodes; thus, the flow near the HV1 electrodes is faster. The direction of the diagonal flow obtained up to 50 ms is changing gradually, from leftward to rightward as measured by the PTV method. Vortexes are observed near the HV1 electrodes up to 55 ms. A complete tilt toward the right part is observed at 80 ms. At 120 ms, it can be considered that a steady state was achieved.

3. Conclusions

Acknowledgements

Author details

References

Kazuo Shimizu\* and Marius Blajan

Shizuoka University, Hamamatsu, Japan

A microplasma actuator is a device with no moving parts that has advantages over conventional actuators. It can be energized at voltages of about 1 kV, and thus requires a smaller size power supply and less electrical insulation compared with macro-plasma actuators. Moreover, due to the use of pulse voltage to energize the actuator, the energy consumption is low.

Dielectric Barrier Discharge Microplasma Actuator for Flow Control

http://dx.doi.org/10.5772/intechopen.75802

19

The use of a microplasma actuator for flow control showed the capability of the actuator to induce flow and also to change the direction of the flow. By using devices with field-effect transistor (FET) switches to energize the actuator, the flow direction was changed from diagonal leftward to diagonal rightward when the duty ratio of applied voltage was changed. This could be useful if the actuator were to be attached to small drones. The characteristics of the flow were investigated both experimentally and using numerical simulations. A numerical simulation code was developed based on the Suzen-Huang model, which calculates the body force from the potential of the external electric field and the potential of the charge density of the plasma, and implements the body force in Navier-Stokes equations. The experimental

The authors would like to thank Professor Hitoki Yoneda from the University of Electro-Communications, Tokyo, Professor Damon Chandler, Mr. Yoshinori Mizuno, Mr. Akihiko Ito and Mr. Daisuke Nonaka from Shizuoka University, Hamamatsu, for the fruitful discussions.

results confirmed the validity of the developed code.

\*Address all correspondence to: shimizu@cjr.shizuoka.ac.jp

[1] Gadelhak M. Flow Control. Cambridge: Cambridge University Press; 2000

Meeting Exhibit; January 9–12, 2006. DOI: 10.2514/6.2006-1203

[2] Roth JR, Din X. Optimization of the aerodynamic plasma actuator as an electrohydrodynamic (EHD) electrical device. In: Proceedings of the 44th AIAA Aerospace Sciences

Figure 17. Flow for duty ratio D = 70%: After 50 ms, the leftward diagonal flow is changing gradually to rightward diagonal flow.

The maximum values of the flow of 0.82 m/s were calculated near the HV1 electrodes. For the diagonal flow, the value was 0.43 m/s. The value obtained by PTV method was 0.42 m/s for the diagonal flow. In the case of experimental results, due to the microplasma light emission near the active electrodes, the flow could not be measured properly; thus, the simulation gave us valuable insight into the microplasma actuator phenomena. The leftward diagonal flow obtained at D = 20% has a smaller angle with the horizontal axis compared with the rightward diagonal flow at D = 70%.
