Script for computing the temporal evolution of the Lacey mixing index for a binary mixture of particles in EDEM
Altair EmployeeOverview
This python script computes the temporal evolution of the Lacey mixing index [1] for a binary mixture of particles in EDEM from the mixture component mass fractions in a dense grid of rectilinear sampling bins as described below.
An example result from the script is shown in Figure 1. A value of zero corresponds to the perfectly segregated state and a value of one to the randomly mixed state. In addition to the graphs shown in Figure 1 the script outputs the results in a comma delimited text file.
Figure 1 a) Dense grid of rectilinear bins for the calculation of the Lacey mixing index and b) the temopral evolution of the Lacey mixing index.
[1] Gu, Z. and Chen, J. J. J. (2015) ‘A probabilistic analysis of some selected mixing indices’, Chemical Engineering Research and Design. Institution of Chemical Engineers, 93(April), pp. 293–303. doi: 10.1016/j.cherd.2014.04.014.
Pre-Requisite
EDEM 2022.3 and EDEMpy 1.2.0
Usage/Installation Instructions
The settings of the script are contained in the Lacey_settings.txt file and defined in Figure 2.
Figure 2 Inputs for the Lacey mixing index script
The python file can be executed directly from the EDEM Analyst using the Run EDEMpy script feature (File->Run EDEMpy script). The EDEM deck files, Lacey.py file and Lacey_settings.txt file should all be in the same folder. The script will find and post-process all EDEM decks in the folder that it is in and output the results in the same location.
Comments
-
Good afternoon.
Can you please extend the above code (which calculates the lacey mixing index) for n species and n timesteps.
I tried but couldn't do it,Thanking you in anticipation.
Thanks & Regard
Akash Pal
0 -
May I ask why several mixing index values are negative?
0 -
Alice Sun said:
May I ask why several mixing index values are negative?
Hi Yurui,
The Lacey index is between 0 and 1 only if the mass or number fractions of the two materials are equal (50:50). Otherwise it is normal to get negative values but not values above 1.
Best regards,
Stefan
0 -
AKASH_22404 said:
Good afternoon.
Can you please extend the above code (which calculates the lacey mixing index) for n species and n timesteps.
I tried but couldn't do it,Thanking you in anticipation.
Thanks & Regard
Akash Pal
Hi Pal,
The Lacey mixing index cannot be generalised to higher order mixtures. Other indices have been proposed for this. For example:
Cho, M., Dutta, P. and Shim, J. (2017) ‘A non-sampling mixing index for multicomponent mixtures’, Powder Technology. Elsevier B.V., 319, pp. 434–444. doi: 10.1016/j.powtec.2017.07.011.
Best regards,
Stefan
0 -
Stefan Pantaleev_21979 said:
Hi Yurui,
The Lacey index is between 0 and 1 only if the mass or number fractions of the two materials are equal (50:50). Otherwise it is normal to get negative values but not values above 1.
Best regards,
Stefan
Thanks very much for your reply.
I use the EDEM 2023 and run this case. Why is the result different from the attached picture? There are several negative values.
Also, Can I run this script using EDEM 2023?
0 -
Alice Sun said:
Thanks very much for your reply.
I use the EDEM 2023 and run this case. Why is the result different from the attached picture? There are several negative values.
Also, Can I run this script using EDEM 2023?
Hi Yurui,
It is also possible to get some negative values or values above 1 if the number if bins or number of particles are low. We updated the script since the images were generated so you may get a different result to the image with the default values. The script should run ok in EDEM 2023.0 but let us know if you get any errors and we can update it.
Best regards,
Stefan
0 -
Stefan Pantaleev_21979 said:
Hi Yurui,
It is also possible to get some negative values or values above 1 if the number if bins or number of particles are low. We updated the script since the images were generated so you may get a different result to the image with the default values. The script should run ok in EDEM 2023.0 but let us know if you get any errors and we can update it.
Best regards,
Stefan
Thank you for your kindness.
Can I use this script directly for other binary mixture case (different mixer)? And where should I revise?
0 -
Alice Sun said:
Thank you for your kindness.
Can I use this script directly for other binary mixture case (different mixer)? And where should I revise?
You can use it for any geometry and any binary mixture. You only need to change the grid dimensions and might need to change the number of bins.
Best regards,
Stefan
0 -
Stefan Pantaleev_21979 said:
You can use it for any geometry and any binary mixture. You only need to change the grid dimensions and might need to change the number of bins.
Best regards,
Stefan
Great! Thanks again
0 -
Hi, Stefan
I am also confused as to what N refers to in this script?
I saw the picture you uploaded about the principle of calculating the Lacey mixing index. I also checked some literature, N should refer to the average number of particles in each bin.
Maybe you can explain why take N is the average number of particles of type 1 in each bin？
0 -
Community Guest said:
Hi, Stefan
I am also confused as to what N refers to in this script?
I saw the picture you uploaded about the principle of calculating the Lacey mixing index. I also checked some literature, N should refer to the average number of particles in each bin.
Maybe you can explain why take N is the average number of particles of type 1 in each bin？
Hi .,
Thank you for bringing this up. There is a typo in the above description. Is should say N is the average number of particles per bin. We will correct this. The calculation in the script is using the correct definition of N.
Best regards,
Stefan
0