HPC performance study of different collision models using the Lattice Boltzmann solver Musubi. (15th April 2023)

• Record Type: Journal Article
• Authors: Spinelli, Gregorio Gerardo; Horstmann, Tobias; Masilamani, Kannan; Soni, Malav Mukesh; Klimach, Harald; Stück, Arthur; Roller, Sabine
• Is Part Of: Computers & fluids. Volume 255(2023)
• Publication Date: 2023-04-15
• Subjects: Lattice Boltzmann method -- Hybrid recursive regularized BGK -- Cumulant -- Turbulence modeling -- Musubi; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2023.105833
• Languages: English
• ISSNs: 0045-7930

Abstract

Abstract: Over the past decades, the lattice Boltzmann method (LBM) has become increasingly popular thanks to its capabilities in the domain of Large-Eddy Simulations (LES). Different collision schemes have been proposed to extend the scope of application to higher Reynolds number flows. This study compares the accuracy and the performance of some of these schemes on a D 3 Q 27 lattice, including the original Multiple Relaxation Times (MRT) model, the Hybrid Recursive Regularized Bhatnagar–Gross–Krook (HRR) operator, as well as the Projected Recursive Regularized Bhatnagar–Gross–Krook (PRR) operator and the parametrized Cumulant collision scheme. For this purpose, the above-mentioned schemes are implemented in the HPC LBM solver Musubi and tested on a well-documented test-case describing the flow past a circular cylinder at a Reynolds number of 3900. Three different subgrid scale (SGS) models are used to account for the unresolved turbulence, i.e. the Smagorinsky model, the Wall-Adapting Local Eddy-viscosity (WALE) model, and the Vreman model. The Cumulant scheme uses an Implicit LES (ILES) subgrid scale model and shows the best agreement with the experimental data followed by MRT with WALE, and HRR with Vreman. The examined collision models are able to capture the second peak at f = 3 f vs of the power spectra density of the y velocity component first discovered in experiments. With respect to performance, the collision models are compared in terms of MLUPs/node and parallel efficiency for a strong scaling analysis. Again the Cumulant scheme outperforms the other collision models even when they are run on the reduced D 3 Q 19 stencil. All the collision schemes show a strong scaling parallel efficiency above 60% on up to 16384 cores in our implementation. Highlights: Comparison of advanced collision models: Cumulant, HRR, MRT and PRR. Comparison of SGS turbulence models: Smagorinsky, Vreman and WALE. Comparison of LES numerical with experimental results of a Re 3900 Cylinder. Computational performance of the collision models: strong scaling and MLUPs/node. Implementation in the opensource solver Musubi.