A Parallel Algorithm for Super-Quadric Discrete Elements Based on the CUDA-GPU Architecture
-
摘要: 大规模离散元的并行计算通常基于理想的球体单元,然而自然界或工业生产中普遍存在的是由非球形颗粒组成的复杂体系,其在不同空间尺度下的动力学行为及力学性质与球形颗粒具有显著差异.基于连续函数包络的超二次曲面单元能有效地构造非球形颗粒的几何形态,并通过非线性Newton迭代算法准确计算单元间的作用力.针对非球形颗粒间接触判断的复杂性及其大规模离散元计算的需求,该文发展了基于CUDA-GPU构架下超二次曲面单元并行算法.该方法在球形颗粒并行计算的基础上,通过核函数建立单元包围盒的粗判断列表及Newton迭代的细判断列表,并优化并行算法和内存访问模式以提高算法的计算效率.为检验超二次曲面并行算法的可靠性,对非球形颗粒的流动过程进行离散元模拟, 并与试验结果进行对比验证.在此基础上,进一步分析了颗粒单元不同长宽比和表面尖锐度对颗粒材料流动特性的影响,为非球形颗粒材料的大规模离散元模拟提供一种有效的数值方法.
-
关键词:
- CUDA-GPU架构 /
- 超二次曲面单元 /
- 非球形颗粒 /
- 离散单元方法 /
- 并行算法
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. -
[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, MLLER 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] HHNER 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, LONARD 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, MRIAUX 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.
点击查看大图
计量
- 文章访问数: 1073
- HTML全文浏览量: 159
- PDF下载量: 333
- 被引次数: 0