The Effect of The Introduction of An Exit Tube On The Separation Efficiency in A Cyclone
The Effect of The Introduction of An Exit Tube On The Separation Efficiency in A Cyclone
Vol. 30, No. 03, pp. 627 - 641, July - September, 2013
(Submitted: February 22, 2012 ; Revised: August 14, 2012 ; Accepted: August 25, 2012)
Abstract - The paper shows the analysis of the effect of introduction of an exit tube into the cylindrical part
of a cyclone on the flow velocity, pressure losses and, above all, the separation efficiency. We made
measurements and carried out CFD simulations for three levels of the exit tube introduction. The ratios of the
depth of the exit tube introduction to the cyclone diameter (Hp/D) were 0.4, 0.475 and 0.89 at three velocity
settings of 8 m.s-1, 13 m.s-1 and 15 m.s-1. The Reynolds Stress Model (RSM), which was compatible with the
experimental results, was used for numerical solutions. The efficiency of cyclone separation was explored on
a sample of oak sawdust. The efficiency of cyclone separation increased with the geometric size of the
particles, the inlet velocity and a deeper introduction of the exit tube.
Keywords: Cyclone; Separation; Computational fluid dynamics; Oak sawdust.
The works of Lapple (1951) and Barth (1956) are A literature survey showed that the vortex length
considered to be the first studies dealing with cyclone can be important for predicting the separation
models. The studies were simple and provided efficiency, especially in short cyclones, but for
acceptable results with a limited range of parameters. longer cyclones the turbulence does not have to
Later, those theories were improved by Smith (1962) reach the cone top. If the length of a cyclonic
and Muschelknautz (1970). In their theories they cylinder increases, the separation efficiency in-
considered more parameters such as surface creases to a certain value and then it begins to
roughness, area of distribution of particles, shape and decrease; it is also similar in an overflow pipe. By
size of particles, pressure losses and others. The lengthening the height of the cyclonic cylinder or by
efficiency of vortex separators depends on the shortening the length of the overflow pipe we can
regime, which is also influenced, apart from the lessen the pressure drop to a certain extent (Zhu and
already mentioned parameters, by inlet velocity, Lee, 1999). The definition and composition of the
temperature and viscosity, which are introduced in pressure drop in the cyclone were analysed in the
the work of Avci and Karagoz (2003). This work paper of Chen and Shi (2006). The pressure drop in
investigated the cyclone flow patterns introducing the cyclone includes spreading losses at the inlet into
different shapes and diameters of exit tube simulated the cyclone and outlet spinning losses caused by the
at diverse flow rates. Lim et al. (2004) examined friction near the walls of the cyclone and the
experimentally the effect of vortex finder shape on dynamic energy of the diffused gas in the overflow
the characteristics of the collection efficiency at pipe. The surface resistance or losses caused by
different flow rates. Elsayed et al. (2010) evaluated friction have a substantial influence on the vortex
the flow patterns and performance of a cyclone with length; with decreasing surface resistivity the length
different vortex finder diameters to examine the of turbulence in the cyclone increases. Numerical
effect of the vortex finder diameter on the studies of separable properties of different concentra-
performance and velocity profile. In the experimen- tions of particles in the cyclone at the inlet were
tal work of Hoffman et al. (1995) it was proved that investigated in the paper of Oian et al. (2007). The
the inlet velocity has a significant influence on the effects of inlet concentration of particles on the
length of turbulence depending on cyclone geometric tangential velocity, particle size, pressure drop and
dimensions. The paper of Lidén and Gudmundsson separation efficiency in the cyclone were deter-
(1996) described the use of semi-empirical polyhedra mined. The impacts of cone dimensions on the
for separation efficiency depending on the operating cyclone performance were examined in the study of
conditions and cyclone dimensions. The length of Xiang et al. (2000).
turbulence in the cyclone increased with increasing The mathematical expression for three-dimen-
inlet velocity. In the work of Martignoni et al. sional flow in the cyclone is complex and has not
(2007), the impact of cyclone geometry was studied been explored in depth. The calculations of physical
by means of creating a symmetrical double inlet and laws are dependent on the attention, time and effort.
a cyclone model with a spiral outlet compared with a Therefore, it was necessary to examine these theories
common cyclone with a tangential inlet. It was on the computer by means of polyhedra and to
verified that the total separation efficiency increases compare the results obtained with experimental data
with decreasing pressure drop in both cyclone or theoretical knowledge. For a better idea of flow
shapes, but in a cyclone with shaped symmetrical and physical actions occurring in the cyclone, CFD
inlets a bigger pressure drop was recorded than in the (Computational Fluid Dynamics) models are utilized
cyclone with the spiral outlet. The results of the at present. The results of numerical calculations and
experimental study of Kim, Lee and Kuhlman (2001) gas flow simulations in the typical Lapple cyclone
provide information about the efficiency of particle are described in the paper of Wang et al. (2006). The
separation for a modified inner surface of the body accuracy of the numerical solution of a flow
of cyclones with spiral, circumferential and vertical simulation in the cyclone was achieved using the
grooves. The estimation of separation efficiency RSM model. The CFD flow simulations completed
dependent on the concentration of particles in axial with properties of heat transfer in the cyclone with a
flow in a cyclone dust collector was published by tangential inlet were published by Karagoz and Kaya
Ogawa (1999). Under the centrifugal force solid (2007). The publication of Griffiths and Boysan
particles slide along the cone surface and fall into a (1996) focused on the comparison of performance of
dust hopper, where some particles are trapped and three small cyclone types with the use of CFD data
some are blown away by secondary flow back to the from experiments and three empirical theories. Gas
cyclone body. flow and a solid component in the cyclone are
described in the paper of Meier and Mori (1998). In where ρ is the density, u is the mean velocity, x is the
this work, the CFD model was based on an Eulerian position, t is the time, p is the pressure, μ is the
approach for both phases with the condition of a 3D dynamic viscosity, u′i u′j is the Reynolds stress tensor,
symmetrical domain, the mathematical model was
compiled for the utilization of the k-ε turbulent g is the acceleration of gravity and Fi is the external
model. body force which arises due to the interaction with
The possibility of using CFD techniques enables dispersed solid particles.
one to make a proposal of more effective separation This equation of motion now represents the
conditions without the need of conducting expensive averaged velocity values (or main flow). Turbulence
experiments. Currently, CFD simulations have rela- is represented by means of the “Reynolds stresses”
_____
tively wide application for relatively precise mathe- ρ u′i u′j . FLUENT describes the effect of main flow
matical approximation of flow, while the arduousness
of its performance and time restraints in solving and the Reynolds stresses by means of the above
concrete physical problems must not be forgotten. mentioned models. The Reynolds Stress Model
includes the calculation of particular Reynolds
stresses by means of a transport differential equation
TURBULENT MODEL DESCRIPTION in the form:
The use of CFD (computational fluid dynamics) ∂ u′i u′j ∂ u′i u′j
+ uk =
for numerical calculations of gas flow in the cyclone ∂t ∂x k
has increasingly been applied. Three models are ( LTM ) ( )
Cij
usually used for the simulation in the cyclone: the
k-ε model, Algebraic Stress Model (ASM) and ∂ ⎡ p′ ∂ ⎤
Reynolds Stress Model (RSM). Wang et al. (2006) −
∂x k ⎣
(
⎢ u′i u′ju′k +
ρ
) ( )
δkju′i + δki u′j − υ
∂x k
( )
u′i u′j ⎥
found that the standard k-ε turbulence model is ⎦
inadequate to simulate the flow with swirling motion Dij ( ) (3)
because it leads to excessive turbulence viscosity and
unrealistic tangential velocity. Therefore, the accu- ⎡ ∂u j ∂u ⎤ p′ ⎡ ∂u ∂u j ⎤
− ⎢ u′i u′k + u′i u′k i ⎥ + ⎢ i + ⎥
racy of the numerical solution can be improved by ⎣ ∂x k ∂x k ⎦ ρ ⎢⎣ ∂x j ∂x i ⎥⎦
using the Reynolds Stress Model (RSM), which was
also utilized in this paper for the simulation of
( )
Pij Φ ij( )
physical actions in the cyclone by using the ⎡ ∂u ∂u i ⎤
FLUENT commercial program. −2 υ ⎢ i ⎥ + Su ,
In turbulent flow, which is accompanied by ⎣ ∂x k ∂x k ⎦
pulsation and fluctuation, the total velocity equals ( )
εij
the sum of the central and fluctuation components of
velocity: where the left two terms are the local time derivative
of stress (LTM) and the convective transport term
u i = ui + u i ' . (1) (Cij), respectively. The right five terms are: Dij is the
stress diffusion term, Pij is the shear production term,
By substituting this expression into the equation Фij is the pressure- strain term, εij is the dissipation
of motion we get the averaged Navier-Stokes term and Su is the source term.
equations of motion in the form: For better expression, the individual members are
approximated for the purpose of closing the set of
∂ ∂ simultaneous equations from this differential equa-
∂t
( ρu i ) +
∂x j
(
ρu i u j ) tion. The Reynolds stresses are consequently substi-
tuted into the equation of momentum transfer.
Generally, the following equations are solved by
∂ ⎛ ⎡ ∂u i ∂u j ⎤ ⎛ 2 ∂u j ⎞ ⎞ using the RSM model:
= ⎜μ ⎢ + ⎥−⎜ μ ⎟⎟ (2)
∂x j ⎜⎝ ⎣⎢ ∂x j ∂x i ⎥⎦ ⎝⎜ 3 ∂x j ⎠⎟ ⎟⎠ 6 transport equations in the form (3),
3 transport equations for the mean component of
∂p ∂
−
∂x i
+ ρg + Fi +
∂x j
(
−ρ u′i u′j ,) velocity,
continuity equation,
Brazilian Journal of Chemical Engineering Vol. 30, No. 03, pp. 627 - 641, July - September, 2013
630 J. Cernecky and K. Plandorova
du iP 3νC D Re p
dt
=
4d 2S
( )
u i − u iP + gi , (4)
dx i
u iP = . (5)
dt
where ReP is the particle Reynolds number defined A set of measurements was carried out for three
as: positions of introduction of an exit tube into the
inner part of the cyclone at a barometric pressure of
d u j − u Pj 1 025 hPa. The exit tube is in the shape of a cylinder
Re P = . (8) with a length of 300 mm and was designed in a way
ν so that it is possible to set the necessary depth of
introduction of the tube in the cyclone. The first
measurements were made at a ratio of Hp/D = 0.89,
CONDITIONS OF THE EXPERIMENT AND Figure 2 (a). The second measurements were performed
CFD SIMULATION at a ratio of Hp/D = 0.475, just below a bottom wall
of the inlet cyclone pipeline. The third measurements
Figure 1 shows a designed cyclone model with a were carried out at a ratio of Hp/D = 0.4, Figure 2 (c).
tangential inlet and Table 1 shows basic cyclone A five-port dynamic probe clamped in a special
dimensions. The cyclone model is placed in an holder was used for these measurements. The probe
experimental overpressure system for the loose was calibrated on a calibration device before the
matter separation of diverse composition and measurements (Polansky and Stech, 2011). A sample
concentration. This device was used for experimental of oak sawdust with a density of 670 kg.m-3 and
purposes with an inlet concentration of particles of dampness of 6.5% was used for finding out the
0.01 kg. N-1.m-3. In terms of air flow it is an air- cyclone separation efficiency.
technical system where the transport air is drawn These experiments were carried out to verify the
from the room to a pipe system and from the numerical model. The inlet air velocity as well as the
cyclone separator transported back to the ambient particle velocity were set with a control damper on
environment. the intake opening of the fan at the values of 8 m.s-1,
13 m.s-1 and 15 m.s-1. Three types of CFD models it is 92,606 CFD cells and for Hp/D = 0.89 it is
were created. They represented three different 111,661 CFD cells. Due the use of the polyhedral
examples of the introduction of an exit tube into the grid, which is of high quality, the grid optimization
cylindrical part of the cyclone. Surfaces were was not done during the computation.
covered with a hexagonal grid of tetra-hybrid type, In order to calculate the particle trajectories the
Figure 3, which was converted to a grid of the poly- following parameters of DPM (Discrete phase
hedral type in the FLUENT program. A computa- modelling) were set: the total flow rate – 0.0029 kg/s,
tional domain for the variant of Hp/D = 0.4 comprises the particle diameter – 1 mm, uniform particle distri-
79,393 CFD cells, for the variant of Hp/D = 0.475 bution and particle material oak.
Brazilian Journal of Chemical Engineering Vol. 30, No. 03, pp. 627 - 641, July - September, 2013
632 J. Cernecky and K. Plandorova
Figure 5: A comparison of numerical and calculated Figure 6: A comparison of numerical and calculated
dynamic pressure. tangential velocities.
of Hp/D = 0.4, 0.475 and 0.89. Figures 7 – 9 show vortex along the cyclone axis and in a closed
the contours of tangential velocity for the introduc- container.
tion of an exit tube at Hp/D = 0.4, 0.475 and 0.89 In all three cases of the introduction of an exit
and for three settings of the inlet velocity of 8 m.s-1, tube, the vortex core is in the shape of a twisted
13 m.s-1 a 15 m.s-1. From the figures it is obvious cylinder and is not axially completely symmetrical,
that, with increasing inlet velocity, the vortex core in especially in the conical part of the cyclone towards
the whole cyclone section is changed. The dominant the container of trapped particles. The axis of the
velocity in the cyclone is a tangential velocity forced vortex is not identical in time with the
component. The value of the tangential velocity geometric cyclone axis, which is in the shape of
equals to zero on the wall and in the centre of field a curve. With increasing inlet velocity the length of
flow. The highest tangential velocity is achieved by turbulence in the cyclone increases.
sucking from the inlet pipeline; then the velocity is Figure 10 shows the contours of tangential velocity
decreased by the vortex gas motion downwards for the introduction at three inlet velocities. From the
along the cyclone wall. Negative values of the figures it can be seen that, with increasing inlet
tangential velocity are achieved in a created, forced velocity, the tangential velocity rises proportionally.
Brazilian Journal of Chemical Engineering Vol. 30, No. 03, pp. 627 - 641, July - September, 2013
634 J. Cernecky and K. Plandorova
Figure 10: Contours of tangential velocity for the ratios of Hp/D = 0.4, 0.475 and 0.89.
Figures 11 – 13 show the contours of axial velocity sample of oak sawdust with a density of 670 kg. m-3
for Hp/D = 0.89 and at the inlet velocity of 8 m.s-1, at three velocities and at three settings of an exit
13 m.s-1 and 15 m.s-1. In the figures it is possible to tube, Figure 15. The particle sizes of the sample used
see that the axial velocity near the cyclone wall were 1 mm, 0.6 mm and 0.125 mm and the separation
achieves negative values and towards the cyclone period was 60 s. From the measurement results it can
centre these velocity values turn into positive ones. be said that the separation efficiency decreased with
Higher values of the axial velocity were recorded in the size of the particles. The next important factor by
the centre of the cylindrical part than in the conical which the separation efficiency can be influenced is
cyclone part. The cyclone hopper has negative values separation velocity and the introduction of an exit
of axial velocity in the whole cross-section. The tube. With increasing inlet velocity the separation
highest axial velocity is achieved near the wall of an efficiency increased for all three settings of an exit
exit tube. tube. At Hp/D = 0.4 the separation was lower because
Figure 14 shows the contours of axial velocity for the particles were dragged from the trajectory of
the introduction at three inlet velocities. swirling motion into the exit tube where they
The separation efficiency was investigated on a escaped from the cyclone as light ashes.
Brazilian Journal of Chemical Engineering Vol. 30, No. 03, pp. 627 - 641, July - September, 2013
636 J. Cernecky and K. Plandorova
Figure 14: Contours of axial velocity for the ratios of Hp/D = 0.4, 0.475 and 0.89.
Figure 15: The efficiency of cyclone separation for oak sawdust at the ratios of Hp/D = 0.4, 0.475 and 0.89.
In Figure 16 there is a simulation of the trajectory velocity at the inlet into the cyclone and the
of oak particles with a size of 1 mm coloured downwards along the conical cyclone part the
according to the velocity size for Hp/D = 0.475 at velocity decreases. The fact that the velocity
the inlet velocity of 15 m.s-1; the time step was 0.2 s decreases is due to the abrasion of particles on the
and the time stop was 10 s. From the figure it is cyclone wall, the action of gravity and also the
possible to see that the particles achieve the highest weight of newly arriving particles in the cyclone.
Brazilian Journal of Chemical Engineering Vol. 30, No. 03, pp. 627 - 641, July - September, 2013
638 J. Cernecky and K. Plandorova
The following figure (Figure 17) shows a graphical for the fractions of 1 mm, 0.6 mm and 0.125 mm, parti-
dependence of the fractional separation by particle size cle sizes that occur mainly in the wood-processing indus-
for Hp/D = 0.475 for the flow velocities of 8, 13 and try. A comparison of the experimental results with the
15 m.s-1. The separation was investigated experimentally outputs of CFD simulations can be seen in the Figure 18.
Figure 18: The efficiency of cyclone separation for diverse flow velocities, particle sizes and ratios of Hp/D.
flow cyclone dust collector. Journal of Thermal Smith, J. L., An experimental study of the vortex in
Science, 8, No. 3, 143 (1999). the cyclone separator. Journal of Basic
Polansky, J. and Stech, J., Calibration of a five-port Engineering, 84, No. 4, 602 (1962).
probe – in Czech. ZCU Plzen (2011). Wang, B., Xu, D. L., Chu, K. W. and Yu, A. B.,
Qian, F. and Wu, Y., Effects of the inlet section Numerical study of gas-solid flow in a cyclone
angle on the separation performance of a cyclone. separator. Applied Mathematical Modelling, 30,
Chemical Engineering Research and Design, 87, No. 11, 1326 (2006).
No. 12, 1567 (2009). Xiang, R., Park, S. H. and Lee, K. W., Effects of
Raoufi, A., Shams, M., Farzaneh, M., and Ebrahimi, cone dimension on cyclone performance. Journal
R., Numerical simulation and optimization of of Aerosol Science, 32, No. 4, 549 (2001).
fluid flow in cyclone vortex finder. Chemical Zhu, I. and Lee, K. W., Experimental study on small
Engineering and Processing, 47, No. 1, 128 cyclones operating at high flowrates. Journal of
(2008). Aerosol Science, 30, No. 10, 1303 (1999).
Brazilian Journal of Chemical Engineering Vol. 30, No. 03, pp. 627 - 641, July - September, 2013