Please use this identifier to cite or link to this item:
https://hdl.handle.net/2440/126390
Citations | ||
Scopus | Web of Science® | Altmetric |
---|---|---|
?
|
?
|
Type: | Journal article |
Title: | A scalable parallel computing SPH framework for predictions of geophysical granular flows |
Author: | Yang, E. Bui, H.H. De Sterck, H. Nguyen, G.D. Bouazza, A. |
Citation: | Computers and Geotechnics, 2020; 121:1-22 |
Publisher: | Elsevier |
Issue Date: | 2020 |
ISSN: | 0266-352X 1873-7633 |
Statement of Responsibility: | Edward Yang, Ha H.Bui, Hans De Sterck, Giang D.Nguyen, Abdelmalek Bouazza |
Abstract: | This paper presents a parallel computing Smoothed Particle Hydrodynamics (SPH) framework for geophysical granular flows scalable on large CPU clusters. The framework is accomplished by adopting a Message Passing Interface (MPI) approach with domain partitioning strategy. The Orthogonal Recursive Bisection (ORB) technique is utilised to subdivide the computational domain. The ORB algorithm is implemented such that any number of MPI processes can be used instead of being limited to powers of two. To avoid global communications in the particle distribution process, a diffusion-based distribution algorithm is implemented and demonstrated to be much faster than global communication approaches when distributing particles to non-neighbouring processes. The proposed parallel scheme achieves 95% weak scaling efficiency and up to 900 times strong scaling speedup on 1024 CPU cores. The parallel scheme enables previously unfeasible simulations to be carried out and here we apply it to the investigation of the granular column collapse experiment under full three-dimensional, axisymmetric conditions for aspect ratios up to 30, not attempted previously using numerical techniques in the literature. Enabled by the parallel scheme, the simulations use up to 11.7 million SPH particles. The investigation is conducted using two popular constitutive models commonly used in modelling of granular flows: the elasto-plastic model with Drucker-Prager yield criterion and the μ (I) rheological model. While very good agreement with experimental data has been reported for both models for small and intermediate aspect ratios, the large-scale simulations conducted for large aspect ratios show that the Drucker-Prager model tends to over-predict final deposit height, and the μ (I) model under-predicts it. Furthermore, due to the capability of the parallel scheme to model the 3D axisymmetric column collapse at higher resolutions, we demonstrate that the elasto-plastic approach is capable of capturing arching effects in the stress profile, whereas the μ (I) model cannot. |
Keywords: | Parallel computing; Message Passing Interface (MPI); Smoothed Particle Hydrodynamics (SPH); granular flows; geophysical flows |
Rights: | © 2020 Elsevier Ltd. All rights reserved. |
DOI: | 10.1016/j.compgeo.2020.103474 |
Grant ID: | http://purl.org/au-research/grants/arc/DP160100775 http://purl.org/au-research/grants/arc/DP170103793 http://purl.org/au-research/grants/arc/DP190102779 |
Published version: | http://dx.doi.org/10.1016/j.compgeo.2020.103474 |
Appears in Collections: | Aurora harvest 8 Geology & Geophysics publications |
Files in This Item:
There are no files associated with this item.
Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.