Challenge

Week 8 - Simulating Cyclone separator with Discrete Phase Modelling

Aim :

1) To perform analysis on cyclone separator and calculate the separation efficiency and pressure drop.

2) To write a few words about any four empirical models used to calculate the cyclone separator efficiency.

3)To perform an analysis on a given cyclone separator model by varying the particle diameter from 1 μm to 5 μm and calculate the separation efficiency in each case and discuss the results.

4)To perform an analysis on a given cyclone separator model by varying the particle velocity from 1 m/sec to 5 m/sec and calculate the separation efficiency and pressure drop in each case and discuss the results.

Theory :

Cyclone separators or simply cyclones are separation devices that use the principle of inertia to remove particulate matter from flue gases. Cyclone separators is one of many air pollution control devices known as precleaners since they generally remove larger pieces of particulate matter. This prevents finer filtration methods from having to deal with large, more abrasive particles later on. In addition, several cyclone separators can operate in parallel, and this system is known as a multicyclone.

It is important to note that cyclones can vary drastically in their size. The size of the cyclone depends largely on how much flue gas must be filtered, thus larger operations tend to need larger cyclones. For example, several different models of one cyclone type can exist, and the sizes can range from a relatively small 1.2-1.5 m tall to around 9 meters .

Cyclone separators work much like a centrifuge, but with a continuous feed of dirty air. In a cyclone separator, dirty flue gas is fed into a chamber. The inside of the chamber creates a spiral vortex, similar to a tornado. This spiral formation and the separation is shown in Figure 2. The lighter components of this gas have less inertia, so it is easier for them to be influenced by the vortex and travel up it. Contrarily, larger components of particulate matter have more inertia and are not as easily influenced by the vortex.

Since these larger particles have difficulty following the high-speed spiral motion of the gas and the vortex, the particles hit the inside walls of the container and drop down into a collection hopper. These chambers are shaped like an upside-down cone to promote the collection of these particles at the bottom of the container. The cleaned flue gas escapes out the top of the chamber.

Most cyclones are built to control and remove particulate matter that is larger than 10 micrometers in diameter. However, there do exist high efficiency cyclones that are designed to be effective on particles as small as 2.5 micrometers. As well, these separators are not effective on extremely large particulate matter. For particulates around 200 micrometers in size, gravity settling chambers or momentum separators are a better option.

Out of all of the particulate-control devices, cyclone separators are among the least expensive.They are often used as a pre-treatment before the flue gas enters more effective pollution control devices. Therefore, cyclone separators can be seen as "rough separators" before the flue gas reaches the fine filtration stages.

1) Emperical methods used to calculate cyclone seperator efficiency .

a) Iozia and Leith Model

Iozia and Leith [4] logistic model is a modified version of Barth Model, which is developed based on force balance.
The model assumes that a particle carried by the vortex endures 165 the influence of two forces: a centrifugal force, Z and a flow resistance, W . The collection efficiency `eta_i` of particle diameter `d_(p)` can be calculated from
`eta_i =1/(1+(d_(pc)+d_p)^beta`
`beta` is an expression for slope parameter derived based on the statistical analysis of experimental data of a cyclone with D = 0.25m given as
`beta = 0.62 - 0.87ln(d_(pc)/100) + 5.21ln((ab)/D^2) +1.05 (ln((ab)/D^2))^2`
b) Li and Wang model
The Li and Wang model includes particle bounce or re-entrainment and turbulent diffusion at the cyclone wall. A two dimensional analytical expression of particle distribution in the cyclone is obtained. Li and Wang model was developed based
on the following assumptions
i) The radial particle velocity and the radial concentration profile are not constant for uncollected particles within the cyclone
ii) Boundary conditions with the consideration of turbulent diffusion coefficient and particle bounce re-entrainment on the cyclone wall are:
c=`c=c_0 at theta=0`
`D_r (delc)/(delr) = (1-alpha)omegac` at r=D/2
iii) The tangential velocity is related to the radius of cyclone by `uR^n=`constant
The resultant expression for collection efficiency of particle of any size is given by
`eta_i = 1-exp{-lambda theta_1}`
where `theta_1 = 2pi(S+L)/a`
c) Koch and Licht Model
Koch and Licht [2] collection theory recognized the inherently turbulent nature of cyclones and the distribution of gas
residence times within the cyclone. Koch and Licht described particle motion in the entry and collection regions with the additional following assumptions
i) The tangential velocity of a particle is equal to the tangential velocity of the gas flow, i.e. there is no slip in the tangential direction between the particle and the gas .
ii) Tangential velocity is related to radius of cyclone by `uR^n =` constant
A force balance and an equation on the particles collection yields the grade efficiency `eta_i`
d) Lapple model
Lapple model was developed based on force balance without considering the flow resistance. Lapple assumed that a particle entering the cyclone is evenly distributed across the inlet opening. The particle that travels from inlet half width to the
wall in the cyclone is collected with 50% efficiency. The semi empirical relationship developed by Lapple to calculate a
50% cut diamter `d_(pc)`
efficiency of collection of any size of particle is given by
Geometry
The STEP file of the model was imported to space claim and required volume was extracted using extract with edges tool .
Extracted Volume
Geometry clean up was thus completed and meshing module was opened .
Mesh
The cartesian meshing method was used to obtain hexa mesh for the model . The element size was set to 5 mm and method element size was set to 2.5 mm .
For case 3,4 & 5 , another type of meshing was used to get a steady solution. Hexa dominant meshing was used with element size of 10mm along with face sizing with element lenth of 5mm .
Named selection was used to name various surfaces
Thus meshing was completed and set up module was opened
SETUP

The general solver setting was set to pressure based with absolute velocity formulation and steady state solution .Gravity was set to 9.81 along -y direction . The working fluid was selected as air and material for seperator was chosen as Aluminium.

The viscous model was set to RNG k-epsilon with swirl dominated flow and standard wall functions.

The Discrete Phase Modelling option was chosen with Anthracite as material , particle of size 1micrometer was chosen and was injected along inlet surface with velocity 3m/s .

Boundary Conditions

Inlet :

Velocity =3 m/s

DPM phase type = Reflect

Outlet heavy :

DPM phase type = Trap

Outlet light :

DPM phase type = Escape

Wall :

DPM phase type = Reflect

The under relaxation factors for solution controls were altered and momentum , turbulent kinetic energy and turbulent dissipation rate discretization was changed to First order upwind .

Report Definitions were requested for area weighted pressure calculation at inlet and outlet .

The solution was initialized using standard option with reference as inlet .

The solution was run for 400 Iterations

Results :

Case 1:

Particle size : 1 `mum`

inlet velocity : 3 m/s

Scaled Residuals

Particle Time Contour

Vortex Region

Pressure Contour

Case 2:

Particle size : 5 `mum`

inlet velocity : 3 m/s

Scaled Residuals

Particle Time Contour

Vortex Region

Pressure Contour

Case 3:

Particle size : 5 `mum`

inlet velocity : 1 m/s

particle velocity : 1m/s

Scaled Residuals

Particle Time Contour

Vortex Region

Pressure Contour

Case 4:

Particle size : 5 `mum`

inlet velocity : 3 m/s

particle velocity : 1m/s

Scaled Residuals

Particle Time Contour

Vortex Region

Pressure Contour

Case 5:

Particle size : 5 `mum`

inlet velocity : 5 m/s

particle velocity : 1m/s

Scaled Residuals

Particle Time Contour

Vortex Region

Pressure Contour

Case no Inlet Velocity (m/s) Particle Size( `mum`) Pressure drop (Pa)

Efiiciency in %

Trapped/Tracked

1 3 1 5.74 33 .59
2 3 5 5.746 38.09
3 1 5 0.653 15.55
4 3 5 5.46 18.5
5 5 5 14.92 20.5

Inference

From the table , we can conclude that

1) As particle size increases , efficiency increases for a given inlet velocity .

2) As inlet velocity increases, pressure drop increases as well as efficiency increase for a given particle size.

Your portfolio, next

Ship work like this
in your own course.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.