Investigation of a Dielectric Barrier Discharge Plasma Actuator to Control Turbulent Boundary Layer Separation

: A numerical investigation was carried out to explore the effects of a dielectric barrier discharge (DBD) plasma actuator on a three-dimensional incompressible, separated ﬂow. The test article selected for the simulations was the National Aeronautical and Space Administration (NASA) wall-mounted hump model. The simulations were run at a Reynolds number of 936,000, based on hump chord length, and a freestream Mach number of 0.1. Hybrid partially averaged Navier–Stokes/large-eddy simulations (PANS/LES) were completed using CALC-LES, a well-validated computational ﬂuid dynamics (CFD) code, developed by Chalmers University of Technology. The baseline code was modiﬁed to simulate the effects of the actuator, which were modeled as source terms in the momentum equation and were assumed to be steady and constant in the span-wise direction. The numerical simulations were carried out for a baseline (no control) case and ﬁve plasma control cases. To optimize the performance of the actuator, the variation of actuator location and voltage frequency was investigated. For the baseline case, a comparison of time-averaged skin friction, the coefﬁcient of pressure, and velocity proﬁles was made of the available experimental results. The results of the baseline case showed good agreement for a hybrid turbulence model, thus strengthening the solver’s ability to predict a three-dimensional separated ﬂow with reasonable accuracy. The results with the plasma actuator turned on showed improved ﬂow characteristics compared to the baseline simulations by reducing the region of separated ﬂow. The actuator placed just downstream of the separation point at an operational frequency of 5kHz completely eliminated the separated ﬂow for our test conditions.


Introduction
A common definition of boundary layer flow control was given by Flatt [1]: "Boundary layer control includes any mechanism or process through which the boundary layer of a fluid flow is caused to behave differently than it normally would were the flow developing naturally along a smooth straight surface". Due to the increase of computational resources and the development of sensor and actuator technology, the numerical and experimental implementation of advanced flow actuation has accelerated rapidly during the past 15 years [2]. The ability to control unwanted fluid flow occurrences is motivated by the desire to improve the aerodynamic performance of air vehicles. For decades, the aerospace industry has actively developed technologies that can eliminate, or at least mitigate, unwanted boundary layer separation [3][4][5]. In fluid dynamics, separated flow may cause an increase in aerodynamic drag, induce mechanical vibration, and decrease lift. These flow features contribute to poor aerodynamic characteristics and structural failures. Thus, by providing effective control and separation mitigation, we can design optimized high-performing air vehicles.
A + e − = A + + 2e − (1) This ionized gas (plasma) moves towards the covered electrode. The ions produce a body force due to the momentum transfer of the flow close to the wall surface. This causes a suction effect above the uncovered electrode and works similarly to a wall jet. The plasma development is strongest in the area close to the electrodes. During the positive half cycle, the upper electrode functions as the cathode, and the dielectric material covers the lower electrode and works as a pseudo-anode. Because of the dielectric material and asymmetric electrode configuration, the plasma formation is relatively low during the negative half cycle. A detailed description of plasma formation and the induced body force mechanisms was described in [11,12]. There are fundamentally two methodologies to numerically model a DBD plasma actuator:

•
First principles based models • Phenomenological or simplified models First principles based modeling approaches are classified as fluid models, kinetic models, or hybrid models. These models are accurate for capturing plasma physics and near wall physics. There were several publications that investigated the use of plasma actuators based on first principles based models [13][14][15][16][17]. On the other hand, phenomenological models focus on the bulk effect of the plasma on the flowfield. These models are popular due to their simplicity of modeling and implementation. There are several validated phenomenological models that are used to model plasmas. Among the most popular are the Shyy model [18], the Suzen and Huang model [19,20], the Roth model [21], and the lumped-element circuit model [22,23].
Following Werner's work, research was conducted showing the effective use of DBD plasma actuators on airfoils, hump models, bluff bodies, compressors, rearward facing steps, and flat plates. Additional studies, such as Suzen et al. [19], studied on the effect of a DBD plasma actuator on PAK-Bblade flow separation and a low-pressure turbine flow with large flow separation in a quiescent environment [20]. Khoshkoo et al. [24] and Gaitonde et al. [25] studied DBD actuators to control flow separation over a NACA 0015 airfoil at a high angle-of-attack, showing the ability of phenomenological DBD actuator models to control low-speed flow. Yu et al. [26] studied a DBD plasma actuator model with a correction factor to examine the complex flow structure. Mushyam et al. [27] used the marker and cell (MAC) based predictor-corrector method to simulate plasma induced electro-hydro-dynamic (EHD) flow and validated their numerical findings with experimental results. Gang et al. [28] conducted experiments of DBD plasma actuators to control external and internal flow separation over a hump and compressor cascade, respectively. A creative stair-shaped DBD plasma actuator configuration was introduced by Rodrigues et al. [29] to increase induced flow velocity and mechanical efficiency with an acceptable actuator durability. Corke et al. [30] showed that DBD plasma actuators can effectively reduce the drag and increase the lift. Rizzetta et al. [31] studied the effects of DBD plasma based control on turbulent separation on a rear facing ramp. Hasan et al. [32], Lu et al. [33], He et al. [34], and others studied the effects of DBD plasma actuators on the NASA hump model. The Shyy model, proposed by Wei Shyy, was adapted for use in this work due to the reasonable accuracy, low computational cost, and simple implementation [24]. The details of this model are presented in the Computational Methodology Section.

Computational Methodology
The computational methodology of the governing equations, numerical solver, plasma actuator model, and computational domain are presented in this section. CALC-LES was used to compute the three-dimensional governing fluid flow equations, known as the Navier-Stokes equations (NSE). The DBD plasma actuator model was numerically implemented in the NSE as source terms based on the linearized Shyy model.

Governing Fluid Flow Equations: Navier-Stokes Equations
The Navier-Stokes equations typically consist of the conservation mass, momentum, and energy. The current work was solely dedicated to capturing the complex physics of the incompressible turbulent flow separation with and without a plasma actuator; however, we neglected the temperature rise due to the plasma as it was considered to be very low and thus ignored [24,25]. In our computations, only the continuity and the momentum equations were considered, where the DBD plasma actuator was modeled as the source term in the momentum equation. The modified set of non-dimensional governing equations for the unsteady, incompressible, plasma induced flow is given by: Equation (2) represents the continuity equation, and the momentum equations are shown in (3)-(5) in the x, y, and z directions, respectively. The density is ρ; P is the pressure; Re is the Reynolds number; and u, v, w are the velocity components in the x, y, and z directions, respectively. The right-hand side of the momentum equation represents the plasma body force source term in three-dimensional space. Equations (2)-(5) represent the complete set of equations used in this modified form of CALC-LES.

DBD Plasma Actuator Model
The phenomenological Shyy model is considered a weakly ionized, non-thermal plasma. This concept simulates the "bulk" motion of fluid flow due to the induced plasma force and is based on the model developed by Roth et al. [21] and Corke et al. [30]. An illustration of the model is shown in Figure 2. This phenomenological model assumes a linearized electric field [11,18,35,36]. Here, the body force generated from the Lorentzian collisions between ions is modeled as a source term in the momentum equation. At a high frequency, electrons travel with high velocity and produce a short duration plasma in which fluid species dissociate and recombine. Due to the greater mass of ions compared to electrons, the plasma body force is induced by momentum transfer of the ions. This is a key concept in the discharge process of the plasma actuator. The electric field lines, shown in Figure 2, roughly display the formation and strength of the electric field due to an asymmetric electric configuration and cause unidirectional force generation. Additionally, the electric field is strongest in the region near the electrodes and decreases in strength further away from the electrodes. Outside of this region, we assumed there was no plasma formation; hence, no force was generated. The force generated by the plasma can be formulated in a qualitative expression. Shyy assumed that the electric field shown was uniform, linear, and parallel. This linearized electric field distribution is shown in Figure 3. The electric field was assumed to decrease in strength as it moved further away from the source. The electric field, Em varies linearly and is expressed as: where E 0 is the maximum electric field and k 1 and k 2 are constants representing the gradient of the electric field intensity along the x and y directions, respectively. Equation (6) shows that the electric field strength, E, varies with position and decreases with distance from the electrodes. The maximum electric field strength is expressed as: d is the distance between the electrodes. The constants k 1 and k 2 are calculated based on the plasma boundary conditions by introducing the breakdown field strength, E b , at locations A and B from Figure 3. The equations for the electric field components are shown in (8) and (9), and the magnitude of the body force is formulated based on the definition of the electrostatic force components in the Lorentz force, shown in Equations (10) and (11).
ρ c is the charge density, and e c is the electron charge; both are constants. Thus, the force components are only a function of electric field strength. As a result, there is a linear increase of the induced body force magnitude with electric field strength. The function δ is used to describe the effectiveness of the plasma-induced body force. Outside the electric field, δ is zero and is unity inside the electric field.
In our case, E cr is also represented by the breakdown electric field strength, E b . The discharge time, ∆t, of the plasma actuator is very small and is a fraction of the time of the duty cycle. The force is produced for a short amount of time, ∆t, of the first half cycle. The force during the second half cycle is negligible due to little or no plasma formation. Moreover, the force is considered as a constant because of the high-voltage frequency of the discharge. The time-averaged force over a complete cycle is expressed as: In Equations (12) and (13), ϑ is the operating frequency and is the reciprocal of the time period, T, whereas, α is the collision efficiency. F tave,x and F tave,y are the time-averaged force components in the x and ydirections, respectively. For the two-dimensional case studies, the time-averaged forces shown in Equations (12) and (13) are added as a source term to the Navier-Stokes equations. Considering, δ = 1, the source term is represented as: The first term in the source vector, S, is the right-hand side of the continuity equation, whereas the second and third components represent the non-dimensional body force due to the plasma actuator. Equations (12) and 13 are used along with the introduction of D c . This term is the non-dimensional plasma force magnitude and is the ratio of the electrical to inertial forces [25,37]. E x and E y are the normalized electric field strengths.

CFD Solver
CALC-LES was used to simulate both the baseline and plasma actuator cases. CALC-LES is an incompressible in-house CFD solver developed by Dr. Lars Davidson of Chalmers University of Technology. It is a FORTRAN based, single block-structured finite volume CFD solver with large eddy simulation (LES), hybrid large eddy simulation (HLES), and Reynolds-averaged Navier-Stokes (RANS) turbulence modeling capabilities [38]. A detailed review of the CALC-LES structure and methodology was explained in [38][39][40][41]. CAL-LES uses the implicit fractional step technique with a multigrid pressure Poisson solver [42]. For spatial discretization of the convective term in the momentum equation, the central difference scheme was blended with an upwinding scheme. Time marching was achieved through the use of the Crank-Nicholson scheme. A hybrid RAN-LES based partially-averaged Navier-Stokes (PANS) model was used to model fine scale turbulence. PANS uses a partial averaging concept, which applies a filtering operation of the fluctuating scales.

Numerical Setup and Computational Domain
The NASA wall-mounted hump model was designed to assess the ability of turbulence models to predict separation and reattachment in an adverse pressure gradient induced separation. There are several excellent publications describing the numerical validation of this model using different numerical schemes [39,[42][43][44][45][46]. The work presented in this paper was the first time that a DBD plasma actuator model was incorporated into CALC-LES to investigate and control separated flow. Moreover, the available experimental values were mostly based on a nominally two-dimensional flow field analysis. Nevertheless, the computational domain and the geometry used in this study were three-dimensional with a spanwise length of 0.    Figure 4 shows the entire domain geometry; however, the modified version of the grid was used to reduce the computational domain, shown in Figure 5. The grid used for this work was the mandatory grid used for hump model computations and was provided by the Advanced Turbulence Simulation for Aerodynamic Application Challenges (ATAAC) project. The grid was fine in area close to the point of incipient separations. The mesh used for our cases had 312 × 120 × 64 cells with Z max = 0.2. The details of the full geometry and initial flow conditions can be found on the NASA CFDVAL2004 and ATAAC project homepage [47]. The maximum heights of the hump and the channel were h/c = 0.128 and H/c = 0.91, respectively.
Due to the three-dimensional wall boundary effects, the top wall was modeled as a slip wall. Inlet flow conditions were adapted from the results of a two-dimensional RANS SST k − ω simulation of Prof. Strelets [42]. The Neumann boundary condition was used for the outflow (right) wall, and a symmetry boundary condition was used in the spanwise direction. A Reynolds number of 936,000, based on chord length, and the inlet velocity of 34.6 m/s were used. As both the velocity and length were normalized by the reference values, the dynamic viscosity was therefore 1/Re.

Results
The result analyses were conducted for two cases: (1) the baseline case with no plasma actuator and (2) the case with the plasma actuator turned on. For the baseline case simulation, subsonic flow over the NASA wall-mounted hump model was simulated with a Mach number of 0.1 and a Reynolds number of 936,000. This baseline case analysis was mainly carried out to observe the turbulent flowfield, capture the separation on the hump surface, and compare the results to the available experimental data. The flow parameters used for the simulations are tabulated below in Table 1.

Baseline Case
In order to investigate the effects of a plasma actuator on a subsonic, incompressible flowfield numerically, it was necessary to understand the baseline flowfield prior to computing results with the plasma actuator turned on. The computational results were compared to the experimental results obtained by Greenblatt et al. [48,49]. As the three-dimensional domain was used for the computation, the upper wall of the domain was modified to compensate for the effects of the side plates. The baseline model of the hump was also studied by several researchers [39,[42][43][44][45][46]. The time step of the computations was 0.002. Time-averaged variables were calculated after 7500 time steps, and sample data were collected after 7500 time steps. The time steps were chosen based on the Strouhal number of 0.2. Figure 6 shown above highlights the turbulent boundary layer flow separation in terms of velocity streamlines and vectors. Figure 6a,b represents the velocity streamlines and vectors, respectively. The streamline figure helps to observe visually the location of the unsteady flow separation on the leeside of the hump due to the adverse pressure gradient. It was observed from the result that the flow separation region was roughly between the non-dimensional length of 0.65 and 1.2. The velocity vectors in Figure 6b magnified the presence of the flow reversal. The instantaneous u-velocity contour shown in Figure 7 also captured the flow separation.   Figure 9. It is noticeable from both figures that the turbulent length scales were larger downstream of the separation region.   For external flow over an object, the resistive drag force consisted of the pressure drag and the drag due to surface roughness, known as the skin friction. The skin friction coefficient and the pressure coefficient are both non-dimensional parameters that represent the pressure and skin friction drag normalized by dynamic pressure and are defined as, Figures 10 and 11 represent the skin friction coefficient and the pressure coefficient, respectively. In both cases, the numerical result was compared with the available experimental result, recorded by Greenblatt [49]. It was observed that the baseline numerical results using CALC-LES predicted the experimental results with reasonable accuracy, thus proving the solver's ability to compute the flowfield using a hybrid turbulence model. The separation and reattachment of the flow near the wall could be estimated from the skin friction where the C f value was zero. From Figure 10, the separation location was identified at x/c = 0.66 for both the numerical and experimental results; whereas the reattachment location for the numerical result was slightly higher than the experimental. The reattachment location was at x/c = 1.1 from the experimental result and x/c = 1.13 from the numerical result. Figure 11 shows the pressure coefficient along the hump wall.  The pressure coefficient, C p , showed the adverse pressure gradient produced on the hump wall, which resulted in flow separation. The pressure increased until the flow reached its location at x/c = 0.66, followed by a sharp decrease of −C p at that position, which indicated the flow separation at location x/c = 0.66. Reattachment at the hump's trailing edge is also found in the figure. The reason for the discrepancy between the numerical and experimental results in the separated region could be due to the end plate effects in the experimental measurements. This was also supported by the results shown by Duda et al. [50], where the numerical results considering the end-plates gave better agreement with the experimental values. Noteworthy is that the numerical simulations for this work were carried out without the consideration of end plates. Figure 12 shows the mean u-velocity component comparisons at different streamwise stations. The stations were selected locations that were comprised of pre-separation, separation, and post-separation regions, located at x/c = 0.65, x/c = 0.80, x/c = 0.90, x/c = 1.1, x/c = 1.2, and x/c = 1.3. The difference between the experimental (orange dotted line) and numerical results (blue solid line) was minimal and therefore showed good agreement. The results were consistent with similar numerical simulations carried out in [19,34,42,44,51]. However, the usage of high-fidelity models could improve the numerical predictions. Locations x/c = 0.80 and x/c = 0.90 were the locations inside the separation region and could be identified from the velocity boundary layer profile seen in Figure 12. The other locations were unaffected by flow separation, in which the first one, i.e, x/c = 0.65, was the location upstream of the separation and x/c = 1.3 was the downstream location of the separation. The Reynolds shear stresses were compared to the experimental values. Figure 13 shows the mean Reynolds shear stresses comparisons at similar streamwise locations that were shown in Figure 12. The numerical results underpredicted the Reynolds stresses compared to the experimental particle image velocimetry (PIV) results. The most noticeable difference was found at location x/c = 0.65, which showed greater shear. The PIV technique inherently had measurement uncertainties close to the wall as seed particles could not get very close to the wall. This was one potential reason for the difference between the numerical and experimental results shown in Figure 13. The error sources of the PIV technique and relative uncertainties associated with it were mentioned in [52]. Smith et al. discussed Reynolds stress uncertainties due to PIV errors [53]. This inherent error is equally problematic for flows suffering separation. It was noted in [44] that the PIV experimental value just upstream of the flow separation is not sufficiently accurate in the near wall region. Figure 13 represents the mean uv normalized stress components and indicates added shear stress accompanied by the separation.

Plasma Actuator Case
The hump model is a modified design of the Glauert airfoil and is used as a benchmark model for the CFD validation in flow separation at the NASA Langley workshop for synthetic jets and turbulence separation control. The baseline results suggested that the separation began at position x/c = 0.66, and so, the DBD plasma actuator was placed at that location. A total of five test cases were investigated for the plasma actuator simulations, based on the applied voltage frequency and the plasma actuator location. Table 2 below displays the test cases. x/c = 0.66 3 kHz 2 x/c = 0.63 3 kHz 3 x/c = 0.70 3 kHz 4 x/c = 0.66 1.5 kHz 5 x/c = 0.66 5 kHz Section 3.2 only describes the results of Test Case 1, whereas, the remaining test cases are discussed in the subsequent subsections. The length and the height of the plasma actuator were denoted by b and a respectively, where the normalized height a/c = 0.018 and normalized length b/c = 0.024. The plasma actuator location was modeled as a body force term in the momentum equation stated in Equation (14). The parametric values of the body force terms described in Table 3 were based on the original Shyy model.
The non-dimensional plasma body force term, D c , was calculated using the formulation stated in Equation (15). The value of D c was calculated based on model geometry and flow conditions. The plasma induced force pushed the existing fluid through the separation zone to reduce or eliminate the flow separation. The magnitude of this force was directly dependent on the plasma size, applied voltage frequency, and electric field strength. Thus, changing these parameters would cause changes in the magnitude of the force. The parameter values needed for the plasma model are illustrated in Table 3. Table 3. Parametric values of the plasma body force modeling.

Parameters Values
Charge density, ρ 1× 10 11 /cm 3 Electron charge, e 1.6 × 10 −19 C Discharge time, ∆t 67 µs Plasma height, a/c 0.018 Plasma length, b/c 0.024 Applied voltage frequency, ϑ 3 kHz Breakdown electric field strength, E b 30 kV/Cm Non-dimensional plasma force parameter, D c 110 Figure 14 shows the comparison between the baseline and plasma cases using velocity streamlines. The plasma induced streamlines are at the left side of Figure 14, whereas the right side of the figure represents the baseline result also shown in Figure 6a. The side-by-side representation of the streamline results indicated the changes in the flow field due to the plasma actuator. There was a substantial reduction in flow separation over the hump surface due to plasma induced force. Nevertheless, the quantitative description of flow separation and flow reattachment was understood from the coefficient of skin friction and the coefficient of pressure displayed in Figures 15 and 16, respectively. The skin friction coefficient gave a precise location of the flow separation and flow reattachment. At the time of flow separation, the C f value decreased to zero, and it rose to zero when the flow reattached. Figure 15 shows the skin friction coefficient in the plasma actuator case and compares the result to the baseline case. The blue line represents the plasma actuator case, and the red line represents the baseline case. A change in the skin friction coefficient is shown in the figure. There was a small delay in separation at the leeside of the hump. The separation started at location x/c = 0.69. The biggest difference was found where the flow reattached. Flow reattachment occurred at location x/c = 0.9. The overall separation region for the plasma enabled flow was from x/c = 0.68 to x/c = 0.90, and the separation region was between location x/c = 0.66 and x/c = 1.13 for the baseline case. Thus, the overall flow separation was reduced by a reasonable margin. The horizontal dotted line is placed at C f equal to zero to determine the streamwise locations of separation and reattachment.   Figure 16 represents the change in pressure coefficient, C p , due to plasma actuator effects, and it compares the plasma result to the baseline result. The pressure dropped by a significant margin at the position where the separation started, i.e., x/c = 0.66. The sudden decrease in C p indicated that the induced plasma body force increased the overall flow velocity compared to the baseline case. This flow remained separated and then reattached downstream of the trailing edge. The plasma force accelerated the fluid near the plasma affected area, which caused a local pressure drop. A significant increase in −Cp was observed at the same location, where a decrease in −Cp was seen for the baseline case. This was due to the pressure increasing suddenly at the separation point and caused the fluid flow near the wall to retard.
The idea behind the plasma actuator was to add momentum to this retarded flow through a plasma induced body force. This was achieved and displayed in Figure 16, where the pressure decreased rapidly to verify the plasma actuator's effectiveness. The coefficient of pressure and the coefficient of skin friction plot shown above sufficiently shows the change in flow distribution, the reduction in flow separation, and the existence of added momentum. As the total mass was conserved, the increase in velocity along the y-axis would provide detailed information of the plasma actuator effect on the velocity field. Figure 17 shows the u-velocity profile along the y-axis at selected axial locations to compare to the baseline values.  , and x/c = 1.1, the flow pattern in the boundary followed in the reverse direction for the separation and recirculation, which is also described in Figure 12. As the plasma actuator was located at the location where the separation began, the actuator assisted the flow to remain attached. Location x/c = 0.80 was still under the separation region, whereas location x/c = 0.90 was free from separation. Figure 17 also holds up the same conclusion with the skin friction coefficient results in Figure 15 regarding the flow reattachment at x/c = 0.90. The further the location was from the flow reattachment, due to sudden flow acceleration, the velocity increased rapidly at the boundary, but slowed down due to relative low momentum flow at the top. This caused the velocity to rise to reach the freestream condition. Figures 18 and 19 show the effect of the plasma induced force by comparing the Q-criterion and the vorticity magnitude, respectively. The positive value of the Q-criterion was the proof of the dominance of rotational components over the strain or stretching components. A push in the flow downstream is visible in Figure 18, and it visually demonstrates the effects of the plasma enabled body force. The blue contours in Figure 19 indicate the reversed flow. It is seen that, for the plasma case, flow reversal, and hence the separation, was significantly smaller than the separation in the baseline case.

Plasma Actuator Location
The results based on the plasma actuator showed that the flow separation and the overall flowfield were changed due to plasma based control. Moreover, a reduction in flow separation was observed for the plasma located at x/c = 0.66. To test the plasma effect without changing its electric field strength or applied voltage frequency, additional analyses were completed by changing the plasma actuator location. The effect of changing the plasma actuator location is described in Figure 21 using the skin friction coefficient. Two specific locations were chosen: (1) upstream of the baseline separation point, i.e., x/c = 0.63, and (2) downstream of the baseline separation point, i.e., x/c = 0.70.
In Figure 21, the skin friction coefficient for three different plasma locations is compared. For the skin friction in both cases with the plasma actuator located at x/c = 0.63 and x/c = 0.66, the line clearly crossed C f = 0, thus demonstrating flow separation and flow reattachment. However, the line for x/c = 0.7 was different than those two positions, so it showed better flow quality in terms of flow separation. This was because the baseline flow separation was substantial at that location, and the plasma actuator positioned at x/c = 0.66 and x/c = 0.63 failed to eliminate the flow separation completely. However, an improvement in reducing the separated flow was evident when the plasma actuator was located downstream of the baseline flow separation location.

Applied Voltage Frequency
The applied voltage frequency used for the DBD plasma actuator was 3 kHz. As at this frequency, the DBD actuator could vary from 1-10 kHz [2], the analysis was done to understand the significance of the applied voltage frequency in the plasma body force. It should be noted that the plasma body force was proportional to the applied voltage frequency, and hence, an increase in frequency would increase the body force. However, this increased induced body force may not guarantee improvement in flow separation. To observe the effects and flow field conditions, skin friction coefficient analysis was done as displayed in Figure 22. Figure 22 indicates that the high frequency created no separation as the yellow line representing 5000 Hz does not cross the C f = 0 line. From the current analysis based on applied voltage frequency variation, it was concluded that the DBD plasma actuator with 5kHz was more effective for flow separation control on the hump model and merited the need for future DBD plasma actuator simulations. As stated earlier, a total of five test cases were run for the plasma actuator simulations, shown in Table 2. The summary of these five cases' results is shown below in Table 4. Test Cases 3 and 5 gave us improved performances for the DBD plasma actuator.

Conclusions
A hybrid RANS-LES (PANS) turbulence model was used to compute the flow over a NASA wall-mounted hump model with and without a DBD plasma actuator. A phenomenological, linearized Shyy model was used over first principles based plasma actuator models for simplicity and accuracy. Comparison of the pressure coefficient, skin friction coefficient, and velocity profiles was made, and good agreement between the numerical and experimental results was found. This was consistent with the results obtained by similar numerical investigations. Hence, it validated the solver's (CALC-LES) ability to reasonably predict flow separation and reattachment in an incompressible flowfield. This work was the first attempt to model a DBD plasma actuator as the source terms in CALC-LES. Here, the baseline code was modified to simulate "bulk" plasma actuator effects. Flow structure comparisons of the baseline case and plasma cases showed that the plasma reduced the separated flow region. The plasma actuator effectively "pushed" the separated fluid downstream and eventually reduced the overall flow separation. Active aerodynamic flow control techniques are still maturing, and optimization (e.g., power, location, frequency) of plasma based flow control devices is necessary for future aerospace applications. The goal of this research was to provide a framework for the implementation of phenomenological plasma based actuators to mitigate the adverse effects of separated flow.

Limitations
The authors did not perform a time independence test for this case study. The time step selected for the simulations was based a Strouhal number of 0.2. The grid used in the work was the mandatory grid from the ATAAC project for the hump model.

Acknowledgments:
We would like to convey our special thanks to Lars Davidson of the Chalmers University of Technology for allowing us to use his solver CALC-LES, which is capable of turbulence LES and RANS-LES hybrid modeling.

Conflicts of Interest:
The authors declare no conflict of interest.

C l
Coefficient of lift C f Coefficient of skin friction C p Coefficient of pressure Re Reynolds number Sr Strouhal number Pr Prandtl number M Mach number