ACU-T: 6106 AcuSolve - EDEM Bidirectional Coupling with Mass Transfer

Prerequisites

This tutorial introduces you to the workflow for setting up and running a basic AcuSolve-EDEM bidirectional coupling simulation with mass transfer. Prior to starting this tutorial, you should have already run through the introductory tutorial, ACU-T: 1000 Basic Flow Set Up and ACU-T: 6100 Particle Separation in a Windshifter using Altair EDEM, 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.

Problem Description

The spray coating update physics models allows you to simulate the coating process in which the mass and volume of a ‘spray’ particle is transferred to a solid particle when they come into contact. In this tutorial, a tablet coating process is simulated, in which spray particles are sprayed onto a row of tablets. The spray particles are modeled as solid particles in EDEM, and their properties are set to that of water. When a spray particle touches the tablet, the volume of spray particle is transferred to the tablet and is registered a custom property called ‘Volume Added’. This custom property can be used as a metric for quantifying the coating received by the tablet.


Figure 1.
In a typical tablet coater, hot air is injected during the coating process. This allows for the solvent to be evaporated and the solid content of the coating gets deposited on the surface of the tablet. By coupling the spray coating model with the heat and mass transfer models in AcuSolve, you can simulate the drying process as well. In this tutorial, the initial temperature of the air is set to 363 K and the walls are considered as adiabatic. The tablet particles are created with an initial temperature of 350 K. The spray particles are injected for the first 0.25 s of the simulation and then the coating is left to evaporate. In addition to the Volume Added (coating volume), the evaporation rate on each particle can also be monitored.


Figure 2.

Part 1 - EDEM Simulation

Start Altair EDEM from the Windows start menu by clicking Start > Altair 2023 > EDEM 2023 .

Open the EDEM Input Deck

  1. From the menu bar, click File > Open.
  2. In the dialog, browse to your problem directory and open the tablet.dem file.
    The geometry and the materials are loaded.


    Figure 3.

Review the Bulk Material and Particle Shape

  1. Right-click in the Creator Tree and select Expand All.
  2. Under Bulk Material, click spray and verify that the properties are set as shown below.


    Figure 4.
  3. Click Properties under spray. Verify that the particle radius is set to 0.0001 m.
  4. Click the tablet bulk material and verify that the properties are set as shown below.


    Figure 5.
  5. Click Properties under tablet and review the tablet shape.


    Figure 6.

    The tablet is made up of four different shapes and the overall dimension of the tablet is 0.004 m x 0.002 m x 0.002 m. The coordinate digit in the figure is slightly different from those in the database.

Set the Physics Models

  1. In the Creator Tree, click Physics.
  2. Set the Interaction to Particle to Particle then click Edit Contact Chain.


    Figure 7.
  3. In the dialog, activate the checkbox for the Spray Coating model. Set the Base and Friction models as shown in the figure below and click OK.


    Figure 8.
  4. In the Creator Tree, select Spray Coating and then click .


    Figure 9.
  5. In the Spray Coating Model Configuration dialog, click , select spray from the list of materials, and then click OK.


    Figure 10.
  6. In the Creator Tree, change the interaction to Particle to Geometry and then click Edit Contact Chain.
  7. In the dialog, set the same parameters as above and then click OK.


    Figure 11.
  8. In the Creator Tree, select Spray Coating and then click .
  9. In the Spray Coating Model Configuration dialog, click , select spray from the list of materials, and then click OK.
  10. In the Creator Tree, change the interaction to Particle Body Force and then click Edit Contact Chain.
    Note: In this simple example, the tablet particles are created such that there are no tablet-tablet and tablet-geometry contacts. Hence the contact models are turned off to save computational time. The spray-tablet contacts are still active and are handled by the spray coating physics model.
  11. In the dialog, activate the checkboxes for Temperature Update and Spray Coating Update and then click OK.


    Figure 12.
  12. In the Creator Tree, select Temperature Update and then click .


    Figure 13.
  13. Set the specific heat capacity of tablet and spray to 1200 and 4000 J/kg-K, respectively, and then click OK.
  14. Save the model.

Define Geometries and Factories

  1. Under Geometries, click box and verify that the type is set to Virtual.
  2. Right-click box and select Add Factory.
  3. Right click New Factory 1 and select Change Factory Type.
  4. Verify that the Factory Type is set to static and set the particle generation parameters as shown in the figure below. Make sure that the Material/Meta-Particle is set to tablet.


    Figure 14.
  5. Under Parameters, set the Position option to cubic and then click .


    Figure 15.
  6. In the dialog, enter the position parameters as shown in the figure below and then click OK.
    This creates a row of tablet particles spaced evenly along the y-axis and located at the center of the box.


    Figure 16.
  7. Similarly, click next to the Temperature setting, set the particle temperature to 350 K, and then click OK.
  8. Under Geometries, click spray and then change the type to Virtual.
  9. Right-click spray and select Add Factory.
  10. Set the particle generation parameters as shown in the figure below. Set the Material to spray.


    Figure 17.
  11. Under Parameters, set the Velocity option to spray and then click .
  12. In the dialog, set the spray velocity parameters as shown in the figure below and then click OK.


    Figure 18.
  13. Similarly, click next to the Temperature setting, set the particle temperature to 350 K, and then click OK.
  14. Save the model.

Define the Simulation Settings

  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 de-activate the Auto Time Step checkbox.
  3. Set the Fixed Time Step to 1e-5 s.
    Note: Generally, a value of 20-40% of the Rayleigh Time Step is recommended as the time step size to ensure stability of the DEM simulation.
  4. Set the Total Time to 0.5 s and the Target Save Interval to 0.01 s.
  5. Set the Cell Size to 50 R min.
  6. Set the Selected Engine to CPU Solver and set the Number of CPU Cores based on availability.


    Figure 19.
  7. Once the simulation settings have been defined, save the model.

Part 2 - 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 20.
    The Open File dialog opens.
  3. Browse to the directory where you saved the model file. Select the HyperMesh file ACU-T6106_tablet.hm and click Open.

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

Set Up Flow

Set the General Simulation Parameters

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


    Figure 22.
    The Setup dialog opens.
  2. Under the Physics models setting, select the Multiphase flow radio button.
  3. Change the Multifluid type to Bidirectional EDEM Coupling.
  4. Set Time step size and Final time to 0.001 and 0.5, respectively. Select Spalart-Allmaras for the Turbulence model.
  5. Set the Pressure scale to Absolute.
  6. Activate Heat transfer and set the Temperature scale to Absolute.
  7. Set the Number of passive species to 1.

    This setting is necessary for the mass transfer feature to be active. The evaporated moisture is tracked using the species transport equation.



    Figure 23.
  8. Click the Solver controls setting. Set the Minimum and Maximum stagger iterations to 2 and 2, respectively.


    Figure 24.
  9. Close the dialog and save the model.

Assign Material Properties

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


    Figure 25.
  2. Verify that Air-EDEM Particle – 2way has been assigned as the material.
    If not assigned, click the box geometry and select Air-EDEM Particle – 2way from the microdialog.
  3. On the guide bar, click to exit the tool.

    Since you are applying an adiabatic condition on the box walls, you don’t need to set any thermal boundary condition explicitly.

Generate the Mesh

  1. From the Mesh ribbon, click the Volume tool.


    Figure 26.
    The Meshing Operations dialog opens.
  2. Set the Mesh size to Average size and change the Maximum element size to 0.005.
  3. Deactivate Curvature-based surface refinement and then click Mesh.


    Figure 27.
  4. Once the meshing process is complete, save the model.

Define Nodal Outputs

Once the meshing is complete, you are automatically taken to the Solution ribbon.

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


    Figure 28.
    The Field Output dialog opens.
  2. Check the box for Write initial conditions.
  3. Uncheck the box for Write results at time step interval
  4. Check the box for Write results at time interval.
  5. Set the Time step interval to 0.01.


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

Submit the Coupled Simulation

  1. Start the coupling server by clicking Coupling Server in EDEM.


    Figure 30.
    Once the Coupling server is activated, the icon changes.


    Figure 31.
  2. Return to HyperMesh CFD.
  3. From the Solution ribbon, click the Run tool.


    Figure 32.
    The Launch AcuSolve dialog opens.
  4. Set the Parallel processing option to Intel MPI.
  5. Set the Number of processors to 6.
  6. Expand Default initial conditions, uncheck Pre-compute flow and set the velocity values to 0. Uncheck Pre-compute Turbulence.
  7. Set the Temperature to 363 K.
  8. Click Run to launch AcuSolve.


    Figure 33.
    Once the AcuSolve simulation is launched, it connects with the EDEM coupling server. A warning message is displayed indicating that the base contact models have not been set. Click Ignore. Since the tablet particles are not in contact with each other or with the geometry, you don’t need any contact models. The spray-tablet contacts are handled by the spray coating model.
  9. In the dialog, right-click on the AcuSolve run and select View log file.
    If the coupling with EDEM is successful, that information is printed in the log file.


    Figure 34.
    Once the simulation is complete, the summary of the run time is printed at the end of the log file.


    Figure 35.

Post-Process the Results with EDEM

  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 select box.
  3. Verify that the Display Mode is set to Filled and set the Opacity to 0.2.


    Figure 36.
  4. In the Analyst Tree, click Particles then click tablet.
  5. Set the coloring to Volume Added.
  6. Activate the Auto Update checkbox for both Min and Max Value.
  7. Activate the Show Legend checkbox.
  8. Click Apply All.


    Figure 37.
  9. On the menu bar, set the time to 0 by clicking:


    Figure 38.
  10. Set the View plane to Default.


    Figure 39.
  11. In the Viewer window, set the Playback Speed to 0.1x and then click to play the particle flow animation.


    Figure 40.

    Observe that the added volume increases in the beginning as the tablets receive the spray coating. Once the spray injection is stopped, the coating starts to dry up and the added volume decreases gradually.

Summary

In this tutorial, you have learned how to set up and run a basic AcuSolve-EDEM bidirectional (two-way) coupling problem with mass transfer. You learned how to create particles with specific position and internal spacing and also learned how to define a spray injection. Once the simulation was completed, you processed the results to view the coating variation over time.