Document Type : Research Paper
Authors
 Seyed Masoud Vahedi ^{1}
 Farzad Parvaz ^{} ^{2}
 Mohsen Khandan Bakavoli ^{1}
 Mohammad Kamali ^{3}
^{1} M.S. Student, Department of Mechanical Engineering, Semnan University, Semnan, Iran
^{2} PhD. Candidate, Department of Mechanical Engineering, Najafabad Branch, Islamic Azad University, Isfahan, Iran
^{3} M.S. Student, Department of Mechanical Engineering, Najafabad Branch, Islamic Azad University, Isfahan, Iran
Abstract
Separation of suspended droplets in a fluid flow has been a great concern for scientists and technologists. In the current study, the effect of the surface roughness on flow field and the performance of a gasoil cyclone is studied numerically. The droplets and the turbulent airflow inside the cyclone are considered to be the discrete and continuous phases respectively. The Reynolds stress model (RSM) is employed to simulate the complex, yet strongly anisotropic, flow inside the cyclone while the EulerianLagrangian approach is selected to track droplet motion. The results are compared to experimental studies; according to the results, the tangential and axial velocities, pressure drop, and Euler number decrease when the surface roughness increases. Moreover, the cyclone efficiency drops when the vortex length decreases as a result of a rise in surface roughness. The differences between the numerical and experimental results become significant at higher flow rates. By calculating the impact energy of droplets and imposing the filmwall condition on the walls, splash does not occur.
Keywords
Main Subjects
1. Introduction
An aerocyclone is a device to which a gasoil mixture is introduced from the inlet channel(s) located near the upper section, and the separation happens in the lower part (Brar et al., 2015). The separation happens as a result of the centrifugal and gravitational forces. Due to their simplicity and ease of installation, aerocyclones are used in a wide range of applications in gas and petroleum industries. The prediction of cyclone performance relies on a good understanding of flow behavior and geometric parameters which influence the flow pattern (Elsayed and Lacor, 2011).
Many researchers have estimated the cyclone efficiency by presenting various mathematical models which are based on physical phenomena happening inside the cyclone (Brar et al., 2015; Baltrenas and Chlebnikovas, 2015). Some neglect the surface friction in their numerical investigations of cyclone aerodynamics (Shukla et al., 2013; Xianget al., 2001; Hosseini et al., 2015). The surface friction of a cyclone depends on roughness characteristics due to corrosion and/or collision and the accumulation of particles on the surface. Therefore, the surface roughness is an important parameter remarkably affecting the cyclone performance. The friction factor was first presented by Barth (1956) who calculated the tangential velocity of the gas in the vortex flow. He also evaluated the pressure drop and the cyclone performance for the Rankine vortex region. The relation extracted by Barth for pressure drop consisted of the wall and outlet loss terms. Barth’s model was extended by Muschelknautz and Krambrock (1970) when they modified the friction coefficient and investigated the friction impact by injecting and adding solid materials to cyclone walls; their model estimated the efficiency more accurately. Meissner and Loffler (1978) derived an empirical expression for the tangential velocity which took the wall friction into account. This model also included two additional friction factors. Karagoz and Avci (2005) presented another model to predict the pressure drop and the efficiency of the cyclone. Unlike other studies, their friction coefficient was a function of the Reynolds number and the surface roughness. They showed that the wall friction reduces the vortex flow, and as a consequence, the pressure drop decreases, which was considered to be constant in all models. Zhou et al. (2017) studied the effect of roughness on the vortex length and the efficiency of gasparticle cyclones. Their results showed that some particles cause surface abrasion, and as a result, they significantly increase surface roughness. Also, they reported that the appropriate cyclone surface roughness should be controlled to obtain the optimal separation performance for an optimized cyclone separator. Kaya et al. (2011) investigated the flow field and efficiency of two types of cyclone using commercial CFD software. There was good agreement between their numerical simulation results and experimental data. Wang et al. (2006) studied the impacts of particle size and gas inlet velocity on the performance and gas flow field of the Lapple cyclone. They indicated that an increase in particle size and volumetric gas flow rate raises the cyclone performance. A particleladen flow was simulated in three different cyclones by Shi and Bayless (2007). They verified their results by comparing them with the experimental data. Raoufi et al. (2008) and Kępa (2010) studied the flow field inside the cyclones with different outlet geometries and concluded that these geometries significantly affect the performance of the aerocyclones.
In the present study, a comprehensive research is conducted to evaluate the effects of roughness, inlet velocity, and droplet diameter on the efficiency of a Stairmand cyclone. The dimensionless distribution of tangential and axial velocities are illustrated for three different materials with three distinct roughness values. The performance of the cyclone was studied at various droplet sizes in a gasdroplet mixture. The RSM model was utilized to simulate the turbulent flow, whereas the EulerianLagrangian statistical approach was employed to investigate the cyclone performance. The length of vortex flow, which is an essential factor in the separation of water droplets, is also calculated. Further, the probability of splashing is considered.
2. Conceptual model
The Stairmand cyclone consists of four parts, namely inlet channel, central cylinder, vortex finder, and the cone section (see Figure 1). The dimensions of these parts are listed in Table 1. The twophase flow of fluid and water droplet are blown through the cyclone inlet channel tangentially. The tangential velocity decreases when the current moves down the cyclone until it reaches the cone section where the decreasing trend of tangential velocity is inverted. In the cone section, where the second vortex formation happens, lighter water droplets are moved by the centrifugal force originated from tangential velocity towards the flow core, and they ascended to the cyclone outlet.
Table 1
Geometrical properties of the cyclone investigated in the present work.
Parameters 

Dimensional (m) 
0.290 
0.145 
0.058 
0.145 
0.145 
0.108 
0.435 
0.1160 
0.435 
0.290 
Nondimensional 
1 
0.5 
0.2 
0.5 
0.5 
0.375 
1.5 
4 
1.5 
1 
Figure 1
A schematic representation of Stairmand cyclone geometry.
3. Governing equations
3.1. Continuous phase
The fluid flow inside a cyclone is threedimensional, turbulent, and isothermal. Since the maximum fluid velocity in the cyclone is at the inlet, the Mach number is equal to 0.06; therefore, the flow can be considered incompressible. The Reynolds averaged NavierStokes and continuity equations can be written as follows (Parvaz et al., 2017; Gao et al., 2013);
(1) 

(2) 

3.2. Turbulence modeling
Since the flow inside the cyclone is a highspeed swirling flow, the RSM model is applied. Streamlines with high curvatures, swirl, rotation, and rapid changes in the strain rate could be simulated in a more rigorous manner with RSM than with oneequation and twoequation turbulent models (Parvaz et al., 2017; ANSYS Fluent User’s Guide 2013; Sun et al., 2017). Thus, the following relations can be used:
(3) 
The first term on the righthand side of the above equation represents turbulent diffusion term which is expressed by:
(4) 
The second term refers to stress production term:
(5) 
The third term is also known as the pressurestrain correlation term:
(6) 
The last term is called the dissipation term:
(7) 
3.3. Discrete phase
The effect of droplets on the change in flow field can be neglected when the assumption of the low volume fraction (less than 10%) is considered for the droplets. The governing equation for droplet motion is given as follows:
(8) 
where, is droplet velocity, and stands for the drag force imposed to droplets. , g, and represent droplet density, gravitational acceleration, and other forces like Brownian and Soffman’s lift forces respectively (Kaya et al., 2011). In the current study, these two latter forces are not considered. For spherical droplets, the drag force can be written by:
(9) 
where, the Reynolds number is calculated by (Parvaz et al., 2017):
(10) 
Four regimes of droplet behavior, namely stick, rebound, spread, and splash regimes, could happen for droplets while the chance of occurrence depends on the impact energy and the cyclone wall temperature (see Figure 2).
Figure 2
Four regimes of droplet behavior (ANSYS Fluent User’s Guide 2013).
Splashing has been studied in the current work at a constant wall temperature of 273 K. Below the boiling temperature of the liquid, the droplets can either stick, spread, or splash, whereas above the boiling temperature, the droplets can either rebound or splash. It should be noted again that the occurrence of these regimes depends on the impact energy and boiling temperature of the liquid. The dimensionless impact energy of the droplets can be defined as:
(11) 
When becomes less than 16, the sticking regime is applied, which means that the velocity of the droplet is the same as the wall velocity. Nevertheless, at wall temperatures higher than the boiling temperature of liquid, the impact energy is less than the critical energy. Splashing happens when the impact energy is higher than a critical energy threshold 57.7 (ANSYS Fluent User’s Guide 2013).
3.4. Wall roughness
The flow resistance, specifically because of surface roughness, causes momentum transport, especially in the boundary layer. Hence, the law of the wall is influenced by the surface roughness, and it can be modified as follows:
(12) 
where, (VonKarman constant) equals 0.41, and B is an empirical constant for smooth walls set at 5.56 (Sun et al., 2017).
For the hydraulically rough flows, the roughness function is expressed in:
(13) 
where, is a constant equal to 0.5 (Sun et al., 2017).
3.5. Vortex length
The outer vortex flow axis changes when its distance from the vortex finder increases. This downward vortex flow is known as the free vortex flow as schematically demonstrated in Figure 3. The free vortex length is calculated using the radial distribution of pressure inside the cyclone. This parameter is rarely considered to investigate the droplets efficiency; however, it has a vital role in cyclone performance estimation (Alahmadi and Nowakowski, 2016; Alexander, 1949; Bryant et al., 1983; Ji et al., 1991).
Figure 3
The free vortex length of a cyclone.
Although the vortex length is almost unpredictable, a number of correlations are presented for it. Alexander (1949) presented a relation for the free vortex length of a glassy cyclone (30.2 to 200 mm).
(14) 
Bryant et al. (1983) and Ji et al. (1991) also reported empirical Equations 15 and 16 respectively:
(15) 

(16) 

3.6. Boundary conditions
After preparing the geometry and grids, the boundary conditions should be specified (Sabziani and Sari, 2015). The air flow enters the cyclones tangentially at a velocity of . The air flow exits from the cyclone through the vortex finder mounted on top of the cyclone. Noslip boundary condition is considered for all the walls. The bottom section is regarded as a trap to study the cyclone performance (Parvaz et al., 2017). The wallfilm model is applied to the walls so as to investigate the splashing. The wall roughness of the three different materials used in producing the cyclone is listed in Table 2.
Table 2
The surface roughness of the substances used in the cyclone (Shames, 1982).
Material 
Glass 
Steel (wrought iron) 
Galvanized iron 
Roughness (mm) 
0.0003 
0.046 
0.15 
4. Solver settings
The detailed information used in the present CFD modeling such as the discretization scheme of pressure, momentum, turbulent kinetic energy, turbulent dissipation rate, the Reynolds stress, and pressure/velocity coupling are tabulated in Table 3.
Table 3
Numerical setting of the current simulation.
Numerical setting 
Pressure discretization 
Pressure velocity coupling 
Momentum discretization 
Turbulent kinetic energy 
Turbulent dissipation rate 
Reynolds stress 
Scheme 
Body force weighted 
SIMPLE 
QUICK 
Secondorder upwind 
Secondorder upwind 
Firstorder upwind 
According to the related literature, the turbulence intensity and turbulence characteristic length are set at 5% and 0.07 times of the inlet width respectively. The discrete phase loading was considered to be per the unit volume of air with a mean diameter of and a geometric standard deviation of 2.08. Using RosinRammler distribution function, the relationship between droplet diameter and the mass fraction of droplets can be obtained (Fathizadeh et al., 2015).
(17) 
where, and represent the characteristic diameter and distribution parameter respectively. The size of the droplets is selected in the range of to . The droplets are injected in the face normal condition at a velocity equal to that of the flowing fluid at the cyclone inlet, and discrete phase model (DPM) setting is defined as a trap at the bottom of the cyclone. The droplets reach the cyclone walls, collide with them, and reflect off them. Under the low droplet loading condition, it is entirely reasonable to assume that the collision of discrete phase droplets with the cyclone walls is perfectly elastic; hence, the coefficients of restitution (COR) in tangential as well as normal directions of the walls are considered to be unity (Fathizadeh et al., 2015).
The residence time of an unsteady simulation is obtained by the following relation (Parvaz et al., 2017):
(18) 
The calculated residence time of the present simulation is equal to 0.73 s.
5. Grid independence study
Hexagonal cells produced via Gambit software are depicted in Figure 4. To prevent variations in the results with grid size, the static pressure drop was used as a criterion to study the mesh independency. Several grids were generated as listed in Table 4. In addition to pressure, the results of studying grid independence are demonstrated for tangential and axial velocities at Z=0.75D in Figure 5. According to Table 3 and the profile of velocities, a geometry containing 765520 cells with an error of 0.96% was selected.
Figure 4
Grid generation.
Table 4
Details of grid independence study.
Number of cells 
585424 
630952 
765520 
Static pressure drop 
339.487 
355.104 
358.522 
Deviation (%) 
 
4.6 
0.96 

Figure 5
Grid independence study of a) the tangential velocity and b) the axial velocity.
The validity range of the logarithmiclaw for mean velocity is 11.225 < y+ < 300, while the laminar stressstrain relationship is valid for y+ < 11.225. The distance of the first grid center from the wall was checked to lie within the range of the logarithmiclaw (11.225 <y+ < 300), which is the default range in FLEUNT (see Figure 6).
By comparing the obtained results with the experimental data of Hoekstra (2000), the present study is validated. The comparison between the numerical results and experimental data on the axial and tangential velocities and pressure drop is displayed in Figures 7 and 8a. Moreover, the results of the dispersed phase are compared to the experimental data of Ehteram et al. (2012) as shown in Figure 8b. It can be seen that the present study simulates the continuous and discrete phases inside the cyclone at an acceptable degree of accuracy.
Glass 
Steel (wrought iron) 
Galvanized iron 
Figure 6
Contour plots of the wall y+ for a) glass, b) steel (wrought iron), and c) galvanized iron.

Figure 7
The comparison of the profiles of a) tangential velocity and b) axial velocity between the current numerical simulation and the experimental data of Hoekstra (2000) in the section Z = 0.75D.

Figure 8
The comparison of the numerical simulation data on a) pressure drop and b) grade efficiency respectively with the works of Hoekstra (2000) and Ehteram et al. (2012).
6. Results and discussion
6.1. Tangential velocity distribution
The tangential velocity has an Sshaped profile as can be seen in Figure 7a. It is divided into two parts. The first part is a forced vortex flow which is mainly located in the central region of the cyclone, and the second part is a free vortex region far from the central axis of the cyclone. In this region, the dissipation increases while the velocity declines. The sum of these two areas is known as Rankine vortex. The effect of roughness on tangential velocity distribution is depicted in Figure 9. As it is obvious, an increase in velocity results in a rise in pressure, which is originated from the swirling flow of gas and its contraction at the entrance of the outlet tube. According to the tangential velocity distributions in Figure 9, the tangential velocity forms both of the regions. Adding roughness to the cyclone walls leads to a rise in the friction between the gas flow and the wall, which reduces the radial distribution of tangential velocity. Therefore, the velocity distribution of glass is located at a higher level than the distributions of steel (wrought iron) and galvanized iron. Although roughness has no noticeable effect on the radial distribution of tangential velocity, its influence on the erosion rate is really important in industry. It is worth noting that increasing roughness shrinks the Rankine vortex region and subsequently reduces the cyclone performance.

Figure 9
Radial distribution of tangential velocity by considering the wall roughness: a) Z = 0.75D and b) Z = 0.625D.
6.2. Axial velocity distribution
The axial velocity has an essential role in the separation of droplets. According to Figure 7b and experimental data, the axial velocity profile is Mshaped, and an outer vortex flow is created due to descending gas flow. Following that, the inner vortex flow is formed after the first vortex flow impacts on the bottom of the cyclone, which transfers small droplets through the vortex finder.
Axial velocity influences both the flow pattern and cyclone performance. The axial velocity distributions at different roughness values are illustrated in Figure 10 at Z = 0.75D and Z = 0.625D. The gas flow has a downward trend near the walls, while it experiences an upward trend in the central part of the cyclone. This phenomenon is known as a reverse flow. Since the outer downward flow rises with an increase in radius, but the upward inner flow increases by a decrease in radius, the maximum axial velocity happens at the central axis of the cyclone. Therefore, the wall roughness has a more significant effect on the upward flow than the downward flow.

Figure 10
Radial distribution of axial velocity by considering the wall roughness: a) Z = 0.75D and b) Z = 0.625D.
6.3. Pressure distribution
The impact of surface roughness on the pressure is also studied herein. The static pressure contours and local static pressure distribution along the cyclone length are presented in Figures 11 and 12. As can be seen, an increase in the roughness leads to a lower pressure drop at the cyclone walls. In other words, the decreased tangential and axial velocities resulted from an increase in the roughness causes the velocity gradients to drop. As a result, it can be expected that the pressure must be reduced, which is supported by the data depicted in Figures 11 and 12.
Pressure (Pa) 

Pressure (Pa) 

Pressure (Pa) 




Galvanized iron 


Steel (wrought iron) 


Glass 


Figure 11
Static pressure at different roughness values.
Figure 12
The impact of roughness on static pressure.
The pressure drop is determined as the difference between the static pressure at the inlet and the one at outlet. A higher pressure drop represents higher energy loss during the separation process of droplets. Hence, the decrease of pressure drop is an essential parameter in designing cyclones. According to the results shown in Figure 13a, increasing the roughness and inlet velocity both decrease the pressure drop, so the pressure drop in a cyclone made of glass is much more than that in a cyclone made of galvanized iron.


Figure 13
The a) pressure drop and b) pressure drop coefficient (or Euler number) at different roughness values and velocities.
On the other hand, it can be seen that the tangential velocity declines at higher surface roughness. Moreover, the energy loss between the dispersed and continuous phases rises, so the pressure drop increases. The variation in the pressure drop because of changing wall roughness broadens at higher velocities. In this situation, the pressure drop coefficient (or Euler number) is defined as:
(19) 
The pressure drop coefficient at various roughness values is delineated in Figure 13b. Increased wall roughness results in decreasing the Euler number at every Reynolds number. Furthermore, it can be inferred from Figure 13b that, for a glass cyclone, an increase in Reynolds number raises Euler number; nonetheless, for a steel (wrought iron) and galvanized iron cyclone, the trend is inverted, that is, an increase in Reynolds number leads to a decrease in pressure drop coefficient.
6.4. Vortex length and the cyclone performance
The vortex length is defined as the axial distance between the outlet tube and the ending vortex flow. It can be seen clearly that the vortex length given in Equations 14 to 16 is defined by three dimensionless geometrical parameters which are not related to the wall roughness. While in the present study, the impact of surface roughness on the vortex length is studied. It can be concluded from Figure 11 that the severe changes in pressure at different wall locations result in a lower vortex length, which affects the cyclone performance from two points of view. It leads to a reduction in pressure drop, which is preferred; however, it influences the separation efficiency reversely. Figure 14 reveals that vortex length decreases as the surface roughness rises.
Figure 14
The impact of surface roughness on vortex length.
In fact, the surface roughness leads to friction between gas and the cyclone wall and decreases the gas swirling, which reduces the vortex length and the cyclone efficiency. According to Figure 15, the cyclone separation efficiency improves when it is made of smooth materials like glass rather than rough materials like galvanized steel. Since increasing the surface roughness causes tangential velocity to decline, it may be concluded that increased surface roughness of the cyclone walls results in lower collection efficiency. Nevertheless, the efficiency of the cyclone enhances at higher velocities.

Figure 15
The efficiency of cyclone at different roughness values at a velocity of a) 10 m/s and b) 16 m/s.
7. Conclusions
The wide range of applications of cyclones in industry has inspired many researchers to study the effects of several parameters on the cyclone performance. In the present study, the effects of three values of wall roughness alongside different droplet sizes and various gas inlet velocities on the flow parameters of cyclones are studied. Hence, the radial distribution of axial and tangential velocities and contour plots of pressure distribution are analyzed and presented. Moreover, vortex length and cyclone separation efficiency are investigated numerically. Also, the splashing phenomenon, as a controversial topic among researchers, is considered. The RSM model and EulerianLagrangian approach are employed to simulate the turbulent flow and droplet motion. Finally, based on the results obtained, the following conclusions can be drawn:
 An increase in the surface roughness of the cyclone wall leads to a decline in tangential velocity as the flow resistance rises.
 The surface roughness affects the pressure drop significantly at a higher inlet velocity, while its impact on pressure drop is negligible at lower inlet velocities.
 An increase in the surface roughness of the cyclone wall leads to a decline in swirling flow, the tangential and axial velocities, the pressure drop, and Euler number.
 An increased surface roughness reduces the vortex length and the cyclone efficiency.
 Splashing phenomenon does not happen on the walls of the cyclone.
Nomenclature
Drag coefficient 

Friction factor 

Gravitational acceleration 

turbulence kinetics energy 

Mean pressure 

Stress generation 

Volumetric flow rate 

Reynolds stress tensor 

Droplet resident time 

xcomponent of velocity 

Fluctuating component of velocity 

Mean velocity 

ycomponent of velocity 

Cyclone volume 

zcomponent of velocity 

Location 

Greek Symbols 

Kinematic viscosity 

Turbulent viscosity 

Density 

Turbulence dissipation 

Dynamic viscosity 

Subscripts 

Droplet 

Turbulence 