Numerical Simulation on Hydraulic Characteristics of Nozzle in Waterjet Propulsion System

: As an important over-current component of the waterjet propulsion system, the main function of a nozzle is to transform the mechanical energy of the propulsion pump into the kinetic energy of the water and eject the water flow to obtain thrust. In this study, the nozzle with different geometry and parameters was simulated based on computational fluid dynamics simulation and experiment. Numerical results show a good agreement with experimental results. The results show that the nozzle with a circular shape outlet shrinks evenly. Under the designed flow rate condition, the velocity uniformity of the circular nozzle is 0.26% and 0.34% higher than that of the elliptical nozzle and the rounded rectangle nozzle, respectively. The pump efficiency of the circular nozzle is 0.31% and 0.14% higher than that of the others. The pressure recovery and hydraulic loss of the circular nozzle are superior. The hydraulic characteristics of the propulsion pump and waterjet propulsion system are optimal when the nozzle area is 30% times the outlet area of the inlet duct. Thus, the shaft power, head, thrust, and system efficiency of the propulsion pump and waterjet propulsion system are maximized. The system efficiency curve decreases rapidly when the outlet area exceeds 30% times the outlet area of the inlet duct. The transition curve forms greatly affect thrust and system efficiency. The transition of the linear contraction shows improved uniformity, and the hydraulic loss is reduced. Furthermore, the hydraulic performance of the nozzle with a linear contraction transition is better than that of others.


Introduction
Waterjet propulsion is a type of special propulsion system that is different from a propeller.This device has been widely used in many high-speed ships because of its high propulsion efficiency, low noise and vibration, and simple transmission structure [1].Waterjet propulsion systems are gradually developed from small and medium to large-sized, highly efficient, low noise, and high-speed mechanisms to meet the requirements of industrial production and military operations [2][3][4][5].A stern-mounted waterjet propulsion system used in commercial applications can be divided into four components, namely, inlet duct, propulsion pump, nozzle, and steering device.Each component of the waterjet propulsion system should function effectively to achieve higher efficiency and reduce energy loss.The inlet duct transfers the bottom water to the propulsion pump.Its performance directly affects the efficiency of the waterjet propulsion systems.The hydraulic loss of the inlet duct also strongly influences the water jet propulsion system.From the perspective of energy loss, Verbeek et al. [6] found that approximately 7-9% of the total power is lost in the inlet duct due to local flow separation and non-uniformity.Jiao et al. [7] used the numerical simulation method to simulate cavitation two-phase flow in the waterjet propulsion pump section and waterjet propulsion system.Ding et al. [8] determined the flow loss of the inlet duct through computational fluid dynamics (CFD).Park et al. [9] observed flow separation and predicted the location of the stagnation point on the lip, in accordance with particle image velocimetry measurement.The flow phenomena occurring within the inlet duct are essential to reduce power loss.The pump is the core component of the waterjet propulsion system, and it delivers the head to produce the jet at the nozzle exit.The main function of the pump is turning the circumrotation power energy into the kinetic energy of water.Waterjet propulsion pumps are of two main types: Mixed-flow pump and axial-flow pump.The efficiency of waterjet propulsion depends on the pump and system efficiency.Insufficient understanding of large deviations between theoretical efficiency and measured efficiency is a major problem in waterjet propulsion design and application.Cao et al. [10] found that the low efficiency was primarily due to the non-uniform inflow of the water-jet pump.Cheng et al. [11] observed the rotating stall region, which is an unstable head curve when many mixed-flow and axial-flow pumps are operated.This region should be generally avoided given the risk of instabilities during pump start and operation.Xia et al. [12] analyzed the rotating stall at a low flow rate and suppressed it with separators.Wang et al. [13,14] simulated the self-priming process of a multistage self-priming centrifugal pump by CFD and optimized the design of a typical multistage centrifugal pump based on energy loss model and CFD.Kim et al. [15] and Etter et al. [16] obtained a number of useful results on the performance optimization of waterjet propulsion pumps through model tests.The thrust of the waterjet propulsion system is obtained by the reaction force of the water flow ejected from the propulsion pump.The main function of the nozzle, which is an important part of the waterjet propulsion system, is to transform the mechanical energy of the propulsion pump into the kinetic energy of the water and eject the water flow to obtain thrust.The gross thrust is used as an alternative because measuring the net thrust of a waterjet installation is cumbersome.The relation between net thrust and gross thrust is not fully understood.Eslamdoost et al. [17] used numerical simulations to investigate this relation.Park et al. [18] conducted a numerical simulation to obtain the complicated viscous flow feature of the waterjet and predict the performance of thrust and torque to obtain waterjet propulsion characteristics.In actual operation, the momentum of the water flow cannot be fully converted into the thrust that propels the ship forward due to the existence of energy loss.The hydraulic loss of the nozzle is closely related to the jet velocity loss because the nozzle is installed immediately after the outlet of the pump.Jian et al. [19] used CFD software Fluent 14.0 for the numerical simulation of four different nozzles and analyzed the effects of geometric and dynamic parameters of nozzles on the momentum thrust of the waterjet propulsion system of autonomous underwater vehicles.Abcand et al. [20], Chin [21], and Jiao et al. [22] performed an overall optimization analysis of the waterjet propulsion system with a nozzle.The steering device can deflect the jet to create steering and reversing forces, and its performance is an important characteristic of the waterjet propulsion [23,24].
In previous decades with the rapid development of computer technology, CFD technology has been widely used in many fields, such as pressure fluctuation and vibration [25][26][27], heat and mass transfer [28,29], flow control [30,31].At the same time, the test method is still one of the effective research methods [32][33][34].In this study, the nozzle with different geometric parameters was simulated based on CFD simulation and experiment.The influences of the nozzles with different geometric parameters on energy loss, efficiency, internal flow characteristics, and hydraulic performance of waterjet propulsion were investigated.The selection of nozzle affects the waterjet propulsion system efficiency.A reasonable nozzle structure is conducive to improving the propulsion system efficiency and reducing energy consumption.

Numerical Model
Calculations are performed using the commercial code ANSYS CFX 14.5.ANSYS CFX uses the element-based finite volume method.The hull boundary layer, inflow velocity, and pressure affect the water flow.Thus, the water flow into the inlet duct is not uniform.The water around the inlet of the waterjet propulsion system should also be included in the computational domain.As shown in Figures 1 and 2, the entire model includes the water body and waterjet propulsion pump system, composed of the inlet duct, propulsion pump, and nozzle.The propulsion pump consists of an impeller with six blades and a guide vane with seven vanes.The rotating speed of the impeller is 700 rpm.D 0 represents the inlet diameter of the impeller.θ represents the dip angle of the inlet duct.
Processes 2019, 7, x FOR PEER REVIEW 3 of 20 Calculations are performed using the commercial code ANSYS CFX 14.5.ANSYS CFX uses the element-based finite volume method.The hull boundary layer, inflow velocity, and pressure affect the water flow.Thus, the water flow into the inlet duct is not uniform.The water around the inlet of the waterjet propulsion system should also be included in the computational domain.As shown in Figures 1 and 2, the entire model includes the water body and waterjet propulsion pump system, composed of the inlet duct, propulsion pump, and nozzle.The propulsion pump consists of an impeller with six blades and a guide vane with seven vanes.The rotating speed of the impeller is 700 rpm.D0 represents the inlet diameter of the impeller.θ represents the dip angle of the inlet duct.

Governing Equations
Reynolds-averaged N-S equation and continuity equation were selected to describe the 3D incompressible viscous flow.The continuity equation is as follows: The momentum conservation equation is as follows: ( ) where ρ represents water density (in m 3 /s), i u , j u represents velocity component of fluid in the i and j directions (in m/s), t represents time (in s), p represents pressure (in Pa), i F represents volume force component in the i direction(in N), μ represents dynamics viscosity coefficient, and i x , j x represents coordinate component.Calculations are performed using the commercial code ANSYS CFX 14.5.ANSYS CFX uses the element-based finite volume method.The hull boundary layer, inflow velocity, and pressure affect the water flow.Thus, the water flow into the inlet duct is not uniform.The water around the inlet of the waterjet propulsion system should also be included in the computational domain.As shown in Figures 1 and 2, the entire model includes the water body and waterjet propulsion pump system, composed of the inlet duct, propulsion pump, and nozzle.The propulsion pump consists of an impeller with six blades and a guide vane with seven vanes.The rotating speed of the impeller is 700 rpm.D0 represents the inlet diameter of the impeller.θ represents the dip angle of the inlet duct.

Governing Equations
Reynolds-averaged N-S equation and continuity equation were selected to describe the 3D incompressible viscous flow.The continuity equation is as follows: The momentum conservation equation is as follows: ( ) where ρ represents water density (in m 3 /s), i u , j u represents velocity component of fluid in the i and j directions (in m/s), t represents time (in s), p represents pressure (in Pa), i F represents volume force component in the i direction(in N), μ represents dynamics viscosity coefficient, and i x , j x represents coordinate component.

Governing Equations
Reynolds-averaged N-S equation and continuity equation were selected to describe the 3D incompressible viscous flow.The continuity equation is as follows: The momentum conservation equation is as follows: where ρ represents water density (in m 3 /s), u i , u j represents velocity component of fluid in the i and j directions (in m/s), t represents time (in s), p represents pressure (in Pa), F i represents volume force component in the i direction(in N), µ represents dynamics viscosity coefficient, and x i , x j represents coordinate component.

Grid Sensitivity Analysis
The computational domain consists of the water body, inlet duct, propulsion pump, and nozzle.The entire computational domain is generated with hexahedral grids and tetrahedral grids by ICEM software.An O-type grid of hexagonal cells was created around the impeller and stator blades by using an extrusion method to ensure good mesh quality in terms of size and skewness.The interior of the impeller and guide vane domain is filled with an unstructured mesh of tetrahedral cells because of the complex topology of the pump.The grid of the calculation domain is shown in Figure 3.
Processes 2019, 7, x FOR PEER REVIEW 4 of 20 The computational domain consists of the water body, inlet duct, propulsion pump, and nozzle.The entire computational domain is generated with hexahedral grids and tetrahedral grids by ICEM software.An O-type grid of hexagonal cells was created around the impeller and stator blades by using an extrusion method to ensure good mesh quality in terms of size and skewness.The interior of the impeller and guide vane domain is filled with an unstructured mesh of tetrahedral cells because of the complex topology of the pump.The grid of the calculation domain is shown in Figure 3.
A grid sensitivity study was conducted to assess the required grid density.Several grid sizes were considered, ranging from a total number of cells of 4.3 × 10 5 up to 2.1 × 10 6 .The standard k-ε is selected as the calculation model in this study.Hence, the y+ values are guaranteed to vary between 30 and 100 to satisfy the computational requirements while changing the grid density.Figure 4 shows the mesh sensitivity analysis of the entire computational domain.The figure reveals that the efficiency and head of the waterjet propulsion pump system increase with the increase in the number of cells.When the number of cells reaches 1.5 million, the head and efficiency remain unchanged with the increase in the number of grids.Theoretically, with the increase in the number and density of grids, the calculation accuracy is generally improved.However, as the number of cells increases, the requirement for computer resources also increases, and the computing speed slows down.Therefore, the final cell number of the entire computational domain is 1,752,299.The y+ value of the propulsion pump is in the range of 30-100 in this study.

Turbulence Model Selecting
The turbulence model was introduced into the numerical calculation to solve the N-S equations.k-ε and k-ω are the commonly used turbulence models in calculation [35][36][37].To select the most A grid sensitivity study was conducted to assess the required grid density.Several grid sizes were considered, ranging from a total number of cells of 4.3 × 10 5 up to 2.1 × 10 6 .The standard k-ε is selected as the calculation model in this study.Hence, the y+ values are guaranteed to vary between 30 and 100 to satisfy the computational requirements while changing the grid density.
Figure 4 shows the mesh sensitivity analysis of the entire computational domain.The figure reveals that the efficiency and head of the waterjet propulsion pump system increase with the increase in the number of cells.When the number of cells reaches 1.5 million, the head and efficiency remain unchanged with the increase in the number of grids.Theoretically, with the increase in the number and density of grids, the calculation accuracy is generally improved.However, as the number of cells increases, the requirement for computer resources also increases, and the computing speed slows down.Therefore, the final cell number of the entire computational domain is 1,752,299.The y+ value of the propulsion pump is in the range of 30-100 in this study.The computational domain consists of the water body, inlet duct, propulsion pump, and nozzle.The entire computational domain is generated with hexahedral grids and tetrahedral grids by ICEM software.An O-type grid of hexagonal cells was created around the impeller and stator blades by using an extrusion method to ensure good mesh quality in terms of size and skewness.The interior of the impeller and guide vane domain is filled with an unstructured mesh of tetrahedral cells because of the complex topology of the pump.The grid of the calculation domain is shown in Figure 3.
A grid sensitivity study was conducted to assess the required grid density.Several grid sizes were considered, ranging from a total number of cells of 4.3 × 10 5 up to 2.1 × 10 6 .The standard k-ε is selected as the calculation model in this study.Hence, the y+ values are guaranteed to vary between 30 and 100 to satisfy the computational requirements while changing the grid density.Figure 4 shows the mesh sensitivity analysis of the entire computational domain.The figure reveals that the efficiency and head of the waterjet propulsion pump system increase with the increase in the number of cells.When the number of cells reaches 1.5 million, the head and efficiency remain unchanged with the increase in the number of grids.Theoretically, with the increase in the number and density of grids, the calculation accuracy is generally improved.However, as the number of cells increases, the requirement for computer resources also increases, and the computing speed slows down.Therefore, the final cell number of the entire computational domain is 1,752,299.The y+ value of the propulsion pump is in the range of 30-100 in this study.

Turbulence Model Selecting
The turbulence model was introduced into the numerical calculation to solve the N-S equations.k-ε and k-ω are the commonly used turbulence models in calculation [35][36][37].To select the most

Turbulence Model Selecting
The turbulence model was introduced into the numerical calculation to solve the N-S equations.k-ε and k-ω are the commonly used turbulence models in calculation [35][36][37].To select the most suitable turbulence model for calculation, this study used several different turbulence models, and the numerical results were compared with the experimental results.Table 1 presents the numerical and experimental results of different turbulence models under the operating conditions of 700 rpm.A slight difference was observed between the calculated values.However, under the condition of the standard k-ε turbulence model, the calculated values are in good agreement with the experimental values.Therefore, the standard k-ε is selected as the calculation model, and the scalable wall-function is used to improve robustness and accuracy.

Boundary Conditions
The inlet of the water body was set as the inlet boundary of the entire computational domain, and normal speed 8 m/s was adopted as the inlet boundary condition.The Reynolds number of the water body inlet is 4.19 × 10 6 .The inflow velocity is equal to ship speed, and nominal turbulence intensities (with a value equal to 5%) are used at the inlet boundary.The outlet of the water body and nozzle were set as the outlet boundary.An average static pressure outlet boundary condition is applied with 1 atm at the water body outlet, and the mass-flow rate was adopted as the outlet boundary of the nozzle.No-slip condition was applied at solid boundaries.The interfaces between the rotational impeller and static diffuser were set as the frozen stage condition.The convergence precision is set to 10 −5 .

Test Rig Set-Up
A test rig is built, as shown in Figure 5, to study the hydraulic performance of the waterjet propulsion pump.The test rig is a closed circulation system that consists of two circulation pipeline systems: Main circulation pipeline and second circulation pipeline.The main circulation pipeline system is designed to ensure the water circulation of the test rig and provide the bottom speed.The second circulation pipeline system is the circulating pipe of the water jet propulsion pump, which is used to test the performance of the waterjet propulsion pump.The propulsion pumping system consists of a propulsion pump, a guide vane, and an inlet duct, as shown in Figure 6.
The propulsion pump with six impeller blades and the guide vane with seven stator vanes are built into the system.The scale model of a waterjet pump is used for ship propulsion.The head, flow rate, torque, and speed were tested to obtain the hydraulic performance of the waterjet propulsion pump.The pump is driven by a DC electromotor at speeds varying from 700 rev/min to 2400 rev/min and equipped with an auxiliary axial pump to regulate the flow rate.The test rig has two electromagnetic flowmeters with an absolute accuracy of ±0.5%.These flowmeters are used to test the flow rate of the main circulation pipeline system and the second circulation pipeline system.The head of the waterjet propulsion pump is measured by the differential pressure transmitter with an absolute accuracy of ±0.2%.

Experimental Verification
As shown in Figure 7, in the numerical calculation, the front section (1-1) of the impeller inlet and the back section of the guide vane outlet (2-2) were considered the pressure measuring sections for the head calculation and analysis.The pump efficiency η, shaft power N, and head H were used to define the hydraulic characteristics of the waterjet propulsion pump system.The calculation formula is as follows: (3)

Experimental Verification
As shown in Figure 7, in the numerical calculation, the front section (1-1) of the impeller inlet and the back section of the guide vane outlet (2-2) were considered the pressure measuring sections for the head calculation and analysis.The pump efficiency η, shaft power N, and head H were used to define the hydraulic characteristics of the waterjet propulsion pump system.The calculation formula is as follows: (3)

Experimental Verification
As shown in Figure 7, in the numerical calculation, the front section (1-1) of the impeller inlet and the back section of the guide vane outlet (2-2) were considered the pressure measuring sections for the head calculation and analysis.

Experimental Verification
As shown in Figure 7, in the numerical calculation, the front section (1-1) of the impeller inlet and the back section of the guide vane outlet (2-2) were considered the pressure measuring sections for the head calculation and analysis.The pump efficiency η, shaft power N, and head H were used to define the hydraulic characteristics of the waterjet propulsion pump system.The calculation formula is as follows: (3) The pump efficiency η, shaft power N, and head H were used to define the hydraulic characteristics of the waterjet propulsion pump system.The calculation formula is as follows:

Flow direction Propulsion pump
where P is the pressure of the section (in Pa), ρ is water density (in kg/m 3 ), g is gravitational acceleration (in m/s 2 ), T is the torque of blades (in N•m), n is the rotating speed of the impeller (in r/min), Q is the flow rate (in m 3 /s), and N is the shaft power (in kW).
The comparative finding of the calculation and the experiment results reveals that the overall trend is similar, as shown in Figure 8.In general, the predicted H and η are in good agreement with the experiment results.Figure 8 shows that the difference between the experiment results and predicted results of H and η is minimum under the design flow-rate condition.In the numerical results, the performance of the head is better than that in the test results.At the design flow rate, the head difference between the test data and numerical data is 1.2%, and the efficiency difference between them is 2.4%.The findings show that the numerical results are reliable.where P is the pressure of the section (in Pa), ρ is water density (in kg/m 3 ), g is gravitational acceleration (in m/s 2 ), T is the torque of blades (in N•m), n is the rotating speed of the impeller (in r/min), Q is the flow rate (in m 3 /s), and N is the shaft power (in kW).
The comparative finding of the calculation and the experiment results reveals that the overall trend is similar, as shown in Figure 8.In general, the predicted H and η are in good agreement with the experiment results.Figure 8 shows that the difference between the experiment results and predicted results of H and η is minimum under the design flow-rate condition.In the numerical results, the performance of the head is better than that in the test results.At the design flow rate, the head difference between the test data and numerical data is 1.2%, and the efficiency difference between them is 2.4%.The findings show that the numerical results are reliable.

The Influence of Nozzle with Different Geometric Parameters on Hydraulic Characteristics of Waterjet Propulsion System
Conducting an in-depth study on the nozzle is necessary to further improve the hydraulic performance and waterjet propulsion efficiency of the waterjet propulsion system.Nine cases are designed to analyze the influence of nozzles with different geometric parameters on the waterjet propulsion system.As shown in Table 2, Cases 1-3 were selected to further investigate the influence of different nozzle outlet shapes on the hydraulic characteristics of the waterjet propulsion system.Case 1 and Cases 4-7 were selected to further investigate the influence of different nozzle outlet areas on the hydraulic characteristics of the waterjet propulsion system.Case 1 and Cases 8-9 were selected to further investigate the influence of different transition curve forms of the nozzle on the hydraulic characteristics of the waterjet propulsion system.

The Influence of Nozzle with Different Geometric Parameters on Hydraulic Characteristics of Waterjet Propulsion System
Conducting an in-depth study on the nozzle is necessary to further improve the hydraulic performance and waterjet propulsion efficiency of the waterjet propulsion system.Nine cases are designed to analyze the influence of nozzles with different geometric parameters on the waterjet propulsion system.As shown in Table 2, Cases 1-3 were selected to further investigate the influence of different nozzle outlet shapes on the hydraulic characteristics of the waterjet propulsion system.Case 1 and Cases 4-7 were selected to further investigate the influence of different nozzle outlet areas on the hydraulic characteristics of the waterjet propulsion system.Case 1 and Cases 8-9 were selected to further investigate the influence of different transition curve forms of the nozzle on the hydraulic characteristics of the waterjet propulsion system.

The Influence of Different Nozzle Outlet Shapes on the Hydraulic Characteristic
The nozzle with a circular shape was mostly used in practical applications, but whether the hydraulic performance of the circular nozzles is optimal remains unclear.Three different nozzle outlet shapes were set up for comparison, as shown in Figure 9, to explore the influence of different nozzle shapes on the hydraulic characteristics of the waterjet propulsion system.The nozzle outlet section shape is changed under the premise of maintaining a constant outlet area of the nozzle.The outlet section of the nozzle in Case 1 is circular, the outlet section of the nozzle in Case 2 is elliptical, and the outlet section of the nozzle in Case 2 is a rounded rectangle.All three outlet sections of different shapes have an area of 30% Ad. a and b denote the long axis and short axis of the elliptical outlet section, respectively.Meanwhile, a and b denote the width and height of the rounded rectangle outlet section.The ratio of a and b is 1.5.The pressure coefficient p C is used to characterize the change of the nozzle outlet pressure.The pressure coefficient p C is a dimensionless value that describes the relative pressure throughout the flow field in fluid dynamics.

Circle shape
Processes 2019, 7, x FOR PEER REVIEW 8 of 20

The Influence of Different Nozzle Outlet Shapes on the Hydraulic Characteristic
The nozzle with a circular shape was mostly used in practical applications, but whether the hydraulic performance of the circular nozzles is optimal remains unclear.Three different nozzle outlet shapes were set up for comparison, as shown in Figure 9, to explore the influence of different nozzle shapes on the hydraulic characteristics of the waterjet propulsion system.The nozzle outlet section shape is changed under the premise of maintaining a constant outlet area of the nozzle.The outlet section of the nozzle in Case 1 is circular, the outlet section of the nozzle in Case 2 is elliptical, and the outlet section of the nozzle in Case 2 is a rounded rectangle.All three outlet sections of different shapes have an area of 30% Ad. a and b denote the long axis and short axis of the elliptical outlet section, respectively.Meanwhile, a and b denote the width and height of the rounded rectangle outlet section.The ratio of a and b is 1.5.The pressure coefficient p C is used to characterize the change of the nozzle outlet pressure.The pressure coefficient p C is a dimensionless value that describes the relative pressure throughout the flow field in fluid dynamics.

Elliptical shape
Processes 2019, 7, x FOR PEER REVIEW 8 of 20

The Influence of Different Nozzle Outlet Shapes on the Hydraulic Characteristic
The nozzle with a circular shape was mostly used in practical applications, but whether the hydraulic performance of the circular nozzles is optimal remains unclear.Three different nozzle outlet shapes were set up for comparison, as shown in Figure 9, to explore the influence of different nozzle shapes on the hydraulic characteristics of the waterjet propulsion system.The nozzle outlet section shape is changed under the premise of maintaining a constant outlet area of the nozzle.The outlet section of the nozzle in Case 1 is circular, the outlet section of the nozzle in Case 2 is elliptical, and the outlet section of the nozzle in Case 2 is a rounded rectangle.All three outlet sections of different shapes have an area of 30% Ad. a and b denote the long axis and short axis of the elliptical outlet section, respectively.Meanwhile, a and b denote the width and height of the rounded rectangle outlet section.The ratio of a and b is 1.5.

The Influence of Different Nozzle Outlet Shapes on the Hydraulic Characteristic
The nozzle with a circular shape was mostly used in practical applications, but whether the hydraulic performance of the circular nozzles is optimal remains unclear.Three different nozzle outlet shapes were set up for comparison, as shown in Figure 9, to explore the influence of different nozzle shapes on the hydraulic characteristics of the waterjet propulsion system.The nozzle outlet section shape is changed under the premise of maintaining a constant outlet area of the nozzle.The outlet section of the nozzle in Case 1 is circular, the outlet section of the nozzle in Case 2 is elliptical, and the outlet section of the nozzle in Case 2 is a rounded rectangle.All three outlet sections of different shapes have an area of 30% A d .a and b denote the long axis and short axis of the elliptical outlet section, respectively.Meanwhile, a and b denote the width and height of the rounded rectangle outlet section.The ratio of a and b is 1.5.

The Influence of Different Nozzle Outlet Shapes on the Hydraulic Characteristic
The nozzle with a circular shape was mostly used in practical applications, but whether the hydraulic performance of the circular nozzles is optimal remains unclear.Three different nozzle outlet shapes were set up for comparison, as shown in Figure 9, to explore the influence of different nozzle shapes on the hydraulic characteristics of the waterjet propulsion system.The nozzle outlet section shape is changed under the premise of maintaining a constant outlet area of the nozzle.The outlet section of the nozzle in Case 1 is circular, the outlet section of the nozzle in Case 2 is elliptical, and the outlet section of the nozzle in Case 2 is a rounded rectangle.All three outlet sections of different shapes have an area of 30% Ad. a and b denote the long axis and short axis of the elliptical outlet section, respectively.Meanwhile, a and b denote the width and height of the rounded rectangle outlet section.The ratio of a and b is 1.5.The pressure coefficient C p is used to characterize the change of the nozzle outlet pressure.The pressure coefficient C p is a dimensionless value that describes the relative pressure throughout the flow field in fluid dynamics.
where p is instantaneous pressure (in Pa), p is average pressure (in Pa), and v out is averaged outlet velocity at the nozzle (in m/s).Figure 10 presents the pressure contours and streamlined diagram of the nozzle outlet section with different nozzle shapes.As shown in Figure 10a, the pressure on the outlet of the circular nozzle is distributed uniformly as a ring.The pressure decreases uniformly from the circumference to the center of the circle.The streamline distribution shows that the flow lines on the outlet section intersect at the center of the circle, and the flow direction of the streamlines is consistent with the rotation direction of the impeller mainly because the guide vanes fail to completely recover the velocity loop.As shown in Figure 10b, the pressure distribution of the elliptical nozzle is relatively uneven.The pressure is symmetrically distributed diagonally.The flow direction of the streamlines is consistent with the direction of the rotation of the impeller.The intersection line of the streamlines at the center of the ellipse is curved.Figure 10c shows four high-pressure zones on the outlet section of the rounded rectangle nozzle.The high-pressure zone is mainly distributed at four rounded corners.The area of the high-pressure zones of the upper and lower side walls is larger than that of the left and right side walls.The flow direction of the streamlines is consistent with the direction of the rotation of the impeller.The intersection line of the streamlines at the center of the rounded rectangle is curved.
Processes 2019, 7, x FOR PEER REVIEW 9 of 21 where p is instantaneous pressure (in Pa), p is average pressure (in Pa), and vout is averaged outlet velocity at the nozzle (in m/s).Figure 10 presents the pressure contours and streamlined diagram of the nozzle outlet section with different nozzle shapes.As shown in Figure 10a, the pressure on the outlet of the circular nozzle is distributed uniformly as a ring.The pressure decreases uniformly from the circumference to the center of the circle.The streamline distribution shows that the flow lines on the outlet section intersect at the center of the circle, and the flow direction of the streamlines is consistent with the rotation direction of the impeller mainly because the guide vanes fail to completely recover the velocity loop.As shown in Figure 10b, the pressure distribution of the elliptical nozzle is relatively uneven.The pressure is symmetrically distributed diagonally.The flow direction of the streamlines is consistent with the direction of the rotation of the impeller.The intersection line of the streamlines at the center of the ellipse is curved.Figure 10c shows four high-pressure zones on the outlet section of the rounded rectangle nozzle.The high-pressure zone is mainly distributed at four rounded corners.The area of the high-pressure zones of the upper and lower side walls is larger than that of the left and right side walls.The flow direction of the streamlines is consistent with the direction of the rotation of the impeller.The intersection line of the streamlines at the center of the rounded rectangle is curved.The uniformity of the velocity distribution at the nozzle outlet is represented by the axial velocity distribution coefficient V .The axial velocity distribution is the best when the axial velocity distribution coefficient is close to 100%.The formula for the axial velocity distribution coefficient V is as follows: where uai is the axial velocity of each element of the calculated section (in m/s), a u is the averaged axial velocity of the calculated section (in m/s), and n is the number of cells of the calculated section.
The velocity uniformity of the three cases is 97.57%, 97.31%, and 97.23%.Different nozzle shapes provide a slight difference in the velocity distribution of the nozzle outlet section.Figure 11 shows the axial velocity distribution curve, and the data point coordinates are normalized.The axial velocity points were uniformly selected for the outlet sections of the three cases along the direction of the vertical and horizontal centerlines.Figure 11a shows the velocity curve along the vertical direction The uniformity of the velocity distribution at the nozzle outlet is represented by the axial velocity distribution coefficient V.The axial velocity distribution is the best when the axial velocity distribution coefficient is close to 100%.The formula for the axial velocity distribution coefficient V is as follows: where u ai is the axial velocity of each element of the calculated section (in m/s), u a is the averaged axial velocity of the calculated section (in m/s), and n is the number of cells of the calculated section.
The velocity uniformity of the three cases is 97.57%, 97.31%, and 97.23%.Different nozzle shapes provide a slight difference in the velocity distribution of the nozzle outlet section.Figure 11 shows the axial velocity distribution curve, and the data point coordinates are normalized.The axial velocity points were uniformly selected for the outlet sections of the three cases along the direction of the vertical and horizontal centerlines.Figure 11a shows the velocity curve along the vertical direction of the center of the outlet section.The vertical direction is the X-axis direction in Figure 10.The high-velocity zones of the three cases are concentrated in the outlet center.The axial velocity drops rapidly close to the upper and lower walls of the nozzle.The three curves show that the axial velocity curve of Case 1 changes more smoothly and evenly.Among the three cases, Case 3 has the greatest range of velocity variation, and the velocity at the center of the nozzle outlet has the largest difference from the velocity at the side wall.Figure 11b shows the velocity curve along the horizontal direction of the center of the outlet section.The horizontal direction is the Y-axis direction in Figure 10.The high-velocity zones of the three cases are concentrated in the outlet center.However, the velocity of the outlet center fluctuates greatly in the horizontal direction.In particular, because the section shrinkage of Cases 2 and 3 is not as uniform as that of Case 1, the velocity changes of Cases 2 and 3 fluctuate intensely.Among the three cases, Case 1 has the greatest range of velocity variation.However, the velocity curve of Case 1 is the most uniform and has a smooth curve.The velocity drops rapidly at the side wall because the velocity near the nozzle wall is very small.The velocity of the three cases in the Y-axis direction is not as uniform as that in the X-axis direction.In particular, the velocity fluctuates greatly near the center in the Y-axis direction.
Processes 2019, 7, x FOR PEER REVIEW 10 of 21 of the center of the outlet section.The vertical direction is the X-axis direction in Figure 10.The highvelocity zones of the three cases are concentrated in the outlet center.The axial velocity drops rapidly close to the upper and lower walls of the nozzle.The three curves show that the axial velocity curve of Case 1 changes more smoothly and evenly.Among the three cases, Case 3 has the greatest range of velocity variation, and the velocity at the center of the nozzle outlet has the largest difference from the velocity at the side wall.Figure 11b shows the velocity curve along the horizontal direction of the center of the outlet section.The horizontal direction is the Y-axis direction in Figure 10.The highvelocity zones of the three cases are concentrated in the outlet center.However, the velocity of the outlet center fluctuates greatly in the horizontal direction.In particular, because the section shrinkage of Cases 2 and 3 is not as uniform as that of Case 1, the velocity changes of Cases 2 and 3 fluctuate intensely.Among the three cases, Case 1 has the greatest range of velocity variation.However, the velocity curve of Case 1 is the most uniform and has a smooth curve.The velocity drops rapidly at the side wall because the velocity near the nozzle wall is very small.The velocity of the three cases in the Y-axis direction is not as uniform as that in the X-axis direction.In particular, the velocity fluctuates greatly near the center in the Y-axis direction.
(a) Vertical direction (b) Horizontal direction The hydraulic performance of the nozzle and guide vane is measured by hydraulic loss.The pressure energy recovery coefficient is introduced to reflect the recovery of static pressure on the nozzle and guide vane.The formulas for the hydraulic loss h  and pressure energy recovery coefficient ξ are as follows: where in P is the total pressure of the inlet section of the nozzle or guide vane (in Pa), and out P is the total pressure of the outlet section of the nozzle or guide vane (in Pa).
Figure 12 shows the relationship between the hydraulic loss and the pressure energy recovery coefficient of the nozzle and guide vane.As shown in Figure 12a, the hydraulic loss of the nozzle is small when the pressure energy recovery coefficient of the nozzle is large.The hydraulic loss and pressure energy recovery coefficients of Cases 1 and 2 are nearly identical when the flow rate is 0.5 to 0.75 times the design flow rate.The hydraulic loss of Case 3 is larger than that of Cases 1 and 2, and the pressure energy recovery coefficient is smaller than that of Cases 1 and 2. The increase in the hydraulic loss of Case 3 is slower when the flow rate is 0.75 to 1 time the design flow rate, initially The hydraulic performance of the nozzle and guide vane is measured by hydraulic loss.The pressure energy recovery coefficient is introduced to reflect the recovery of static pressure on the nozzle and guide vane.The formulas for the hydraulic loss h and pressure energy recovery coefficient ξ are as follows: ξ = P out P in × 100% (10) where P in is the total pressure of the inlet section of the nozzle or guide vane (in Pa), and P out is the total pressure of the outlet section of the nozzle or guide vane (in Pa).
Figure 12 shows the relationship between the hydraulic loss and the pressure energy recovery coefficient of the nozzle and guide vane.As shown in Figure 12a, the hydraulic loss of the nozzle is small when the pressure energy recovery coefficient of the nozzle is large.The hydraulic loss and pressure energy recovery coefficients of Cases 1 and 2 are nearly identical when the flow rate is 0.5 to 0.75 times the design flow rate.The hydraulic loss of Case 3 is larger than that of Cases 1 and 2, and the pressure energy recovery coefficient is smaller than that of Cases 1 and 2. The increase in the hydraulic loss of Case 3 is slower when the flow rate is 0.75 to 1 time the design flow rate, initially larger than that of Cases 1 and 2, and then gradually smaller than that of Cases 1 and 2. The pressure energy recovery coefficient of Case 3 is smaller than that of Cases 1 and 2, and then larger than that of Cases 1 and 2. Under the design flow rate condition, the pressure energy recovery coefficient of Case 3 is the largest, and the hydraulic loss is the smallest.The pressure energy recovery coefficient of Case 2 is minimum.At this time, the pressure recovery coefficients of Cases 1, 2, and 3 are 98.26%, 98.21%, and 98.32%, respectively.When the flow rate is greater than the designed flow rate, the pressure recovery performance of Case 3 is better than that of Cases 1 and 2. Thus, the hydraulic loss of Case 3 is smaller than that of Cases 1 and 2. When the flow rate is greater than 1.25 times the designed flow rate, the hydraulic loss of the three cases decreases more severely with the increase in the flow rate.When the flow rate is 1.5 times the designed flow rate, the pressure energy recovery coefficient of Case 3 is 0.26% and 0.43% higher than that of Cases 1 and 2, respectively.As shown in Figure 12b, the relationship between the pressure energy recovery coefficient and hydraulic loss of the guide vane is opposite to that of the nozzle.In addition, the relationship the pressure energy recovery coefficient (and the hydraulic loss) and the flow rate of the guide vanes is not a monotonic function.The pressure energy recovery coefficient of the guide vane initially increases and then decreases with the increase in the flow rate.The hydraulic loss of the guide vane initially decreases and then increases with the increase in the flow rate.When the flow rate is 1.25 times the designed flow rate, the hydraulic loss of the guide vanes of the three cases reaches the maximum, and the pressure energy recovery coefficient is the smallest.Figure 12a,b shows that when the pressure recovery of the guide vane is worse, the pressure recovery of the nozzle is better.Under the conditions of small flow rate and designed flow rate, the difference between the pressure recovery coefficient and hydraulic loss of nozzles and guide vanes in the three cases is extremely small.Under large flow rate conditions, the pressure energy recovery coefficient and hydraulic loss of Case 3 are optimal, and those of Case 2 are the worst.
Processes 2019, 7, x FOR PEER REVIEW 11 of 21 larger than that of Cases 1 and 2, and then gradually smaller than that of Cases 1 and 2. The pressure energy recovery coefficient of Case 3 is smaller than that of Cases 1 and 2, and then larger than that of Cases 1 and 2. Under the design flow rate condition, the pressure energy recovery coefficient of Case 3 is the largest, and the hydraulic loss is the smallest.The pressure energy recovery coefficient of Case 2 is minimum.At this time, the pressure recovery coefficients of Cases 1, 2, and 3 are 98.26%, 98.21%, and 98.32%, respectively.When the flow rate is greater than the designed flow rate, the pressure recovery performance of Case 3 is better than that of Cases 1 and 2. Thus, the hydraulic loss of Case 3 is smaller than that of Cases 1 and 2. When the flow rate is greater than 1.25 times the designed flow rate, the hydraulic loss of the three cases decreases more severely with the increase in the flow rate.When the flow rate is 1.5 times the designed flow rate, the pressure energy recovery coefficient of Case 3 is 0.26% and 0.43% higher than that of Cases 1 and 2, respectively.As shown in Figure 12b, the relationship between the pressure energy recovery coefficient and hydraulic loss of the guide vane is opposite to that of the nozzle.In addition, the relationship between the pressure energy recovery coefficient (and the hydraulic loss) and the flow rate of the guide vanes is not a monotonic function.The pressure energy recovery coefficient of the guide vane initially increases and then decreases with the increase in the flow rate.The hydraulic loss of the guide vane initially decreases and then increases with the increase in the flow rate.When the flow rate is 1.25 times the designed flow rate, the hydraulic loss of the guide vanes of the three cases reaches the maximum, and the pressure energy recovery coefficient is the smallest.Figure 12a,b shows that when the pressure recovery of the guide vane is worse, the pressure recovery of the nozzle is better.Under the conditions of small flow rate and designed flow rate, the difference between the pressure recovery coefficient and hydraulic loss of nozzles and guide vanes in the three cases is extremely small.Under large flow rate conditions, the pressure energy recovery coefficient and hydraulic loss of Case 3 are optimal, and those of Case 2 are the worst.Figure 13 shows the distributions of the pressure contours of the longitudinal section of the nozzles.Four longitudinal sections are intercepted equidistantly in the nozzle.As illustrated in the figure, seven high-pressure zones due to the influence of guide vane are shown in Plane 4. Compared with the three results, the change in the shape of the nozzle affects the pressure distribution at the nozzle.In Case 1, as the section approaches the outlet of the nozzle, the pressure distribution in the section becomes increasingly uniform.The pressure distribution in Case 3 is the most uneven.In the three cases, the pressure distribution near the outlet is similar to the shape of the outlet section.Figure 13 shows the distributions of the pressure contours of the longitudinal section of the nozzles.Four longitudinal sections are intercepted equidistantly in the nozzle.As illustrated in the figure, seven high-pressure zones due to the influence of guide vane are shown in Plane 4. Compared with the three results, the change in the shape of the nozzle affects the pressure distribution at the nozzle.In Case 1, as the section approaches the outlet of the nozzle, the pressure distribution in the section becomes increasingly uniform.The pressure distribution in Case 3 is the most uneven.In the three cases, the pressure distribution near the outlet is similar to the shape of the outlet section.

The Influence of Different Nozzle Outlet Areas on the Hydraulic Characteristic
Figure 15 shows the 3D diagram of the circular nozzle with different outlet areas.The hydraulic characteristics of the waterjet propulsion pump system are better when the circular nozzle is used.The hydraulic characteristics of different outlet areas of the nozzle were studied under the condition of the circular nozzle.The hydraulic characteristics of the waterjet propulsion systems with five different nozzle areas were studied.The nozzle area varies from 10% to 50% times the outlet area of the inlet duct.

The Influence of Different Nozzle Outlet Areas on the Hydraulic Characteristic
Figure 15 shows the 3D diagram of the circular nozzle with different outlet areas.The hydraulic characteristics of the waterjet propulsion pump system are better when the circular nozzle is used.The hydraulic characteristics of different outlet areas of the nozzle were studied under the condition of the circular nozzle.The hydraulic characteristics of the waterjet propulsion systems with five different nozzle areas were studied.The nozzle area varies from 10% to 50% times the outlet area of the inlet duct.

The Influence of Different Nozzle Outlet Areas on the Hydraulic Characteristic
Figure 15 shows the 3D diagram of the circular nozzle with different outlet areas.The hydraulic characteristics of the waterjet propulsion pump system are better when the circular nozzle is used.The hydraulic characteristics of different outlet areas of the nozzle were studied under the condition of the circular nozzle.The hydraulic characteristics of the waterjet propulsion systems with five different nozzle areas were studied.The nozzle area varies from 10% to 50% times the outlet area of the inlet duct.
Figure 16 shows the pump efficiency, head, shaft power, and IVR (Inlet Velocity Ratio) curves of the waterjet propulsion pump system with different outlet areas of the nozzle.IVR refers to the ratio of the ship speed to the averaged axial outflow velocity at the duct outlet.The averaged axial outflow velocity at the duct outlet is an important parameter to describe the flow phenomena in the inlet, where the speed is changed from the ship speed to the pump velocity.IVR is defined as follows: 11) where v duct is the averaged axial outflow velocity at the duct outlet (in m/s), and v s is the ship speed (in m/s).Figure 16 shows the pump efficiency, head, shaft power, and IVR (Inlet Velocity Ratio) curves of the waterjet propulsion pump system with different outlet areas of the nozzle.IVR refers to the ratio of the ship speed to the averaged axial outflow velocity at the duct outlet.The averaged axial outflow velocity at the duct outlet is an important parameter to describe the flow phenomena in the inlet, where the speed is changed from the ship speed to the pump velocity.IVR is defined as follows: where duct v is the averaged axial outflow velocity at the duct outlet (in m/s), and s v is the ship speed (in m/s).
Figure 16a presents the relationship between efficiency and head with the nozzle area.The head decreases and efficiency initially increases and then decreases with the increase in the nozzle area.The head decreases slowly when the nozzle outlet area changes from 10% Ad to 30% Ad.From the curve, the head drops from 1.14 Hd to 1.0 Hd.Meanwhile, the efficiency curve increases, and the efficiency rises from 0.53 ηd to 1.0 ηd.When the nozzle area exceeds 30% Ad, the head curve drops rapidly, and the efficiency curve begins to decline.When the nozzle outlet area changes from 30% Ad to 50% Ad, the head drops from 1.0 Hd to 0.36 Hd, and the efficiency is reduced from 1.0 ηd to 0.67 ηd.The best efficiency point is obtained when the outlet area is near 30% Ad. Figure 16b shows the relationship between power and IVR with the nozzle area.With the increase in the nozzle outlet area, IVR increases, and the power initially increases and then decreases.When the nozzle outlet area changes from 10% Ad to 30% Ad, IVR rises from 0.32 IVRd to 1.0 IVRd, and the shaft power rises from 0.69 Nd to 1.0 Nd.When the nozzle outlet area changes from 30% Ad to 50% Ad, the IVR curve rises slowly from 1.0 IVRd to 1.37 IVRd.When the nozzle outlet area exceeds 30%Ad, the power curve no longer rises and begins to fall.The maximum power point is obtained when the outlet area is near 30%Ad.Therefore, the maximum efficiency and shaft power occur when the nozzle outlet area is 30% Ad, and the hydraulic characteristics of the waterjet propulsion pump system are better.Figure 16a presents the relationship between efficiency and head with the nozzle area.The head decreases and efficiency initially increases and then decreases with the increase in the nozzle area.The head decreases slowly when the nozzle outlet area changes from 10% A d to 30% A d .From the curve, the head drops from 1.14 H d to 1.0 H d .Meanwhile, the efficiency curve increases, and the efficiency rises from 0.53 η d to 1.0 η d .When the nozzle area exceeds 30% A d , the head curve drops rapidly, and the efficiency curve begins to decline.When the nozzle outlet area changes from 30% A d to 50% A d , the head drops from 1.0 H d to 0.36 H d , and the efficiency is reduced from 1.0 η d to 0.67 η d .The best efficiency point is obtained when the outlet area is near 30% A d .Figure 16b shows the relationship between power and IVR with the nozzle area.With the increase in the nozzle outlet area, IVR increases, and the power initially increases and then decreases.When the nozzle outlet area changes from 10% A d to 30% A d , IVR rises from 0.32 IVR d to 1.0 IVR d , and the shaft power rises from 0.69 N d to 1.0 N d .When the nozzle outlet area changes from 30% A d to 50% A d , the IVR curve rises slowly from 1.0 IVR d to 1.37 IVR d .When the nozzle outlet area exceeds 30%A d , the power curve no longer rises and begins to fall.The maximum power point is obtained when the outlet area is near 30%A d .Therefore, the maximum efficiency and shaft power occur when the nozzle outlet area is 30% A d , and the hydraulic characteristics of the waterjet propulsion pump system are better.Figure 16 shows the pump efficiency, head, shaft power, and IVR (Inlet Velocity Ratio) curves of the waterjet propulsion pump system with different outlet areas of the nozzle.IVR refers to the ratio of the ship speed to the averaged axial outflow velocity at the duct outlet.The averaged axial outflow velocity at the duct outlet is an important parameter to describe the flow phenomena in the inlet, where the speed is changed from the ship speed to the pump velocity.IVR is defined as follows: where duct v is the averaged axial outflow velocity at the duct outlet (in m/s), and s v is the ship speed (in m/s).
Figure 16a presents the relationship between efficiency and head with the nozzle area.The head decreases and efficiency initially increases and then decreases with the increase in the nozzle area.The head decreases slowly when the nozzle outlet area changes from 10% Ad to 30% Ad.From the curve, the head drops from 1.14 Hd to 1.0 Hd.Meanwhile, the efficiency curve increases, and the efficiency rises from 0.53 ηd to 1.0 ηd.When the nozzle area exceeds 30% Ad, the head curve drops rapidly, and the efficiency curve begins to decline.When the nozzle outlet area changes from 30% Ad to 50% Ad, the head drops from 1.0 Hd to 0.36 Hd, and the efficiency is reduced from 1.0 ηd to 0.67 ηd.The best efficiency point is obtained when the outlet area is near 30% Ad. Figure 16b shows the relationship between power and IVR with the nozzle area.With the increase in the nozzle outlet area, IVR increases, and the power initially increases and then decreases.When the nozzle outlet area changes from 10% Ad to 30% Ad, IVR rises from 0.32 IVRd to 1.0 IVRd, and the shaft power rises from 0.69 Nd to 1.0 Nd.When the nozzle outlet area changes from 30% Ad to 50% Ad, the IVR curve rises slowly from 1.0 IVRd to 1.37 IVRd.When the nozzle outlet area exceeds 30%Ad, the power curve no longer rises and begins to fall.The maximum power point is obtained when the outlet area is near 30%Ad.Therefore, the maximum efficiency and shaft power occur when the nozzle outlet area is 30% Ad, and the hydraulic characteristics of the waterjet propulsion pump system are better.The purpose of the waterjet propulsion system is to produce thrust to propel a ship.To study the influence of the nozzle outlet area on the waterjet propulsion system, this study considers two important factors, namely, thrust and system efficiency, in addition to the performance curves.
The thrust of the waterjet propulsion system is defined, considering the influence of the boundary, as follows: where A is the outlet area of the nozzle, v in is the mass-averaged ingested velocity at the duct inlet, and α is the effect coefficient of the boundary layer, which is 0.95 according to the reference documentation.
The efficiency of the waterjet propulsion system η c denotes the ratio of the output power of the system to that of the propulsion pump.Considering the pipeline loss coefficient, the formula for the η c is as follows: where K 1 is the pipeline loss coefficient in the range of 0.40-0.50.In this study, K 1 is 0.45, k' is the ratio of v out to v s vs, and β is the effect coefficient of the boundary layer.
Figure 17 shows the system efficiency and thrust curves of the waterjet propulsion system with different nozzle outlet areas.The thrust and system efficiency of the waterjet propulsion system initially increase and then decrease with the increase in the outlet area.When the nozzle outlet area changes from 10% A d to 30% A d , the thrust rises from 0.27 F d to 1.0 F d , and the system efficiency rises from 0.89 η cd to 1.0 η cd .When the nozzle outlet area exceeds 30% A d , the thrust and system efficiency curves no longer rise and begin to fall.When the nozzle outlet area changes from 30% A d to 50% A d , the thrust reduces from 1.0 F d to 0.51 F d , and the system efficiency declines from 1.0 η cd to 0.19 η cd .The maximum system efficiency and thrust are obtained when the outlet area is close to 30% A d .The nozzle area greatly influences the thrust.When the outlet area exceeds 30% A d , the system efficiency curve decreases rapidly.When the nozzle outlet area is extremely large, the nozzle diffuses the water flow.Thus, the axial velocity of the nozzle outlet v out is less than the velocity of the inlet section of the inlet duct v in , and the thrust value is negative.
shaft power and IVR of the waterjet propulsion pump system under the nozzle outlet area of 30% Ad, respectively.
The purpose of the waterjet propulsion system is to produce thrust to propel a ship.To study the influence of the nozzle outlet area on the waterjet propulsion system, this study considers two important factors, namely, thrust and system efficiency, in addition to the performance curves.
The thrust of the waterjet propulsion system is defined, considering the influence of the boundary, as follows: where A is the outlet area of the nozzle, in v is the mass-averaged ingested velocity at the duct inlet, and α is the effect coefficient of the boundary layer, which is 0.95 according to the reference documentation.
The efficiency of the waterjet propulsion system c η denotes the ratio of the output power of the system to that of the propulsion pump.Considering the pipeline loss coefficient, the formula for the c η is as follows: where K1 is the pipeline loss coefficient in the range of 0.40-0.50.In this study, K1 is 0.45, k' is the ratio of out v to s v vs, and β is the effect coefficient of the boundary layer.
Figure 17 shows the system efficiency and thrust curves of the waterjet propulsion system with different nozzle outlet areas.The thrust and system efficiency of the waterjet propulsion system initially increase and then decrease with the increase in the outlet area.When the nozzle outlet area changes from 10% Ad to 30% Ad, the thrust rises from 0.27 Fd to 1.0 Fd, and the system efficiency rises from 0.89 ηcd to 1.0 ηcd.When the nozzle outlet area exceeds 30% Ad, the thrust and system efficiency curves no longer rise and begin to fall.When the nozzle outlet area changes from 30% Ad to 50% Ad, the thrust reduces from 1.0 Fd to 0.51 Fd, and the system efficiency declines from 1.0 ηcd to 0.19 ηcd.The maximum system efficiency and thrust are obtained when the outlet area is close to 30% Ad.The nozzle area greatly influences the thrust.When the outlet area exceeds 30% Ad, the system efficiency curve decreases rapidly.When the nozzle outlet area is extremely large, the nozzle diffuses the water flow.Thus, the axial velocity of the nozzle outlet vout is less than the velocity of the inlet section of the inlet duct vin, and the thrust value is negative.

The Effect of Different Transition Curve Forms on the Hydraulic Characteristic
Figure 18 shows the transition curve forms of the nozzle.The first component is the nozzle, and the second component is the guide vane outlet.The total length of the nozzle and outlet section of the guide vane L', guide vane outlet diameter D g , nozzle outlet diameter D 1 , and shrinkage arc radius of the guide vane outlet section R 1 are constant.In Cases 1, 8, and 9, the transition curve forms of the nozzle were changed.In Case 1, the linear contraction is selected as the transition curve form.The shrinkage angle θ of the nozzle with linear contraction is 13 • .In Case 8, the combined contraction of the curve and the straight line is selected as the transition curve form.The curve segment consists of arc segments with radii R2 and R3.In Case 9, the arc contraction is selected as the transition curve form.The geometric parameters of the nozzle are shown in Table 3.

The Effect of Different Transition Curve Forms on the Hydraulic Characteristic
Figure 18 shows the transition curve forms of the nozzle.The first component is the nozzle, and the second component is the guide vane outlet.The total length of the nozzle and outlet section of the guide vane L', guide vane outlet diameter Dg, nozzle outlet diameter D1, and shrinkage arc radius of the guide vane outlet section R1 are constant.In Cases 1, 8, and 9, the transition curve forms of the nozzle were changed.In Case 1, the linear contraction is selected as the transition curve form.The shrinkage angle ' θ of the nozzle with linear contraction is 13°.In Case 8, the combined contraction of the curve and the straight line is selected as the transition curve form.The curve segment consists of arc segments with radii R2 and R3.In Case 9, the arc contraction is selected as the transition curve form.The geometric parameters of the nozzle are shown in Table 3.   19 shows the ratio of the hydraulic performance of the three types of transition curves to that of Case 1.The thrust value is the axial thrust from the inlet duct to the nozzle section.The figure indicates that the system efficiency and thrust of the system corresponding to Case 1 are the greatest for the three types of nozzles, but their values are similar to those of Case 8.The main reason is that the hydraulic losses of the two types of nozzles are relatively small.The hydraulic losses of Case 9 are larger than those of the others given the excessive contraction at the nozzle.Figure 19 shows the ratio of the hydraulic performance of the three types of transition curves to that of Case 1.The thrust value is the axial thrust from the inlet duct to the nozzle section.The figure indicates that the system efficiency and thrust of the system corresponding to Case 1 are the greatest for the three types of nozzles, but their values are similar to those of Case 8.The main reason is that the hydraulic losses of the two types of nozzles are relatively small.The hydraulic losses of Case 9 are larger than those of the others given the excessive contraction at the nozzle.

The Effect of Different Transition Curve Forms on the Hydraulic Characteristic
Figure 18 shows the transition curve forms of the nozzle.The first component is the nozzle, and the second component is the guide vane outlet.The total length of the nozzle and outlet section of the guide vane L', guide vane outlet diameter Dg, nozzle outlet diameter D1, and shrinkage arc radius of the guide vane outlet section R1 are constant.In Cases 1, 8, and 9, the transition curve forms of the nozzle were changed.In Case 1, the linear contraction is selected as the transition curve form.The shrinkage angle ' θ of the nozzle with linear contraction is 13°.In Case 8, the combined contraction of the curve and the straight line is selected as the transition curve form.The curve segment consists of arc segments with radii R2 and R3.In Case 9, the arc contraction is selected as the transition curve form.The geometric parameters of the nozzle are shown in Table 3.   19 shows the ratio of the hydraulic performance of the three types of transition curves to that of Case 1.The thrust value is the axial thrust from the inlet duct to the nozzle section.The figure indicates that the system efficiency and thrust of the system corresponding to Case 1 are the greatest for the three types of nozzles, but their values are similar to those of Case 8.The main reason is that the hydraulic losses of the two types of nozzles are relatively small.The hydraulic losses of Case 9 are larger than those of the others given the excessive contraction at the nozzle.The ideal nozzle outflow condition is that the outlet water flow angle is perpendicular to the nozzle outlet section.Thus, the mechanical energy of the propulsion pump can be maximally converted into the kinetic energy of the water.The outlet water flow cannot be completely perpendicular to the nozzle outlet section due to the tangential velocity of the nozzle outlet.Thus, the weighted-velocity average swirl angle θ is used to measure the outflow conditions of the nozzle.The outlet water flow angle perpendicular to the nozzle outlet section is better when θ is closer to 90 • .The formula for velocity-weighted average swirl angle θ is as follows: (17) where u ti is the tangential velocity of each element of the calculated section (in m/s), and u ai is the axial velocity of each element of the calculated section (in m/s).
Table 4 shows the hydraulic performance of the nozzle with different contraction curves.v max denotes the maximum axial outlet velocity of the nozzle, v denotes the averaged outlet axial velocity of the nozzle, V denotes the axial velocity distribution coefficient, and θ denotes the velocity-weighted average swirl angle.
As shown in Table 4, the maximum axial velocity values of the nozzle outlet section of the three cases are basically the same.Therefore, the three cases have good shrinkage effects.The comparative result of the average axial velocity of the three cases shows that the average axial velocity of Case 1 is much larger than that of Cases 2 and 3.The average axial velocity of Case 1 is 1.025 and 1.021 times that of Cases 2 and 3, respectively.The comparative result of the axial velocity distribution coefficient of the three cases shows that the axial velocity distribution coefficient of Case 1 is the best.The axial velocity distribution uniformity of Case 1 is 3.33% and 1.5% more than that of Cases 2 and 3.The velocity-weighted average swirl angle of Case 2 is 80.1 • , which is 1.37 • and 5.03 • higher than that of Cases 1 and 3, respectively.The velocity-weighted average swirl angle of Case 2 is the best mainly due to the rectification effect of the straight section before the outlet of the nozzle in Case 2. Figure 20 shows the axial velocity contours and streamlines of the nozzle outlet section with different contraction curves.The figure shows that the axial velocity of the three cases is distributed in a ring-like manner, and the high-speed zone is concentrated at the center of the circle.The area of the high-speed zone in Cases 2 and 3 is larger, but the axial velocity distribution uniformity of Case 1 is better.The streamline distribution indicates that Case 2 is more affected by the straight section of the outlet of the nozzle, and the flow lines on the outlet section are more deflected.

Conclusions
In this study, the hydraulic characteristics of the nozzle of the waterjet propulsion system were studied by a combination of numerical simulation and experimental verification.Different grid values and turbulence models were used for calculation and analysis to determine the appropriate numerical settings and ensure the accuracy of the calculation.The model test was used to verify the external characteristics of the waterjet propulsion pump, and the numerical results are in good agreement with the test results.Therefore, the numerical results are reliable.Nine different cases were established to systematically study the hydraulic characteristics of the nozzle in accordance with the different nozzle outlet shapes, nozzle outlet areas, and nozzle transition curves.The following results are obtained through calculation and analysis: (1) The pressure coefficient Cp, axial velocity distribution coefficient V , hydraulic loss h  , and pressure energy recovery coefficient ξ were introduced to analyze the hydraulic characteristics of the waterjet propulsion system with different nozzle outlet shapes.The outlet section pressure and streamline the distribution of the circular nozzle case is superior to those in the other two cases.The pressure on the outlet of the circular nozzle is distributed uniformly as a ring, and the flow lines on the outlet section intersect at the center of the circle.Under the design flow-rate condition, the velocity uniformity of the circular nozzle is 0.26% and 0.34% higher than that of the elliptical nozzle and the rounded rectangle nozzle, respectively.The pump efficiency of the circular nozzle is 0.31% and 0.14% higher than that of the others.At this point, the pressure recovery and hydraulic loss of the circular nozzle are superior.Under large flow-rate conditions, the pressure recovery and hydraulic loss of the rounded rectangle nozzle are superior to those in the other two cases.The external characteristic curves show that the circular nozzle case is better.The hydraulic performance of the elliptical nozzle is lower because the elliptical section is more irregular and shrinks more sharply than the other two sections.
(2) The comprehensive hydraulic characteristics are analyzed in accordance with the different nozzle outlet areas, waterjet propulsion pumps, and waterjet propulsion systems.The analysis of the hydraulic performance of the waterjet propulsion pump reveals that when the nozzle outlet area is close to 30% Ad, the waterjet propulsion system has the highest efficiency and the best hydraulic performance.When the nozzle area exceeds 30% Ad, the head curve descends faster and the IVR curve ascends slowly.The analysis of the hydraulic performance of the waterjet propulsion system reveals that when the outlet area is close to 30% Ad, the thrust and system efficiency reach the maximum value.When the outlet area exceeds 30% Ad, the system efficiency curve decreases rapidly.
(3) The comparison among different transition curve forms of the nozzle shows that the transition curve forms greatly affect the thrust and system efficiency.The linear contraction of Case 1 slightly affects the hydraulic performance of the nozzle mainly because the transition of the linear contraction is more uniform with less hydraulic loss.The average axial velocity and axial velocity

Conclusions
In this study, the hydraulic characteristics of the nozzle of the waterjet propulsion system were studied by a combination of numerical simulation and experimental verification.Different grid values and turbulence models were used for calculation and analysis to determine the appropriate numerical settings and ensure the accuracy of the calculation.The model test was used to verify the external characteristics of the waterjet propulsion pump, and the numerical results are in good agreement with the test results.Therefore, the numerical results are reliable.Nine different cases were established to systematically study the hydraulic characteristics of the nozzle in accordance with the different nozzle outlet shapes, nozzle outlet areas, and nozzle transition curves.The following results are obtained through calculation and analysis: (1) The pressure coefficient Cp, axial velocity distribution coefficient V, hydraulic loss h, and pressure energy recovery coefficient ξ were introduced to analyze the hydraulic characteristics of the waterjet propulsion system with different nozzle outlet shapes.The outlet section pressure and streamline the distribution of the circular nozzle case is superior to those in the other two cases.The pressure on the outlet of the circular nozzle is distributed uniformly as a ring, and the flow lines on the outlet section intersect at the center of the circle.Under the design flow-rate condition, the velocity uniformity of the circular nozzle is 0.26% and 0.34% higher than that of the elliptical nozzle and the rounded rectangle nozzle, respectively.The pump efficiency of the circular nozzle is 0.31% and 0.14% higher than that of the others.At this point, the pressure recovery and hydraulic loss of the circular nozzle are superior.Under large flow-rate conditions, the pressure recovery and hydraulic loss of the rounded rectangle nozzle are superior to those in the other two cases.The external characteristic curves show that the circular nozzle case is better.The hydraulic performance of the elliptical nozzle is lower because the elliptical section is more irregular and shrinks more sharply than the other two sections.
(2) The comprehensive hydraulic characteristics are analyzed in accordance with the different nozzle outlet areas, waterjet propulsion pumps, and waterjet propulsion systems.The analysis of the hydraulic performance of the waterjet propulsion pump reveals that when the nozzle outlet area is close to 30% A d , the waterjet propulsion system has the highest efficiency and the best hydraulic performance.When the nozzle area exceeds 30% A d , the head curve descends faster and the IVR curve ascends slowly.The analysis of the hydraulic performance of the waterjet propulsion system reveals that when the outlet area is close to 30% A d , the thrust and system efficiency reach the maximum value.When the outlet area exceeds 30% A d , the system efficiency curve decreases rapidly.
(3) The comparison among different transition curve forms of the nozzle shows that the transition curve forms greatly affect the thrust and system efficiency.The linear contraction of Case 1 slightly affects the hydraulic performance of the nozzle mainly because the transition of the linear contraction is more uniform with less hydraulic loss.The average axial velocity and axial velocity distribution Development of Jiangsu Higher Education Institutions (PAPD), Jiangsu Province Postdoctoral Science Foundation (Grant No. 1701189B), Jiangsu Province Science Foundation for Youths (Grant No. BK20170507).

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 3 .
Figure 3. Grid of the computational domain.

Figure 3 .
Figure 3. Grid of the computational domain.
(a) Surface mesh of impeller and guide vane (b) Inlet duct

Figure 3 .
Figure 3. Grid of the computational domain.

Figure 8 .
Figure 8.Comparison of calculation and experiment results.Qd represents the design flow rate.HBEP and ηd represent the head and efficiency of the waterjet propulsion pump under the designed flow rate condition, respectively.

Figure 8 .
Figure 8.Comparison of calculation and experiment results.Q d represents the design flow rate.H BEP and η d represent the head and efficiency of the waterjet propulsion pump under the designed flow rate condition, respectively.
Ad represents outlet area of the inlet duct.

Figure 9 .
Figure 9. Three-dimensional diagram of nozzles with different shapes.
Ad represents outlet area of the inlet duct.

Figure 9 .
Figure 9. Three-dimensional diagram of nozzles with different shapes.
D1 = 0.55 D0 Ad represents outlet area of the inlet duct.

Figure 9 .
Figure 9. Three-dimensional diagram of nozzles with different shapes.
Processes 2019, 7, x FOR PEER REVIEW 8 of 20 D1 = 0.55 D0 Ad represents outlet area of the inlet duct.

Figure 9 .
Figure 9. Three-dimensional diagram of nozzles with different shapes.

Figure 9 .
Figure 9. Three-dimensional diagram of nozzles with different shapes.

Figure 10 .
Figure 10.Pressure contours and streamline with different nozzle shapes.

Figure 10 .
Figure 10.Pressure contours and streamline with different nozzle shapes.

Figure 12 .
Figure 12.Hydraulic loss and pressure energy recovery coefficient.

Figure 12 .
Figure 12.Hydraulic loss and pressure energy recovery coefficient.

Figure 13 .
Figure 13.Distributions of pressure contours of the longitudinal section of nozzles.

Figure 14
Figure 14 shows the comparison of the hydraulic performance of different nozzle shapes.The figure shows that the trend of the head and efficiency curves of the three cases is consistent.The head initially decreases and then increases with the increase in the flow rate, and then continues to decrease.The head increases with the flow rate when the flow rate is 0.5 to 0.75 times the designed flow rate.The efficiency initially increases and then decreases with the increase in the flow rate.It has the highest efficiency point.The corresponding flow rate of the high-efficiency area is (1.0-1.2) Qd.When the flow reaches (0.42-1.2) Qd, the flow and head curves of the three cases are almost coincident.Each performance value is higher near the design flow range.The head and efficiency are lower when the flow rate is larger than the designed flow rate.The head and efficiency of Cases 2 and 3 are lower than those of Case 1.When the flow rate is larger than 1.25 times the designed flow rate, the head and efficiency curves of Case 2 drop rapidly because the elliptical nozzle is more irregular and shrinks more sharply than those in the other two cases.

Figure 14 .
Figure 14.Hydraulic performance curve of pump system with different nozzle shapes.

Figure 13 .
Figure 13.Distributions of pressure contours of the longitudinal section of nozzles.

Figure 14 21 (Figure 13 .
Figure 14 shows the comparison of the hydraulic performance of different nozzle shapes.The figure shows that the trend of the head and efficiency curves of the three cases is consistent.The head initially decreases and then increases with the increase in the flow rate, and then continues to decrease.The head increases with the flow rate when the flow rate is 0.5 to 0.75 times the designed flow rate.The efficiency initially increases and then decreases with the increase in the flow rate.It has the highest efficiency point.The corresponding flow rate of the high-efficiency area is (1.0-1.2) Q d .When the flow reaches (0.42-1.2) Q d , the flow and head curves of the three cases are almost coincident.Each performance value is higher near the design flow range.The head and efficiency are lower when the flow rate is larger than the designed flow rate.The head and efficiency of Cases 2 and 3 are lower than those of Case 1.When the flow rate is larger than 1.25 times the designed flow rate, the head and efficiency curves of Case 2 drop rapidly because the elliptical nozzle is more irregular and shrinks more sharply than those in the other two cases.

Figure 14
Figure 14 shows the comparison of the hydraulic performance of different nozzle shapes.The figure shows that the trend of the head and efficiency curves of the three cases is consistent.The head initially decreases and then increases with the increase in the flow rate, and then continues to decrease.The head increases with the flow rate when the flow rate is 0.5 to 0.75 times the designed flow rate.The efficiency initially increases and then decreases with the increase in the flow rate.It has the highest efficiency point.The corresponding flow rate of the high-efficiency area is (1.0-1.2) Qd.When the flow reaches (0.42-1.2) Qd, the flow and head curves of the three cases are almost coincident.Each performance value is higher near the design flow range.The head and efficiency are lower when the flow rate is larger than the designed flow rate.The head and efficiency of Cases 2 and 3 are lower than those of Case 1.When the flow rate is larger than 1.25 times the designed flow rate, the head and efficiency curves of Case 2 drop rapidly because the elliptical nozzle is more irregular and shrinks more sharply than those in the other two cases.

Figure 14 .
Figure 14.Hydraulic performance curve of pump system with different nozzle shapes.

Figure 14 .
Figure 14.Hydraulic performance curve of pump system with different nozzle shapes.

Figure 15 .
Figure 15.3D diagram of the circular nozzle with different outlet areas.

Figure 15 .
Figure 15.3D diagram of the circular nozzle with different outlet areas.

Figure 15 .
Figure 15.3D diagram of the circular nozzle with different outlet areas.

Figure 16 .Figure 16 .
Figure 16.Hydraulic performance curve of waterjet propulsion pump system with different outlet areas.*AYrepresents the impeller outlet area.ηd , Hd, Nd, and IVRd represent the pump efficiency, head, Figure 16.Hydraulic performance curve of waterjet propulsion pump system with different outlet areas.* A Y represents the impeller outlet area.η d, H d , N d , and IVR d represent the pump efficiency, head, shaft power and IVR of the waterjet propulsion pump system under the nozzle outlet area of 30% A d , respectively.

Figure 17 .
Figure 17.Hydraulic performance curve of waterjet propulsion system with different outlet areas.

Figure 18 .
Figure 18.Transition curve forms of the nozzles.

Figure 19 .
Figure 19.Hydraulic performance curve of waterjet propulsion pump system with different contraction curves.

Figure 18 .
Figure 18.Transition curve forms of the nozzles.

Figure 17 .
Figure 17.Hydraulic performance curve of waterjet propulsion system with different outlet areas.

Figure 18 .
Figure 18.Transition curve forms of the nozzles.

Figure 19 .
Figure 19.Hydraulic performance curve of waterjet propulsion pump system with different contraction curves.

Figure 19 .
Figure 19.Hydraulic performance curve of waterjet propulsion pump system with different contraction curves.

9 Figure 20 .
Figure 20.Axial velocity contours and streamlines of the nozzle outlet section with different contraction curves.

Figure 20 .
Figure 20.Axial velocity contours and streamlines of the nozzle outlet section with different contraction curves.

Table 1 .
Numerical and experimental results with different turbulent models.td represents the test value of head under design flow rate condition.η td represents the test value of pump efficiency under design flow rate condition. H

Table 3 .
Geometric parameters of the nozzles with different contraction curves.

Table 3 .
Geometric parameters of the nozzles with different contraction curves.

Table 3 .
Geometric parameters of the nozzles with different contraction curves.

Table 4 .
Hydraulic performance of the nozzle with different contraction curves.
R 2 first transition arc radius of nozzle in case 8, m R 3 second transition arc radius of nozzle in case 8, m R 4 transition arc radius of nozzle in case 9, m θ s weighted-velocity average swirl angle, • u ti tangential velocity of each element of the calculated section, m/s v max maximum axial outlet velocity of the nozzle, m/s v averaged outlet axial velocity of the nozzle, m/s