<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJFD</journal-id><journal-title-group><journal-title>Open Journal of Fluid Dynamics</journal-title></journal-title-group><issn pub-type="epub">2165-3852</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojfd.2017.71002</article-id><article-id pub-id-type="publisher-id">OJFD-73490</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  A Study on the Half-Ducted Axial Flow Fan Designed by a Diagonal Flow Fan Design Method
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Atsushi</surname><given-names>Kaji</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yoichi</surname><given-names>Kinoue</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Norimasa</surname><given-names>Shiomi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Toshiaki</surname><given-names>Setoguchi</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Institute of Ocean Energy, Saga University, Saga, Japan</addr-line></aff><aff id="aff1"><addr-line>Department of Mechanical Engineering, Saga University, Saga, Japan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>kinoue@me.saga-u.ac.jp(YK)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>27</day><month>12</month><year>2016</year></pub-date><volume>07</volume><issue>01</issue><fpage>15</fpage><lpage>25</lpage><history><date date-type="received"><day>December</day>	<month>20,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>January</month>	<year>13,</year>	</date><date date-type="accepted"><day>January</day>	<month>16,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  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 values of the meridional velocity and the tangential velocity are larger than the design values at the tip side of the radial distribution.
 
</p></abstract><kwd-group><kwd>Axial Flow Fan</kwd><kwd> Half-Ducted Fan</kwd><kwd> Internal Flow</kwd><kwd> Streamline Inclination</kwd><kwd> Numerical Analysis</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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 non- ducted 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 [<xref ref-type="bibr" rid="scirp.73490-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.73490-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.73490-ref3">3</xref>] . 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 [<xref ref-type="bibr" rid="scirp.73490-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.73490-ref5">5</xref>] , a half-ducted axial flow fan has much more complicated and three-dimensional internal flow than a full- ducted 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 [<xref ref-type="bibr" rid="scirp.73490-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.73490-ref5">5</xref>] .</p></sec><sec id="s2"><title>2. Design Method</title><p>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 [<xref ref-type="bibr" rid="scirp.73490-ref6">6</xref>] , 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 [<xref ref-type="bibr" rid="scirp.73490-ref7">7</xref>] . 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.</p><disp-formula id="scirp.73490-formula50"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x2.png"  xlink:type="simple"/></disp-formula><p>When the compressibility of the fluid is ignored,</p><disp-formula id="scirp.73490-formula51"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x3.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.73490-formula52"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x4.png"  xlink:type="simple"/></disp-formula><p>where,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x5.png" xlink:type="simple"/></inline-formula>: meridional velocity [m/s],<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x6.png" xlink:type="simple"/></inline-formula>: tangential velocity [m/s],<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x7.png" xlink:type="simple"/></inline-formula>: total pressure [Pa],<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x8.png" xlink:type="simple"/></inline-formula>: fluid density [kg/m<sup>3</sup>],<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x9.png" xlink:type="simple"/></inline-formula>: radius [m],<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x10.png" xlink:type="simple"/></inline-formula>: curvature radius of meridional line [m],<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x11.png" xlink:type="simple"/></inline-formula>: angle between meridional normal line and</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x12.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x13.png" xlink:type="simple"/></inline-formula>: angle between radial direction and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x14.png" xlink:type="simple"/></inline-formula></p><p>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x15.png" xlink:type="simple"/></inline-formula>: inclination angle of revolving stream surface [deg].</p><p>The solution of Equation (1) is given as</p><disp-formula id="scirp.73490-formula53"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x16.png"  xlink:type="simple"/></disp-formula><p>where, C is an integral constant, which satisfies the following equation.</p><disp-formula id="scirp.73490-formula54"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x17.png"  xlink:type="simple"/></disp-formula><p>where, G: designed flow rate [kg/s], K<sub>B</sub>: end-wall blockage coefficient. K<sub>B</sub> was presumed to be 0.96. The total pressure rise is presumed to be able to calculate by the Euler equation as following.</p><disp-formula id="scirp.73490-formula55"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x18.png"  xlink:type="simple"/></disp-formula><p>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. <xref ref-type="table" rid="table1">Table 1</xref> 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.</p><p>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.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Design parameter</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Tip diameter [mm]</th><th align="center" valign="middle" >200</th></tr></thead><tr><td align="center" valign="middle" >Hub-tip ratio</td><td align="center" valign="middle" >0.6</td></tr><tr><td align="center" valign="middle" >Blade number</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >Flow coefficient</td><td align="center" valign="middle" >0.264</td></tr><tr><td align="center" valign="middle" >Pressure-rise coefficient</td><td align="center" valign="middle" >0.336</td></tr><tr><td align="center" valign="middle" >Rotation speed [min<sup>−1</sup>]</td><td align="center" valign="middle" >3000</td></tr><tr><td align="center" valign="middle" >Pressure rise [Pa]</td><td align="center" valign="middle" >200</td></tr><tr><td align="center" valign="middle" >Efficiency</td><td align="center" valign="middle" >0.6</td></tr><tr><td align="center" valign="middle" >Blockage coefficient</td><td align="center" valign="middle" >0.96</td></tr></tbody></table></table-wrap><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x19.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x20.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x21.png" xlink:type="simple"/></inline-formula> (7)</p><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x22.png" xlink:type="simple"/></inline-formula> is speed on the rotor tip. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x23.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x24.png" xlink:type="simple"/></inline-formula> are the diameter of the rotor tip and hub respectively. Q is the volume flow rate, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x25.png" xlink:type="simple"/></inline-formula>is the static pressure rise, ρ is atmospheric density, T is the rotor torque and ω is the rotating speed of the rotor.</p><p>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. Three- dimensional 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 [<xref ref-type="bibr" rid="scirp.73490-ref8">8</xref>] .</p></sec><sec id="s3"><title>3. Design Result</title><p>The streamlines on the meridional plane of half-ducted fan were shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. 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 [<xref ref-type="bibr" rid="scirp.73490-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.73490-ref5">5</xref>] 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 <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p></sec><sec id="s4"><title>4. Numerical Method and Result</title><p>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 <xref ref-type="table" rid="table2">Table 2</xref> and <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p>The computation domain was extended 300 mm upstream from the leading edge and 500 mm downstream from the trailing edge, respectively. <xref ref-type="fig" rid="fig3">Figure 3</xref> shows</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Design result of meridional streamline</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x26.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Design result of NACA65 blade shape</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x27.png"/></fig><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Simulation conditions</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Software: Ansys CFX Release 16.0</th></tr></thead><tr><td align="center" valign="middle" >RANS (Reynolds Averaged Navier-Stokes Simulation)</td></tr><tr><td align="center" valign="middle" >Air Temperature: 25˚C (Density: 1.185 kg/m<sup>3</sup>)</td></tr><tr><td align="center" valign="middle" >Turbulence model: SST k-ω model</td></tr><tr><td align="center" valign="middle" >Relative Coordinate System (Shaft Speed: 3000 min<sup>−1</sup>)</td></tr><tr><td align="center" valign="middle" >Mass Flow Inlet (0.1971 kg/s, Uniform)</td></tr><tr><td align="center" valign="middle" >Pressure Static Outlet (1 atm, Uniform)</td></tr><tr><td align="center" valign="middle" >Wall Boundary: Dirichlet boundary, No Slip (Counter Rotating Wall in Shroud)</td></tr></tbody></table></table-wrap><p>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.</p><p>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.</p><p><xref ref-type="fig" rid="fig4">Figure 4</xref> 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</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Simulation layout</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x28.png"/></fig><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Designed and simulated meridional streamline.</title></caption><fig id ="fig4_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x30.png"/></fig><fig id ="fig4_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x29.png"/></fig></fig-group><fig-group id="fig5"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Meridional velocity.</title></caption><fig id ="fig5_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x32.png"/></fig><fig id ="fig5_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x31.png"/></fig></fig-group><fig-group id="fig6"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Tangential velocity.</title></caption><fig id ="fig6_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x34.png"/></fig><fig id ="fig6_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x33.png"/></fig></fig-group><fig-group id="fig7"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Radial velocity.</title></caption><fig id ="fig7_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x36.png"/></fig><fig id ="fig7_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x35.png"/></fig></fig-group><p>tangential velocities are normalized by the rotational velocity at the blade tip<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x37.png" xlink:type="simple"/></inline-formula>. In addition, the radius r is normalized by the values of the hub <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x38.png" xlink:type="simple"/></inline-formula> and casing<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x39.png" xlink:type="simple"/></inline-formula>. The positive direction of the meridional, tangential and radial velocity indicate the main flow direction, rotational direction and radial direction, respectively. The negative directions of them mean the reverse directions. In this method, the radial inclination of meridional streamline was taken into account for the half-ducted fan design. Especially, at rotor inlet the negative radial velocity at tip side is also important value as well as the positive radial velocity at hub in <xref ref-type="fig" rid="fig7">Figure 7</xref>. For the case at rotor inlet, at the flow rate coefficient of f = 0.264, the CFD values of the meridional, tangential and radial velocity are near the design value. The meridional, tangential and radial velocity come close so well between design and CFD. The vortex design of the constant tangential velocity 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 f = 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 <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x40.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x41.png" xlink:type="simple"/></inline-formula> 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.</p></sec><sec id="s5"><title>5. Experiment</title><p>The designed and calculated half-ducted rotor with four NACA65 blades was fabricated by a three-dimensional printer shown in <xref ref-type="fig" rid="fig8">Figure 8</xref> 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 <xref ref-type="fig" rid="fig9">Figure 9</xref> with an 880 mm &#215; 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<sup>−1</sup> 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</p><fig-group id="fig8"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Fabricated rotor.</title></caption><fig id ="fig8_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x43.png"/></fig><fig id ="fig8_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x42.png"/></fig></fig-group><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Experimental apparatus</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x44.png"/></fig><p>rotating speed of the booster fan. Axial torque coefficient and efficiency are defined as:</p><disp-formula id="scirp.73490-formula56"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x45.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.73490-formula57"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2320364x46.png"  xlink:type="simple"/></disp-formula><p>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 <xref ref-type="fig" rid="fig1">Figure 1</xref>0. 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<sup>−1</sup> in this paper. Calculated and measured pressure rise, torque and efficiency is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>0. For the comparison between the design values and the experimental values at the design flow rate coefficient of f = 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</p><fig-group id="fig10"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> Stage performance.</title></caption><fig id ="fig10_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x48.png"/></fig><fig id ="fig10_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x47.png"/></fig><fig id ="fig10_3"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2320364x49.png"/></fig></fig-group><p>flow rate coefficient of f = 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.</p></sec><sec id="s6"><title>6. Conclusions</title><p>The half-ducted axial flow fan was designed by a diagonal flow fan design method.</p><p>Half-ducted fans have been designed, numerically analyzed, manufactured and tested to confirm the effect of this design method for improving design method on the overall performance and the three-dimensional flow field in design. The conclusions are summarized as follows.</p><p>1) For the comparison between the design values and the experimental values at the design flow rate coefficient of f = 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.</p><p>2) For the case at rotor outlet, at the flow rate coefficient of f = 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 <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x50.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2320364x51.png" xlink:type="simple"/></inline-formula> at rotor outlet 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.</p><p>3) For the comparison between the experimental values and the CFD values at the design flow rate coefficient of f = 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.</p></sec><sec id="s7"><title>Acknowledgements</title><p>The authors would like to acknowledge the financial support of Saga University for this research project.</p></sec><sec id="s8"><title>Cite this paper</title><p>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. http://dx.doi.org/10.4236/ojfd.2017.71002</p></sec></body><back><ref-list><title>References</title><ref id="scirp.73490-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Zangeneh, M. (1996) Inverse Design of Centrifugal Compressor Vaned Diffusers in Inlet Shear Flows. ASME Transactions on Journal of Turbomachinery, 118, 385-393. https://doi.org/10.1115/1.2836653</mixed-citation></ref><ref id="scirp.73490-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Zangeneh, M., Goto, A. and Harada, H. (1998) On the Design Criteria for Suppression of Secondary Flows in Centrifugal and Mixed Flow Impellers. ASME Transactions on Journal of Turbomachinery, 120, 723-735.  
https://doi.org/10.1115/1.2841783</mixed-citation></ref><ref id="scirp.73490-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Goto, A. (2016) Historical Perspective on Fluid Machinery Flow Optimization in an Industry. International Journal of Fluid Machinery and Systems, 9, 75-84.  
 https://doi.org/10.5293/IJFMS.2016.9.1.075</mixed-citation></ref><ref id="scirp.73490-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Shiomi, N., Kinoue, Y. and Setoguchi, T. (2012) Three Dimensional Velocity Fields at Rotor Outlet of a Semi-Opened Propeller Fan. Turbomachinery, 40, 218-225. (In Japanese)</mixed-citation></ref><ref id="scirp.73490-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Shiomi, N., Kinoue, Y. and Setoguchi, T. (2012) Experimental Study on Flow Fields with Vortex in a Semi-Opened Propeller Fan. Turbomachinery, 40, 688-696. (In Japanese)</mixed-citation></ref><ref id="scirp.73490-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Liu, P., Oka, Y., Kinoue, Y., Shiomi, N., Setoguchi, T. and Jin, Y. (2013) Design of Half-Ducted Axial Flow Fan Considering Radial Inflow and Outflow. Open Journal Fluid Dynamics, 3, 1-8.  https://doi.org/10.4236/ojfd.2013.32A001</mixed-citation></ref><ref id="scirp.73490-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Novak, R.A. (1967) Streamline Curvature Computing Procedures for Fluid-Flow Problems. ASME Transactions on Journal of Engineering for Power, 89, 478-490. 
 https://doi.org/10.1115/1.3616716</mixed-citation></ref><ref id="scirp.73490-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Inoue, M., Ikui, T., Kamada, Y. and Tashiro, M. (1980) A Quasi Three-Dimensional Design of Diagonal Flow Impellers by Use of Cascade Data. Proceeding of 10th Symposium of IAHR, Tokyo, 403-414.</mixed-citation></ref></ref-list></back></article>