Next Article in Journal
Massive Outbreak of Aurelia coerulea in Geoje Bay, Korea
Next Article in Special Issue
Investigation and Improvement of Centrifugal Slurry Pump Wear Characteristics via CFD-DEM Coupling
Previous Article in Journal
A Practical, Adaptive, and Scalable Real-Time Control Approach for Stormwater Storage Systems
Previous Article in Special Issue
Study on Transient Flow Characteristics of Pump Turbines during No-Load Condition in Turbine Mode Startup
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Influence of Particles and Blade Tip Clearance on the Wear Characteristics of a Submersible Sewage Pump

1
Research Center of Fluid Machinery Engineering and Technology, Jiangsu University, Zhenjiang 212013, China
2
Key Laboratory of Fluid Machinery and Engineering, Xihua University, Chengdu 610039, China
3
Fluid Machinery of Wenling Research Institute, Jiangsu University, Wenling 318000, China
4
Taizhou Longda Submersible Pump Co., Ltd., Taizhou 225500, China
5
Leo Group Co., Ltd., Wenling 317500, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(19), 2845; https://doi.org/10.3390/w16192845
Submission received: 7 September 2024 / Revised: 29 September 2024 / Accepted: 30 September 2024 / Published: 7 October 2024
(This article belongs to the Special Issue Hydrodynamic Science Experiments and Simulations)

Abstract

:
A submersible sewage pump is designed for conveying solid–liquid two-phase media containing sewage, waste, and fiber components, through its small and compact design and its excellent anti-winding and anti-clogging capabilities. In this paper, the computational fluid dynamics–discrete element method (CFD-DEM) coupling model is used to study the influence of different conveying conditions and particle parameters on the wear of the flow components in a submersible sewage pump. At the same time, the energy balance equation is used to explore the influence mechanism of different tip clearance sizes on the internal flow pattern, wear, and energy conversion mechanism of the pump. This study demonstrates that increasing the particle volume fraction decreases the inlet particle velocity and intensifies wear in critical areas. When enlarging the tip clearance thickness from 0.4 mm to 1.0 mm, the leakage vortex formation at the inlet is enhanced, leading to increased wear rates in terms of the blade and volute. Consequently, the total energy loss and turbulent kinetic energy generation increased by 3.57% and 2.25%, respectively, while the local loss coefficient in regard to the impeller channel cross-section increased significantly. The findings in this study offer essential knowledge for enhancing the performance and ensuring the stable operation of pumps under solid–liquid two-phase flow conditions.

1. Introduction

Submersible sewage pumps are extensively utilized in urban sewage management, agricultural irrigation, and industrial water treatment, due to their compact design and superior anti-winding capabilities. These pumps effectively transport sewage media that contain solid particles and fiber materials. Due to velocity disparities between solid particles and the fluid medium, particle collisions with pump components induce wear, adversely affecting the hydraulic performance of the submersible sewage pump, which may lead to shutdown [1].
To date, researchers globally have conducted comprehensive investigations into the complexities of solid–liquid two-phase flow in pumps. For instance, Luo et al. [2] conducted an analysis of the energy conversion characteristics of solid–liquid two-phase flow in an in-line hybrid pump, using an integrated analysis method, based on entropy production theory and relative energy rotor enthalpy. Their findings revealed that the peak values of turbulent kinetic energy (TKE) and the entropy production rate were localized near the tongue, and the flow rate exerted a more significant influence on the TKE than the particle concentration. Tane et al. [3] employed the RANS–DEM method to analyze the motion characteristics of cylindrical particles in a slurry pump, revealing that the average velocity of cylindrical particles decreases with increasing particle density and concentration. Wang et al. [4,5] utilized the renormalization group k-ε model and the dense discrete phase model in Fluent to examine the characteristics of solid–liquid two-phase flow and wear in an open impeller centrifugal pump that handles dense, fine particles. Additionally, the CFD–DEM method was employed to examine the distribution of flow in terms of binary mixed particles within the pump, as well as to investigate particle–fluid interactions at varying volume concentrations. It was found that the linked reaction in terms of the fluid field on the particles decreased as the particle concentration increased, and the time taken to transport small particles decreased with the concentration. Zheng [6] developed a CFD–DEM coupling method, incorporating a fiber model, to analyze how the flow and fiber length affect fiber movement and blockage, revealing that fiber retention locations coincide with vortex sites, notably at the impeller inlet, blade tip, mid-blade, and pump tongue.
Wang et al. [7] studied the space distribution and motion characteristics of particles with varying diameters using the CFD–DEM coupling algorithm, indicating that increasing the particle diameter reduces stratification phenomena, while enhancing the vortex intensity and scale within the guide vane. Cao et al. [8] simulated the flow in pumps featuring diverse particle sizes and shapes by utilizing the CFD–DEM approach in conjunction with the Archard model, finding that larger particles and reduced sphericity resulted in heightened average wear in overflow areas. Wang et al. [9] used a dense discrete phase model to analyze particle–fluid flow and wall erosion in a slurry pump, suggesting that the volute experiences the most erosion and is generally less sensitive to the particle concentration. Guo [10] employed the Rosin–Rammler function to model the sand particle size distribution and investigated key factors influencing component wear in terms of a multiphase pump blade. They identified that the particle Reynolds number primarily affects wear on the pressure side of the blade, while the sand concentration plays a crucial role on the suction side. Peng et al. [11] examined impeller wear in a centrifugal slurry pump, using the Euler–Euler multiphase flow model. It was observed that increasing the particle concentration under identical conditions reduces the pump head and efficiency, while gradually increasing the particle slip velocity at the blade outlet. Cheng et al. [12] employed the Fluent and EDEM coupling approach to study solid–liquid two-phase flow within a pump, observing that high flow rates lead to particle aggregation in the flow channel, forming a long wear strip at the trailing edge of the blade on the pressure surface, with minimal wear on the suction side. Based on the four-way coupled Euler–Lagrangian method, Wang et al. [13] implemented a semi-open centrifugal pump to emulate the behavior of solid–liquid flow with varying particle sizes. They observed that due to effective particle tracking and suppression of the tip leakage vortex, particles readily impacted the trailing edge of the suction surface and the blade hub. As the particle size decreased, more particles flowed into the tip clearance area.
To mitigate friction between the dynamic and static components of the submersible sewage pump, a specific tip clearance is maintained between the blade tip and the pump shell in open impeller submersible sewage pumps. However, the tip clearance induces leakage flow during pump operation, resulting in a tip leakage vortex that impacts the hydraulic performance of the pump. Therefore, investigating the effect of the tip clearance size on the internal flow dynamics of the pump is crucial for optimizing its performance. Yang et al. [14] studied the combined effects of diverse two-phase flow patterns and different tip gap sizes on pump performance by using the multi-factor variance analysis method. The findings indicate that pump performance significantly deteriorates as the clearance size increases, and the sensitivity of the pump efficiency to the change in the tip clearance size depends on the flow rate. Ji et al. [15] carried out energy performance and pressure tests to investigate how tip clearance affects pressure pulsations in mixed-flow pumps, identifying low-frequency pulsations with concentrated energy at the impeller inlet as the primary interference frequency. Increasing the tip clearance expanded the high-amplitude region of the wavelet spectrum towards lower frequencies. Lin et al. [16] explored the fluctuation characteristics of the tip clearance leakage flow and pressure pulsation in response to different flow conditions. The observations indicate that with higher flow rates, the amplitude of the pressure fluctuations at the tip clearance increases initially and then declines, while the effect of the tip clearance size on the tip clearance flow rate becomes less significant.
Yang et al. [17] discussed the effects of tip clearance design on both the energy characteristics and flow dynamics, through a combination of numerical calculations and experimental verification, and found that the complicated mixing action of leakage and the mainstream flow is the main source of energy dissipation in a semi-open impeller. Guo [18] investigated the influence of four different tip clearance thicknesses in inducers on the high-speed performance of centrifugal pumps under various flow conditions. The findings indicate that increased tip clearance, particularly under high-flow conditions, enhances both the efficiency and pump head. Peng et al. [19] examined the impact of the tip gap between the impeller and blade tip on the hydraulic performance of a single-stage centrifugal pump, through experiments and numerical simulations. They observed that an increase in the tip clearance significantly reduced the impact losses at the tongue and gradually mitigated interference from water exiting the tip clearance at the impeller inlet. Li [20] investigated the impact of leakage flow through the tip clearance on the rotating stall and observed that the critical stall point shifted to higher flow rates under a smaller tip clearance, consequently reducing the stable operation scope of mixed-flow pumps. Han [21] discovered that increasing the tip clearance in submersible electric pumps decreases the energy performance. Specifically, at lower flow rates, increasing the tip clearance proportionally decreases the energy performance. Under higher flow conditions, a greater tip clearance results in significant leakage increases.
In conclusion, extensive research has been undertaken on solid–liquid two-phase flow and tip clearance within this field, yet the dynamics within pumps and the mechanisms governing phase interactions are still not fully understood. Furthermore, the interaction between particles and the tip clearance size under multiphase solid–liquid conditions, and their combined effects on the pump’s energy characteristics, have not been thoroughly investigated. Therefore, there is an urgent need for further exploration of the flow characteristics in terms of solid–liquid two-phase systems and their interaction mechanisms with tip clearance. In this study, an open impeller submersible sewage pump is taken as the research object and the CFD–DEM coupling model is deployed to explore the wear properties of the internal flow components under different transportation flow rates and particle parameters. Additionally, three different tip clearance sizes are considered to comprehensively analyze the tip leakage vortex structure, internal flow field patterns, and wear characteristics, across different configurations. Utilizing the energy balance equation, the relationship between the tip clearance and the energy conversion in the pump is deeply explored. The research in this paper can not only enrich the theory on solid–liquid two-phase flow within pump systems, but also serves as a reference for designing and optimizing the tip clearance of open impeller sewage pumps, possessing significant theoretical and practical engineering value.

2. Geometric Model and Parameters

This study focuses on an open impeller submersible sewage pump, characterized by a low specific speed, with a specific speed of ns = 24.8, a design flow rate of Qd = 10 m3/h, a design head of Hd = 28 m, and a rated speed of n = 2980 r/min. Figure 1 illustrates the two-dimensional and three-dimensional configuration of the pump. To mitigate the effects of inlet and outlet reflux, the inlet and outlet pipes are lengthened to five times their diameter. The key geometric parameters of the model pump are listed in Table 1.

3. Numerical Modeling

3.1. Particle Model

To ascertain the impact of the solid–liquid medium properties on the internal wear and flow dynamics in the submersible sewage pump, various solid particle configurations were employed within the solid–liquid two-phase flow, as detailed in Table 2. Three variables, namely the flow conditions, particle volume fraction, and particle shape, were considered in the study. Figure 2 illustrates the three particle shapes: spherical particles with a diameter of 1 mm, two spherical particles piled up to form a short cylinder (length 2.5 mm), and four spherical particles piled up to form a long cylindrical particle (length 5 mm).

3.2. Numerical Simulation Method

To address potential issues arising from high local particle volume fractions during two-phase flow simulations, the following fluid control equations are implemented:
t α ρ f + α ρ f v f = 0
ρ f t ( α v f ) + α ρ j v f v f = ( α P ) + ( α μ ) v f + ( v f ) T 2 3 v f l + ρ f g + F
where ρf stands for fluid density; νf represents the fluid velocity vector; t is the time variable; P indicates pressure; g denotes gravitational acceleration; α is the volume fraction of the liquid phase; and F represents other additional forces (including the interaction force between discrete phases, etc.).
During the simulation of the fluid flow within the submersible sewage pump, the finite volume method and SST k-ω turbulence model are applied, and the rotation effect of the impeller is realized by applying the sliding grid method. In addition, the particles are discretized, according to the Euler–Lagrange method, to calculate the movement of the particle phase. In this study, the particle volume fraction is below 10%, which means that it belongs to a sparse solid–liquid two-phase flow. The solid particles are described by the following governing equation:
m s d v s d t = F drag + F g + F vm + F p + F saff + F Magn + F c
I s d ω s d t = T c + T f
I s = 0.1 m s d s 2
where Fdrag represents the drag force on the particles; Fg denotes the resultant force of gravity and buoyancy acting on the particles; Fvm signifies the additional mass force on the particles; Fp denotes the pressure gradient force on the particles; Fsaff accounts for the force generated by the difference in the vertical flow velocity between the particles; FMagn represents the force on the particles due to rotation; Fc denotes the contact force between the particles and between the particles and the wall; Is represents the moment of inertia of the particle; ωs represents the angular velocity of the particle; Tc and Tf [22], respectively, represent the contact torque of the particle and the torque resulting from the fluid phase; and ms indicates the mass of the particle, whereas ds refers to its size.
Particles are primarily influenced by gravity and drag force during movement [23]. The Di Felice drag model, which considers the fluid volume fraction and porosity, is utilized to compute the drag force. The formula is given below:
F drag = 1 8 π d s 2 ρ f C D v f + v f v s v f + v f v s α f 1 α
α = 3.7 0.65 e x p 1.5 log 10 Re s , α 2 2
where Res,α is the particle Reynolds number under the action of the fluid volume fraction, which can be expressed as:
R e s , α = ρ f d s α v f + v f v s μ f
The combined force of gravity and buoyancy Fg can be expressed as follows:
F g = m s ρ s ρ f ρ s g
where ρs represents the particle density.
The pressure gradient force arises from the pressure differential induced by the accelerated fluid motion over the surface of the particles, resulting from the particle movement within the flow field with changing pressures. It can be expressed as follows:
F p = V s p
The additional mass force is the force required to push the fluid to accelerate when the particle velocity is greater than the fluid velocity, which can be expressed as:
F vm = ρ f 2 ρ s v s v f + v f d v s d t
when the ratio of the fluid density to the particle density exceeds 0.1, the effects of the virtual mass force and pressure gradient force on the particle motion become significant [24]. In this paper, the fluid density ρf = 998.2 kg/m3 and particle density ρs = 2650 kg/m3 ratio is greater than 0.1, thus necessitating the virtual mass force and pressure gradient force in regard to the coupling interface.
The Fsaff is the vertical force on the particle due to the different fluid velocities across the upper and lower surfaces [25], which is expressed as:
F saff = 1.615 d s 2 ( μ f ρ f ) 1 / 2 C saff | w c | 1 / 2 ( v f v s ) w c
where wc is the vorticity of the fluid; Csaff denotes the slip-shear lift coefficient, which correlates closely with the particle Reynolds number, and can be expressed as:
C saff   = ( 1 0.3314 β ) e x p R e s 10 + 0.3314 β ,   R e s 40 0.0524 β R e s ,   R e s > 40
In the formula, the shear Reynolds number is expressed as Res = ρfds2|wc|/μf and the coefficient is β = Res/2 Res,α.
The FMagn is the force applied to the particle due to rotation, calculated according to the Kutta–Joukowski theorem [26]. The expression is:
F Magn = 1 2 A p C Magn ρ f | V | | Ω | ( V × Ω )
where Ap is the projection area of the particles; Ω signifies the angular velocity of the fluid in relation to the particles; V stands for the linear velocity of the particle as compared to the fluid; and CMagn is the rotational lift coefficient, predominantly impacted by the rotational Reynolds number Res,r and the particle Reynolds number Res,α. It is given by:
C Magn = 0.45 + ( R e s , r R e s , α 0.45 ) e ( 0.057 π ) R e s , r 0.4 R e s , α 0.3 )
R e s , r = ρ f d s 2 ω c μ f
Fc denotes the force arising from mutual contact between the particles and can be formulated as:
F c = F c , n + F c , t
where Fc,n refers to the normal force exerted on the particle after collision and Fc,t represents the tangential force on the particle after collision, which can be calculated using the following expression:
F c , n = 1.5 Y 1 μ 2 R 1 , 2 u n 3
F c , t = 12 ( G 1.2 R 1.2 u n ) u q 3
where μ signifies Poisson’s ratio, Y defines the equivalent elastic modulus, R1,2 represents the equivalent contact radius, un represents the normal displacement, G1,2 represents the equivalent shear modulus, and uq represents the tangential displacement.
The Archard wear model [27] is a wear model grounded in contact mechanics between the particles and the walls, used to characterize solid material interactions on frictional surfaces. It assumes that wear occurs due to the relative movement between small particles or surface protrusions, which break or detach upon contact, resulting in material loss. The semi-empirical formula is given by:
Q v = K H v F n d s
where Qv represents the amount of wear, where the unit is mm3; Hv represents the hardness of the material surface, where the unit is N∙mm−2; K denotes the wear constant; Fn denotes the point load, in N; and ds represents the sliding distance, where the unit is mm.
The wear amount is represented by the wear depth h, which can be expressed as:
h = Q v A
where A indicates the contact region between the wall and the particle, where the unit is mm2.

3.3. Boundary Condition

The Fluent–EDEM coupling calculation simulates solid–liquid two-phase flow in pumps by treating the liquid phase continuously and the solid phase discretely. The Fluent 2022 R1 software computes the liquid flow field, while EDEM simulates particle motion. Figure 3 depicts the coupling mechanism used to transfer momentum between the solid and liquid phases. As demonstrated in our previous research [28], the Fluent–EDEM coupling calculation can effectively simulate the solid–liquid two-phase flow in the pump.
The Fluent simulation employs the SST k-ω turbulence model, with an inlet velocity of 1.0514 m/s, under the design conditions, while the outlet boundary condition is configured as free outflow (outflow). In the EDEM, the pump’s inlet surface is configured as a particle factory, generating a certain number of particles per second based on the computational parameters. The contact model between the particles and between the particle walls uses the Hertz–Mindlin (no-slip) model, and the wear rate of the wall in terms of the flow-through parts is calculated using the Archard wear model. The time step for the calculations is set to 1 × 10−6 s in the EDEM and 1 × 10−4 s in Fluent. Both simulations use a residual convergence criterion of 1 × 10−6. Table 3 presents the key parameters for the EDEM and the Fluent simulation.

3.4. Grid Division and Independence Analysis

The grid quantity and quality substantially affect both the speed and accuracy of numerical simulations of submersible sewage pumps. To optimize the computational resources, Fluent meshing was employed to mesh the components, including the inlet pipe, tip clearance, impeller, rear chamber, volute, and outlet pipe of a submersible sewage pump, featuring a 0.7 mm tip clearance, in clear water conditions. Five sets of polyhedral grids were examined for grid independence, all exhibiting a grid quality above 0.3. Figure 4 shows the details of the polyhedral grid generation.
The variations in the pump head and efficiency across different grid configurations were analyzed to determine the optimal mesh density for a submersible sewage pump. Based on the findings shown in Table 4 and Figure 5, the fluctuations in the efficiency and pump head remain within 1% when the grid count exceeds 652,543. Consequently, grid scheme group 3, comprising 652,543 grids, was selected as the most suitable configuration.

4. Analysis of the Influence of Particle Parameters on Pump Performance

4.1. The Influence of the Flow Rate on Flow Patterns and Wear Characteristics

The wear characteristics of the impeller and the volute of a submersible sewage pump, when subjected to long cylindrical particles at a volume fraction of 1%, are comprehensively investigated under varying flow conditions. Figure 6 illustrates the velocity vector distribution within the pump across different flow rates. It has been discerned that at lower flow rates in terms of solid–liquid two-phase flow, a significant low-speed region emerges on the suction side of the impeller near the tongue, leading to flow separation and the formation of small vortices that can potentially cause channel blockage. Conversely, the pressure side exhibits relatively stable flow conditions, with no distinct low-speed zone. Under conditions in terms of the designed flow rate and high flow rates, the fluid velocity within the impeller channels increases, causing a gradual weakening and migration of flow separation-induced vortices towards the blade’s pressure surface. The highest velocity area appears near the second section of the volute, and the peak decreases with increasing flow rates, while the velocity at the outlet of the volute gradually decreases with the increase in the flow rate.
Figure 7 depicts the TKE distribution within the pump across various solid–liquid two-phase flow conditions. The maximum turbulent kinetic energy is observed at the volute tongue and the middle section of the impeller channel. Under the design conditions, the high turbulent kinetic energy area is significantly reduced at lower and higher flow rates. The impeller channel exhibits a more uniform distribution in terms of the TKE, indicating more regular fluid movement. Near the volute outlet, the turbulent kinetic energy gradually increases with the flow rate, leading to a more disordered distribution. This phenomenon results from higher flow rates increasing the fluid velocity at the volute outlet, intensifying the turbulent kinetic energy in that region.
Figure 8 demonstrates the distribution pattern of the particles within the pump under varying flow rates. It can be observed that the inlet of the impeller is in a low-speed zone and the velocity of the particles is low. In this region, the impact of the particles on the leading edge of the blade, along with their subsequent rebound and further collisions, will be of particular significance. With the increase in the flow rate, the particle velocities increase, augmenting their travel distance per unit time and, thereby, diminishing particle aggregation towards the impeller inlet, resulting in enhanced dispersion. Consequently, the time required for particle concentration within the pump to stabilize is reduced. In addition, with the particle volume fraction held as constant, an increase in the flow rate leads to an increase in both the number of particles within the pump and the transport load on the impeller. The quantity of particles conveyed through each flow channel of the impeller is significantly increased, resulting in serious wear to the flow components in this area. Concurrently, particle aggregation intensifies along the volute wall surface with increasing flow rates, predicting a gradual increase in the wear rate. At the outlet pipe of the submersible sewage pump, significant particle aggregation is only observed at low flow rates. This phenomenon is primarily attributed to the inadequate fluid velocity within the outlet pipe at low flow rates, resulting in insufficient transportation capacity for solid particles.
Figure 9 illustrates the distribution of the wear rates for the blade and volute in the submersible sewage pump for various flow rates. From the diagram, it can be seen that the area with the largest wear rate on the blade is mainly distributed on the leading edge and tail of the blade. With a constant solid volume fraction, the leading-edge area exhibits progressively higher wear rates as the flow rate increases. This trend is primarily attributed to elevated particle velocities and increased particle counts associated with higher flow rates, leading to intensified wear. The rate of wear in the volute grows with the flow rate from section I to section V, while it decreases initially and then increases again from section V to section VIII. This is due to the small flow rate and the insufficient capacity of the submersible sewage pump to transport particles, resulting in the deposition of particles from section V to section VIII. The wear rate of the volute decreases as the flow rate increases, owing to the enhanced particle conveyance capacity of the pump. However, the excessive flow rate increases the particle velocity and quantity, accelerating volute wear and, consequently, increasing the wear rate.

4.2. The Influence of the Particle Volume Fraction on the Pump Wear Characteristics

To investigate the characteristics of the particle transportation by a submersible sewage pump for various particle volume fractions and particle distributions, the wear analyses were conducted at 1%, 2%, and 3% particle volume fractions. The flow rate was maintained in terms of the design condition and the particles were long cylinders in terms of their shape. Figure 10 illustrates the particle distribution within the submersible sewage pump at varying particle volume fractions. The number of particles transported per unit of time exhibited a linear increase with an increasing particle volume fraction. Collisions both among the particles and between the particles and the leading edge of the blade were considerably intensified, leading to a reduced particle velocity at the inlet and increased momentum loss. Therefore, a substantial number of particles accumulated at the impeller inlet, causing a blockage effect. Particles entering the impeller passage experienced significant collisions, predominantly at the leading edge of the blade’s pressure side and the trailing edge of the suction surface. The amount of particles transported through each flow channel of the impeller increases substantially with the increasing concentration. The primary accumulation zone for the particles extends from section V to section VIII, along the volute wall. Higher concentrations result in more pronounced particle aggregation within this region.
The wear rates of the blades and volute in the submersible sewage pump are shown in Figure 11, as influenced by varying particle volume fractions. The leading edge located on the pressure surface, the trailing edge on the suction surface, and the V–VIII region of the volute, are the main wear areas. The leading edge of the blade is where the initial particle collisions occur, thus causing heightened wear at the front edge of the blade’s pressure side and the trailing edge on the suction surface. An increase in particle volume fraction expands the wear area and significantly increases the wear depth, necessitating timely replacement of the impeller due to continuous operation under high volume concentrations. Figure 11b illustrates that as the concentration increases, the volute’s high-wear area experiences deeper wear, leading to substantial material removal from the wall and, thereby, reducing the volute’s lifespan considerably. Moreover, the tongue barrier intensifies particle collisions and increases the collision frequency. However, the volute exhibits a lower wear rate. This is mainly due to particles colliding with the tongue, which prevents direct impact on the volute wall and, thereby, reduces the wear on the volute tongue to some extent.

4.3. The Influence of Particle Shape on the Wear Characteristics

In practical solid–liquid fluid transportation applications, submersible sewage pumps are frequently employed to convey particles of varying shapes. This study simulates operational scenarios involving a submersible sewage pump, with a designed flow particle volume fraction of 1%, to investigate how different particle shapes affect the wear characteristics of the flow components. Additionally, the volume concentration of the particles generated per unit of time for different particle shapes is the same.
Figure 12 illustrates the particle distribution in the submersible sewage pump for various particle shapes. It reveals that at the same particle volume fraction, the pump transports different quantities of particles depending on their shapes. The order of the shapes, from the highest to the lowest number of particles, is spherical, mixed, short cylindrical, and long cylindrical. Upon observing the particle distribution, it was noted that the velocity distribution among the particles of different shapes in the submersible sewage pump is nearly uniform. However, the uniformity of the particle distribution within the impeller channel correlates notably with the particle shape. Spherical particles exhibit the most even distribution in the impeller channel, followed by mixed particles and then short cylindrical particles. Long cylindrical particles display the most chaotic distribution, consistent with their lower distribution value among the different particle shapes. The reason is that, on the one hand, the contact area with the fluid is expanded due to the increase in the surface area of the particles, which intensifies the interaction with the fluid and elevates the resistance of the particles. On the other hand, the equivalent diameter of spherical particles is smaller. Compared with long and short cylindrical particles and mixed-shape particles, the Stokes number (St number) of spherical particles is lower, and the drag force is dominant. Consequently, spherical particles demonstrate superior flowability and more uniform distribution and movement within the fluid.
The wear rate distribution in submersible sewage pumps for different particle shapes is shown in Figure 13. As evidenced by Figure 13a, the high wear rate area distribution in the submersible sewage pump when transporting the above four shapes of particles is consistent, primarily distributed along the pressure leading edge and suction surface trailing edge of the impeller blade. High wear is observed along the leading edge of the blade, extending rearward from the bottom. The wear rate significantly declines in the upper portion of the leading edge, while the center of the blade shows irregular, localized wear patterns.
The primary cause lies in the lower velocity of the particles at the bottom of the blade’s leading edge, leading to prolonged accumulation and severe wear in this region. In the middle of the blade, particles acquire higher kinetic energy, resulting in more random impacts and uneven wear patterns. At the trailing edge of the blade, the particle velocity peaks, significantly intensifying the blade impact and causing extensive high wear rates. As the particle surface area increases, the contact area between the particles and the fluid expands, leading to a higher wear rate on the blade. As illustrated in Figure 13a, the area with a high wear rate, for elongated cylindrical particles, is the most extensive, whereas for spherical particles, the high wear rate area is the smallest. From Figure 13b, it is evident that the impact of the particle shape on the wear rate distribution in the volute mirrors that of the impeller. Long cylindrical particles demonstrate the highest wear rate and deepest wear depth on the volute, followed by short cylindrical particles with a lower wear rate. Spherical particles exhibit the lowest wear rate, while mixed particles display intermediate wear rates between short cylindrical and spherical particles.

5. Analysis of the Influence of the Tip Clearance Size on Performance

To prevent friction between the dynamic and static components of open-type submersible sewage pumps, a specially designed tip clearance gap induces leakage flow of the working medium under a pressure differential, forming a tip leakage vortex that reduces the hydraulic performance. Thus, the following research investigates the effect of varying tip clearance thicknesses as the sole variable, employing long cylindrical particles in terms of the design condition, with a particle volume fraction of 1%, to simulate the influence of the tip clearance thickness change on the solid–liquid two-phase flow dynamics and energy dissipation by the pump. Various tip clearance parameter schemes are presented in Table 5.

5.1. The Influence of Different Tip Clearances on Tip Leakage Vortex Structure

Figure 14 presents the color coding for the TKE that represents the distribution of the tip leakage vortices in a submersible sewage pump for varying tip clearances. As depicted, increasing the gap thickness, enlarges the leakage vortex at the tip clearance inlet and markedly expands the region of irregular small-scale vortices along the blade. Conversely, the TKE distribution within the vortex gradually diminishes with the increasing tip clearance thickness. This is attributed to the fact that enlarging the tip clearance thickness allows more fluid to reach the volute directly through the clearance, diminishing the driving effect of the blade on the fluid and leading to increased resistance as the fluid moves past the blade and, consequently, more turbulence. Additionally, the fluid velocity decreases, exacerbating energy dissipation. However, the reduced fluid velocity concurrently lowers the turbulent kinetic energy within the vortex.
Figure 15 examines the analysis of the tip leakage flow velocities for different tip clearances, to assess the impact on the performance of the submersible sewage pump. The tip leakage flow mainly travels from the pressure side to the suction side of the impeller, with its velocity affected by varying tip clearance thicknesses. Furthermore, under a constant tip clearance thickness, the tip leakage flow velocity increases progressively from the leading edge to the trailing edge of the blade, attaining its highest value at the trailing edge of the pressure surface, as illustrated. As the tip clearance thickness increases, the velocity distribution at identical positions along the tip clearance diminishes gradually, which aligns with the previously discussed TKE distribution analysis. This trend primarily stems from the augmented tip clearance thickness enabling greater fluid passage, thereby reducing the driving effect of the blade on the fluid.

5.2. The Influence of Different Tip Clearances on the Flow Pattern

The velocity and TKE distributions are used to describe the fluid motion and energy distribution within the pump, directly representing the internal flow pattern of the submersible sewage pump. Figure 16 presents the velocity distribution contour for various tip clearances. With the increasing tip clearance thickness, the maximum fluid velocity within the impeller and volute of the submersible sewage pump decreases gradually. The variation in the tip clearance significantly influences the fluid velocity within the impeller passages, particularly along the suction surface of the blades. A progressive increase in the fluid velocity at the suction surface of the blade is noted with an increase in the tip clearance. Moreover, under the same tip clearance, the maximum velocity region of the impeller in the submersible sewage pump is situated at the trailing edge of the pressure surface, with the greatest velocity region in the volute appearing in section II; the minimum velocity region of the impeller is observed at the impeller inlet, and a similar region is identified at the volute outlet.
Figure 17 depicts the distribution of the TKE within the impeller and volute in the submersible sewage pump for various tip clearances. The diagram shows that as the tip clearance enlarges, the TKE within the impeller in the submersible sewage pump gradually increases, especially in the region of the impeller flow channel adjacent to the volute tongue. This indicates that turbulence generated at the impeller in the submersible sewage pump increases with a greater tip clearance, resulting in increased energy loss and reduced operational stability. Meanwhile, the TKE value in the V to VIII section of the volute shows a gradual decrease with the increasing tip clearance.

5.3. The Influence of Different Tip Clearances on the Wear Characteristics

Figure 18 and Figure 19 depict the distribution of the particles and wear rate across the primary flow components for different tip clearances in terms of the pump. It is apparent from the figures that as the tip clearance increases, particle aggregation becomes more pronounced, accompanied by a gradual reduction in the particles within the impeller flow channel. Figure 19 reveals a consistent wear rate distribution trend on the impeller blades and volute in the submersible sewage pump for various tip clearances. The highest wear rates on the blades predominantly occur at the leading edge of the pressure region and the trailing edge of the suction region, while for the volute, they are concentrated in sections I and IV. Furthermore, as the tip clearance increases incrementally, the wear rates on the impeller blades and volute both exhibit a gradual increase, with the most notable increase observed at the trailing edge of the suction surface of the blade. The observed phenomenon can be explained by the alteration in the flow pattern and velocity distribution of the solid–liquid two-phase fluid between the impeller and volute caused by the increased tip clearance. This change induces greater eddy currents and turbulence, thereby intensifying the particle impact and surface friction on the impeller and volute surfaces, consequently elevating the wear rate. Moreover, the increased tip clearance facilitates greater direct contact between the particles and pump surfaces, further augmenting the wear severity.

5.4. The Influence of Different Tip Clearances on the Energy Loss Mechanism

During the design and operation of open-type submersible sewage pumps, the tip clearance thickness is closely linked to the pump’s energy conversion efficiency, directly influencing internal flow dynamics and energy conversion processes. Energy losses occur within the primary components, such as the impeller and volute, during the flow process of submersible sewage pumps. Using energy balance equation theory, the alterations in the energy loss for the flow components are evaluated, highlighting the relationship between the tip clearance thickness and the energy conversion process inside the pump and providing a reference for selecting the most suitable tip clearance in engineering practice.
Lu et al. [29] used the unsteady Reynolds-Averaged Navier–Stokes method to analyze local, unsteady losses using the energy equation. It was observed that, during non-rotating stall conditions, changes in the average flow kinetic energy had negligible influence on local energy losses within the pump; instead, the majority of losses were converted into turbulent kinetic energy. Building upon this finding and referencing Wilhelm et al.’s work [30], a simplified equation can be derived:
p ¯ ρ g + z + 1 2 g v f   i ¯ 2 ρ g v f   j ¯ x j = x j ( v f   i ¯ ) ( ρ v f   i v f   j ¯ ) p 1 + μ x j ( v f   i ¯ ) ( D ij ¯ ) p 2 ( ρ v f   i v f   j ¯ ) v f   i ¯ x j p 3 μ ( D ij ¯ ) v f   i ¯ x j p 4
D ij = 2 S ij = v fi x j + v f   j x i
ρ v f   i v f   j ¯ = μ t ( v f   i ¯ x j + v f   j ¯ x i ) 2 3 δ ij ρ k
where ρ represents the fluid density; g is the acceleration of gravity; z signifies the height of the fluid mass; μ stands for the fluid’s dynamic viscosity, in units of N·s/m2; νfi denotes the fluid velocity component in the direction of i; v fi ¯ is the average velocity component in the direction of i and v fi is the fluctuating component of the velocity in the direction of i, both with units of m/s; ρ v f   i v f   j ¯ represents the Reynolds stress, which can be expressed by the average velocity gradient and turbulent viscosity μt; and μt serves as a parameter in the SST k-ω turbulence model.
By adding up the local energy losses, the total energy dissipation within the pump is determined. The energy balance equation for the pump is formulated as:
L = g Q Δ H = V x j v f   i ¯ ρ v f   i v f   j ¯ d V L 1 V μ v f   i ¯ D ij ¯ x j d V L 2       + V ( ρ v f   i v f   j ¯ ) v f   i ¯ x j d V L 3 + V μ D ij ¯ v f   i ¯ x j d V L 4
where Q represents the mass flow, where the unit is kg/s; L represents the total energy loss term, where the unit is W; L1 and L2 denote the diffusion of kinetic energy through Reynolds stress and viscous stress, respectively, known as the Reynolds stress transfer term and viscosity transfer term. L3 represents the conversion of average kinetic energy into turbulent kinetic energy, termed turbulent kinetic energy generation. The TKE is eventually converted into heat through viscous dissipation by small-scale turbulence. L4 denotes the viscous dissipation of the average kinetic energy, referred to as the viscous dissipation term [31].
Illustrated in Figure 20 is the distribution of various energy loss types within the submersible sewage pump for different tip clearances. It can be observed from the diagram that under different tip clearance thicknesses, the turbulent kinetic energy generation term L3 has the highest value, while the viscous dissipation term L4 has the lowest value in terms of the two-phase flow condition in the submersible sewage pump. The value of L3 significantly exceeds that of the other terms, aligning with the transient flow characteristics observed in terms of the clear water conditions described earlier. This underscores that the generation of irregular turbulence under solid–liquid two-phase flow conditions constitutes the primary energy loss mechanism during the operation of submersible sewage pumps. The influence on the energy loss in the submersible sewage pump is primarily attributed to the turbulent kinetic energy generation term L3. As the tip clearance thickness gradually increases, the diagram illustrates a corresponding increase in the total energy loss term L within the submersible sewage pump, consistent with the trend observed in L3. The analysis indicates that higher fluid velocity differentials passing through the blade clearance are a result of the increase in the tip clearance, thereby elevating the turbulent kinetic energy generation term L3 and, consequently, increasing the total energy loss L. This also shows that to investigate the connection between the tip clearance change and energy conversion in the submersible sewage pump, it is necessary to further analyze the turbulent kinetic energy generation term L3 for each flow component.
Figure 21 depicts the variation in the turbulent kinetic energy generation term L3 among the principal flow components of the submersible sewage pump under solid–liquid two-phase flow conditions and different tip clearances. The diagram reveals that as the tip clearance thickness increases, the L3 values in regard to the tip clearance (BTC) and impeller (IMP) progressively increase, whereas the L3 value in terms of the volute (VOL) decreases gradually. The L3 value in the IMP exhibits the most significant change, increasing by 36.8 W. Conversely, L3 in the VOL decreases by 34.7 W for different tip clearances. The change in L3 for the BTC is less pronounced compared to the IMP and VOL, increasing by only 11.7 W. These results highlight that the variations in the tip clearance thickness primarily affect L3 in the impeller and volute. The slight increase in L3 in the BTC with the increasing tip clearance thickness is primarily attributed to enhanced leakage flow, intensifying the non-uniform fluid velocity distribution at the tip clearance and, consequently, augmenting L3. Due to the minor volume in terms of the tip clearance compared to the impeller and volute, changes in L3 are relatively insignificant. Additionally, an augmentation in the tip clearance directly enlarges the gap between the blade and the volute, and more fluid flows directly into the volute without the pressure and growth of the impeller, resulting in a more dramatic change in the fluid velocity gradient in the impeller, an irregular eddy current and turbulence increase, resulting in a considerable increase in the L3 value and a sharp increase in the energy loss. As the tip clearance thickness increases within the volute, L3 exhibits a gradual decrease, signifying enhanced fluid flow stability, reduced turbulence distribution, and diminished energy loss within this component.

5.5. The Influence of Different Tip Clearances on Local Energy Loss

To analyze the energy loss component along the cross-section, a local loss coefficient li is introduced to quantify the magnitude of the term pi (Equation (22)) within specific regions of the pump body. Its definition is as follows:
l i = p i ( ρ u 2 3 ) / α D
where α is a constant with a value of 5.556; u2 is defined as the circumferential velocity of the exit to the impeller, where the unit is m/s; D represents the impeller diameter, measured in meters; l1 signifies the local loss coefficient for the Reynolds stress transfer term; l2 denotes the local loss coefficient for the viscosity transfer term; l3 indicates the local loss coefficient for the turbulent kinetic energy generation term; and l4 represents the local loss coefficient for the viscous dissipation term.
Figure 22 depicts the distribution of the local loss coefficient l3 for the turbulent kinetic energy generation term in the central region of the impeller and volute for different tip clearances. The regions of elevated l3 are predominantly observed in the impeller inlet area, volute tongue area, and the junction of the trailing edge of the blade and the volute. The change in l3 in the central flow channel area between the two blades near the volute tongue is more obvious, which indicates that the increasing tip clearance leads to a gradual diffusion of the turbulence from the impeller inlet to the center of the flow channel. This results in increased flow instability within the impeller, accompanied by higher operational energy losses.
The impeller and tip clearance in the submersible sewage pump have a cross-section along the x-axis and y-axis directions. The positional distribution of the four sections (P1 to P4) is illustrated in Figure 23.
The distribution of the local loss coefficient l3 related to the turbulent kinetic energy generation term is shown in Figure 24 across the four sections (P1 to P4), highlighting the intensity of turbulent kinetic energy generation in the local area. According to the diagram, under different tip clearance thicknesses, the high l3 regions in sections P1 through to P4 are primarily concentrated around the edge of the impeller blade. This situation is particularly evident in terms of the distribution of l3 at the tip clearance, where a sharp decrease in l3 is observed in regions not in contact with the blade. This reduction occurs due to the absence of blade contact, mitigating the flow separation phenomena and, subsequently, reducing the turbulent energy. With the increasing tip clearance thickness, the region with elevated l3 values at the tip clearance steadily expands towards the midsection of the impeller passage. The most noticeable diffusion changes are observed in the X1 to X4 and Y1 to Y4 regions, as represented in the figure. It is apparent that in these regions of the impeller, the l3 value increases gradually with the thicker tip clearance, accompanied by an expanding diffusion range. This demonstrates that as the tip clearance increases, the leakage flow increases, impacting the fluid velocity and pressure distribution. Concurrently, the turbulence intensity within the impeller flow channel shows an increasing trend, consistent with the earlier findings. Moreover, under an identical tip clearance thickness, the l3 distribution is notably more intricate in the X3, Y2, and Y3 regions, with higher l3 values compared to the other five regions marked by red circles. The main reason is that the cross-sections of P2 and P3 are adjacent to the volute tongue, and the impeller rotation interacts with the tongue to produce a dynamic and static coherence effect. As a result of the interference effect, the flow properties in the central region of the impeller channel undergo considerable changes, forming more complex turbulence.

6. Conclusions

In the present study, the CFD–DEM coupling methodology is employed to investigate how varying particle parameters and tip clearance thicknesses affect the performance of an open impeller submersible sewage pump, operating under solid–liquid two-phase flow conditions. Utilizing the energy balance equation, the tip leakage vortex structure, the internal flow pattern, and wear characteristics of the pump for various tip clearance thicknesses are compared and analyzed, including the effect on energy conversion within the pump. The following key findings are summarized:
  • At a low flow rate, a substantial low-speed region develops on the suction surface of the impeller, leading to the formation of minor vortices. The highest velocity region is located near the second section of the volute, and this peak velocity diminishes gradually with an increasing flow rate. Furthermore, higher flow rates lead to more dispersed particle distribution, increasing the wear rates, particularly at the leading edge of the blade. Conversely, the wear rates impacting the volute decrease;
  • The relationship between the particle parameters and pump wear characteristics is compared and analyzed. It is observed that increasing the particle volume fraction leads to decreased inlet particle velocity. Meanwhile, the extent of the wear on the blade pressure leading edge, the suction surface trailing edge, and the volute area become larger. The wear rate is highest for long cylindrical particles, with spherical particles exhibiting the lowest wear rate. A low particle velocity is observed at the leading edge of the blade, which prolongs the particle contact time and causes significant wear in that area. The velocity of the particles at the trailing edge of the blade reaches the maximum, the impact effect is significantly enhanced, and a large area related to a high wear rate area appears;
  • Through comparing the tip clearance thickness of 0.4 mm, 0.7 mm, and 1.0 mm, changes in the internal flow and wear rates of the flow passage components were analyzed. Observations indicate that a greater tip clearance thickness intensifies the leakage vortex at the tip clearance inlet. The tip leakage flow velocity escalates from the leading to the trailing edge of the blade, peaking at the trailing edge of the pressure surface, which also makes the TKE on the impeller gradually increase. Moreover, an increased tip clearance exacerbates particle aggregation, leading to a gradual reduction in the particle concentration within the impeller channel, and the wear rate of the impeller and volute gradually increase;
  • By applying the energy balance equation, the effect of varying tip clearance thicknesses on the energy dissipation in the pump is assessed under solid–liquid two-phase flow conditions. The results indicate that the turbulent kinetic energy generation term L3 exhibits the highest magnitude, while the viscous dissipation term L4 exhibits the smallest proportion, highlighting turbulence as the primary source of energy dissipation in submersible sewage pump operations. Simultaneously, regions with elevated values in terms of the local loss coefficient l3 from turbulent kinetic energy generation are primarily concentrated at the impeller inlet and blade edges. Under the same tip clearance thickness, dynamic and static interference results in higher l3 values in regions X3, Y2, and Y3 near the volute tongue, indicating more complex flow dynamics.

Author Contributions

Conceptualization, G.P. and J.Y.; methodology, L.M. and Y.L.; software, J.Y., Y.L. and Z.W.; validation, G.J.; formal analysis, J.Y. and S.H.; investigation, L.M.; resources, H.C.; data curation, J.Y. and L.M.; writing—original draft preparation, J.Y. and Z.W.; writing—review and editing, J.Y. and L.M.; supervision, G.P.; project administration, H.C.; funding acquisition, G.P. All authors have read and agreed to the published version of the manuscript.

Funding

The work was sponsored by the National Natural Science Foundation of China (Grant No. 52409114, No. 52079058), the China Postdoctoral Science Foundation (Grant No. 2024M751272), the Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering (Grant No. sklhse-2024-E-02), the Natural Science Foundation of Jiangsu Province (Grant No. BK20230011), and the Taizhou Science and Technology Support Plan (Grant No. TG202336).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

Author Guangjie Peng was employed by the company Taizhou Longda Submersible Pump Co., Ltd. and author Hao Chang was employed by the company Leo Group Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Nomenclature

English Symbols
b1Impeller outlet width, mm
b2Volute inlet width, mm
b3Tip clearance width, mm
b4Rear cover plate width, mm
φWrapping angle, °
ZNumber of blades
nsSpecific speed
HdDesign head, m
QdDesign flow rate, m3/h
ρfThe fluid density, kg/m3
ρsThe particle density, kg/m3
νfThe fluid velocity vector
gGravity acceleration, m/s2
ωsAngular velocity of the particle, rad/s
IsThe moment of inertia of the particle, kg·m2
msParticle mass, kg
CsaffSlip-shear lift coefficient
wcThe vorticity of the fluid
VThe linear velocity of the particle relative to the fluid, m/s
ΩAngular velocity of the fluid relative to the particles, rad/s
Res,rThe rotational Reynolds number
Res,αThe particle Reynolds number
μPoisson’s ratio
YThe equivalent elastic modulus
QvThe amount of wear, mm3
HvThe hardness of the material surface, N/mm2
KThe wear constant
v fi ¯ The average velocity component in the direction of i, m/s
v fi The fluctuating component of the velocity in the direction of i, m/s
LThe total energy loss term, W
L1The Reynolds stress transfer term, W
l1The local loss coefficient for the Reynolds stress transfer term
L2The viscosity transfer term, W
l2The local loss coefficient for the viscosity transfer term
L3The turbulent kinetic energy generation term, W
l3The local loss coefficient for the turbulent kinetic energy generation term
L4The viscous dissipation term, W
l4The local loss coefficient for the viscous dissipation term
Acronyms
CFDComputational fluid dynamics
DEMDiscrete element method
RANSReynolds-averaged Navier–Stokes
SSTShear stress transport
TKETurbulent kinetic energy
IMPImpeller
BTCThe blade tip clearance
VOLVolute

References

  1. Wu, B.; Wang, X.-L.; Liu, H.; Xu, H.-L. Numerical simulation and analysis of solid-liquid two-phase three-dimensional unsteady flow in centrifugal slurry pump. J. Central S. Univ. 2015, 22, 3008–3016. [Google Scholar] [CrossRef]
  2. Luo, X.; Cao, S.; Pei, Y.; Xiao, X.; Wang, L.; Xue, Z.; Yang, C. Energy conversion characteristics of an in-line mixing pump in solid-liquid two-phase dynamic mixing. Arab. J. Sci. Eng. 2023, 48, 16953–16971. [Google Scholar] [CrossRef]
  3. Tan, M.; Shao, C.; Wu, X.; Liu, H.; Pan, B. Numerical study on cylindrical particle solid-liquid two-phase flow in a slurry pump. Comput. Part. Mech. 2024, 11, 921–933. [Google Scholar] [CrossRef]
  4. Wang, Y.; Tao, R.; Han, C.; Li, W.; He, T.; Zhu, Z. Numerical study on flow and wear characteristics of dense fine particle solid-liquid two-phase flow in centrifugal pump. AIP Adv. 2022, 12, 045109. [Google Scholar] [CrossRef]
  5. Wang, Y.; Zu, Y.; He, T.; Deng, X.; Li, Y. Effect of the volume concentration of binary mixed particles on the flow and wear characteristics of centrifugal pumps. Processes 2024, 12, 1185. [Google Scholar] [CrossRef]
  6. Zheng, S.; Yang, C.; Yan, C.; Chai, M.; Sun, Z. Study on fiber clogging mechanism in sewage pump based on CFD-DEM simulation. Energies 2022, 15, 1598. [Google Scholar] [CrossRef]
  7. Wang, R.; Guan, Y.; Jin, X.; Tang, Z.; Zhu, Z.; Su, X. Impact of particle sizes on flow characteristics of slurry pump for deep-sea mining. Shock Vib. 2021, 2021, 6684944. [Google Scholar] [CrossRef]
  8. Cao, W.; Wang, H.; Tang, J. Study on wear characteristics of a guide vane centrifugal pump based on CFD-DEM. J. Mar. Sci. Eng. 2024, 12, 593. [Google Scholar] [CrossRef]
  9. Wang, H.; Tan, Z.; Kuang, S.; Yu, A. Numerical modeling and analysis of particle-fluid flow and wall erosion in centrifugal slurry pumps under different solid concentrations. Powder Technol. 2022, 410, 117861. [Google Scholar] [CrossRef]
  10. Guo, X.; Shi, G.; Xiao, Y.; Ye, X. A study of the relationship between sand movement and flow field distribution and wear causes in a multiphase pump. J. Mar. Sci. Eng. 2024, 12, 1203. [Google Scholar] [CrossRef]
  11. Peng, G.; Chen, Q.; Bai, L.; Hu, Z.; Zhou, L.; Huang, X. Wear mechanism investigation in a centrifugal slurry pump impeller by numerical simulation and experiments. Eng. Fail. Anal. 2021, 128, 105637. [Google Scholar] [CrossRef]
  12. Cheng, W.; Fan, H.; Cheng, W.; Shao, C. Investigation on wear induced by solid-liquid two-phase flow in a centrifugal pump based on EDEM-Fluent coupling method. Flow Meas. Instrum. 2024, 96, 102542. [Google Scholar] [CrossRef]
  13. Wang, K.; Liu, H.; Wang, L.; Guo, P.; Wang, Y.; Yang, J. Effect of particle size on vortex structure and erosion behavior of semi-open centrifugal pump. Energy 2024, 293, 130576. [Google Scholar] [CrossRef]
  14. Yang, Y.; Wang, H.; Hu, Q.; Ji, L.; He, Z.; Shi, W.; Song, X.; Zhou, L. Two-phase flow investigation of sewage pumps with different tip clearance via computational fluid dynamics and multi-factor ANOVA. Eng. Appl. Comput. Fluid Mech. 2024, 18, 2322514. [Google Scholar] [CrossRef]
  15. Ji, L.; Li, W.; Shi, W.; Zhou, L.; Agarwal, R. Experimental study of pressure pulsation in a mixed-flow pump with different tip clearances based on wavelet analysis. Shock Vib. 2020, 2020, 9041686. [Google Scholar] [CrossRef]
  16. Lin, Z.; Yang, F.; Guo, J.; Jian, H.; Sun, S.; Jin, X. Leakage flow characteristics in blade tip of shaft tubular pump. J. Mar. Sci. Eng. 2023, 11, 1139. [Google Scholar] [CrossRef]
  17. Yang, Y.; Wang, H.; Chen, X.; Li, J.; Jiao, W.; Ji, L.; He, Z.; Shi, W.; Zhou, L. Leakage vortices and energy characteristics of semi-open sewage pump with various blade tip alterations. Phys. Fluids 2024, 36, 035170. [Google Scholar] [CrossRef]
  18. Guo, X.; Jiang, C.; Qian, H.; Zhu, Z. The influence of tip clearance on the performance of a high-speed inducer centrifugal pump under different flow rates conditions. Processes 2023, 11, 239. [Google Scholar] [CrossRef]
  19. Peng, G.; Hong, S.; Chang, H.; Fan, F.; Zhang, Y.; Shi, P. Numerical and experimental research on the influence of clearance between impeller and cover on the pump performance. Mechanika 2022, 28, 67–72. [Google Scholar] [CrossRef]
  20. Li, W.; Ji, L.; Li, E.; Zhou, L.; Agarwal, R.K. Effect of tip clearance on rotating stall in a mixed-flow pump. J. Turbomach. 2021, 143, 091013. [Google Scholar] [CrossRef]
  21. Han, Y.; Bai, L.; Du, D.; Shi, W.; Zhou, L. Effects of tip clearance on energy performance of three-stage electrical submersible pump. Geoenergy Sci. Eng. 2023, 226, 211696. [Google Scholar] [CrossRef]
  22. Rubinow, S.I.; Keller, J.B. The transverse force on a spinning sphere moving in a viscous fluid. J. Fluid Mech. 1961, 11, 447–459. [Google Scholar] [CrossRef]
  23. Xu, L.; Zhang, Q.; Zheng, J.; Zhao, Y. Numerical prediction of erosion in elbow based on CFD-DEM simulation. Powder Technol. 2016, 302, 236–246. [Google Scholar] [CrossRef]
  24. Fluent Inc. Ansys Fluent 22 User’s Guide; Fluent Inc.: Lebanon, NH, USA, 2022. [Google Scholar]
  25. Saffman, P.G. The lift on a small sphere in a slow shear flow. J. Fluid Mech. 1965, 22, 385–400. [Google Scholar] [CrossRef]
  26. Crowe, C.T. Multiphase Flows with Droplets and Particles, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2012. [Google Scholar]
  27. Archard, J.F. Contact and rubbing of flat surfaces. J. Appl. Phys. 1953, 24, 981–988. [Google Scholar] [CrossRef]
  28. Hong, S.; Peng, G.; Chang, H.; Yu, D.; Ji, G.; Ma, L.; Wang, Z. Numerical study of erosion characteristics in U-shaped elbow and slurry pump. Phys. Fluids 2024, 36, 073343. [Google Scholar] [CrossRef]
  29. Lu, G.; Zuo, Z.; Liu, D.; Liu, S. Energy balance and local unsteady loss analysis of flows in a low specific speed model pump-turbine in the positive slope region on the pump performance curve. Energies 2019, 12, 1829. [Google Scholar] [CrossRef]
  30. Wilhelm, S.; Balarac, G.; Métais, O.; Ségoufin, C. Analysis of head losses in a turbine draft tube by means of 3D unsteady simulations. Flow Turbul. Combust. 2016, 97, 1255–1280. [Google Scholar] [CrossRef]
  31. Yang, G.; Shen, X.; Shi, L.; Meng, J.; Luo, W.; Zhang, D.; van Esch, B.P.M. Unsteady numerical investigations of the effect of guide vane openings on the hydrodynamic characteristics under stall conditions in a pump-turbine pump mode. Energy Convers. Manag. 2023, 293, 117499. [Google Scholar] [CrossRef]
Figure 1. Models of the submersible sewage pump.
Figure 1. Models of the submersible sewage pump.
Water 16 02845 g001
Figure 2. The particle shapes.
Figure 2. The particle shapes.
Water 16 02845 g002
Figure 3. The diagram of coupling.
Figure 3. The diagram of coupling.
Water 16 02845 g003
Figure 4. The diagram of polyhedral meshing of submersible sewage pump.
Figure 4. The diagram of polyhedral meshing of submersible sewage pump.
Water 16 02845 g004
Figure 5. A comparison of the pump head and efficiency for different grid groups.
Figure 5. A comparison of the pump head and efficiency for different grid groups.
Water 16 02845 g005
Figure 6. Velocity vector distribution for different flow rates.
Figure 6. Velocity vector distribution for different flow rates.
Water 16 02845 g006
Figure 7. Distribution of TKE in solid–liquid two-phase flow conditions for different flow rates.
Figure 7. Distribution of TKE in solid–liquid two-phase flow conditions for different flow rates.
Water 16 02845 g007
Figure 8. The particle distribution in the submersible sewage pump for different flow rates.
Figure 8. The particle distribution in the submersible sewage pump for different flow rates.
Water 16 02845 g008
Figure 9. The wear rate distribution in the submersible sewage pump for different flow rates.
Figure 9. The wear rate distribution in the submersible sewage pump for different flow rates.
Water 16 02845 g009
Figure 10. The particle distribution for different particle volume fractions.
Figure 10. The particle distribution for different particle volume fractions.
Water 16 02845 g010
Figure 11. The wear rate distribution on the flow components for different particle volume fractions.
Figure 11. The wear rate distribution on the flow components for different particle volume fractions.
Water 16 02845 g011aWater 16 02845 g011b
Figure 12. The particle distribution for different particle shapes.
Figure 12. The particle distribution for different particle shapes.
Water 16 02845 g012
Figure 13. The wear rate distribution for different particle shapes.
Figure 13. The wear rate distribution for different particle shapes.
Water 16 02845 g013
Figure 14. The distribution of the tip leakage vortex structure for different tip clearances.
Figure 14. The distribution of the tip leakage vortex structure for different tip clearances.
Water 16 02845 g014
Figure 15. The distribution of the tip leakage flow velocity for different tip clearances.
Figure 15. The distribution of the tip leakage flow velocity for different tip clearances.
Water 16 02845 g015
Figure 16. The distribution of velocity for different tip clearances.
Figure 16. The distribution of velocity for different tip clearances.
Water 16 02845 g016
Figure 17. The TKE distribution in the pump for different tip clearances.
Figure 17. The TKE distribution in the pump for different tip clearances.
Water 16 02845 g017
Figure 18. The distribution of the particles for different tip clearances.
Figure 18. The distribution of the particles for different tip clearances.
Water 16 02845 g018
Figure 19. The wear rate distribution on the flow components for different tip clearances.
Figure 19. The wear rate distribution on the flow components for different tip clearances.
Water 16 02845 g019aWater 16 02845 g019b
Figure 20. The energy loss distribution for different tip clearances.
Figure 20. The energy loss distribution for different tip clearances.
Water 16 02845 g020
Figure 21. The L3 distribution of the flow parts for different blade tip clearances.
Figure 21. The L3 distribution of the flow parts for different blade tip clearances.
Water 16 02845 g021
Figure 22. The l3 distribution in the submersible sewage pump for different blade tip clearances.
Figure 22. The l3 distribution in the submersible sewage pump for different blade tip clearances.
Water 16 02845 g022
Figure 23. The cross-section diagram of blade tip clearance and impeller.
Figure 23. The cross-section diagram of blade tip clearance and impeller.
Water 16 02845 g023
Figure 24. The l3 distribution in blade tip clearance and impeller sections for different clearance widths.
Figure 24. The l3 distribution in blade tip clearance and impeller sections for different clearance widths.
Water 16 02845 g024
Table 1. The geometrical parameters of the submersible sewage pump.
Table 1. The geometrical parameters of the submersible sewage pump.
Geometrical ParametersValueGeometrical ParametersValue
Inlet pipe diameter Di/mm58Outlet pipe diameter Do/mm50
Impeller inlet diameter D1/mm71Wrapping angle φ/°194.57
Impeller outlet diameter D2/mm180Base circle diameter of volute D3/mm180
Impeller outlet width b1/mm7.62Volute inlet width b2/mm26
Number of blades Z3Diameter of volute outlet D4/mm50
Blade inlet setting angle β1/°16.3Tip clearance width b3/mm0.7
Blade outlet setting angle β2/°13.15Rear cover plate width b4/mm0.5
Table 2. Solid–liquid two-phase flow research plan.
Table 2. Solid–liquid two-phase flow research plan.
VariablesFlow Conditions (Q/Qd)Particle Volume Fractions (%)Particle Shapes
Flow conditions0.6, 1.0, 1.41Long cylindrical particle
Particle volume fractions1.01, 2, 3Long cylindrical particle
Particle shapes1.01Spherical, short cylinder particle, long
cylindrical particle, mixed particle
Table 3. Key parameter settings for the submersible sewage pump in Fluent and EDEM.
Table 3. Key parameter settings for the submersible sewage pump in Fluent and EDEM.
Fluid propertiesParametersUnitValue
Densitykg/m3998.2
Import speedm/s1.0514
Free outflow//
Particle propertiesDensitykg/m32650
Poisson’s ratio/0.3
Young’s modulusPa1.9 × 109
Restitution coefficient/0.4
Static friction coefficient/0.7
Wall propertiesDensitykg/m37825
Poisson’s ratio/0.285
Young’s modulusPa2 × 1011
Restitution coefficient/0.5
Static friction coefficient/0.4
Table 4. Grid independence analysis.
Table 4. Grid independence analysis.
SchemeGrid NumberHead (m)Efficiency (%)
140295731.4234.85
252882531.1434.23
365254330.7133.95
473454230.6833.89
587551130.6733.82
Table 5. Simulation scheme for different blade tip clearance parameters.
Table 5. Simulation scheme for different blade tip clearance parameters.
SchemeFlow Rate (Q/Qd)Particle Volume Fraction (%)Particle ShapeBlade Tip Clearance (mm)
11.01Long cylindrical particles0.4
21.01Long cylindrical particles0.7
31.01Long cylindrical particles1.0
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

Peng, G.; Yang, J.; Ma, L.; Wang, Z.; Chang, H.; Hong, S.; Ji, G.; Lou, Y. Research on the Influence of Particles and Blade Tip Clearance on the Wear Characteristics of a Submersible Sewage Pump. Water 2024, 16, 2845. https://doi.org/10.3390/w16192845

AMA Style

Peng G, Yang J, Ma L, Wang Z, Chang H, Hong S, Ji G, Lou Y. Research on the Influence of Particles and Blade Tip Clearance on the Wear Characteristics of a Submersible Sewage Pump. Water. 2024; 16(19):2845. https://doi.org/10.3390/w16192845

Chicago/Turabian Style

Peng, Guangjie, Jinhua Yang, Lie Ma, Zengqiang Wang, Hao Chang, Shiming Hong, Guangchao Ji, and Yuan Lou. 2024. "Research on the Influence of Particles and Blade Tip Clearance on the Wear Characteristics of a Submersible Sewage Pump" Water 16, no. 19: 2845. https://doi.org/10.3390/w16192845

APA Style

Peng, G., Yang, J., Ma, L., Wang, Z., Chang, H., Hong, S., Ji, G., & Lou, Y. (2024). Research on the Influence of Particles and Blade Tip Clearance on the Wear Characteristics of a Submersible Sewage Pump. Water, 16(19), 2845. https://doi.org/10.3390/w16192845

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