Please wait a minute...
Frontiers of Chemical Science and Engineering

ISSN 2095-0179

ISSN 2095-0187(Online)

CN 11-5981/TQ

Postal Subscription Code 80-969

2018 Impact Factor: 2.809

Front. Chem. Sci. Eng.    2022, Vol. 16 Issue (6) : 874-885    https://doi.org/10.1007/s11705-021-2131-1
RESEARCH ARTICLE
A “Sequential Design of Simulations” approach for exploiting and calibrating discrete element simulations of cohesive powders
Xizhong Chen1,2,3, Chunlei Pei1,4(), James A. Elliott1()
1. Department of Materials Science and Metallurgy, University of Cambridge, Cambridge CB3 0FS, UK
2. Process and Chemical Engineering, School of Engineering, University College Cork, Cork T12 K8AF, Ireland
3. Department of Chemical and Biological Engineering, The University of Sheffield, Sheffield S10 2TN, UK
4. Key Laboratory for Green Chemical Technology of Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
 Download: PDF(2266 KB)   HTML
 Export: BibTeX | EndNote | Reference Manager | ProCite | RefWorks
Abstract

The flow behaviours of cohesive particles in the ring shear test were simulated and examined using discrete element method guided by a design of experiments methodology. A full factorial design was used as a screening design to reveal the effects of material properties of partcles. An augmented design extending the screening design to a response surface design was constructed to establish the relations between macroscopic shear stresses and particle properties. It is found that the powder flow in the shear cell can be classified into four regimes. Shear stress is found to be sensitive to particle friction coefficient, surface energy and Young’s modulus. A considerable fluctuation of shear stress is observed in high friction and low cohesion regime. In high cohesion regime, Young’s modulus appears to have a more significant effect on the shear stress at the point of incipient flow than the shear stress during the pre-shear process. The predictions from response surface designs were validated and compared with shear stresses measured from the Schulze ring shear test. It is found that simulations and experiments showed excellent agreement under a variety of consolidation conditions, which verifies the advantages and feasibility of using the proposed “Sequential Design of Simulations” approach.

Keywords discrete element method      cohesive materials      parameter calibration      ring shear cell      design of experiments     
Corresponding Author(s): Chunlei Pei,James A. Elliott   
Online First Date: 01 March 2022    Issue Date: 28 June 2022
 Cite this article:   
Xizhong Chen,Chunlei Pei,James A. Elliott. A “Sequential Design of Simulations” approach for exploiting and calibrating discrete element simulations of cohesive powders[J]. Front. Chem. Sci. Eng., 2022, 16(6): 874-885.
 URL:  
https://academic.hep.com.cn/fcse/EN/10.1007/s11705-021-2131-1
https://academic.hep.com.cn/fcse/EN/Y2022/V16/I6/874
Fig.1  A SDoS approach for studying and calibrating DEM simulations.
Property Value
Particle density/(kg·m–3) 1600
Number of particles ~80000
Particle radius/μm 125 (30%), 75 (61%), 60 (9%)
Surface energy/(J·m–2) 0.0001–0.1
Young’s modulus/GPa 0.13–1.3
Poisson’s ratio 0.1–0.5
Particle-particle friction coefficient 0.1–0.8
Particle-particle restitution coefficient 0.3–0.9
Wall Young’s modulus/GPa 1.3
Wall Poisson’s ratio 0.3
Particle-wall friction coefficient 0.1
Particle-wall restitution coefficient 0.6
Tab.1  Simulation parameters used in ring shear test
Pattern Case ID Young’s modulus/GPa Poisson ratio Friction coefficient Restitution coefficient Surface energy/(J·m–2)
++++− 1 1.3 0.5 0.8 0.9 0.0001
+−−++ 2 1.3 0.1 0.1 0.9 0.1
++−++ 3 1.3 0.5 0.1 0.9 0.1
++−−− 4 1.3 0.5 0.1 0.3 0.0001
+−+−− 5 1.3 0.1 0.8 0.3 0.0001
−++−− 6 0.13 0.5 0.8 0.3 0.0001
−+++− 7 0.13 0.5 0.8 0.9 0.0001
−−+−+ 8 0.13 0.1 0.8 0.3 0.1
−−−+− 9 0.13 0.1 0.1 0.9 0.0001
−+−−− 10 0.13 0.5 0.1 0.3 0.0001
+−−−− 11 1.3 0.1 0.1 0.3 0.0001
+−−+− 12 1.3 0.1 0.1 0.9 0.0001
++−−+ 13 1.3 0.5 0.1 0.3 0.1
+−+−+ 14 1.3 0.1 0.8 0.3 0.1
−−−++ 15 0.13 0.1 0.1 0.9 0.1
−−−−+ 16 0.13 0.1 0.1 0.3 0.1
+−++− 17 1.3 0.1 0.8 0.9 0.0001
+−+++ 18 1.3 0.1 0.8 0.9 0.1
−++++ 19 0.13 0.5 0.8 0.9 0.1
−+−+− 20 0.13 0.5 0.1 0.9 0.0001
+−−−+ 21 1.3 0.1 0.1 0.3 0.1
+++++ 22 1.3 0.5 0.8 0.9 0.1
++−+− 23 1.3 0.5 0.1 0.9 0.0001
−+−−+ 24 0.13 0.5 0.1 0.3 0.1
+++−+ 25 1.3 0.5 0.8 0.3 0.1
−−−−− 26 0.13 0.1 0.1 0.3 0.0001
−−++− 27 0.13 0.1 0.8 0.9 0.0001
−−+−− 28 0.13 0.1 0.8 0.3 0.0001
−+−++ 29 0.13 0.5 0.1 0.9 0.1
+++−− 30 1.3 0.5 0.8 0.3 0.0001
−−+++ 31 0.13 0.1 0.8 0.9 0.1
−++−+ 32 0.13 0.5 0.8 0.3 0.1
Tab.2  Design table of the full factorial design used in this study
Fig.2  The simulated shear strain-stress results by the full factorial design.
Effect p-Value LogWorth
Friction coefficient 4.2e−26 25.37
Surface energy 2.23e−9 8.65
Young's modulus 7.13e−2 1.14
Restitution coefficient 4.04e−1 0.39
Poisson ratio 9.55e−1 0.02
Tab.3  The effect summary for pre-shear shear stress predicted by the full factorial design
Fig.3  The simulated pre-shear shear stresses by the full factorial design.
Effect p-Value LogWorth
Friction coefficient 9.0e−22 14.85
Surface energy 1.4e−15 21.05
Young’s modulus 8.07e−4 3.09
Restitution coefficient 7.42e−1 0.13
Poisson ratio 7.5e−1 0.12
Tab.4  The effect summary for shear stress at the incipient flow predicted by the full factorial design
Fig.4  The simulated shear stress at the incipient flow by the full factorial design.
Fig.5  Four different flow regimes extracted from the full factorial design. (a) Low particle friction and low cohesion; (b) high particle cohesion and low cohesion; (c) low particle friction and high cohesion; (d) high particle friction and high cohesion.
Pattern Case ID Friction coefficient Surface energy/(J·m–2)
1 0.1 0.0001
0a 2 0.45 0.0001
+− 3 0.8 0.0001
a0 4 0.1 0.05005
00 5 0.45 0.05005
A0 7 0.8 0.05005
−+ 8 0.1 0.1
0A 9 0.45 0.1
++ 10 0.8 0.1
Tab.5  Design table of the response surface design used in this study
Fig.6  The surface plots predicted by the response surface design (a) response surface of pre-shear stress by varying particle friction and surface energy (b) response surface of shear stress at incipient flow by varying particle friction and surface energy. The black mesh planes represent the experimental measurement values.
Fig.7  The actual and predicted stress fitted by the response surface design (a) shear stress at pre-shear (b) shear stress at incipient flow. The line of fit is solid red and the confidence bands (95%) are shaded red. The dashed horizontal blue line is the mean stress of the simulation cases.
Fig.8  The prediction profiler from the fitting of response surface design results.
Fig.9  (a) Validation of the pre-shear shear stress and shear stress at the incipient flow predicted from response surface design; (b) test of the validity of the calibrated material properties for the predictions of shear stress at the incipient flow under different normal stresses.
Fig.10  Comparisons of the shear stresses from experimental measurement and DEM simulation using calibrated material properties under a 5 kPa pre-consolidation normal stress.
1 M Sakai. How should the discrete element method be applied in industrial systems?: a review. Kona Powder and Particle Journal, 2016, 33: 169–178
https://doi.org/10.14356/kona.2016023
2 Y Guo, J S Curtis. Discrete element method simulations for complex granular flows. Annual Review of Fluid Mechanics, 2015, 47(1): 21–46
https://doi.org/10.1146/annurev-fluid-010814-014644
3 H Zhu, Z Zhou, R Yang, A Yu. Discrete particle simulation of particulate systems: a review of major applications and findings. Chemical Engineering Science, 2008, 63(23): 5728–5770
https://doi.org/10.1016/j.ces.2008.08.006
4 A Jarray, M Habibi, B J Scheper, H Shi, S Luding. Mixing of bidisperse cohesive granular materials in food processes. International Journal of Food Engineering, 2019, 5(3): 195–199
https://doi.org/10.18178/ijfe.5.3.195-199
5 C J Coetzee. Review: calibration of the discrete element method. Powder Technology, 2017, 310: 104–142
https://doi.org/10.1016/j.powtec.2017.01.015
6 M Marigo, E H Stitt. Discrete element method (DEM) for industrial applications: comments on calibration and validation for the modelling of cylindrical pellets. Kona Powder and Particle Journal, 2015, 32(0): 236–252
https://doi.org/10.14356/kona.2015016
7 T Roessler, A Katterfeld. Scaling of the angle of repose test and its influence on the calibration of DEM parameters using upscaled particles. Powder Technology, 2018, 330: 58–66
https://doi.org/10.1016/j.powtec.2018.01.044
8 M Combarros, H Feise, H Zetzener, A Kwade. Segregation of particulate solids: experiments and DEM simulations. Particuology, 2014, 12: 25–32
https://doi.org/10.1016/j.partic.2013.04.005
9 Z Yan, S Wilkinson, E Stitt, M Marigo. Discrete element modelling (DEM) input parameters: understanding their impact on model predictions using statistical analysis. Computational Particle Mechanics, 2015, 2(3): 283–299
https://doi.org/10.1007/s40571-015-0056-5
10 M Alizadeh, M Asachi, M Ghadiri, A Bayly, A Hassanpour. A methodology for calibration of DEM input parameters in simulation of segregation of powder mixtures, a special focus on adhesion. Powder Technology, 2018, 339: 789–800
https://doi.org/10.1016/j.powtec.2018.08.028
11 C Wensrich, A Katterfeld. Rolling friction as a technique for modelling particle shape in DEM. Powder Technology, 2012, 217: 409–417
https://doi.org/10.1016/j.powtec.2011.10.057
12 T A Simons, R Weiler, S Strege, S Bensmann, M Schilling, A Kwade. A ring shear tester as calibration experiment for DEM simulations in agitated mixers—a sensitivity study. Procedia Engineering, 2015, 102: 741–748
https://doi.org/10.1016/j.proeng.2015.01.178
13 X Bednarek, S Martin, A Ndiaye, V Peres, O Bonnefoy. Calibration of DEM parameters on shear test experiments using Kriging method. Powders and Grains 2017. Les Ulis: EDP Sciences, 2017, 15016
14 D Schulze. Powders and Bulk Solids: Behavior, Characterization, Storage and Flow. 1st ed. Berlin Heidelberg: Springer, 2008, 35–74
15 H Huang, E Tutumluer. Discrete element modeling for fouled railroad ballast. Construction & Building Materials, 2011, 25(8): 3306–3312
https://doi.org/10.1016/j.conbuildmat.2011.03.019
16 M Rackl, K J Hanley. A methodical calibration procedure for discrete element models. Powder Technology, 2017, 307: 73–83
https://doi.org/10.1016/j.powtec.2016.11.048
17 S D Figard. The basics of experimental design for multivariate analysis. In: SAS Global Forum. Buckinghamshire: SAS Institute Inc., 2009, 284
18 K J Hanley, C O’Sullivan, J C Oliveira, K Cronin, E P Byrne. Application of Taguchi methods to DEM calibration of bonded agglomerates. Powder Technology, 2011, 210(3): 230–240
https://doi.org/10.1016/j.powtec.2011.03.023
19 M Johnstone, J Ooi. Calibration of DEM models using rotating drum and confined compression measurements. In: 6th World Congress on Particle Technology. Nuremberg, 2010
20 C Kloss, C Goniva, A Hager, S Amberger, S Pirker. Models, algorithms and validation for opensource DEM and CFD-DEM. Progress in Computational Fluid Dynamics, an International Journal, 2012, 12(2-3): 140–152
21 X Chen, J A Elliott. On the scaling law of JKR contact model for coarse-grained cohesive particles. Chemical Engineering Science, 2020, 227: 115906
https://doi.org/10.1016/j.ces.2020.115906
22 K Johnson, K Kendall, A Roberts. Surface energy and the contact of elastic solids. In: Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. London: The Royal Society, 1971, 301–313
23 J M P Version. 15. Buckinghamshire: SAS Institute Inc., 2019
24 M Ajmal, T Roessler, C Richter, A Katterfeld. Calibration of cohesive DEM parameters under rapid flow conditions and low consolidation stresses. Powder Technology, 2020, 374: 22–32
https://doi.org/10.1016/j.powtec.2020.07.017
25 J Pachón-Morales, H Do, J Colin, F Puel, P Perre, D Schott. DEM modelling for flow of cohesive lignocellulosic biomass powders: model calibration using bulk tests. Advanced Powder Technology, 2019, 30(4): 732–750
https://doi.org/10.1016/j.apt.2019.01.003
26 ASTM International. Standard Test Method for Bulk Solids Using Schulze Ring Shear Tester. West Conshohocken, PA: ASTM International, 2016
27 M Proust. Design of Experiments Guide. 1st ed. Cary, NC: SAS Institute Inc., 2010, 298–300
28 G E Box, W G Hunter, J S Hunter. Statistics for Experimenters. 1st ed. New York: John Wiley and Sons, 1978, 173–200
29 P A Johnson, H Savage, M Knuth, J Gomberg, C Marone. Effects of acoustic waves on stick-slip in granular media and implications for earthquakes. Nature, 2008, 451(7174): 57–60
https://doi.org/10.1038/nature06440
30 O Dorostkar, J Carmeliet. Potential energy as metric for understanding stick-slip dynamics in sheared granular fault gouge: a coupled CFD-DEM study. Rock Mechanics and Rock Engineering, 2018, 51(10): 3281–3294
https://doi.org/10.1007/s00603-018-1457-6
31 C K Lieou, E G Daub, R A Guyer, R E Ecke, C Marone, P A Johnson. Simulating stick‐slip failure in a sheared granular layer using a physics‐based constitutive model. Journal of Geophysical Research. Solid Earth, 2017, 122(1): 295–307
https://doi.org/10.1002/2016JB013627
32 D N M Klaumünzer. Stick-slip shear banding in metallic glasses. Dissertation for the Doctoral Degree. Zurich: Eidgenössische Technische Hochschule Zürich, 2012, 25–56
33 R Cabiscol, J Finke, H Zetzener, A Kwade. Characterization of mechanical property distributions on tablet surfaces. Pharmaceutics, 2018, 10(4): 184
https://doi.org/10.3390/pharmaceutics10040184
34 Y Guo, K Buettner, V Lane, C Wassgren, W Ketterhagen, B Hancock, J Curtis. Computational and experimental studies of flexible fiber flows in a normal‐stress‐fixed shear cell. AIChE Journal. American Institute of Chemical Engineers, 2019, 65(1): 64–74
https://doi.org/10.1002/aic.16397
35 A Aigner, S Schneiderbauer, C Kloss, S Pirker. Determining the coefficient of friction by shear tester simulation. In: Proceedings of the 3rd International Conference on Particle-Based Methods. Barcelona: CIMNE, 2013, 335–342
36 W R Ketterhagen, J S Curtis, C R Wassgren, B C Hancock. Predicting the flow mode from hoppers using the discrete element method. Powder Technology, 2009, 195(1): 1–10
https://doi.org/10.1016/j.powtec.2009.05.002
[1] ZHAO Yongzhi, JIANG Maoqiang, CHENG Yi. Particle-scale simulation of fluidized bed with immersed tubes[J]. Front. Chem. Sci. Eng., 2008, 2(3): 341-345.
[2] ZHAO Yongzhi, JIANG Maoqiang, CHENG Yi, JIN Yong. Numerical simulation of two-dimensional spouted bed with draft plates by discrete element method[J]. Front. Chem. Sci. Eng., 2008, 2(1): 5-9.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed