Efficient computation of particle-fluid and particle-particle interactions in compressible flow

TL;DR

This paper introduces an efficient, high-order accurate MPI+MPI hybrid algorithm for simulating dense, compressible particle-laden flows with particle collisions, suitable for complex geometries and large-scale HPC systems.

Contribution

It presents a novel particle collision and projection operator that enables scalable, high-order simulations of dense particle-fluid flows on high-performance computing architectures.

Findings

Achieves exact binary particle collision calculations on arbitrary core counts.

Demonstrates excellent scaling on state-of-the-art HPC systems.

Validates accuracy and efficiency through extensive benchmark tests.

Abstract

Particle collisions are the primary mechanism of inter-particle momentum and energy exchange for dense particle-laden flow. Accurate approximation of this collision operator in four-way coupled Euler-Lagrange approaches remains challenging due to the associated computational cost. Adopting a deterministic collision model and a hard-sphere approach eases time step constraints but imposes non-locality on distributed memory architectures, necessitating the inclusion of collision partners from each grid element in the vicinity. Retaining high-order accuracy and parallel efficiency also ties into the correct and compact treatment of the particle-fluid coupling, where adequate kernels are required to effectively project the work of the particles to the Eulerian grid. In this work, we present an efficient particle collision and projection operator based on an MPI+MPI hybrid approach to enable time-resolved and high-order accurate simulations of compressible, four-way coupled particle-laden flows at dense concentrations. A distinct feature of the proposed particle collision algorithm is the efficient calculation of exact binary inter-particle collisions on arbitrary core counts. Combining the particle operator with a hybrid discretization operator based on a high-order discontinuous Galerkin method and a localized low-order finite volume operator allows an accurate treatment of highly compressible particle-laden flows. The approach is extensively validated against a range of benchmark problems. Contrary to literature, the scaling properties are demonstrated on state-of-the-art high performance computing systems. Finally, the proposed algorithm is compatible with unstructured, curved high-order grids which permits the handling of complex geometries as is emphasized by application of the framework to large-scale application cases.