Aerodynamic Investigation of Datum and Slotted Blade Proﬁles under Di ﬀ erent Mach Number Conditions

: Mach number e ﬀ ects on loss and loading are evaluated in both the datum and slotted compressor proﬁles under a wide range of incidences based on two-dimensional (2D) computational ﬂuid dynamic (CFD) simulations. First, total pressure loss and loading abilities are compared. Then, three kinds of deﬁcit thickness are deﬁned and evaluated, and a correlation is made between the loading and the momentum deﬁcit thickness at the proﬁle trailing edge. Finally, the nondimensionalized destruction of mean mechanical energy and dissipation function are employed to analyze the loss mechanism. The slotted proﬁle broadens the low loss range towards the positive incidence range. The slotted proﬁle allows a higher di ﬀ usion factor (DF) than the datum proﬁle. It is hard to distinguish failure simply based on the DF values, whereas the Zweifel loading coe ﬃ cient connects well with the low momentum deﬁcit in the proﬁle trailing edge. The peak of the V-shaped distributions in the Ψ - θ def plot could better suggest the design condition and determine the correct operating range despite the occurrence of bulk separation. The slotted proﬁle gains the ability of the boundary layer ﬂow near the suction surface to resist the adverse pressure gradient, hence, a reduced shear thickness and a uniformed downstream ﬂow ﬁeld is obtained.


Introduction
Air transport will contribute 1.5 billion tonnes of CO 2 emissions annually by 2025 [1]. Hence, it is critical to study aerodynamics with the view to creating more energy-efficient aero-engine. However, the development of more energy-efficient aero-engine requires less cost, lighter weight, and fewer emissions. As the key component of aero-engine, the compressor system play a vital role to the performance. A nominal 1% drop in compressor efficiency (if first stage) would result in about three percent increase in the specific fuel consumption (SFC), due to the changes in both pressure ratio and efficiency [2]. For the development of higher loading, higher efficiency compressors, improving the thrust-to-weight ratio is a permanent pursuit. Usually, increasing the total pressure ratio per stage of compressors is required, resulting in high blade loading. The loss may significantly increase with higher loading.
Designs for high loading are considerably complex but of significant importance. During recent decades, many attempts have been made to obtain high loading for compressors. Wennerstrom [3] made an earlier survey on the approaches taken to achieve a high loading design in axial-flow compressors up to the 1970s. To increase diffusion limits, efforts were associated with some form of boundary layer control. Generally passive approaches were favoured such as vortex generators [4], slots [5][6][7] and performance. The effective operating range is reduced with an increase in Mach number. Hence, analyses of high loading blade profiles with a variation of Mach number play an important role in modern compressor design.
The objective of this paper is to evaluate the effects of Mach number on the loss effect and loading capacities for the datum and slotted blade profiles and explore the underlying loss mechanism. Because the profile loss refers to the loss generated in the blade boundary layer that is located relatively far away from the endwall region, the profile loss could be evaluated in a two-dimensional (2D) way. Performance for the datum profile and the slotted profile with double slots are elaborately compared based on experimentally validated numerical simulations. Firstly, the Mach number effects on the performance and the deficit thickness are evaluated. Then, associations are made for the deficit thickness to the loss and loading capacities. Sources of loss and the physical mechanisms responsible are explored.

Blade Profile
Numerical investigations are conducted on the mid-height section of a highly loaded stator in a modern compressor. In our previous research, a linear compressor cascade using the same blade profile has been numerically and experimentally studied on corner separation control for both the datum blade and the blade with blade end slots [29][30][31]. Table 1 lists the geometrical parameters of the datum profile. The high turning that approaches 53 degrees is representative of the high loading design. The objective of the previous studies aimed to validate the effect of the blade end slots to suppress the three-dimensional (3D) corner separation [29,30]. The subsequent work compared the 3D blading with blade end slots and full-span slots [31], the purpose of which also lay on the control of the corner separation flow. The current research aims at investigating the Mach number effects on the 2D profile properties; thus, the performance for both datum and slotted profiles corresponding to the previous study are evaluated. The slotted profile is cut off from the slotted part of the end-slotted blade that was experimentally tested previously [29,30]. Figure 1 illustrates the geometry for the cross section of the end-slotted blade. The dashed line depicts the datum blade, and the regions covered with section line represent the slotted profile. More details for the slotting design are available in reference [30]. The design philosophy of the slot is quite different from the prior works. The outlet of the slot is placed prior to the separation (considering the endwall flow). The slot is designed to be contracted with its throat located approaching the slot outlet. Small ratios of the slot throat width to the slot exit curvature that promote the Coanda effect are chosen to enhance the attachment of the flow at the slot exit region. Additionally, the employment of befitting axial overlap is also favoured for better guiding of the slot-out-jet. Although the slotted profile is designed for the blade end region, it is still worth evaluating its profile loss and determining the underlying mechanisms.

Computational Details
Steady CFD simulations are conducted using Fluent software, which has been experimentally validated in our previous studies [31,38]. The baseline blade investigated is cut off from a middle stage stator of the compressor rear stage, the design Mach number of which is 0.5884. To investigate Mach number effects, Mach numbers of 0.4, 0.5, 0.5884, 0.7 and 0.8, with a wide incidence range covering −10° to 8°, are considered. The SST turbulence model is known as a relatively robust and accurate two-equation eddyviscosity turbulence model to predict aeronautical flows [39]. To investigate the 3D corner separation in a linear cascade with this same blade profile, the SST model has been modified using helicity to account for energy backscatter in our previous study [38]. Validations shows that the modified SST-Helicity model could significantly improve predicting accuracy for the core separation region [31,38], and it also demonstrated that the model could predict 2D flow (such as mid-span section) well. The SST-Helicity model degenerate into SST model when predicting the 2D flow, thus the SST model is used in the current research to predict the 2D flow field of both datum and slotted profile.
In the simulation, single blade passage is employed due to the periodicity, and the translational periodic boundary condition is used in the pitchwise direction. The variation in Mach number is achieved by automatching the target mass flow rate; thus, 3% in span height is chosen with the translational periodic boundary condition used in the spanwise to make the span height infinite. The structured multi-block O4H grid is used for both datum and slotted profiles, and modifications are only made in the O-mesh region to make the simulation more comparative. Three individual O-mesh blocks are inserted into the original O-mesh zone for the three aerofoil elements in the slotted part. The distribution of the 2D grid section remains the same as for the investigation in reference [31]. The final grid counts is approximately 197 thousand for the datum profile, and 298 thousand for the slotted profile. For the grid of the slotted profile, refinements are made in both the slotted region and the region downstream of the slot outlet.
The finite volume technique is used, and both convective and diffusive fluxes, are discretized by the second-order upwind scheme. The turbulence intensity is specified to be the same as the experiments with a value of 2%. An ideal gas is assumed, and the molecular viscosity is calculated based on Sutherland's viscosity law.

Total Pressure Loss
There are many approaches to estimate loss for turbomachinery flow. The total pressure loss coefficient is the most commonly used one. It is defined as

Computational Details
Steady CFD simulations are conducted using Fluent software, which has been experimentally validated in our previous studies [31,38]. The baseline blade investigated is cut off from a middle stage stator of the compressor rear stage, the design Mach number of which is 0.5884. To investigate Mach number effects, Mach numbers of 0.4, 0.5, 0.5884, 0.7 and 0.8, with a wide incidence range covering −10 • to 8 • , are considered. The SST turbulence model is known as a relatively robust and accurate two-equation eddy-viscosity turbulence model to predict aeronautical flows [39]. To investigate the 3D corner separation in a linear cascade with this same blade profile, the SST model has been modified using helicity to account for energy backscatter in our previous study [38]. Validations shows that the modified SST-Helicity model could significantly improve predicting accuracy for the core separation region [31,38], and it also demonstrated that the model could predict 2D flow (such as mid-span section) well. The SST-Helicity model degenerate into SST model when predicting the 2D flow, thus the SST model is used in the current research to predict the 2D flow field of both datum and slotted profile.
In the simulation, single blade passage is employed due to the periodicity, and the translational periodic boundary condition is used in the pitchwise direction. The variation in Mach number is achieved by automatching the target mass flow rate; thus, 3%c in span height is chosen with the translational periodic boundary condition used in the spanwise to make the span height infinite. The structured multi-block O4H grid is used for both datum and slotted profiles, and modifications are only made in the O-mesh region to make the simulation more comparative. Three individual O-mesh blocks are inserted into the original O-mesh zone for the three aerofoil elements in the slotted part. The distribution of the 2D grid section remains the same as for the investigation in reference [31]. The final grid counts is approximately 197 thousand for the datum profile, and 298 thousand for the slotted profile. For the grid of the slotted profile, refinements are made in both the slotted region and the region downstream of the slot outlet.
The finite volume technique is used, and both convective and diffusive fluxes, are discretized by the second-order upwind scheme. The turbulence intensity is specified to be the same as the experiments with a value of 2%. An ideal gas is assumed, and the molecular viscosity is calculated based on Sutherland's viscosity law.

Total Pressure Loss
There are many approaches to estimate loss for turbomachinery flow. The total pressure loss coefficient is the most commonly used one. It is defined as where the subscript 0 represents the stagnation pressure, and the subscript 1 stands for the inlet value of the computational domain.
To depict the effects of Mach number on loss, Figure 2 contours the mass-averaged total pressure loss coefficient at the position 60%c length downstream of the profile trailing edge. Regions in blue depict that relatively low loss is obtained by the corresponding operating conditions. Regions in red indicate the high-risk operating conditions that should be avoided in the routine design. Figure 2a reflects that the relatively low loss range is wide under low Mach number conditions for the datum profile. With an increase of Mach number the range for the blue regions is reduced sharply. The high-loss region (filled in red) occupies a small negative incidence under Mach numbers above 0.7. With an increase in Mach number, the total pressure loss coefficient increases more rapidly under negative incidence angles compared to that under positive incidence angles. Compared with Figure 2a, Figure 2b reflects that the employment of slots alters the low loss operating range. Clearly, the relatively low loss range is broadened in all of the Mach number conditions tested, and the low loss region moves towards the positive incidence range. Low loss under positive incidence conditions contributes to achieving better aerodynamic performance. The lowest loss is obtained ranging from the incidence angle of 0 • and 2 • , covering Mach numbers of 0.5 to 0.7. Contrary to the datum profile, the increased Mach number makes the total pressure loss coefficient increase more progressively in negative incidence angles than in positive incidence angles for the slotted profile. The high-risk region (filled in red) nearly remains, whereas the white region (with a moderate level of loss) is significantly reduced, particularly under positive incidence angle conditions. Energies 2020, 2, x FOR PEER REVIEW 5 of 22 where the subscript 0 represents the stagnation pressure, and the subscript 1 stands for the inlet value of the computational domain.
To depict the effects of Mach number on loss, Figure 2 contours the mass-averaged total pressure loss coefficient at the position 60% length downstream of the profile trailing edge. Regions in blue depict that relatively low loss is obtained by the corresponding operating conditions. Regions in red indicate the high-risk operating conditions that should be avoided in the routine design. Figure 2a reflects that the relatively low loss range is wide under low Mach number conditions for the datum profile. With an increase of Mach number the range for the blue regions is reduced sharply. The highloss region (filled in red) occupies a small negative incidence under Mach numbers above 0.7. With an increase in Mach number, the total pressure loss coefficient increases more rapidly under negative incidence angles compared to that under positive incidence angles. Compared with Figure 2a, Figure  2b reflects that the employment of slots alters the low loss operating range. Clearly, the relatively low loss range is broadened in all of the Mach number conditions tested, and the low loss region moves towards the positive incidence range. Low loss under positive incidence conditions contributes to achieving better aerodynamic performance. The lowest loss is obtained ranging from the incidence angle of 0° and 2°, covering Mach numbers of 0.5 to 0.7. Contrary to the datum profile, the increased Mach number makes the total pressure loss coefficient increase more progressively in negative incidence angles than in positive incidence angles for the slotted profile. The high-risk region (filled in red) nearly remains, whereas the white region (with a moderate level of loss) is significantly reduced, particularly under positive incidence angle conditions.

Loading Abilities
The de Haller number and local diffusion factor ( ) are the common diffusion parameters used in the compressor routine design. To obtain the de Haller number, the blade inlet and outlet velocities are required, and the demands the velocity distribution of the blade suction surface. Although both the de Haller number and measures the pressure rise ability, they fail in

Loading Abilities
The de Haller number and local diffusion factor (D loc ) are the common diffusion parameters used in the compressor routine design. To obtain the de Haller number, the blade inlet and outlet velocities are required, and the D loc demands the velocity distribution of the blade suction surface. Although both the de Haller number and D loc measures the pressure rise ability, they fail in adequately exposing the aerodynamic difficulties of achieving it. Lieblein [40] developed the most commonly used measurement for loading, termed diffusion factor (DF), in 1953, by evaluating the adverse pressure gradient experienced by the blade suction surface boundary layer. The DF is correlated with D loc , and straightforwardly associated with velocity triangles and solidity. All these loading parameters read as V 1 and V 2 are the average velocities into and out of a blade row, ∆V θ is the change in whirl velocity in the row, and σ is the solidity of the profile row.
All these parameters depict the diffusion achieved by the aerofoil. However, aerodynamic failure could arise under high DF conditions. D loc more directly reflects the real working capacity of the blade, although one drawback for resolving D loc is that the velocity on the wall is zero under the no-slip wall boundary conditions with the viscosity considered. The static pressure distribution can be employed as a replacement, whereas, the 'suction peak' may disappear under large positive incidence angles. Thus, another quantity, the Zweifel loading coefficient, is employed to evaluate the real loading that the blade experienced. The Zweifel loading coefficient is defined as p ps stands for the pressure along the profile pressure surface, and p ss stands for the pressure along the profile suction surface. p 01 is the total pressure into the blade row, and p 1 is the static pressure into the blade row.
Here, the DF is employed to evaluate the diffusion capacity of the configuration, and the Zweifel loading coefficient is used to measure the real loading that suffered by the profile. Figure 3 compares the DF and Zweifel loading coefficient for the datum and slotted profiles varying with incidence angles under different Mach number conditions. The high-risk operating conditions are excluded. Theoretically, the loading that was experienced by the profile is directly associated with the diffusion capacity; hence, a linear relationship is presumed between the DF and Zweifel loading coefficient. However, there appears to be a unimodal behaviour between DF and Zweifel loading coefficient under the same Mach number conditions. The Zweifel loading coefficient peaks are obtained at the incidence angle of -4 • for the datum profile when the Mach number is less than 0.7; and it improves to −1.69 • for the slotted profile. The slotted profile obtains extensively higher diffusion abilities than the datum profile when the DF is over 0.5, and the 'aerodynamic failure' occurs relatively late for the slotted profile. At extremely high positive incidences, the 'aerodynamic failure' under different Mach number conditions shows a tendency to converge.
Obviously, the slotted profile could achieve higher DF than the datum profile when the same Zweifel loading is experienced in the tested operating conditions with the exception of small negative incidence angles under Mach numbers of 0.7. The Zweifel loading coefficient falls steeply for the datum profile (shown in full lines) under positive incidence angles, whereas it reduces gradually for the slotted profile (shown in dashed lines). The highest DF obtained in the extremely high incidence angle is concentrated in the region about 0.65 for the datum profile, and it improves to the range of 0.7 to 0.75 for the slotted profile. The slotted profile acts more effectively under the Mach number of 0.4, 0.5 and 0.5884. When the Mach number increases to 0.7, the loading abilities are reduced for both datum and slotted profiles. The slotted profile shows worse DF or Zweifel loading coefficient than the datum profile in the incidence range below −6 • under Mach numbers of 0.7, whereas it dramatically reverses in the incidence range above −1.69 • . There exists a loading limit for the slotted profile under Mach numbers of 0.7. When the DF for the slotted profile is over 0.75, it does not increase any longer with the further increase of incidence angles. Under Mach numbers of 0.8, a similar loading limit appears for the slotted profile, but it improves to DF = 0.8. The drop of Zweifel loading coefficient under positive incidence angles limits the increasing rate of DF, and the reduction occurs more gradually for the slotted profile. In general, the slotted profile enables a higher DF to be obtained when experiencing the same blade loading.
Energies 2020, 2, x FOR PEER REVIEW 7 of 22 loading coefficient under positive incidence angles limits the increasing rate of DF, and the reduction occurs more gradually for the slotted profile. In general, the slotted profile enables a higher DF to be obtained when experiencing the same blade loading.

Velocity Distribution at the Trailing Edge
In real turbomachinery flows, the limited space as well as the extremely strong shear, high adverse pressure gradient, and unsteadiness, make it difficult to generate a pure boundary layer as the plate flows. Because separation is a key characteristics of compressor flow, the flow mechanism for a compressor blade is far more complex than the conclusions drawn by the plate boundary layer flow investigations. Figure 4 plots the velocity distribution at the profile trailing edge varying with Mach number at some typical incidence angles. The velocity along the line perpendicular to the profile trailing edge is extracted. Figure 4a provides the datum results and Figure 4b gives the slotted results. Velocity distribution on the suction side (SS) is more sensitive to the variation of incidences and incoming flow Mach number. In contrast, the velocity distribution on the pressure side (PS) nearly remains the same shape despite the velocity magnitude for the main flow.

Velocity Distribution at the Trailing Edge
In real turbomachinery flows, the limited space as well as the extremely strong shear, high adverse pressure gradient, and unsteadiness, make it difficult to generate a pure boundary layer as the plate flows. Because separation is a key characteristics of compressor flow, the flow mechanism for a compressor blade is far more complex than the conclusions drawn by the plate boundary layer flow investigations. Figure 4 plots the velocity distribution at the profile trailing edge varying with Mach number at some typical incidence angles. The velocity along the line perpendicular to the profile trailing edge is extracted. Figure 4a provides the datum results and Figure 4b gives the slotted results. Velocity distribution on the suction side (SS) is more sensitive to the variation of incidences and incoming flow Mach number. In contrast, the velocity distribution on the pressure side (PS) nearly remains the same shape despite the velocity magnitude for the main flow.  Seen from Figure 4a, a C-shape velocity distribution is obtained for the datum profile on the SS with the appearance of negative velocity magnitude near the wall bound. For the flow near the wall bound of the profile, a friction force is experienced by the effect of wall shear, whereas the main flow provides the driving force to the flow near the wall bound on the other side. The combination of the driving force from the main flow and the friction force from wall shear makes the boundary layer accumulated more thickly in the downstream development. For the plate, the velocity within the boundary layer is always positive, and the reverse flow cannot occur. In contrast, the existence of the curve makes the boundary layer properties of the blade profile far more complex than the plate flow. The diffusion tunnel determined by the blade profile generates a differential pressure force opposite to the main flow. The wall friction force and the differential pressure force lead to the flow near the wall to decelerate. When the velocity of the boundary layer fluid reduces to near zero, the wall friction force reduces, whereas the differential pressure force remains. The adverse pressure gradient makes the flow near the wall reverse. The reversed flow meets with the downstream flow, then moves away from the boundary wall, leading to the lift off of the boundary layer flow. Thus, the reverse flow at the trailing edge results in negative velocity distribution on the SS. The incidence angle of -4º is extracted as the typical negative incidence. At this angle, a small reverse flow region appears under incoming Mach numbers of 0.4, 0.5, 0.5884, and 0.7, whereas bulk separation occurs when the incoming flow Mach number increases to 0.8. With an increase of incidence angles, the flow trends to more easily separate. The boundary bound of the separated region thickens at the incidence angle of 0°, and the reverse flow region enlarges for the Mach number of 0.8. At the incidence angle of 6°, bulk reverse flow occurs for the datum profile under all Mach numbers tested. Seen from Figure 4a, a C-shape velocity distribution is obtained for the datum profile on the SS with the appearance of negative velocity magnitude near the wall bound. For the flow near the wall bound of the profile, a friction force is experienced by the effect of wall shear, whereas the main flow provides the driving force to the flow near the wall bound on the other side. The combination of the driving force from the main flow and the friction force from wall shear makes the boundary layer accumulated more thickly in the downstream development. For the plate, the velocity within the boundary layer is always positive, and the reverse flow cannot occur. In contrast, the existence of the curve makes the boundary layer properties of the blade profile far more complex than the plate flow. The diffusion tunnel determined by the blade profile generates a differential pressure force opposite to the main flow. The wall friction force and the differential pressure force lead to the flow near the wall to decelerate. When the velocity of the boundary layer fluid reduces to near zero, the wall friction force reduces, whereas the differential pressure force remains. The adverse pressure gradient makes the flow near the wall reverse. The reversed flow meets with the downstream flow, then moves away from the boundary wall, leading to the lift off of the boundary layer flow. Thus, the reverse flow at the trailing edge results in negative velocity distribution on the SS. The incidence angle of −4º is extracted as the typical negative incidence. At this angle, a small reverse flow region appears under incoming Mach numbers of 0.4, 0.5, 0.5884, and 0.7, whereas bulk separation occurs when the incoming flow Mach number increases to 0.8. With an increase of incidence angles, the flow trends to more easily separate. The boundary bound of the separated region thickens at the incidence angle of 0 • , and the reverse flow region enlarges for the Mach number of 0.8. At the incidence angle of 6 • , bulk reverse flow occurs for the datum profile under all Mach numbers tested.
As reflected by Figure 4b, the slotted profile alters the velocity distribution at the trailing edge and produces a more uniform exit flow field. Compared to Figure 4a, the proportion of the main flow region is significantly enlarged by the slotted profile, with lower main flow velocities obtained under relatively larger incidence angles. The slotted profile splits the total turning into multi-part, and the rear elements possess a refreshed boundary layer. This layer contributes to diminishing the accumulated thick boundary layer, and the newly developed boundary layer possesses better ability to resist the adverse pressure gradient. The flow that goes into the slot inlet can be regarded as being similar to the incoming flow, and the slot outlets share approximate static pressure with the local fluid outside the slot outlet region. Thus, the slot-out-jet is self-adaptive. The velocity of the jet is approaching or slightly larger than the local main flow velocity in the blade tunnel near the slot outlet. Additionally, the mass flow rate of the jet flow is also automatically determined by the pressure differential. Because the fluid within the slots suffers less loss than the local fluid in the blade passage near the slot outlet, the self-adaptive jet has a higher momentum than the local low momentum fluid. Figure 4b reflects that the slot-out-jet shares a similar velocity magnitude with the main flow. These self-adaptive jet could reenergize the local low momentum flow, gaining its abilities to resist the adverse pressure gradient, and bringing it downstream, subsequently preventing the accumulation of the low momentum fluid at the blade aft. Compared with the datum profile, the pitchwise extent of the low momentum fluid is extensively reduced. The C-shaped velocity distribution lifts off from the wall due to the jet, and the reverse flow is avoided. The employment of the slots impairs the wall friction that acted on the low momentum boundary flow downstream of the slot outlet and acts as a down flow traction. This traction goes against the differential pressure force in the SS, thus delaying the occurrence of separation.
Additionally, due to the fact that the jet out of the slot is self-adaptive, exaggerative mixing loss would be avoided in comparison to the traditional active jet with fixed mass flow rate.

Deficit Thickness at the Trailing Edge
There are three kinds of quantitative estimations for the plate boundary layer flow, namely the boundary layer thickness, the boundary layer displacement thickness, and the boundary layer momentum thickness. The boundary layer thickness is obtained as the boundary that shares 0.99 times the velocity as the main flow. The boundary layer displacement thickness is defined in terms of volumetric flow rate, and the boundary layer momentum thickness is defined in terms of momentum flux. All these thickness concepts are defined for a thin boundary layer on a plate without curvature. There is no occurrence of reverse flow within or near the boundary layer. For the curved surface, such as the profile investigated in this study, the boundary layer structures are more complex, and the definition of three kinds of boundary layer thickness is no longer suitable for a separated boundary layer with reverse flow under an adverse pressure gradient. To quantify the separation effect within the blade tunnel, three kinds of deficit thickness are defined referencing to the boundary layer thickness for a plate. Similarly, the displacement thickness, δ * de f , is used to evaluate the effect of separation on blockage. This parameter is also used by Gbadebo [41] to estimate the separation extent of the 3D corner separation in a linear compressor cascade. Analogically, the momentum deficit thickness, θ de f , and the kinetic energy deficit thickness, δ de f , are defined, to estimate the separation effect on momentum and kinetic energy, as well as their uniformity. These three kinds of deficit thickness within the blade tunnel take the following form: The curvature and adverse pressure gradient effect makes it difficult to determine the main flow. The velocity distribution at the trailing edge (shown in Figure 4) reflects that the flow is maintained unseparated in the middle region of the blade tunnel. Thus, the velocity in the middle region of the blade tunnel is chosen as the free stream and represented by the subscript f s. The resolution of three kinds of deficit thickness is based on the free stream velocity magnitude obtained along the pitchwise. Considering that the separation in the PS is thin, the blockage and loss is mainly caused by the deficit in the SS. The deficit in the SS are resolved in the profile trailing edge, with y standing for the distance from the wall. The lower deficit thickness represents less deficit and more uniform distribution. Figure 5 compares the displacement thickness at the profile trailing edge. Figure 5 shows that below a Mach number of 0.7 (including 0.7), the displacement thickness for the datum profile is low under a negative incidence angle of −4 • . The displacement thickness increases with the incidence, and rises dramatically at an incidence angle of 6 • , indicating the occurrence of bulk separation. At each incidence angle, the displacement thickness increases with Mach number. Under a Mach number of 0.8, the blockage at the negative incidence angle of −4 • is also increased. Compared with datum profile, the slotted profile effectively decreases the blockage under all the conditions listed. The blockage in the SS is extensively reduced, particularly under a positive incidence angle of 6 • .
Energies 2020, 2, x FOR PEER REVIEW 10 of 22 The curvature and adverse pressure gradient effect makes it difficult to determine the main flow. The velocity distribution at the trailing edge (shown in Figure 4) reflects that the flow is maintained unseparated in the middle region of the blade tunnel. Thus, the velocity in the middle region of the blade tunnel is chosen as the free stream and represented by the subscript . The resolution of three kinds of deficit thickness is based on the free stream velocity magnitude obtained along the pitchwise. Considering that the separation in the PS is thin, the blockage and loss is mainly caused by the deficit in the SS. The deficit in the SS are resolved in the profile trailing edge, with standing for the distance from the wall. The lower deficit thickness represents less deficit and more uniform distribution.
Figure5 compares the displacement thickness at the profile trailing edge. Figure 5 shows that below a Mach number of 0.7 (including 0.7), the displacement thickness for the datum profile is low under a negative incidence angle of −4°. The displacement thickness increases with the incidence, and rises dramatically at an incidence angle of 6°, indicating the occurrence of bulk separation. At each incidence angle, the displacement thickness increases with Mach number. Under a Mach number of 0.8, the blockage at the negative incidence angle of −4° is also increased. Compared with datum profile, the slotted profile effectively decreases the blockage under all the conditions listed. The blockage in the SS is extensively reduced, particularly under a positive incidence angle of 6°. As shown in Figure 6, the momentum deficit thickness shares a similar variation trend to the displacement thickness. At an incidence angle of −4°, the momentum deficit thickness nearly remains for the slotted profile in comparison to the datum profile. With an increase in incidence angle, the slotted profile results in much lower momentum deficit thickness, depicting that the momentum deficit caused by the separation in SS is effectively reduced by the slotted profile, and the flow uniformity in the SS is enhanced. As shown in Figure 6, the momentum deficit thickness shares a similar variation trend to the displacement thickness. At an incidence angle of −4 • , the momentum deficit thickness nearly remains for the slotted profile in comparison to the datum profile. With an increase in incidence angle, the slotted profile results in much lower momentum deficit thickness, depicting that the momentum deficit caused by the separation in SS is effectively reduced by the slotted profile, and the flow uniformity in the SS is enhanced. Figure 7 reflects that the slotted profile slightly increases the kinetic energy deficit thickness at a negative incidence angle of −4 • below a Mach number of 0.7. The positive effect of the slot is shown at 0 • and 6 • incidence angles. The improvement of slotted profile on kinetic deficit thickness at positive incidence angle is much better than the deterioration under negative incidence angle.  Figure 7 reflects that the slotted profile slightly increases the kinetic energy deficit thickness at a negative incidence angle of −4° below a Mach number of 0.7. The positive effect of the slot is shown at 0° and 6° incidence angles. The improvement of slotted profile on kinetic deficit thickness at positive incidence angle is much better than the deterioration under negative incidence angle. Traditionally, it is difficult for a passive flow control method to achieve general applicability in a wide operating range. Most of the time, one optimized scheme functions well under a certain operating condition, whereas it deteriorates under another operating conditions. Contrary to the majority of passive flow control schemes, the proposed slotting scheme has been demonstrated to be self-adaptive regardless of the variation in Mach number and incidence angle. It suppresses the separation and improves the performance at positive incidence, and does not deteriorate much in negative incidence.

DF Varying with Momentum Deficit Thickness
Lieblein [40] demonstrated an almost unique correlation between profile loss and local diffusion through a measurement of the adverse pressure gradient applied to the suction surface boundary layer, known as the diffusion factor (DF). Over nearly seven decades, DF has been widely used in guiding the design of compressors. It is based on the principle that the profile losses are uniquely correlated with the local diffusion ability. In the current research, one question arises regarding  Figure 7 reflects that the slotted profile slightly increases the kinetic energy deficit thickness at a negative incidence angle of −4° below a Mach number of 0.7. The positive effect of the slot is shown at 0° and 6° incidence angles. The improvement of slotted profile on kinetic deficit thickness at positive incidence angle is much better than the deterioration under negative incidence angle. Traditionally, it is difficult for a passive flow control method to achieve general applicability in a wide operating range. Most of the time, one optimized scheme functions well under a certain operating condition, whereas it deteriorates under another operating conditions. Contrary to the majority of passive flow control schemes, the proposed slotting scheme has been demonstrated to be self-adaptive regardless of the variation in Mach number and incidence angle. It suppresses the separation and improves the performance at positive incidence, and does not deteriorate much in negative incidence.

DF Varying with Momentum Deficit Thickness
Lieblein [40] demonstrated an almost unique correlation between profile loss and local diffusion through a measurement of the adverse pressure gradient applied to the suction surface boundary layer, known as the diffusion factor (DF). Over nearly seven decades, DF has been widely used in guiding the design of compressors. It is based on the principle that the profile losses are uniquely correlated with the local diffusion ability. In the current research, one question arises regarding Traditionally, it is difficult for a passive flow control method to achieve general applicability in a wide operating range. Most of the time, one optimized scheme functions well under a certain operating condition, whereas it deteriorates under another operating conditions. Contrary to the majority of passive flow control schemes, the proposed slotting scheme has been demonstrated to be self-adaptive regardless of the variation in Mach number and incidence angle. It suppresses the separation and improves the performance at positive incidence, and does not deteriorate much in negative incidence.

DF Varying with Momentum Deficit Thickness
Lieblein [40] demonstrated an almost unique correlation between profile loss and local diffusion through a measurement of the adverse pressure gradient applied to the suction surface boundary layer, known as the diffusion factor (DF). Over nearly seven decades, DF has been widely used in guiding the design of compressors. It is based on the principle that the profile losses are uniquely correlated with the local diffusion ability. In the current research, one question arises regarding whether the classical DF is still valid for modern high loading profiles with variation of Mach number. Another question is how the slotted profile impacts the trend of DF. To illustrate these points, the momentum deficit thickness in SS of the profile trailing edge varying with DF is depicted in Figure 8. whether the classical DF is still valid for modern high loading profiles with variation of Mach number. Another question is how the slotted profile impacts the trend of DF. To illustrate these points, the momentum deficit thickness in SS of the profile trailing edge varying with DF is depicted in Figure 8.  Figure 8 reflects that for the datum profile, the momentum deficit thickness increases with incidence angle except for the extremely small negative incidence angles of −10° and −8°. Above the DF value of 0.6, a sharp increase in momentum deficit thickness appears for the datum profile, and the DF keeps increasing slowly regardless of the bulk increase in momentum deficit. The variation of momentum deficit thickness to the DF under different Mach number conditions collapses onto an almost unique curve above the DF value of 0.6 for the datum profile. The DF fails to indicate the rapid increase in momentum deficit. For the slotted profile, the critical DF value, above which the sharp increase in momentum deficit occurs, moves from 0.6 to 0.7. Similar to the datum profile, the DF keeps rising slowly with the sharp increase in momentum deficit thickness. However, it is difficult to distinguish failure simply based on the DF value for a given geometry.
As a simplified empirical parameter, the DF is derived based on the boundary layer theory. However, the classical boundary layer theory is established for the plate boundary layer without streamwise pressure gradient. The separation in the SS near the trailing edge of the curved profile caused by the adverse pressure gradient within the compressors makes the flow complicated. Additionally, the derivation of DF is also based on the assumption that the distribution of the SS velocity must occur for the maximum velocity peak near the leading edge to the trailing edge velocity. However, the maximum velocity peak may disappear under large positive incidence angles. All of these factors make the DF only suitable for the flow around a profile without separation and under small incidence angles. When large reverse flow occurs under a positive incidence angle, the diffusion ability drops, whereas the slowly increased DF could not forecast the aerodynamic failure. It should be noted that the implementation of DF is not arbitrary for all conditions, and the higher DF value does not indicate a higher diffusion ability.
Comprehensively, the DF fails to adequately pick up the variation in momentum deficit of the profile as was its original intention with increasing incidence angle under high loadings. Thus, the large DF is not equal to the high loading that could be achieved. However, before the sharp increase of loss, the DF value could be considered in evaluating the loading, even for the slotted profile.  Figure 8 reflects that for the datum profile, the momentum deficit thickness increases with incidence angle except for the extremely small negative incidence angles of −10 • and −8 • . Above the DF value of 0.6, a sharp increase in momentum deficit thickness appears for the datum profile, and the DF keeps increasing slowly regardless of the bulk increase in momentum deficit. The variation of momentum deficit thickness to the DF under different Mach number conditions collapses onto an almost unique curve above the DF value of 0.6 for the datum profile. The DF fails to indicate the rapid increase in momentum deficit. For the slotted profile, the critical DF value, above which the sharp increase in momentum deficit occurs, moves from 0.6 to 0.7. Similar to the datum profile, the DF keeps rising slowly with the sharp increase in momentum deficit thickness. However, it is difficult to distinguish failure simply based on the DF value for a given geometry.
As a simplified empirical parameter, the DF is derived based on the boundary layer theory. However, the classical boundary layer theory is established for the plate boundary layer without streamwise pressure gradient. The separation in the SS near the trailing edge of the curved profile caused by the adverse pressure gradient within the compressors makes the flow complicated. Additionally, the derivation of DF is also based on the assumption that the distribution of the SS velocity must occur for the maximum velocity peak near the leading edge to the trailing edge velocity. However, the maximum velocity peak may disappear under large positive incidence angles. All of these factors make the DF only suitable for the flow around a profile without separation and under small incidence angles. When large reverse flow occurs under a positive incidence angle, the diffusion ability drops, whereas the slowly increased DF could not forecast the aerodynamic failure. It should be noted that the implementation of DF is not arbitrary for all conditions, and the higher DF value does not indicate a higher diffusion ability.
Comprehensively, the DF fails to adequately pick up the variation in momentum deficit of the profile as was its original intention with increasing incidence angle under high loadings. Thus, the large DF is not equal to the high loading that could be achieved. However, before the sharp increase of loss, the DF value could be considered in evaluating the loading, even for the slotted profile.

Zweifel Loading Coefficient Varying with Momentum Deficit Thickness
Because the DF fails to show a better correlation with the momentum deficit, another loading parameter, the Zweifel loading coefficient, is used to correct the loss to loading. Figure 9 depicts the momentum deficit thickness in SS of the profile trailing edge varying with the Zweifel loading coefficient. To improve clarity, the region that is concentrated with the majority of points is amplified. momentum deficit thickness in SS of the profile trailing edge varying with the Zweifel loading coefficient. To improve clarity, the region that is concentrated with the majority of points is amplified. Seen from Figure 9, a V-shaped distribution is obtained under each Mach number condition. The opening extent of the V-shape for the datum profile is larger than the slotted profile under the same Mach number conditions. A large Zweifel loading coefficient connects well with the low momentum deficit in the profile trailing edge. With an increase of Mach number, the peak of the V-shape moves towards a thicker momentum deficit thickness and lower Zweifel loading both for datum and slotted profiles. Under different Mach number conditions, the branches of the V-shape are nearly parallel to each other, indicating that the increasing rate of the momentum deficit in the SS to the Zweifel loading coefficient of the datum profile is undisturbed by the variation in Mach number. The largest Zweifel loading coefficient does not meet with the lowest momentum deficit thickness, and the peak of the V-shape indicates the operating condition between the lowest momentum deficit and the largest Zweifel loading. Although the drop in Zweifel loading lags behind the increase in momentum deficit, it could correctly forecast the aerodynamic failure in diffusion ability caused by the high momentum deficit. Comprehensively, the peak of the V-shape could better suggest the design condition, under which the profile could achieve better performance both in loss and loading ability. A region that is near the V-shape peak could be chosen to extract the correct operating range. The light grey in the background of Figure 9 contours the region with high performance-the darker the colour, the better the performance. Thus, theplot could assist in finding the design point and determine the preferred operating range for a particular profile. This plot is not influenced by the incidence range and the occurrence of bulk separation under positive incidence angles. The recommended design incidence angle (refer to the peak) for the slotted profile is larger than the datum profile, and the effective operating range suggested by the plot is extensively broadened for the slotted profile. It accords well with the total pressure loss contours shown in Figure 2. Hence, this judgement is suitable for slotted profile.

The Underlying Loss Mechanism
The viscous shear stress dominates in the energy transform and transport process. Equation (9) describes changes in kinetic energy per unit mass for a fluid particle [42]. It depicts the averaged turbulence motion. Seen from Figure 9, a V-shaped distribution is obtained under each Mach number condition. The opening extent of the V-shape for the datum profile is larger than the slotted profile under the same Mach number conditions. A large Zweifel loading coefficient connects well with the low momentum deficit in the profile trailing edge. With an increase of Mach number, the peak of the V-shape moves towards a thicker momentum deficit thickness and lower Zweifel loading both for datum and slotted profiles. Under different Mach number conditions, the branches of the V-shape are nearly parallel to each other, indicating that the increasing rate of the momentum deficit in the SS to the Zweifel loading coefficient of the datum profile is undisturbed by the variation in Mach number.
The largest Zweifel loading coefficient does not meet with the lowest momentum deficit thickness, and the peak of the V-shape indicates the operating condition between the lowest momentum deficit and the largest Zweifel loading. Although the drop in Zweifel loading lags behind the increase in momentum deficit, it could correctly forecast the aerodynamic failure in diffusion ability caused by the high momentum deficit. Comprehensively, the peak of the V-shape could better suggest the design condition, under which the profile could achieve better performance both in loss and loading ability. A region that is near the V-shape peak could be chosen to extract the correct operating range. The light grey in the background of Figure 9 contours the region with high performance-the darker the colour, the better the performance. Thus, the Ψ-θ de f plot could assist in finding the design point and determine the preferred operating range for a particular profile. This plot is not influenced by the incidence range and the occurrence of bulk separation under positive incidence angles. The recommended design incidence angle (refer to the peak) for the slotted profile is larger than the datum profile, and the effective operating range suggested by the plot is extensively broadened for the slotted profile. It accords well with the total pressure loss contours shown in Figure 2. Hence, this judgement is suitable for slotted profile.

The Underlying Loss Mechanism
The viscous shear stress dominates in the energy transform and transport process. Equation (9) describes changes in kinetic energy per unit mass for a fluid particle [42]. It depicts the averaged turbulence motion. d dt Energies 2020, 13, 1673 14 of 22 The first term on the left-hand side is the transport of k, where k represents the mean kinetic energy, and it takes k = 1 2 u 2 i . The first term, ∂ − p ρ u j /∂x j , on the right-hand side expresses the work on the fluid per unit of mass by pressure during unit time, and the second term, ∂ 2νu i S ij /∂x j , represents the transport of mechanical energy to the fluid motion by the viscous force of the averaged flow. 2νS ij S ij is the viscous dissipation for the averaged flow. Another two terms that includes the fluctuation reflect the influence of fluctuation to the averaged flow. The term −∂ u i u j u i /∂x j represents the transport of mechanical energy by Reynolds stress, and the term u i u j S ij is the production of the mean kinetic energy.
The last two terms exhibit the destruction of mechanical energy due to turbulent shear work and the shear work by the mean strain. There are two kinds of force that resists the deformation in fluid motion. One of them is the molecular viscosity, and the other is the Reynolds force caused by turbulence fluctuation. The work done on the viscous force that resists the deformation transforms the mechanical energy into kinetic energy for thermal motion of the molecule. Thus, it is regarded as a dissipation term, and it contributes to the energy transformation between mechanical and internal energies. The work done on the Reynolds stress that resists the mean flow deformation transforms the mean flow mechanical energy into kinetic energy for turbulence fluctuation. Namely it is the production of k. Hence, the nondimensionalized destruction of mean mechanical energy for the Reynolds-averaged flow takes the following form: For the flow resolved by Reynolds-averaged Navier-Stokes (RANS) methods with SST turbulence model, the destruction of the mean mechanical energy that is modelled could be reduced as: where U e is the characteristic velocity of the flowfield. Here U e is converted by the incoming flow Mach number and temperature conditions. ν t is the turbulence kinetic viscosity, and ν is the molecular kinetic viscosity. The second term of the c d also terms as dissipation function, φ d . Hence, the nondimensionalized dissipation function is defined as: The shear work done by the molecular viscosity irreversibly dissipates the mechanical energy of the fluid motion as heat. The dissipation term remains positive, and it acts as a source term rather than a transport term. Additionally, this term is closely associated with entropy production. Here, both the nondimensionalized destruction coefficient and the nondimensionalized dissipation function are employed to analyze the control mechanism. High destruction exists where the shear strain dominates. The attached boundary layer in the SS has high destruction coefficient. For the datum profile, the high destruction coefficient region covers a certain axial chord length in the suction surface, then the high destruction bound lifts off from the suction surface, whereas the high destruction coefficient region nearly remains attached in the suction surface for the slotted profile. Thus, the relatively high destruction region is reduced by the slotted High destruction exists where the shear strain dominates. The attached boundary layer in the SS has high destruction coefficient. For the datum profile, the high destruction coefficient region covers a certain axial chord length in the suction surface, then the high destruction bound lifts off from the suction surface, whereas the high destruction coefficient region nearly remains attached in the suction surface for the slotted profile. Thus, the relatively high destruction region is reduced by the slotted profile, particularly when premature lift off of high destruction bound occurs for the datum profile. A high destruction region appears for the datum profile behind the profile trailing edge, whereas it is diminished for the slotted profile.

Destruction Coefficient
To quantitatively extract the lift off location, and evaluate the Mach number effects on destruction coefficient near the wall, the nondimensionalized destruction coefficient distribution along the profile suction surface is compared in Figure 11. To quantitatively extract the lift off location, and evaluate the Mach number effects on destruction coefficient near the wall, the nondimensionalized destruction coefficient distribution along the profile suction surface is compared in Figure 11.
The nondimensionalized destruction coefficient distribution along the profile suction surface is closely associated with the boundary layer development in the profile SS. Figure 11a shows that a valley appears at approximately a 12% axial chord length for the datum profile regardless of the variation in Mach number. This observation indicates the transformation of a favourable pressure gradient to an adverse pressure gradient. The occurrence of a destruction coefficient that approaches zero and then remains in the suction surface for the datum profile; indicates the lift off of the high destruction bound. The location of the initial lift off point is remarkably influenced by incidence angles. The increasing incidence angle promotes the premature lift off of the high destruction bound. Besides, with an increase of Mach number, the initial lift off point moves upstream along the suction surface.
The premature lift off makes it easy to form bulkier reverse flow region in the profile SS, because the high destruction coefficient bound comes into the main flow region earlier. A peak exists for the suction surface distribution of destruction coefficient downstream of the valley value. Either the increased incidence angle or the increased Mach number contributes to the decreased destruction coefficient peak value. However, the boundary layer with a lower destruction coefficient has a weaker ability to resist the adverse pressure gradient. Thus, the separation occurs more easily for the datum profile under larger incidence angles or larger Mach number conditions. The nondimensionalized destruction coefficient distribution along the profile suction surface is closely associated with the boundary layer development in the profile SS. Figure 11a shows that a valley appears at approximately a 12% axial chord length for the datum profile regardless of the variation in Mach number. This observation indicates the transformation of a favourable pressure gradient to an adverse pressure gradient. The occurrence of a destruction coefficient that approaches zero and then remains in the suction surface for the datum profile; indicates the lift off of the high destruction bound. The location of the initial lift off point is remarkably influenced by incidence angles. The increasing incidence angle promotes the premature lift off of the high destruction bound. Besides, with an increase of Mach number, the initial lift off point moves upstream along the suction surface. The premature lift off makes it easy to form bulkier reverse flow region in the profile SS, because the high destruction coefficient bound comes into the main flow region earlier. A peak exists for the suction surface distribution of destruction coefficient downstream of the valley value. Either the increased incidence angle or the increased Mach number contributes to the decreased destruction coefficient peak value. However, the boundary layer with a lower destruction coefficient has a weaker ability to resist the adverse pressure gradient. Thus, the separation occurs more easily for the datum profile under larger incidence angles or larger Mach number conditions. Figure 11b reflects that the destruction coefficient valley appears slightly downward for the slotted profile in the first element of the profile. Then, the destruction coefficient increases until the end of the first profile element. Another destruction coefficient valley is generated in the second profile element at the exit of the front slot, and a destruction coefficient peak is also generated in the second profile element. These destruction coefficient peaks have larger destruction coefficient values than the peak values for the datum profile under the same operating conditions, indicating that the jet from the slot outlet enhances the attachment of boundary layer flow along the suction surface of the second profile element. The valley does not appear in the last profile element due to the strong adverse pressure gradient, whereas two higher peak values are obtained in the last profile element. One peak is located at the exit of the rear slot, and the other peak appears where a balance is obtained between the slot-out-jet and the local main flow. The attached flow in the suction side of the last profile element is enhanced more effective; thus, the slotted profile gains the ability of the boundary layer flow near the suction surface to resist the adverse pressure gradient, and the lift off of the high destruction coefficient bound is eliminated. Figure 11. Nondimensionalized destruction coefficient distribution along profile suction surface. Figure 11b reflects that the destruction coefficient valley appears slightly downward for the slotted profile in the first element of the profile. Then, the destruction coefficient increases until the end of the first profile element. Another destruction coefficient valley is generated in the second profile element at the exit of the front slot, and a destruction coefficient peak is also generated in the second profile element. These destruction coefficient peaks have larger destruction coefficient values than the peak values for the datum profile under the same operating conditions, indicating that the jet from the slot outlet enhances the attachment of boundary layer flow along the suction surface of the second profile element. The valley does not appear in the last profile element due to the strong adverse pressure gradient, whereas two higher peak values are obtained in the last profile element. One peak is located at the exit of the rear slot, and the other peak appears where a balance is obtained between the slot-out-jet and the local main flow. The attached flow in the suction side of the last profile element is enhanced more effective; thus, the slotted profile gains the ability of the boundary layer flow near the suction surface to resist the adverse pressure gradient, and the lift off of the high destruction coefficient bound is eliminated.

Dissipation Function
The logarithm of the nondimensionalized dissipation function is shown in Figure 12. High dissipation indicates a strong shear region. For the datum profile, the high shear bound that lifted off from the suction surface processes a high magnitude of dissipation, and the wake also processes a high dissipation value. The slotted profile suppresses the lift off of the high shear bound, and also results in a reduced shear thickness. Hence, a more uniformed downstream flow field is obtained by the slotted profile.
To evaluate the influence of the dissipation function, the pitchwise distribution of a nondimensionalized dissipation function in 20%c downstream of the profile trailing edge is compared for the datum and slotted profiles. It should be noted that the scale for the vertical coordinate is different. To evaluate the influence of the dissipation function, the pitchwise distribution of a nondimensionalized dissipation function in 20% downstream of the profile trailing edge is compared for the datum and slotted profiles. It should be noted that the scale for the vertical coordinate is different.
The dissipation in the PS is small for both datum and slotted profiles. As shown in Figure 13a, two dissipation peaks appears for the datum profile in the SS. The higher peak exposes the dissipation caused by the wake, and the other peak is due to the high destruction bound, namely the shear between the separated flow on the aft of the profile SS and the main flow generates high dissipation. The increase in incidence angle facilitates a wider range of high dissipation range. With an increase of Mach number, peak values of dissipation reduce, and the high dissipation range over the pitchwise direction is enlarged. Because the dissipation finally dissipates as heat, it more effectively suppresses high dissipation, thus reducing the loss. Figure 13b reflects that nearly two thirds of dissipation extent is reduced by the slotted profile. Apparently, a narrow wake is obtained with lower loss. The high dissipation peak in the SS for the slotted profile is increased to three in accordance to two for the datum profile. Referring to the dissipation contours shown in Figure 12, the high dissipation peak in The dissipation in the PS is small for both datum and slotted profiles. As shown in Figure 13a, two dissipation peaks appears for the datum profile in the SS. The higher peak exposes the dissipation caused by the wake, and the other peak is due to the high destruction bound, namely the shear between the separated flow on the aft of the profile SS and the main flow generates high dissipation. The increase in incidence angle facilitates a wider range of high dissipation range. With an increase of Mach number, peak values of dissipation reduce, and the high dissipation range over the pitchwise direction is enlarged. Because the dissipation finally dissipates as heat, it more effectively suppresses high dissipation, thus reducing the loss. Figure 13b reflects that nearly two thirds of dissipation extent is reduced by the slotted profile. Apparently, a narrow wake is obtained with lower loss. The high dissipation peak in the SS for the slotted profile is increased to three in accordance to two for the datum profile. Referring to the dissipation contours shown in Figure 12, the high dissipation peak in the detached line is caused by the wake of each profile element. Due to the injection of the slot-out-jet, wakes from the front and second profile elements do not lift off and follow well with the profile suction surface. Dissipation for the wake of the front and second profile element is lower than for the high destruction bound of the datum profile. This is due to the mitigation of a local adverse pressure gradient by the slot-out-jet. Additionally, the slot-out-jet driven by the pressure difference is self-adaptive; thus, high mixing loss is avoided in the jet area. The dissipation in the pressure surface near the trailing edge is impaired by the slotted profile; hence, the dissipation formed in the profile trailing edge is reduced. Under an incidence angle of −4 • , a small extent of separation occurs for the datum profile; thus, the entire relatively high dissipation range for the slotted profile is slightly wider than the datum profile, but the dissipation magnitude for the slotted profile is smaller. Under an incidence angle of 0 • , both the high dissipation range and the dissipation magnitude are reduced through the implementation of slots. The high dissipation range increases to nearly half of the pitch wide for the datum profile under the incidence angle of 6 • , whereas it is reduced to no more than a quarter pitch wide for the slotted profile, and the peak value is also halved.
Energies 2020, 2, x FOR PEER REVIEW 18 of 22 the detached line is caused by the wake of each profile element. Due to the injection of the slot-outjet, wakes from the front and second profile elements do not lift off and follow well with the profile suction surface. Dissipation for the wake of the front and second profile element is lower than for the high destruction bound of the datum profile. This is due to the mitigation of a local adverse pressure gradient by the slot-out-jet. Additionally, the slot-out-jet driven by the pressure difference is selfadaptive; thus, high mixing loss is avoided in the jet area. The dissipation in the pressure surface near the trailing edge is impaired by the slotted profile; hence, the dissipation formed in the profile trailing edge is reduced. Under an incidence angle of −4°, a small extent of separation occurs for the datum profile; thus, the entire relatively high dissipation range for the slotted profile is slightly wider than the datum profile, but the dissipation magnitude for the slotted profile is smaller. Under an incidence angle of 0°, both the high dissipation range and the dissipation magnitude are reduced through the implementation of slots. The high dissipation range increases to nearly half of the pitch wide for the datum profile under the incidence angle of 6°, whereas it is reduced to no more than a quarter pitch wide for the slotted profile, and the peak value is also halved.

Conclusions
The Mach number effects on loss and loading are evaluated both for the datum and slotted profiles under a wide incidence range, and the underlying loss mechanism is explored based on experimentally validated CFD simulations. Conclusions are drawn as follows: • The mass-averaged total pressure loss contours showed that the low loss range for the datum profile reduced sharply with increasing Mach number. The employment of slots alters the low loss operating range towards a positive incidence range and broadens the low loss range in all Mach number conditions tested.

•
The DF and Zweifel loading coefficient are compared in evaluating the loading. The reduced Zweifel loading coefficient limits the increase rate of DF, and the slotted profile allows a higher DF to be obtained than the datum profile when the same blade loading is suffered. • Three kinds of deficit thickness are defined and evaluated for the flow in the profile trailing edge. All the deficit thicknesses in the SS are much higher than that in the PS. The slotted profile contributes to diminishing the accumulated thick boundary layer at the profile trailing edge. As a kind of passive flow control method, it suppresses the separation and improves the performance at positive incidence, and does not deteriorate much in negative incidence.

•
The correction of loading to the momentum deficit thickness demonstrates that it is difficult to distinguish failure simply based on the DF values for a given geometry, particularly under high loadings. The Zweifel loading coefficient connects well with the low momentum deficit in the profile trailing edge. V-shaped distributions appear for each Mach number condition in the Ψ-θ de f plot. Although the drop in Zweifel loading lags behind the increase in momentum deficit, the peak of the V-shape could better suggest the design condition, under which the profile could achieve better performance both in loss and loading ability, and a region near the V-shape peak could better indicate the correct operating range. Additionally, this judgement is not limited by the incidence range and the occurrence of bulk separation under positive incidences.

•
The nondimensionalized destruction of mean mechanical energy and dissipation function are employed to analyze the loss mechanism. Either the increased incidence angle or the increased Mach number results in decreased coefficient peak values for the datum profile in the suction surface. Because as a majority part of the destruction of mean mechanical energy, the Reynolds stress that resists the mean flow deformation transforms the mean flow kinetic energy into kinetic energy for turbulence fluctuation, the boundary layer with a lower destruction coefficient has less ability to resist the adverse pressure gradient; thus, the separation occurs more easily for the datum profile under larger incidence or high Mach number conditions. The slotted profile gains the ability of the boundary layer flow near the suction surface to resist the adverse pressure gradient, and the lift off of the high destruction coefficient bound is eliminated. A reduced shear thickness and a uniformed downstream flow field is obtained by the slotted profile.