A Parallel Algorithm for Super-Quadric Discrete Elements Based on the CUDA-GPU Architecture

WANG Siqiang; JI Shunying

doi:10.21656/1000-0887.390267

Volume 40 Issue 7

Turn off MathJax

Article Contents

Article Navigation > Applied Mathematics and Mechanics > 2019 > 40(7): 751-767

WANG Siqiang, JI Shunying. A Parallel Algorithm for Super-Quadric Discrete Elements Based on the CUDA-GPU Architecture[J]. Applied Mathematics and Mechanics, 2019, 40(7): 751-767. doi: 10.21656/1000-0887.390267

Citation:

PDF( 5367 KB)

A Parallel Algorithm for Super-Quadric Discrete Elements Based on the CUDA-GPU Architecture

doi: 10.21656/1000-0887.390267

State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology, Dalian, Liaoning 116024, P.R.China

Funds: The National Key R&D Program of China(2016YCF1401505); The National Natural Science Foundation of China(11572067;11772085)

Received Date: 2018-10-15
Rev Recd Date: 2018-10-30
Publish Date: 2019-07-01

Abstract

Abstract

The traditional parallel computing of large-scale discrete elements was suitable for spherical particles. However, in natural fields or industrial applications, the granular systems commonly comprise non-spherical particles. Meanwhile, the dynamic behaviors and mechanical properties of non-spherical particles are strongly different from those of spherical particles at different spatial scales. Super-quadric elements based on the continuous function envelope were used to effectively describe the geometric shapes of irregular particles, and accurately calculate the contact forces between elements with the non-linear Newtonian method. In view of the complexity of the contact detection between non-spherical particles and the large-scale computational requirements of the discrete element method, a CUDA-GPU parallel algorithm was developed for super-quadric elements. Based on the parallel calculation of spherical particles, the rough contact list of the element envelope and the accurate contact list of the Newtonian method were established with the kernel function. Meanwhile, the parallel algorithm and the memory access mode were optimized to improve the computation efficiency. To examine the reliability of the parallel algorithm, the flow process of non-spherical particles was simulated with the discrete element method and compared with the experimental results. Furthermore, the influences of the aspect ratio and the sharpness parameter of elements on the flow characteristics of non-spherical particles were analyzed. This study provides an effective numerical method for large-scale simulation of non-spherical granular materials.
- CUDA-GPU architecture,
- super-quadric element,
- non-spherical particle,
- discrete element method (DEM),
- parallel algorithm

FullText(HTML)

References(52)

References

[1]	季顺迎. 非均匀颗粒材料的类固-液相变行为及本构方程[J]. 力学学报, 2007,39(2): 223-237.(JI Shunying. The quasi-solid-liquid phase transition of non-uniform granular materials and their constitutive equation[J]. Chinese Journal of Theoretical and Applied Mechanics,2007,39(2): 223-237.(in Chinese))
[2]	徐泳, 孙其诚, 张凌, 等. 颗粒离散元法研究进展[J]. 力学进展, 2003,33(2): 251-260.(XU Yong, SUN Qicheng, ZHANG Ling, et al. Advances in discrete element methods for particulate materials[J]. Advances in Mechanics,2003,33(2): 251-260.(in Chinese))
[3]	CUNDALL P A, STRACK O D L. A discrete numerical mode for granular assemblies[J]. Géotech-nique,1979,29(1): 47-65.
[4]	ZHU H P, ZHOU Z Y, YANG R Y, et al. Discrete particle simulation of particulate systems: a review of major applications and findings[J]. Chemical Engineering Science,2008,62(13): 3378-3396.
[5]	常晓林, 马刚, 周伟, 等. 颗粒形状及粒间摩擦角对堆石体宏观力学行为的影响[J]. 岩土工程学报, 2012,34(4): 646-653.(CHANG Xiaolin, MA Gang, ZHOU Wei, et al. Influences of particle shape and inter-particle friction angle on macroscopic response of rockfill[J]. Chinese Journal of Geotechnical Engineering,2012,34(4): 646-653.(in Chinese))
[6]	严颖, 赵金凤, 季顺迎. 块石含量和空间分布对土石混合体抗剪强度影响的离散元分析[J]. 工程力学, 2017,34(6): 146-156.(YAN Ying, ZHAO Jinfeng, JI Shunying. Discrete element analysis of the influence of rock content and rock spatial distribution on shear strength of rock-soil mixtures[J]. Engineering Mechanics,2017,34(6): 146-156.(in Chinese))
[7]	崔泽群, 陈友川, 赵永志, 等. 基于超二次曲面的非球形离散单元模型研究[J]. 计算力学学报, 2013,30(6): 854-859.(CUI Zequn, CHEN Youchuan, ZHAO Yongzhi, et al. Study of discrete element model for non-sphere particles base on super-quadric[J]. Chinese Journal of Computational Mechanics,2013,30(6): 854-859.(in Chinese))
[8]	LU G, THIRD J R, MLLER C R. Discrete element models for non-spherical particle systems: from theoretical developments to applications[J]. Chemical Engineering Science,2015,127: 425-465.
[9]	ZHONG W, YU A, LIU X, et al. DEM/CFD-DEM modelling of non-spherical particulate systems: theoretical developments and applications[J]. Powder Technology,2016,302: 108-152.
[10]	严颖, 季顺迎. 颗粒形态对离散介质剪切强度的影响[J]. 岩土力学, 2009,30(S1): 225-230.(YAN Ying, JI Shunying. Effects of particle shape on shear strength of discrete media[J]. Rock and Soil Mechanics,2009,30(S1): 225-230.(in Chinese))
[11]	GUI N, YANG X, TU J, et al. Effect of roundness on the discharge flow of granular particles[J]. Powder Technology,2017,314: 140-147.
[12]	LIU S D, ZHOU Z Y, ZOU R P, et al. Flow characteristics and discharge rate of ellipsoidal particles in a flat bottom hopper[J]. Powder Technology,2014,253: 70-79.
[13]	ZHENG Q J, ZHOU Z Y, YU A B. Contact forces between viscoelastic ellipsoidal particles[J]. Powder Technology,2013,248: 25-33.
[14]	CLEARY P W, SINNOTT M D, MORRISON R D, et al. Analysis of cone crusher performance with changes in material properties and operating conditions using DEM[J]. Minerals Engineering,2017,100: 49-70.
[15]	CLEARY P W, HILTON J E, SINNOTTO M D. Modelling of industrial particle and multiphase flows[J]. Powder Technology,2017,314: 232-252.
[16]	刘璐, 龙雪, 季顺迎. 基于扩展多面体的离散单元法及其作用于圆桩的冰载荷计算[J]. 力学学报, 2015,47(6): 1046-1057.(LIU Lu, LONG Xue, JI Shunying. Dilated polyhedra based discrete element method and its application of ice load on cylindrical pile[J]. Chinese Journal of Theoretical and Applied Mechanics,2015,47(6): 1046-1057.(in Chinese))
[17]	孙珊珊, 严颖, 赵春发, 等. 往复荷载下铁路道砟沉降特性的扩展多面体离散元分析[J]. 铁道学报, 2015,37(11): 89-95.(SUN Shanshan, YAN Ying, ZHAO Chunfa, et al. Dilated polyhedral discrete element analysis of settlement characteristics of railway ballast under cyclic loading[J]. Journal of the China Railway Society,2015,37(11): 89-95.(in Chinese))
[18]	许文祥, 孙洪广, 陈文, 等. 软物质系颗粒材料组成、微结构与传输性能之间关联建模综述[J]. 物理学报, 2016,65(17): 75-98.(XU Wenxiang, SUN Hongguang, CHEN Wen, et al. A review of correlative modeling for transport properties, microstructures, and compositions of granular materials in soft matter[J]. Acta Physica Sinica,2016,65(17): 75-98.(in Chinese))
[19]	SU D, YAN W M. 3D characterization of general-shape sand particles using microfocus X-ray computed tomography and spherical harmonic functions, and particle regeneration using multivariate random vector[J]. Powder Technology,2018,323: 8-23.
[20]	KAWAMOTO R, AND E, VIGGIANI G, et al. All you need is shape: predicting shear banding in sand with LS-DEM[J]. Journal of the Mechanics and Physics of Solids,2018,111: 375-392.
[21]	LIM K W, KAWAMOTO R, AND E, et al. Multiscale characterization and modeling of granular materials through a computational mechanics avatar: a case study with experiment[J]. Acta Geotechnica,2015,11(2): 243-253.
[22]	王嗣强, 季顺迎. 基于超二次曲面的颗粒材料缓冲性能离散元分析[J]. 物理学报, 2018,67(9): 094501. DOI: 10.7498/aps.67.20172549.(WANG Siqiang, JI Shunying. Discrete element analysis of buffering capacity of non-spherical granular materials based on super-quadric method[J]. Acta Physica Sinica,2018,67(9): 094501. DOI: 10.7498/aps.67.20172549.(in Chinese))
[23]	DELANEY G W, CLEARY P W. The packing properties of superellipsoids[J]. Europhysics Letters,2010,89(3): 34002. DOI: 10.1209/0295-5075/89/34002.
[24]	MA H, ZHAO Y. Modelling of the flow of ellipsoidal particles in a horizontal rotating drum based on DEM simulation[J]. Chemical Engineering Science,2017,172: 636-651.
[25]	MA H, ZHAO Y. Investigating the flow of rod-like particles in a horizontal rotating drum using DEM simulation[J]. Granular Matter,2018,20(3): 20-41.
[26]	HHNER D, WIRTZ S, SCHERER V. A study on the influence of particle shape on the mechanical interactions of granular media in a hopper using the discrete element method[J]. Powder Technology,2015,278: 286-305.
[27]	AMRITKAR A, DEB S, TAFTI D. Efficient parallel CFD-DEM simulations using OpenMP[J]. Journal of Computational Physics,2014,256: 501-519.
[28]	CHEN J, MATUTTIS H G. Optimization and OpenMP parallelization of a discrete element code for convex polyhedra on multi-core machines[J]. International Journal of Modern Physics C,2013,24(2): 1350001. DOI: 10.1142/S0129183113500010.
[29]	LEMIEUX M, LONARD G, DOUCET J, et al. Large-scale numerical investigation of solids mixing in a V-blender using the discrete element method[J]. Powder Technology,2008,181(2): 205-216.
[30]	YAN B, REGUEIRO R A. A comprehensive study of MPI parallelism in three-dimensional discrete element method (DEM) simulation of complex-shaped granular particles[J]. Computational Particle Mechanics,2018,5(4): 553-577.
[31]	LIU H, TAFTI D K, LI T. Hybrid parallelism in MFIX CFD-DEM using OpenMP[J]. Powder Technology,2014,259: 22-29.
[32]	BERGER R, KLOSS C, KOHLMEYER A, et al. Hybrid parallelization of the LIGGGHTS open-source DEM code[J]. Powder Technology,2015,278: 234-247.
[33]	付晓东, 盛谦, 张勇慧. 基于OpenMP的非连续变形分析并行计算方法[J]. 岩土力学, 2014,35(8): 2401-2407.(FU Xiaodong, SHENG Qian, ZHANG Yonghui. Parallel computing method for discontinuous deformation analysis using OpenMP[J]. Rock and Soil Mechanics,2014,35(8): 2401-2407.(in Chinese))
[34]	严成增, 郑宏, 孙冠华, 等. 基于OpenMP的二维有限元-离散元并行分析方法[J]. 岩土力学, 2014,9(35): 2717-2724.(YAN Chengzeng, ZHENG Hong, SUN Guanhua, et al. Parallel analysis of two-dimensional finite-discrete element method based on OpenMP[J]. Rock and Soil Mechanics,2014,9(35): 2717-2724.(in Chinese))
[35]	WASHIZAWA T, NAKAHARA Y. Parallel computing of discrete element method on GPU[J]. Applied Mathematics,2013,4(1): 242-247.
[36]	XU J, QI H, FANG X, et al. Quasi-real-time simulation of rotating drum using discrete element method with parallel GPU computing[J]. Particuology,2011, 9(4): 446-450.
[37]	SEO I S, KIM J H, SHIN J H, et al. Particle behaviors of printing system using GPU-based discrete element method[J]. Journal of Mechanical Science and Technology,2014,28(12): 5083-5087.
[38]	LONG X, JI S, WANG Y. Validation of microparameters in discrete element modeling of sea ice failure process[J]. Particulate Science and Technology,2018,37(5): 546-555. DOI: 10.1080/02726351.2017.1404515.
[39]	QI J, LI K, JIANG H, et al. GPU-accelerated DEM implementation with CUDA[J]. International Journal of Computational Science and Engineering,2015,11(3): 330-337.
[40]	HE Y, EVANS T J, YU A B, et al. A GPU-based DEM for modelling large scale powder compaction with wide size distributions[J]. Powder Technology,2018,333: 219-228.
[41]	狄少丞, 季顺迎. 海冰与自升式海洋平台相互作用GPU离散元模拟[J]. 力学学报, 2014,46(4): 561-571.(DI Shaocheng, JI Shunying. GPU-based discrete element modelling of interaction between sea ice and jack-up platform structure[J]. Chinese Journal of Theoretical and Applied Mechanics,2014,46(4): 561-571.(in Chinese))
[42]	GAN J Q, ZHOU Z Y, YU A B. A GPU-based DEM approach for modelling of particulate systems[J]. Powder Technology,2016,301: 1172-1182.
[43]	GOVENDER N, WILKE D N, PIZETTE P, et al. A study of shape non-uniformity and poly-dispersity in hopper discharge of spherical and polyhedral particle systems using the Blaze-DEM GPU code[J]. Applied Mathematics and Computation,2018,319: 318-336.
[44]	GOVENDER N, WILKE D N, KOK S. Collision detection of convex polyhedra on the NVIDIA GPU architecture for the discrete element method[J]. Applied Mathematics and Computation,2015,267: 810-829.
[45]	SOLTANBEIGI B, PODLOZHNYUK A, PAPANICOLOPULOS S A, et al. DEM study of mechanical characteristics of multi-spherical and superquadric particles at micro and macro scales[J]. Powder Technology,2018,329: 288-303.
[46]	AH B. Superquadrics and angle-preserving transformations[J]. IEEE Computer Graphics and Applications,1981, 1(1): 11-23.
[47]	NISHIURA D, SAKAGUCHI H. Parallel-vector algorithms for particle simulations on shared-memory multiprocessors[J]. Journal of Computational Physics,2011,230: 1923-1938.
[48]	PORTAL R, DIAS J, DE SOUSA L. Contact detection between convex superquadric surfaces[J]. Archive of Mechanical Engineering,2010,57(2): 165-186.
[49]	PODLOZHNYUK A, PIRKER S, KLOSS C. Efficient implementation of superquadric particles in discrete element method within an open-source framework[J]. Computational Particle Mechanics,2016,4(1): 101-118.
[50]	ZHOU Z Y, ZOU R P, PINSON D, et al. Dynamic simulation of the packing of ellipsoidal particles[J]. Industrial and Engineering Chemistry Research,2011,50(16): 9787-9798.
[51]	GOLDMAN R. Curvature formulas for implicit curves and surfaces[J]. Computer Aided Geometric Design,2005,22(7): 632-658.
[52]	OWEN P J, CLEARY P W, MRIAUX C. Quasi-static fall of planar granular columns: comparison of 2D and 3D discrete element modelling with laboratory experiments[J]. Geomechanics and Geoengineering,2009,4(1): 55-77.