Next Article in Journal
Evaluation of Soil–Foundation–Structure Interaction for Large Diameter Monopile Foundation Focusing on Lateral Cyclic Loading
Next Article in Special Issue
Analysis of Flow and Runner Dynamic Response Characteristics under Pump Conditions of Variable-Speed Pump Turbine
Previous Article in Journal
Research on Scale Effect of Resistance Components for Full-Formed Ship Based on Large-Scale Model Towing Test
Previous Article in Special Issue
Influence of Guide Vane Opening on the Runaway Stability of a Pump-Turbine Used for Hydropower and Ocean Power
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Impeller Optimization Method for the High Specific Speed Mixed-Flow Reactor Coolant Pump Applied to Marine Nuclear Power

1
National Research Center of Pumps, Jiangsu University, Zhenjiang 212013, China
2
School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China
3
Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(7), 1301; https://doi.org/10.3390/jmse11071301
Submission received: 20 May 2023 / Revised: 23 June 2023 / Accepted: 23 June 2023 / Published: 26 June 2023
(This article belongs to the Special Issue Advancement in the Multiphase Flow in Fluid Machinery)

Abstract

:
The reactor coolant pump (RCP) is the only rotating equipment in the primary circuit system of a nuclear power plant and the “heart” of the nuclear reactor. The L formula is defined, and the L/himp is introduced to study the influence of impeller blade type on the performance of the RCP. Twenty groups of models are designed, the concept of arc height ratio is proposed from the perspective of himp and L, and the distribution of internal entropy production within the impeller of the RCP under different Ls and himps of the impeller blade type is analyzed. The results show that when himp remains un-changed and L increases, the low-pressure area at the inlet of the impeller expands while the high-pressure area at the outlet decreases under the design flow or large flow conditions. The smoother blade profile reduces the occurrence of secondary flow phenomena and makes the RCP pressure distribution more uniform. Under design flow and large flow conditions, smaller L/himp and higher himp lead to higher efficiency and head performance. However, higher efficiency and lower head performance can be achieved under small flow conditions with larger L/himp and lower himp.

1. Introduction

The Pressurized Water Reactor (PWR) is the most technologically mature marine nuclear reactor at present, which can be applied to floating power stations, nuclear submarines, nuclear-powered icebreakers, and marine resources development. It is an important measure to solve future offshore energy problems. The Reactor Coolant Pump (RCP) is the “heart” of a nuclear reactor and the only rotating equipment in the primary circuit system of a nuclear power plant [1]. Its primary function is to ensure the circulation of the nuclear primary coolant and remove heat from the core during regular operation. As the primary energy-consuming equipment in nuclear power plants, the RCP’s continuous and stable operation has a critical impact on the safety and stability of the entire plant [2]. Therefore, developing an efficient hydraulic model for the RCP is essential.
In recent years, scholars have conducted extensive analysis of the high-efficiency hydraulic model and internal flow field of pumps. The study of impellers mainly focuses on the in-depth exploration of geometric parameters such as blade thickness distribution, pack angle, and blade inlet position, and finding the optimal values of these parameters using genetic algorithms. Yang Minguan et al. [3] discovered that an impeller with a maximum thickness of about 1/3 of the blade area from the inlet can achieve relatively high hydraulic efficiency, and that an impeller with a rounded blade working face has the best hydraulic performance [4]. Qiu W et al. [5] found that within a certain range, the larger the impeller pack angle, the better the inertia performance of the corresponding pump; the smaller the number of impeller blades, the better the inertia performance of the corresponding pump. Sambhrant et al. [6] found that at the impeller entrance, the inclined blade position was more suitable than the trapezoidal blade, because the paradigm iso effect force distribution is smaller in the inlet inclined blade. Wang Xiuli et al. [7] discovered that the dominant frequency of the impeller remains unchanged with changes in flow or inlet diameter of the short blade. The high frequency of the suction side of the blade gradually decreases, while the high frequency of the pressure face gradually increases with an increase in the inlet diameter of the short blade. Al-Obaidi [8] found that the number of impeller blades has a high impact on pressure, shear stress, magnitude velocity, axial velocity, radial velocity, tangential velocity, and average pressure, and he also [9] studied the behaviors of the flow field and pressure fluctuations in both time and frequency domains in an axial flow pump via the changing of various impeller blade angles. Li H [10] found that different impeller trailing edges can affect the efficiency of the pump; the thinner the trailing edge, the higher the efficiency. Wang H [11] variously tested slot structure geometric parameter combinations to explore this relationship: slot position p, slot width b1, slot deflection angle β, and slot depth h with (3–4) levels were selected for each factor on an L16 orthogonal test table, and the results show that b1 and h are the major factors influencing pump performance under small and rated flow conditions. Xie Rong et al. [12] utilized an optimized Latin hypercube experimental design method to analyze the impact of impeller geometry and its distance from the guide vane inlet on the hydraulic performance of an RCP. Lu Yeming [13,14] employed a united optimal design technology to evaluate the effects of crucial structures (impeller and vane blades) on pump performance. They determined that swirl velocity and installation angle were key design variables for the impeller blades and vane blades, respectively. D X Ye [15] proposed an optimization approach based on the Kriging model and genetic algorithm, utilizing an experimental design with a Latin Hypercube to create 16 design cases containing three main parameters of blade: inlet angle β1, pack angle ϕ, and outlet blade angle β2.
Most previous studies have focused on analyzing individual parameters of the impeller or using algorithms to analyze the optimal values of a few parameters but have yet to consider the coupling relationships between parameters to analyze the influence of impeller geometry on hydraulic performance. This study focuses on a mixed-flow RCP with ns = 461, which utilizes the design concept of axial flow pumps for impeller blade-type design. The L formula is defined, and the L/himp is introduced to study the influence of impeller blade type on the performance of the RCP. Twenty groups of models are designed according to the orthogonal test method, the concept of arc height ratio is proposed from the perspective of himp and L, and the distribution of internal entropy production within the impeller of the RCP under different Ls and himps of the impeller blade type is analyzed. By studying the flow loss of the RCP, the study derives the effect of himp and L on the hydraulic performance of the RCP.

2. Materials and Methods

This study focuses on the mode pump of a mixed-flow RCP with a scaling ratio of 2.866 and a specific speed of 461. The fluid medium density in the pump is 997.561 kg/m3, and the dynamic viscosity of the fluid medium is 8.8871 × 10−4 Pa·s. Table 1 illustrates the scaled model pump parameters. The impeller and guide vane are built using CFturbo software, while a class of spherical snail casing is built using three-dimensional software CREO. Figure 1 depicts the fluid calculation domain of the RCP.
Mixed flow pumps have a specific speed ranging from 300 to 500, while axial flow pumps range from 500 to 1200. The specific speed of the RCP examined in this study is 461, which can be either an axial or mixed flow solution at this parameter. Axial flow pumps operate based mainly on the basic theory of lift forces in fluid mechanics. When the centrifugal impeller runs at high speed, the fluid in the casing is primarily subject to the radial lift force of the impeller to perform work. The design method of traditional mixed flow pumps is basically the same as that of centrifugal pumps. It does not take into account the influence of the radial lift forces on the fluid in the impeller, while the fluid in the impeller of mixed flow pumps is subject to both centrifugal and radial lift forces.
The geometry parameters of axial flow pump impellers mainly consist of chord length, outlet diameter, pitch, blade overlap coefficient, inlet and outlet angles, and hub ratio. This study incorporates the design philosophy of axial flow pumps and axial flow impeller design by treating the blades as made up of infinite arcs. The vanes are unfolded to examine the impact of impeller blades on the hydraulic performance of RCP concerning arc length (chord length and blade overlap coefficient), blade height (outlet diameter and hub ratio), and the interaction between them.
The K1 factor is used to represent the arc length on the horizontal plane, while the K2 coefficient is introduced to correct errors resulting from the angle of repose and to make the arc length on the horizontal plane more accurate. The K2 coefficient accomplishes this by adding an import/export width difference parameter, thereby converting a planar arc length into a spatial arc length, as illustrated in Figure 2. The arc length is calculated using the following formula:
L = 1 2 L H + L S = [ [ K 1 D 1 + ( 1 K 1 ) D 2 ] × ϕ cos [ K 2 β 1 + ( 1 K 2 ) β 2 ] ] 2 + [ 1 2 ( b 1 b 2 ) ] 2 ,
h i m p 1 2 h L + h T = 1 2 b 1 + b 2 ,
where φ is the pack angle of the impeller blade; K1 and K2 are the correction factors; K1 is 1.9274, K2 is 0.5246; D1, D2 are the inlet and outlet diameters of the impeller; β1, and β2 are the inlet and outlet placement angles of the impeller; b1 and b2 are the inlet and outlet widths of impeller; LH is the hub length; LS is the shroud length; himp is the blade height; hL is the leading edge length; hT is the trailing edge length; and L/himp is the ratio of arc length to blade height.
From Formula (1), we can see that the L and himp are mainly determined by φ, D1, D2, β1, and β2. β1 and β2 are mainly determined by the velocity triangle of the blade inlet and outlet and are thus kept constant; b1 is mainly determined by the impeller inlet diameter and has a set value. Therefore, for mixed flow impeller blade types, the L and himp are mainly controlled by φ and b2.
This study utilized an orthogonal test design with 20 different combinations of pack angles ranging from 105° to 125° and outlet widths of 95, 99, 103, and 107. The aim was to investigate how the L and himp affect the performance of RCPs. The design scheme for the orthogonal tests is provided in Table 2.
In this study, entropy production theory is utilized to analyze the distribution of hydraulic losses and the extreme value region of flow losses within the RCP. As the medium being transported is water, it is less affected by temperature changes, so the heat transfer effect can be ignored. Only the specific entropy generation rate Spro,DT caused by the dissipation effect is considered. The internal flow dissipation loss of the RCP is much greater than the wall loss. Therefore, the specific entropy generation rate caused by wall friction is ignored, and the specific entropy generation rate in a turbulent flow Spro,DT is mainly composed of the specific entropy generation rate caused by the time-averaged movement and the fluctuating terms. It can be expressed as
S pro , DT = S p r o , D ¯ + S p r o , D .
The expression for the entropy generation rate caused by the time-averaged movement is
S pro , D ¯ = μ T ¯ 2 u ¯ x 2 + v ¯ y 2 + w ¯ z 2 + u ¯ y + v ¯ x 2 + u ¯ z + w ¯ x 2 + v ¯ z + w ¯ y 2 .
The expression for the entropy generation rate caused by the fluctuation velocity is
S pro , D = μ T ¯ 2 u x 2 + v y 2 + w z 2 + u y + v x 2 + u z + w x 2 + v z + w y 2 .
However, Spro,D is still an unknown term but may be related to the turbulence model [16]. According to Kock’s entropy production theory [17], the entropy generation rate caused by the fluctuation velocity is defined as
S p r o , D = ρ ε T ¯ .
The total entropy production of the calculated domain is the volume integral of the dissipation specific entropy generation rate. SD can be expressed as
S D = V S p r o , D ¯ d V + V S p r o , D d V ,
where μ is the dynamic viscosity; ε is the turbulent dissipation rate; u, v, w are the components of the velocity of the mass in a right-angle coordinate system; and T is the temperature of the fluid mass.
This study utilised the STAR-CCM+ software’s built-in mesh generation module was utilized to generate the mesh. The surface mesh generator model uses the surface remesher and automatic surface repair, whereas the pyramidal mesh generator uses the polyhedral mesh generator. The boundary layer mesh generator selected the prism layer mesher and locally refined the blades of the impeller and guide vane. The mesh of the RCP is displayed in Figure 3. The independence of the fluid calculation domain’s grid for the RCP was examined. Schemes 1 to 5 were established from sparse to dense in proportion of the number of cells for each scheme, as shown in Table 3. When the number of cells reached 1,855,000, the change in RCP efficiency was less than 0.05%. Considering computational resources and the computational accuracy of entropy production, Scheme 3 was selected for numerical simulation.
The fluid medium used in the calculations was water at a temperature of 25 °C. It was assumed to be incompressible with a density of 997.561 kg/m3 and a dynamic viscosity of 8.8871 × 10−4 Pa∙s. The reference pressure was set at 101,325 Pa. The inlet boundary condition was set as a mass flow inlet while the outlet boundary condition was set as a pressure outlet. Interfaces were created for the contact surfaces of each adjacent component and were set as internal interfaces. The realizable k-ε two-layer turbulence model was utilized as the turbulence model with a computational convergence accuracy of 10−4.
Figure 4a displays the test bench photograph, which includes equipment such as a pressure transmitter, force sensor, torque meter, DN350 electromagnetic flow meter, and pressure sensor. The accuracy of the pressure transmitter is ±0.1%, the accuracy of the torque meter is ±0.02%, the accuracy of the DN350 electromagnetic flow meter is ±0.5%. By using the pump unit comprehensive error, the overall performance of the pump unit system can be evaluated, and the efficiency comprehensive error of the pump unit system can be obtained based on the systematic error and random error of the pump unit system. The accuracy error of the pump unit efficiency in this experiment is ±1.48%, which meets the precision requirements of GB33216-89B level and ISO/DIS5198A level, indicating that the experimental results are reliable. Figure 4b shows a comparison between the experiments and simulations. As the simulation and test results exhibit a consistent trend with an error of less than 5%, the simulation data are deemed reliable.

3. Results

The external characteristics of the 20 orthogonal design models are listed in Table 4. It is evident from the table that the efficiency under large flow conditions aligns with the trend of head change, which is primarily influenced by the L and himp. There is little variation under small flow conditions, while under the design flow condition, the efficiency (Ef) changes slightly but there is a significant change in head. The 5th, 6th, 10th, 11th, 15th, and 16th models exhibit an extensive variation range, with all parameters exhibiting an increasing trend except for the efficiency under the design flow condition. The head also rises when the L is kept constant and the himp is increased. On the other hand, when the himp remains unchanged and the L is changed, the head (H) drops. Each group of models maintains a fixed himp and varies the φ design within the range of 105° to 125° in five-degree increments. Comparing the H of the 1st, 6th, 11th, and 16th models, which have different himps but the same L, shows that all have a consistent upward trend, and the efficiency curve also shows a constant upward trend.
A correlation analysis was conducted to investigate the relationship between the geometric parameters of the impeller, including L, himp, L/himp and S, with efficiency and head performance. Figure 5 indicates that S has a weak correlation with efficiency and head performance, while L and himp have a strong correlation with efficiency and head performance. Under both 1.0Qd and 1.2Qd conditions, the impact of each factor on head and efficiency is similar. Specifically, L, L/himp and S are all negatively correlated with head and efficiency, whereas himp is positively correlated with head and efficiency. Moreover, the contribution of each geometric parameter to head and efficiency can be ranked as follows: L/himp > himp > L > S.
Under the 0.8Qd condition, each factor has opposite effects on efficiency and head performance. Specifically, L, L/himp, and S positively affect efficiency, while the head is negatively correlated with them. The contribution of each geometric parameter to efficiency can be ranked as follows: L/himp > L > himp > S; and for head performance, the order is L/himp > himp > L > S. Based on the correlation analysis, it can be concluded that smaller L/himp and higher himp lead to higher efficiency and head performance under design flow and large flow conditions. However, larger L/himp and lower himp under small flow conditions result in higher efficiency and lower head performance.
The head and efficiency indicators of twenty impeller models were comprehensively analyzed under different flow conditions, and four impeller models with different L/himp, L, and himp were selected to analyze these parameters’ effects on the RCP’s work and influx energy conversion mechanism. The selected models were Model 6, 8, 16, and 18, as shown in Figure 6, and named Model I, II, III, and IV, respectively. By comparing these models, analyzing constant himp, increasing L, and larger L/himp, it was found that the efficiency increased at both small flow conditions and design flow conditions, while the efficiency at large flow conditions tended to decrease significantly. When the L remained unchanged, and the himp decreased, the efficiency increased under small flow conditions but decreased under design flow conditions and large flow conditions.

4. Discussion

Figure 7 presents the axial velocity streamline of the RCP under different conditions. A significant number of low-velocity vortex groups can be observed inside the spherical volute, mainly distributed downstream of the volute (A) and the inner part of the diffusion segment (B). Under small flow conditions, the flow in the pump casing is more uniform, with fewer vortex groups and higher flow velocities in the impeller and guide vane flow paths. Models I and III have shorter L and larger velocity gradients within the impeller, resulting in significantly larger high-speed vortex groups that have a significant impact on the inlet side of the guide vane. The fluid outflow after the flow state of the guide vane is worse, and it has a more significant impact on the inner wall of the volute, forming a pronounced secondary flow and vortex groups in the volute.
Model IV impeller exhibits significantly smaller vortex groups at high velocities in the impeller flow path with more uniform velocity gradient changes. As a result, the fluid follows a better streamline after entering the guide vane and volute, leading to smaller vortex groups in the A and B zones. Therefore, under small flow conditions, the impeller is slender (with a large L/himp), and it can perform more adequate work.
Under the 1.0Qd condition, the Model I and II impellers exhibit obvious high-velocity vortex groups, and there is a high-velocity zone and boundary layer separation in the guide vane. However, the flow in the volute is significantly better, with an obvious low-velocity vortex zone in the B zone. The flow in the Model IV impeller is the best, but there is a large vortex in the volute, and there is backflow phenomenon at the guide vane outlet. At 1.2Qd, models I, III, and IV all have good flow patterns within the impeller, but there is reflux at the exit of the guide vane channel, a large number of vortex groups in the volute, and an obvious high-speed jet in the B zone. Model II has an obvious vortex group in the impeller, with high velocity backflow in the guide vane channel, but it performs relatively better inside the volute. At 0.8Qd and 1.0Qd conditions, Model IV has the best flow pattern with uniform velocity gradient changes within the impeller, resulting in the highest efficiency. However, at 1.2Qd conditions, Model IV has a large number of vortex groups inside the volute, leading to the lowest efficiency.
Based on the above analysis, it is evident that the L, himp, and L/himp of the impeller have a significant impact on the flow inside the RCP. A detailed analysis of dissipation losses within the impeller, guide vane, and volute flow channel of the RCP is necessary to assess the impeller’s influence on the flow inside the RCP and to gain a more comprehensive understanding of the work and energy conversion mechanism of the RCP impeller.
Figure 8 and Figure 9 illustrate the entropy production distribution on the impeller’s working surface and suction surface. As depicted in Figure 8, hydraulic losses on the blade’s working surface are primarily concentrated at the blade outlet and outer edge of the blade. As the flow rate increases, the fluid’s impact on the blade intensifies, leading to boundary layer separation and the generation of wake flow.
Consequently, there is an increase in hydraulic losses at the blade outlet and outer edge. Figure 8a shows that at small flow conditions, the entropy production on the working surfaces of different impellers is concentrated near the blade outlet, specifically in the area close to the front cover. Moreover, the entropy production distribution on the suction surface is significantly greater than on the working surface. Model I and III blades have higher entropy production on the outlet edge’s inner (C) side due to fluid decoupling. In contrast, models II and IV exhibit less loss in the C zone and higher loss at the outlet close to the rear cover. As the L increases and the himp decreases, the flow channel narrows, causing the blades to be more constrained to the fluid.
Comparing Figure 8 and Figure 9, the hydraulic losses on the impeller mainly occur on the suction side of the impeller. Upon examination of Figure 9’s entropy distribution diagram for the impeller’s suction surface, it becomes apparent that hydraulic losses at the outer edge of the blade (D) increase as the flow rate increases. Moreover, when the L increases, hydraulic losses in the D zone increase due to hydraulic impact and friction losses. As the L remains unchanged but the himp increases, the hydraulic loss at the inlet of the impeller’s suction side decreases. Model II, with its smaller L and himp, has a greater blade curvature, resulting in a larger hydraulic loss at the inlet. At small flow conditions, secondary flow occurs at the blade inlet, creating a sizable area of hydraulic loss at the inlet due to the narrow flow path of Model II. At large flow conditions, the range of hydraulic losses in the D zone expands due to the rapid impingement and friction of the fluid.
Figure 10 displays the entropy distribution of the axial section of the RCP, providing a visual representation of the internal hydraulic losses in the RCP. The area of the entropy generation rate distribution reflects the efficiency of the RCP, with larger areas indicating higher losses and lower efficiency. Figure 10 shows that hydraulic losses within the RCP are mainly located at the impeller outlet, guide vanes, inner flow channel and E zone.
At 0.8Qd flow condition, the energy losses are mostly concentrated in the guide vane channel. At 1.2Qd flow condition, the energy losses concentrate at the E zone. Under 1.0Qd condition, the overall loss is small, but the loss in the guide vane channel of Model II is larger. Based on the aforementioned velocity flow field analysis, fluid is likely to experience turbulence when flowing along the guide vane channel, which can cause phenomena such as impact. These effects can lead to higher pressure areas in the local guide vane region.
For models I, II, III, and IV, the loss in the E zone increases with increasing L or decreasing himp. Model I has smaller losses in the E zone, but greater losses at the impeller outlet and the inlet of the guide vane. Model IV has the smallest entropy distribution area while Model I has the largest entropy distribution area under small flow conditions.
Due to the increased flow velocity at the impeller outlet of Model I, rotor-stator interaction occurs, resulting in larger losses at both the guide vane outlet and within the guide vane’s flow channel. To address this issue, increasing the arc length can result in more uniform flow in the pump, lower flow velocity at the impeller outlet, and reduced losses. By increasing the L, the flow in the pump becomes more uniform, flow velocity at the impeller outlet decreases and losses are reduced. When the L is unchanged and the himp is reduced, the flow velocity in the impeller channel decreases. This reduces the influence of static and dynamic blade interference, reducing guide vane channel losses. The highest efficiency occurs at an arc height ratio of 2.427, while the lowest efficiency occurs at an arc height ratio of 2.161.
Under large flow conditions, models II and IV experience greater losses due to their smaller impeller himp. When the fluid flows out of the guide vane and into the outlet section, mutual extrusion can occur, resulting in chaotic flow and causing higher losses. Additionally, losses increase at the corner where the volute connects to the exit section as the L/himp increases. Model I, with L/himp of 2.161, is the most efficient, while Model IV, with L/himp of 2.427, is the least efficient.
To further investigate the internal flow losses of the RCP under various flow conditions, the percentage of losses in each component of the RCP and the entropy production of each component was analyzed. Figure 11 and Figure 12 show that the flow loss in the volute is the primary contributor, and as the flow rate increases, the loss in the volute also increases. At 1.2Qd condition, the entropy production in the volute accounts for as much as 76%.
When the flow rate changes from 0.8Qd to 1.0Qd, the proportion of entropy production of Model I and II impellers remains unchanged, while the proportion of entropy production of Model III and IV impellers increases. Although the proportion of Model I and II remains unchanged according to Figure 12, the entropy generation rate of Model I impeller decreases, while that of the Model II impeller increases. This indicates that Model I is suitable for operation under 1.0Qd conditions, while Model II, with its reduced himp and narrower flow channel, is suitable for operation under low flow conditions.
The guide vane entropy production of Model I, III, and IV decreases greatly, indicating that the impeller and guide vane coupling is better. The fluid entering the guide vane can better fit the linear flow of the guide vane. On the other hand, the guide vane entropy production of Model II increases because the short and narrow impeller blade shape leads to a higher velocity of fluid entering the guide vane, causing impact losses to the guide vane.
As the flow rate expands from 1.0Qd to 1.2Qd, the impeller entropy generation rate decreases, and the actual entropy production value also decreases. The entropy production of the guide vane is smaller under large flow conditions. Still, the entropy production of Model II drops significantly due to turbulent flow occurring in the volute which increases losses.

5. Conclusions

The CFD method combined with entropy production theory was used to analyze the effects of variations in L, himp, and L/himp on the hydraulic performance of the model under different conditions of the RCP. Based on the analysis, the following conclusions were drawn:
  • The different impeller types have a significant impact on the hydraulic performance of the RCP. The impeller blade with L/himp of 2.427 exhibits the highest efficiency at 0.8Qd and 1.0Qd conditions, while the impeller blade with L/himp of 2.161 has the highest efficiency at the 1.2Qd condition. Conversely, the impeller blade with L/himp of 2.427 has the lowest efficiency at the 1.2Qd condition. These findings suggest that selecting the appropriate impeller design is critical for optimizing the hydraulic performance of RCPs under different conditions.
  • Increasing the L of the impeller leads to a smoother blade profile and reduces the occurrence of secondary flow phenomena. Additionally, under low flow conditions, increasing the L while decreasing the himp narrows the flow channel, allowing the fluid to better utilize the blade shape to do work. This optimization can lead to improved hydraulic performance and efficiency of the RCP.
  • The dissipation loss of the different models varies greatly as the flow rate increases. When the flow rate changes from 0.8Qd to 1.0Qd, the entropy production of the guide vane decreases significantly for Models I, III, and IV, while the entropy production of the guide vane for Model II increases. Under large flow conditions, the entropy production of the guide vane for each model is lower, and the entropy production of the guide vane for Model II decreases significantly. However, the loss in the volute increases, indicating that Model II is better suited for operation at low flow conditions. These results highlight the importance of selecting the appropriate model based on the conditions to optimize the hydraulic performance of the RCP.

Author Contributions

Conceptualization, Q.F. and Y.L.; methodology, W.Z.; software, Y.Z.; validation, R.Z. and Q.F.; formal analysis, Y.Z.; data curation, W.Z.; writing—original draft preparation, Y.Z.; writing—review and editing, Y.Z.; visualization, Q.F.; supervision, R.Z.; project administration, R.Z.; funding acquisition, Q.F. All authors have read and agreed to the published version of the manuscript.

Funding

Joint Funds of the National Natural Science Foundation of China (U20A20292); National Key Research and Development Program of China (2018YFB0606105); Natural Science Foundation of Jiangsu Province (BK20210771); National Natural Science Foundation of China (51906085); China Postdoctoral Science Foundation Funded Project (2021M701847 and 2022M711768).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhong, Y.; Liu, Y.; Zhou, W.; Xia, D. Investigation on methodology of coastdown analysis and design for reactor coolant pump. Nucl. Power Eng. 2017, 38, 84–88. [Google Scholar]
  2. Cai, L.; Zhang, L. Brief introduction of a pressurized water reactor main pump in a nuclear power plant. Water Pump Technol. 2007, 176, 1–5+9. [Google Scholar]
  3. Yang, M.; Lu, S.; Gao, B. Influence of blade thickness on energy performance of mixed flow nuclear main pump. Fluid Mach. 2015, 43, 28–32. [Google Scholar]
  4. Yang, M.; Zhang, Y.; Zhou, Z.; Wang, Z.; Gao, B. Influence of different blade trailing edges on energy performance of impellers in a mixed-flow nuclear main pump. J. Harbin Eng. Univ. 2017, 38, 230–234. [Google Scholar]
  5. Qiu, W.; Zhu, R.; Wang, X. Influence of optimization of impeller geometric parameters on coast-down performance of nuclear reactor coolant pump. J. Drain. Irrig. Mach. Eng. 2019, 37, 1019–1024. [Google Scholar]
  6. Srivastava, S.; Roy, A.K.; Kumar, K. Design of a mixed flow pump impeller blade and its validation using stress analysis. Procedia Mater. Sci. 2014, 6, 417–424. [Google Scholar] [CrossRef] [Green Version]
  7. Xie, R.; Hao, M.; Jin, W.; Wang, X.; Cheng, Z.; Lu, Z. Hydraulic Model Design and Optimization of Nuclear Reactor Coolant Pump Based on Approximation Model. J. Eng. Thermophys. 2016, 37, 1427–1431. [Google Scholar]
  8. Lu, Y.; Wang, X.; Liu, H.; Li, Y. Investigation of the effects of the impeller blades and vane blades on the CAP1400 nuclear coolant pump’s performances with a united optimal design technology. Prog. Nucl. Energy 2020, 126, 103426. [Google Scholar] [CrossRef]
  9. Lu, Y.M.; Wang, X.F.; Wang, W.; Zhou, F.M. Application of the modified inverse design method in the optimization of the runner blade of a mixed-flow pump. Chin. J. Mech. Eng. 2018, 31, 105. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, X.; Lu, Y.; Zhu, R.; Fu, Q.; Yu, H.; Chen, Y. Study on bidirectional fluid-solid coupling characteristics of reactor coolant pump under steady-state condition. Nucl. Eng. Technol. 2019, 51, 1842–1852. [Google Scholar] [CrossRef]
  11. Ye, D.X.; Lai, X.D.; Zhang, X.; Wang, Q.L. Hydraulic optimization of impeller blade in reactor coolant pump base on Kriging model. IOP Conf. Ser. Earth Environ. Sci. 2019, 240, 032051. [Google Scholar] [CrossRef]
  12. Wang, H.; Long, B.; Wang, C.; Han, C.; Li, L. Effects of the impeller blade with a slot structure on the centrifugal pump performance. Energies 2020, 13, 1628. [Google Scholar] [CrossRef]
  13. Al-Obaidi, A.R. Investigation of the influence of various numbers of impeller blades on internal flow field analysis and the pressure pulsation of an axial pump based on transient flow behavior. Heat Transf. 2020, 49, 2000–2024. [Google Scholar] [CrossRef]
  14. Al-Obaidi, A.R. Analysis of the effect of various impeller blade angles on characteristic of the axial pump with pressure fluctuations based on time-and frequency-domain investigations. Iran. J. Sci. Technol. Trans. Mech. Eng. 2021, 45, 441–459. [Google Scholar] [CrossRef]
  15. Li, H.; Chen, Y.; Yang, Y.; Wang, S.; Bai, L.; Zhou, L. CFD simulation of centrifugal pump with different impeller blade trailing edges. J. Mar. Sci. Eng. 2023, 11, 402. [Google Scholar] [CrossRef]
  16. Zhou, L.; Hang, J.; Bai, L.; Krzemianowski, Z.; El-Emam, M.A.; Yasser, E.; Agarwal, R. Application of entropy production theory for energy losses and other investigation in pumps and turbines: A review. Appl. Energy 2022, 318, 119211. [Google Scholar] [CrossRef]
  17. Kock, F.; Herwig, H. Entropy production calculation for turbulent shear flows and their implementation in CFD codes. Int. J. Heat Fluid Flow 2005, 26, 672–680. [Google Scholar] [CrossRef]
Figure 1. Fluid domain of the RCP.
Figure 1. Fluid domain of the RCP.
Jmse 11 01301 g001
Figure 2. Blade unfolding diagram.
Figure 2. Blade unfolding diagram.
Jmse 11 01301 g002
Figure 3. Grid division of the RCP.
Figure 3. Grid division of the RCP.
Jmse 11 01301 g003
Figure 4. RCP test bench and comparison curve.
Figure 4. RCP test bench and comparison curve.
Jmse 11 01301 g004
Figure 5. Influence of different geometric parameters on the hydraulic characteristics.
Figure 5. Influence of different geometric parameters on the hydraulic characteristics.
Jmse 11 01301 g005
Figure 6. Characteristics of the four model impeller blade shapes.
Figure 6. Characteristics of the four model impeller blade shapes.
Jmse 11 01301 g006
Figure 7. Axial velocity flowline of RCP under different flow conditions.
Figure 7. Axial velocity flowline of RCP under different flow conditions.
Jmse 11 01301 g007
Figure 8. Entropy production distribution on impeller blade working surface.
Figure 8. Entropy production distribution on impeller blade working surface.
Jmse 11 01301 g008aJmse 11 01301 g008b
Figure 9. Entropy production distribution on impeller blade suction surface.
Figure 9. Entropy production distribution on impeller blade suction surface.
Jmse 11 01301 g009
Figure 10. Entropy production distribution of RCP under three flow rates.
Figure 10. Entropy production distribution of RCP under three flow rates.
Jmse 11 01301 g010aJmse 11 01301 g010b
Figure 11. Entropy production ratio of each component of the RCP under the three flow rates.
Figure 11. Entropy production ratio of each component of the RCP under the three flow rates.
Jmse 11 01301 g011
Figure 12. Entropy production of each component of the RCP under the three flow rates.
Figure 12. Entropy production of each component of the RCP under the three flow rates.
Jmse 11 01301 g012
Table 1. Main parameters of the RCP.
Table 1. Main parameters of the RCP.
ParameterValueParameterValue
Rated flow, Qd1080 m3/hImpeller inlet diameter, D1230 mm
Rated head, Hd12 mImpeller outlet diameter, D2278 mm
Rated speed, n1485 r/minNumber of impeller blades, Zi4
Impeller outlet width, b299 mmNumber of guide vane blades, Zd11
Table 2. Impeller designs with different arc lengths and blade heights.
Table 2. Impeller designs with different arc lengths and blade heights.
Sampleφb2LhimpL/himpL× himp(S)
110595244.08107.52.271 26,239
211095254.67107.52.369 27,377
311595265.32107.52.468 28,521
412095276.00107.52.567 29,671
512595286.73107.52.667 30,824
610599244.46109.52.232 26,768
711099255.09109.52.330 27,933
811599265.79109.52.427 29,104
912099276.52109.52.525 30,279
1012599287.29109.52.624 31,459
11105103244.86111.52.196 27,302
12110103255.55111.52.292 28,494
13115103266.28111.52.388 29,691
14120103277.06111.52.485 30,893
15125103287.88111.52.582 32,099
16105107245.29113.52.161 27,840
17110107256.03113.52.256 29,059
18115107266.81113.52.351 30,283
19120107277.63113.52.446 31,511
20125107288.50113.52.542 32,745
Table 3. Grid independence test.
Table 3. Grid independence test.
SchemeCell NumberEfficiency/%
Scheme 198022782.1847
Scheme 2152524982.7779
Scheme 3185513882. 8299
Scheme 4212476382.8722
Scheme 5245640982.9176
Table 4. Table of design characteristics.
Table 4. Table of design characteristics.
0.8Qd1.0Qd1.2Qd
ModelEf/%H/mEf/%H/mEf/%H/m
179.19815.24181.14712.87861.6488.154
281.45115.28681.42912.66260.1027.644
381.82115.10080.83812.23157.0006.939
482.15314.94478.49811.57553.1556.227
583.60814.79577.24411.07649.6885.545
680.52415.76281.65013.30067.2739.184
781.15115.64481.82513.09064.5348.549
881.80815.20282.66012.96456.8926.926
981.82615.23280.75912.28656.6036.955
1082.09714.96878.68411.68053.8376.345
1179.07915.69781.72613.47671.57110.083
1279.28715.65682.17313.45468.5929.381
1380.60215.65481.74913.14467.1138.915
1481.44015.52382.56612.93565.1728.366
1581.77315.25281.24612.42757.6587.109
1678.16915.77082.05213.97374.93410.898
1778.61215.69182.51613.74571.81410.157
1879.80715.66182.63713.52170.3559.660
1980.68915.63883.13213.40769.4439.256
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fu, Q.; Zhao, Y.; Lu, Y.; Zhao, W.; Zhu, R. An Impeller Optimization Method for the High Specific Speed Mixed-Flow Reactor Coolant Pump Applied to Marine Nuclear Power. J. Mar. Sci. Eng. 2023, 11, 1301. https://doi.org/10.3390/jmse11071301

AMA Style

Fu Q, Zhao Y, Lu Y, Zhao W, Zhu R. An Impeller Optimization Method for the High Specific Speed Mixed-Flow Reactor Coolant Pump Applied to Marine Nuclear Power. Journal of Marine Science and Engineering. 2023; 11(7):1301. https://doi.org/10.3390/jmse11071301

Chicago/Turabian Style

Fu, Qiang, Yun Zhao, Yonggang Lu, Weiqiang Zhao, and Rongsheng Zhu. 2023. "An Impeller Optimization Method for the High Specific Speed Mixed-Flow Reactor Coolant Pump Applied to Marine Nuclear Power" Journal of Marine Science and Engineering 11, no. 7: 1301. https://doi.org/10.3390/jmse11071301

APA Style

Fu, Q., Zhao, Y., Lu, Y., Zhao, W., & Zhu, R. (2023). An Impeller Optimization Method for the High Specific Speed Mixed-Flow Reactor Coolant Pump Applied to Marine Nuclear Power. Journal of Marine Science and Engineering, 11(7), 1301. https://doi.org/10.3390/jmse11071301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop