Numerical Prediction of Hydromechanical Behaviour of Controllable Pitch Propeller

The research described in this paper was carried out to predict hydrodynamic and frictional forces of controllable pitch propeller (CPP) that bring about fretting problems in a blade bearing.The governing equations are Reynolds-averagedNavier-Stokes (RANS) and are solved by OpenFOAM solver for hydrodynamic forces behind the ship’s wake. Frictional forces are calculated by practical mechanical formulae. Different advance velocities with constant rotational speed for blades are used to achieve hydrodynamic coefficients in open water and the wake behind the propeller. Results are compared at four different pitches. Detailed numerical results of 3D modelling of the propeller, hydrodynamic characteristics, and probability of the fretting motion in the propeller are presented. Results show that the probability of the fretting movement is related to the pitch.


Introduction
The main goal of any propulsion system is to obtain the required speed of a vessel.One way to increase the control, manoeuvrability, and efficiency is to utilize controllable pitch propellers instead of constant pitch ones.Specific characteristics of controllable propellers permit the blades to rotate around the propeller axle.Blades rotate with complex mechanism which acts mechanically or hydraulically.The ability of the CPP is to adjust the pitch in order to avoid the motor to be overloaded and lead to a decrease in shaft speed as well as to increase the motor duration.During recent years numerous investigations have been carried out on propeller designs.Consequently different methods have been presented.
Because of different application of propellers numerous layouts of propeller have emerged.Today CPP is used especially on vessels which operate at various operation conditions such as containers, fishing vessels, and towing vessels.For the first time in 1840 CPP was introduced.At that time practical usage of CPP was in gas turbines.In twentieth century when diesel was introduced, usage of CPP increased suddenly.At that time pitch control was simple and mechanical.Dutch Schooner was the first one who applied CPP for diesels.In 1934, Escher Wyss used hydraulic mechanism in CPP which can be used for larger vessels.In 1937, Kamewa propellers were introduced.The hub mechanism of the mentioned propellers was based on Kaplan turbine design [1,2].
In recent years numerous researches have been conducted on mechanism and the controlling system of CPP and on the ways to decrease the friction resistance and vibrations of the blades.Tendency to increase both hydrodynamic forces in CPP and usage of hydraulic systems for pitch adjustment leads to more investigations on CPP.For example, in 2005, Bakker studied different control strategies.He introduced mechanical systems of CPP propellers in detail [3].In 2006, Dallinga calculated the forces acting on the CPP in different sea conditions [4].In 2009, Stuart et al. calculated the forces acting on blades under cavitation condition [5].A wake model for the prediction of propeller performance at low advance ratios was investigated by Tian and Kinnas [6].In 2009, Godjevac et al. [7] focused on the forces acting in a controllable pitch propeller, mechanism, and the manifestation of vibrating motion (fretting) in a blade bearing caused by seaway.
Fretting is a special wear process which happens in contact surfaces which have a relative oscillatory motion of small amplitude, typically smaller than 1 mm.Usually, the motion is caused by vibrations of the machinery.The fretting wear can occur at the CPP of the contact between blades and hub due  to variable hydrodynamic loading acting on the blades.The critical moment for the commencement of fretting motion is when the spindle torque conquers the friction torque.As the thrust and torque of a blade change during one revolution, the Coulomb friction in the blade bearing alters as well during one revolution.Looking at the Coulomb friction in the blade bearing and spindle torque of the corresponding blade, it is possible to illustrate the occurrence of fretting motion.Figure 1 shows the CPP configurations and its coordinate system.Blade bearing is the bearing that connects the blade foot, blade carrier, and hub.Blade bearing has the radial and axial parts.Radial part of blade bearing has the inner side (closer to x-axis) and outer side [8].
The -, -, and -axes are used to explain the forces and moments acting in a blade bearing.It also shows the radial and axial part of the blade bearing.The blade foot and the blade carrier are fastened with bolts.The critical moment for the beginning of fretting motion is when the spindle torque overcomes the friction torque.As the thrust and torque of a blade change during one revolution, the friction force in the blade bearing changes as well during one revolution.Looking at the friction force in the blade bearing and spindle torque of the corresponding blade, it is possible to illustrate the happening of fretting motion.
The efficiency of the propulsion system is strongly dependent on propeller performance, thrust force, torque of propeller, and its efficiency.Therefore this study focuses on the hydrodynamic performance and prediction of fretting motion in special CPP.The purpose of this study is to identify the effect of changes in pitch on propeller fretting motion.
Therefore, a simple method has been investigated to graph the propeller hydrodynamic coefficients with respect to advance coefficient and fretting occurring during one revolution of the blades.This investigation considers the situation that leads to fretting motion.When the hydrodynamic forces change during rotation of the blades, the friction in the blade bearings will change accordingly.By calculation of the friction in the blade bearing and spindle torque of the corresponding blade, it is possible to describe the occurrence of fretting motion.

Three-Dimension Modelling of the Investigated Propeller
B-series propeller is a famous series which is used in almost all practical usages [1].Pitch-diameter ratio in these series is Stationary Zone.This zone is cylindrical and includes the boss, propeller, and moving zone.We considered upstream length as 2, downstream length is 10, and the zone diameter is 3, as shown in Figure 2 where  is the diameter of the propeller.

Numerical Procedure
To investigate the hydrodynamic performance and flow fields around the propeller, Reynolds Navier-Stokes equations (RANS) for incompressible flow are applied.The governing equations are as follows.Conservation of mass (quantity equation): conservation of momentum: where   = body force.
Unstructured triangle cells are used to describe the control volume.By using the momentum, conservation of mass and energy equation around the propeller and the mathematical condition imposed at solid and fluid boundaries.Governing equations are solved.We should consider the time and cost in calculation.We used different number of meshes to reach its independency for the final results of propeller characteristics (  ,   , and   ).Size of meshes near the hub, blade edges, and tip of the blades is smaller than other parts.The domain was meshed by tetrahedral cells as shown in Figure 3. Instead of rotating the propeller, fluid region around the propeller is rotated.
To ensure the grid independency, different grids numbers have been chosen [10].It is found that with the number 880000 grids all parameters such as thrust coefficient are converged.The important values  + are less than 50 for all advance velocity ratios () when   is converged.Gird number, thrust coefficient, and its  + values are shown in Table 1.

Hydrodynamic Performance of the Investigated Propeller
The most important characteristics of the marine propellers are the hydrodynamic performances that are operated in where  and  are the thrust and torque of the propeller, respectively. is the advance velocity ratio and, defined by  =   /, means advance velocity (  ) divided by rotational speed () and diameter ().
Here, we selected that the rotational speed of the propeller is constant and equal to 120 rpm and diameter is 0.508 [m].The advance velocity is changed in order to calculate the propeller performance.The hydrodynamics characteristics of the propeller are performed and compared at four pitchdiameter ratios, that is, / = 0.5, 0.8, 1.1, and 1.4.Figures 4, 5, 6, and 7 show the comparison of the hydrodynamic coefficients.General trends of the results are good but there are some discrepancies at off-design conditions.The velocity and pressure contours around the propeller in  = 1 are shown in Figures 8 and 9, respectively.As observed in Figure 8, the pressure shows at face side of the propeller.The pressure is different on each blade due to nonuniform wake flow.

Friction Forces
CPP mechanism is much complicated because the propeller should utilize the engine power to obtain the optimum efficiency at all conditions.Propeller operates behind the ship hull where the flow is nonuniform and oscillating flow.Thus, the resultant velocity will change during one cycle and as a result the thrust will be oscillated.For this reason the propeller pitch should be adjusted based on the required thrust loading.Therefore, the blade pitch changes at any time during one revolution.The blade seat on the hub is made by a bearing that supports forces in both axial and radial directions, as shown in Figure 10.
The friction forces due to controlling the blade on the hub bearing are divided into two axial and radial directions.The forces and moments have been evaluated by the procedure proposed [7].The resulting forces in the axial part of the blade are defined through   where  hd, and  hd, are the components of the hydrodynamic forces, means  hd, = (Thrust),  hd, = (Torque)/  .  means the radius of centroid of torque. ce, is the component of centrifugal forces.The axial torque due to the axial force is given by where   is the distance between the point where the axial force is supposed to be applied and the hub friction seat of the propeller blade root.  is the nondimensional frictional coefficient of the blade and hub materials.These parameters are given by the following formulae [2]: ,   =  (1 + 0.0477 sin 2  + 0.544sin 4  −1.051sin 6  + 0.6982sin 8 ) .Advance coefficient (J) where  hd, and  hd, are the components of the hydrodynamic moments and  ce, is the component of centrifugal moments.The hydrodynamic moments acting on the centre of hydrodynamic (CH) of a CPP blade are calculated by where ( CH ,  CH ,  CH ) is the hydrodynamic force centre.Finally, the friction torque in the radial part of the blade bearing is expressed as [11]  Fr,rad = 4  ⋅  ⋅   (cos Φ 0 +  ⋅ arcsin ()) , where  is a coefficient that describes the ratio of the radial force and the maximum force resulting from the bending moment and Φ 0 is the angle of circumference around blade bearing.Total friction torque of the blade bearing is obtained by  Fr = √ 2 Fr, +  2 Fr,rad .The fretting motion will happen when the blade spindle torque becomes higher than the friction torque, which means  hd, >  Fr .

Results and Discussion
Fretting motion starts when the sum of the forces acting on the body overcomes the friction force.As mentioned, it may occur when the blade spindle torque becomes bigger than the friction torque.In order to investigate the start of the fretting motion in a CPP, the relation between friction torque and the spindle torque during one revolution was accomplished.
In this study, we focused on the effect of different pitches on fretting behaviour of the CPP.Propeller is working in behind the ship hull where the wake flow which is into the propeller is nonuniform and oscillating flow.Wake factor (1 − ) is defined by   /  , where   is advance velocity behind the ship and   is the speed of the ship.Due to oscillating wake filed, the hydrodynamic forces will change in the wake field behind the propeller.The investigated wake field (1 − ) during one revolution of the blades shows different radius of the blades in Figure 11.The results of the spindle torque coefficient for different pitches for one revolution are presented in Figure 12.The mean values of the hydrodynamic Forces are presented in nondimensional values in Table 2.
behaviour of the propeller at different pitches for the investigated propeller are presented.It shows that for the current propeller, when pitch ratio increases, the possibility of the fretting motion increases.
In the near future, an investigation will be carried out if the geometry of this or other propellers is made available.Also, effect of the ocean wave may be studied on the fretting problem.More accurate formulae can be achieved to determine the frictional forces on the blade bearing.

Figure 3 :
Figure 3: Final mesh for the calculation domain.

( 6 )
is the radius of blade bearing.The angle  depends on the axial force, friction coefficient, radius of blade bearing and
open-water flow.The hydrodynamic open-water characteristics are expressed as follows: