Investigation of the effects of a fan-shaped Mewis duct before a propeller on propulsion performance

To define the guidelines for the design of a Mewis duct for a small bulk cargo ship, a numerical study was carried out to investigate the effects of a fan-shaped Mewis duct on propeller performance of a 38,000 t bulk cargo ship with a three-blade propeller. Calculations were performed using the STAR-CCM+ software based on the solution of the Reynolds averaged Navier–Stokes (RANS) equation. Computations were carried out for a wide range of locations from 0.11D to 0.21D upstream from the propeller disc, a radius of trailing edge from 0.65R to 1.05R, and angle of attack of the inner fins, where from 0° to 10° of the fan-shaped Mewis duct, where D is the diameter of the propeller and R is the radius of the propeller. With the installation of the duct, the axial mean wake fraction on the propeller disc decreased, while the tangential mean wake fraction increased significantly. As the installation position was far away from the propeller, the axial mean wake fraction varied in a wave pattern, while the tangential mean wake fraction increased. With the increase of the radius and the angle of attack of the inner fins of the duct, the axial mean wake fraction first decreased and then increased, while the tangential mean wake fraction increased gradually. Correspondingly, the propeller force coefficients and efficiencies with the duct were all more than those without the duct and changed with the variations of the duct parameters. When the propeller and duct were taken as a propulsion system, the efficiency of the system was less than that of only the propeller. At the same time, with variations of the location, radius, and angle of attack of the inner fins of the duct, the efficiency of the system first increased and then decreased gradually, which indicates that there was an optimum design value for all of the duct parameters. Moreover, only when the design of these duct parameters are reasonable can energy savings be achieved; otherwise, the opposite effect is produced. Eventually, with the installation of the fan-shaped Mewis duct, the wake harmonic degree on the propeller disc can be significantly improved.


Introduction
At present, marine propellers are still the most commonly used propulsion device. Energy savings and low-carbonemission output are currently the two inevitable trends in the shipbuilding industry. Meanwhile, the aims of ship energy savings and emission reductions have promoted the rapid development of the research in energy-saving devices (ESDs), which can be installed before the propeller, at the propeller station, and after the propeller [1]. ESDs can reduce energy loss by producing additional thrust, improving the inflow on the propeller plane, or recovering the rotation energy of the tail of the propeller.
The propeller works in the wake field at the stern of the ship. Due to the influence of the wake, the actual efficiency of the propeller is different from that of being in open water. A high wake zone is usually found in the upper half of the propeller disc which results in the decrease of propulsion efficiency, the deterioration of cavitation performance, and the increase of exciting force levels [2]. Many additional devices have been designed to improve the hydrodynamic performance and to reduce cavitation and noise. A lot of research has been done on the working mechanisms and effectiveness of these devices [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18].
For the equalization of the propeller inflow and improvement of the propulsion efficiency, one of the energy-saving devices used widely in ships is the pre-duct such as the wake equalizing duct (WED) and the Mewis duct, while Schneekluth [19] introduced the WED first. In fact, the WED was first fitted onto a large ship in March 1984, while model tests were carried out beforehand. WEDs aim to improve the overall propulsion efficiency of the ship by improving the inflow of the propellers.
To explore the special features of rudder balls and WEDs, Dayuan [20] carried out model tank tests for the 65,655 DWT bulk carriers and satisfactory results were obtained. He believed that the simultaneous use of rudder balls and WEDs was a practical solution for energy saving in fully laden vessels.
To investigate the scale effects of the ducts, Friesch and Johannsen [21] performed resistance and self-propulsion tests on a large single screw tanker ship with and without a WED in the large hydrodynamics and cavitation tunnel of the Hamburg Ship Model Basin (HYKAT). They reported that the effectiveness of the energy savings of a WED would decrease when increasing the Reynolds number further and disappear totally at full scale. They thought that the measurements should be understood as a possible explanation for the discrepancy between predicted gain and reality.
To study the effect of a partial WED on the powering characteristics of the subject vessel, Korkut [22] carried out resistance, self-propulsion, and flow visualization measurements for a river going general cargo ship. To improve the stern flow characteristics of the original hull form, the stern form of the vessel was optimized and two different hull forms were generated. He found that a partial WED with optimized stern forms would result in less resistance and effective power requirement than those of the original hull form, which indicates that the duct can improve the powering characteristics.
Çelik [23] numerically investigated the effect of a WED on the propulsion performance of a chemical tanker. The propulsion characteristics, path lines over the after body, and viscous resistance forces were analyzed for six WED cases combined with the ship and propeller at various ship speeds. He thought that the ship's viscous resistance forces and propeller propulsive characteristics were not changed with ducts, but the WEDs generated a positive thrust, which resulted in a maximum increase of the propulsive efficiency by 9.7%. He concluded that a well-designed WED can improve the propulsion characteristics of a ship considerably, while a poorly designed WED can cause negative effects.
Mewis [24][25][26] developed the Mewis duct for smaller container ships and bulk carriers by combining a WED with pre-swirl fins. Studies show that a Mewis duct can lead to an energy savings between 7 and 9%.
Hans-Jürgen Heinke [27] carried out calculations by computational fluid dynamic (CFD) in model and full scale for a typical container ship with a WED to investigate the scale effect on the flow around the WED and on the inflow to the propeller. The calculation results showed that the change of the propulsion coefficients thrust deduction fraction, wake fraction, and hull efficiency of the ships with WEDs can be predicted with reasonable accuracy using the ITTC 1978 propulsion prognosis method. The use of a WED leads to a reduction of the propeller induced pressure fluctuations in the first and second harmonic orders.
Carlton [1] pointed out that a WED can help to establish a more uniform inflow into the propeller to reduce the amount of separation over the after body of the vessel by accelerating the flow in the upper region of the propeller disc, and by attempting to minimize the tangential velocity components in the wake field.
Feizi Chekab and Ghadimi [28] investigated noise reduction methods for submerged marine propellers by modifying the wake inflow with a WED and a Mewis duct. He reported that the efficiency of the propulsion system can be improved, and the propeller's cavitation and vibrations could be reduced with the modification of wake inflow.
Go and Yoon [29] aimed at providing the general information for the duct design regardless of the ship type, and considered a wide range of duct diameters and angles of the attack. They numerically investigated the effects of duct diameters and angle of attack on the hydrodynamic forces and efficiency of a marine propeller by carefully analyzing the related flow fields. He imposed a uniform inflow condition without considering the type of ship to investigate the net effect of a duct on propeller performance. He found that when the drag of the duct excluded the thrust coefficient, torque coefficient, and propeller efficiency of the propeller all augments with the increasing angle of attack of the duct and all coefficients of the propeller with the duct were larger than those of the propeller without the duct. In addition, it was found that a duct with a high angle of attack provides faster inflow into the propeller than in the case without a duct.
According to our review of the references, most researchers that have studied the working mechanisms or the effects of a pre-duct believe that a pre-duct can help improve the propulsion performance of a vessel. However, the effect of energy savings of a pre-duct usually depends on the hull form of the stern. For a Mewis duct, in addition to the hull form, its several characteristic parameters have great influence on its performance including the installation position, radius, angle of attack of the inner fins, profile of the inner fins and outer duct, and so on. To achieve better energy-saving effects, these parameters should be optimized. However, in the past, the design of a pre-duct was largely a matter of experience based on model tests, because it is difficult to find guidelines for the design of these parameters of the duct independent of hull form. Thus, the designers had to test and verify their designs in actual use, which often makes the applications of pre-ducts questionable.
In the present study, for a small bulk cargo ship, the effects of a fan-shaped Mewis duct with different characteristic parameters were investigated to find the effect it had on the propulsion performance of the propeller. The fan-shaped Mewis duct was developed by the Shanghai Merchant Ship Design and Research Institute. The choice of characteristic parameters of a Mewis duct depends on the ship type, and it is futile to find guidelines applicable to all ship types. Therefore, the conclusions from the present study are only applicable to small bulk cargo ships. Our research mainly focused on the relationship between propulsion performance of propeller and characteristic parameters of the Mewis duct. If this relationship is known, the optimized design of the Mewis duct can be created more efficiently for a small bulk cargo ship.
For a 38,000-t bulk cargo ship, a wide range of duct locations, duct radius, and angle of attack of the inner fins were considered and the effects of these parameters on the hydrodynamic forces, efficiency, and wake harmonic degree were numerically investigated by carefully analyzing the wake fractions of the related flow fields. The relationship between the propulsion performance of propeller and the characteristic parameters of the Mewis duct was needed, and it is believed that the results will be helpful to the design of the Mewis duct for small bulk cargo ships.

Governing equations
The three-dimensional incompressible steady-state flow problem is governed by continuity and the Navier-Stokes equations [29], which can be described as where x i are Cartesian coordinates,u i are the velocity components, ρ is the water density, p is the pressure, µ is the dynamic viscosity, and − u � i u � j is the Reynolds stress term closed by the SST k-ω turbulence model with the standard wall function. This turbulence model was recommended by the Specialist Committee on CFD in Marine Hydrodynamics . These equations are discretized into algebraic equations based on the finite-volume method. The commercial code STAR-CCM+ was used for solving the algebraic equations. The convection term of the momentum equation was discretized by QUICK, while the diffusion term was discretized by the second-order central difference scheme. The SIMPLEC algorithm was used for the velocity-pressure coupling and the moving reference frame (MRF) scheme was adopted for the rotation of the propeller. A convergence criterion of 10 −6 was used, and further details can be found in the STAR CCM+ manuals [37].

Boundary conditions
A fixed Cartesian coordinate system with respect to the ship was used. As shown in Fig. 1, the computation domain was defined by the boundaries as follows. The inflow boundary, with a 1.5 L WL upstream distance from the ship's midsection, was set as the velocity-inlet boundary condition, where L WL is the length of water line of the ship. The outflow boundary, with a 3.5 L WL downstream distance from the ship's midsection, was set as the undisturbed pressure-outlet boundary condition. The outer boundary, with a 1.5 L WL distance from the central axis of the model, was set as the undisturbed velocity boundary condition. The hull surface and duct were set as no-slip boundary conditions. The propeller blades were set as rotating non-slip surfaces. in Table 1, while the hull's outline is shown in Fig. 2.

Geometry
The principal parameters of the propeller are shown in Table 2. The propeller with duct is shown in Fig. 3 and the back view of the fan-shaped Mewis duct is shown in Fig. 4.
The fan-shaped Mewis duct consists of an outer duct and four inner guide fins, which was designed using NACA0015. The angle of attack of the outer duct was 15° and a chord length of 0.28D was used, where D is the diameter of the propeller. The radius of the leading edge was 1.0R, where R is the radius of the propeller. The chord length of the inner fins was 0.18D. The inner fins were installed on the shaft of the propeller and connected with the outer duct. As shown in Fig. 4, β 1 , β 2 , β 3 , and β 4 were the intersection angles between the guide fins and the longitudinal centerline of the hub, where β 1 and β 4 were 72°, β 2 was 28°, and β 3 was 42°. α 1 , α 2 , α 3, and α 4 were the angle of attack of the inner fins of the duct.

Grid system and validation
The Cartesian cut cell method was used to generate a computational mesh to cover the global fixed and moving block. The generation of the present grid system follows the same procedure as in Ref. [2].  Locally refined and hybrid grids were deployed in the computational domain grid division to guarantee the computational precision and operation speed requirements. In the stern area, including the propellers and duct, unstructured grids were adopted, while the rest used structured grids. The grids are refined to ensure that the aspect ratios of the numerical cells were less than 300 at the middle of the hull and 100 at the bow. The computational grid is shown in Fig. 5.
The grid verification study was performed with three different grids (coarse, medium, and fine) to determine the effect of grid systems on the numerical solution, as shown in Table 3. The results of verification and validation are listed in Tables 4, 5 and 6, using total resistance coefficient C T , thrust coefficient, K T , torque coefficient, K Q , and propeller efficiency, η 0 [38].
As shown in Table 5, the grid converges ratios R G of the four parameters is all less than 1, which suggests that the computational grids exhibit monotonic convergence. A comparison between the calculated results according to different grid systems with the experimental data is listed in Table 6. It can be seen that even if the number of the three sets of grids varies greatly, the errors of K T , K Q , and η 0 are always controlled within 5%, which shows that the solutions depend on the grid systems weakly. According to the error of efficiency, η 0 , we believe that the medium grid system is the most suitable one to guarantee the calculation accuracy and control computational time.

Results and discussion
Assuming that the thrust, torque, and rotational speed of the propeller is T 0 , Q 0, and n, respectively, and the ship's total resistance is R T0 in the cases without the fan-shaped Mewis duct, the efficiency of the propeller without a duct is When moving with constant forward speed V m and rotational speeds n, assuming that propeller thrust, propeller torque, and the ship's total resistance are T, Q, and R T, respectively, the efficiency of the propeller with a duct is With regard to the propeller and the fan-shaped Mewis duct as a propulsion system, the thrust of the system T E is    In addition, the efficiency η E can be expressed as To predict speed performance, it is necessary to obtain a self-propulsion point of the ship for the given V m . The self-propulsion factors are obtained using the revised ITTC'78 method [39]. The towing force, F D , that is difference in force between the viscous ship resistance and the propeller thrust [39,40].
The total resistance and propeller hydrodynamic force of the ship model with the design speed of 1.2731 m/s were calculated at three different propeller revolution speeds near the self-propulsion point. The results are shown in Table 7, according to which it can be obtained that the propeller revolution speed at the self-propulsion point is 9.11 rps. The total resistance, the propeller hydrodynamic force, and the mean wake fraction on the propeller disc are recalculated at the self-propulsion point, as shown in Table 8.
Under the condition that the design speed was 1.2731 m/s and the propeller speed was 9.11 rps, the effects of duct location, duct radius of the trailing edge, and the angle of attack of the inner fins on the propulsion performance were analyzed.

Duct location
As shown in Fig. 6, the duct location refers to the longitudinal upstream distance between the trailing edge of the duct and the center of the propeller. The angle of attack of the inner fins was taken as 10°, while the radius was 0.85R. The duct locations used were: 0.11D, 0.135D, 0.16D, 0.185D, and 0.21D.
The axial mean nominal wake fraction and tangential wake fraction on the propeller disc according to the duct location are shown in Fig. 7. It can be seen from Fig. 7a that the axial mean wake fraction decreases when the fan-shaped Mewis duct was installed, which shows that the presence of the duct accelerates the wake flow of the ship. As can be seen from Fig. 7b, when the duct was installed the mean tangential wake fraction increased significantly, and when the duct location increased the mean tangential average wake fraction was augmented. The negative value represents that the flow and the propeller rotation are in opposite directions. The effect of the duct on the propeller performance was caused by the effects of the axial wake and tangential wake simultaneously. The decrease of the axial flow increased the inflow speed of the propeller, while the decrease of the tangential flow increased the relative revolution speed of the propeller. The propeller thrust increased, while the real  propeller speed did not change, thereby increasing the efficiency of the propeller. The hydrodynamic performance of the single propeller and the propulsion system consisting of the propeller and duct according to the duct location are plotted in Fig. 8. From Fig. 8, the following conclusions were drawn. The thrust, torque, and efficiency of the propeller with the duct were significantly larger than that without the duct, which means that the duct was favorable to improve the propeller performance. With the installation position of the duct further away from the propeller disc, the thrust, torque, and efficiency of the propeller fluctuated. The thrust and torque had the same changing trend, contrary to the efficiency. This shows that the different installation positions had an impact on the performance of the propeller.
As can be seen from Fig. 8a, when considering the propeller and duct as a propulsion system, the thrust of the system was smaller than that of the single propeller itself. This shows that the propeller efficiency had been improved with the installation of the duct, but additional resistance was also generated. Therefore, the energy-saving effect depends on whether the increase of propeller efficiency can outweigh the adverse effects caused by the additional resistance. Figure 8c shows that the efficiency of the propulsion system was significantly lower than that of the single propeller. With the change of the duct location, the efficiency of the propulsion system first increased and then decreased, reaching the maximum at the duct location of 0.16D, and only when the duct location was approximately between 0.13D and 0.18D was the energy savings achieved.
As discussed above, the wake flow field varied along with the parameters of the duct. Therefore, a poor design and installation of the duct will cause turbulent changes to the wake flow field, which leads to the increase of noise and exciting forces of the propeller. This is the main reason leading to the vibration of the ship's stern. Therefore, it was necessary to analyze the effect of the duct on the wake harmonic degree on the propeller disc. In this paper, the distribution of the axial velocity of the wake field was analyzed by Fourier transformations to analyze the wake harmonic degree at the radius 0.75R of the propeller disc when the duct parameters changed. Figure 9 shows the numerical results of the axial wake harmonic degree at the radius 0.75R of the propeller disc according to duct location, where smaller amplitudes indicate a more uniform wake flow on the propeller disc. The first eight terms of the amplitude are listed. The main components of the axial wake amplitude were the first three terms. After installing the duct, the axial wake flow was more uniform than that without the duct, and the first two terms of its amplitudes decreased significantly. The uniformity of the axial wake flow on the propeller disc surface was  optimal when the duct location was 0.16D, which is very beneficial in the reduction of propeller noise and excited forces.

Duct radius of the trailing edge
As shown in Fig. 10, the duct radius in this paper refers to the radius of the trailing edge of the duct. The angle of attack of the inner fins was taken as 10° and the duct location of 0.16D. The duct locations were 0.65R, 0.75R, 0.85R, 0.95R, and 1.05R, respectively. The axial mean nominal wake fraction and tangential wake fraction on the propeller disc with respect to the duct radius are plotted in Fig. 11. Figure 12 shows the hydrodynamic performance of the single propeller and the propulsion system consisting of the propeller and duct according to the duct radius. Figure 11a shows that the axial mean wake fraction on the propeller disc varied greatly when the duct radius was different. The axial mean wake fraction decreased rapidly when the duct radius was less than 0.85R while increasing rapidly when the duct radius was greater than 0.85R. When the duct radius was greater than the radius of the propeller, the wake fraction was greater than that without the duct. This shows that the outer duct plays a dominant role in the rectification and acceleration of the wake field of the ship when the duct radius is small. However, when the duct radius increases, the effect of the rectification and acceleration of the wake field of the outer duct decreases until it disappears. Meanwhile, the inner fins play a dominant role in blocking the the axial inflow of the propeller. As can be seen from Fig. 11b, the mean tangential wake fraction increased along with the increase of the duct radius, which is because the inner fins can affect larger regions of the inflow of the propeller as their length increases. Figure 12 shows the hydrodynamic performance of the single propeller and the propulsion system consisting of the propeller and duct according to the duct radius. As shown in Fig. 12, the thrust and torque of the propeller increased with the increasing duct radius, while the efficiency increased first and then decreased reaching the maximum at the duct radius of 0.85R. As shown in Fig. 12a, the total thrust of the propeller and the duct was greater than that of the single propeller when the duct radius was smaller than 0.75R and less than that of the single propeller when the duct radius was larger than 0.75R. This suggests that additional thrust was produced on the fan-shaped Mewis duct when the duct radius was small. In fact, additional thrust was produced on the outer duct, and additional resistance was caused by the inner fins. When the duct radius was small, the thrust was greater than the resistance; thus, additional thrust was generated on the fan-shaped Mewis duct. With the increase of the duct radius, the resistance increased more rapidly than the thrust because of the increase of the inner fins length, which results in the additional resistance on the fan-shaped Mewis duct. Figure 12c shows that the efficiency of the propeller and duct increased with the increasing duct radius, and then decreased rapidly when the duct radius was greater than 0.75R. When the duct radius was greater than 0.85R, the efficiency was lower than that of propeller without the duct, which indicates that the duct radius should be less than 0.85R. Figure 13 shows numerical results of the axial wake harmonic degree at the radius 0.75R of the propeller disc according to the duct radius. The first eight terms of the amplitude are listed. As can be seen from Fig. 13 Fig. 11 Wake for different duct radius: a axial wake fraction and b tangential wake fraction increase of the duct radius, the first term of the axial wake amplitude varied slightly, while the second and third terms of the axial wake amplitudes decreased gradually.

Angle of attack of the inner fins of the duct
The duct radius was taken as 0.75R and the duct location of 0.16D; the angle of attack of the inner fins was 0°, 5°, 10°,   15°, and 20°, respectively, as shown in Fig. 14. Figure 15a shows that the axial mean wake fraction on the propeller disc varied greatly when the angle of attack of inner fins was different. The axial mean wake fraction decreased rapidly when the angle of attack of the inner fins was less than 10° while increasing rapidly when the angle of attack of the inner fins was greater than 10°. This shows that the outer duct played a dominant role in the rectification and acceleration of the wake field of the ship when the angle of attack was small. However, when the angle of attack of the inner fins increased, the axial component of the wake flow decreased rapidly with the water flowing through the inner fins, thus reducing the axial inflow velocity of the propeller. It can be seen from Fig. 15b that the mean tangential wake fraction increased along with the increase of the angle of attack of the inner fins, which is because the greater the angle of attack of the inner fins is, the greater the tangential component of the wake flow will be.
The hydrodynamic performance of the single propeller and the propulsion system consisting of the propeller and duct according to the angle of attack of the inner fins is plotted in Fig. 16. From Fig. 16, the following conclusions are drawn. The angle of attack of the inner fins had great influence on the hydrodynamic performance of the propeller. The thrust, torque, and efficiency of the propeller increased with the increase of the angle of attack of the inner fins. It can be seen from Fig. 15 that when the angle of attack of the inner fins was less than 10°, the axial mean wake fraction decreased, while the tangential mean wake fraction increased; thus, the propeller efficiency increases gradually. When the angle of attack of the inner fins was greater than 10°, the mean axial wake fraction and the mean tangential wake fraction were all growing, and the effect of tangential flow apparently played a leading role. Thus, the propeller efficiency still increased with the increase of the angle of attack of the inner fins. As seen from Fig. 16a, the total thrust of the propeller and duct was less than that of the single propeller. The total thrust of the propeller and duct first increased and then decreased with the increase of the angle of attack of the inner fins, and reached the maximum when the angle of attack of the inner fins was 10°. This is because when the angle of attack of inner fins was greater than 10°, larger tangential wake flow not only made the relative rotation speed of the propeller higher, improving the propeller efficiency, but also resulted in greater additional resistance. This offset the gains of the improvement of the propeller efficiency. It can be seen from Fig. 16c that the efficiency of the propulsion system consisting of the propeller and the fan-shaped Mewis duct reached the highest at the 10° angle of attack for the inner fins. Moreover, only when the blade angle of attack was approximately between 15° and 12° was the energy-saving effect achieved. Figure 17 shows the numerical results of the axial wake harmonic degree at the radius 0.75R of the propeller disc according to the angle of attack of the inner fins. The first eight terms of the amplitude are listed. From Fig. 17, the following conclusions were drawn. When the angle of attack of the inner fins was small, the axial wake with the duct was more uniform than that without duct, and the first two terms of the axial wake amplitudes decreased remarkably. When the angle of attack of the inner fins was greater than 15°, the first term of the axial wake amplitude with the duct increased approaching that of without the duct, which shows that too large of an angle of attack of the inner fins was unfavorable to the uniformity of the axial wake. After installing the duct, the axial wake flow was more uniform than that without the duct, and the first two terms of its amplitudes decreased significantly. It can be understood that when the angle of attack of the inner fins was 10°, the axial wake flow at the propeller disc was most uniform.

Conclusions
The effect of a fan-shaped Mewis duct on propeller performance was numerically investigated by considering a wide range of locations from 0.11D to 0.21D, radius from 0.65R to 1.05R, and angle of attack of the inner fins from 0° to 10° of the fan-shaped Mewis duct for a small bulk cargo ship with a three-blade propeller. The profile of the outer duct and inner fins were all designed using NACA0015. The angle of attack of the outer duct was 15° with a chord length of 0.28D, while the chord length of the inner fins was 0.18D. The design speed of the ship model was 1.2731 m/s. The mean wake fractions, propulsion characteristics and wake harmonic degree were analyzed carefully. The main conclusions of the present study are summarized as follows.
The existence of the fan-shaped Mewis duct reduced the axial wake and increased the tangential wake on the propeller plane. With the decrease of the axial flow, the inflow speed of the propeller was increased, and with the increase   Fig. 17 Axial wake harmonic degree at r/R = 0.7 on propeller disc according to angle of attack of inner fins of the duct of the tangential wake, the relative rotation speed of the propeller was increased. This resulted in an increased propeller thrust. However, the real rotation speed of the propeller did not change; thus, it improved the efficiency of the propeller. As the installation position became further away from the propeller, the axial mean wake fraction varied in a wave pattern, while the tangential mean wake fraction decreased gradually. When not considering the additional resistance produced by the fan-shaped Mewis duct, with the increasing distance between the fan-shaped Mewis duct and the propeller, the thrust and torque of the propeller showed the same increasing patterns with that of the mean axial wake fraction. However, the efficiency varies in the opposite way. It was observed that the tangential flow was so small that the axial wake is a more important variable. In addition, the hydrodynamic forces and efficiency of the propeller itself were greatly improved because of the installation of the fan-shaped Mewis duct. All the coefficients of the propeller with the duct installed were larger than that of the propeller without the duct. However, if regarding the propeller and duct as a propulsion system, the thrust and efficiency of this system were significantly smaller than that of propeller itself. In addition, only when the fan-shaped Mewis duct was mounted approximately between 0.13D and 0.18D of the distance from the propeller disc, the efficiency of the propulsion system was higher than that in the case without installing the fan-shaped Mewis duct. The maximum efficiency was achieved near 0.16D.
With the increase of the radius of the fan-shaped Mewis duct, the axial mean wake fraction first decreased and then increased, while the tangential mean wake fraction increased quickly. Both thrust and torque of the propeller increased, while the radius of the fan-shaped Mewis duct increased. The efficiency of the propeller increased first and then decreased. The propeller performance was affected by the axial wake and the tangential wake simultaneously, but the tangential wake was a more important variable. The maximum efficiency of the propulsion system consisting of the propeller and the fan-shaped Mewis duct was achieved when the duct radius was equal to 0.75R. Moreover, only when the duct radius was less than 0.85R did the efficiency of the propulsion system reach higher values than that of the propeller without the fan-shaped Mewis duct.
With the increase of the angle of attack of the inner fins of the fan-shaped Mewis duct, the axial mean wake fraction first decreased and then increased, while the tangential mean wake fraction increased quickly. The thrust, torque, and efficiency of the propeller all increased, while the angle of attack of the inner fins of the fan-shaped Mewis duct increased. Thus, the tangential wake was more important than the axial wake. The efficiency of the propulsion system consisting of the propeller and the fan-shaped Mewis duct increased gradually along with the increase of the angle of attack of the inner fins, and reached the maximum near 10°. When the angle of attack of the inner fins was greater than 10°, the efficiency dropped rapidly. When the angle of attack of inner fins was approximately between 5° and 12°, the efficiency of the propulsion system was higher than that of the propeller without the fan-shaped Mewis duct. This is because the tangential wake increased along with the increase of the angle of attack of the inner fins, which was equivalent to increasing the relative revolution speed of the propeller thus increasing the propeller thrust. However, when the angle of attack of the inner fins was greater than 10°, the resistance of the fan-shaped Mewis duct increased faster than the propeller thrust.
The axial wake harmonic degree on the propeller disc had been significantly improved with the installation of the fan-shaped Mewis duct, which means that the propeller noise and excited forces were effectively controlled. There were also optimal parameters of the fan-shaped Mewis duct corresponding to the best axial wake harmonic degree, but the parameters were sometimes inconsistent with that corresponding to the best energy efficiency.
The research in this study showed that the thrust of the propeller was increased with the installation of the fanshaped Mewis duct, but the ship's total resistance was also increased. Therefore, to achieve energy savings, it is necessary to optimize the parameters of the fan-shaped Mewis duct to ensure that the gains from the increase of the propeller thrust is greater than the adverse effects caused by the increase in resistance.
Consequently, it is believed that the main findings given in this study are helpful in understanding the relationship between that of the propulsion performance of the propeller and the parameters of the Mewis duct, which is very beneficial to the design of these parameters. In the future, more parameters such as the angle of attack and the profile of the outer duct and inner fins should be considered.