Experimental Investigation of Propeller Wake Velocity Field to Determine the Major Factors Affecting Propeller Wake Wash

The propeller jet from a ship has a significant component directed upwards towards the free surface of the water, which can be used for ice management. This paper describes a comprehensive laboratory experiment where the operational factors affecting a propeller wake velocity field were investigated. The experiment was conducted using a steady wake field to investigate the characteristics of the axial velocity of the fluid in the wake and the corresponding variability downstream of the propeller. The axial velocities and the variability recorded were time-averaged. Propeller rotational speed was found to be the most significant factor, followed by propeller inclination. The experimental results also provide some idea about the change of the patterns of the mean axial velocity distribution against the factors considered for the test throughout the effective wake field, as well as the relationships to predict the axial velocity for known factors.


Introduction
Propeller wake wash can be used as a means of clearing pack ice and small ice masses from the vicinity of offshore installations and shipping channels.In a full-scale study of various pack ice management techniques performed offshore Newfoundland, propeller wake wash was found to have an appreciable degree of success [1].In more recent studies of ice management techniques, the wake of an azimuth thruster was shown to be useful for ice breaking [2,3], in addition to clearing ice from designated areas.Propeller wake wash was found to be an effective means of breaking ice ridges, and pushing away medium sized ice floes and even glacial ice.This technique enables ice management without having any physical interaction with the ice.
The interaction between an ice mass and a propeller jet depends on the jet velocity at a particular location downstream of the propeller, which depends on the propeller shaft rotational speed (n), propeller inclination (θ), and propeller submergence from the fluid surface (H).The application of propeller wake wash as an ice management tool has been investigated at a model tank scale by Ferrieri [4], where change in ice concentration was investigated as a function of the factors affecting propeller wake wash.Bastin [5] developed a simple mathematical model of propeller wake wash based on Ferrieri's experimental results.
A detailed model-scale experimental study on a propeller wake velocity field in the vicinity of the propeller was performed in the large circulating water channel of CNR-INSEAN using LDV (Laser Doppler Velocimetry) measurement techniques [6].Islam [7] investigated the kinematics of propeller wake wash for ice management by using podded propulsors (both puller and pusher arrangements) at different shaft speeds, pod configurations, and propeller inclinations.
Brewster [8] reported that a propeller jet velocity decays exponentially as it moves downstream from the propeller plane or efflux plane.Lam et al. [9,10] reported that the tangential velocity is already significantly diminished a distance x/D p = 3.68 downstream from the propeller.McGarvey [11] and Brewster [8] reported that the radial velocity decays by about 80% within x/D p = 0.30 from the efflux plane.In the present study, only the axial component of fluid velocity is considered, which is expected to be significant throughout the wake field and important for ice management.

Methodology
The experiment was designed using the center composite design (CCD) technique of response surface methodology for three factors (propeller-rotational speed, inclination and submergence) tested at five levels.A total of 18 combinations of the three factors was considered, and measurements were taken at 168 locations in the wake field, thereby requiring 3024 measurements for each response variable.The length of the wake field covered in the experiment was x/D p = 30.5, which was divided into three zones: near field (x/D p ≤ 3.5), intermediate field (3.5 ≤ x/D p ≤ 15.5) and far field (15.5 ≤ x/D p ≤ 30.5).The response variable of interest was the mean axial flow velocity.

Design of Experiment
The required minimum number of factorial combinations for the experimental design was calculated as = {2 3 + (2 × 3) + 1} = 15.There are eight factorial runs, six axial point runs and one center point runs.Three extra center point runs (combinations of factors) were added to the experiment to ensure an accurate estimation of experimental error [12].Design Expert 10.0 software [13] was used to obtain the combinations of factors and to analyze the results.The values of factors used at five levels during the experiment are listed in Table 1.The lowest propeller rotational speed was 6.0 rps, which was determined on the basis of the formula for Reynolds numbers (Propeller Reynolds number 'Re prop ', and Flow Reynolds number 'Re flow ') to meet the requirement for minimum Reynolds number, such that the viscous scaling effect on the experimental results can be neglected [14].

Experimental Set-Up
The experiment was done in a tow tank (Length × Breadth × Depth = 58.27× 4.57 × 3.04 m).The arrangement of equipment used in the experiment is shown in Figure 1.The major components were the open boat propulsion unit, the EMCM (Electro-magnetic Current Meter) sensors, and the data acquisition system.
A B-series propeller of diameter 250 mm was located at the center of the tank.The wake survey was carried out in the bollard pull condition.A stock open boat propulsion unit was used to drive the propeller.The measurements were taken in terms of the center line of the tank at half the width of the wake field, where the propeller was always positioned along the center line of tank.

Scaling Effect
The variables of the prediction model were non-dimensionalized using 'Buckingham-π theory'.The viscous scaling effects of the experimental results were checked using the empirical formula developed by Verhey [14], who suggested that the scaling effect due to viscosity is negligible, if the Reynolds number of the propeller (Reprop) and the Reynolds number of the flow (Reflow) are greater than 7 × 10 4 and 3 × 10 3 , respectively, for a propeller wake.Verhey used the formulae: Reflow = ν and Reprop = , where V0 = efflux velocity of the propeller jet, and V = 1.59nD √C [15], Ct is the thrust coefficient of model propeller, DP = Propeller diameter, n = propeller rotational speed (rps), ν = kinematic viscosity of fluid (water) = 1.0023 × 10 −6 m 2 /s at 20 °C; Lm = length term dependent on blade area ratio (β), number of blades of propeller (N), diameter of hub (Dh) and DP.Lm is defined [16] as: L = β.D .π. 2N 1 − . According to Verhey's formulas, the Reynolds numbers for the lowest shaft rotational speed (6 rps) of the model propeller were Reflow = 2.87 × 10 5 and Reprop = 8.54 × 10 4 .

Open Boat Propulsion System
The open boat propulsion system integrates all the equipment that was used to generate the propeller wake under a particular propulsion condition.A total of 18 experimental conditions were used to develop the prediction models.The unit consisted of the propeller, the open boat propulsion unit, and the frame to locate the open boat unit on the carriage of the tow tank.The propeller was a B4-55 series of diameter Ø = 250 mm.The frame had a mechanism of adjusting the inclination of the propeller.The power delivered was adjustable which was synchronized to the shaft rotational speed of propeller.Table 2 illustrates the properties of the prototype propeller.

EMCM Equipment
An EMCM (Electro-magnetic Current Meter) system was used to measure the wake flow velocity during the experiment.This system included an EMCM package (EMCM sensors and a built-in mounting to maintain the orientation of sensors) to measure the axial flow velocity, a frame to mount and operate the sensors, a linear displacement transducer to measure the transverse

Scaling Effect
The variables of the prediction model were non-dimensionalized using 'Buckingham-π theory'.The viscous scaling effects of the experimental results were checked using the empirical formula developed by Verhey [14], who suggested that the scaling effect due to viscosity is negligible, if the Reynolds number of the propeller (Re prop ) and the Reynolds number of the flow (Re flow ) are greater than 7 × 10 4 and 3 × 10 3 , respectively, for a propeller wake.Verhey used the formulae: Re flow = V 0 D P ν and Re prop = nL m D P ν , where V 0 = efflux velocity of the propeller jet, and V 0 = 1.59nDP √ C t [15], C t is the thrust coefficient of model propeller, D P = Propeller diameter, n = propeller rotational speed (rps), ν = kinematic viscosity of fluid (water) = 1.0023 × 10 −6 m 2 /s at 20 • C; L m = length term dependent on blade area ratio (β), number of blades of propeller (N), diameter of hub (D h ) and D P .
L m is defined [16] as . According to Verhey's formulas, the Reynolds numbers for the lowest shaft rotational speed (6 rps) of the model propeller were Re flow = 2.87 × 10 5 and Re prop = 8.54 × 10 4 .

Open Boat Propulsion System
The open boat propulsion system integrates all the equipment that was used to generate the propeller wake under a particular propulsion condition.A total of 18 experimental conditions were used to develop the prediction models.The unit consisted of the propeller, the open boat propulsion unit, and the frame to locate the open boat unit on the carriage of the tow tank.The propeller was a B4-55 series of diameter Ø = 250 mm.The frame had a mechanism of adjusting the inclination of the propeller.The power delivered was adjustable which was synchronized to the shaft rotational speed of propeller.Table 2 illustrates the properties of the prototype propeller.

EMCM Equipment
An EMCM (Electro-magnetic Current Meter) system was used to measure the wake flow velocity during the experiment.This system included an EMCM package (EMCM sensors and a built-in mounting to maintain the orientation of sensors) to measure the axial flow velocity, a frame to mount and operate the sensors, a linear displacement transducer to measure the transverse displacement of the sensors, a platform to support the system equipment, and connecting wires.The tank platform was installed across the breadth of the tank (4.57m), on the top of which the EMCM equipment with the supporting frame, and the data acquisition system were installed.The EMCM sensors could slide on the frame with respect to a datum.The datum was marked at the center of the wake field (y = 0).The readings along the transverse direction of the propeller wake (along y/D p ) were collected by sliding the sensors following a pre-marked scale on the frame.The EMCM system was connected with the data acquisition system, where all the data were recorded and stored immediately after the capture of the reading.In the experiment, each measurement was taken for a period of at least 30 s at 4 Hz.

Results and Discussion
The objective of this study was to determine the influence of the three input factors (Table 1) affecting the propeller wake wash.A total of 18 propeller conditions (experimental runs) comprising different combinations of the three factors were used to create an experimental plan with the Design Expert 10.0 software for a Center Composite Design (CCD) of the Response Surface Methodology.The resulting data were used to develop prediction equations for the mean axial velocity of flow in the wake.Flow velocity measurements were taken longitudinally and transversely covering half-width of the wake field for three water depths (d = 0.25D p , 0.55D p and 1.05D p ).A decision was made to assume longitudinal symmetry about the mid plane and thus only to survey half the width of wake field based on empirical analyses performed by earlier researchers [8,9,11,14,[16][17][18][19][20], who reported the axisymmetric nature of the wake velocity distribution for the far zone downstream of propeller wake.This allowed a more detailed survey of the flow within the limitations of time and resources.

Individual Effect of Factors
Initially the influences of individual factors on the mean axial velocity of flow were assessed by applying an analysis approach known as 'One Factor at a Time (OFAT)'.In OFAT, the configurations of factors were selected such that only one factor was varied at a time, keeping the other two factors unchanged.Three cases were considered illustrating the individual effect of each factor.The combinations of factors in each case are shown in Table 3.The above cases illustrate the effect of each factor on the mean axial flow velocity along the normalized dimension x/D p (longitudinally downstream of propeller).The transverse location for these data was along the center of the wake field (y/D p = 0), where the tank wall boundary effect was the minimum.

Individual Effect of Factors on Mean Axial Flow Velocity
Case #1: The individual effect of propeller shaft rotational speed (n) on the mean axial velocity of flow (V a ) along x/D p at the center of the wake field (y/D p = 0) is illustrated in Figure 2a-c for the three different depths of measurements: 0.25D p , 0.55D p and 1.05D p , respectively.These illustrate that the mean axial velocity along the centerline generally increases with 'n'.
The above cases illustrate the effect of each factor on the mean axial flow velocity along the normalized dimension x/Dp (longitudinally downstream of propeller).The transverse location for these data was along the center of the wake field (y/Dp = 0), where the tank wall boundary effect was the minimum.

Individual Effect of Factors on Mean Axial Flow Velocity
Case #1: The individual effect of propeller shaft rotational speed (n) on the mean axial velocity of flow (Va) along x/Dp at the center of the wake field (y/Dp = 0) is illustrated in Figure 2a-c for the three different depths of measurements: 0.25Dp, 0.55Dp and 1.05Dp, respectively.These illustrate that the mean axial velocity along the centerline generally increases with 'n'.Case #2: The individual effect of propeller inclination (θ) on the mean axial velocity of flow (Va) along x/Dp at the center of the wake field (y/Dp = 0) is illustrated in Figure 3a-c for the three depths of measurements: 0.25Dp, 0.55Dp and 1.05Dp, respectively.These illustrate that the characteristics of the mean axial velocity along x/Dp with 'θ' depend on the depth.In general, at shallow depth, the mean axial velocity increases with 'θ', and the effect of 'θ' reduces with increasing depth.The effect of propeller inclination on the mean axial velocity is most pronounced near the surface and negligible at the lowest depth.Moreover, the effect of 'θ' is not as strong as 'n'.This phenomenon becomes clearer in 3D surface plots in the following section.Case #2: The individual effect of propeller inclination (θ) on the mean axial velocity of flow (V a ) along x/D p at the center of the wake field (y/D p = 0) is illustrated in Figure 3a-c for the three depths of measurements: 0.25D p , 0.55D p and 1.05D p , respectively.These illustrate that the characteristics of the mean axial velocity along x/D p with 'θ' depend on the depth.In general, at shallow depth, the mean axial velocity increases with 'θ', and the effect of 'θ' reduces with increasing depth.The effect of propeller inclination on the mean axial velocity is most pronounced near the surface and negligible at the lowest depth.Moreover, the effect of 'θ' is not as strong as 'n'.This phenomenon becomes clearer in 3D surface plots in the following section.
Case #3: The individual effect of propeller submergence (H) on the mean axial flow velocity (V a ) along x/D p at the center of the wake field (y/D p = 0) is illustrated in Figure 4a-c for the three different depths of measurements: 0.25D p , 0.55D p and 1.05D p , respectively.These plots illustrate that the characteristics of the mean axial velocity along x/D p with 'H' also largely depend on the depth.In general, at shallow depth, the axial velocity increases with decreasing 'H', and at deeper depth the mean axial velocity is not significantly affected in the downstream and shows mixed effects in the region near the propeller.Case #3: The individual effect of propeller submergence (H) on the mean axial flow velocity (Va) along x/Dp at the center of the wake field (y/Dp = 0) is illustrated in Figure 4a-c for the three different depths of measurements: 0.25Dp, 0.55Dp and 1.05Dp, respectively.These plots illustrate that the characteristics of the mean axial velocity along x/Dp with 'H' also largely depend on the depth.In general, at shallow depth, the axial velocity increases with decreasing 'H', and at deeper depth the mean axial velocity is not significantly affected in the downstream and shows mixed effects in the region near the propeller.
(a) (b) Case #3: The individual effect of propeller submergence (H) on the mean axial flow velocity (Va) along x/Dp at the center of the wake field (y/Dp = 0) is illustrated in Figure 4a-c for the three different depths of measurements: 0.25Dp, 0.55Dp and 1.05Dp, respectively.These plots illustrate that the characteristics of the mean axial velocity along x/Dp with 'H' also largely depend on the depth.In general, at shallow depth, the axial velocity increases with decreasing 'H', and at deeper depth the mean axial velocity is not significantly affected in the downstream and shows mixed effects in the region near the propeller.Given that the focus of this study is on the use of wake flow in ice management, the wake flow at the surface is more significant in terms of the effect on floating ice.Thus, the surface flow characteristics were analyzed in more detail.
Case #1: The individual effect of propeller shaft rotational speed (n) on the mean axial velocity of flow along y/Dp at various longitudinal locations (x/Dp) is illustrated in Figure 5, for the depth of measurement closest to the water surface, which was d = 0.25Dp.The two longitudinal locations closest to the propeller (at x/Dp = 0.50 and 1.50) are ignored due to the confused flow in those regions.The readings taken at the locations close to the tank wall were affected by the boundary condition and so are ignored for developing the curves within the plots (Figures 5-7).Therefore, the curves illustrating the distribution of the mean axial velocities along y/Dp are the approximate pattern of the velocity distribution along y/DP.Figure 5 (below) illustrates that the mean axial velocity along y/Dp also increases with 'n'.In addition, the mean axial velocity along transverse locations is the highest at x/DP = 7.50 for all three 'n' considered in the study, where θ = 4.5° and H = 350 mm.
Case #2: The individual effect of propeller inclination (θ) on the mean axial velocity of flow (Va) along y/Dp at various longitudinal locations (x/Dp) is illustrated in Figure 6, for the depth closest to the water surface.Figure 6 shows that, the mean axial velocity along y/Dp generally increases with 'θ', but the changes are different at different locations.Also, the changes are not as significant as those due to 'n'.The locations of the highest velocities also depend on 'θ'.At θ = 0°, the most effective zone (the overall mean axial velocities in the zone are higher than other zones) is approximately at x/Dp = 11.50.At θ = 4.5°, the most effective zone is approximately at x/Dp = 7.50.At θ = 9°, the most effective zone is approximately from x/Dp = 11.5 to 15.5.
Case #3: The individual effect of propeller submergence (H) on the mean axial velocity (Va) along y/Dp at various longitudinal locations (x/Dp) is illustrated in Figure 7 for the depth of measurement closest to the water surface (at 0.25Dp).The effect of 'H' on the mean axial velocity is opposite to the other factors, particularly in the near field (x/Dp ≤ 3.5) and intermediate field (3.5 ≤ x/Dp ≤ 15.5) zones.With a decrease of 'H', the mean axial velocity increases at any location in the near field and the intermediate field zones, whereas the effect of 'H' is inconsistent in the far field zone.

Individual Effect of Factors on the Mean Axial Velocity of Flow along y/D p at 0.25D p
Given that the focus of this study is on the use of wake flow in ice management, the wake flow at the surface is more significant in terms of the effect on floating ice.Thus, the surface flow characteristics were analyzed in more detail.
Case #1: The individual effect of propeller shaft rotational speed (n) on the mean axial velocity of flow V a D P ν along y/D p at various longitudinal locations (x/D p ) is illustrated in Figure 5, for the depth of measurement closest to the water surface, which was d = 0.25D p .The two longitudinal locations closest to the propeller (at x/D p = 0.50 and 1.50) are ignored due to the confused flow in those regions.The readings taken at the locations close to the tank wall were affected by the boundary condition and so are ignored for developing the curves within the plots (Figures 5-7).Therefore, the curves illustrating the distribution of the mean axial velocities along y/D p are the approximate pattern of the velocity distribution along y/D P .
Figure 5 (below) illustrates that the mean axial velocity along y/D p also increases with 'n'.In addition, the mean axial velocity along transverse locations is the highest at x/D P = 7.50 for all three 'n' considered in the study, where θ = 4.5 • and H = 350 mm.
Case #2: The individual effect of propeller inclination (θ) on the mean axial velocity of flow (V a ) along y/D p at various longitudinal locations (x/D p ) is illustrated in Figure 6, for the depth closest to the water surface.Figure 6 shows that, the mean axial velocity along y/D p generally increases with 'θ', but the changes are different at different locations.Also, the changes are not as significant as those due to 'n'.The locations of the highest velocities also depend on 'θ'.At θ = 0 • , the most effective zone (the overall mean axial velocities in the zone are higher than other zones) is approximately at x/D p = 11.50.At θ = 4.5 • , the most effective zone is approximately at x/D p = 7.50.At θ = 9 • , the most effective zone is approximately from x/D p = 11.5 to 15.5.
Case #3: The individual effect of propeller submergence (H) on the mean axial velocity (V a ) along y/D p at various longitudinal locations (x/D p ) is illustrated in Figure 7 for the depth of measurement closest to the water surface (at 0.25D p ).The effect of 'H' on the mean axial velocity is opposite to the other factors, particularly in the near field (x/D p ≤ 3.5) and intermediate field (3.5 ≤ x/D p ≤ 15.5) zones.With a decrease of 'H', the mean axial velocity increases at any location in the near field and the intermediate field zones, whereas the effect of 'H' is inconsistent in the far field zone.
The effects of each of the three factors on the response variable are illustrated explicitly from the study, and the interaction effect of the factors can be shown through a 3D plot or a regression equation.The effects of each of the three factors on the response variable are illustrated explicitly from the study, and the interaction effect of the factors can be shown through a 3D plot or a regression equation.

Interaction Effect of Factors
The interaction effect of multiple factors can be represented by regression equations or 3D contour/surface plots, which show the relationship between factors and the effect on the response variable for each location.In this section, 3D surface plots showing the two-factorial-interaction (2FI) effect of the most influential factors (filtered by the Design Expert software) are presented for locations along the center of the wake field at the depth of measurement closest to the fluid surface (d = 0.25Dp).The locations that are selected for the investigation are at x/Dp = 7.5, 11.5 and 15.5 at d = 0.25Dp, where the overall mean axial velocity was higher than other zones.Additional plots are provided in [21], which is the basis of this paper.Figure 8a-c show the change of the non-dimensionalized mean axial velocity with the change of both 'n' and 'θ' along the center of the wake field at x/Dp = 7.5, 11.5 and 15.5 at d = 0.25Dp respectively.
These plots show that, at x/Dp = 7.5 the velocity increases with 'n' and 'θ', but the effect of 'n' is larger than 'θ'.The effect of 'θ' gradually diminishes and the effect of 'n' decreases with x/Dp.Also it is noticeable that the effect of 'n' in comparison to the effect of 'θ' is greater at all three locations, and it is expected that this trend will be applicable for all locations downstream in the wake field.
In Figure 8a-c, all the parameters are non-dimensionalized.The typical criteria used by the Design Expert 10.0 software for ANOVA testing are: (i) Normality: whether the distributions of the residuals are normal.
(ii) Constant variance: the variance of data in groups should be the same.
(iii) Independence: all samples are independent of one another.
(iv) Randomness: all samples are random samples from their respective populations.

Interaction Effect of Factors
The interaction effect of multiple factors can be represented by regression equations or 3D contour/surface plots, which show the relationship between factors and the effect on the response variable for each location.In this section, 3D surface plots showing the two-factorial-interaction (2FI) effect of the most influential factors (filtered by the Design Expert software) are presented for locations along the center of the wake field at the depth of measurement closest to the fluid surface (d = 0.25D p ).The locations that are selected for the investigation are at x/D p = 7.5, 11.5 and 15.5 at d = 0.25D p , where the overall mean axial velocity was higher than other zones.Additional plots are provided in [21], which is the basis of this paper.Figure 8a-c show the change of the non-dimensionalized mean axial velocity with the change of both 'n' and 'θ' along the center of the wake field at x/D p = 7.5, 11.5 and 15.5 at d = 0.25D p respectively.These plots show that, at x/D p = 7.5 the velocity increases with 'n' and 'θ', but the effect of 'n' is larger than 'θ'.The effect of 'θ' gradually diminishes and the effect of 'n' decreases with x/D p .Also it is noticeable that the effect of 'n' in comparison to the effect of 'θ' is greater at all three locations, and it is expected that this trend will be applicable for all locations downstream in the wake field.
In Figure 8a-c, all the parameters are non-dimensionalized.The typical criteria used by the Design Expert 10.0 software for ANOVA testing are: (i).Normality: whether the distributions of the residuals are normal.
(ii).Constant variance: the variance of data in groups should be the same.(iii).Independence: all samples are independent of one another.(iv).Randomness: all samples are random samples from their respective populations.

Axial Velocity Prediction Equations
The functional relationship between the non-dimensionalized response variable (non-dimensionalized by mean axial velocity) and the non-dimensionalized input factors can be expressed as follows: Buckingham-π theory was used to develop the non-dimensionalized functional relationship among the factors: propeller rotational speed 'n', inclination of propeller 'θ', depth of submergence of propeller 'H', and the response variable of the study, which is mean axial velocity of fluid 'V a '.
Assuming: and H D p = Z; the common form of the quadratic equation representing the relationship among the factors and the mean axial velocity can be written as:  Design Expert 10.0 software was used to generate the coefficients of Equation (1) for known propeller shaft rotational speed (n), propeller inclination (θ), propeller submergence (H) and corresponding mean axial flow velocity for all 168 locations.Then, the relationships among the 10 coefficients in terms of x D p , y D p and d D p were determined.To increase the accuracy, the regression equations for predicting the coefficients were developed separately for each zone.
There are 10 unknown coefficients in Equation ( 1), each of which can be derived by using the regression equations listed in Tables 4-6 for each of the three zones of interest.Stepwise regression analysis was used to model the prediction equations to determine the coefficients for up to cubic order, to maintain good prediction capacities of the equations.Hence, for a given location in the 3D wake field, the coefficients of Equation ( 1) can be determined, and for a given propeller condition (where n, θ, and H are known), the mean axial flow velocity can be estimated.

Conclusions
The present experimental investigation on a propeller wake velocity (axial) field for the major factors affecting propeller wake wash performance leads to the following concluding remarks: -Among the three factors, propeller rotational speed is the most influential, followed by propeller inclination angle and propeller depth of submergence.The propeller speed 'n' and inclination

Figure 2 .
Figure 2. (a) Mean axial flow velocity along x/Dp at y/Dp = 0 and d = 0.25Dp for different 'n'; (b) Mean axial flow velocity along x/Dp at y/Dp = 0 and d = 0.55Dp for different 'n'; (c) Mean axial flow velocity along x/Dp at y/Dp = 0 and d = 1.05Dp for different 'n'.

Figure 2 .
Figure 2. (a) Mean axial flow velocity along x/D p at y/D p = 0 and d = 0.25D p for different 'n'; (b) Mean axial flow velocity along x/D p at y/D p = 0 and d = 0.55D p for different 'n'; (c) Mean axial flow velocity along x/D p at y/D p = 0 and d = 1.05D p for different 'n'.

Figure 3 .Figure 3 .
Figure 3. (a) Mean axial flow velocity along x/D p at y/D p = 0 and d = 0.25D p for different 'θ'; (b) mean axial flow velocity along x/D p at y/D p = 0 and d = 0.55D p for different 'θ'; (c) mean axial flow velocity along x/D p at y/D p = 0 and d = 1.05D p for different 'θ'.

Figure 4 .
Figure 4. (a) Mean axial flow velocity along x/D p at y/D p = 0 and at d = 0.25D p for different 'H'; (b) mean axial flow velocity along x/D p at y/D p = 0 and at d = 0.55D p for different 'H'; (c) mean axial flow velocity along x/D p at y/D p = 0 and at d = 1.05D p for different 'H'.

Figure 5 .
Figure 5.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities ν

Figure 6 .
Figure 6.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities along y/Dp with the change of 'θ'.

Figure 5 .
Figure 5.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities V a D P ν along y/D p with the change of 'n'.

Figure 5 .
Figure 5.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities ν

Figure 6 .
Figure 6.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities along y/Dp with the change of 'θ'.Figure 6.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities

Figure 6 .
The approximate patterns of the distribution of the non-dimensionalized mean axial velocities V a D P ν along y/Dp with the change of 'θ'.

Figure 7 .
Figure 7.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities along y/Dp with the change of 'H'.

Figure 7 .
Figure 7.The approximate patterns of the distribution of the non-dimensionalized mean axial velocities V a D P ν along y/Dp with the change of 'H'.

Figure 8 .
Figure 8.(a) 3D surface plot showing the interaction effect of 'n' and 'θ' on the mean axial flow velocity at the point location of (x/D p , y/D p , d/D p ) = (7.5, 0.0, 0.25); (b) 3D surface plot showing the interaction effect of 'n' and 'θ' on the mean axial flow velocity at the point location of (x/D p , y/D p , d/D p ) = (11.5,0.0, 0.25); (c) 3D surface plot showing the interaction effect of 'n' and 'θ' on the mean axial flow velocity at the point location of (x/D p , y/D p , d/D p ) = (15.5,0.0, 0.25).

There are 10
unknown coefficients in the relationship among the factors and the response, values which depend on the location x D the corresponding coefficients were also developed for each zone (near field, intermediate field and far field), which are shown in Tables4-6.The kinematic viscosity of water (ν) is considered as 1.0023 × 10 −6 m 2 /s to calculate the coefficients.

Table 1 .
Factors used at different levels of the experiment.

Table 2 .
Properties of the model propeller.

Table 2 .
Properties of the model propeller.

Table 3 .
Combinations of factors under each case to assess the individual effect of factors.

: Influence of Propeller Depth of Submergence 'H' on the Mean Axial Velocity of Flow
* Values for the Center Point Run are the average of the four center point runs (combinations of factors) considered in the experiment (n = 9 rps, θ = 4.5 • , H = 350 mm).