3.1. Unsteady Exciting Force
In the present paper, three conditions of
J = 0.6, 0.8217 (design condition), 1.0 are considered and simulated.
Figure 9 shows the time histories of unsteady
KT and
KQ under
J = 0.8217.
It can be seen that similar to traditional propeller system, the unsteady thrust and torque also fluctuate close to a certain value slightly. However, the situation becomes more complicated because of the pod, which is downstream of the propeller and will interferes with the propeller.
Figure 10 shows the distributions of pressure coefficient on the pressure side under
J = 0.8217 with time step
t0,
t0 + 5Δ
t,
t0 + 10Δ
t, and
t0 + 15Δ
t.
The rotation of the propeller is right hand. It is evident that compared to the blades which are far away from the strut, the blades near the strut have a higher region of the pressure. With the rotation of the propeller, the high-pressure region on the blade that gradually away from the strut disappears slowly. However, a new high-pressure region slowly appears on the blade, which is gradually approaching the strut. Moreover, the largest region of high pressure is on the blade that is sweeping the strut. The wake flow field of the propeller is inhomogeneous and periodic. When the water flows impact the strut periodically, the strut will be subjected to periodic exciting force. Due to the interaction of forces, the propeller will also get exciting force. Consequently, the pressure side on the blade that is sweeping the strut will generate a high-pressure region. This phenomenon is called rotor–stator interaction. The smaller distance between the rotor and the stator, the stronger effect of rotor–stator interaction will be generated.
It can also be observed that when the blades sweep the fin in the bottom, there seems no special region of high pressure appears on the pressure side of the blade. The reason is that the distance between the propeller and fin in the bottom is larger than it between the propeller and the strut, therefore the effect of rotor–stator interaction has been decreased significantly. Owing to rotor–stator interaction, the podded propulsor will obtain higher thrust than traditional propeller system. Besides, compared to the normal propeller which often been mounted behind the ship directly, the installation position of the podded propulsor is far away from the wake flow field of the ship, which will lead to a more uniform inflow condition. Combined with the characteristic of high thrust, the efficiency of the podded propulsor will be higher than other traditional propellers.
Figure 11 shows the frequency domain result of the thrust coefficient under
J = 0.6, 0.8217, 1.0. It is obvious that the discrete frequency is
fBPF. According to the normal law of the open water performance of the propeller, it is definite that the thrust provided by the propeller will decrease with the increase of the advance coefficient. However, in the frequency domain, the energy contribution of the thrust is variant under different conditions, as is shown in
Figure 11. Under three conditions, although the main frequencies are all blade passage frequencies, the amplitudes are different; there is the maximum amplitude of the thrust in the blade passage frequency under
J = 1.0, followed by 0.8 and 0.6. Obviously, the energy is more centralized under a higher-advance coefficient, while some energy is dispersed elsewhere under low-advance coefficient. More specifically, it is in the lower-frequency band with a characteristic frequency 0.6
fR approximately, which is marked by red circle box in
Figure 11. This means that there are low-frequency interference signals in the flow field at low advance coefficient, which may be induced by the periodic shedding vortex generated by the complex flow field under low conditions. Under higher conditions which are closer to the design condition, the flow field will be more uniform, thus the low-frequency signals may be covered by the main frequency.
The time history and frequency domain result of the lateral force of the pod is shown in
Figure 12. Similar to the thrust of the propeller, the lateral force of the pod is also vibrating. Meanwhile, all the values have the same sign, which means that the pod is always pushed to one side.
Figure 13 shows the constrained streamlines and the distribution of pressure coefficients on X = ±0.35D planes under
J = 0.8217 at time step
t0 and
t0 + 15Δ
t. The view is from the top of the pod to the bottom. It is obvious that although the strut and fin at the bottom are designed by some symmetrical airfoils, the pressure distribution is asymmetric. This is evidently attributed to the nonuniform outflow by the propeller. On the X = ±0.35D planes, the local flow field around the strut and the fin is obviously affected by the lateral outflow, which is induced by the motion of the propeller. The leading edge of the strut being impacted by the water flow will firstly generate a high-pressure region, while the other side will generate a low-pressure regio, because of the faster velocity induced by the larger distance at the same time. This apparent misalignment presents in both the strut and the fin simultaneously.
What is more, the directivity of the pressure difference is complicated, the details about the pressure coefficient distributions on the airfoil section lines of the strut and the fin (shown in
Figure 13b,d) are shown in
Figure 14. It should be noted that the view observed from the pod is from back to front, thus the side names are defined correspondingly. It is significant that the direction of the local lateral forces suffered by the strut and the fin are inconsistent. Both point to one side firstly, and then point to another side. The turning points in the section lines of the strut and fin are located at 0.37 and 0.7 times the chord length, respectively. In addition, the pressure coefficient fluctuates at the back of the strut, while it does not fluctuate at the same location of the fin. The reason may be that the strut receives the impact of the propeller wake flow earlier than the fin, and the flow filed has a higher pulsation level.
As a comparison, the pressure coefficient distributions on the same location in
Figure 14 at another time step
t0 are plotted, as shown in
Figure 15. In the aspect of the strut, both of the two time steps have similar distribution and turning points. However, the pressure fluctuation is different in the region marked by the blue circle box, where alternating phenomena presents. In other words, with the motion of the propeller, pressure changes periodically in the region mentioned above may lead to vortex shedding. As for the fin, there seems no obvious difference between the two time steps, which is to say that the influence of the wake flow field of the propeller on the flow around the fin is not very obvious. At another location, where on X = ±0.45D planes, with the time steps mentioned above, the phenomenon of pressure fluctuation becomes more violent, as shown in
Figure 16. The pressure on the pod not only fluctuates at different time steps, but also on the left and right sides at the same time step. As for the fin, there are also pressure fluctuations at different time steps, and the location of the turning points are different, which is almost at the same location in X = ±0.35D planes. Although the distributions of local pressure in the strut and the fin are complex, the overall lateral force points to the right side of the pod, which is identified in
Figure 12. Since the lateral force is obtained by integrating the local pressure with the area, and the area of the strut is much larger than the fin. In addition, the lateral force of the pod has the same discrete frequencies with the propeller, which is the result of the pod acting on the propeller wake flow field.
3.2. Characteristics of Pressure Pulsation
In order to investigate the characteristics of the podded propulsor, some monitors of the pressure pulsation have been set, as is shown in
Figure 17.
To investigate the pressure coefficient distributions in the outside and rear of the propeller, some cylindrical surfaces are defined with different diameters 0.7D, 0.9D, 1.2D, 1.5D. As shown in
Figure 18, for a clearer display, the cylindrical surfaces are extended by the phase angle θ, which begins at positive direction of X, then towards to the positive direction of Y. It should also be pointed out that all the surfaces do not intersect with the propeller region, which is why there are some blank regions on surfaces with diameters of 0.7D and 0.9D. Compared with
Figure 13, it can be seen that the region with a phase angle between 180° and 360° is on the left side of the pod, since there are high pressure regions on the certain side of the strut and the fin, as shown in
Figure 18a,e,i. In three different conditions, the pressure coefficient distributions are similar. The pressure pulsation is developed rotationally during the evolution from the region behind the strut towards the left side of the pod, while it is behind the fin towards to the right side of the pod. However, the levels of the pressure coefficient are obviously different. Under the condition
J = 0.6, the levels are the highest, followed by 0.8217 and 1.0.
Figure 19 shows a diagram of a velocity triangle to present the relationship between the advance coefficient and attack angle of propeller airfoil. Under the low-advance coefficient condition, the inflow velocity is lower, and the circumferential velocity is constant on a certain radius, thus the tangent angle between the above two velocities is smaller. Then, the pitch angle is also a certain constant, which means the attack angle, calculated by the pitch angle minus the tangent angle (ignore the influence of induced velocity), is larger than the high advance coefficient condition. According to the theory of airfoils, pressure coefficients increase gradually within a certain range of positive attack angle. Thus, the levels of the pressure field with low-advance coefficients are higher than when under a high advance coefficient.
As shown in
Figure 18a,c, some typical pressure coefficient distribution regions have been named and marked by red rectangular boxes. Other locations and conditions have similar typical regions, but they are not marked in
Figure 18. The region β, which is distributed on both sides of the pod, is induced by the propeller wake flow directly. The region α and γ are behind the strut and the fin, respectively, they are wake flow fields of the strut and the fin. However, these two regions are not similar to normal flow around airfoil, they are not symmetric. The characteristics of the region α and γ are closer to the region β. Obviously, on the surface with diameter 0.7D, due to the influence of the propeller wake flow, which is rotational and pulsating, the characteristics of wake flow fields of the strut and the fin have changed, and the structures of distribution tend to be similar to the propeller wake flow field. With the increase of the surface diameter, the pressure pulsation becomes more violent, as shown in
Figure 18b. The pressure coefficient distribution on this surface is caused by the wake flow generated by the tip of the propeller. According to
Figure 10, it can be seen that the pressure gradient at the tip of the propeller is very evident, which will lead to the higher velocity pulsation, thus the flow from here to the back will have a higher pulsation level. On the surface with diameter 1.2D, the typical distributions are changed. The part beyond the surfaces with diameter 0.7D and 0.9D is the outside of the propeller, where the region δ is located. This periodic phenomenon of alternating positive and negative pressures is caused by the water sucking and drainage of the propeller. In addition, regions α, β, and γ are no longer similar to the regions on the surface with diameters of 0.7D and 0.9D, and become more symmetric, which are close to the characteristics of the symmetric airfoil around flow. On the surface with the largest diameter 1.5D, the distributions are similar to the surface with a diameter of 1.2D, but the levels of pressure coefficients are decreased. According to the analysis above, it can be found that the propeller wake flow has a great impact on the flow field near the pod. However, such an influence is limited and bounded. As shown in
Figure 18, the influence mainly exists in the regions within the propeller diameter. In other regions, the flow field presents similar characteristics to normal flow around the airfoil.
To obtain the attenuation law of pressure pulsation behind the podded propulsor, the rotating frequency and blade passage frequency data under
J = 0.8217 at monitors (see
Figure 17) L1~12, M1~12, N1~12, are extracted, as shown in
Figure 20. It can be seen that the amplitudes of the pressure coefficients are extremely outstanding at the direction of 0° and 180°, where the places behind the fin and the strut are, respectively. Furthermore, the amplitudes at the direction of the strut are larger than the fin, which means that the influence of the strut to the flow field is greater than the fin. The amplitudes of pressure coefficients at other directions are comparatively close. This means that the pod has a significant impact on the propeller wake flow field, which is reflected in the pressure distributions. In addition, the energy contributions in the rotating frequency and blade passage frequency are different at various locations. With the increase of the distance between the monitors and the pod, the energy proportion in the rotating frequency gradually increases, while the blade passage frequency gradually decreases. The overall pressure energy at these monitors must be decreased, it means that the energy in the blade passage frequency decays faster than the energy in the rotating frequency.
3.3. Vortex Structure and Evolution
To investigate the structure of the vortex that is generated from the flow field of the podded propulsor, a method named
Q criterion is utilized to calculate the vortex core regions and present them as iso-surfaces.
Q criterion can be described as [
25]:
where Ω is the rotation rate tensor and
S is the strain rate tensor. ||·||
F means the Frobenius norm, which is the square root of the sum of the squares of all tensor elements. If
Qc > 0, the rotating force will be the predominant role and the vortex will be generated in the flow field. What is more, the contours on the iso-surfaces are obtained by the z-vorticity, which is defined as [
26,
27]:
where
vx and
vy are the velocity components. The reason why vorticity at the direction of
Z is used is that the flow direction is
Z, hence z-vorticity is the best value which can reflect the characteristics of the vortex. The positive and negative values of z-vorticity indicate a different swirl direction. In addition, the vorticity magnitude
ω is used to present the shape of the vortex, which is defined as:
where
ωX and
ωY have similar definitions as
ωZ.
Figure 21 shows the vortex structures of the podded propulsor, which are visualized by iso-surfaces with
Qc = 500 1/s
2 and colored by z-vorticity with value range −100 1/s to 100 1/s. Overall, the vortex structures under three different advance coefficients are similar, but it is stronger under
J = 0.6 which can be sure according to the larger area of high vorticity region shown in
Figure 22. As is discussed in
Section 3.2, with the lower advance coefficient, the larger attack angle will be induced, and the phenomenon of flow separation will be more obvious, thus the vortex will be stronger.
Unlike a traditional single propeller system, the vortex of the podded propulsor has more characteristics. Some typical vortex structures have been marked by black rectangular boxes, as shown in
Figure 21. The tip vortex (TV) is the first typical structure that occurs in the flow field, the root vortex (RV) also appears at the same time. Both are induced by the differential pressure between the pressure side and the suction side. According to the scale calculation in
Figure 21c, the distance between two vortex circles correspond to the pitches at 0.9R and 0.5R, respectively, which are consistent with their generating positions on the blades. The tip vortex can be observed under all the three conditions, but a sspecial type of structure that occurs under
J = 0.6, is called the leapfrogging vortex (LV), shown in
Figure 21a,b and
Figure 22b. This kind of vortex has been discussed extensively in wind turbines [
28]. Felli et al. [
29] also found it in underwater propellers in their experiments. When the adjacent tip vortices pass within a distance, the mutual inductance effect appears, so that they will experience the influence of each other. This effect causes the adjacent tip vortices to roll up around each other, which lead to the leapfrogging phenomenon. With the advance coefficient decreases, the leapfrogging vortex occurs earlier. It is only observed in
Figure 21a,b and
Figure 22b, and close to the propeller. The influence of the pod on the propeller wake flow is not violent in this region, hence the leapfrogging vortex under low advance coefficient will appear. However, with the evolution of the propeller wake flow, the influence of the pod increases, which leads to more complicated vortex structures and cannot find leapfrogging vortex anymore. As shown in
Figure 22, only on Z = 0.35D plane under
J = 0.6 can it be observed, because the leapfrogging vortex has an obvious pitch that is different to the pitch of the tip vortex.
As for the root vortex, under
J = 0.6 and 0.8217, it is wound around the pod, and develops towards the back with a spiral shape. However, it will not develop to the rear of the pod and dissipate quickly due to the effect of the viscous force. As shown in
Figure 22, the root vortex can be found on Z = 0.2D to 1D planes under
J = 0.6 and 0.8217. However, under
J = 1.0, the root vortex can hardly be observed under iso-surface with
Qc = 500 1/s
2, as shown in
Figure 21e and
Figure 22i, which means that the root vortex here is much weaker than the tip vortex.
With the evolution of the tip vortex, it impacts the leading edge of the strut, and be separated to left part and right part. Besides, it can also be observed that there are some cracked and detached vortex (DV), as is shown in
Figure 21c,d and
Figure 22b. They are not connected to the tip vortex, and also has a certain pitch, but it does not match to any section pitch of the propeller blade, which can be observed in
Figure 22b clearly. It may be induced by the energy separation because of the impact.
With the evolution continuing, the tip vortex flows around and through the pod and the fin early or late, which lead to the appearance of the strut vortex (SV) and the fin vortex (FV).
Figure 23 shows the distributions of z-vorticity on different planes. During this evolution process, part of the energy of the tip vortex is used to generate two new vortex structures mentioned above, thus changing its original evolution state. As shown in
Figure 23, the tip vortex attaches the trailing edge of the strut and fin gradually, and the bound vortex occurs on the plane X = −0.35D and 0.45D. It has been known that there is phenomenon of the pressure fluctuation located in where the bound vortices appear, as shown in
Figure 15 and
Figure 16. Then, the bound vortices develop toward to back, and the strut vortex and the fin vortex are generated. This means that the guess mentioned in
Section 3.1 is right: pressure fluctuation is one of the reasons of the vortex generated.
Under all the three conditions, the left part of the tip vortex offsets upward, the blank region generated by offset is occupied by the fin vortex, shown in
Figure 21a,c,e. Correspondingly, the right part of the tip vortex offsets downward, and the region is occupied by the strut vortex, shown in
Figure 21b,d,f. This phenomenon is different from the traditional single propeller system, in which the tip vortex should be parallel or contracted and develop toward the back. Although the two vortex structures mentioned above occupy the region generated by the offset of the different parts of the tip vortex, they are connected with another part of the tip vortex, respectively, as shown in
Figure 22h,l. However, as shown in
Figure 22d, the structure behind the pod is different. The strut vortex seems no connected with the left part of the tip vortex, but it is connected with the hub vortex (HV). Besides, under
J = 0.6, it is observed that there is a merged vortex (MV), which is also induced by the mutual inductance effect. Both of the two phenomena mentioned above can prove the instability of the flow field under a low-advance coefficient.
At the rear of the pod, the hub vortex occurs, and two parts of the tip vortex offset continuously. Because of the lack of grids refined in the location where the hub vortex occurs, there are no more details about the hub vortex and its mutual inductance effect with other vortex structures.