B-spline-based material point method with dynamic load balancing technique for large-scale simulation

作     者：Hidano, Soma Pan, Shaoyuan Yoshida, Keina Nomura, Reika Miki, Yohei Kawai, Masatoshi Moriguchi, Shuji Nakajima, Kengo Terada, Kenjiro 

作者机构：Tohoku Univ Grad Sch Engn Dept Civil & Environm Engn 6-6-06 AzaaobaAramakiAoba Ku Sendai 9808572 Japan Tohoku Univ Int Res Inst Disaster Sci 468-1 AzaaobaAramakiAoba Ku Sendai 9808572 Japan Univ Tokyo Informat Technol Ctr 5-1-5 Kashiwanoha Kashiwa Chiba 2778589 Japan Nagoya Univ Informat Technol Ctr Furo ChoChikusa Ku Nagoya Aichi 4648603 Japan 

出 版 物：《ENGINEERING WITH COMPUTERS》 (Eng Comput)

年 卷 期：2025年

页      面：1-18页

核心收录：

学科分类：08[工学] 0802[工学-机械工程] 0812[工学-计算机科学与技术（可授工学、理学学位）] 

基　　金：Japan Society for the Promotion of Science [23KJ0124, 22H00507, 19H01094] JSPS KAKENHI [JPMJSP2114, jh220019, jh230070, jh240040] JST SPRING 

主　　题：Material point method B-spline basis functions Domain decomposition Dynamic load balancing Large-scale landslide simulation 

摘      要：In this study, a dynamic load-balancing (DLB) technique based on the sampling method is developed for MPMs using higher-order B-spline basis functions for parallel MPI calculations based on domain decomposition, in order to perform large-scale, long-duration landslide simulations in realistic computation time. Higher-order B-spline basis functions use a range of influence across cells compared to general basis functions, but this DLB technique dynamically adjusts the size of the computational domain according to the material point distribution, so that the material points are almost equally distributed across all cores. This allows the load bias between cores to be mitigated and the advantages of parallel computation to be fully exploited. Specifically, the novel contribution of this study is that the domain decomposition allows for proper communication between control points, even if the physical regions of the cores are staggered or non-adjacent, and even if the area of influence of B-spline basis functions spans multiple subdomains at this time. In numerical examples, the quasi-3D benchmark solid column collapse problem is computed for multiple core configurations to verify the effectiveness of the DLB method in terms of scalability and parallelization efficiency. The simulation of the full 3D column collapse problem also illustrates the applicability of the proposed DLB method to large-scale disaster simulations. Finally, to demonstrate the promise and capability of the DLB technique in the MPM algorithm, a full-scale size landslide disaster simulation is carried out to illustrate that it can withstand some practical size calculations.