Investigation on the Flow Behavior of Side Channel Pumps Based on Vortex Identification

The momentum flow exchange between the impeller and side channel produces highly turbulent flows in side channel pumps. The turbulent flows feature complex patterns of vortex structures that are partly responsible for the dissipation of energy losses and unsteady pressure pulsations. The concept of turbulent flows in side channel pumps requires a reliable vortex identification criterion to capture and predict the effects of the vortex structures on the performance. For this reason, the current study presents the application of the new Ω-criterion to a side channel pump model in comparison with other traditional methods such as Q and λ2 criteria. The 3D flow fields of the pump were obtained through unsteady Reynolds-averaged Navier-Stokes (RANS) simulations. Comparative studies showed that the Ω-criterion identifies the vortex of different intensities with a standard threshold, Ω=0.52. The Q and λ2 criteria required different thresholds to capture vortex of different intensities thus leads to subjective errors. Comparing the Ω-criterion intensity on different planes with the entropy losses and pressure pulsation, the longitudinal vortex plays an important role in the momentum exchange development which increases the head performance of the pump. However, the rate of exchange is impeded by the axial and radial vortices restricted in the impeller. Therefore, the impeller generates the highest entropy loss and pressure pulsation intensities which lower the output efficiency. Finally, the findings provide a fundamental background to the morphology of the vortex structures in the turbulent flows which can be dependent upon for efficiency improvement of side channel pumps.


Introduction
Description of the evolution of vortex structures explains the flow theory behind turbulent flows. Thus, vortex structure identification has been extensively identified and studied not only to model the turbulence but reveal turbulence production and dissipation in turbomachines. To date, there exist no general identification method of vortex structures and the available methods are case dependent [1]. Most of the commonly used identification methods in turbomachinery depend on a key parameter called the velocity gradient tensor. The velocity gradient tensor ∇ � V can be defined as: where P, Q and R are the invariants of the velocity gradient tensor. The eigenvalues, λ of ∇ � V is determined by calculating the cubic characteristic polynomial, det ∇ − → V − I = 0 , here I is the identity matrix.
Additionally, ∇ � V can be decomposed into a symmetric part, A and antisymmetric part, B.
The popular Q-criterion proposed by Hunt et al. [2] identifies vortex as a fluid region with a positive second invariant Q of the velocity gradient tensor ∇ � V for incompressible flows. Q-criterion can be expressed as the difference between the magnitude of the rate of shear (1) 3 + P 2 + Q + R = 0, Furthermore, Jeong and Hussain [3] postulated the λ 2 -criterion which also depends on ∇ � V . This method seems to reduce the pressure but removes the effects of unsteady strain and viscosity by taking the gradient of the incompressible Navier-Stokes equation. The equation decomposes into two parts A and B representing symmetric part and antisymmetric part respectively.
Concerning the λ 2 -criterion, vortex exists in the domain where there are two negative eigenvalues of A 2 + B 2 i.e. 1 ≥ 2 ≥ 3 , 2 < 0 . The relation between negative eigenvalues of A 2 + B 2 and Q can be expressed as: This method depends on the minimum value of λ 2 extracted to represent the vortex core position. The vortex core lines revealed that using this method usually are not continuous because they are composed of several vortex structure segments.
To capture the vortex structures in turbulent flows, the Q and λ 2 criteria have been commonly used. Wang et al. [4] analyzed the instantaneous vortical structures depicted by the vortex core region using the λ 2 -criterion for different turbulence models in side channel pump. They concluded that the evolution of most vortical structures in the circumferential direction are responsible for the process of fluid exchange. Minakov et al. [5] showed the typical vortex structures in the operation of Francis turbines behind the runner using the Q-criterion. It is revealed that large vortices are found behind the runner and are the major source of the pressure pulsations. Asim et al. [6] depicted the variations in the downstream of the tongue region of a centrifugal pump. They indicated that the positive Q-criterion values correspond to higher vorticity magnitude than the shear strain-rate. Furthermore, it was pointed out that the vortical structures formed downstream of the tongue region, and their strength decreases as the distance from the tongue increases.
Other notable literature on the application of these methods in different turbomachines can be found in Refs. [7][8][9][10][11][12]. However, these methods can be challenging in the identification of the actual vortical rotational motion due to their dependence on the threshold value selected. In other words, inappropriate thresholds misidentify vortex structures which do not produce accurate information on the inner flow dynamics. To solve this problem, Liu et al. [13] recently put forward a method of vortex identification based on the decomposing of ∇ � V into a vortical part and non-vortical part and introduced a parameter called Ω. As stated in Eq. (2), Mathematically, Ω is given as: To avoid undefined cases and for practical purposes, a small value of ε = δ(b − a) max is added to the denominator to Eq. (6). With this criterion, Ω = 0 depicts no flow rotation while Ω → 1 reveals that the fluid does rigid body rotation. A standard threshold of Ω = 0.52 clearly identifies the vortex structures. The Ω-method has been rarely used to identify vortex structures in turbomachines. Zhang et al. [14][15][16] are among the few that have employed the Ω-criterion to capture the vortex structures in the complex flow fields in turbine machines. They concluded that the advances of the Ω-criterion in vortex identification of cavitation are well elaborated and analyzed. Furthermore, authors [1,[17][18][19] have also applied the Ω-criterion in other engineering applications. All vortical structures (weak and strong) were identified under a single standard threshold, Ω=0.52.
For purposes of this study, a side channel pump model with a convex blade is considered. Side channel pumps are becoming popular in the field of fluid machinery due to their ability to supply fluids at higher heads compared to centrifugal and axial pumps at the same tip speed [20,21]. Numerous studies, remarkably Authors [20,[22][23][24][25][26] have been conducted on the complex flow patterns and performance optimization of side channel pumps. Until now, the flow theory underlying the operation of side channel pumps remains a challenge due to the inability to fully identify and reveal the evolution mechanism of vortex structures in the turbulent flows. The Reynolds number of the investigated pump based on the mean velocity of the fluid, v 1 and inlet pipe diameter, D 1 under design conditions is determined by Re = ρv 1 D 1 /µ = 10 5 . Here, v 1 and D 1 represent the velocity and diameter of the inlet suction pipe respectively. This depicts and confirms the turbulent nature of the flow characterized in side channel pumps [4]. Therefore, this paper mainly presents comparative studies on the robustness and reliability of the newly proposed Ω-criterion against two traditional methods (Q and λ 2 ) to accurately capture the vortex structures in the highly turbulent flows linked to the investigated pump model. Besides, the pressure pulsation and entropy loss distributions characterizing the flows in the pump are comprehensively studied on different planes concerning the vortex structures evolved. The main contribution of this paper is the novel application of Ω-criterion to predict the local flow theory to enable the efficient design of side channel pumps.

Pump Model Description
Wei et al. [27] first demonstrated the necessity of the convex blade at the outer radius. The unsteady pressure pulsating within the pump significantly reduced with the outer radius convex blade. However, in this study, the side channel pump model replicated for this study is similar to that used by Wang et al. [4] however, a convex blade is introduced at the inner radius. The convex blade used in the present study has thickness, b 2 =2 mm and positioned at a height, h=7 mm. The pump operated at a rotational speed of 1500 r/min. The primary design specifications of the investigated side channel pump model are enlisted in Table 1. Figure 1 is a description of the studied side channel pump model showing the impeller, convex blade and side channel.

Turbulence Model Flow Solver
Numerical simulations were conducted using the commercial CFX software. The theoretical background of the flow motion is governed by the Navier-Stokes equations. The equations comprise the conservation of mass and momentum equations. Assuming the fluid used (water at 25ºC) is constant and incompressible, and then the equations are given as: where u i and u j represent the absolute and relative velocity components respectively. The pressure and density of the fluid are denoted by p and ρ , respectively. µ defines the kinematic viscosity measured in m 2 /s. −u ′ j u ′ i is the Reynold stress tensor and defined as: where k is the turbulent kinetic energy, and δ ij represent the Kronecker delta.
In the present study, the Navier-Stokes equations were closed by the shear stress transport (SST) k − ω turbulence model to solve the turbulence scales generated in the investigated pump [22,28]. This turbulence model was developed by Menter [29,30]. Owing to the robust ability of the SST k − ω model to combine the k − ω model in the wall boundary region and k − ε model in the main stream flow region. The formulations of the SST k − ω model are provided as: The coefficients of the models of the k − ω and k − ε models applied are provided in Table 2.
Total pressure of 1kPa was used as the inlet boundary condition. The flow direction was normal to the (14) boundary condition. Opening cylindrical velocity terms were specified at the outlet boundary condition for each operating condition. The results of the steady-state calculations failed to provide accurate flow characteristics and did not agree with the experimental measurements. Therefore, the impeller-side channel interface was updated to transient rotor-stator for the unsteady-state calculation with a time step, t = 1.111 × 10 −4 s . The total time for the unsteady calculations was 2.0 × 10 −1 s for five impeller revolutions. High-resolution advection and the second-order backward Euler transient schemes were recommended to solve the convection and diffusion terms. To ensure numerical accuracy, the residual targets of the continuity, momentum, k, ω and ε equations were set at 10 −4 .

Mesh Domains and Independence Study
As a requirement for numerical simulations, computational mesh cells were generated using ANSYS ICEM. From Refs. [31][32][33][34], the hexahedral structured meshing technique is applied to create cells for the impeller and side channel computational domains. The mesh of the convex blade and clearances (both axial and radial) are created together with the impeller (see Figure 2(b) and (c)). Similarly, the outlet pipe and side channel were meshed as a single domain (see Figure 2(d)) to avoid creating several interfaces which in turn reduces the numerical deviations. Mesh refinement was done at the boundary where the flow exchanges between the impeller and side channel.
To ensure that the simulated flow in the pump is independent of the number of nodes, five different sets of meshes were created with different number of nodes. Table 3 provides a mesh independence study. The head is determined using: Mesh 3 with a suitable topology and number of nodes was chosen and used for the numerical calculations due to the minimal effect of the number of nodes on the head of the pump. Since the mesh size is related to the Y+ value, the Y+ is an important quantity to determine the non-dimension distance from the wall to the first mesh  node. The averaged Y+ values of the impeller and side channel are 13 and 11 respectively. The maximum Y+ value is less than 65 for the entire pump, and it is acceptable for the industrial applications. The Y+ distribution on the impeller and side channel walls are presented in Figure 3 using the SST k − ω model. Most regions reveal Y+ value below 10 except for some portions of the side channel. Thus, the boundary details of all physical surfaces are set as to automatic non-slip wall function with smooth wall roughness.

Validation of Numerical Model
The reliability of the numerical simulations was verified by comparing them with the head observed in a closed experimental set-up. Figure 4 shows the head coefficient of the investigated side channel pump. The mean flow coefficient,ϕ and head coefficient,ψ were evaluated using the following expressions respectively.   For the measured operating conditions, the flow coefficient is in a range of φ=0.00249 to 0.016. Both results from the experiments and numerical approach show a common decreasing tendency from ψ = 0.0032774 to 0.001484. A good agreement is established with a maximum relative error of 12.3% at φ = 0.0128 over the experimental measurements. This deviation is partly due to the highly turbulent patterns of the flows and the losses such as mechanical and disk frictional losses that were not considered in the simulation. Therefore, the numerical calculation model used in this study is suitable to predict the complex vortex structures accompanying the flow in the investigated pump.

Determination of Appropriate δ Value for Ω-Criterion in Side Channel Pumps
It should be noted that the arithmetic average of the transient dataset obtained from the unsteady simulations are used to generate the results for detailed analyses. As pointed out earlier, a small positive value, ε is introduced to Eq. (6) to achieve the validity of the equation as presented in Eq. (19).
where ε = δ(b − a) max . According to Dong et al. [35], the value of ε is case dependent due to the computational noises associated with different flows. Importantly, ε is introduced only to remove the noise but not a threshold to vary the iso-surfaces of the vortex structures. Different values of δ in the range of 10 −9 ≤ δ ≤ 10 −2 are chosen to determine the influence of ε in Ω-criterion applications in side channel pumps. Figure 5 examines the dependency of δ on the value of ε on inlet, impeller and side channel domains of the pump. Here, the area here is the proportion of the iso-surfaces of the vortex structures formed in each domain. Generally, the flow characteristics in the impeller account for the largest number of vortex structures generated in the entire pump. When lgδ ≤ 10 −7 , there is minimal or no change in the value of the area for all domains thus no significant effect on the ε value is observed after δ=10 −7 . The deviation between the area of the impeller and side channel is 24.98%. The error defined by Eq. (20) at δ=10 −7 is the same for the impeller and side channel domains at 0.014%. The error recorded for the inlet pipe is 0.47%.
The iso-surfaces of the vortex structures formed under different δ values are presented in Figure 6. The strength of the vortex structures was colored by the velocity. Ambiguous vortex structures are revealed in the impeller passage when δ=10 −3 . This observation is counter to Figure 5 where the impeller records the largest area of the vortex structure iso-surfaces. On the other hand, no distinct differences can be recognized between the structures generated at δ=10 −5 and δ=10 −7 . Thus, from both qualitative and quantitative analyses, the small positive parameter is obtained using ε=10 −7 (b−a) max for the Ω-criterion in this case.

Comparison of Ω against Q and λ 2 Criteria
The operation of side channel pumps leads to various types of vortex structures in the internal flow field. Therefore, the iso-surfaces of the aforementioned vortex criteria are critically applied to accurately identify the vortex patterns in the entire investigated pump domains.  between the impeller and side channel [23]. Notably, this is an important observation since the strength of the exchange induces the several vortex structures produced. With Q-criterion, as the threshold value decreases (from Q= 1×10 −5 s −2 to Q= 5×10 −6 s −2 ), only weak vortex structures are identified. The λ 2 -criterion renders a similar tendency as the threshold value decreases from λ 2 = −5×10 −5 s −2 to λ 2 = −1×10 −4 s −2 . The tendency emanates from Eq. (5), such that when Q is positive, the eigenvalue λ 2 is usually, and a counter observation will be seen if Q is negative [36]. This makes it more challenging to select the suitable threshold using Q and λ 2 -criteria.
To avoid obscuring some of the vortex structures in the pump, it is essential to choose a standard thresholdbased criterion which is not dependent on the experience of the user. Both strong and weak vortex structures were captured by Ω-criterion when the threshold value was varied from Ω=0.52 to 0.64. No significant deviations were observed when the threshold changes because vortex is defined as a region where vorticity overtakes deformation. Thus, Ω> 0.52 perfectly captures all vortex structures in the entire pump domains. This confirms that amongst the three criteria, Ω demonstrates less sensitivity to threshold selection. The strength of the velocity distribution is common to the three criteria. From the earlier observations, Q and λ 2 -criteria have confirmed to be threshold-value dependent. Therefore, suitable thresholds of Q=1×10 −5 s −2 and λ 2 =−1×10 −5 s −2 were chosen to deeply reveal the vortex structures evolved in the impeller passage of the pump. The selection of these threshold values was guided by the results of Ω-criterion with a standard threshold of Ω= 0.52. From the inflow, similar trends are seen until the mid-section of the pump for the three criteria (see Figure 10). It is expected that water occupies the impeller during the flow exchange between the impeller and side channel. For example, at interrupter where there is no channel the impeller passage is fully occupied as identified by the three criteria. The Q-criterion depicts larger empty areas at the inner radius of the blade around the mid-section and closer to the outflow as labeled as X1 and X2 at the same region for the three criteria. A common observation is realized with Ω-criterion but λ 2 reveals otherwise. The inner radius is occupied with the flow at mid-section under λ 2 -criteria. It can be assumed that the λ 2 -criteria misidentifies the vortex structures at that point (X2) compared to Ω-criterion. Qualitatively, the iso-surfaces produced from the three criteria in the impeller are in close accordance. However, the orientation of some of the vortex structures illustrates some discrepancies therefore the vortex area under each criterion is plotted as supplied in Figure 11. Guided by the physical definition of Ω-criterion, the area under Ω-criterion is the smallest while Q and λ 2 -criteria recorded the highest. The area covered under each criterion is expressed as Ω, Q and λ 2 for the values 0.06262, 0.06278 and 0.06262 respectively, using the thresholds stated in the paragraph above. The Ω-criterion proves to be reliable and avoids subjective errors during the selection of the threshold.

Vortex Structures Intensity under Ω-criterion
For detailed insight into the vortex evolution and mechanism, different planes (see Figure 12) are set from the inflow to outflow positions in different directions. Planes I to VI, planes VII to IX and planes X to XII are set in the longitudinal, axial and radial directions respectively. Following previous works [23,26,37], a high-performance side channel pump ought to be characterized by strong and orderly organized longitudinal vortices to reduce the flow losses in the impeller passage caused by axial and radial vortices restricted in the impeller. The momentum transfer or exchange of energy of the flow from the impeller to the side channel is responsible for the positive operation (longitudinal vortex)of these kinds of pumps. Typically, the flow is expected to be guided by the outer radius of the impeller into the side channel and re-enter at the inner radius as portrayed (see Figure 12(c)). Side channel pumps mainly dominate three vortex structure types (longitudinal, axial and radial).
The Ω-criterion is non-dimensional and measures from 0 to 1. According to Liu et al. [13], a vortex is evolved in fluid flow when the threshold is greater than 0.52. Therefore, values of Ω greater than 0.8 identify the evolution mechanism and development of a vortex core. In general, the higher the intensity of longitudinal vortex, the more energy the impeller transfers to the water and more work is done. Figure 13 shows the vortex intensity distribution coupled with the velocity vector and streamline plots in the longitudinal direction. There is a number of high vortex intensity regions discovered in the impeller, side channel and as well as the interface. The regions highlighted in black represent the symmetrical longitudinal vortex developed close to the interface. The longitudinal vortex structures promote flow exchange between the impeller and side channel. Furthermore, the axial and radial vortex structures that account for the high vortex intensities are restricted in the impeller passage. These regions are highlighted in red and considered to have negative influence by restraining the longitudinal vortex mechanism. These are in good relation with the convergence of the streamlines.
On plane I which is the interrupter, all the vortex generated are restricted to the impeller due to the absence of the channel. These negatively influenced vortex structures are likely to be transferred to the inflow section of the pump causing undesirable flow patterns. On plane II, the water is received from the inlet pipe, and the flow begins to develop the longitudinal vortex structures as observed in the velocity vector plots due to the weak centrifugal force received. Several vortex structures recognized on this plane are of high intensities greater than 0.65 and constrained in the impeller. This leads to weak longitudinal vortex with interspersed distortions.
As the flow approaches plane III, strong and weak longitudinal vortex structures evolve close to the interface. A vortex core is formed and promotes the flow exchange however some vortex structures restricted close to the impeller impede the longitudinal flow exchange. The velocity vectors and streamlines are in close accordance with the vortex plot. The velocity vectors show uniform and expected flow patterns at plane IV (mid-section). At the mid-section, the flow possesses a significant amount of centrifugal force to enhance the exchange, but the sudden change in geometrical structures restricts some vortex structures in the impeller. The exchange enhances on plane V with some portions of weak vortex structures in the impeller. Incomplete longitudinal vortex structures are generated close to the interface. At the outflow section VI, the flow leaves the pump with high intensity accounting for the high head built at the outlet pipe.
It is worth noticing that uncharacteristic flow patterns are highlighted in purple on the velocity vector plot on plane VI; instead of the flow to exit at the outer radius, some portions of the flow enter at the outer radius. For plane II, the flow leaves the impeller at the inner radius instead of the outer radius. These uncharacteristic phenomena lead to violent collisions resulting to a vortex restricted in the impeller. From the above discussions, the intensity of the longitudinal vortex in the longitudinal direction is a little higher but some vortex structures restricted in the impeller relatively generate certain energy losses which affect the generation of longitudinal vortex between the impeller and side channel. The results supplied by the Ω-criterion confirm the physical meaning of vortex as proposed by Liu et al. [13]. As discussed earlier, the vortex structures centered in the impeller are answerable to the low efficiency of side channel pumps because they retard the beneficial effects of the longitudinal vortex. To define the axial vortex in  Figure 12 is close to the impeller and side channel interface; Plane VIII is set at the middle of the impeller while Plane IX is set at the bottom of the impeller. The velocity vectors and streamlines in the axial planes are also studied per the axial vortex intensity distribution (see Figure 14). The intensity or strength of the axial vortex generally decreases downwards the height of the impeller. On plane VII, high vortex intensities of about Ω= 0.6 are identified at the outer radius of each blade pitch except the interrupter. From the velocity vector, secondary and reverse flows are established at the outer radius confirming the results in the axial vortex intensity plots. Sets of axial vortex structures evolve and scatter around the central and inner radiuses of the blade.
On plane VIII, the intensity of the axial vortex formed increases and widely scatters around the inner, central and outer radiuses of the blade. The central radius from visual observation show strong intensity distribution of axial vortex structures. In the middle of the impeller, the velocity vector and streamline plots show that the flow does not gain its full flow characteristics as it moves towards the interface for the exchange with the side channel. Possibly, this could be responsible for the high intensity of axial vortex structures generated. Weak and small amounts of axial vortex structures evolve at the bottom of the impeller close to the pump casing. Noticeably, an axial vortex core develops at central radius on plane VIII. The core increases in size and strength at plane VIII and diminishes on plane IX. Again, the axial vortex structures are concentrated on the pressure side of the blade on all the axial planes. Figure 15 shows the vortex intensity together with the velocity vectors and streamlines on the radial planes at the mid-section (about 150° offset from inflow) of the impeller. The radial vortex structures are localized in the impeller by the sides of the blades thereby opposing the momentum exchange developed by the positive

Energy Loss Distribution
For side channel pumps and other hydraulic processes, irreversible losses are associated with their operations. The mechanical energy is always transformed into internal energy owing to the viscosity and highly turbulent flow characteristics. Furthermore, shear can also cause energy dissipation and the energy dissipation partly depends on the negative influences of some vortex structures generated. Therefore, to study the relationship between the energy dissipation and the vortex structures identified by Ω-criterion, the local entropy production method [38,39] provides an effective approach to investigate the negative effects emanating from the highly unsteady flow characteristics in the investigated pump. The turbulent dissipation accounts for the main losses in side channel pumps compared with direct dissipation. Moreover, the impeller generates the highest entropy loss compared to the inlet pipe, side channel and outlet pipe. The mathematical formulations of the entropy production method and their applications in side channel pumps are well explained in Refs. [24,31] . The coefficient of entropy is defined as: where D ′ represents the averaged turbulent dissipation, U 2 is the circumferential velocity of impeller outlet, V is volumetric flow rate, and d 2 is impeller outer diameter. As said earlier, the flow enters the pump from the inlet in the circumferential direction on plane II. As the impeller rotates, the flow at the upper part of the inner radius of blade leaks to enter the clearance and rotates along with the impeller. The flow then rotates in the radial direction with the centrifugal force gained. Besides, as the rotational kinetic energy of the impeller on the flow is transformed into pressure head at the outlet, there is backflow from the outflow to the inflow which does not aid in the energy conversion process. From the velocity vector plot, vortex structures along the flow direction are developed in the side channel which induces the flow at the upperpart of the inner radius of the blade. The intensity of the vortex increases to form a vortex core as indicated on the longitudinal vortex plot. This creates high entropy loss, ξ D ′ ≥ 11.0 at the upper part of the inner radius of the blade as shown in Figure 16. Moreover, the collision as realized on the velocity vector plot in Figure 13 on plane II causes a large vortex which is not characteristic of the pump. Thus, the boundary of the side channel generates high energy loss, ξ D ′ ≥ 11.0 . The (21) Figure 15 Distributions of vortex intensity and velocity vector plots on radial planes energy loss gradually reduces from plane III in both side channel and impeller. On planes IV and V, the channel records low entropy coefficient of less 7.0 but the losses at the upper part of the inner radius is still evident. Generally, losses produced in the longitudinal plane are mainly located at the inner radius of the blade which affects the output efficiency of side channel pumps. Particularly, the losses located in the impeller are consistent with the vortex structures restricted in the impeller. Planes I to VI show good agreements with the longitudinal vortex plot.
On the other hand, the outer radius generates a low loss, ξ D ′ ≤ 7.0. Figure 17 presents the entropy loss distribution in the axial planes as referenced in Figure 12. The losses generated in the axial planes are concentrated at the inner radius of both suction and pressure sides of the blade. As registered on the velocity vector and streamline plots in the axial plane, rotating vortex structures are formed as the flow enters the channel predominately at the inner radius. This weakens the energy output by the impeller to the flow making it distract the flow exchange. Thus, the losses close to the impeller-side channel interface (plane VII) are relatively higher than plane VIII. At the bottom of the impeller, there are high flow leakages that occur at the axial clearance thus efficient longitudinal vortex cannot develop fully to promote the energy transfer process. Therefore, plane IX records the highest ξ D ′ at the inner radius and extends to the central radius as well.
The flow enters the impeller from the inner radius and as shown in the velocity plots in the longitudinal direction, portions of the flow leave at the inner radius causing high vortex cores. The cores partly contribute to the high loss that occurs at the inner radius. Plane X shows the highest entropy loss as shown in Figure 18. As the flow travels to the outer radius of the blade along the radial direction, the degree of the distribution decreases as the plane is moved from XI to XII. The vortex structures identified in the radial direction show a larger deviation with the entropy loss distribution. The clearance at the bottom of the impeller contributes maximum entropy loss especially on planes X and XI.

Pressure Pulsation Distribution
As suggested by Zhang et al. [26], the dynamic pressure pulsations in side channel pumps are mainly influenced by the flow exchange between the impeller and side channel. Similarly, the flow exchange brings about the vortex evolution and considering the complex flow characteristics in the investigated pump, and it is important to analyze the pressure pulsations in different planes to the vortex structures identified. The pressure pulsation is monitored for 1800 cycles and the results of the last 24 cycles are examined in this study. The pressure pulsation coefficient is given by the root-mean-square method as: where N, p i and p denote the sample number, the pressure at each time-step, and arithematic average of pressure respectively.
Generally speaking, Figure 19 reveals that the contour of pressure pulsation coefficient; C ′ p in the impeller is higher than the side channel on the longitudinal plane. The impeller reported C ′ p > 0.7 whiles the side channel also reported C ′ p ≤ 0.7 . With the development of the vortex structures confined in the impeller, it is found that high C ′ p are obvious at the upperparts and roots of the inner radius of the blade on planes I to VI. The intensity of the pressure pulsation at the root of the impeller at the inner radius even reaches the maximum value at all planes in the longitudinal direction. Therefore, the vortex cores developed in the impeller generates significant pulsations. Figure 20 shows that the distribution of the pressure pulsations intensity on the axial planes in the impeller is high and irregular due to flow turbulence from the vortex structures evolved. Most areas on plane VII reveal C ′ p > 0.7 . The average C ′ p decreases from planes VII to plane VIII but increases on plane IX. Plane IX records maximum C ′ p values at the inner radius of the blade from the inflow to outflow. This could be attributed to the vortex structures generated close to the bottom of the impeller which is partly due to the flow leakage into the axial clearance. By and large, the outer radius showed high C  Figure 18 Distributions of entropy loss on radial planes