Abstract
This study uses a combination of computational fluid dynamics (CFD) and experimental research to explore the influence of the length of the inlet duct on the hydraulic performance of the water jet propulsion pump device. By extending the oblique, straight, and long pipe section of the inlet duct, six sets of schemes were designed. The research results show that as the length of the inlet duct increases, the head, thrust, and inlet duct efficiency of the water jet propulsion pump first increase and then decrease, reaching the maximum value in scheme 4, but it has little effect on the efficiency of the entire water jet propulsion system. In addition, after extending the obliquely straight long pipe section of the inlet flow channel, the velocity distribution of the outlet section of the flow channel is gradually uniform, and the absolute value of the vorticity first increases and decreases, but the uniformity of the velocity distribution and the weighted average angle change little. This research enriches the research theory of the inlet duct and also provides a reference for the selection of inlet duct length and its performance optimization.
1. Introduction
Although traditional propeller propulsion has high propulsion efficiency in low and medium speeds, and the propeller propulsion efficiency of displacement ships is about 65% [1], the water jet propulsion pump device has high efficiency and low noise when the ship is traveling at high speed, and it is widely used on high-speed ships [2]. In the past decades, with the rapid development of computer technology, CFD technology has been widely used in many fields, among which a lot of research has been done on the hydraulic performance of pumps in fluid machinery [3–7]. In order to ensure the accuracy of numerical simulation, a large number of studies through the combination of experimental and numerical simulation methods were compared to obtain more reliable research results [8–12]. All kinds of mathematical methods have been widely used, combined with numerical simulation to solve the practical engineering problems [13–16]. The thrust of water jet propulsion is obtained by the reaction force of the water jet from the propulsion pump. One of the main components of the water inlet duct is to transfer the water flow from the bottom of the ship to the propulsion pump. Therefore, the performance of the water inlet duct will affect the performance of the entire propulsion system.
Regarding the parameter optimization of the inlet duct, Jiao et al. [17, 18] analyzed the relationship between the performance of the inlet duct and the ship’s speed and obtained the geometric parameter values under the best optimized condition. Huang et al. [19] improved the NSGA-II and TOPSIS technologies and optimized the performance of the inlet duct with the unevenness and verticality of the outlet and hydraulic efficiency as the optimization goals. Gong et al. [20] designed a model predictive controller for steering system control. Regarding the research on the internal flow characteristics of the inlet duct, Park et al. [21] analyzed the velocity vector and streamline numerical results of the inlet duct as well as the pressure distribution, showed the eddy caused by the separation along the corner of the side wall, and explained the flow phenomenon inside the inlet channel. Ding and Wang [22] used computational fluid dynamics to determine the flow loss of the inlet channel, thereby further improving the prediction accuracy of the propulsion performance of the entire water jet propulsion system. Xu et al. [23] used HSV to show the evolution of TLC and solved the problem of tip leakage cavitation. Cao et al. [24] carried out numerical simulations on the water jet propulsion pump device and revealed the internal relationship between the unevenness of the internal flow of the inlet duct and the performance deviation of the water jet propulsion pump device and essentially studied the performance deviation between the uniform and nonuniform suction flow. Regarding the innovation of the research method of the inlet flow path, Huang and Luo [25] improved the NSGA-II and TOPSIS technologies and optimized the performance of the inlet flow channel with the unevenness and verticality of the outlet and hydraulic efficiency as the optimization goals. Lindau et al. [26] used the RANS method to conduct an in-depth study on the cavitation performance of water jet axial flow pumps. Gong et al. [27] used PIV to measure the velocity distribution. Kimball and Taylor [28], Kerwin [29], etc., based on the lifting surface theory of propeller design and combined with RANS numerical calculations, carried out research on the design method of water jet propulsion pumps and obtained certain law. Gong et al. [30] studied the interaction between the water jet propulsion system and the hull and found that the different design of the inlet passage would lead to different interactions between them, leading to the performance differences among the four water jet propulsion systems. Predecessors have carried out a lot of research work on the inlet duct and the entire water jet propulsion pump device, but there are few aspects of the analysis and optimization of the inlet duct parameters.
In view of this, the performance optimization of the inlet duct is of great significance to the application of the entire water jet propulsion pump device. In order to understand the influence of the length of the inlet duct on the external characteristics of the inlet duct, the uniformity of the outlet of the inlet duct, and the internal characteristics of the inlet duct, numerical simulation and qualitative and quantitative analyses will be carried out in this article. This is conducive to the application of the water jet propulsion pump device on ships and also expands the theoretical analysis content of the inlet duct.
2. Numerical Simulation Methods
2.1. Computational Domain
Because the performance of the inlet channel is inseparable from the hull structure and working conditions, and the performance of the inlet channel will also affect the hydraulic performance of the propulsion pump, the control body should be around the propulsion pump and the water inlet at the bottom of the stern when performing performance analysis on the inlet duct. Figure 1 shows the entire water jet propulsion system, which includes the inlet duct, propulsion pump (mixed flow pump), and nozzles. According to the literature [31], the length, width, and height of the flow field control body are, respectively, 30D, 10D, and 8D (D is the diameter of the outlet of the inlet duct, D = 223.4 mm).

2.2. Mesh Generation and Independence Analysis
Use ICEM software to complete the structured grid of the inlet duct, the propulsion pump, and the water body; the schemes with different grid numbers are shown in Table 1. As the number of grids increases, the accuracy of simulation results gradually improves, but at the same time the demand for computing resources also increases. Domestic and foreign institutions have made clear guidance and definitions for CFD credibility analysis, and most of them use Richardson extrapolation to evaluate errors caused by the calculation domain. Whether it is a high-order or low-order format, the numerical calculation results of the grid will tend to be close to the exact solution [32].
In order to obtain the most suitable number of grids, the method of changing the grid size of the water inlet channel is used to verify the grid independence; that is, as the number of grids increases, the efficiency of the water jet propulsion pump device is calculated until the calculation result is within the allowable range.
As shown in Figure 2, when the grid number of the inlet flow channel is 1.26 million, 1.35 million, and 1.48 million, the efficiency of the pump device hardly changes. The calculation results show that when the number of grids exceeds 1.35 million, the change of pump device efficiency is 0.018% and 0.004%, respectively, which can be regarded as the result has nothing to do with the grid. Therefore, the number of grids of the inlet runner is 1.35 million (as shown in Figure 3).


2.3. Boundary Conditions
This article is a steady calculation. IVR = 0.45 (IVR is the ratio of the ship speed to the averaged axial outflow velocity at the passage outlet). In this paper, the inlet boundary of the water body is set as the speed inlet, and the speed is 8 m/s. The nominal turbulence intensity (value 5%) is used, and the hydraulic diameter is 993 mm. Both the outlet of the water body and the outlet of the nozzle are set to average static pressure, the value is one atmospheric pressure, and the impeller speed is 700 r/min. The frame change/mixing model of the two interfaces related to the impeller is set to frozen rotor. The frame change/mixing model of the other interfaces is set to frozen rotor. The turbulence model selects the standard k–ε turbulence model [33], applies the first-order upwind scheme, and sets the convergence accuracy to 10−5 (as shown in Figure 4).

2.4. Research Case
In order to analyze the influence of the length of the inlet duct on itself and the entire water jet propulsion system, six schemes were designed (as shown in Table 2). Taking option 4 as the original schemes, the length of the inclined straight pipe section is L = 1.12D (D is the diameter of the outlet section of the inlet duct), the inclination angle of the inlet duct is 35°, and the inlet duct height H = Lcos35° + h (h is the vertical distance between the upper end of the inclined straight pipe section and the lower end of the inlet duct exit section) (as shown in Figure 5). When the length L of the inclined straight pipe section is changed, the length of the entire inlet duct will change accordingly. Because the angle is constant, the height of the runner also changes accordingly (as shown in Table 2). Through the analysis of these six schemes, certain rules can be drawn.

2.5. Analysis Parameters
2.5.1. External Characteristic
Efficiency is a comprehensive performance index for evaluating water jet propulsion system, which can well reflect the pros and cons of water jet propulsion system. Therefore, formulas (1)–(3) and formula (4) will be used to calculate the inlet duct efficiency and the efficiency of the entire system, respectively, and will be used as the analysis index to obtain the graph of the change rule of efficiency with the length of the inlet duct, where represents the efficiency of the inlet runner in the fourth scheme.
The formula for calculating the efficiency of the inlet duct is
In the formula, , , and are the inlet duct efficiency, the total energy of the runner inlet, and the total energy of the runner exit, respectively.
Water Jet Propulsion System Efficiency. The physical meaning of system propulsion efficiency is the ratio of the output power of the system to the output power of the propulsion pump [34], and considering the pipeline loss coefficient K1, since (A is the nozzle area), it is brought into
In the formula, is the axial flow velocity of the nozzle; is the speed of the ship; α is the influence coefficient of the boundary layer, generally 0.95; and K1 is the pipeline coefficient, generally 0.45.
2.5.2. Analysis of Flow Uniformity
In order to quantitatively analyze the flow uniformity of the outlet part of the inlet duct (as shown in Figure 6), the objective function velocity distribution uniformity and the weighted average angle θ are hereby selected as the judgment basis. (θ refers to the weighted average of the speed of each node on the selected section).

Speed distribution uniformity is represented aswhere is the flow velocity distribution uniformity of the outlet section of the inlet duct; is the average axial velocity of the inlet section of the impeller chamber; is the axial velocity of each unit at the inlet section of the impeller chamber; and m is the number of units divided by the section during numerical calculation of the flow field.
Weighted average angle of exit velocity distribution [35] iswhere is the weighted average angle of velocity at the exit section of the runner; is the axial velocity of each unit at the inlet section of the impeller chamber; and is the lateral velocity of each unit at the inlet section of the pump.
2.5.3. Thrust Characteristics
The water jet propulsion system relies on the water jet from the nozzle to obtain the reverse thrust to propel the ship. Therefore, the thrust of the water jet propulsion system is a very important indicator. In an ideal state, the outlet flow from the nozzle should not rotate at all, only the axial high-speed flow part, so the velocity distribution of the nozzle is closely related to the thrust of the water jet [36]. Taking the computational domain in this article as the controlling body, according to the momentum theorem, we can get
To this end, the following formula is used to calculate the thrust of the device, and the relative increment is calculated on the basis of Scheme 4.
Thrust of the whole device [37] iswhere A is the cross-sectional area of the nozzle; is the fluid density; is the axial flow velocity of the nozzle; is the flow velocity at the inlet of the inlet duct; is the speed of the ship; and is the boundary layer influence coefficient, generally 0.95.
3. Experimental Verification
In this paper, the model in the literature [38] is selected. In order to verify the reliability of the numerical simulation, the same boundary conditions as in the literature are set to calculate the efficiency of the water jet propulsion pump under the same working conditions. The selected test equipment is shown in Figure 7. A test facility is built which consists of two closed-loop circuits. One circuit is used to provide navigation speed for the propulsion pumping system. The other is applied to measuring performance of the pumping system. The propulsion pumping system consists of a propulsion pump, a guide vane, and an inlet passage. The pump is driven by a DC electric motor via a system. A frequency controller is used to regulate the shaft speed. The flow rate of the pump is measured by an electromagnetic flow meter with an accuracy of ±0.5%. The system uncertainty of the instrument measuring the head is ±0.1% (the uncertainty of the entire experiment includes system uncertainty and random uncertainty).

Figure 8 is a performance comparison between test results and experimental results. The graph uses a series of flow points as the abscissa and head and efficiency as the ordinate to compare the numerical simulation results with the experimental data. It can be seen from the figure that the head decreases with the increase of the flow, and the efficiency increases with the increase of the flow. The experimental results of the water jet propulsion device and the overall trend of the numerical results are consistent, so the numerical simulation results in this article are credible.

4. Performance Analysis
In order to fully analyze the influence of the length of the inlet duct on the inlet duct and the entire water jet propulsion system, the following will present a quantitative and qualitative analysis from three aspects: the external characteristics of the inlet duct, the uniformity of the outlet flow, and the internal flow of the runner. A certain rule is found to facilitate the application of runners at different lengths.
4.1. External Characteristic Analysis
Figure 9 shows the variation of the head H and system efficiency with the length of the inlet flow channel. H0 and ηs0 are the head and system efficiency of the original scheme. It can be seen from the figure that the efficiency of the system does not change much with the length of the inlet duct, and the variation range is within 0.5%, indicating that the entire water jet propulsion pump device is little affected by the inlet duct length. The head of the water-jet propulsion pump varies greatly by the length of the inlet duct. The head of Case 3 is the lowest, which is reduced by 2.7%, and the head of Case 4 is the highest. When the inlet duct is short (schemes 1, 2, and 3), the net head is lower than that of scheme 4 due to the poor flow state of the inlet flow passage; when the flow passage is long (schemes 5 and 6), the net head will also be low due to the backflow. Overall, the runner’s length of Case 4 is the most suitable.

Figure 10 shows the inlet duct efficiency with the length of the inlet duct. ηd0 is the inlet duct efficiency of the original scheme. It can be seen from the figure that the efficiency of the inlet duct increases first and then decreases as the length of the flow channel increases. Case 6 has the lowest efficiency of the inlet runner, which is reduced by 3.1%. Although the channel is short (schemes 1, 2, and 3), the loss along the channel is low, but the flow pattern is not good, which consumes a lot of energy, so the channel efficiency is not high; when the channel is long (schemes 5 and 6), the channel efficiency is also low because of the increase of the loss along the channel and the backflow loss. Therefore, Case 4 has the highest efficiency of the inlet duct, which is consistent with the original efficiency.

4.2. Analysis of Flow Uniformity at the Outlet of the Inlet Duct
Figure 11 is a cloud diagram of the velocity distribution of the inlet pipe and the length of the inlet pipe. The results show that there is a low-speed zone above the pump shaft and a high-speed zone below the pump shaft. As the length of the intake pipe increases, the area of the low-speed area above the drive shaft continuously decreases to almost zero, and the area of the high-speed area below the pump shaft also decreases to the middle, so the medium-speed distribution area continues to increase. This shows that as the length of the intake pipe increases, the velocity distribution moves closer to the middle value, the velocity gradient decreases, and the velocity uniformity of the entire water area increases.

Figure 12 is a graph showing the variation of the vorticity at the exit of the inlet duct with the length of the inlet duct. The figure shows that when the length of the inlet duct is short, the vorticity distribution in the area below the pump shaft is symmetrical on both sides. There is a negative vortex in the upper left area of the pump shaft. As the length of the inlet duct increases, the negative vortex area gradually decreases, and the absolute value of the vortex decreases, indicating that the strength of the vortex decreases with the increase of the inlet duct length. The positive vortex is distributed on the upper right side of the pump shaft. The closer to the upper side, the greater the vortex intensity. In particular, the vortex intensity at the upper right of Case 4 increases sharply, but the overall trend is that the vortex value decreases with the increase in the length of the inlet duct. Taken together, the vortex intensity of the outlet section of the inlet duct decreases with the length of the inlet duct.

Table 3 shows the velocity distribution uniformity Vu and the weighted average angle θ of the outlet section of the inlet duct change with the inlet duct length. It can be seen from the table that the uniformity of velocity distribution is the largest when the length of the inlet duct is 1.2L, reaching 65.66%, and the uniformity of velocity distribution is the smallest when the length of the inlet duct is 0.4L, indicating that the proper extension of the inlet duct length can improve the uniformity of the velocity distribution degrees, but the weighted average angle has little effect on the length of the inlet duct.
4.3. Analysis of Flow Characteristics in the Inlet Duct
Figure 13 is a static pressure streamline diagram of the section of the inlet duct with different inlet duct lengths. The figure shows that as the length of the inlet duct increases, the static pressure at the slope surface, inlet elbow, and lip corner of the inlet duct changes significantly. Specifically, the high-pressure area at the corner of the lip gradually decreases, and the static pressure at the slope surface and the inlet elbow does not change much.

The vortex inside the inlet duct occurs in the area above the inlet duct close to the drive shaft, which shows that the drive shaft hinders the flow of water. After Case 4, the vortex area obviously moves to the direction of the water inlet of the flow channel and leaves the vicinity of the drive shaft. This shows that when the length of the inlet duct increases to a certain extent, because the return area gradually increases, it is difficult for the propulsion pump device of the same speed to drive the water flow in the inlet flow passage into the pump body, so the fourth option is the best.
4.4. Thrust Characteristics
Figure 14 shows the change rule of thrust along with the extension of the inclined straight pipe section of the inlet duct. Ft0 is the thrust of the original scheme. With the extension of the inlet duct, the thrust of the entire device first increases and then decreases. The thrust of the fourth solution is the peak value. When the water inlet channel is shortened by 0.2 times or extended by 0.2 times, the thrust is reduced by 2.0%, indicating that a certain length of the inlet channel is beneficial to the propulsion of the ship by the water jet propulsion pump device. The thrust depends on the exit speed. When the flow passage is short (schemes 1, 2, and 3), the inlet velocity of the impeller chamber decreases due to the vortex and other bad flow patterns, so the outlet velocity of the shorter flow passage is smaller under the same speed of the impeller; when the flow passage is long (schemes 5 and 6), the increase of the loss along the way and the occurrence of the backflow phenomenon also leads to the decrease of the inlet velocity of the impeller chamber. In this article, the thrust reaches the highest in Case 4.

5. Conclusion
Through a combination of numerical simulation and experimental verification, the influence of the change in the length of the inlet duct on itself and the entire water jet propulsion system is studied. The research conclusions are as follows:(1)As the length of the inlet duct increases, both the head and the inlet duct efficiency firstly increase and then decrease. When L = 1.12D, they reached the peak value.(2)Extending the length of the inlet duct is beneficial to improve the vortex condition and improve the uniformity of the water flow in the pipeline.(3)The thrust firstly increases and then falls as the length of the inlet duct increases. When L = 1.12D, they reached the highest thrust value.
To sum up, among the six schemes studied in this article, scheme four has the highest comprehensive performance, and the length of the inclined straight pipe section is L = 1.12D.
Data Availability
The curve data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.