Research Article  Open Access
Study on Flow Characteristics in Volute of Centrifugal Pump Based on Dynamic Mode Decomposition
Abstract
To investigate the unsteady flow characteristics and their influence mechanism in the volute of centrifugal pump, the Reynolds timeaveraged NS equation, RNG kε turbulence model, and structured grid technique are used to numerically analyze the transient flowfield characteristics inside the centrifugal pump volute. Based on the quantified parameters of flow field in the volute of centrifugal pump, the velocity mode contours and oscillation characteristics of the midspan section of the volute of centrifugal pump are obtained by dynamic mode decomposition (DMD) for the nominal and low flowrate condition. The research shows that the firstorder average flow mode extracted by DMD is the dominant flow structure in the flow field of the volute. The secondorder and thirdorder modes are the most important oscillation modes causing unsteady flow in the volute, and the characteristic frequency of the two modes is consistent with the blade passing frequency and the 2x blade passing frequency obtained by the fast Fourier transform (FFT). By reconstructing the internal flow field of the volute with the blade passing frequency for the nominal flowrate condition, the periodic variation of the unsteady flow structure in the volute under this frequency is visually reproduced, which provides some ideas for the study of the unsteady structure in the internal flow field of centrifugal pumps.
1. Introduction
The internal flow of the centrifugal pump is a complex threedimensional unsteady turbulent flow, often accompanied by flow separation, cavitation, hydraulic vibration, etc., which affect the stability of flow in the centrifugal pump. Among them, the rotorstator interaction will produce periodic pressure fluctuation, which will lead to the intensification of pump vibration, noise enhancement, and performance decay, seriously affecting the safe and stable operation of the unit. Especially in the tongue area, since the fluid flow from the impeller outlet impacts the volute tongue, the rotorstator interaction at this place should be the strongest, which is one of the key areas for the study of unsteady flow in centrifugal pumps. For instance, Kelder [1] studied the flow in the volute of a low specificspeed pump both experimentally and numerically near its design point and found that near the design point of this pump, the core flow behaves like a potential flow, provided that no boundary layer separation occurs. González [2] captured the unsteady flow effects inside a centrifugal pump due to the impellervolute interaction with a numerical simulation and provided the modeling possibilities for the prediction of the dynamic interaction between the flow at the impeller exit and the volute tongue. Keller [3] investigated the unsteady flow in a centrifugal pump with vaneless volute based on PIV measurements. The obtained data show that the fluiddynamic bladetongue interaction is dominated by highvorticity sheets (positive and negative) being shed from the impeller channels, especially from the blade trailing edges, and their impingement on the tongue tip with subsequent. Tan [4] elucidated the detailed flow field and cavitation effect in the centrifugal pump volute at partial load condition. The results showed that the maximum amplitudes of pressure fluctuations in volute at serious cavitation condition are twice that at noncavitation condition because of the violent disturbances caused by cavitation shedding and explosion. Hamed [5] investigated reducing the radial force, the effects of concentric volute, and multivolute geometry on the head, efficiency, and radial force of a low speed centrifugal pump at offdesign conditions. The results showed that the triplevolute is the most appropriate volute geometry at offdesign conditions. Wang [6] carried out the transient analysis to investigate the dynamic stress and vibration of volute casing for a large doublesuction centrifugal pump by using the transient fluidstructure interaction theory. The results revealed that, for all operating conditions, the maximum stress located at the volute tongue region, whereas the maximum vibration displacement happened close to the shaft hole region, and the blade passing frequency and its harmonics were dominant in the variations of dynamic stress and vibration displacement. Liu [7] investigated the pressure fluctuation intensity and vortex characteristic of a mixedflow pump as turbine at pump mode with a tip clearance. Hao [8] studied the transient cavitating flows of a mixedflow PAT (pump as turbine) at pump mode experimentally and numerically. Tan [9, 10] investigated the role of blade rotational angle in the energy performance and pressure fluctuation and the influence of original and Tshape blade end on performance of a mixedflow pump through an experimental measurement and numerical simulation, respectively.
Dynamic mode decomposition (DMD) is a datadriven algorithm for extracting dynamic information from unsteady experimental measurements or numerical simulations. For such simulations and experiments provide largescale data, it is necessary to understand essential phenomena from the data provided. So it can be used to analyze the main characteristics of complex unsteady flows or to establish loworder flowfield dynamics models. The DMD method was first proposed by Schmid [11–13] in 2008, and then the flow field data obtained by numerical simulation and experiment were analyzed by DMD method. Seena [14] employed the dynamic mode decomposition to analyze the unsteadiness in extracted modes without the explicit knowledge of evolution operator of the data. Hemati [15] presented two algorithms and compared the application of the two algorithms on cylinder wake data collected from both direct numerical simulations and PIV experiments. Tu [16] demonstrated the utility of this approach by presenting novel sampling strategies that increase computational efficiency and mitigated the effects of noise, respectively. Pan [17] used a linear combination of a sinusoidal unstable wave and its highorder harmonics as a prototype to take an error analysis of DMD algorithm. The result shows that the superimposition of finer structures with less energy dominance might damage the estimation accuracy of the primary structures’ growth rate. Diana [18] proposed a new vectorfiltering criterion for dynamic modes selection that is able to extract dynamically relevant flow features from dynamic mode decomposition of timeresolved experimental or numerical data. Kunihiko [19] described herein some of the dominant techniques for accomplishing these modal decompositions and analyses that have seen a surge of activity in recent decades and presented a brief overview of several of the wellestablished techniques and clearly lay the framework of these methods using familiar linear algebra. Liu [20, 21] employed the dynamic mode decomposition (DMD) and Proper Orthogonal Decomposition (POD) to analyze the coherent structure of cavitating flow around ALE 15 hydrofoil and investigated the complicated transient characteristics of gasliquid twophase flow in a multiphase pump under 10% and 20% inlet GVFs by dynamic mode decomposition (DMD).
2. Numerical Model and Computational Method
2.1. Geometry of the Pump
The object of this paper is a conventional singlestage singlesuction centrifugal pump. The main parameters of the pump are listed in Table 1. The wrap angle α in Table 1 is defined as the angle between the line of the blade inlet edge and the center of the impeller and the line of the blade outlet edge and the center of the impeller. The threedimensional water body (impeller and volute) required for the calculation of the model is generated by threedimensional design software Pro/E 5.0. According to the design parameters, the number of blades is 6, and the volute is a semihelical double volute. It should be pointed out that compared with a single volute, the presence of the baffle in the double volute divides the volute into two flow paths, and the volute tongue and the baffle entrance are symmetrical about the center of rotation. Highspeed fluid from impeller passages still causes unsteady flow in volute when it impacts the inlet of baffle. Therefore, the existence of baffle is equivalent to the presence of a second tongue in double volute. So, this area also needs important attention in the study of internal flow in volute.

2.2. Comparison of Numerical Simulation and Experiment in the Centrifugal Pump
For the entire computational domain (composed of suction pipe, impeller, volute, and outlet pipe), structured hexahedral mesh is generated using AnsysICEM mesh generation tool. The overall calculation domain is shown in Figure 1. After the grid independence test, the total number of grids was 4.5 million. To calculate the unsteady flow structures in the centrifugal pump, the commercial CFD code AnsysFluent is used in the present paper. The RNG kε twoequation model is used under the condition of satisfying the accuracy. The discrete method of the governing equation is the finite volume method. The coupling between the pressure and velocity is solved using the simple scheme, and secondorder upwind scheme is used for spatial discretization. The uniform velocity at the pump inlet is set as the mass flow inlet boundary condition, and the outflow at the pump outlet is set as the outlet boundary condition. The fluid medium is the normal temperature water; and the boundary condition is the noslip solid wall. To guarantee the accuracy of the result, the convergence precision of the continuity residual, xvelocity residual, yvelocity residual, zvelocity residual, k residual, and epsilon residual for steady simulation is set as 10^{−5}.
As shown in Figure 2, the performance test of the centrifugal pump was tested on a closed test stand. And Figure 3 shows the head/flowrate relationship and the efficiency/flowrate relationship curves obtained from both numerical simulation and experiment of the centrifugal pump. The results show that as the flow rate increases gradually, the head of the centrifugal pump gradually decreases, and the efficiency first increases and then decreases. For the different flowrate conditions at the centrifugal pump, both the numerically calculated head and efficiency values are higher than the experimental values, and the errors for the nominal flowrate condition are 5.8% and 3.7%, respectively. The errors between numerical simulation and experiment are due to the fact that the mechanical losses caused by bearings, mechanical seals, etc. are often ignored in numerical calculation. From the comparison of numerical simulation and experimental values of head and efficiency of the centrifugal pump, the numerical simulation method was used to predict the performance of the centrifugal pump having good accuracy of the reference value. In addition, it can be seen from Figure 3 that in the efficiency/flowrate curve and head/flowrate curve, the maximum difference appears at the nominal flowrate condition, but the values obtained by numerical calculation are in good agreement with the experimental values for the low flowrate condition. Therefore, it can be explained that the centrifugal pump is in a state of offdesign condition when it works. The main reason for this phenomenon is that the design parameters of impeller are unreasonable, such as the impeller outlet diameter and the impeller outlet width being too small.
3. Dynamic Mode Decomposition
Dynamic mode decomposition is a method for extracting the mode in a flow based on flowfield snapshots, so that the flow structure can be accurately described. For the linear flow, the DMD method can extract the modes that can characterize the global flow stability. For the nonlinear flow, the DMD method can describe the flow structure in which the observations (such as velocity and pressure) dominate. The mathematical derivation process of it is introduced below.
Suppose there is a set of observation data matrices that vary with time:
where N is the total number of snapshots of the flow field, and the column vector v_{i} is the data of the i th snapshot. It is assumed that the snapshot data of the flow field at two adjacent moments can be represented by a linear transformation matrix A:
As the acquired snapshot data increases, we can further assume that the vector formed by the snapshot data eventually tends to be linearized. Thus, the last flowfield snapshot can be represented as a linear combination of all previous snapshots:The matrix form isIn formula (4), ; r is residual vector.
From formula (1), The new relation can be obtained by combining formula (3) as follows: In formula (6),Obviously, when the residual vector r is small, the eigenvalue of the matrix S is approximated by the eigenvalue of the system matrix A. Therefore, the matrix S is a lowdimensional approximation of the system matrix A, and its eigenvalues can represent the main eigenvalues of the system matrix A. Matrix S is usually obtained by QR decomposition of matrix ; if matrix can not guarantee full rank, its QR decomposition is not unique and may lead to the irreversibility of upper triangular matrix R, which will ultimately lead to the impossibility of finding matrix S. Perform singular value decomposition on matrix :
From formula (5) and formula (7):
From the previous analysis, we can see that the eigenvalues of matrix can represent the main eigenvalues of system matrix A. Then the eigenvalues and eigenvectors of matrix are obtained:
Finally, the DMD modes can be obtained as follows:
It should be pointed out that the eigenvalue μ_{i} contains the information of the mode Φ_{i}. When the eigenvalues are expressed in the complex plane, the modes on the unit circle of the complex plane are relatively steady, while the modes whose eigenvalues are not on the unit circle are unsteady.
The corresponding frequency ω_{i} and growth rate of the modes are defined as follows.According to the similar matrix obtained by the above DMD decomposition, the evolution of the flow field can be further estimated. By singular value decomposition (7), a highdimensional system v_{i} can be mapped to a subspace z_{i}.
By simple transformation there areLet be the column vector of matrix W; N has the eigenvalues of the matrix in its diagonal. Then the eigendecomposition can be expressed as
So, from the previous derivation, a snapshot at any time instant i can be approximated as
From the DMD mode definition formula (10),The mode amplitude α is represented aswhere α_{i} denotes the amplitude of the ith mode, which represents the mode contribution to the initial snapshotv_{1}.
Substituting formula (17) and (18) into formula (16), the flow field at any time instant can be predicted as
Then its snapshot sequences can be expressed asThe dynamic mode decomposition of the velocity field in the midspan section of the centrifugal pump volute can be performed by using the DMD calculation formula derived above. Firstly, we need to edit the formula of DMD calculation in MATLAB and then export the coordinates of grid points in the midspan section and the velocity values at the grid points in CFD Post. From the definition above, it can be seen that the set of velocity values at different grid points at the same time is a column vector of the whole data matrix. The final calculation data matrix is made up of the set of velocity values at different times. Next, the dynamic mode decomposition of the velocity matrix is then performed using MATLAB. Finally, the calculated modal values are edited into the acceptable data format of CFD Post and imported into CFD Post to make the corresponding velocity contours.
4. Results and Discussions
4.1. Unsteady Pressure Fluctuation Analysis inside the Volute
In the unsteady calculation, the spectrum analysis is performed to obtain the characteristic frequency of the unsteady flow inside the centrifugal pump. The dimensionless pressure coefficient C_{p} is introduced to describe the pressure fluctuation characteristics of each monitoring point. The expression is as follows:
where P_{i} is the static pressure value of the monitoring point at a certain time, ; P_{ave} is the average value of static pressure in one cycle, ; u_{2} is the circumferential velocity at the impeller outlet, m/s. The monitoring points P_{1}, P_{2}, P_{3}, and P_{4} are arranged inside the volute as shown in Figure 4. P_{1} and P_{3} are set to monitor pressure fluctuation characteristics of fluid flow in volute tongue and baffle entrance region; P_{2} is set to monitor the pressure fluctuation characteristics of fluid flow at the volute inlet except for the tongue and the baffle entrance region; and P4 is set to monitor the pressure fluctuation characteristics at the volute outlet. By comparing the pressure fluctuation characteristics at these four points, the main factors causing unsteady flow in volute are explored.
Figure 5 shows the spectrum analysis results of the pressure fluctuation at the monitoring points inside the centrifugal pump for the low flowrate (0.5) condition. Obviously, the main frequency of the pressure fluctuation in the centrifugal pump is the blade passing frequency (f_{1}=171Hz), and the pressure amplitude at 2x blade passing frequency (2f_{1}=342HZ) is also prominent. Figure 6 shows the spectrum analysis results of the pressure fluctuation at the monitoring points inside the centrifugal pump under the low flowrate (0.5) condition. For the nominal flowrate condition, the blade passing frequency of the centrifugal pump is f_{1}=169.59Hz and the 2x blade passing frequency is 2f_{1}=339.17HZ. Comparing the results of spectrum analysis for the two operating conditions, it can be seen that the dominant pressure fluctuation main frequency and the high amplitude pressure fluctuation frequency are basically consistent. However, compared with the nominal flowrate condition, it produces a certain low frequency fluctuation at less than the blade passing frequency in the low flowrate condition.
4.2. Modal Analysis inside the Volute at the Low FlowRate Condition
Based on the previous CFD unsteady calculation results, DMD analysis was performed on 120 flowfield speed snapshots with time interval Δt=0.0015s, and the DMD results for the two operating conditions were obtained. Figure 7 shows the distribution of the DMD mode eigenvalue in the complex plane of the unsteady flow field in the midspan section of the volute of the centrifugal pump at low flowrate (0.5) condition. In this figure, the real part of the eigenvalue is the horizontal axis and the imaginary part is the vertical axis. All the eigenvalues are distributed near the unit circle, and some eigenvalues are distributed on the unit circle, which means the corresponding modes are neutral stability. Figure 8 shows the relationship between the corresponding frequency of the DMD modes and the correlation coefficient. The correlation coefficient is used to measure the influence of each mode on the original flow field. In [22], C_{j} is defined as the correlation coefficient.where α_{j} is the modal amplitude and μ_{j} is the eigenvalue and is the mode Frobenius norm. After sorting the DMD modes according to the correlation coefficient, the first six modes with the highest correlation coefficient can be clearly seen from Figure 8. Since the corresponding eigenvalue of the mode with the largest correlation coefficient is a real number, its frequency is zero. Other eigenvalues are complex numbers, and conjugate pairs appear, which can be observed from the distribution of eigenvalues. Since the parameters describing the flowfield information in each mode are their corresponding real parts, modes corresponding to a pair of conjugate eigenvalues are the same mode. The frequencies of the second and thirdorder modes are consistent with the blade passing frequency and 2x blade passing frequency of the centrifugal pump obtained by the FFT in Figure 5, which fully demonstrates that the pressure fluctuation frequency in the original flow field is objectively present, and the capture of the characteristic frequency by the DMD is very accurate. From Figures 7 and 8, it can be clearly seen that the highenergy mode eigenvalues are distributed on the unit circle, so the corresponding flow is relatively steady, which belongs to the main flow structure in the flow field, while the flow corresponding to the lowenergy mode eigenvalues which are not distributed on the unit circle is unsteady, and it is not the main flow structure in the flow field.
Figure 9 is the first sixthorder mode velocity contour of the midspan section of the volute for the low flowrate condition (0.5). Figure 9(a) is the firstorder zerofrequency velocity mode, and the eigenvalue corresponding to this mode is real number, so the frequency of this mode is zero; that is, the flow structure is timeaverage and is the average flow mode of the velocity field, which shows the dominant flow structure in the original flow field. Actually, the average flow mode can be regarded as the basic structure of the velocity field inside the volute, and the original flow field inside the volute can be regarded as being formed by superimposing oscillation modes of different frequencies on this basic structure. It can be seen from the velocity contour that due to the existence of the volute tongue and the volute baffle, two highspeed fluid regions are distributed at the entrance of the tongue and the baffle, which affects the circumferential distribution of the speed in the volute inlet. Due to the existence of the baffle, the velocity field of the basic flow structure inside the volute is symmetrically distributed with respect to the center of rotation, which also changes the flow state inside the volute to some extent.
(a) Firstorder mode
(b) Secondorder mode
(c) Thirdorder mode
(d) Fourthorder mode
(e) Fifthorder mode
(f) Sixthorder mode
Figures 9(b) and 9(c) are velocity contour of the secondorder and thirdorder modes of the midspan section of the volute for the low flowrate condition (0.5). From the comparison of correlation coefficients, the second and thirdorder velocity modes are the two most important modes besides the firstorder zerofrequency modes. Therefore, the twoorder velocity modes are also the main oscillating modes causing unsteady flow in the volute. From the velocity contour in the second and thirdorder modes, it can be seen that there are 6 and 12 periodic high and low velocity clusters in the circumference of the volute inlet, respectively. However, due to the influence of the tongue and the baffle, the high and low velocity fluid regions alternately distributed in the two places are weakened to a certain extent. The frequencies of the twoorder velocity modes correspond exactly to the blade passing frequency and 2x blade passing frequency obtained in the FFT, which indicates that the DMD captures the influence of the centrifugal pump rotorstator interaction flow structure. Figures 9(d), 9(e), and 9(f) are the velocity contour of the fourth, fifth, and sixth modes of the midspan section of the volute, respectively. The fourthorder mode captures the unstable flow structure caused by three high velocity fluid clusters downstream of the baffle inlet. It can be seen from the fifth order mode velocity contour that high and low velocity fluid clusters alternately distribute at the inlet of the tongue and the baffle. And high and low velocity fluid clusters alternately distribute in a direction perpendicular to the line connecting the tongue and the partition in the sixth order mode velocity contour. The distribution indicates that there are vortex structures with opposite rotation directions, reflecting the unsteady characteristics of highorder flow inside the volute.
4.3. Modal Analysis inside the Volute at the Nominal FlowRate Condition
Figures 10 and 11 show the distribution of the DMD mode eigenvalues and correlation coefficients of the unsteady flow field in the midspan section of the volute of the centrifugal pump for the nominal flowrate condition (). In the distribution of eigenvalues, it can be seen that the eigenvalues of the two modes are in the neutral stable range of the unit circle, but because of their low correlation coefficients, the influence on the original flow field is weak.
Figure 12 is the first sixorder velocity mode contour of the midspan section of the volute for the nominal flowrate condition (), where Figure 12(a) is still the firstorder zerofrequency mode, which is also the basic structure of the velocity field inside the volute. The velocity contour of the firstorder mode is similar to the zerofrequency velocity contour in the low flowrate condition. A highspeed fluid region is evenly distributed in the circumferential direction at the volute inlet, and the flow of fluid is still symmetrical about the center of the axis of rotation in the two channels of the volute. The basic flow structure in the volute is still affected by the presence of the tongue and the baffle, but the degree of influence is greatly reduced compared to the low flowrate condition.
(a) Firstorder mode
(b) Secondorder mode
(c) Thirdorder mode
(d) Fourth order mode
(e) Fifth order mode
(f) Sixth order mode
Figures 12(b) and 12(c) are the second and thirdorder velocity mode contour of the midspan section of the volute for the nominal flowrate condition (), which are similar to the corresponding modes for the low flowrate condition. The two modes are still the main oscillation mode for causing unsteady flow in the volute. Similarly, the frequencies of the two modes are the same as the blade passing frequency and 2x blade passing frequency of the centrifugal pump, which indicates that whether it is in a low or nominal flowrate condition, the rotorstator interaction of the impeller and the volute is the main cause of the unsteady flow in the volute of the centrifugal pump. However, unlike the low flowrate condition, the periodic velocity clusters in the two velocity modes are almost not affected by the tongue and baffle for the nominal flowrate condition.
Figures 12(d), 12(e), and 12(f) are the fourth, fifth, and sixthorder velocity mode contour of the midspan section of the volute for the nominal flowrate condition (). In these order velocity mode contours, there are 8, 10, and 7 periodic velocity clusters in the inlet of the volute, respectively. It can be seen from the velocity contours that the influence range of the high and low velocity fluid clusters in these three modes on the flow field is much smaller than that in the second and thirdorder modes in the volute, which is also consistent with the distribution of correlation coefficients. In addition, compared with similar modes for low flowrate condition, the distribution of unsteady fluid clusters is still not affected by the tongue and baffle in the latter three modes for the nominal flowrate condition. Combining with the second and thirdorder modes, it shows that the existence of the tongue and baffle can restrain the unsteady flow structure inside the volute when the flow rate of the centrifugal pump decreases. It should also be pointed out that except the second and thirdorder modes which cause unsteady flow in the volute, the other three modes are all distributed between the blade passing frequency and 2x blade passing frequency for the nominal flowrate condition, while the frequency of the fourthorder mode with a high correlation coefficient is less than the blade passing frequency for the low flowrate condition (0.5), which indicates that the unsteady flow in the volute is more prone to occur at low frequencies when the flow rate decreases.
4.4. Reconstruction of Flow Field in Volute Based on DMD
In order to further observe the effect of dynamic mode decomposition on the extraction of flowfield characteristics in volute of centrifugal pump, a reduced order model of unsteady flow field in volute was established based on formula (19), and the flow field was reconstructed by using the obtained dynamic mode decomposition method. The unsteady flow field is reconstructed by the first ten modes with the highest correlation coefficient obtained by DMD. Figures 13 and 14 show the comparison of the reconstructed velocity contour and the original flow field velocity contour at T/2 for the low flowrate condition and the nominal flowrate condition. From the comparison, it can be seen that the reconstructed results have a high degree of reduction and identification for the flow structure in the flow field.
(a) Original flow field at T/2
(b) Reconstructed flow field at T/2
(a) Original flow field at T/2
(b) Reconstructed flow field at T/2
Since the firstorder mode is an average flow mode, it does not change during the entire rotating period of the centrifugal pump. In order to further study the unsteady flow structure in a certain mode, the mode can be superimposed on the firstorder average flow mode to observe the oscillation law of the mode. The secondorder mode caused by rotorstator interaction is superimposed on the average flow mode to reconstruct the flow field of a single mode, so as to observe the variation of the unsteady structure with time in the mode.
Figure 15 shows the reconstructed velocity contour of the mode corresponding to the blade passing frequency at different times. It can be seen from the velocity contour that six highspeed fluid clusters are uniformly distributed in the circumferential direction of the volute inlet, and at these different times, the six highspeed fluid clusters always occupy a dominant position. The motion process of the six highspeed fluid clusters is consistent with the rotary scanning process of the blade, which indicates that the unsteady flow structure in the volute caused by rotorstator interaction is not fixed but is synchronized with the rotation of the blade.
5. Conclusions
In this paper, the internal flow of the centrifugal pump is numerically simulated, and the velocity field of the midspan section of the centrifugal pump volute is analyzed by dynamic mode decomposition for the low flowrate condition and the nominal flowrate condition. The characteristic frequency of the flow field and the flow mode of different frequencies are extracted, and the reduced order analysis of the original flow field is realized. In the extracted flow modes of different frequencies, the firstorder mode is the flow structures occupying the dominant position of the original flow field for the low flowrate condition and the nominal flowrate condition. The second and thirdorder modes are the main oscillation modes of the original flow field, and the two modes’ characteristic frequency is consistent with the blade passing frequency and 2x blade passing frequency obtained by FFT, which proves that the rotorstator interaction between the impeller and the volute is the main reason for the unsteady flow in the volute. Next, comparing the mode velocity contours of each order for the condition of low flow rate and nominal flow rate, when the flow rate of the centrifugal pump is reduced, the existence of the tongue and the diaphragm will have a certain inhibitory effect on the unstable flow structure inside the volute, but it is more prone to occur in low frequency unsteady flow. Finally, the flow field of a single mode is reconstructed by superimposing the secondorder mode to the firstorder mode under the rated condition. And the results show that six highspeed fluid clusters dominated at the volute inlet rotate periodically with the blade. Generally speaking, the combination of DMD method and CFD unsteady calculation method is beneficial to understand the flow mechanism in the volute of centrifugal pump by extracting the coherent structure of unsteady flow.
Nomenclature
:  Nominal rate flow (m^{3}/h) 
:  Pump head for nominal rate flow (m) 
:  Nominal rotational speed (r/min) 
∆t:  Time step (s) 
Z:  Number of blades 
D_{1}:  Diameter of impeller inlet (m) 
D_{2}:  Diameter of impeller outlet (m) 
α:  Blade wrap angle 
β_{2}:  Blade outlet angle 
D_{3}:  Base circle diameter of volute (m) 
:  Pressure coefficient 
:  Static pressure value of the monitoring point at a certain time (Pa) 
:  Average value of static pressure in one cycle (Pa) 
:  Circumferential velocity at the impeller outlet (m/s) 
:  Observation data matrix 
N:  Total number of snapshots 
:  Data of the i th snapshot 
A:  Linear transformation matrix 
r:  Residual vector 
S:  Companion matrix 
U, V:  Unitary matrix 
Σ:  Diagonal matrix 
:  Eigenvalues of matrix S 
:  Eigenvectors of matrix S 
:  Corresponding frequency 
:  Growth rate 
Ф_{i}:  DMD modes 
:  A subspace 
W:  Matrix with as the column vector 
N:  Matrix with as the column vector 
:  Amplitude of the i th mode 
:  Correlation coefficient 
:  Mode Frobenius norm. 
DMD:  Dynamic mode decomposition 
FFT:  Fast Fourier transformation 
CFD:  Computational fluid dynamics. 
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors have approved the final version of the manuscript submitted. There is no financial or personal interest.
Acknowledgments
The authors are grateful to the financial support from the National Natural Science Foundation of China (Research Project No. 51866009).
References
 J. D. H. Kelder, R. J. H. Dijkers, B. P. M. Van Esch, and N. P. Kruyt, “Experimental and theoretical study of the flow in the volute of a low specificspeed pump,” Fluid Dynamics Research, vol. 28, no. 4, pp. 267–280, 2001. View at: Publisher Site  Google Scholar
 J. González, J. Fernández, E. Blanco, and C. Santolaria, “Numerical simulation of the dynamic effects due to impellervolute interaction in a centrifugal pump,” Journal of Fluids Engineering, vol. 124, no. 2, pp. 348–355, 2002. View at: Publisher Site  Google Scholar
 J. Keller, E. Blanco, R. Barrio, and J. Parrondo, “PIV measurements of the unsteady flow structures in a volute centrifugal pump at a high flow rate,” Experiments in Fluids, vol. 55, no. 10, p. 1820, 2014. View at: Google Scholar
 L. Tan, B. Zhu, Y. Wang, S. Cao, and S. Gui, “Numerical study on characteristics of unsteady flow in a centrifugal pump volute at partial load condition,” Engineering Computations (Swansea, Wales), vol. 32, no. 6, pp. 1549–1566, 2015. View at: Publisher Site  Google Scholar
 H. Alemi, S. A. Nourbakhsh, M. Raisee, and A. F. Najafi, “Development of new "multivolute casing" geometries for radial force reduction in centrifugal pumps,” Engineering Applications of Computational Fluid Mechanics, vol. 9, no. 1, pp. 1–11, 2015. View at: Publisher Site  Google Scholar
 F. J. Wang, L. X. Qu, L. Y. He, and J. Y. Gao, “Evaluation of flowinduced dynamic stress and vibration of volute casing for a largescale doublesuction centrifugal pump,” Mathematical Problems in Engineering, vol. 2013, Article ID 764812, 9 pages, 2013. View at: Publisher Site  Google Scholar
 Y. Liu and L. Tan, “Tip clearance on pressure fluctuation intensity and vortex characteristic of a mixed flow pump as turbine at pump mode,” Journal of Renewable Energy, vol. 129, pp. 606–615, 2018. View at: Publisher Site  Google Scholar
 Y. Hao and L. Tan, “Symmetrical and unsymmetrical tip clearances on cavitation performance and radial force of a mixed flow pump as turbine at pump mode,” Journal of Renewable Energy, vol. 127, pp. 368–376, 2018. View at: Publisher Site  Google Scholar
 T. Lei, Y. Zhiyi, X. Yun, L. Yabin, and C. Shuliang, “Role of blade rotational angle on energy performance and pressure fluctuation of a mixedflow pump,” Proceedings of the Institution of Mechanical Engineers, Part A: Journal of Power and Energy, vol. 231, no. 3, pp. 227–238, 2017. View at: Publisher Site  Google Scholar
 T. Lei, X. Zhifeng, L. Yabin, H. Yue, and X. Yun, “Influence of Tshape tip clearance on performance of a mixedflow pump,” Proceedings of the Institution of Mechanical Engineers, Part A: Journal of Power and Energy, vol. 232, no. 4, pp. 386–396, 2018. View at: Publisher Site  Google Scholar
 P. J. Schmid and J. Sesterhenn, “Dynamic mode decomposition of numerical and experimental data,” in Proceedings of the SixtyFirst Annual Meeting of the APS Division of Fliuid Mechanics, 2008. View at: Google Scholar
 P. J. Schmid, “Dynamic mode decomposition of numerical and experimental data,” Journal of Fluid Mechanics, vol. 656, pp. 5–28, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 P. J. Schmid, L. Li, M. P. Juniper, and O. Pust, “Applications of the dynamic mode decomposition,” Theoretical and Computational Fluid Dynamics, vol. 25, no. 14, pp. 249–259, 2011. View at: Publisher Site  Google Scholar
 A. Seena and H. J. Sung, “Dynamic mode decomposition of turbulent cavity flows for selfsustained oscillations,” International Journal of Heat and Fluid Flow, vol. 32, no. 6, pp. 1098–1110, 2011. View at: Publisher Site  Google Scholar
 M. S. Hemati, M. O. Williams, and C. W. Rowley, “Dynamic mode decomposition for large and streaming datasets,” Physics of Fluids, vol. 26, no. 11, 2014. View at: Google Scholar
 J. H. Tu, C. W. Rowley, D. . Luchtenburg, S. L. Brunton, and J. N. Kutz, “On dynamic mode decomposition: theory and applications,” Journal of Computational Dynamics, vol. 1, no. 2, pp. 391–421, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 C. Pan, D. Xue, and J. Wang, “On the accuracy of dynamic mode decomposition in estimating instability of wave packet,” Experiments in Fluids, vol. 56, no. 8, p. 164, 2015. View at: Google Scholar
 D. A. Bistrian and I. M. Navon, “The method of dynamic mode decomposition in shallow water and a swirling flow problem,” International Journal for Numerical Methods in Fluids, vol. 83, no. 1, pp. 73–89, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 K. Taira, S. L. Brunton, S. T. Dawson et al., “Modal analysis of fluid flows: an overview,” AIAA Journal, 2017. View at: Publisher Site  Google Scholar
 M. Liu, L. Tan, and S. Cao, “Dynamic mode decomposition of cavitating flow around ALE 15 hydrofoil,” Journal of Renewable Energy, vol. 139, pp. 214–227, 2019. View at: Publisher Site  Google Scholar
 M. Liu, L. Tan, and S. Cao, “Dynamic mode decomposition of gasliquid flow in a rotodynamic multiphase pump,” Journal of Renewable Energy, vol. 139, pp. 1159–1175, 2019. View at: Publisher Site  Google Scholar
 J. Kou and W. Zhang, “An improved criterion to select dominant modes from dynamic mode decomposition,” European Journal of Mechanics  B/Fluids, vol. 62, pp. 109–129, 2017. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Yibin Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.