Abstract

An axial fan in the general industry needs to secure stable characteristics in a wide operating range through an antistalling process, concerning its high demands, energy consumption, safety, and lifetime. As a “simple” and “original” method to control the stall of an axial fan, two-dimensional plates so-named anti-stall fin (ASF) were attached inside the casing and toward the shaft. The sensitivity analysis for design variables was conducted through the 2k full factorial design method, and the optimization was performed using the response surface method (RSM). In the case of an axial fan without the ASF, the backflow gradually increased in the spanwise direction toward the hub and the streamwise direction to the inlet passage (upstream), as the flowrate decreased, and the rotating stall was developed through the quickening-growing-setting period. However, an ASF-attached axial fan stably recovered performance degradation in the stalling flow rates and allowed to form a negative slope to 0.5Φd. This study was mainly based on numerical simulations, while the performance before/after application of the ASF was validated through experimental tests.

1 Introduction

Axial fans are in high demand for air conditioning systems such as vehicles, aircraft, home appliances, electronic devices, and underground spaces, with their structural simplicity. Especially for applications in subways or tunnels, where demands cannot be even, the axial fans need to secure stable characteristics in a wide operating range through an antistalling process. Here, “stall” is referred to as a phenomenon that can be frequently contained in fluid machinery and is due to an increase in incidence angle at the low flow rates. Figure 1 shows the inlet velocity triangle under the ideal condition, where u, c, w, i, and β denote circumferential velocity, absolute velocity, relative velocity, incidence angle, and flow angle, and subscripts 1, m, and b denote blade inlet, meridional component, and blade, respectively. If the incidence angle increases beyond a specific limit, the fluid machinery can no longer guarantee stable operation, i.e., in the stall range, the performance curve (flowrate-pressure; Q-P) exhibits a positive slope with a decrease in pressure [1,2] and is generally characterized by noise [3,4], vibration [5,6], and system instability [79]. In our field, the stall stands as a profound problem.

Fig. 1
Velocity triangle and incidence angle at blade inlet (i1=β1b−β1; difference between the inlet blade angle and the inlet flow angle): Design flowrate (green) and Lower flowrate (red)
Fig. 1
Velocity triangle and incidence angle at blade inlet (i1=β1b−β1; difference between the inlet blade angle and the inlet flow angle): Design flowrate (green) and Lower flowrate (red)
Close modal

Some previous studies [1012] experimentally and numerically focused on the internal flow field to visualize the inlet flow pattern in the stalling flow rates, and “rotating stall” was pointed in common. The rotating stall could be defined as a vortex's cluster which was following the blade's rotational direction in the inlet passage (front part of the blade) and was facing the downstream: it could be confirmed when “backflow” was premised from the blade tip; here, the backflow should be located at the spans higher than the rotating stall's location and could obtain a circumferential velocity component. This flow pattern could cause performance degradation in the stalling flow rates and be closely related to noise, vibration, and system instability. In particular, the backflow became more assertive at lower flow rates [1317]. Here, in the case of centrifugal fluid machinery, the rotating stall may be observed only between the blade passages rather than the inlet passage; this study relates to an axial fan.

The researchers have been trying to control the stall for decades; the main focus was to disrupt the formation of backflow or rotating stalls near the shroud. Lim et al. [18] configured the sliding rollers as a part of the casing treatment and provided them outside the casing to be arranged as an annular row so that grooves (or slots) were formed at the moment of stalling condition. The groove, which had been proposed since the 1950s, allowed the high-pressure air at the rear of the blade to be transferred back to the front of the blade; however, it could increase the clearance between the casing and blade tip, which could affect the performance near the design flowrate [19]. Thus, they had to drive the sliding rollers again to obliterate the grooves in case of no stall. Gentry et al. [20] made the grooves permanent with minimized depth and length. Jung et al. [21] perforated the casing at the front and rear of the blade and constructed a seesaw-type link to control the perforation. The perforation could be opened to release the high-pressure air at the front of the blade. McKelvey [22] provided an annular chamber and padded it outside the casing so that the high-pressure air at the rear of the blade could be retransferred to the front of the blade; this recirculation method was reviewed with Bianchi et al. [23] and was similar to the proposal of Dygert and Bushnell [24]. Khaleghi [25] and Salunkhe et al. [26] introduced an air injection method to supply high pressure to the low-pressure region. In terms of blades, Jung et al. [27] suggested a recessed cavity inside the blade tip; this can be understood as the reverse concept against the grooves above. The inlet guide vane could be considered as an additional system so that the absolute flow angle at the blade inlet could be adjusted to secure the stable operation [28]. A bump-shaped obstacle was added on the blade surface [29] or on the hub surface [30] to affect the flow patterns locally, and a sinusoidal-like blade leading edge (LE) was applied to enhance the lift coefficient [31]. These deep efforts eventually paid off in the antistalling performance; however, each of them might face major or minor disadvantages in the case by case:

  • operating devices and systems; cost and time

  • complicated design

  • installation space and maintenance

  • performance degradation (or change) from design specification

The highly inspired focus of this study was a “simple” and “original” method to compensate for the above disadvantages based on the antistalling process of an axial fan. Thus, two-dimensional plates were attached inside the casing and toward the shaft. Since the plates should be fixed and semipermanent as a stator, it could be a kind of passive control for the stall, and it has been named as “anti-stall fin (ASF)” by ourselves. First, the sensitivity analysis for design variables was conducted through the 2k full factorial design method, and the optimization was performed using the response surface method (RSM). The optimized ASF was directly applied to an axial fan and was compared to the non-ASF state. The analysis was focused on the internal flow field and the fast Fourier transform (FFT) results: it was mainly based on the numerical simulations; the performance before/after application of the ASF was validated through experimental tests. If performance degradation in the stalling flow rates is improved without any instability, an efficient operation can be secured through expansion of the operating range (stall margin). Additionally, noise reduction can be expected through suppression of the unfavorable flow patterns near the blade LE. Meanwhile, the antistalling method with ASF of this study was recognized for its originality from a patent application [32].

2 Subject of Study

2.1 Axial Fan Unit.

The axial fan for applying the ASF was being used in the general industry, as shown in Fig. 2; it is a prototype of this study. The assembly was made up of a pair of 10 blades and 11 guide vanes. A bell mouth and hub cap were attached as a front part. The blades had three-dimensional geometry based on the NACA airfoil series at each span. The guide vanes exhibited two-dimensional geometry, which means the same sectional configuration to the spanwise direction. The tip clearance (δt) occupied a range of 0.99–1.00 based on hub-to-shroud normalized span (r*, 0 to 1). Table 1 lists the design specifications and parameters; here, subscript d could imply the design point, but omitted, e.g., Φd= 0.285. The specific speed (Ns), flow coefficient (Φ), and pressure coefficient (Ψ) were defined as
Ns=ωQ(P/ρ)34
(1)
Φ=cm2u2
(2)
Ψ=P12ρu22
(3)

where ω, Q, P, ρ, and subscript 2 denote the angular velocity, volume flowrate, total pressure, air density, and blade outlet, respectively and u2 was assigned for the blade tip. Meanwhile, the axial fan showed a peak point at 0.8Φd on the Q-P (or Φ-Ψ) curve and contained a positive slope in the lower flow rates than 0.8Φd [33].

Table 1

Design specifications and parameters of the axial fan

Specification or parameterValue (unit)
Specific speed (Ns)7.77
Flow coefficient (Φ)0.285
Pressure coefficient (Ψ)0.117
Rotational speed (N)1,470 (rpm)
Hub ratio (rh/rs)0.44
Tip clearance ratio [δt/(rsrh)]0.01
Soliditya (C/S)0.769 (hub), 0.155 (shroud)
Setting angleb49.7 (hub), 23.1 (shroud) (degree)
Number of blades and guide vanes10 and 11
Airfoil seriesNACA 3512
Specification or parameterValue (unit)
Specific speed (Ns)7.77
Flow coefficient (Φ)0.285
Pressure coefficient (Ψ)0.117
Rotational speed (N)1,470 (rpm)
Hub ratio (rh/rs)0.44
Tip clearance ratio [δt/(rsrh)]0.01
Soliditya (C/S)0.769 (hub), 0.155 (shroud)
Setting angleb49.7 (hub), 23.1 (shroud) (degree)
Number of blades and guide vanes10 and 11
Airfoil seriesNACA 3512
a

Blade chord length (C) divided by inlet blade-to-blade pitch (S).

b

Tangential definition.

Fig. 2
Typical assembly of an axial fan
Fig. 2
Typical assembly of an axial fan
Close modal
Table 2

Details of design parameters for the anti-stall fin (ASF)

ParameterValue (unit)Remarks
Radial length (lr/D2)0.01–0.035variable, center: 0.0225
Axial length (la/D2)0.05–0.15variable, center: 0.1
Axial gap (δ/D2)0.01–0.04variable, center: 0.025
Number of fins13
Ellipse ratio1 (semicircle)both leading and trailing edge
ParameterValue (unit)Remarks
Radial length (lr/D2)0.01–0.035variable, center: 0.0225
Axial length (la/D2)0.05–0.15variable, center: 0.1
Axial gap (δ/D2)0.01–0.04variable, center: 0.025
Number of fins13
Ellipse ratio1 (semicircle)both leading and trailing edge

2.2 Anti-Stall Fin.

The ASF was designed to be attached inside the casing as a stator so that no additional space was required for installation (see Fig. 3). Since the ASF's axial directionality (angle; βASF) inevitably causes the absolute flow angle at the blade inlet and leads to the decrease (or change) in performance even near the design flowrate, the ASF was designed so that its surface was not visible in front of an axial fan, i.e., the ASF exhibited two-dimensional geometry just toward the shaft. The basic principle was to prevent the development of backflow and redirected the circumferential velocity component to the axial direction. The features that could be considered from the above concept are:

Fig. 3
Schematic drawing of anti-stall fin (ASF) on meridional plane
Fig. 3
Schematic drawing of anti-stall fin (ASF) on meridional plane
Close modal
  • no operating devices and systems

  • simple configuration

  • immediacy (on-site welding or fastening); semi-permanent

  • guaranteed performance based on design specification

  • regardless of material (iron, rubber, plastic, etc.)

where any material, including a flexible one, would be available under the premise that the ASF and casing are completely attached. Nevertheless, before applying the ASF, a previous consideration on the meridional plane is strongly recommended to prevent contact between the ASF's TE and rotor's LE; the same is true in the case of high-pressure fans (or fluid machinery) that may have arguments with thrust.

Table 2 lists the design parameters of ASF. From an empirical perspective, the axial gap (δ) was selected as the most critical variable; it was based on the fact that the backflow from the blade LE mainly causes instability at low flow rates. In addition, since the backflow occupies a wider region inside the flow passage as the flowrate decreases, the radial and axial length (lr and la) were considered as the variables: the above three parameters were selected as the main variables for ASF (see Fig. 3) and were normalized to the fan diameter (D2). The number of fins (13) was not a variable; the number of fins was considered to be more than the number of blades (10) to sufficiently control the unstable flow associated with the backflow from each blade LE and to avoid overlap with the number of blades and guide vanes (11). The ellipse ratio was selected as 1 (semicircle) for both LE and TE, and the angular distribution (θASF or βASF) should not be considered with two-dimensional geometry. The thickness was designed to be constant in terms of general application; it was the same dimension as that of the guide vane in this study. Besides the main variables (lr, la, δ), the parameters that can affect the function of the ASF may be the number of fins and the tangential angle (θASF); these may be addressed in further studies.

3 Basis of Analysis Method

3.1 Experimental Test.

Most of all, the experimental process and facility of this study fully complied with the international standard, ANSI/AMCA 210-07 [34]. Figure 4(a) shows a photograph of each axial fan connected to the test facility: the case of “none” represents a typical assembly as shown in Fig. 2; in the case of “ASF attached,” only the ivory-colored duct was replaced with the ASF. The rotor (including blades and hub cap) was manufactured using five-axis machining technology. The 0.5 hp (0.373 kW/h) class motor was selected with three-phase electric conversion and four poles and was placed inside the hub duct.

Fig. 4
Experimental test facility: (a) Photograph of axial fan unit for each case of none and ASF attached and (b) Schematic diagram of outlet chamber setup
Fig. 4
Experimental test facility: (a) Photograph of axial fan unit for each case of none and ASF attached and (b) Schematic diagram of outlet chamber setup
Close modal

Regarding further details on the experimental test, the outlet chamber setup was adopted as shown in Fig. 4(b). A straight duct that has twice the axial length for the fan diameter was connected between the fan outlet and chamber inlet. The flow settling means had secured the proper porosity [35]. Relative humidity, barometric pressure, and dry-bulb temperature were measured to calculate the density; here, the density and rotational speed were respectively converted as the same values with the numerical setup. The flowrate was adjusted with nozzles and was calculated from the differential pressure (ΔPs). The flowrate range that could not be measured with a combination of nozzles was handled with a servoblower. The pressure and rotational speed were measured with pressure manometers and a laser tachometer (or a stroboscope). Actually, all measuring instruments in this test facility have been subject to annual certification so that the uncertainty for pressure manometers, stroboscope, and dry-bulb temperature detectors was 0.001–0.005 kPa for the range of 0–1.33 kPa, 0.1–1 rpm for the range of 40–35,000 rpm, and 0.07 °C for the range of 0–60 °C, respectively. The data acquisition system was incorporated with an averaging function for all outputs (measured data) for 6 s. The average measurements during the test for dry-bulb temperature at the front of the nozzles inside the chamber, dry-bulb temperature at the rear of the nozzles inside the chamber, dry-bulb temperature in the laboratory, barometric pressure in the laboratory, and relative humidity in the laboratory were 24.24 °C, 24.08 °C, 24.92 °C, 101.1 kPa, and 46.01%, respectively.

3.2 Numerical Simulation

3.2.1 Governing Equation.

The Reynolds-averaged Navier–Stokes equations were employed in the three-dimensional flow field and were discretized using the finite volume method, and the time variation terms were added for the transient state (unsteady-Reynolds-averaged Navier–Stokes); the conservation of energy was ignored because this study stands for the isothermal (25 °C) condition. The conservation of mass (Eq. (4.1)) and momentum (Eq. (4.2)) could be stated as follows:
ρt+(ρUi)xi=0
(4.1)
(ρUi)t=Pxi+xj[μ(Uixj+Ujxi23Urxrδij)]+ρFi
(4.2)

where t, U (could be substituted as V or W), x (could be substituted as y or z), and Fi denote the time, velocity, coordinate, and body force, respectively, and the terms in square brackets denote the viscous stress tensor (τij). Here, there was no change in density over time; the maximum Mach number at the blade tip was estimated as 0.09 at 25 °C (subsonic flow; Mach number < 0.3). Meanwhile, a high-resolution discretization method was adopted based on the second-order approximation, and the root-mean-square (RMS) residuals were kept within 1.0×104 and 1.0×105.

3.2.2 Turbulence Model and Intensity.

The k-ω-based shear stress transport standard (SST Std) model is known to be quite suitable for rotating machinery: it had been developed to provide accurate predictions in adverse pressure gradients, especially for onset and amount of the flow separation; however, it included transport effects to the eddy-viscosity formulation [36], and flow separation from smooth surfaces could be exaggerated under the influence of adverse pressure gradients. To enhance turbulence levels in the separating shear layers emanating from walls, a modified SST Std model was suggested and was referred to as “shear stress transport reattachment modification (SST RM)” [3739]. The SST RM model considered an additional source term for k-equation to secure the ratio of turbulence production to dissipation, which might be greatly exceeded in large flow separation; it means that the SST RM model is more suitable when focusing on the range of low flow rates as in this study. Meanwhile, the SST Std and RM models had little effect on weak separation such near the design flowrate (0.8–1.2Φd) [40]. As further presented in Sec. 4.2, the SST RM model predicted relatively accurate performance curves at 0.5–1.2Φd, including the stalling flow rates. Obviously, there was a deviation in the stalling flow rates; it was not enough to be of concern, while the numerical accuracy may be enhanced with an adjusted ratio (turbulence production to dissipation).

The turbulence intensity (Tu) was approximately 4.84% at the inlet boundary
Tu=13[(UU¯)2+(VV¯)2+(WW¯)2](U¯2+V¯2+W¯2)=23k(U¯2+V¯2+W¯2)
(5.1)
k=12[(UU¯)2+(VV¯)2+(WW¯)2]
(5.2)

where U, V, W, U¯, V¯, W¯, and k denote the instantaneous and averaged velocity of x, y, z-direction in the orthogonal coordinate system, and turbulence kinetic energy, respectively. According to a general recommendation, it was appropriate that the level of turbulence intensity was selected as “medium,” which corresponds to the range of 1–5%; this range is recommended for flow in not-so-complex devices such as large pipes, fans, wind tunnels, or ventilation flows [34].

3.2.3 Computational Domain and Setup.

Figure 5(a) depicts the whole flow passage: the inlet passage included the ASFs and was sufficiently extended to account for rotating and recirculating flow under the stalling condition; the rotating passage included the blades, and the counter-rotating condition was given to the shroud wall; the outlet passage included the guide vanes. Here, the bell mouse and hub cap could not be considered because their effects were insignificant compared to the straightly extended passage [41]. The atmospheric pressure and mass flowrate were given to the inlet and outlet boundary, respectively, and the working fluid was air at 25 °C. The wall function was selected as automatic, and the boundary walls were treated as smooth and nonslip conditions. The stage (mixing-plane) method with the physical time-step (0.001 s) was applied to each interface of the steady-state simulation, while transient rotor–stator method was given to the transient one. The total time for one revolution was approximately 0.0408 s, and transient data were obtained for each time-step (every 3 deg) with a maximum loop of 5; i.e., 5 iterations per time-step.

Fig. 5
(a) Computational domain and grid system and (b) results of grid test based on grid convergence index method (GCIfine21= 0.000297 and 0.000006 for pressure and shaft power coefficient)
Fig. 5
(a) Computational domain and grid system and (b) results of grid test based on grid convergence index method (GCIfine21= 0.000297 and 0.000006 for pressure and shaft power coefficient)
Close modal

3.2.4 Grid System.

The grid test was conducted as shown in Fig. 5(b); it was at the design flowrate for the case of none. A grid refinement technique, which was established by Roache [42], was employed to quantify the grid convergence. As a result, the convergence corresponding to the N1 set was evaluated as 0.000297 and 0.000006 for the pressure and shaft power coefficient (λ)
λ=L0.5ρA2u23
(6)

where L and A2 denote shaft power and fan outlet area (πrs2πrh2). The evaluated values above (0.000297 and 0.000006) were considerably lower than the self-proposed criteria [43]; the numerical results were hardly affected with the N1 set, and the grid system was applied with the same topology corresponding to N1 set.

3.2.5 Other Information.

A commercial software, ansyscfx 19.1, was basically used for simulations. minitab 17 (Minitab Inc., Pennsylvania, PA) and origin 2019b (OriginLab Corp., Northampton, MA) were utilized to extract the optimized point with the RSM and perform the FFT from the fluctuation signal. The workstation has the specifications as follows: Intel® Xeon® CPU E5-2680 v2; clocked at 2.80 GHz with dual processor; random access memory with 80 GB; 64-bit operating system; parallel computations. The computational time for a set of steady-state and a revolution of transient simulations was approximately 26 and 30 h, respectively.

4 Results

4.1 Design and Optimization of Anti-Stall Fin.

The effects of main variables (lr, la, δ) on antistalling performance were analyzed from the 2k full factorial design method, one of the traditional methods for the design of experiments; this method basically enables the sensitivity analysis for each variable and understands the interactions between each variable to arrive at the optimization [44]. The objective function (a) of this study was selected to negatively maximize the gradient from the flowrate point where the positive slope begins to be formed in the Q-P (or Φ-Ψ) curve for the case of none, i.e., such that a positive slope cannot be formed
a=Ψx2Ψx1Φx2Φx1(x2=x1+101)
(7.1)

where subscripts x1 and x2 denote the multiple of normalized flowrate based on the design flowrate; according to the descriptions below, x1 and x2 would be designated as 0.7 and 0.8, respectively.

As the first step for designing ASF, the low-high levels, which were derived from a feasibility study, were assigned to the reference set (center; see Table 2). The total run time (number of sets) based on the 2k full factorial design method was 8. Including the set of none and reference, each set was subjected to the steady-state simulations in the flowrate range, 0.5–1.0Φd. Table 3 lists the dimensions of the main variables for each set, total pressure rise at each flowrate point, and the values of the objective function: 0 and n in the “-” column mean the set of center and none; the objective function should be recognized as negatively maximizing the gradient in 0.7–0.8Φd. From the sensitivity analysis (Fig. 6(a)) for the total pressure rise, it was difficult to confirm the effect of the main variables in 0.8–1.0Φd; this suggests that ASF had little effect near the design flowrate where backflow or stall is not generally expected; at 0.8, 0.9, and 1.0Φd, each mean of pressure coefficient decreased only by 1.25, 0.20, and 0.08% compared to the set of none. However, at 0.7Φd, the higher rise was obtained with the longer radial length and the narrower axial gap; the ASF played some role in 0.7–0.8Φd; the axial gap was the most sensitive. In 0.5–0.6Φd, the radial and axial length showed more sensitive effects.

Fig. 6
Main effect plot for (a) total pressure and (b) objective function (gradient): Radial length (lr); Axial length (la); and Axial gap (δ)
Fig. 6
Main effect plot for (a) total pressure and (b) objective function (gradient): Radial length (lr); Axial length (la); and Axial gap (δ)
Close modal
Table 3

Numerical results (pressure coefficient, Ψ and objective function, a) for each set of 2k full factorial design (1–8), center (0), and none (n)

Ψ
lr/D2la/D2δ/D20.5Φd0.6Φd0.7Φd0.8Φd0.9Φd1.0Φda
10.010.050.010.1120.1320.1480.1450.1370.1170.102
20.0350.050.010.1280.1660.1560.1440.1360.1170.437
30.010.150.010.1260.1340.1500.1450.1370.1170.185
40.0350.150.010.1720.1650.1560.1440.1360.1170.428
50.010.050.040.1130.1350.1430.1430.1370.1170.004
60.0350.050.040.1190.1460.1430.1430.1360.1170.004
70.010.150.040.1240.1310.1430.1430.1370.1170.003
80.0350.150.040.1550.1460.1440.1430.1360.1170.033
00.02250.10.0250.1250.1570.1480.1420.1370.1170.200
n0.1100.1290.1420.1450.1370.1170.107
Ψ
lr/D2la/D2δ/D20.5Φd0.6Φd0.7Φd0.8Φd0.9Φd1.0Φda
10.010.050.010.1120.1320.1480.1450.1370.1170.102
20.0350.050.010.1280.1660.1560.1440.1360.1170.437
30.010.150.010.1260.1340.1500.1450.1370.1170.185
40.0350.150.010.1720.1650.1560.1440.1360.1170.428
50.010.050.040.1130.1350.1430.1430.1370.1170.004
60.0350.050.040.1190.1460.1430.1430.1360.1170.004
70.010.150.040.1240.1310.1430.1430.1370.1170.003
80.0350.150.040.1550.1460.1440.1430.1360.1170.033
00.02250.10.0250.1250.1570.1480.1420.1370.1170.200
n0.1100.1290.1420.1450.1370.1170.107
Subsequently, Fig. 6(b) presents the sensitivity analysis for the objective function. First of all, the axial gap should be less than 4% (0.04) of the fan diameter to avoid a positive slope in 0.7–0.8Φd; it was the most sensitive variable. As the radial length increased, the effect gradually became insignificant, and the axial length had no notable effect. Although the low-high levels for the radial and axial length were selected within the range predicting negative slope, it could be expected to be positive near the lower range. The sensitivity for each variable was 0.140, 0.569, and 0.078 for radial length, axial length, and axial gap, respectively, where coefficient of sensitivity was 0.05; i.e., the variables were sensitive in order of axial gap > radial length > axial length. Finally, each dimension of the optimized ASF was obtained through an RSM-based regression equation which derived from the correlation between the main variables and the objective function
a=0.1531.4lr̃0.741lã+0.58δ̃+330lr̃2+13.8lr̃lã+360.4lr̃δ̃+6.2lãδ̃
(7.2)

where “∼” imply normalization to D2, and a second-order regression came from the RSM for 2k sets including the center point; the center point implies that the response surface can be curved (nonlinear) rather than linear, and in general RSM, it is preferable to examine second-order regression rather than first-order one, which can reflect more detailed interactions. Here, the optimized dimensions for lr̃, lã, and δ̃ were 0.035, 0.15, and 0.01, respectively, and the gradient (a; objective function) was estimated as 0.442 with R2= 99.34.

4.2 Performance Characteristic.

The optimized ASF was attached to the axial fan and was compared with the case of none. Figure 7(a) shows Q-P (or ΦΨ) curves for each case of none and ASF attached; the transient data represent the average (solid-square symbol) and maximum-minimum value (bar symbol) of time-traced total pressure fluctuations in the rotor's last revolution (120 data per set) exhibiting an adequately repeated pattern, and the description would be included in Sec. 4.4. In the case of ASF attached, the positive slope contained in the case of none was completely reversed to become negative. The ASF-attached axial fan stably recovered pressure degradation in the stalling flow rates and allowed to form a negative slope to 0.5Φd. The gradient in 0.7–0.8Φd, which was selected as an objective function, was 0.428, and a deviation from the estimated value was only 0.014. Moreover, the ASF had little effect near the design flowrate. Here, transmitted data on each interface could be more detailed with the transient simulation because the time-step was more dense (see Sec. 3.2.3); this can be a main cause for the small deviation between the steady-state and transient results. Figure 7(b) depicts the shaft power coefficient (λ) in the same flowrate range. For each shaft power (L) of experimental and numerical data, it was substituted as 3VIcosθ and Tω, respectively. Accordingly, the experimental data had to take motor efficiency into account. Here, motor efficiency might not be constant depending on the load (flowrate or rotational speed), and this fluctuation is generally not easy to be guaranteed as the motor was small size. Nevertheless, the overall tendency seemed to be almost the same as Fig. 7(a). It could additionally be considered that the motor efficiency was approximately 71% (datasheet by HIGEN Motors Inc.) in this experimental condition. Table 4 lists the increase and decrease in the total pressure rise, shaft power, and total efficiency that could be obtained from the ASF. The increase in pressure could compensate for the increase in shaft power in the stalling flow rates, sufficiently.

Fig. 7
Performance curve for each case of none and ASF attached: (a) Flow rate-Total pressure and (b) Flow rate-Shaft power
Fig. 7
Performance curve for each case of none and ASF attached: (a) Flow rate-Total pressure and (b) Flow rate-Shaft power
Close modal
Table 4

Increase and decrease rates of pressure coefficient, shaft power coefficient, and total efficiency (based on the case of none, %)

0.5Φd0.6Φd0.7Φd0.8Φd0.9Φd1.0Φd
ΨExp.+28.6+28.3+11.21.82.21.6
CFD (steady)+56.6+28.5+9.71.10.40.2
CFD (transient)+48.2+30.0+10.41.30.60.5
λExp.+12.4+13.4+3.20.7+0.3+1.3
CFD (steady)+36.0+22.8+8.00.60.10
Ψ/λExp.+14.3+13.1+7.81.12.62.8
CFD (steady)+15.1+4.6+1.50.50.30.2
0.5Φd0.6Φd0.7Φd0.8Φd0.9Φd1.0Φd
ΨExp.+28.6+28.3+11.21.82.21.6
CFD (steady)+56.6+28.5+9.71.10.40.2
CFD (transient)+48.2+30.0+10.41.30.60.5
λExp.+12.4+13.4+3.20.7+0.3+1.3
CFD (steady)+36.0+22.8+8.00.60.10
Ψ/λExp.+14.3+13.1+7.81.12.62.8
CFD (steady)+15.1+4.6+1.50.50.30.2

Obviously, the ASF caused a slight degradation in 0.8–1.0Φd, and it could be more stressed at 0.8Φd from the CFD results, i.e., the ASF might be detrimental to performance near the peak point where the positive slope began to be formed in the case of none; since the peak point implies the minimum limit for the stable operation (or incidence angle) of an axial fan, it seemed reasonable that the performance was affected by ASF near this point. Nevertheless, it was not at a level of concern compared to the huge earnings on stall margin. The degradation in 0.9–1.0Φd seemed to be the blockage effect concerning the thickness of ASF, and it may be advantageous to design the thickness as thin as possible to minimize the blockage effect.

Meanwhile, especially for the case of none, which obtained the pressure fluctuations with the positive slope, the expanded uncertainty for measured pressure during the experimental tests was the greatest at 0.3Φd for the whole flowrate range, but it was only ±0.55%. In the stalling flow rates (0.4, 0.5, 0.6, and 0.7Φd), it was ±0.43, 0.24, 0.24 and 0.20%, respectively. Here, the evaluation of extended uncertainty was performed in a previous study [33].

4.3 Internal Flow Field (Backflow and Rotating Stall).

First, Fig. 8 depicts circumferentially averaged turbulence kinetic energy (k) contours and circumferentially projected three-dimensional streamlines on the meridional plane: the x-axis represents the axial coordinate normalized to the fan diameter (D2); the y-axis represents the normalized span (r*); the white-dotted perpendicular line near the blade LE indicates a mark for descriptions below. In the case of none, a backflow region was identified near the shroud in 0.5–0.7Φd, where the positive slope was formed, and it gradually increased in the spanwise direction toward the hub and the streamwise direction to the inlet passage (upstream) as the flowrate decreased. On the other hand, in the case of ASF attached, most of the backflow was suppressed. Here, the maximum length of backflow for each direction was quantified as presented in Fig. 9(a). At 0.5, 0.6, and 0.7Φd, the axial length (ζa) decreased by 72.8, 69.4, and 58.7% with ASF, and the radial length (ζr) decreased by 62.9, 75.7, and 73.8% with ASF, respectively. Figure 9(b) shows the maximum length of backflow only for the radial direction under the concept of normalized span. Meanwhile, in Figs. 9(a) and 9(b), it was confirmed that the ASF caused some backflow at 0.8Φd; this could be the reason why the case of ASF attached showed a lower pressure rise at 0.8Φd than that of none.

Fig. 8
Turbulence kinetic energy contour and streamline averaged and projected on the meridional plane for each case of (a) none and (b) ASF attached
Fig. 8
Turbulence kinetic energy contour and streamline averaged and projected on the meridional plane for each case of (a) none and (b) ASF attached
Close modal
Fig. 9
Maximum length of backflow for (a) each direction and (b) radial direction only: Axial length (ζa) and Radial length (ζr)
Fig. 9
Maximum length of backflow for (a) each direction and (b) radial direction only: Axial length (ζa) and Radial length (ζr)
Close modal
Fig. 10
Particle tracking for each case of (a) none and (b) ASF attached at 0.5Φd
Fig. 10
Particle tracking for each case of (a) none and (b) ASF attached at 0.5Φd
Close modal
Fig. 11
Circumferential (left) and axial (right) velocity contour, Q-criterion (left), and streamlines (left and right) on sectional plane for case of none: (a) 0.7Φd, (b) 0.6Φd, and (c) 0.5Φd
Fig. 11
Circumferential (left) and axial (right) velocity contour, Q-criterion (left), and streamlines (left and right) on sectional plane for case of none: (a) 0.7Φd, (b) 0.6Φd, and (c) 0.5Φd
Close modal

The stalling and antistalling mechanisms were visualized with Fig. 10; a single particle was seeded near the shroud to track over the whole flow passage. Under the stalling condition, the particle entering the blade (2) was distorted while receiving the same directionality as the blade's rotational direction. In short course, it entered the blade at a span much lower than the shroud; this was because of the occupancy of backflow, which was shown in Fig. 8(a) and would be described later. The particle could not face downstream and headed to the inlet passage as a backflow (3). Subsequently, it showed circling movement along the shroud wall (4); at this time, the circling direction was the same as the blade's rotational direction. The circling particle tried to enter the blade again; however, it became a backflow (5) as in pattern (3). The particle repeated circling movement along the shroud wall (6–7). The ASF blocked this unfavorable flow pattern. The particle entering the blade could be straight without distortion, and the backflow lost its circumferential velocity component, close to zero. The ASF caught the circling movement of backflow and forced it back into the blade.

Fig. 12
Circumferential (left) and axial (right) velocity contour, Q-criterion (left), and streamlines (left and right) on sectional plane for case of ASF attached: (a) 0.7Φd, (b) 0.6Φd, and (c) 0.5Φd
Fig. 12
Circumferential (left) and axial (right) velocity contour, Q-criterion (left), and streamlines (left and right) on sectional plane for case of ASF attached: (a) 0.7Φd, (b) 0.6Φd, and (c) 0.5Φd
Close modal

The flow pattern near the blade LE was further observed from the transient analysis in Figs. 11 and 12: the sectional plane is located at the coordinates corresponding to the white-dotted perpendicular line in Fig. 8; inside each box for a time-step, the left and right one show the circumferential and axial velocity contour; the streamlines are the same for both left and right one; the left one includes the inlet passage with transparency, the blades (black), and a vortex identification method with isosurface (Q-criterion) [45]; the right one contains only an arrow pointing to the center of a certain vortex core identified from the streamlines; the circumferential velocity contour indicates a higher value as it becomes stronger against the blade's rotational direction, and the axial velocity contour indicates a higher value as the component toward the downstream becomes stronger. First, in the case of none at 0.7Φd, the backflow and rotating stall did not fully occupy the shroud wall in the flow passage. There were two cores of rotating stall in accordance with the blade's rotational direction. As expected from Fig. 7, backflow and rotating stall near the blade LE appeared to be in the quickening period at 0.7Φd. At 0.6Φd, backflow and rotating stall entered the growing period. It was confirmed that two cores were strongly developed to be the farthest from each other, and another two cores repeated formation and extinction with weak intensity, i.e., stall cores were rotating with an annular symmetry (every 90 deg); however, the intensity and rotating speed seemed uneven. At 0.5Φd, it could be understood as the setting period. Although the backflow was distributed thicker from the shroud wall, the imaginary line connecting the center of stall core was nearly pentagonal. Moreover, the rotating speed seemed to be uniform. As a qualitative analysis, the rotating speed for each type of rotating stall was approximately 1/3 or less to the blade's rotational speed, and it gradually increased as the flowrate decreased, i.e., at 0.7, 0.6, and 0.5Φd, 85/360  0.23, 103/360  0.28, and 112.5/360  0.31, respectively. Meanwhile, all the stall cores mentioned above had their swirling directionality even while rotating in the flow passage, which was counterclockwise based on the captured figure. On the other hand, in the case of ASF attached, only small eddies due to the filtered flow were confirmed on each pressure side of the ASF. These eddies became progressively active at a lower flowrate; however, the ASF still controlled them. The directionality inside each eddy was quite complicated, and it did not rotate along the circumferential direction. The ASF remarkably stabilized the velocity contour near the blade LE, which was distorted with the backflow and rotating stall in the case of none.

4.4 Pressure Fluctuation and Fast Fourier Transform.

The fluctuations with respect to the pressure rise (differential pressure between the inlet and outlet) in Fig. 7(a) were further indicated in Fig. 13; the y-axis represents the difference of maximum-minimum level (ΔΨ*, bar symbol in Fig. 7(a)). In the case of none, the pressure fluctuation increased sharply in the flowrate range lower than 0.8Φd (0.5–0.7Φd). However, it was slightly decreased at 0.5Φd: at the flowrate point where backflow would be actively developed, the magnitude of pressure fluctuation could be subsided [46]; the tendency in Fig. 11 could be referred in this case. On the other hand, the pressure fluctuations were thoroughly suppressed with the ASF. The ASF eventually secured stability in the whole flow passage of the axial fan.

Fig. 13
Difference of maximum-minimum level of total pressure fluctuation on performance curve
Fig. 13
Difference of maximum-minimum level of total pressure fluctuation on performance curve
Close modal

The stability near the blade LE was evaluated from the FFT results performed with total pressure fluctuation at the four monitoring points (P1–4) in Fig. 5(a): the points were located every 90 deg slightly inward from the shroud wall (0.001D2 far), and the axial coordinates were the same as the white-dotted perpendicular line in Fig. 8; the data extracted here is time-traced total pressure during approximately one revolution of rotating stall's core, here, since stalls are generally unstable factors, the characteristics of each revolution may be different, but data from one revolution was adequate to detect stalls; although the static wall pressure can evaluate the stability of the internal flow, it is considered that the characteristics of the rotating stall cannot be clearly captured because it does not include dynamic pressure. Figure 14 contains the FFT results at each monitoring point. Since the blade passing frequency (BPF) of this axial fan is 245 Hz, it could be confirmed that the ASF completely suppressed instability in the stalling flow rates (0.5–0.7Φd). Even the small eddies in Fig. 12 could not contain any meaningful peaks, and this could be additional proof that the eddies were not rotating. On the other hand, the case of none included quite complicated peaks at lower frequencies than the first BPF. In particular, based on the qualitative analysis of the rotating speed and the number of cores for a rotating stall, its passing frequencies could be inferred as 0.02, 0.11, and 0.15fn at 0.7, 0.6, and 0.5Φd, respectively, where fn denote the BPF. This inference was based on the following two factors: first, the rotating speed was uniform for each core at 0.7, 0.6, and 0.5Φd; second, the number of cores was 1, 4, and 5 at 0.7, 0.6, and 0.5Φd, respectively; e.g., at 0.5Φd, the stall core's rotating speed was estimated as 112.5/360  0.31 times to the blade's rotational speed, then the stall core's passing frequency could be [5×(1470 × 0.31)]/60  37.98 Hz, thus 0.15fn (37.98/245) was indicated in Fig. 14(a). Here, the number of cores at 0.7Φd (Fig. 11(a)) appeared to be two; however, it was not annularly symmetric and was adjacent to each other. Thus, it should be understood as a combined component in terms of the frequency domain. Accordingly, as shown in Fig. 14(a), the unfavorable peaks appeared at 0.02, 0.07–0.12, and 0.15fn for 0.7, 0.6, and 0.5Φd, respectively. It can prove that these peaks were directly originated with the rotating stall. The reason that the peaks at 0.6Φd were confirmed at rather sporadic frequencies for each monitoring point is because each core's rotating speed was not uniformly maintained. At 0.5Φd, the pattern for each monitoring point was almost identical, meaning that the rotating stall had a fairly uniform characteristic. Meanwhile, the magnitude of fluctuations was the highest at 0.7Φd, and those of 0.5Φd and 0.6Φd were similar.

Fig. 14
FFT results at each monitoring point for each case of (a) none and (b) ASF attached
Fig. 14
FFT results at each monitoring point for each case of (a) none and (b) ASF attached
Close modal

5 Conclusion

The ASF was designed and optimized using the 2k full factorial design and response surface method. In terms of the objective function, the axial gap was the most sensitive variable and should be less than 4% of the fan diameter to avoid a positive slope. The radial and axial length could be designed at least 1 and 5% to the fan diameter. The optimized dimensions for radial length, axial length, and axial gap were 3.5, 15, and 1%, respectively, based on the fan diameter. An ASF-attached axial fan stably recovered performance degradation in the stalling flow rates and allowed to form a negative slope to 0.5Φd. From the transient simulations, the total pressure in 0.5–0.7Φd increased by +48.2, +30.0, and +10.4%, respectively. Moreover, the ASF had little effect on the performance near the design flowrate. In the case of an axial fan without the ASF, the backflow gradually increased in the spanwise direction toward the hub and the streamwise direction to the inlet passage (upstream), as the flowrate decreased. The rotating stall was developed through the process of quickening-growing-setting period and was rotating faster, as the flowrate decreased. The ASF blocked the circling backflow near the blade LE and forced it back into the blade. The backflow lost its circumferential velocity component, close to zero. The flow entering the blade could be straight without distortion. The rotating stall could not be formed with the ASF. The ASF resolved the instability inside the flow passage, which was contained in the stalling flow rates. In 0.5–0.7Φd, the total pressure fluctuation was reduced by 92.6, 95.6, and 67.2%, respectively. It did not allow instability even near the blade LE.

As an additional consideration, the optimized dimensions for each variable were limited to each low-high level. From the results of this study, if the variables were considered at a wider low-high level, the radial and axial length could be longer than 3.5 and 15%, and the axial gap could be narrower than 1% to the fan diameter, respectively. However, the sufficiently secured radial and axial length had little effect on the objective function. From a practical point of view, the authors tried to avoid presenting results for the range where the effect of radial and axial length on the objective function was insignificant. Since the axial gap could not be zero, 1% to the fan diameter was selected as the low level, which could be an intuitive value with a nonzero. The authors expect the selection for each low-high level to be understood as an implicit constraint.

Acknowledgment

The antistalling method of this study was recognized for its originality from a patent application (1071231, 10-2022-0099333) and the patent contained most of the results.

Funding Data

  • Korea Institute of Energy Technology Evaluation and Planning (KETEP), Ministry of Trade, Industry and Energy (MOTIE) (Award No. 2021202080026D; Funder ID: 10.13039/501100003052).

Conflict of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Data Availability Statement

The data that support the findings of this study are available from the first and corresponding author, Y.-I.K. and Y.-S.C., upon reasonable request.

Nomenclature

a =

objective function (gradient on Q-P curve)

A2 =

fan outlet area (πrs2πrh2), m2

ANSI =

american national standards institute

ASF =

anti-stall fin

BPF =

blade passing frequency, 245 hertz

c =

absolute velocity, m/s

C =

blade chord length, m

CFD =

computational fluid dynamics

D2 =

fan diameter, m

Exp. =

experimental test

Fi =

body force, kg·m/s2

fn =

blade passing frequency, Hz

FFT =

fast Fourier transform

GCI =

grid convergence index

i =

incidence angle, deg

I =

electric current, A

k =

turbulence kinetic energy, m2/s2

la =

axial length of ASF, m

lr =

radial length of ASF, m

lhexa =

maximum length of a hexahedral grid cell, m

ltetra =

maximum length of a tetrahedral grid cell, m

L =

shaft power, W

LE =

leading edge

n =

number of grid cells

N =

rotational speed, rpm

Ns =

specific speed (type number)

NACA =

national advisory committee for aeronautics

P =

total pressure, Pa

Q =

volume flow rate, m3/s or invariant of velocity gradient tensor, 1/s2

r =

fan radius, m

r* =

hub-to-shroud normalized span

R2 =

coefficient of determination

RMS =

root-means-square

RSM =

response surface method

S =

inlet blade-to-blade pitch, m

SST RM =

shear stress transport reattachment modification

SST Std =

shear stress transport standard

T =

torque, N·m

Tu =

turbulence intensity

t =

time, s

TE =

trailing edge

u (vθ) =

circumferential velocity, m/s

U (could be substituted

as V or W) =

velocity, m/s

V =

voltage, V

va =

axial velocity, m/s

vθ (u) =

circumferential velocity, m/s

w =

relative velocity, m/s

x (could be substituted as y or z) =

coordinate, m

z =

axial coordinate, m

αgrid =

grid ratio

β =

flow angle, deg

βASF =

axial directionality (angle) of ASF, deg

ΔPs =

differential pressure for nozzles, Pa

ΔΨ* =

difference of maximum-minimum level for pressure fluctuations

δ =

axial gap of ASF, m

δILE =

axial distance between upstream interface and blade LE, m

δTEI =

axial distance between blade TE and downstream interface, m

δt =

tip clearance, m

ζa =

axial length of backflow, m

ζr =

radial length of backflow, m

θASF =

tangential directionality (angle) of ASF, deg

λ =

shaft power coefficient

ρ =

air density, kg/m3

τij =

viscous stress tensor, N/m2

Φ =

flow coefficient

Ψ =

pressure coefficient

ω =

angular velocity, rad/s or turbulence eddy frequency, 1/s=Hz

subscript b =

blade

subscript h =

hub

subscript m =

meridional component

subscript s =

shroud or static

subscript x1 =

multiple of normalized flow rate based on the design flow rate, 0.7

subscript x2 =

multiple of normalized flow rate based on the design flow rate, 0.8

subscript 1 =

blade inlet

subscript 2 =

blade outlet

Appendix

The inlet passage was composed of a tetrahedral type, and the rotating and outlet passages were filled with a hexahedral type (see the enlarged window in Fig. 5(a)). Here, different grid types (tetra and hexa) were in contact at the interface between the inlet-rotating passage. Additionally, in this study, since ASF would approach the blade LE, a relatively deeper consideration was required to secure the grid quality, i.e., as shown in the enlarged window of Fig. 5(a), the grid ratio was difficult to be properly secured because it was not evenly spread out: the grid system in this study was under the automatic generation; a local refinement for grid system interferes with consistency and is not efficient for repetitive tasks. Therefore, it was necessary to secure the optimal condition by adjusting the grid size or axial coordinates of the interface. Figure 15 shows the total pressure rise for each grid ratio (αgrid) between the inlet-rotating passage
αgrid=1lhexaltetra,(ltetra>lhexa)
(A1)
αgrid=ltetralhexa1,(ltetra<lhexa)
(A2)
αgrid=0,(ltetra=lhexa)
(A3)

where ltetra and lhexa denote the maximum length that one tetrahedral and hexahedral cell can have in each direction (circumferential, streamwise, or spanwise). The grid system was reviewed with a total of nine sets from G1-I1 to G7-I6: the review was performed at the design flowrate for the case of none; the abbreviation G and I denote grid and interface, respectively; the apostrophe implies it had only adjusted in spanwise direction; the axial distance for upstream interface-blade LE (δILE) and blade TE-downstream interface (δTEI) is listed in Table 5; the total pressure rise was normalized based on G7-I6. From the results, G7-I6 was predicted as the highest rise among the reviewed grid systems, which showed a deviation of only +0.47% compared to the experimental test. At the interface, the grid ratio in each direction (circumferential, streamwise, and spanwise) appeared to have a minimal effect when it was within approximately 0.8 (1:5 based on the number of cells, n): this was consistent with a general recommendation [37]; the grid ratios outside the recommendation led to low predictions. For example, the increase from G2 s to G3 s and from G5-I4 to G6-I5 could be obtained because δILE in shroud and hub were adjusted to improve the grid ratio in each streamwise and circumferential direction. On the other hand, the increase from G3 s to G4-I3 seemed to be due to the adjustment of δTEI in shroud; however, it could not be analyzed as a specific reason because Fig. 15 contained interactions between the grid size and axial coordinates of the interface. Among the sets that maintained the proper grid ratio, the interface condition that allowed the ASF to be attached in the inlet passage was I5, thus the interface corresponding to I5 was combined with the grid system, satisfying the grid ratio within 0.8. Meanwhile, the results in Fig. 15 are just a process and reference, not absolute.

Table 5

Axial distance of upstream interface-blade LE (δILE) and blade TE-downstream interface (δTEI) normalized by fan diameter

δILE/D2δTEI/D2
HubShroudHubShroud
I10.0030.0030.0070.048
I20.0030.0080.0070.048
I30.0030.0080.0070.008
I40.0030.0080.0070.032
I50.0060.0070.0230.048
I60.0130.0240.0070.048
δILE/D2δTEI/D2
HubShroudHubShroud
I10.0030.0030.0070.048
I20.0030.0080.0070.048
I30.0030.0080.0070.008
I40.0030.0080.0070.032
I50.0060.0070.0230.048
I60.0130.0240.0070.048
Fig. 15
Total pressure rise for each grid ratio in circumferential, streamwise, and spanwise direction
Fig. 15
Total pressure rise for each grid ratio in circumferential, streamwise, and spanwise direction
Close modal

References

1.
Bleier
,
F. P.
,
1998
,
Fan Handbook
,
McGraw-Hill
,
New York
.
2.
Eck
,
B.
,
1973
,
Fans
,
Pergamon Press Inc
.,
New York
.
3.
Čudina
,
M.
,
1995
, “
Flow Instabilities and Noise Generation in Vane Axial Fans Due to Rotating Stall and Surge
,”
SAE Trans.: J. Engines
,
104
(
3
), p.
952104
.10.4271/952104
4.
Zhang
,
L.
,
Liang
,
S.
, and
Hu
,
C.
,
2014
, “
Flow and Noise Characteristics of Centrifugal Fan Under Different Stall Conditions
,”
Math. Probl. Eng.
,
2014
(
1
), pp.
1
9
.10.1155/2014/403541
5.
Rendu
,
Q.
,
Vahdati
,
M.
, and
Salles
,
L.
,
2019
, “
Radial Decomposition of Blade Vibration to Identify a Stall Flutter Source in a Transonic Fan
,”
ASME J. Turbomach.
,
141
(
10
), p.
101011
.10.1115/1.4044484
6.
Zhao
,
X.
,
Zhou
,
Q.
,
Yang
,
S.
, and
Li
,
H.
,
2019
, “
Rotating Stall Induced Non-Synchronous Blade Vibration Analysis for an Unshrouded Industrial Centrifugal Compressor
,”
Sensors
,
19
(
22
), p.
4995
.10.3390/s19224995
7.
Shahriyari
,
M. J.
,
Khaleghi
,
H.
, and
Heinrich
,
M.
,
2019
, “
A Model for Stall and Surge in Low-Speed Contra-Rotating Fans
,”
ASME J. Eng. Gas Turbines Power
,
141
(
8
), p.
081009
.10.1115/1.4043251
8.
Sheard
,
A. G.
, and
Corsini
,
A.
,
2012
, “
The Mechanical Impact of Aerodynamic Stall on Tunnel Ventilation Fans
,”
Int. J. Rotating Mach.
,
2012
(
2
), pp.
1
12
.10.1155/2012/402763
9.
Zhang
,
L.
,
Zhang
,
Q.
, and
Hu
,
C.
,
2013
, “
Fluid Dynamics Characteristics During Stall Inception Process in a Centrifugal Fan
,”
Res. J. Appl. Sci., Eng. Technol.
,
6
(
6
), pp.
1059
1064
.10.19026/rjaset.6.4013
10.
Fike
,
M.
,
Bombek
,
G.
,
Hriberšek
,
M.
, and
Hribernik
,
A.
,
2014
, “
Visualisation of Rotating Stall in an Axial Flow Fan
,”
Exp. Therm. Fluid Sci.
,
53
, pp.
269
276
.10.1016/j.expthermflusci.2013.12.020
11.
Liang
,
F.
,
Xie
,
Z.
,
Xia
,
Z.
, and
Zhou
,
M.
,
2020
, “
Aeroelastic Simulation of the First 1.5-Stage Aeroengine Fan at Rotating Stall
,”
Chin. J. Aeronaut.
,
33
(
2
), pp.
529
549
.10.1016/j.cja.2019.05.004
12.
Zhang
,
L.
,
Zhang
,
L.
,
Zhang
,
Q.
,
Jiang
,
K.
,
Tie
,
Y.
, and
Wang
,
S.
,
2018
, “
Effects of the Second-Stage of Rotor With Single Abnormal Blade Angle on Rotating Stall of a Two-Stage Variable Pitch Axial Fan
,”
Energies
,
11
(
12
), p.
3293
.10.3390/en11123293
13.
Lee
,
S. G.
,
Kim
,
Y. I.
,
Yang
,
H. M.
,
Kim
,
S.
,
Lee
,
S. Y.
,
Lee
,
K. Y.
,
Yang
,
S. H.
, and
Choi
,
Y. S.
,
2020
, “
Numerical Study on the Effect of the Tip Clearance of a 100 HP Axial Fan on the Aerodynamic Performance and Unsteady Stall Characteristics
,”
J. Mech. Sci. Technol.
,
34
(
12
), pp.
5117
5137
.10.1007/s12206-020-1115-2
14.
Nishioka
,
T.
,
2013
, “
Inception Mechanism and Suppression of Rotating Stall in an Axial-Flow Fan
,”
IOP Conf. Ser.: Mater. Sci. Eng.
,
52
(
1
), p.
012002
.10.1088/1757-899X/52/1/012002
15.
Zhang
,
X.
,
Jia
,
X.
, and
Jiang
,
X.
,
2021
, “
The Effect of Load Control on the Performance of Contra-Rotating Fans
,”
Processes
,
9
(
7
), p.
1227
.10.3390/pr9071227
16.
Wang
,
J.
, and
Kruyt
,
N. P.
,
2020
, “
Computational Fluid Dynamics Simulations of Aerodynamic Performance of Low-Pressure Axial Fans With Small Hub-to-Tip Diameter Ratio
,”
ASME J. Fluids Eng.
,
142
(
9
), p.
091202
.10.1115/1.4047120
17.
Gad-el-Hak
,
M.
, and
Bushnell
,
D. M.
,
1991
, “
Separation Control: Review
,”
ASME J. Fluids Eng.
,
113
(
1
), pp.
5
30
.10.1115/1.2926497
18.
Lim
,
B. J.
,
Park
,
T. C.
, and
Kang
,
Y. S.
,
2014
,
An Apparatus for Stall Suppression of an Axial Compressor by Casing Treatment
,
MI Patent & Law Firm, Patent No. KR 1387391 B1
.
19.
Tan
,
C. S.
,
Day
,
I.
,
Morris
,
S.
, and
Wadia
,
A.
,
2010
, “
Spike-Type Compressor Stall Inception, Detection, and Control
,”
Annu. Rev. Fluid Mech.
,
42
(
1
), pp.
275
300
.10.1146/annurev-fluid-121108-145603
20.
Gentry
,
T.
,
Reynolds
,
B. D.
, and
Nolcheff
,
N.
,
2020
,
Grooved Shroud Casing Treatment for High Pressure Compressor in a Turbine Engine
,
Honeywell International Inc
., Patent No. US 20,200,182,259 A1.
21.
Jung
,
J. H.
, and
Kim
,
K. H.
,
2014
,
Structure for Preventing Stall of Axial Fan
, Patent No. KR 1422765 B1.
22.
McKelvey
,
W. G.
,
2001
,
Axial Fan Skip-Stall
,
Allied Signal Inc.
, Patent No. US 6,302,640 B1.
23.
Bianchi
,
S.
,
Corsini
,
A.
,
Sheard
,
A. G.
, and
Tortora
,
C.
,
2013
, “
A Critical Review of Stall Control Techniques in Industrial Fans
,”
ISRN Mech. Eng.
,
2013
(
2
), pp.
1
18
.10.1155/2013/526192
24.
Dygert
,
R. K.
, and
Bushnell
,
P. R.
,
2018
,
Stall Margin Enhancement of Axial Fan With Rotating Shroud
,
Carrier Corporation
, Patent No. US 9,885,368 B2.
25.
Khaleghi
,
H.
,
2015
, “
Stall Inception and Control in a Transonic Fan, Part B: Stall Control by Discrete Endwall Injection
,”
Aerosp. Sci. Technol.
,
41
, pp.
151
157
.10.1016/j.ast.2014.12.022
26.
Salunkhe
,
P. B.
,
Joseph
,
J.
, and
Pradeep
,
A. M.
,
2011
, “
Active Feedback Control of Stall in an Axial Flow Fan Under Dynamic Inflow Distortion
,”
Exp. Therm. Fluid Sci.
,
35
(
6
), pp.
1135
1142
.10.1016/j.expthermflusci.2011.03.008
27.
Jung
,
Y. J.
,
Jeon
,
H. S.
,
Jung
,
Y. H.
,
Lee
,
K. J.
, and
Choi
,
M. S.
,
2016
, “
Effects of Recessed Blade Tips on Stall Margin in a Transonic Axial Compressor
,”
Aerosp. Sci. Technol.
,
54
, pp.
41
48
.10.1016/j.ast.2016.04.009
28.
Cyrus
,
V.
,
Trávníček
,
Z.
,
Wurst
,
P.
, and
Kordík
,
J.
,
2013
, “
Flow Control in Axial Fan Inlet Guide Vanes by Synthetic Jets
,”
EPJ Web Conf.
,
45
, p.
01021
.10.1051/epjconf/20134501021
29.
Shrestha
,
B.
, and
Bhattarai
,
B.
,
2021
, “
Numerical Simulation Studies on Stall Suppression of a NACA0015 Airfoil
,”
J. Adv. Coll. Eng. Manage.
,
6
, pp.
23
31
.10.3126/jacem.v6i0.38276
30.
Peng
,
G.
,
Du
,
J.
,
Chang
,
H.
,
Chen
,
Q.
,
Li
,
J.
,
Pan
,
C.
, and
Yang
,
Y.
,
2021
, “
Numerical and Experimental Analysis of Influence of Impeller Structures on Slurry Pump Performance
,”
Front. Energy Res.
,
9
, p.
762159
.10.3389/fenrg.2021.762159
31.
Corsini
,
A.
,
Delibra
,
G.
, and
Sheard
,
A. G.
,
2013
, “
On the Role of Leading-Edge Bumps in the Control of Stall Onset in Axial Fan Blades
,”
ASME J. Fluids Eng.
,
135
(
8
), p.
081104
.10.1115/1.4024115
32.
Choi
,
Y. S.
,
Kim
,
Y. I.
,
Lee
,
K. Y.
, and
Yang
,
H. M.
,
2022
,
Anti-Stall Structure and AXIAL FAN Including Same
(1071231, 10-2022-0099333),
Eruum & Leeon Intellectual Property Law Firm
.
33.
Kim
,
Y. I.
,
Lee
,
S. Y.
,
Lee
,
K. Y.
,
Yang
,
S. H.
, and
Choi
,
Y. S.
,
2020
, “
Numerical Investigation of Performance and Flow Characteristics of a Tunnel Ventilation Axial Fan With Thickness Profile Treatments of NACA Airfoil
,”
Energies
,
13
(
21
), p.
5831
.10.3390/en13215831
34.
Air Movement and Control Association International, Inc.
,
2008
,
ANSI/AMCA 210-07: Laboratory Methods of Testing Fans for Certified Aerodynamic Performance Rating
,
Air Movement and Control Association International, Inc. and the American Society of Heating, Refrigerating, and Air Conditioning Engineers
, Atlanta, GA.
35.
Choi
,
Y. S.
,
Kim
,
D. S.
, and
Yoon
,
J. Y.
,
2005
, “
Effects of Flow Settling Means on the Performance of Fan Tester
,”
KSFM J. Fluid Mach.
,
8
(
5
), pp.
29
34
.10.5293/KFMA.2005.8.5.029
36.
Menter
,
F.
,
1994
, “
Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications
,”
AIAA J.
,
32
(
8
), pp.
1598
1605
.10.2514/3.12149
37.
ANSYS Inc.
,
2016
,
ANSYS User Manual
,
ANSYS Inc
., Canonsburg, PA.
38.
Brown
,
G. J.
,
Fletcher
,
D. F.
,
Leggoe
,
J. W.
, and
Whyte
,
D. S.
,
2018
, “
Investigation of Turbulence Model Selection on the Predicted Flow Behaviour in an Industrial Crystalliser—RANS and URANS Approaches
,”
Chem. Eng. Res. Des.
,
140
, pp.
205
220
.10.1016/j.cherd.2018.10.007
39.
Liu
,
B.
,
An
,
G.
, and
Yu
,
X.
,
2015
, “
Assessment of Curvature Correction and Reattachment Modification Into the Shear Stress Transport Model Within the Subsonic Axial Compressor Simulations
,”
Proc. Inst. Mech. Eng., Part A
,
229
(
8
), pp.
910
927
.10.1177/0957650915595168
40.
Kim
,
Y. I.
,
Choi
,
Y. U.
,
Jeong
,
C. Y.
,
Lee
,
K. Y.
, and
Choi
,
Y. S.
,
2021
, “
Effect of Two- and Three-Dimensionally Designed Guide Vanes With Different Camber Length on Static Pressure Recovery of a Wall-Mounted Axial Fan
,”
Processes
,
9
(
9
), p.
1595
.10.3390/pr9091595
41.
Lee
,
S. G.
,
Lee
,
K. Y.
,
Yang
,
S. H.
, and
Choi
,
Y. S.
,
2019
, “
A Study on Performance Characteristics of an Axial Fan With a Geometrical Parameters of Inlet Hub Cap
,”
KSFM J. Fluid Mach.
,
22
(
4
), pp.
5
12
.10.5293/kfma.2019.22.4.005
42.
Roache
,
P. J.
,
1994
, “
Perspective: A Method for Uniform Reporting of Grid Refinement Studies
,”
ASME J. Fluids Eng.
,
116
(
3
), pp.
405
413
.10.1115/1.2910291
43.
Celik
,
I. B.
,
Ghia
,
U.
,
Roache
,
P. J.
,
Freitas
,
C. J.
,
Coleman
,
H.
, and
Raad
,
P. E.
,
2008
, “
Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications
,”
ASME J. Fluids Eng.
,
130
(
7
), p.
078001
.10.1115/1.2960953
44.
Wang
,
G. G.
,
Dong
,
Z.
, and
Aitchison
,
P.
,
2001
, “
Adaptive Response Surface Method – a Global Optimization Scheme for Approximation-Based Design Problems
,”
Eng. Optim.
,
33
(
6
), pp.
707
733
.10.1080/03052150108940940
45.
Hunt
,
J. C. R.
,
Wray
,
A. A.
, and
Moin
,
P.
,
1988
, “
Eddies, Stream, and Convergence Zones in Turbulent Flows
,” Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1988 Summer Program,
NASA Ames Research Center, Mountain View, CA, pp. 193–208
.
46.
Zhou
,
B.
,
He
,
X.
,
Yang
,
H.
,
Zhu
,
Z.
,
Wei
,
Y.
, and
Zhang
,
Y.
,
2020
, “
Unsteady Flow Characteristics of Rotating Stall and Surging in a Backward Centrifugal Fan at Low Flow-Rate Conditions
,”
Processes
,
8
(
7
), p.
872
.10.3390/pr8070872