A Lagrangian Inertial Centroidal Voronoi Particle method for dynamic load balancing in particle-based simulations
arXiv:1802.04539 · doi:10.1016/j.cpc.2019.01.011
Abstract
In this paper we develop a Lagrangian Inertial Centroidal Voronoi Particle (LICVP) method to extend the original CVP method \cite{fu2017physics} to dynamic load balancing in particle-based simulations. Two new concepts are proposed to address the additional problems encountered in repartitioning the system. First, a background velocity is introduced to transport Voronoi particle according to the local fluid field, which facilitates data reuse and lower data redistribution cost during rebalancing. Second, in order to handle problems with skew-aligned computational load and large void space, we develop an inertial-based partitioning strategy, where the inertial matrix is utilized to characterize the load distribution, and to confine the motion of Voronoi particles dynamically adapting to the physical simulation. Intensive numerical tests in fluid dynamics simulations reveal that the underlying LICVP method improves the incremental property remarkably without sacrifices on other objectives, i.e. the inter-processor communication is optimized simultaneously, and the repartitioning procedure is highly efficient.
arXiv comment31 pages, 12 figures
References in corpus (5)
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- Smoothed Particle Hydrodynamics and Magnetohydrodynamics
- Inviscid SPH
- Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver
- Examining the Accuracy of Astrophysical Disk Simulations With a Generalized Hydrodynamical Test Problem
Cited by in corpus (4)
- A Block-based Adaptive Particle Refinement SPH Method for Fluid-Structure Interaction Problems
- A Consistent Parallel Isotropic Unstructured Mesh Generation Method based on Multi-phase SPH
- A Feature-aware SPH for Isotropic Unstructured Mesh Generation
- Multi-level adaptive particle refinement method with large refinement scale ratio and new free-surface detection algorithm for complex fluid-structure interaction problems