Next Article in Journal
Instrumentalization of a Model for the Evaluation of the Level of Satisfaction of Graduates under an E-Learning Methodology: A Case Analysis Oriented to Postgraduate Studies in the Environmental Field
Next Article in Special Issue
Dielectric Elastomer Generator for Electromechanical Energy Conversion: A Mini Review
Previous Article in Journal
Gender Diversity in Research and Innovation Projects: The Proportion of Women in the Context of Higher Education
Previous Article in Special Issue
Electret Nanogenerators for Self-Powered, Flexible Electronic Pianos
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of the Hydrodynamic Unsteady Flow Inside a Centrifugal Fan and Its Downstream Pipe Using Detached Eddy Simulation

1
College of Engineering, Zhejiang Normal University, Jinhua 321004, China
2
Key Laboratory of Urban Rail Transit Intelligent Operation and Maintenance Technology & Equipment of Zhejiang Province, Jinhua 321005, China
3
School of Mechatronic Engineering, Quzhou College of Technology, Quzhou 324000, China
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(9), 5113; https://doi.org/10.3390/su13095113
Submission received: 12 April 2021 / Revised: 28 April 2021 / Accepted: 29 April 2021 / Published: 2 May 2021
(This article belongs to the Special Issue Environmental Energy Harvesting and Storage)

Abstract

:
The detailed unsteady turbulent flow inside a centrifugal fan and its downstream pipe was studied using detached eddy simulation (DES) at three flowrates, namely, the best efficiency point (BEP), 0.75BEP, and 1.49BEP. Both the mean and fluctuating flow fields were analyzed on the basis of the root-mean-square value as the indication of fluctuating intensity. Results showed that the pressure fluctuation had the minimum value at BEP, but the velocity fluctuation increased with the flowrate. Most regions inside the centrifugal fan underwent large pressure fluctuation with the magnitude of about 10~20% of pref = 0.5 ρu22, where u2 is the blade velocity at the impeller outlet. The pressure fluctuation had a maximum value at the impeller side of the tongue tip rather than the stagnation point, and it decreased rapidly along the outlet pipe with magnitude about 1% of pref after distance of five pipe diameters. The spectra of hydrodynamic pressure showed conspicuous spikes at the blade passing frequency (BPF) in the volute but not in the downstream pipe. At the downstream pipe entrance, pressure fluctuation spectra agreed with experimental results, showing that hydrodynamic pressure fluctuations were dominant; however, the experimental data showed a much slower decreasing rate due to the acoustic fluctuations.

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.

2. Materials and Methods

2.1. Description of the Centrifugal Fan

The main dimensions of the studying centrifugal fan are shown in Table 1. The impeller with the outlet diameter d = 400 mm had Zb = 12 forward-curved blades, and a vaneless rotating diffuser with the diameter d2 = 460 mm. The rotational speed was N = 2900 rpm; therefore, the rotational frequency was RF = N/60 = 48.3 Hz and the blade passing frequency was BPF = RF × Zb = 580 Hz. The blade velocity at the impeller outlet was u2 = 2π N/60 × d/2 = 60.73 m/s, and the reference dynamic pressure was pref = 0.5 ρu22 = 2259.5 Pa.
The airflow leaving from the impeller gathered in the volute housing and then entered the downstream pipe. The fan outlet had a rectangular section and was connected by a transitional pipe piece (length of 0.245 m), which gradually changed from the rectangular section to a circular section with a diameter of 190 mm. A pipe of 7 m long was connected to the circular end of the transitional pipe piece, and thus the pressure fluctuations in the downstream pipe could be investigated.
On the basis of standard ISO 5801:2007 industrial fans performance testing using standardized airways, we tested the fan performance in a standard test rig that had the uncertainty within 2%. The non-dimensional fan performance with the pressure rise and efficiency is shown in Figure 1. The flow coefficient was determined by φ = Q / ( u 2 π d 2 2 / 4 ) , where Q is volume flowrate. The pressure coefficient was found by ψ = Δ p t o t / p r e f ψ = Δ p t o t / p r e f , where Δ p t o t is the total pressure rises. The BEP flow rate falls in 1000~1200 m3/h. In this study, 1116 m3/h was taken as the BEP flowrate, i.e., φ = 0.0406.

2.2. The Computational Mesh and CFD Simulation

The meshing procedure was carried out by dividing the flow domain into four parts, namely, the inlet part, the impeller part (rotating unit), the volute, and the downstream pipe. Structured mesh in ANSYS ICEM-CFD was built to represent the four parts and was converted into unstructured computational meshes. The four mesh parts were then read into ANSYS FLUENT and combined to represent to the whole flow region. Interfaces between neighboring parts were used to interchange flow field data in different parts.
The computational mesh is shown in Figure 2. In the near-wall regions, very fine grid (at least 15 layers) was applied (shown in Figure 2), resulting in many cells. Several grid spaces were considered in successive refinements to check the influence of the grid density on the pressure and efficiency performance of the centrifugal fan, and finally a mesh with 8,975,832 cells was used in the unsteady simulation. The mesh statistics are shown in Table 2, in total, nearly 9 million cells were used. The CFD model without the downstream pipe was used in the study of Lu et al. [25], where LES was carried out. The mesh scheme was found to be not fine enough to performing LES, because according to the recommendation of Versteeg and Malalasekera [26], the non-dimensional wall distance of the first layer in the mesh y+ ≤ 1. However, in the simulation the condition, y+ ≤ 5 can be satisfied on most wall surfaces, which meets the requirement of the kω turbulent model.
The internal flows of the centrifugal fan can be regarded as incompressible because the Mach number with the circumferential velocity at the blade outlet u2/c = 0.179 (<0.3). Therefore, the incompressible N-S equations were used for computational procedure in this study. With the incompressible assumption, only the hydrodynamic pressure is obtained without the acoustic pressure. The compressible N-S equations can certainly simulate the sound waves as well as turbulence. However, it is quite challenging to numerically resolve the acoustic waves because acoustic intensity is oftentimes five to six orders smaller than the mean flow, and a numerical scheme must have extremely low numerical noise in order to accurately obtain the parameters of acoustic waves [27].
The detached eddy simulation was employed to model the turbulence. The mathematical model can be found in the Appendix A. The numerical simulation was realized using a commercial CFD package ANSYS FLUENT. The impeller domain was set to be rotating with a speed of 2900 RPM to include the centrifugal force source. The remaining three parts, i.e., the inlet part, the volute, and the downstream pipe, were stationary. According to the 0.75BEP, BEP, and 1.49BEP flow rates, the velocity was applied at the inlet with the magnitudes of 9.45, 12.6, and 18.73 m/s, respectively. Medium turbulent level with the turbulent intensity of 5% and the hydraulic diameters of 0.177 were set as the turbulence parameters. The pressure outlet boundary condition with the averaged static pressure of 0 Pa was applied at the pipe outlet surface. The pressure–velocity coupling was calculated using the SIMPLEC algorithm [26]. Second-order, upwind discretizations were used for convection terms, and bounded central differencing schemes were used for diffusion terms. For the time-dependent term scheme, we applied an implicit second-order scheme. The convergence criteria for the residuals of the continuity equation, the momentum equations, and the turbulence equations were set as 1.0 × 10−5 during iteration.
A steady-state simulation was initially carried out with the SST kω as the turbulence model, and the results served as the initial conditions for the following unsteady calculation. In the unsteady-state computation, the sliding mesh technique was applied to the interfaces in order to allow unsteady interactions between the impeller and the volute. A complete impeller revolution was divided into 2048 time-steps, i.e., the time step was equal to 1.010237 × 10−5 s. In the previous study [28], one impeller revolution was divided into 512 time steps in the URANS simulation with the standard kε equations as the turbulence modelling. Here, for DES, a much finer computational mesh was adopted, and therefore a smaller time step was needed to ensure the stability of the numerical calculation. In the study of [28], only the internal flow in the centrifugal fan was calculated, and with the acoustic analogy theory the aerodynamic noise was predicted, showing that the volute casing noise was predominant. However, in the present study, the unsteady flow in the fan as well as in the long downstream pipe was simulated in order to study the development and transportation of the pressure fluctuations.
The target convergence criteria for the residuals of the governing equations were 1.0 × 10−5 during iterations, and a max number of 30 iterations per time step. A representative section of convergence history at BEP flowrate is shown in Figure 3, and good convergence was found. In the unsteady-state simulation, the results of this simulation became statistically stable after several impeller revolutions. Then, the data sampling for time statistics in the ANSYS FLUENT solver was activated and the temporal pressure data of monitoring points were saved at every time step. A total of 16,384 time-steps (8 impeller revolutions) were sampled. The calculation of each flowrate took about 30 days in a cluster with 4 computer nodes. Each node had two Intel Xeon Gold 6132 processors, and one processor has 14 cores; therefore, 112 cores were used during the simulation in total.

3. Results and Discussion

From the time-step-averaged flow results, we obtained the values of total pressure increase Δ p t o t by the difference of the mass-averaged total pressure at the fan outlet and inlet. The torque of the impeller T 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 η = Q Δ p t o t / ( T 2 π N / 60 ) = 0.767 . 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
Q = 1 2 | Ω | 2 | S | 2
where Ω and S are the second invariant measure of the vorticity and strain rate tensors
Ω = 1 2 u ( u ) T , S = 1 2 u + ( u ) T
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]
Q * = Q / Q max
where Q max is the maximum Q value. Figure 6d shows the iso-surface of Q * = 0.01 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:
φ RMS = i = 1 n φ i φ ¯ 2 / n
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 5D 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 10D, 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.

Author Contributions

J.-C.C.: writing—original draft, writing—review and editing, funding acquisition, investigation, methodology; H.-J.C.: investigation; V.B.: writing—original draft, investigation. Y.-H.G.: supervision, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number 51976201) and Zhejiang Provincial Natural Science Foundation of China (grant number LY18E060006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that supports the findings of this study are available within the article.

Acknowledgments

The authors kindly acknowledge the support of Nanfang Ventilator Co., Ltd. (Foshan, Guangdong Province, China).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A Mathematical Model and Numerical Method

Appendix A.1. Governing Equations

For incompressible flow, the continuity equation and the momentum equation read
u = 0
u t + ( u ) u = 1 ρ p ¯ + ν 2 u
where u is the flow velocity vector; p is the pressure; and ρ and ν are the density and kinematic viscosity of the fluid, respectively. The Poisson equation of the pressure can be obtained by performing the divergence of the momentum Equation (A2) and considering the divergence free velocity condition (A1), leading to
2 p = ρ ( u ) u
The Poisson equation of pressure is elliptic with the velocity field of the sources, and the pressure is referred to as the hydrodynamic pressure, which decays quickly away from their source of generation. Unlike the hydrodynamic pressure fluctuation governed by the Poisson Equation (A3), the acoustic equation is
1 c 2 2 p t 2 - 2 p = 0
where c is the speed of sound, and p is the acoustic pressure. This is the partial differential equation of hyperbolic type, and the acoustic pressure fluctuations pass through the pipe without attenuation if damping is not taken into account. For low Mach number flow, the famous acoustic analogy theory by Lighthill [31] reformulates the continuity equation and momentum equation into a nonhomogeneous wave equation:
2 ρ t 2 c 2 2 ρ x i 2 = 2 T i j x i x j
With the Lighthill’s tensor,
T i j = ρ v i v j + p p ρ ρ c 2 δ i j σ i j
where the acoustic density variation is ρ = c 2 p , the sub means the reference value of ambient condition, δ i j is the Kronecker delta, and σ i j is viscous stress tensor. Thus, aerodynamic noise generated by the unsteady flow in the centrifugal fan can propagate through the downstream pipe as the acoustic waves causing acoustic pressure fluctuations.
Computational aeroacoustics is not the same as CFD, and it faces a different set of computational challenges, as were pointed out by Tam [27]. In the present study, the acoustic pressure was ruled out due to the incompressible flow model.

Appendix A.2. Turbulence Model

It can be concluded that the internal flow is turbulent in nature on the basis of the characteristic Reynolds number, Re = d u 2 / ν = 1.57 × 10 6 . In turbulence modelling, the flow is usually separated into the mean and fluctuating components as ϕ = ϕ ¯ + ϕ . With this decomposition, the governing N-S equations read as follows:
u ¯ = 0
u ¯ t + ( u ¯ u ¯ ) = 1 ρ p ¯ + ν 2 u ¯ 1 ρ τ R
where τ R is the turbulent or Reynolds stress tensor τ R = ρ u u ¯ . The overbar should be interpreted as the ensemble average procedure in RANS approach and as the spatial filtering operation in LES. The N-S equations are closed by boundary conditions: either a Dirichlet condition u ¯ or a Neumann condition σ n = t ¯ is specified, where σ and n are the stress tensor and unit outward normal vector, respectively, and t ¯ is the specified tensional force. With the help of the constitutive model of the fluid, the directional derivatives of u ¯ are specified in this manner.
An improved version of detached eddy simulation, namely, delayed detached eddy simulation (DDES), was used in this study to calculate the turbulence. In DES, the Reynolds stress is as follows by the Boussinesq approximation:
τ R = ρ u u ¯ = 2 μ T S ¯ 2 3 ρ k I
where S ¯ is the mean strain rate tensor, μT is the turbulent eddy viscosity, k is the turbulent kinetic energy, and I is the identity tensor.
In the context of RANS simulations, the SST kω turbulence model [26] was chosen for its advantages in simulating the near-wall flows. The SST kω turbulence mode employs the kω model in the boundary layer flows and switches to the kε model in free shear streams. For the kω model, the turbulent kinetic energy k and turbulent-specific dissipation ω are descripted by the transport equations [22]:
k t + ( u ¯ k ) = P k + ( ν + ν T / σ k ) k k 3 / l D D E S
ω t + ( u ¯ ω ) = ( ν + ν T / σ ω ) ω + α P k / ν T β ω 2 + 2 ( 1 F 1 ) σ ω 2 ( k ω ) / ω
where α , β , σ k , σ ω , σ ω 2 are model constants, the kinematic eddy viscosity ν T = μ T / ρ . The F1 is the first blending function for an automatic switching between the kω and kε models, and P k is the function of turbulent kinetic energy production:
P k = min ( τ R : u ¯ , 10 β k ω ) .
The turbulent eddy viscosity μ T is bounded by
μ T = ρ a 1 k max ( a 1 ω , S F 2 ) ,
where a 1 = β , S = 2 ( S ¯ i j S ¯ i j ) is the second invariant of the strain rate tensor, and F2 is a second blending function with the expression as
F 2 = tanh max 2 k β * ω y , 500 v y 2 ω 2 .
The DDES turbulent length scale in Equation (A10) can be written as follows:
l D D E S = l R A N S - f d max ( 0 , l R A N S l L E S )
where l L E S = C D E S h max , l R A N S = k / ( C μ ω ) , hmax is the maximum edge length of the cell., and
C D E S = C D E S 1 F 1 + C D E S 2 ( 1 F 2 )
The model constants of the above equations can be found in [22].

Appendix A.3. Numerical Method

The finite volume method is employed to discretize the integral form of the mass and momentum equations or every control volume VP at point P as follows:
f u ¯ f n f S f = 0
t u ¯ P V P + f u ¯ f u m n f u ¯ f S f = f ν f n f ( u ¯ ) f S f 1 ρ ( p ) P V P
where the subscript P represents the volume values, and f means the face values. The VP is the cell volume of each control volume, Sf is the face area, and n f is the unit outward normal vector of the face. The transport equations of the turbulent kinetic energy k and turbulent specific dissipation ω are discretized in the same manner as in Equation (A18) when RANS simulation is performed instead of LES. At the centroid of each cell, the variable values are calculated directly, and for variable values at the control volume surface, suitable interpolation is used with the center values. Finally, suitable quadrature formulae are applied to approximate the surface and volume integrals to a second-order accuracy in this study.

References

  1. Kim, K.Y.; Samad, A.; Benini, E. Design Optimization of Fluid Machinery: Applying Computational Fluid Dynamics and Numerical Optimization; John Wiley & Sons, Inc.: Singapore, 2019. [Google Scholar]
  2. Kjork, A.; Lofdahl, L. Hot-Wire Measurements inside a Centrifugal Fan Impeller. J. Fluids Eng. 1989, 111, 363–368. [Google Scholar] [CrossRef]
  3. Velarde-Suarez, S.; Ballesteros-Tajadura, R.; Santolaria-Morros, C.; Gonzalez-Perez, J. Unsteady Flow Pattern Characteristics Downstream of a Forward-Curved Blades Centrifugal Fan. J. Fluids Eng. 2001, 123, 265–270. [Google Scholar] [CrossRef]
  4. Zhang, H.; Yang, C.; Wang, W.; Chen, J.; Qi, M. Investigation on the Casing Static Pressure Distribution and Stall Behaviors in a Centrifugal Compressor with Volute. Int. J. Mech. Sci. 2019, 160, 318–331. [Google Scholar] [CrossRef]
  5. Ni, D.; Zhang, N.; Gao, B.; Li, Z.; Yang, M. Dynamic Measurements on Unsteady Pressure Pulsations and Flow Distributions in a Nuclear Reactor Coolant Pump. Energy 2020, 198, 117305. [Google Scholar] [CrossRef]
  6. Ballesteros-Tajadura, R.; Velarde-Suárez, S.; Hurtado-Cruz, J.P.; Santolaria-Morros, C. Numerical Calculation of Pressure Fluctuations in the Volute of a Centrifugal Fan. J. Fluids Eng. 2006, 128, 359–369. [Google Scholar] [CrossRef]
  7. Younsi, M.; Djerrada, A.; Belamri, T.; Menter, F. Application of the Sas Turbulence Model to Predict the Unsteady Flow Field Behaviour in a Forward Centrifugal Fan. Int. J. Comput. Fluid Dyn. 2008, 22, 639–648. [Google Scholar] [CrossRef]
  8. Spence, R.; Amaral-Teixeira, J. A Cfd Parametric Study of Geometrical Variations on the Pressure Pulsations and Performance Characteristics of a Centrifugal Pump. Comput. Fluids 2009, 38, 1243–1257. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, J.; Meng, L.; Zhou, L.; Luo, Y.; Wang, Z. Unsteady Internal Flow Field Simulations in a Double Suction Centrifugal Fan. Eng. Comput. 2013, 30, 345–356. [Google Scholar] [CrossRef]
  10. Gao, Z.; Zhu, W.; Lu, L.; Deng, J.; Zhang, J.; Wuang, F. Numerical and Experimental Study of Unsteady Flow in a Large Centrifugal Pump with Stay Vanes. J. Fluids Eng. 2014, 136, 1–10. [Google Scholar] [CrossRef]
  11. Yao, Z.; Yang, M.; Xiao, R.; Wang, F. Influence of Wall Roughness on the Static Performance and Pressure Fluctuation Characteristics of a Double-Suction Centrifugal Pump. Proc. Inst. Mech. Eng. Part A J. Power Energy 2018, 232, 826–840. [Google Scholar] [CrossRef]
  12. Menter, F.R.; Kuntz, M.; Langtry, R. (Eds.) Ten Years of Industrial Experience with the Sst Turbulence Model; Turbulence Heat & Mass Transfer 4; Begell House Inc.: Redding, CT, USA, 2003. [Google Scholar]
  13. Bozorgasareh, H.; Khalesi, J.; Jafari, M.; Gazori, H.O. Performance Improvement of Mixed-Flow Centrifugal Pumps with New Impeller Shrouds: Numerical and Experimental Investigations. Renew. Energy 2021, 163, 635–648. [Google Scholar] [CrossRef]
  14. Guala, M.; Hommema, S.E.; Adrian, R.J. Large-Scale and Very-Large-Scale Motions in Turbulent Pipe Flow. J. Fluid Mech. 2006, 554, 521–542. [Google Scholar] [CrossRef]
  15. Discetti, S.; Bellani, G.; Örlü, R.; Serpieri, J.; Vila, C.S.; Raiola, M.; Zheng, X.; Mascotelli, L.; Talamelli, A.; Ianiro, A. Characterization of Very-Large-Scale Motions in High-Re Pipe Flows. Exp. Fluid Sci. 2019, 104, 1–8. [Google Scholar] [CrossRef]
  16. Tan, L.; Zhu, B.; Wang, Y.; Cao, S.; Liang, K. Turbulent Flow Simulation Using Large Eddy Simulation Combined with Characteristic-Based Split Scheme. Comput. Fluids 2014, 94, 161–172. [Google Scholar] [CrossRef]
  17. Dutta, P.; Saha, S.K. Nityananda Nandi, and Nairit Pal. Numerical Study on Flow Separation in 90° Pipe Bend under High Reynolds Number by K-Ε Modelling. Eng. Sci. Technol. Int. J. 2016, 19, 904–910. [Google Scholar]
  18. Röhrig, R.; Jakirlić, S.; Tropea, C. Comparative Computational Study of Turbulent Flow in a 90° Pipe Elbow. Int. J. Heat Fluid Flow 2015, 55, 120–131. [Google Scholar] [CrossRef]
  19. Zimoń, M.J.; Skillen, A.; Liu, W.; Sawko, R. Surrogate Modelling of Transient Stratified Flow in a U-Shaped Pipe. Nucl. Eng. Des. 2020, 368, 110787. [Google Scholar] [CrossRef]
  20. Spalart, P.R. Detached-Eddy Simulation. Annu. Rev. Fluid Mech. 2009, 41, 181–202. [Google Scholar] [CrossRef]
  21. Nguyen, V.; Nguyen, H.H. Detached Eddy Simulations of Flow Induced Vibrations of Circular Cylinders at High Reynolds Numbers. J. Fluids Struct. 2016, 63, 103–119. [Google Scholar] [CrossRef]
  22. Gritskevich, M.S.; Garbaruk, A.V.; Schütze, J.; Menter, F.R. Development of Ddes and Iddes Formulations for the K-Ω Shear Stress Transport Model. Flow Turbul. Combust. 2012, 88, 431–449. [Google Scholar] [CrossRef]
  23. Zhang, N.; Liu, X.; Gao, B.; Xia, B. Ddes Analysis of the Unsteady Wake Flow and Its Evolution of a Centrifugal Pump. Renew. Energy 2019, 141, 570–582. [Google Scholar] [CrossRef]
  24. Zhang, N.; Jiang, J.; Gao, B.; Liu, X.; Ni, D. Numerical Analysis of the Vortical Structure and Its Unsteady Evolution of a Centrifugal Pump. Renew. Energy 2020, 155, 748–760. [Google Scholar] [CrossRef]
  25. Lu, F.A.; Qi, D.T.; Wang, X.J.; Zhou, Z.; Zhou, H.H. A Numerical Optimization on the Vibroacoustics of a Centrifugal Fan Volute. J. Sound Vib. 2012, 331, 2365–2385. [Google Scholar] [CrossRef]
  26. Versteeg, H.K.; Malalasekera, W. An Introduction to Computational Fluid Dynamics the Finite Volute Method, 2nd ed.; Pearson Education Limited: Essex, UK, 2007. [Google Scholar]
  27. Tam, C.K.W. Recent Advances in Computational Aeroacoustics. Fluid Dyn. Res. 2006, 38, 591–615. [Google Scholar] [CrossRef]
  28. Cai, J.C.; Qi, D.T.; Lu, F.A. Investigation of Fan Casing Aerodynamic Noise at the Blade Passing Frequency. In Proceedings of the ASME 2009 Fluids Engineering Division Summer Meeting, Vail, CO, USA, 2–6 August 2009. [Google Scholar]
  29. Cai, J.; Qi, D.; Lu, F.; Wen, X. Study of Tonal Fan Noise Reduction by Modification of the Volute Cutoff. Acta Acust. United Acust. 2010, 96, 1115–1124. [Google Scholar] [CrossRef]
  30. Wei, Y.; Zhu, L.; Zhang, W.; Wang, Z. Numerical and Experimental Investigations on the Flow and Noise Characteristics in a Centrifugal Fan with Step Tongue Volutes. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2019, 234, 2979–2993. [Google Scholar] [CrossRef]
  31. Glegg, S.; Devenport, W. Aeroacoustics of Low Mach Number Flows Fundamentals, Analysis, and Measurement; Academic Press: London, UK, 2017. [Google Scholar]
Figure 1. Fan performance: (a) pressure rise; (b) efficiency.
Figure 1. Fan performance: (a) pressure rise; (b) efficiency.
Sustainability 13 05113 g001
Figure 2. Fan internal flow region mesh: (a) the whole model; (b) the impeller; (c) snapshot of mesh at the blade trailing edge.
Figure 2. Fan internal flow region mesh: (a) the whole model; (b) the impeller; (c) snapshot of mesh at the blade trailing edge.
Sustainability 13 05113 g002
Figure 3. The iteration convergence history at BEP in the unsteady simulation.
Figure 3. The iteration convergence history at BEP in the unsteady simulation.
Sustainability 13 05113 g003
Figure 4. The mean pressure contour at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Figure 4. The mean pressure contour at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Sustainability 13 05113 g004
Figure 5. The mean velocity distribution at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Figure 5. The mean velocity distribution at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Sustainability 13 05113 g005
Figure 6. (a) Turbulent kinetic energy k (m2/s2), (b) turbulent specific dissipation ω (1/s), (c) vorticity in the z-direction at the middle span of the volute, and (d) the vortices identified by Q-criterion (Q* = 0.01) at 1.49BEP flowrate.
Figure 6. (a) Turbulent kinetic energy k (m2/s2), (b) turbulent specific dissipation ω (1/s), (c) vorticity in the z-direction at the middle span of the volute, and (d) the vortices identified by Q-criterion (Q* = 0.01) at 1.49BEP flowrate.
Sustainability 13 05113 g006
Figure 7. (a) Pressure and (b) velocity characteristics along the downstream pipe.
Figure 7. (a) Pressure and (b) velocity characteristics along the downstream pipe.
Sustainability 13 05113 g007
Figure 8. The volume-averaged RMS pressure in the computational domains.
Figure 8. The volume-averaged RMS pressure in the computational domains.
Sustainability 13 05113 g008
Figure 9. The RMS pressure contour at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Figure 9. The RMS pressure contour at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Sustainability 13 05113 g009
Figure 10. The volume-averaged RMS velocity in the computational domains.
Figure 10. The volume-averaged RMS velocity in the computational domains.
Sustainability 13 05113 g010
Figure 11. The RMS velocity distribution at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Figure 11. The RMS velocity distribution at the middle span of the volute: (a) 0.75BEP; (b) BEP; (c) 1.49BEP.
Sustainability 13 05113 g011
Figure 12. The (a) RMS pressure and (b) RMS velocity along the pipe center line.
Figure 12. The (a) RMS pressure and (b) RMS velocity along the pipe center line.
Sustainability 13 05113 g012
Figure 13. (a) The monitoring points and (b,c) the spectra of pressure fluctuations of 1.49BEP flowrate.
Figure 13. (a) The monitoring points and (b,c) the spectra of pressure fluctuations of 1.49BEP flowrate.
Sustainability 13 05113 g013aSustainability 13 05113 g013b
Table 1. The main dimensions of the centrifugal fan.
Table 1. The main dimensions of the centrifugal fan.
ParametersValue
Impeller outlet diameter, d2 (m)0.400
Impeller outlet width, b2 (m)0.036
Diffusor outlet t diameter, d3 (m)0.460
Impeller inlet diameter, d1 (m)0.155
Impeller inlet width, b1 (m)0.700
Blade outlet angle, β2 (°)126
Blade inlet angle, β1 (°)38
Radius to tongue (m)0.246
Number of blades, Zb12
Rotational speed, N (rpm)2900
Table 2. Mesh statistics of the fan.
Table 2. Mesh statistics of the fan.
NumberInlet PartImpellerVoluteDownstream Pipe
Cell1,003,6842,119,3082,947,7942,905,046
Node974,5792,024,4842,887,5262,869,296
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cai, J.-C.; Chen, H.-J.; Brazhenko, V.; Gu, Y.-H. Study of the Hydrodynamic Unsteady Flow Inside a Centrifugal Fan and Its Downstream Pipe Using Detached Eddy Simulation. Sustainability 2021, 13, 5113. https://doi.org/10.3390/su13095113

AMA Style

Cai J-C, Chen H-J, Brazhenko V, Gu Y-H. Study of the Hydrodynamic Unsteady Flow Inside a Centrifugal Fan and Its Downstream Pipe Using Detached Eddy Simulation. Sustainability. 2021; 13(9):5113. https://doi.org/10.3390/su13095113

Chicago/Turabian Style

Cai, Jian-Cheng, Hao-Jie Chen, Volodymyr Brazhenko, and Yi-Hong Gu. 2021. "Study of the Hydrodynamic Unsteady Flow Inside a Centrifugal Fan and Its Downstream Pipe Using Detached Eddy Simulation" Sustainability 13, no. 9: 5113. https://doi.org/10.3390/su13095113

APA Style

Cai, J. -C., Chen, H. -J., Brazhenko, V., & Gu, Y. -H. (2021). Study of the Hydrodynamic Unsteady Flow Inside a Centrifugal Fan and Its Downstream Pipe Using Detached Eddy Simulation. Sustainability, 13(9), 5113. https://doi.org/10.3390/su13095113

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

Article Metrics

Back to TopTop