Effect of the Reynolds Number and Clearance Flow on the Aerodynamic Characteristics of a New Variable Inlet Guide Vane

Recently, a new type of low-loss variable inlet guide vane (VIGV) was proposed for improving a compressor’s performance under off-design conditions. To provide more information for applications, this work investigated the effect of the Reynolds number and clearance flow on the aerodynamic characteristics of this new type of VIGV. The performance and flow field of two representative airfoils with different chord Reynolds numbers were studied with the widely used commercial software ANSYS CFX after validation was completed. Calculations indicate that, with the decrease in the Reynolds number Rec, the airfoil loss coefficient ω and deviation δ first increase slightly and then entered a high growth rate in a low range of Rec. Afterwards, a detailed boundary-layer analysis was conducted to reveal the flow mechanism for the airfoil performance degradation with a low Reynolds number. For the design point, it is the appearance and extension of the separation region on the rear portion; for the maximum incidence point, it is the increase in the length and height of the separation region on the former portion. The three-dimensional VIGV research confirms the Reynolds number effect on airfoils. Furthermore, the clearance leakage flow forms a strong stream-wise vortex by injection into the mainflow, resulting in a high total-pressure loss and under-turning in the endwall region, which shows the potential benefits of seal treatment.


Introduction
Variable inlet guide vanes (VIGVs) are widely used in turbo-machinery to provide the desired flow angle into the downstream rotor [1][2][3], which can improve the compressor efficiency and stall margin under off-design conditions. As the stagger angle increases, the VIGV enters into a high incidence condition, and the loss coefficient starts to increase rapidly, which is detrimental to the compressor's efficiency. Therefore, some researchers proposed a new type of VIGV featuring a "dual-peak" surface velocity distribution under a high-incidence condition, which shows a wide low-loss operation range in twoand three-dimensional environments [4,5]. Related research has indicated that the chord Reynolds number and clearance flow can have an important influence on the compressor blade's aerodynamic performance [6,7]. Therefore, to reveal the influence pattern for better applications of this new type of VIGV, this paper studies the effect of the Reynolds number and clearance flow on these "dual-peak" VIGV aerodynamic characteristics.
The chord Reynolds number Re c has been known from early research to be one of the factors that have an important impact on the airfoils' aerodynamic performance [6]. In the 1950s, low-speed experiments on the compressor blade 10C4/40 P40, Lighthill blade, and NACA 65 series airfoil indicated that the airfoil loss coefficient shows a steep increase after the Reynolds number Re c decreases to a certain value, and in doing so confirmed the existence of critical Reynolds numbers [8][9][10]. However, it appears that no constant critical Reynolds number can hold for different airfoils [6], and the degradation in aerodynamic performance seems to be more remarkable for airfoils with a larger laminar boundary-layer at the design point (e.g., the Lighthill blade) [9]. Qualitative flow analysis indicates that a steep increase in the loss coefficient due to a low Reynolds number is related to the inception of boundary-layer separation [6]. The airfoil with a larger laminar region at the design point shows more serious performance degradation, due to the larger laminar separation under a low Reynolds number [6]. In addition, for a certain airfoil, the higher free-stream turbulence intensity leads to a lower critical Reynolds number [6,9], which is due to the earlier boundary-layer transition at a higher turbulence intensity.
Recent research has provided more flow details for the Reynolds number effect. In the 2000s, the high-speed experiments on a controlled diffusion airfoil conducted by Schreiber et al. [11] showed that, at a low turbulence intensity, a laminar separation bubble shortly after the maximum velocity (near 35-40% of the chord) on the suction surface could be observed for all Reynolds numbers in the tested range ( 0.7 ∼ 3.0 × 10 6 ) [11]. However, they also found that the transition point moves to the acceleration front portion at a high turbulence intensity with a high Reynolds number [11]. In 2004, Sonoda and Schreiber et al. optimized the stator airfoil's performance with a low Reynolds number by using the evolution strategy (ES) and multi-objective genetic algorithm (MOGA). The two obtained advanced compressor airfoils featured very early boundary-layer transitions (induced by an adverse pressure gradient around the leading edge), which reduced the separation flow with a low Reynolds number and obtained excellent aerodynamic performance [12,13]. This research showed the advantage of airfoil flow customization design through geometry optimization. In 2013, an experimental and numerical study (via ANSYS CFX) on a symmetric VIGV airfoil with Re c = 1.0 ∼ 7.0 × 10 5 by Handel and Barthmes indicated a critical Reynolds number of 2.0 × 10 5 at β s = 90 o , 80 o , and 70 o because of the inception of a separation region [14,15]. The reason for the over-turning phenomenon with a low Reynolds number was revealed and explained by the increase in the effective camber due to the boundary-layer thickening [14,15]. To further improve the low-Reynolds number airfoil's performance, in 2019, Yang et al. used local dynamic surfaces as an active flow control method to suppress the laminar flow separation on the suction surface in a turbine cascade under a low Reynolds number [16]. In 2020, a novel passive control method was proposed by Arif et al. [17], which involved setting a short elastic panel flush mounted on the suction surface of an NACA 0012 airfoil under a low Reynolds number and showed the effectiveness of tonal noise reduction by the absorption of the energy of natural instabilities that evolve in the laminar boundary layer. In 2021, research on the surface roughness effect in a subsonic compressor airfoil conducted by Wang et al. [18] via simulations indicated that the surface roughness determined the loss generation process under a low Reynolds number by affecting the laminar separation bubble structure and the near-wall turbulence level.
The Reynolds number's effect is different due to the difference of the aerodynamic design feature of airfoils. Most of the aforementioned works focused on the compressor airfoil or wing airfoil. Therefore, to provide more information on using this "dual-peak" cambered VIGV, this paper discusses the effect of the Reynolds number Re c and the clearance flow. Numerical calculations were successfully applied in aerodynamic field research [15,18,19], and, with proper setups, a good match with experimental results was obtained, providing a convenient method for systematic research. In this paper, firstly, the effect of the Reynolds number on the performances and boundary-layer behaviors of two "dual-peak" VIGV airfoils is explored using CFX with the RANS method after verification with experimental data. To reveal the related flow mechanism, a boundary-layer analysis is conducted on the transition location, inception, and re-attachment locations of separation to explain the airfoil performance variation. The three-dimensional research confirms the Reynolds number effect observed in plane airfoil research. The loss coefficient of the small-chord VIGV (halved Reynolds number) in the main blade part is significantly higher than that of the normal-chord version at the high stagger angle points. In addition, the leakage flow will form strong stream-wise vortex, generating a high total-pressure loss in the hub and tip regions, which shows the potential benefit from seal treatment.

Research Samples
Two representative "dual-peak" VIGV airfoils (parameters given in Table 1) that show a minimum loss coefficient [5] were chosen to be researched here. The parameter definitions and geometry are shown in Figure 1a-c. The research range of chord Reynolds number (the last line of Table 1) is designed to cover the typical range that a VIGV can encounter in an application environment. The chord Reynolds number Re c is defined by The ρ, U 1 , and µ 1 in Equation (1) are density, velocity, and dynamic viscosity of inflow, respectively. In this paper, the variation of chord Reynolds number Re c is achieved by scaling the airfoil chord as shown in Figure 1d.    Figure 2b. To evaluate the airfoil performance, the total-pressure loss coefficient and deviation at the design point (∆β 2 = 0 • ), high incidence (∆β 2 = 24 • ), and maximum incidence (∆β 2 = 30 • ) were used. The axial-velocity-density ratio was set to 1.0 for the two-dimensional research.

Numerical Method and Verification
To obtain the aerodynamic performance of the airfoils, steady simulations based on the Reynolds-Averaged Navier-Stokes (RANS) method were carried out by a widely used commercial solver, ANSYS CFX. For the setups, the high-resolution advection scheme, the SST turbulence model [20], and the γ − Re θ transition model [21] were used. The physical timescale was used and set to 0.0005 s. An RMS of the residuals below 10 −6 was used as the numerical convergence criteria.
The mesh for simulations had an O4H multi-block topology ( Figure 3), generated with the IGG software package [22]. The inlet and outlet plane were located at x = −40 and 100 mm, respectively. The LE plane of all airfoils were fixed at x = 0. The streamsurface thickness was set to 0.5 mm with five equally spaced nodes. To determine the proper grid size, a grid independence test was conducted (Grids 1-4 in Table 2). For all grids, the blade surface wall cell height was set to 0.0008 mm to yield an average value of y + < 1. The wall cell height expansion ratio was set to 1.10 to increase the grid density around the blade. As shown in Figure 4, in general, with the increase in mesh nodes, the variation of profile loss coefficient ω, outflow angle β 2 , and surface isentropic Mach number Ma is gradually became smaller. The very close results between Grid 3 and 4 indicate that the mesh convergence was reached by Grid 3 (Figure 4).   As a further verification, a symmetric VIGV cascade (has a maximum relative thickness of 15% and solidity of about 1.0) with abundant measurement data from tests by Handel et al. [14] was simulated at a different Reynolds number Re c by using Grid 4 and comparing it with the experimental results. As shown in Figure 5, a good match between the simulated and experimental surface isentropic Mach number was observed. The expansion of the suction surface laminar separation bubble with the reduction in Reynolds number Re c (from 300 to 100 k) at the design point were well captured ( Figure 5a). The increasing trend of the LE separation plateau and the reduction trend of the spike Mach number (from 500 to 200 k) at the high incidence point were also predicted by the simulation, but the magnitude of the pressure plateau and suction spike was slightly underestimated because of the slightly lower stagger angle relative to the experiment (the simulation at γ = 30 • was found with a fully detached flow at low Reynolds numbers of 100 k and 200 k, so γ = 28 • was used in the simulation [15]).
For aerodynamic parameters, the variation trend of relative loss coefficient ω/ω D by simulation was consistent with the experimental results (Figure 5b,d), but some underestimation at conditions with the appearance of flow separation was observed. The variation trend of the outflow angle β 2 at γ = 30 • was also reflected by the simulation (the value is approached if corrected by 2 deg, the difference of stagger angle). In general, the numerical method with the above setups can predict crucial phenomena such as the expansion of the laminar separation bubble and the variation trend of the loss coefficient with acceptable accuracy, which can be used in the following numerical calculations.

Loss Coefficient and Deviation
As shown in Figures 6 and 7, the variation of the Reynolds number Re c has significant impact on airfoil performance. With the reduction in the Reynolds number, the airfoil loss coefficient ω and deviation δ first increase slightly and then enters a high growth rate in the low range of Re c . For Case 1 (Ma 1 = 0.45) and Case 2 (Ma 1 = 0.60), the Reynolds number, corresponding to the rapid performance degradation, is about Re c = 1.5 × 10 5 for the design point and 2.0 × 10 5 for high incidence and max incidence (marked by the arrows).

PEER REVIEW
7 of 17 ted (Grid 4) and experimental results [14] for the effect of chord Reynolds number Re on the symmetric he design (γ = 0°) condition and (c,d) for the high incidence condition (γ = 20°).

Loss Coefficient and Deviation
As shown in Figures 6 and 7, the variation of the Reynolds number has significant impact on airfoil performance. With the reduction in the Reynolds number, the airfoil loss coefficient and deviation first increase slightly and then enters a high growth rate in the low range of . For Case 1 ( = 0.45) and Case 2 ( = 0.60), the Reynolds number, corresponding to the rapid performance degradation, is about = 1.5 × 10 for the design point and 2.0 × 10 for high incidence and max incidence (marked by the arrows).

Flow Mechanism: Pattern of Transition and Separation
To reveal the flow mechanism behind the effect of chord Reynolds number Re c , this section will focus on the boundary-layer behaviors and flow field with emphasis on the key flow phenomena. In the analysis, the skin friction coefficient C f was used to determine the location and scale of boundary-layer separation, which is calculated by C f = µ(∂u/∂y) wall / 0.5ρ 1 u 2 1 inlet main flow (2) The relative momentum thickness δ 2 /c measures the momentum loss of boundarylayer and can be used to evaluate the generation of total-pressure loss [7], which is defined by The design point (∆β 2 = 0 • ) and maximum incidence point (∆β 2 = 30 • ) are two important operation points. Figure 8 shows the evolution of the suction surface boundary layer with a change in the chord Reynolds number Re c . Figure 9 compares the important boundary-layer parameter distribution of Case 1 at 4 Reynolds number (SP1 to SP4). As shown in Figure 8, the effect of the Reynolds number Re c on boundary-layer development at the design point can be divided into two phases: the laminar transition phase and the laminar separation phase.
I. Laminar transition phase. In the high Reynolds number range (Re c > 2.6 × 10 5 for Case 1, and > 3.2 × 10 5 for Case 2), the suction surface boundary layer keeps a laminar state from the LE stagnation point to the initial point of transition x Tr.Inp. . With the reduction in the Reynolds number Re c , the initial point x Tr.Inp. and the accomplished point x Tr.Acp. of the transition region both move towards TE nearly linearly. Therefore, the length of the turbulent boundary-layer region L Tur. decreases 17.2% in chord for Case 1 and 12.8% in chord for Case 2, which tends to reduce the loss level. However, on the other side, according to the boundary-layer equation [7] (which can also be observed by comparing SP1 and SP2 in Figure 9c), the reduction in Reynolds number Re c will lead to a thicker momentum thickness δ 2 /c of the laminar boundary-layer portion, namely, increasing the loss level. The antagonism between the two factors (the length of the turbulent portion and the thickness of the laminar portion) explained why the airfoil performance is not sensitive to the Reynolds number in this phase.  II. Laminar separation phase. Once the chord Reynolds number Re c reduces to a certain value (1.9 × 10 5 for Case 1, and 2.4 × 10 5 for Case 2), the laminar separation occurs and becomes the dominant factor affecting the loss level. As shown in Figure 8, with the reduction in the Reynolds number Re c , the laminar separation inception point x Inp. moves upstream towards LE; meanwhile, the re-attachment point x Ret. moves downstream towards TE, leading to a rapid increase in the laminar separation length L LAS. and the momentum thickness δ 2 /c (Figures 8 and 9c). This explains the high growth rate of the loss coefficient in the low Reynolds number range ( Figure 6).

•
Maximum incidence (∆β 2 = 30 • ) For the maximum incidence (∆β 2 = 30 • ), Case 1 (Ma 1 = 0.45) is found to have a complex separation pattern with 3 phases: LE separation, LE and secondary separation, and secondary separation (Figure 10a): LE separation. In the high Reynolds number range (Re c > 2.7 × 10 5 ), a laminar separation region is formed in the LE region due to the adverse pressure gradient behind the LE suction spike. Although the reduction in the Reynolds number Re c in this phase is considerable (4.5 to 2.7 × 10 5 ), the length of the LE separation L LES only shows a slight increase of 2.4% in chord length, so the growth rate of the loss coefficient is low in this phase (Figure 10a). II.
LE and secondary separation. In the intermediate Reynolds number range (1.3 × 10 5 < Re c < 2.7 × 10 5 ), a secondary separation region appears at about 15% chord location, on the suction surface (Figure 10a). The boundary-layer separation region identified by the laminar separation criteria, i.e., shape factor H b > 3.5 (Figure 12a), fits well with the range identified by the skin friction coefficient C f (Figure 12b), which indicates that it is a laminar separation. In this phase, with the reduction of Re c , the LE laminar separation region decreases 1.4% in chord length; meanwhile, the secondary laminar separation region extends 6.3% in chord length, so the airfoil loss coefficient keep increasing with a rate higher than phase I. III.
Secondary separation. In the low Reynolds number range (Re c < 1.3 × 10 5 ), it is found the LE laminar separation disappears. This is because the relative thickness of boundary layer (δ 1 /c) will be larger at low Reynolds number and this changed the effective aerodynamic shape of the airfoil (Figure 11a), which results in lower peak Mach number Ma p (Figure 10b). The decreased peak Mach number reduces the strength of subsequent adverse pressure gradient, and then the LE laminar separation disappears. As shown in Figure 10a, with the reduction of Reynolds number Re c from 1.3 to 0.3 × 10 5 , the separation inception point x SES.Inp. remains nearly constant, but the re-attachment point x SES.Ret. moves towards downstream from x = 0.21 to x = 0.35, leading to a rapid extension of the separation region. In addition, the strength and height of the boundary-layer velocity loss is also significantly increased (comparing SP3 with SP4 in Figure 11b). These two factors (the increase in the length and height of the laminar separation region) causes the high growth of the loss coefficient in this phase (Figure 10a).
For Case 2, with a higher Mach number Ma 1D = 0.60, the flow separation pattern has only one type: LE separation (Figure 10c). There is an "insensitive range" (Re c > 3.4 × 10 5 ) where the LE laminar separation length L LES is almost unaffected by the Reynolds number (Figure 10c). However, once the Reynolds number Re c < 3.4 × 10 5 , the growth rate of the loss coefficient substantially increases due to the rapid increase in length and height of LE laminar separation region, which is same to Case 1.

Methodology
Compared with airfoil flow, the three-dimensional flow in the guide vane is more complex because of the flow property variation in the radial direction, interactions between adjacent stream surfaces, and the endwall effect. This section focuses on the effect of the Reynolds number and the clearance flow on the VIGV in a three-dimensional environment.
The dual-peak VIGV with an optimized LE portion, given in previous paper [5], was used as the blueprint to obtain the small-chord version by scaling the airfoils from root to tip with a constant ratio of 0.5 (Figure 13a,b). The numerical calculations in this part are conducted by the software ANSYS CFX with a steady and unsteady RANS method. The high-resolution advection scheme, the SST turbulence model [20], and the γ − Re θ transition model [21] are used. The shroud and hub are the same for the normal-and small-chord VIGVs. To consider the influence of leakage flow, the tip and hub clearances are designed to a gap of 0.5 mm at a minimum state (Figure 13c). Pure axial flow with a total pressure of 101.325 kpa and a total temperature of 288.15 k is used as the inlet boundary condition. The mass flow rate is specified as the outlet boundary condition. For the unsteady RANS calculation (URANS), the numerical converged RANS result is used as the initial values to start the calculation. The same simulation mesh is used in the RANS and URANS calculations. The mesh has an O4H topology with a total node of 4971.9 k ( Figure 14). For one layer, the mesh has 377 stream-wise nodes (233 on the profile surface) and 81 pitch-wise nodes (53 across the O-block). To research the Reynolds number and clearance flow effect on the three-dimensional VIGV performance in the entire operation range, five operating points from a design condition (P1) to a high stagger angle condition (P5) are calculated (Table 3).    Figure 15 shows the operation schedule, Reynolds number, the radial distribution of the inlet Mach number Ma 1 , and the outflow angle β 2 of the normal (×1.0) and small-chord versions (×0.5) at each operation point (same as Schedule 1 in [5]). The Reynolds number Re c at the mid-span (R = 0.5) of the normal-chord VIGV is 4.2 × 10 5 at the design point and varied to 2.6 × 10 5 at the high stagger angle point P5. With a chord scale of 0.5, the Reynolds number Re c of the small-chord version is halved (Figure 15b). The inlet Mach number Ma 1 of the two versions are very close with a difference within 0.01 at most of the blade span (Figure 15c). The slightly higher inlet Mach number Ma 1 in the root part ( R = 0 ∼ 0.1) of the small-chord VIGV is caused by the contraction of the flow path because its inflow measurement plane (x = −10 mm) is 10 mm towards the downstream relative to the normal-chord version.

Reynolds Number Effect in Three-Dimensional VIGV
As shown in Figure 16a-e, the small-chord version (lower Reynolds number) has a higher loss coefficient ω in the main blade part ( R = 0.1 ∼ 0.85), especially at the high stagger angle points (P4 and P5). Specifically, the averaged loss coefficient of the small-chord VIGV in the main part ( R = 0.1 ∼ 0.85) is 44.3% and 80.5% higher than the normal-chord version at P4 (∆γ = 30 deg) and P5 (∆γ = 35 deg), respectively. The unsteady RANS (URANS) calculation confirms the effects of the chord Reynolds number reduction and clearance flow on the total-pressure loss level (Figure 16d-f). However, the specific loss coefficient distribution along the blade span is different from the RANS method. For example, the magnitude of the loss coefficient bulges (around R = 0.2 and 0.4) are reduced in the URANS calculation (Figure 16e), which reflects that the radial flow mixing is stronger in the URANS calculation.
For the flow guidance ability, the small-and normal-chord versions are almost the same in the mainflow region ( R = 0.1 ∼ 0.85) except at the operation point P5 (Figure 15d). The small-chord version shows a 1.6 • over-turning relative to the normal-chord version around span height R = 0.2. Flow visualization indicates this flow over-turning is induced by the increase in effective camber due to the boundary-layer thickening around the former suction surface.
However, in the near-wall region (hub region R = 0 ∼ 0.05; tip region R = 0.95 ∼ 1), the loss level of the normal-chord version is higher than the small-chord version (the enlarged drawings in Figure 16a-e), which is because of the stronger leakage flow due to the larger clearance at LE and TE. The net effect of the mainflow and endwall regions are given in Figure 16f. The passage-averaged total pressure recovery coefficient σ of the small-chord VIGV shows a higher value with 0.12 and 0.36% at operating points P1 and P5, relative to the normal-chord version, which indicates that the reduction in the Reynolds number can cause considerable degradation in three-dimensional VIGV performance, especially at the high stagger angle.

Impact of Endwall Leakage Flow
To reveal the impact of the leakage flow, Figure 16 displays the calculated results of the normal-chord VIGV without clearance. For all simulated operation points, the influence range of the clearance flow is the endwall region, namely, R = 0 ∼ 0.1 at the hub and 0.85 ∼ 1.0 at the tip ( Figure 16). The clearance VIGV has a much higher loss coefficient the endwall region (Figure 16b-e), which leads to a 0.2 percentage reduction in total-pressure recovery coefficient σ at the high stagger angle points (P4 and P5). Flow visualization shows that the entropy production in the leakage flow vortex forming and mixing with the mainflow leads to the high total-pressure loss in the corner region (Figures 17 and 18). Furthermore, the total-pressure loss decomposition (to (1) the profile loss, (2) the endwall loss due to friction, (3) the corner loss due to leakage flow and the intersection of the vane and the wall boundary layer), by the method of Meyer and Engel [23], indicates that the proportion of corner loss increases from the design point to the high stagger angle point and is significantly larger with the existence of leakage flow. From operation point P2 to P5, the ratio of corner loss to total loss is about 4.3-6.7% in the VIGV without clearance and increases to 14.9-19.8% in the VIGV with clearance. As shown in Figure 18a, the leakage flow, driven by the blade surface pressure difference, shows a rapid increase in tangential velocity v at the edge of the pressure surface (pressure potential energy into kinetic energy) and injects into the suction side mainflow, generating a strong, anti-clockwise axial vortex (Figure 18b). Because of the decreased pressure difference from LE to TE, the strength of the clearance flow gradually reduces towards TE (Figure 18a,c). These results indicate the potential benefit in the loss reduction from seal treatment applications and that it should be placed in the former portion of the VIGV.  For the outflow angle, the VIGV with and without clearance are totally different near the endwall region (marked by the dashed line in Figure 15d). The VIGV with clearance shows insufficient turning (in other words, a large deviation) near the endwall region ( ∆H = 0.5% ∼ 5% for hub; ∆H = 3% ∼ 6% for tip) because of the loss of the tangential velocity v caused by the clearance leakage flow (Figure 19a). For the VIGV without clearance, the tangential velocity v in the endwall region is close to the design value (Figure 19a), but the axial velocity w is rapidly decreased (Figure 19b) because of the endwall boundary layer, which leads to the over-turning of the outflow (about 3% blade height for hub and 4% for shroud). For the surface pressure distribution, two major differences exist. First, in the VIGV with clearance, the LE suction peak collapses (evaluated by the Ma p , decreased by about 20%) compared with the VIGV without clearance. This can be illustrated by considering the clearance as a communication passage that unloads the pressure difference around the leading-edge region (Figure 18a-c), resulting in the collapse of the suction spike. Second, in the VIGV with clearance, a low-pressure region appears on the former suction surface (peak at x/c = 15%, Figure 19b), which is induced by the accumulation of the high velocity leakage flow.

Conclusions
This paper investigates the effect of the Reynolds number and the clearance flow on the aerodynamic characteristics of a new type of VIGV. Key conclusions are summarized as follows:

•
With the reduction in Reynolds number Re c , the airfoil loss coefficient ω and deviation δ first increase slightly and then enter a high growth rate in the low range. For the two researched airfoils, the initial Reynolds number to the rapid performance degradation is 1.5 × 10 5 for the design point and 2.0 × 10 5 for the high incidence and max incidence.

•
For VIGV airfoils at the design point, in the high Reynolds number range, the antagonism between the reduction in the turbulent portion length and the thickening of the laminar portion makes the airfoil performance insensitive to the Reynolds number.
Once the Reynolds number reduces to a certain value, the laminar separation bubble appears on the rear portion of the suction surface and extends rapidly because of the downstream migration of the re-attachment point, which leads to the steep increase in loss coefficient. • For VIGV airfoils at maximum incidence, the flow separation pattern is not the same for the two researched airfoils. However, for both airfoils, the rapid increase in the loss coefficient at a low Reynolds number is caused by the increase in the length and height of the laminar separation region.
The three-dimensional research confirms the Reynolds number effect. The average loss coefficient of the small-chord VIGV (Reynolds number halved) in the main part ( R = 0.1 ∼ 0.85) is 44.3% and 80.5% higher than the normal-chord version at the high stagger angle points P4 (∆γ = 30 deg) and P5 (∆γ = 35 deg), respectively. In the VIGV with clearance, the entropy production in the leakage flow vortex mixing with the mainflow leads to the high loss level in the endwall region ( R = 0 ∼ 0.1 and 0.85 ∼ 1.0). The strength of the clearance flow gradually decreases towards the trailing edge. This indicates the potential benefit in loss reduction from seal treatment applications and that it should be placed in the former portion. Data Availability Statement: The data of the researched VIGV airfoils geometry are available from the corresponding author upon request.

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