Cyclone Seprator

1. Create Case

After opening SimFlow, we will now create a new case cyclone_separator

  1. Go to New panel

  2. Provide name cyclone_separator

  3. Click Create Case

cs 1 create case

2. Import Geometry

After creating case Download Geometry: Cyclone

  1. Click Import Geometry

  2. Select geometry file cyclone.stl

  3. Click Open

cs 2 import geometry

3. Display Geometry

After importing geometry, it will appear in the 3D window

  1. Click Fit View to zoom out the geometry

cs 3 geometry view 1

4. Scale Geometry (I)

Geometry was created in millimeters, we will scale it to meters

  1. Click Options next to the cyclone geometry

  2. Select Scale from the options list

cs 4 scale geometry 1

5. Scale Geometry (II)

Scale Geometry window will appear on the left side

  1. Check Uniform scaling option

  2. Set Scaling Factor to 1e-03

  3. Click Scale

cs 5 scale geometry 2

6. Fit View Geometry

After rescaling geometry, we should zoom in to see the geometry properly

  1. Click Fit View to zoom in the geometry

cs 3 geometry view 1

7. Meshing Parameters - Cyclone

Now we will set meshing parameters for cyclone geometry

  1. Go to Hex Meshing panel

  2. Enable meshing on the cyclone geometry

cs 7 mesh geometry

8. Base Mesh

We will define the base mesh now

  1. Go to Base tab

  2. Define initial mesh extends
    Min \(\sf{[m]}\)-0.1-0.1-0.4499
    Max \(\sf{[m]}\)0.14990.1010.5499
    (these dimensions have been set, so that inlet, bottom and top faces of the geometry fall slightly outside the box)

  3. Define the number of divisions
    Division2520100

cs 8 base mesh

9. Base Mesh Boundaries

To be able to define different conditions on each boundary of the domain we need to assign an individual name to each side of the base mesh

  1. Change the following boundary names accordingly
    X+ inlet
    Z- bottom
    Z+ top

  2. Define boundary types accordingly
    X- wall
    X+ patch
    Y- wall
    Y+ wall
    Z- wall
    Z+ patch

cs 9 base mesh boundaries

10. Material Point

Make sure the material point is located inside the geometry

  1. Go to Point tab

  2. Make sure that the material point is set accordingly
    Material Point000

cs 10 material point

11. Start Meshing

Now everything is set up, we can begin the meshing process

  1. Go to Mesh tab

  2. Click Mesh to start the meshing process

cs 11 meshing

12. Mesh

After meshing process is finished the mesh should appear in the 3D graphics window

  1. Click Fit View

cs 12 mesh

13. Setup Solver - MPPIC

We will use MMPICFoam solver. This is a Lagrangian solver based on the Discrete Particle Modelling(DPM) method with Multiphase Particle-in-Cell (MPPIC) collision handling.

  1. Go to Setup panel

  2. Select Transient filter

  3. Select Lagrangian model

  4. Pick MPPIC (MMPICFoam) solver

  5. Select solver

cs 13 setup solver

14. Discrete Phase - Properties

We will now define properties of the discrete phase

  1. Go to Discrete Phase panel

  2. Set density of the discrete phase
    \(\rho_0\) \(\sf{[\frac{kg}{m^2}]}\)2500

  3. Set packing factor
    \(\alpha_{packed}\) \(\sf{[-]}\)0.6

cs 14 discrete phase properties

15. Discrete Phase - Injection

We will now define the injection of the discrete phase through the inlet boundary. Mass flow rate of the particles is equal to the mass flow rate of the air that will be defined later in the Boundary Conditions panel.

  1. Go to Injection tab

  2. Create Boundary Injector

  3. 4, 5. Set following parameters accordingly
    Total Mass \(\sf{[kg]}\)1.225
    SOI \(\sf{[s]}\)1
    Duration \(\sf{[s]}\)10
    Parcels Per Second10000
    Boundaryinlet
    \(U_0\) \(\sf{[m/s]}\)-300

cs 15 discrete phase injection 1

16. Discrete Phase - Distribution

We will now define a normal distribution of the discrete phase size

  1. Go to Distribution tab

  2. Set Distribution to Normal

  3. Set distribution parameters accordingly
    Min \(\sf{[m]}\)4e-05
    Max \(\sf{[m]}\)3.6e-04
    \(\mu\) \(\sf{[m]}\)5e-05
    \(\sigma\) \(\sf{[m]}\)2e-04

cs 16 discrete phase distribution

17. Discrete Phase - Models

We will now define the particle drag model

  1. Go to Models tab

  2. Select Ergun-Wen-Yu model

cs 17 discrete phase models

18. Discrete Phase - Solution

We will now define methods of computing and interpolating an average of the Lagrangian phase

  1. Go to Solution tab

  2. Set Averaging method to Dual

  3. Expand Source Terms options

  4. Set Semi-Implicit discretization of the U term

cs 18 discrete phase solution

19. Turbulence

For turbulence modeling, we will use the LES model

  1. Go to Turbulence panel

  2. Select LES turbulence modeling

  3. Select \(k \; Equation\) model

cs 19 turbulence

20. Transport Properties - Fluid

Now we will define the transport properties of fluid material

  1. Go to Transport Properties panel

  2. Click Material Database

  3. Select air material

  4. Click Apply

cs 20 transport properties fluid

21. Solution - PIMPLE

To increase the stability of the simulation we will increase the number of pressure corrector iterations

  1. Go to Solution panel

  2. Select the PIMPLE tab

  3. Increase the number of Correctors to 2

cs 21 solution pimple

22. Boundary Conditions - Bottom (Particles)

Now we will set a bottom boundary to be transmissive for particles

  1. Go to Boundary Conditions panel

  2. Select bottom boundary

  3. Select Particles tab

  4. Change particle interaction to Escape

cs 22 boundary conditions bottom particles

23. Boundary Conditions - Cyclone (Particles)

  1. Select cyclone boundary

  2. Set following values accordingly
    e \(\sf{[-]}\)0.97
    \(\mu\) \(\sf{[-]}\)0.09

cs 23 boundary conditions cyclone particles

24. Boundary Conditions - Inlet (Particles)

We will now set boundary conditions on inlet boundary

  1. Select inlet boundary

  2. Change boundary condition for inlet to Velocity Inlet

  3. Change particle interaction to Rebound

cs 24 boundary conditions inlet particles

25. Boundary Conditions - Inlet (Flow)

  1. Switch to Flow tab

  2. Set velocity at inlet
    Reference Value \(\sf{[m/s]}\)20

cs 25 boundary conditions inlet flow

26. Run - Time Controls

  1. Go to Run panel

  2. Set Simulation Time [s] to 5

  3. Set time step \(\Delta t [s]\) to 1e-04

cs 26 run time controls

27. Run - Output

  1. Go to Output tab

  2. Set Write Control Interval [s] to 0.1 seconds
    (it will force solver to write results on the hard drive every 0.1 seconds of the simulation)

cs 27 run output

28. Run - CPU

This simulation will require high CPU usage and might take several hours to finish the calculation depending on the machine used. To speed up the calculation process increase the number of CPUs basing on your PC capability. We recommend using at least 4 cores for this tutorial. The free version allows you to use only 2 processors in parallel mode. To get the full version, you can use the contact form to Request 30-day Trial

Estimated computation time for 2 processors: 4 hours

  1. Switch to CPU tab

  2. Use parallel mode

  3. Increase the Number of processors

  4. Click Run Simulation button

cs 27 run cpu

29. Postprocessing - Open ParaView

Open ParaView software to display results

  1. Go to Postprocessing panel

  2. Start ParaView

cs 27 postprocessing open paraview

30. ParaView - Import Results

After opening the ParaView we will import result from the simulation. We will import fluid region and parcels cloud separately to be able to assign them different display properties.

  1. Select cyclone_separator.foam

  2. Click Apply to import results

  3. Click Last Frame to select the latest result set

  4. Refresh the results

  5. Uncheck lagrangian/kinematicCloud region

  6. Click Apply

cs 30 paraview import results

31. ParaView - Streamlines

We will now plot streamlines for the fluid region and set the coloring variable to velocity

  1. Select cyclone_separator.foam

  2. Click Stream Tracer button to add streamlines

  3. Click Apply

  4. Change the coloring variable to U.bulk

cs 31 paraview streamlines

32. ParaView - Display Geometry

To show the transparent geometry together with the streamlines we will follow below steps:

  1. Select cyclone_separator_foam

  2. Click on the eye next to cyclone_separator_foam

  3. Select Solid Color from the list

  4. Change the opacity in a properties tab to 0.3

cs 32 paraview geometry opacity

33. ParaView - Display Particles (I)

In order to display particles you can import the same case into ParaView once again

  1. Select Open from the top menu

  2. Make sure you are in folder with cyclone_separator_foam case

  3. Select cyclone_separator.foam from file

  4. Click OK

cs 33 paraview import geometry 1

34. ParaView - Display Particles (II)

Now we will select to display particles only

  1. Click Apply

  2. Go to Mesh regions, uncheck all regions and check only lagrangian/kinematicCloud

  3. Click once again Apply to confirm

  4. Change the coloring variable to the age

Note that after first applying, lagrangian/kinematicCloud will appear on the list of mesh regions.
cs 34 paraview display particles

35. ParaView - Results

The results are displayed in the graphics window.

  1. Play with an animation buttons to track the results of analysis

Note that this tutorial is meant only to demonstrate the capabilities of the software and not to solve the problem in the best possible way. Therefore, some assumptions are taken to keep case setup time and computational time low. In particular, to refine the model, one could in the first place consider refining the grid and choosing a more suitable drag model for the particles (e.g. Ergun-Wen-Yu model). Subsequently, it is worth considering enabling isotropic particle packing, isotropic particle timescale, and stochastic isotropy model.
cs 35 paraview results