New practical discrete non-spherical N-body method: Validation with the Brazil nut effect

Zeng, Xia; Wen, Tongge; Yu, Yang; Cheng, Bin Kang; Qiao, Dong

doi:10.1016/j.icarus.2022.115201

Cited by 8 publications

(5 citation statements)

References 46 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…Then, the algorithm queries the triangular facets enveloped by the leaf BVs for overlap. The corresponding algorithm is consistent with the previous CPU code and is omitted here for simplicity (see Wen et al 2023;Zeng et al 2022 for more details). Finally, the contact force and torque of a pair of intersected triangular facets, i.e., active contact elements, can be obtained through Equation (2).…”

Section: Contact Detection By Using Lbvhmentioning

confidence: 97%

“…To further check the acceleration of PolyDEM, the single-core CPU-based code is used to simulate the accretion of the 10400-particle assembly. The implementation of the CPU-based code follows the same algorithms described in Zeng et al (2022). The CPU simulations are performed using Intel(R) Xeon(R) Gold 6138 CPU @2.0 GHz.…”

Section: Simulation Of Self-gravitating Aggregationmentioning

confidence: 99%

“…In our recent study, a CPU-based code suite has been developed to solve the dynamical evolution of the nonspherical granular system (Zeng et al 2022). The particle shapes have no restrictions where the force-based contact model was adopted.…”

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Accelerating Polyhedral Discrete Element Method with CUDA

Wen,

Zeng

2023

Self Cite

View full text Add to dashboard Cite

This paper presents an efficient CUDA-based implementation of a nonspherical discrete element method where irregular particles are described by using polyhedrons. Two strategies are employed to exploit the parallelism of the numerical method. One is to perform contact detection based on the contact pair level instead of the traditional particle level. The second is to reduce the computational burden of each kernel function by allocating thread blocks reasonably. Contact detection between potential contact pairs is the most complicated, time-consuming, and essential process for the polyhedral discrete element method. The linear bounding volume hierarchies are introduced to fix this issue. The hierarchies of the bounding volume tree are organized in a spatially coherent way. Such a structure can minimize branch divergence and is very suitable for parallel implementation with GPU. Two numerical examples are presented to show the performance of the code. It is found from the scenario of two sphere collision that improving the mesh resolution of polyhedral particles can reduce the computational error while slowing down the computational speed correspondingly. A trade-off must be made between accuracy and efficiency. The other example of self-gravitating aggregation demonstrates the code is convergent, stable, and highly efficient. Particularly, with a mainstream GPU, the proposed method easily performs hundreds of times faster than the serial CPU code that does the same function.

show abstract

Section: Contact Detection By Using Lbvhmentioning

confidence: 97%

Section: Simulation Of Self-gravitating Aggregationmentioning

confidence: 99%

See 1 more Smart Citation

Accelerating Polyhedral Discrete Element Method with CUDA

Wen,

Zeng

2023

Self Cite

View full text Add to dashboard Cite

show abstract

“…航天动力学领域分析中常用的三角形网格多面体法由Werner [58] 于1994年提出, 相较于多面体数值积分近似求解方法随着计算机技术的发展, 离散元法近年来在小天体结构 [62] 和地质演化 [63] 的动力学分析中应用越来越图 5 (网络版彩图)小天体示意图及着陆参考坐标系每个颗粒的演化以及整个系统的运行特征. 离散元法由早期的硬球离散元 [64] 发展至目前的软球离散元方法 [65,66] , 且已发展出与多面体方法相结合的不规则离散元方法 [67] , 未来有望出现多面体颗粒与球形颗粒混 .…”

Section: 多面体法、离散元法和有限单元法均可用于小天unclassified

Recent development of landing dynamics over small celestial bodies

ZENG,

WEN,

2024

Sci. Sin.-Tech.

Self Cite

View full text Add to dashboard Cite

“…Furthermore, we expect to find more voids in systems of irregularly shaped grains due to the increased shear strength of aggregates versus their spherical counterparts, leading to a more porous equilibrium state when packing under uniform gravity. The BNE has been investigated several times in the past with pkdgrav (Matsumura et al 2014;Maurel et al 2016;Chujo et al 2018), and recent studies have shown how ellipsoidal shapes influence the process of boulder stranding (Zeng et al 2022). However, no investigation has yet been published using the varied geometries that we can model with our approach.…”

Section: Brazil Nut Effectmentioning

confidence: 99%

An Efficient Numerical Approach to Modeling the Effects of Particle Shape on Rubble-pile Dynamics

Marohnic,

DeMartini,

Richardson

et al. 2023

Planet. Sci. J.

View full text Add to dashboard Cite

We present an approach for the inclusion of nonspherical constituents in high-resolution N-body discrete element method (DEM) simulations. We use aggregates composed of bonded spheres to model nonspherical components. Though the method may be applied more generally, we detail our implementation in the existing N-body code pkdgrav. It has long been acknowledged that nonspherical grains confer additional shear strength and resistance to flow when compared with spheres. As a result, we expect that rubble-pile asteroids will also exhibit these properties and may behave differently than comparable rubble piles composed of idealized spheres. Since spherical particles avoid some significant technical challenges, most DEM gravity codes have used only spherical particles or have been confined to relatively low resolutions. We also discuss the work that has gone into improving performance with nonspherical grains, building on pkdgrav's existing leading-edge computational efficiency among DEM gravity codes. This allows for the addition of nonspherical shapes while maintaining the efficiencies afforded by pkdgrav's tree implementation and parallelization. As a test, we simulated the gravitational collapse of 25,000 nonspherical bodies in parallel. In this case, the efficiency improvements allowed for an increase in speed by nearly a factor of 3 when compared with the naive implementation. Without these enhancements, large runs with nonspherical components would remain prohibitively expensive. Finally, we present the results of several small-scale tests: spin-up due to the YORP effect, tidal encounters, and the Brazil nut effect. In all cases, we find that the inclusion of nonspherical constituents has a measurable impact on simulation outcomes.

show abstract

New practical discrete non-spherical N-body method: Validation with the Brazil nut effect

Cited by 8 publications

References 46 publications

Accelerating Polyhedral Discrete Element Method with CUDA

Accelerating Polyhedral Discrete Element Method with CUDA

Recent development of landing dynamics over small celestial bodies

An Efficient Numerical Approach to Modeling the Effects of Particle Shape on Rubble-pile Dynamics

Contact Info

Product

Resources

About