ACU-T: 6100 Particle Separation in a Windshifter using Altair EDEM

This tutorial introduces you to the workflow for setting up and running a basic sequential coupling (one-way steady) simulation using AcuSolve and EDEM. Prior to starting this tutorial, you should have already run through the introductory tutorial, ACU-T: 1000 Basic Flow Set Up, and have a basic understanding of HyperMesh CFD, AcuSolve, and EDEM. To run this simulation, you will need access to a licensed version of HyperMesh CFD, AcuSolve, and EDEM.

Before you begin, copy the file(s) used in this tutorial to your working directory.
Note: This tutorial does not cover the steps related to geometry cleanup and mesh settings.

Problem Description

The problem to be solved is shown schematically in Figure 1. It is a windshifter model in which air enters the domain from the inlet at the bottom at a very high velocity (20 m/sec) and exits through the outlet. When particles are introduced into the domain through the particle inlet, the lighter particles get carried away with the air and the heavier particles exit through the opening at the bottom. Since the concentration of the particles is sparse, the effect of the particles on the fluid field is not considered for this simulation. Hence, one-way coupling between AcuSolve and EDEM is used to simulate the separation of particles; only the effect of the fluid forces on the particles is considered.


Figure 1.
The model consists of a cylindrical pipe with a 45-degree bend. The radius of the pipe is 0.25 m and the particle inlet is located midway through the length of the pipe. To simplify things, the particle inlet is considered as a wall for the CFD simulation. The workflow for an AcuSolve-EDEM sequential coupling simulation is shown below.


Figure 2.

The tutorial consists of two parts:

  1. AcuSolve simulation
  2. EDEM simulation

The setup for the AcuSolve simulation is done using HyperMesh CFD. Once the AcuSolve simulation is complete, the velocity nodal data is exported into CGNS format using the AcuTrans utility available in AcuSolve. This field data will be imported into EDEM and will be used to calculate the drag forces on the particles due to the fluid.

Two different bulk materials used in the EDEM simulation and their properties are listed below:

Name Density (kg/m3) Particle radius (m) Average weight of individual particle (kg) Rate of generation (particles per sec)
Heavy particle 900 0.03 0.004 100
Light particle 100 0.03 0.0004 100

The fluid drag forces on the particles are calculated using the Schiller-Naumann drag model. The drag force on the particle is calculated using the relative velocity of the particle with respect to the local fluid velocity. This particle body force will used to update the location of the particles and then the new drag forces are calculated. This loop is repeated until the end of the simulation.

Part 1 - AcuSolve Simulation

Start HyperMesh CFD and Open the HyperMesh Database

  1. Start HyperMesh CFD from the Windows Start menu by clicking Start > Altair <version> > HyperMesh CFD.
  2. From the Home tools, Files tool group, click the Open Model tool.


    Figure 3.
    The Open File dialog opens.
  3. Browse to the directory where you saved the model file. Select the HyperMesh file ACU-T6100_windshifter.hm and click Open.
  4. Click File > Save As.
  5. Save the database as windshifter_sequential in the same directory as the other input files.
    This will be the working directory and all the files related to the simulation will be stored in this location.

Validate the Geometry

The Validate tool scans through the entire model, performs checks on the surfaces and solids, and flags any defects in the geometry, such as free edges, closed shells, intersections, duplicates, and slivers.

To focus on the physics part of the simulation, this tutorial input file contains geometry which has already been validated. Observe that a blue check mark appears on the top-left corner of the Validate icon on the Geometry ribbon. This indicates that the geometry is valid, and you can go to the flow set up.


Figure 4.

Set Up Flow

Set the General Simulation Parameters

  1. From the Flow ribbon, click the Physics tool.


    Figure 5.
    The Setup dialog opens.
  2. Under the Physics models setting:
    1. Verify that the Incompressible option is selected under Single phase flow.
    2. Activate the Export geometry to EDEM check-box.
    3. Set Time marching to Steady.
    4. Select Spalart-Allmaras as the Turbulence model.


    Figure 6.
  3. Close the dialog and save the model.

Assign Material Properties

  1. From the Flow ribbon, click the Material tool.


    Figure 7.
  2. Verify that Air has been assigned as the material.
  3. On the guide bar, click to exit the tool.

Define Flow Boundary Conditions

  1. From the Flow ribbon, Profiled tool group, click the Profiled Inlet tool.


    Figure 8.
  2. Click on the inlet face highlighted in the figure below and enter a value of 20 m/s for the Average velocity in the microdialog.


    Figure 9.
  3. On the guide bar, click to execute the command and exit the tool.
  4. Click the Outlet tool.


    Figure 10.
  5. Select the face highlighted in the figure below and then click on the guide bar.


    Figure 11.
  6. Click the No Slip tool.


    Figure 12.
  7. Select the face highlighted in the figure below.


    Figure 13.
  8. In the Boundaries legend, double-click Wall and rename it to Particle_inlet.
    You will use this surface as a reference geometry to create a particle factory in EDEM. Hence, you are placing this surface in a different geometry group.
  9. Click on the guide bar.
  10. Save the model.

Generate the Mesh

To focus on the solver setup, the mesh settings are predefined in the input file given to you.
  1. From the Mesh ribbon, click the Volume tool.


    Figure 14.
    The Meshing Operations dialog opens.


    Figure 15.
  2. Click Mesh.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.
    Tip: Right-click on the mesh job and select View log file to view a summary of the meshing process.
  3. Save the model.

Run AcuSolve

  1. From the Solution ribbon, click the Run tool.


    Figure 16.
    The Launch AcuSolve dialog opens.
  2. Set the Parallel processing option to Intel MPI.
  3. Optional: Set the number of processors to 4 or 8 based on availability.
  4. Click Run to launch AcuSolve.


    Figure 17.

    Once the AcuSolve simulation is launched, HyperMesh CFD also exports the EDEM files which contain the geometry. The surfaces are organized based on the CFD setup.

Export Velocity Field Data

  1. From the Solution ribbon, click the Convert tool.


    Figure 18.
    The AcuSolve Data Translator dialog opens.
  2. In the dialog, click Browse, navigate to the AcuSolve problem directory, and open the windshifter_sequential.1.Log file.
  3. Set the Output format to CGNS.
  4. Select velocity from the list of variables.
  5. Click Convert to export the velocity nodal data in CGNS format.


    Figure 19.

    Once the process is successful, the same information is printed in the Output field and a CGNS file is created in the AcuSolve problem directory.

Post-Process the Results with HM-CFD Post

  1. Once the solution is completed, navigate to the Post ribbon.
  2. From the menu bar, click File > Open > Results.
  3. Select the AcuSolve .log file in your problem directory to load the results for post-processing.
    The solid and all the surfaces are loaded in the Post Browser.
  4. Select the Left face on the view cube to align the model to the x-z plane.


    Figure 20.
  5. In the Post Browser, turn off the display of the boundary surfaces by clicking on the icon next to Flow Boundaries.


    Figure 21.
  6. Click the Slice Planes tool.


    Figure 22.
  7. In the modeling window, click on the plane that is parallel to the screen (the x-z plane).
  8. In the slice plane microdialog, click to create the slice plane.
  9. In the display properties microdialog, set the display to velocity and activate the Legend toggle.
  10. Click and set the Colormap name to Rainbow Uniform.


    Figure 23.
  11. Click on the guide bar.


    Figure 24.

Part 2 - EDEM Simulation

Start Altair EDEM from the Windows start menu by clicking Start > Altair 2023 > EDEM 2023 . The user interface in EDEM is divided into three tabs: Creator, Simulator and Analyst. The Creator is used to setup and initialize your model. It is where you import particles and geometries and define the other model parameters. The Simulator is where you configure and control the EDEM simulation engine, and where you can observe the progress of your simulation. The Analyst is the post-processor used to analyze and visualize the results of your simulation.

The steps required to setup and run a basic CFD-DEM one-way coupling are explained below. Refer to the EDEM help documentation for more details.

Open the EDEM Input Deck

As mentioned earlier, when the AcuSolve simulation was launched, HyperMesh CFD created a set of EDEM files in the problem directory. You will open that EDEM input deck and setup the DEM simulation.

  1. In the Creator tab in EDEM, go to File > Open.
  2. In the dialog, browse to the AcuSolve problem directory and open the windshifter_sequential.dem file located in the EDEM folder.
    The geometry is loaded.
  3. Click the Environment tab under the Creator Tree, uncheck Auto Update from Geometry, and then check the box again to fit the geometry within the boundary.


    Figure 25.


    Figure 26.

Define the Bulk Materials and Equipment Material

In this step, you will define the material models for the heavy and light bulk material and the equipment material.

  1. In the Creator Tree, right-click Bulk Material and select Add Bulk Material.
  2. Rename the material to Heavy.
  3. In the Creator Tree, set the Solids Density property to 900 kg/m3.
    You will use the default values for other properties for this tutorial.


    Figure 27.
  4. Click below Interaction to define the interaction properties for collisions among the heavy particles. In the dialog, click OK.
  5. In the Creator Tree, right click Heavy and select Add Shape from Library > Dual Sphere Shape.
  6. Rename the particle to Heavy particle.
  7. Under Heavy particle, click Properties.
  8. In the Heavy particle Spheres panel, set the Physical Radius of both the spheres to 0.03 m and press Enter.


    Figure 28.
  9. In the Creator Tree, click Calculate Properties.


    Figure 29.
  10. In the Creator Tree, right-click Bulk Material and select Add Bulk Material.
  11. Rename the material to Light.
  12. In the Creator Tree, set the Solids Density property to 100 kg/m3.
    You will use the default values for other properties for this tutorial.


    Figure 30.
  13. Click below Interaction to define the interaction properties for collisions among the heavy particles. In the dialog, select Heavy and then click OK.
  14. Click again to define the interaction properties for collisions among the light particles. In the dialog, select Light and then click OK.
  15. In the Creator Tree, right click on Light and select Add Shape from Library > Dual Sphere Shape.
  16. Rename the particle to Light particle.
  17. Under Light particle, click Properties.
  18. In the Light particle Spheres panel, set the Physical Radius of both the spheres to 0.03 m and press Enter.


    Figure 31.
  19. In the Creator Tree, click Calculate Properties.


    Figure 32.
  20. In the Creator tree, right-click Equipment Material and select Add Equipment Material. Rename it to Steel.
  21. Set the Density to 7800 kg/m3.
  22. Click below Interaction to define the interaction properties for collisions among the heavy particles. In the dialog, select Heavy and then click OK.
  23. Click again to define the interaction properties for collisions among the light particles. In the dialog, select Light and then click OK.
  24. Save the model.

Create the Particle Factory

  1. Expand Geometries under the Creator Tree tab. Next, right-click the Particle inlet surface group and select Copy Geometry > Single Copy.
  2. Rename the new geometry section to Particle_factory.
  3. Right-click the Default Wall geometry section in the Creator Tree and select Merge Geometry(s).
  4. In the Merge Geometry dialog, select the Particle_inlet and then click OK.
  5. Click the Default Wall geometry section, set the Type to Physical, and the Material to Steel (if not set already).


    Figure 33.
  6. In the Creator Tree, click the Inlet section and change the Type to Virtual.
  7. Similarly, change the Type to Virtual for the Outlet and Particle_factory sections.
  8. Under Particle_factory, click Transform. Set the X-Position to 0.035 m


    Figure 34.
    Note: Set the Opacity value to 0.2 to see the transformed surface location inside the pipe geometry.

    This is done to make sure that the particles are generated inside the fluid domain.

Define the Particle Factory

Now that the bulk material, geometry sections, and equipment materials are defined, you need to create a particle factory to generate the particles. You will create one factory for each bulk material.

  1. In the Creator Tree, right-click Particle_factory and select Add Factory.
  2. Rename the new factory to Heavy factory.
  3. Set the particle generation parameters as shown in the figure below.


    Figure 35.
  4. Click besides Velocity, set the X-velocity to 1 m/s, and then click OK.
  5. Repeat steps 1-4 to create another factory named Light factory using the same parameters but with Light as the Material.

Define the Physics and Import CFD Field Data

In this step, you will define the physics models for particle collisions and the particle body force.

  1. In the Creator Tree, click Physics.
  2. Click the Interaction drop-down menu and select Particle Body Force.


    Figure 36.
  3. Click Edit Contact Chain.
  4. In the dialog, select the Schiller and Naumann Drag model and then click OK.
  5. In the Physics tab, click Schiller and Naumann Drag and then click .
  6. In the Schiller and Naumann Drag Parameter Values dialog, enter the values as shown in the figure below and then click OK.


    Figure 37.
  7. In the Physics tab, click the Field Data Manager.
  8. In the Field Data Manager dialog, click then browse to the AcuSolve problem directory where you saved the CGNS file exported using the CFD field data. Open windshifter_sequential_run1.cgns.
  9. Once the data import is complete, double-click Velocity and rename it to velocity. (The letter ‘v’ should be in lower case). Then, close the dialog.


    Figure 38.
  10. Save the EDEM deck.

Define the Environment

In this step, you will define the extents of the domain for the EDEM simulation and the direction of gravitational acceleration.

  1. In the Creator Tree, click Environment.
  2. Activate the checkbox for Auto Update from Geometry (if not already selected).
    When a moving particle touches the bounding faces of the domain (environment), it will be removed from the simulation.
  3. Activate Gravity and set the z-value to -9.81 m/s2.
  4. Save the EDEM deck.

Define the Simulation Settings and Run the Simulation

  1. Click in the top-left corner to go to the EDEM Simulator tab.
  2. In the Simulator Settings tab, set the Time Integration scheme to Euler and activate the Auto Time Step checkbox (if not set already).
  3. Set the Total Time to 1 s and the Target Save Interval to 0.01 s.
  4. Set the Cell Size to 4 R min.
    Generally, a value in the range of 3-6 Rmin is recommended as the optimum cell size. The cell size in EDEM doesn’t have any impact on the accuracy of the simulation and affects only the run time.
  5. Set the Selected Engine to CPU Solver and set the Number of CPU Cores based on availability.


    Figure 39.
  6. Once the simulation settings have been defined, click to start the EDEM simulation.

Analyze the Results

  1. Once the EDEM simulation is complete, click in the top-left corner to go to the EDEM Analyst tab.
  2. In the Analyst Tree, expand Display > Geometries and then click Default Wall.
  3. Verify that the Display Mode is set to Filled and set the Opacity to 0.2.


    Figure 40.
  4. In the Analyst Tree, expand Particles and click Heavy particle.
  5. Change the display color to Magenta.


    Figure 41.
  6. Click Light particle and set the display color to Green.


    Figure 42.
  7. On the menu bar, set the time to 0 by clicking:


    Figure 43.
  8. Set the View plane to + Y.


    Figure 44.
  9. In the Viewer window, set the Playback Speed to 0.1x then click on the play icon to play the particle flow animation.


    Figure 45.


    Figure 46.

    Observe that the lighter particles (green) get carried by the fluid and escape the domain through the outlet at the top and the heavier particles (magenta) stay inside the domain for a longer time while some of them fall through the bottom of the pipe.

Summary

In tutorial you learned how to setup a basic AcuSolve-EDEM sequential (one-way steady) coupling problem. In the first part, you set up and solved a steady state flow simulation using AcuSolve and then exported the CFD field data using AcuTrans. Next, you set up the EDEM model and imported the CFD field data using the Field Data Manager. Once the EDEM simulation was completed, you learned how to view results and create animations in EDEM Analyst.