1. Introduction
Fans, as typical fluid machinery devices, are widely employed for industrial and civilian usage; for transportation of working fluids in industry; for ventilation in buildings, tunnels, and cars; and for cooling in electric and electronic devices [
1]. Since fans are ubiquitous, they consume a huge amount of electricity, which is usually generated by fossil fuel. Improvement of fan efficiency can reduce the energy consumption so that it is helpful to environmental sustainability. Detailed investigation of the internal unsteady flows of fans is the first step to improve the fan performance and reduce fan noise.
Fans can be classified as either axial, centrifugal, or mixed flow types. Centrifugal fans can achieve high pressure rise in a short axial distance compared with axial fans due to the fact that the work done by the impeller is through the centrifugal flow. For centrifugal fans, the internal flow of the impeller can be very complex due to the influence of the centrifugal and Coriolis forces as well as the pressure gradients. The impeller rotation causes pressure fluctuations, which can lead to exciting the structures, such as the vibrations of the volute casing and the pipeline, as well as aerodynamic noise producing. Therefore, it is important and necessary to study the internal flow fields including both the mean flow and the fluctuating flow for long-term safe and quiet operation.
The flow fluctuations inside a centrifugal fan are mainly induced by the periodic interaction between the rotating impeller and the stationary casing, especially the volute tongue. The blade passing frequency (BPF) fluctuation induced by blades sweeping the volute tongue is usually the predominant component in the spectra. There are also broadband fluctuations due to the turbulent flow fluctuations in the inlet stream and boundary layers over structure surfaces.
In an earlier experimental study, Kjork and Lofdahl [
2] carried out the measurements of the three mean velocity components and the Reynolds stresses values inside a centrifugal fan impeller. They found that the mean flow can be characterized as an attached flow with almost linearly distributed velocity profiles, and for the turbulent flow, turbulent stresses with relatively low values predominate in the center region of the channel. Velarde-Suarez et al. [
3] experimentally investigated the flow at exit radial locations of a centrifugal fan with forward-curved blades using hot wire anemometry techniques and found a strong flow asymmetry with considerable changes in both magnitude and direction along the different circumferential positions, especially in circumferential positions closer to the volute tongue. Recently, Zhang et al. [
4] investigated the casing static pressure distribution of a centrifugal compressor ranging from the choke to stall operating conditions in order to find the relationship between the casing pressure evolution and stall behavior under different rotational speeds. Ni et al. [
5] studied the unsteady pressure and velocity fluctuations using laser doppler velocimeter in pump and discussed the relation between the unsteady pulsation and flow structure.
Computational fluid dynamics (CFD) has been a powerful tool to simulate the unsteady 3D flowfields of turbomachinery. Ballesteros-Tajadura et al. [
6] carried out a simulation of the complete 3D unsteady flow inside an centrifugal fan, and analyzed the pressure fluctuations in some locations over the volute wall, showing an important BPF peak. Younsi et al. [
7] solved the unsteady Reynolds-averaged Navier–Stokes equations (URANS) with scale adaptive simulation (SAS) as the turbulence modelling to study the overall performances of a centrifugal fan and the wall pressure fluctuations upon the volute casing surface. Spence and Amaral-Teixeira [
8] performed a CFD parametric study of the effect of geometrical variations on the pressure pulsations and performance characteristics of a centrifugal pump, finding that the cutwater gap and vane arrangement had the greatest influence across the various monitored locations and the flow range. Jing et al. [
9] performed simulation using the URANS approach with the shear stress transport (SST)
k−
ω turbulence model to study the interactions between the non-uniform impeller flow and the fixed volute through which the significant pressure fluctuations arise. Gao et al. [
10] studied the unsteady flow in a large centrifugal pump with stay vanes on the basis of numerical and experimental approaches. Yao et al. [
11] used detached eddy simulations (DES) [
12] to investigate the influence of wall roughness on the static performance and pressure fluctuations of a double-suction centrifugal pump. Bozorgasareh et al. [
13] designed a new impeller configuration with innovative shrouds, studying its influence on the pressure head and efficiency of the centrifugal pump from both experimental and numerical approaches.
There are only a few works studying the unsteady flow fields in the downstream pipe of a centrifugal pump or fan, although in many applications, pumps and fans are installed in the pipeline. At the same time, many papers concerned pipe flow exclusively, and the 90° pipe bends in particular can be found. Guala et al. [
14] and Discetti et al. [
15] studied the large-scale and very large scale motions in turbulent pipe flow using particle image velocimetry (PIV) measurement. Tan et al. [
16] performed large eddy simulation (LES) combined with a characteristic-based split scheme to study the flow in two circular sections of 90° pipe bends, observing an additional pair of vortexes near the curved section of the pipe inner side. Dutta et al. [
17] numerically studied flow separation in 90° pipe bend with high Reynolds number by
k−
ε turbulence modelling. Röhrig et al. [
18] made a comparative computational study of turbulent flow in a 90° pipe elbow, demonstrating the superiority of LES over the RANS approaches; however, this was at the cost of significant increase in computational effort. Zimoń et al. [
19] computed the mixed convection flow through a u-shaped bend with the transient thermal boundary conditions applied at the inlet.
Regarding this fact, in this work, the full 3D unsteady turbulent flow at three flow rates inside a centrifugal fan and its downstream pipe was simulated though CFD calculations by DES. The detailed time-averaged and fluctuating flow fields were analyzed, especially the characteristics of the hydrodynamic pressure fluctuations in the downstream pipe at different flowrates. DES approach was employed to address the challenge of internal unsteady flow with high Reynolds number and massive flow separations. DES combines LES and RANS, spurred by the belief that each alone was powerless to solve the problem at hand [
20].
Although LES showed good performance in capturing unsteady flow features, the computational cost is much higher than in industrial widely adopted unsteady RANS models [
21]. For the no-slip condition on solid walls, the LES-filtered N-S equations are integrated to the wall, which requires fine cells with near-wall cell points
y+ ≤ 1. Moreover, the resolution needs in the outer region of the boundary layer are very high, i.e., with at the least 20 points per thickness
δ in each direction [
20]. In the present study, the centrifugal fan had many wall surfaces, which leads to enormous computational cost of pure LES; therefore, the DES, specifically the delayed DES (DDES) [
22], was employed here. Recently, Zhang et al. [
23,
24] applied DDES on the basis of SST
k−
ω as RANS approach to study the unsteady flow in a centrifugal pump, and good agreement with the PIV experimental measurement results was found.
3. Results and Discussion
From the time-step-averaged flow results, we obtained the values of total pressure increase
by the difference of the mass-averaged total pressure at the fan outlet and inlet. The torque of the impeller
was received through the integration of moment in z-direction of the pressure and friction force over the impeller walls. The fan efficiency was calculated on the basis of the equation
. The simulation results of fan performance were compared with the experimental results as shown in
Figure 1. Good agreement can be found showing the validity of the application of DDES in the centrifugal fan. The numerical results showed slightly better fan performance. It may have due to the fact that other flow loss such as internal leakage had not been resolved well in the flow simulation.
3.1. Mean Flow Field
The mean pressure contour on the middle span plane of the volute is shown in
Figure 4, where the non-dimensional pressure was obtained by dividing current value of pressure on
pref. During the simulation, the obtained results were saved on different final time steps, which led to different pitch angles of the impeller, but it did not affect the general conclusions of the flow fields with different flowrates. As can be seen from the figure, static pressure rose from the impeller inlet along the impeller and further increased in the downstream pipe where part of the dynamic pressure was converted into static pressure. The pressure value on the blade pressure side was higher than on the suction side. The highest static pressure value was observed around the tongue tip, where the stagnation occurred. According to the velocity triangles in turbomachinery, the outlet flow angle increased with large flow rate. As a result, the stagnation point shifted from the outlet pipe side to the impeller side with the increasing of the flowrate. The impeller side of the volute tongue suffered a low-pressure region, where the air was sucked back into the volute with the increase of the velocity.
The averaged velocity magnitude distribution on the middle span of the volute is shown in
Figure 5. As we can see from the figure in the tongue region, the flow was split at the tongue tip, resulting in the stagnation point therein at the nominal flow rate. The maximum absolute velocity was located at the blade trailing edge regions, which can reach 1.5 times of the circumferential velocity at the blade trailing edge
u2 = 60.73 m/s. Forward-curved blade can lead to large absolute velocity value at the impeller outlet according to the velocity triangle of a rotor in turbomachinery theory.
At the impeller outlet, the velocity on the blade pressure side was larger than its value on the suction side. The circumferential non-uniform patterns of the velocity and pressure at the impeller outlet generated the regular pressure fluctuations on the volute casing surface, especially on the volute tongue surface that was close to the impeller. When a blade passed by, a pulsation occurred, and thus the pressure fluctuations had the BPF as the fundamental frequency.
At 0.75BEP flowrate, the flow leaving the impeller had a small flow angle. As a result, a vortex region existed at the lower part of the transitional pipe, as can be seen in
Figure 5a. At 1.49BEP flowrate, the flow angle increased, and the flow tended to attack the lower part of the outlet pipe. Thus, a vortex region formed at the upper part of the transitional pipe (see
Figure 5c).
The turbulent kinetic energy
k, the turbulent specific dissipation
ω, and the vorticity were similar for the three flowrates. Therefore, only the 1.49BEP case is shown in
Figure 6, which was taken as the example for explanation. Large turbulent kinetic energy exited in the places of the volute tongue region and the impeller outlet. The lower part of the transitional section of the downstream pipe also suffered intensive turbulence. Checking the streamlines in
Figure 5c we found that vortices formed in that region due to the divergent shape of the transitional pipe.
As we see from
Figure 6b large turbulent specific dissipation values
ω were observed near wall regions where large velocity gradients existed. In the blade channel, the turbulent-specific dissipation was also quite large. The reason can be found by analyzing z-direction vorticity
ωz in
Figure 6c Large negative values of
ωz were observed around the blade trailing edges, on the pressure side of the beginning part of blades, suction side of the blade middle part, while large values of positive value
ωz existed in the blade channels. In large vorticity regions, fluid particles had a large velocity difference resulting in large turbulent energy dissipation.
The 3D vortical structure can be examined using the Q-criterion, which is represented as
where
and
are the second invariant measure of the vorticity and strain rate tensors
A positive
Q value indicates that the rotating effect is stronger than the strain rate effect. The non-dimensional value of
Q is calculated by [
24]
where
is the maximum
Q value.
Figure 6d shows the iso-surface of
with the
prms as the contour color. In order to see the vortices more clearly, we took the view angle from the hub side because of the vortex sheet on the shroud side. Here, we can see the large vortical structures on the blade suction sides and in most blade channels. Vortices also appeared around the impeller wheel exit, where airflow discharged into the volute, causing expansion. At the beginning of the spiral volute, intensive vortices existed, especially in the tongue region, where interactions between the rotating impeller and the stationary volute was the strongest.
The distributions of the pressure and velocity along the central line of the downstream pipe are shown in
Figure 7, in which the distance refers to the distance from the inlet section of the transitional pipe. The downstream pipe had an inner diameter of
D = 190 mm, and its wall thickness was 5 mm; therefore, the outer radius was
D = 200 mm. As one can see, the mean pressure increased through the transitional pipe, owing to its diffuser shape, and then decreased gradually due to the flow loss along the pipe. The velocity decreased through the transitional pipe part and reached a constant value where the pipe had a uniform circular section with a diameter of 190 mm.
3.2. Fluctuating Flow Field
The root-mean-square (RMS) value was taken as the indicator of fluctuating strength. The RMS value of time history data of any quantity
is defined as follows:
where
means the average value. Note that with this equation, the RMS value indicates only the fluctuations, since the mean value is subtracted. The volume-averaged RMS pressure of the four computational domains is shown in
Figure 8. The pressure fluctuation was the smallest at BEP. At 0.75BEP, flow rate pressure fluctuations increased slightly, while at 1.49BEP, flow rate the pressure fluctuations increased noticeably, especially in the impeller domain. In the downstream pipe, the pressure fluctuation was much smaller in comparison with the other parts, which denotes that hydrodynamic pressure fluctuations decreased quickly in the downstream pipe.
The RMS pressure distribution on the middle span of the volute is shown in
Figure 9 Inside the impeller, the RMS pressure was around 10–20% of the reference dynamic pressure
pref. The pressure fluctuation was more intensive on the blade suction side in comparison with the pressure side for forward curved blades. In this case, flow separation can occur more easily on the blade suction side, especially in the large flowrate condition when the flow attacks the blade pressure side. Flow separation can lead to flow fluctuations.
The maximum pressure fluctuations occurred at the volute tongue, wherein the strongest interaction and unsteadiness were induced by the impeller rotation relative to the fixed volute. The RMS pressure value was able to reach above 30% of
pref at 1.49BEP flowrate. The location of the maximum pressure fluctuations was not at the stagnation point but at the impeller side of the tongue where the air flow was suctioned back to the volute. In consequence, the tongue region was the main aerodynamic sound source of centrifugal fans. Modifications of the volute tongue can reduce the fan noise effectively, for example, through the usage of inclined-tongue [
29] and step tongue [
30].
The volume-averaged RMS value distribution of velocity is shown in
Figure 10 with logarithmic scale usage. The velocity in the impeller was the relative velocity observed from the rotating impeller, and its fluctuation was quite large in comparison with the absolute velocity fluctuations in the stationary parts. As we can see from the figure inside the centrifugal fan as well as the downstream pipe, the intensity of velocity fluctuation increased with the flowrate. The velocity fluctuation was sufficiently small in the inlet part. It is explained by the fact that the airflow suffered fluctuations when it went through the rotating impeller, and the fluctuations were carried to downstream with little influence on the inlet part at the upstream. In the inlet part, the velocity fluctuation was the smallest at BEP because the airflow angle was consistent with the blade inlet angle, causing the weakest incidence.
The RMS value distribution of velocity magnitude on the middle span of the volute is shown in
Figure 11. As can be seen, the impeller suffered great velocity fluctuations during the rotation, and its magnitude was approximately in the same order with
u2. Unlike the pressure fluctuation, which had the minimum value at BEP, velocity fluctuations increased with the flowrate. Intensive fluctuations were found at the impeller outlet and in the volute tongue region, as was expected. The magnitude was around 20% of
u2 for BEP and 0.75BEP flowrates; however, for 1.49BEP flowrate, it could be over 30% of
u2 in the volute tongue region. For low flow rates, large velocity fluctuation was detected at the lower part of the transitional pipe, and according to
Figure 5a,b, we can see that the flow separates in that region due to the diffuser configuration, hence vortices appeared. For large flow rate, large velocity fluctuation exited at the upper part of the transitional pipe because the flow leaving the impeller had a larger flow angle and tended to attack the low part of the outlet pipe, causing a wake region at the upper part of the transitional pipe, as explained in
Section 3.1.
After the flow entered the downstream pipe, it became smoother in a way that the pressure fluctuations decreased sharply, as can be seen in
Figure 12, where the RMS pressure and velocity values along the pipe centerline are shown. The pressure fluctuation intensity may have increased slightly at the inlet of the transitional pipe, which can be ascribed to the enhancement of the local flow disturbance due to the divergent shape, as shown in
Figure 5. Then, the pressure fluctuations decreased sharply after the transitional section of the downstream pipe. For the 0.75BEP and BEP flowrates, the pressure fluctuation decreased to 0.5%
pref after 5D, and about 1%
pref for the 1.49BEP flowrate. The decaying trend of the velocity fluctuations demonstrated a similar pattern with that of the pressure fluctuations, except that in the transitional pipe section, the velocity fluctuation decreased slower than the pressure fluctuations. This means that local changes of flow passages had more influences on velocity fields than pressure fields.
Experiment measurements of the pressure fluctuations in the downstream pipe were carried out. The downstream pipe was made of PVC material. It was drilled easily with small holes of 3 mm diameter along the axial direction. The Kulite miniature pressure transducer with the model XTL-140M was used to measure the dynamic pressure. Its tap diameter was 2.54 mm, and it was inserted into the pipe with the depth of 12.5 mm. The analog voltage signal was acquired by Simcenter SCADAS XS with the 24-bit A/D converter. The measurement was repeated eight times, and the standard deviation of the results was found to within 5%.
The measured pressure fluctuation at the 1.49BEP flowrate is also shown in
Figure 12a. Here, it should be pointed out that the numerical results were along the centerline while the experimental results were received along the line, which was 12.5 mm away from the wall surface. We can see from the figure that numerical pressure fluctuation had the same order of magnitude in comparison with the numerical result. Moreover, the trend was also similar. However, the experimental data showed higher pressure fluctuation intensity through the whole downstream pipe. This may be ascribed to the acoustic signals in the experimental data, which were completely ignored in numerical simulations with incompressible flow models.
From the study, we can conclude that the hydrodynamic fluctuation faded quickly as soon as the flow left the centrifugal fan, and acoustic fluctuations were dominant thereafter.
3.3. Spectra of Pressure Fluctuations
Several monitoring points were set inside the flow domain to record the pressure time history data (see
Figure 13a). Vlt1~Vlt7 were seven points on the casing wall along the volute profile at the mid-span of the volute: Vlt1 was taken at the tongue tip, with every 60° in the circumferential direction setting for Vlt2~Vlt6, and Vlt7 was 30° away from Vlt6. Pp00~Pp10 were taken along the centerline of the downstream pipe: Pp00 was at the fan outlet, Pp0 was at the transitional pipe outlet, Pp1 was 0.2 m from the Pp0, and Pp5 and Pp10 were respectively 1 and 2 m from Pp0.
The 1.49BEP flowrate was selected as an example to show the spectra of the pressure fluctuations of monitoring points (see
Figure 13b,c). For the pressure fluctuations in the volute, the discrete components of BPF and its second and third harmonics were noticeable in the spectra. The pressure fluctuation amplitudes of the monitoring points in the volute were larger than those in the downstream pipe, and the volute tongue had the most intensive flow fluctuations, as we already expected. At the volute tongue, the BPF component was able to reach a value of 330 Pa, which was nearly 15% of
pref, and this component decreased quickly, which can be described in mathematical way according to the elliptical Poisson equation of pressure (see Equation (A3) in the
Appendix A). After 5
D distance in the downstream pipe (
D is the pipe diameter), the BPF component reduced to 2.9 Pa, i.e., less than 0.2% of
pref. At the beginning of the downstream pipe, low-frequency hydrodynamic pressure fluctuations were noticeable, but after distance was equal to 10
D, the pressure fluctuations decreased to being lower than 0.1% of
pref, similar to the BPF component.
From the experimental results shown in
Figure 13c, we can see that at the beginning of the downstream pipe, they agreed with each other to a certain extent, showing that hydrodynamic pressure fluctuations were dominant. The hydrodynamic pressure fluctuations, especially the low-frequency components, reduced quickly along the pipe. However, the pressure fluctuations decreased much more slowly than the hydrodynamic pressure fluctuations because most fluctuations were related to acoustic pressures. From the measured pressure spectra, we can see that at Pp5 and Pp10, the spectra of experimental results were similar; thus, the overall intensity of pressure fluctuations was almost the same at these two points, as shown in
Figure 11a.
4. Conclusions
The unsteady hydrodynamic flow fields of a centrifugal fan with a long downstream pipe, at BEP, 0.75BEP, and 1.49BEP flowrates, were analyzed in detail using CFD with the delayed detached eddy simulation (DDES) as the turbulence modeling. The turbulence and vorticity characteristics are discussed. The flowfields revealed the characteristics of impeller–volute interactions, which generated the pressure and velocity fluctuations with the blade-passing frequency (BPF) as the predominant component.
The time-averaged flow fields clearly showed that the pressure increased gradually through the impeller, further recovered from the velocity in the volute, and gradually decreased along the downstream pipe due to friction. Owing to the forward-curved blades, at the impellor exit, the flow had very large velocity, and the maximum velocity near the trailing edge was able to reach 1.5 u2, where u2 is the blade velocity at the impeller outlet. The root-mean-square value pressure distributions showed that most places inside the centrifugal fan underwent large pressure fluctuations with a magnitude of about 10% of the reference dynamic pressure pref = 0.5 ρu22; the maximum value located at the tongue tip was able to reach 30% of pref. Unlike the pressure fluctuation, which had a minimum value at BEP, velocity fluctuations increased with the flowrate. The velocity fluctuation was intensive at the impeller outlet, and its magnitude was around 20% of u2 in the volute tongue region, but for 1.49BEP flowrate, it can be over 30% of u2 in the tongue region.
The pressure fluctuation magnitude decreased rapidly along the outlet pipe: after 5D (D is the outlet pipe diameter), the magnitude was about 0.5% pref for the 0.75BEP and BEP flowrates, and 1% pref for the 1.49BEP flowrate. The discrete components at BPF and its second and third harmonics were prominent in the pressure fluctuation spectra at the monitoring points in the volute. The BPF component had a maximum value of 15% of pref in the tongue region at the 1.49BEP flowrate, and it decreased sharply along the downstream pipe with an amplitude less than 0.2% of pref after 5D distance.
The spectra of the hydrodynamic pressure were compared with the measured pressure fluctuations at 1.49BEP flowrate. At the beginning of the downstream pipe, they agreed with each other to a certain extent, showing that hydrodynamic pressure fluctuations are dominant wherein, however, the experimental pressure fluctuations decreased more slowly than the hydrodynamic pressure due to the acoustic pressure waves in the downstream pipe.
This study focused on the internal flow of the centrifugal fan as well as its downstream pipe, and it is expected to contribute to the deep understanding of the hydrodynamic unsteady flow. Future study will include acoustic fluctuations that will be dominant in places away from disturbance such as in the downstream pipe.