A Study on the Half-Ducted Axial Flow Fan Designed by a Diagonal Flow Fan Design Method

A study on the half-ducted axial flow fan designed by a diagonal flow fan design method was conducted. The rotor which has NACA65 blades was designed, calculated numerically, manufactured and tested experimentally. As a result of the design and CFD, the meridional streamline and three distributions of the meridional, tangential and radial velocity at inlet and outlet go well as designed values of the half ducted fan. On the other hand, the values of the meridional velocity and the tangential velocity are little smaller than the design values at the hub side of the radial distribution. The improvement of the design is prospected for this point, that is, the approach between the design value and the actual flow is prospected if the tangential velocity is assigned small at hub and is assigned large at the tip so as to accord the actual flow in the vortex design of the rotor blade. Then the designed half-ducted rotor with four NACA65 blades was fabricated by a three-dimensional printer and tested in the wind tunnel in order to validate the half-ducted design method. For the comparison between the design values and the experimental values at the design flow rate coefficient of φ = 0.264, the experimental values of the pressure rise coefficient ψ and the efficiency η are rather small than the design values, while the experimental value of the torque coefficient τ is almost the same as the design value. However, the experimental value of approximately 0.45 of the maximum efficiency is comparably large value considering for the limitation of the situation of half-ducted. For the comparison between the experimental values and the CFD values at φ = 0.264, the CFD values are almost the same values as the experimental values for all the values of ψ, τ and η. In addition, the tendencies of the CFD values when the flow rate coefficient changes are almost similar as the experimental tendencies, though the flow rate coefficient for the CFD values when ψ or η takes the peak value shifts toward larger flow rate. For the case at rotor outlet at φ = 0.264, two How to cite this paper: Kaji, A., Kinoue, Y., Shiomi, N. and Setoguchi, T. (2017) A Study on the Half-Ducted Axial Flow Fan Designed by a Diagonal Flow Fan Design Method. Open Journal of Fluid Dynamics, 7, 15-25. https://doi.org/10.4236/ojfd.2017.71002 Received: December 20, 2016 Accepted: January 13, 2017 Published: January 16, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access

values at the design flow rate coefficient of φ = 0.264, the experimental values of the pressure rise coefficient ψ and the efficiency η are rather small than the design values, while the experimental value of the torque coefficient τ is almost the same as the design value.However, the experimental value of approximately 0.45 of the maximum efficiency is comparably large value considering for the limitation of the situation of half-ducted.For the comparison between the experimental values and the CFD values at φ = 0.264, the CFD values are almost the same values as the experimental values for all the values of ψ, τ and η.In addition, the tendencies of the CFD values when the flow rate coefficient changes are almost similar as the experimental tendencies, though the flow rate coefficient for the CFD values when ψ or η takes the peak value shifts toward larger flow rate.For the case at rotor outlet at φ = 0.264, two values of the meridional velocity and the tangential velocity are larger than the design values at the tip side of the radial distribution.

Introduction
A lot of axial fans, which are small size and high efficiency, are used in our daily life, such as a power unit cooling fan of personal computer, a room ventilation fan, and a radiator fan in car engine room.In most applications of fans nonducted types are used for the space limitation or cost saving.In comparison with the conventional full-ducted fan mainly used in industrial applications, they are called as half-ducted type, semi-opened type and opened type, respectively.Among these types, a half-ducted axial flow fan is focused in this paper.The corresponding design methods for half-ducted types of propeller fans are not studied by many researchers.Recently, many of fans are designed by the inverse design method [1] [2] [3].A half-ducted axial flow fan is usually designed provided that both the inflow and the outflow are full-ducted as if it were in the straight pipe.However, in our previous result [4] [5], a half-ducted axial flow fan has much more complicated and three-dimensional internal flow than a fullducted fan.All conventional design methods are not taking the radial velocity component into account.In this paper, the radial inclinations of meridional streamline were considered for the improvement of the design method by using the diagonal flow fan design method.The half-ducted axial flow fans are designed by the controlled vortex design by specifying the constant tangential velocity both at inlet and outlet of the rotor.Conventional half-ducted fan is usually designed to prescribe almost uniform inflow and outflow as if it was in the straight pipe.However, many axial flow fans are not used in the straight pipe, such as the application of using in the ventilation and cooling systems without pipe.Thus it is important to take the real flow situation into account in design.Therefore, half-ducted axial flow fan was designed to compare with the traditional design of ducted ones by specifying the flow angles according to the previous experimental results of authors [4] [5].

Design Method
The quasi three-dimensional flow theory was applied to investigate the flow of the axial flow fans.The meridional flow and the revolutional flow between blades were calculated by the method of streamline curvature.Based on the theory, the meridional flow was calculated by adopting the radial balance equations [6], while the calculation of the blade to blade flow was obtained by 2D cascade data with the correction by a potential flow theory so as to consider the axial flow velocity change and the inclination of the flow surface [7].In calculation of meridional flow, the radial balance equation was evaluated.Pressure gradient from hub to casing occurs to be balanced against centrifugal force by swirled blades.The radial direction distribution of state quantity such as flow velocity is dependent on such the balancing condition, which is called radial equilibrium condition.The following radial balance equation was evaluated at the quasi-orthogonal direction q on meridional plane.
( ) ( ) When the compressibility of the fluid is ignored, ( ) ( ) where, m C : meridional velocity [m/s], C θ : tangential velocity [m/s], t P : total The solution of Equation ( 1) is given as where, C is an integral constant, which satisfies the following equation.

cos d
where, G: designed flow rate [kg/s], K B : end-wall blockage coefficient.K B was presumed to be 0.96.The total pressure rise is presumed to be able to calculate by the Euler equation as following.
( ) The controlled vortex design method has been applied for half-ducted and ducted propeller fan design by specifying the constant tangential velocity at both inlet and outlet of the fan rotor.Therefore, the meridional velocity and the tangential velocity can be obtained so that the calculation of meridional flow is finished.Table 1 shows the designed parameters of ducted fan and half-ducted fan.All these parameters are specified the same value for both of the designed fans.4 where t U is speed on the rotor tip.t D and h D are the diameter of the rotor tip and hub respectively.Q is the volume flow rate, s p ∆ is the static pressure rise, ρ is atmospheric density, T is the rotor torque and ω is the rotating speed of the rotor.
The blade profile on the revolutional plane was selected by referring to the diagram of NACA65 carpet.In order to use the cascade data, the through flow in cylindrical or cone surface must be considered on an averaged stream surface.The flow on revolving stream surface is projected to a plane.Thus, the NACA65 blade with quadrilateral blade on the meridional plane was adopted.Threedimensional effect of inclination and thickness variation of stream surface on flow field is substituted by the distributions of vortex and divergence on the potential theory [8].

Design Result
The streamlines on the meridional plane of half-ducted fan were shown in Figure 1.The eleven streamlines from hub to tip were designed so that flow rates between each streamline became equal.In this case, the inclinations of inflow and outflow are given based on the experimental data [4] [5] so as to consider the effect of radial velocity component on the internal flow field.The flow angles of the streamlines are given 63, 16, 38 degrees respectively on tip streamline at inlet, on tip streamline at outlet, on hub streamline at outlet.The blade shapes on top view were obtained showing in Figure 2.

Numerical Method and Result
The analysis of the three-dimensional internal flow fields of half ducted fan are conducted comparing to the numeral computation results in the commercial software (Ansys CFX Release 16.0).Simulation conditions and layouts are shown in Table 2 and Figure 3.
The computation domain was extended 300 mm upstream from the leading edge and 500 mm downstream from the trailing edge, respectively.Figure 3 shows  the perspective view of the calculated region for the full-pitch calculation.The unstructured, tetrahedral grids were used.The total number of the grids is approximately 5,550,000.The tip clearance of 2 mm is set for the numerical calculation.
The turbulence model was SST.The use of a k-ω formulation in the inner parts of the boundary layer makes the model directly usable all the way down to the wall through the viscous sub-layer.The SST formulation also switches to a k-ε behaviour in the free-stream and thereby avoids the common k-ω problem that the model is too sensitive to the inlet free-stream turbulence properties.Especially, the influence of stall by separated flow at low flow rate can be better simulated in SST turbulence model.Therefore SST turbulence model was adopted.from hub to tip was adopted.At rotor inlet, tangential velocity is almost zero, which goes as designed.So it can be argued that applying this design method considering inclination of meridional streamline to half ducted fan is effective.
For the case at rotor outlet, at the flow rate coefficient of φ = 0.264, two values of the meridional velocity and the tangential velocity are larger than the design values at the tip side of the radial distribution.On the other hand, the values of 2 m C and 2 C θ are smaller than the design values at the hub side of the radial distribution.The improvement of the design is prospected for this point, that is, the approach between the design value and the actual flow is prospected if the tangential velocity is assigned small at the hub and is assigned large at the tip so as to accord the actual flow in the vortex design of the rotor blade.

Experiment
The designed and calculated half-ducted rotor with four NACA65 blades was fabricated by a three-dimensional printer shown in Figure 8 and tested in the wind tunnel in order to validate the half-ducted design method.The rotor has been tested in experiment equipment system in Figure 9 with an 880 mm × 880 mm square cross section and a length of 3 m.The rotor speed can arbitrarily set on the range of 0 to 3100 min −1 driven by the direct-current electric motor.In order to get smaller static pressure rise, the centrifugal fans were set in the outlet of the system.The flow rate of the system can be well controlled by altering the  ( ) The static pressure of the propeller fan is obtained by averaging the data from 16 static pressure taps setting at downstream of rotor as show in Figure 10.The flow rate is calculated by the differential pressure between the upstream and downstream of the nozzle where eight pressure taps are set respectively.The higher flow rate, corresponding to larger differential pressure, can be obtained by elevating the rotating speed of two centrifugal booster fans, and the rotor speed is set by 3000 min −1 in this paper.Calculated and measured pressure rise, torque and efficiency is shown in Figure 10.For the comparison between the design values and the experimental values at the design flow rate coefficient of φ = 0.264, the experimental values of the pressure rise coefficient ψ and the efficiency η are rather small than the design values, while the experimental value of the torque coefficient is almost the same as the design value.However, the experimental value of approximately 0.45 of the maximum efficiency is comparably large value considering for the limitation of the situation of half-ducted.For the comparison between the experimental values and the CFD values at the design pressure [Pa], ρ : fluid density [kg/m 3 ], r : radius [m], m r : curvature radius of meridional line [m], ε : angle between meridional normal line and ϕ : inclination angle of revolving stream surface [deg].
The hub ratio has to be set a little larger to avoid the flow focus on the blade tip.Designed pressure rise and efficiency are set 200 [Pa] and 60[%] which are comparably large value than the conventional one for the superior target.The flow rate and pressure rise are represented with no dimensional form of flow coefficient and pressure rise coefficient which are defined as follows.

Figure 4 Figure 3 .
Figure 4 shows the designed and simulated streamlines.The designed and simulated streamlines are similar tendency.Figures 5-7 present the distributions of circumferentially averaged meridional, tangential and radial velocity at rotor inlet and rotor outlet in terms of numerical results and design data, respectively for half-ducted fan along the radial direction.In the figures, the meridional and