Week 8 - Simulating Cyclone separator with Discrete Phase Modelling

Aim:- To perform analysis on cyclone separator and calculate the separation efficiency and pressure drop.

Objective:-

1. To write a few words about any four empirical models used to calculate the cyclone separator efficiency.
2. 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. Discuss the results. [Use both the velocity's as 3m/sec.]
3. 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. Discuss the results. [Use particle diameter size as 5 μm for all cases & keep flow velocity same as particle velocity]

Boundary Conditions:-

Inlet and wall:- Reflect

Outlet(Top):- Escape

Outlet(Bottom/dustbin) - Trap

What is cyclone separator?

• Cyclonic separation is a method of removing particulates from an air, gas or liquid stream, without the use of filters, through vortex separation.
• When removing particulate matter from liquid, a hydrocyclone is used; while from gas, a gas cyclone is used.
• Rotational effects and gravity are used to separate mixtures of solids and fluids. The method can also be used to separate fine droplets of liquid from a gaseous stream.

Working of a Cyclone Seperator :

• A high speed rotating (air)flow is established within a cylindrical or conical container called a cyclone.
• Air flows in a helical pattern, beginning at the top (wide end) of the cyclone and ending at the bottom (narrow) end before exiting the cyclone in a straight stream through the center of the cyclone and out the top.
• Larger (denser) particles in the rotating stream have too much inertia to follow the tight curve of the stream, and thus strike the outside wall, then fall to the bottom of the cyclone where they can be removed.
• In a conical system, as the rotating flow moves towards the narrow end of the cyclone, the rotational radius of the stream is reduced, thus separating smaller and smaller particles.
• The cyclone geometry, together with volumetric flow rate, defines the cut point of the cyclone.
• This is the size of particle that will be removed from the stream with a 50% efficiency.
• Particles larger than the cut point will be removed with a greater efficiency, and smaller particles with a lower efficiency as they separate with more difficulty or can be subject to re-entrainment when the air vortex reverses direction to move in direction of the outlet.

Empirical Models:

1. Iozia and Leith Model:

• This model is based on force balance.
• The model assumes that a particle carried by the vortex endures the influence of two forces: a centrifugal force, Z and a flow resistance, W.
• The collection efficiency of a particle diameter dpi can be calculated as,

• β is an expression for slope parameter derived based on thestatistical analysis of experimental data of a cyclone with D = 0.25 m given as,

• dpc is given by,

• where core length, zc, and core diameter, dc, are given as,

• Values of H, S, B etc can be taken from diagram,

2. Li and Wang Model:

• This 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:

• The tangential velocity is related to the radius of cyclone by: = constant.
• The concenntration distribution in cyclone is given as:

• The resultant expression of the collection efficiency for particle of any size is given as:

3. Koch and Licht Model:

• This 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:

-- 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.

-- The tangential velocity is related to the radius of cyclone by: = constant.

• A force balance and an equation on the particles collection yields the grade efficiency,

• G is a factor related to the configuration of the cyclone, n isrelated to the vortex and τ is the relaxation term.

4. Lapple Model:

• Lapple [1] 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 diameter, dpc, is:

• where Ne is the number of revolutions,

• The efficiency of collection of any size of particle is given by,

Simulation Procedure :

The simmulation was carried out in two parts each having 5 cases.

Part 1: With varying particle diameter.

Part 2: With varying the velocity.

Creating Geometry:

• Load the provided model into spaceclaim.
• Go to volume extract option and select edge selection.
• Select the inlet edge, outlet edge (top and bottom both) and hit the green tick mark to create the volume.
• Make sure to supress the solid body and keep just the extracked volume.
• The geometry is ready for meshing.
• The above procedure for geometry creation is same for all the cases.

Meshing:

• Set the global element size to 10 mm. (free choice)
• Right click on mesh and select 'methods' option.
• Select 'Cartesian' for methods and generate the mesh.
• This is the cutcell method which generates hexahedral mesh.

Setting up the Physics:

• Viscous Model: k-epsilon RNG Swirl Dominanted Flow, this gives a better approximaiton for flows with swirls.

• Select Discrete Phase model.
• Select 'Interaction with Continuous Phase'. This allows the interaction between the continuous phase and the particles.
• Continuous phase is the fluid and the discrete phase is the particle to be separated (in this case).

• Select 'Injections'.
• Set 'Injection type' to surface and select inlet.
• Assign 3 m/s in x-direction.
• One can use the 'Inject using Face Normal Direction', this assures that the velocity is normal to the surface selected.

Boundary Conditions:-

Inlet and wall:- Reflect

Outlet(Top):- Escape

Outlet(Bottom/dustbin) - Trap

• Navigate to 'Methods' and set the following parameters to:

• Conduct a Standard Initialization.
• Compute from inlet.
• Give the number of iteraion and hit Run.

PART 1 : To perform an analysis on a given cyclone separator model by varying the particle diameter from 1 μm to 5 μm

Keeping the velocity constant as 3 m/s.

• Part 1 , Case 1 = 1 μm
• Part 1 , Case 2 = 2 μm
• Part 1 , Case 3 = 3 μm
• Part 1 , Case 4 = 4 μm
• Part 1 , Case 5 = 5 μm

Results for part 1 :

Formula : Separation efficiency = (particles trapped/particles tracked)*100.

CASE 1 :1μm

Residual plot :

• Seperation efficiency = (61/90)*100 = 67.77 %.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

Case 2 : 2μm

Residual plot :

• Separation effiency : (62/91)*100 = 68.131%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

Case 3 : 3μm

Residual plot :

• Separation effiency = (62/91)*100 = 68.131%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

Case 4 : 4μm

Residual plot :

• Separation efficiency = (64/91)*100 = 70.32%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

Case 5 : 5μm

Residual plot :

• Separation efficiency = (59/91)*100 = 64.83%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

Conclusion:

• Seperation efficiency is higher for 4 micrometer.
• The results do not show consistency with respect to change in diameter.
• Therefore, furthur mesh refinement is necessary.

Part 2 : Varying the inlet velocity from 1m/s to 5m/s and keeping particle diameter constant.

• Everything is similar to Part 1, the only change is in the setup of discrete phase change tab.
• Instead of varying the particle diameter, we will vary the x-velocity.
• Note that in boundaries the velocity must be same as particle velocity. ie. x-velocity.

Keeping the diameter constant as 5e-06.

• Part 2, Case 1 = 1 m/s
• Part 2, Case 2 = 2 m/s
• Part 2, Case 3 = 3 m/s
• Part 2, Case 4 = 4 m/s
• Part 2, Case 5 = 5 m/s

Results for part 2 :

Formula = Separation Efficiency = (Particle Traped/Particle Tracked)*100.

case 1 : 1 m/s.

Residual plot :

• Separatipn Efficiency = (68/91)*100 = 74%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

case 2 : 2 m/s.

Residual plot :

• Separatipn Efficiency = (66/91)*100 = 72.52%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

case 3 : 3 m/s.

Residual plot :

• Separatipn Efficiency = (59/91)*100 = 64.83%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

case 4 : 4 m/s.

Residual plot :

• Separatipn Efficiency = (54/91)*100 = 59.34%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

case 5 : 5 m/s.

Residual plot :

• Separatipn Efficiency = (52/91)*100 = 57.14%.

Anthracite particle time :

Pressure fluent PT for anthracite :

Velocity fluent PT for anthracite :

Conclusion :

• The highest separation efficiency is achieved for 1m/s case.
• The mesh needs to be more finer so that we can have more particles.
• The pressure drop increases with increase in velocity.

*______________________________________*_________________________________________*________________________________________________*_________________________________*