Linear stability and isotropy properties of athermal regularized lattice Boltzmann methods

Lightweight lattice Boltzmann

A regularized version of the lattice Boltzmann method for efficient simulation of soft materials is introduced. Unlike standard approaches, this method reconstructs the distribution functions from available hydrodynamic variables (density, momentum, and pressure tensor) without storing the full set of discrete populations. This scheme shows significantly lower memory requirements and data access costs. A series of benchmark tests of relevance to soft matter, such as collisions of fluid droplets, is discussed to validate the method. The results can be of particular interest for high-performance simulations of soft matter systems on future exascale computers.

Discrete-velocity Boltzmann model: Regularization and linear stability

A discrete-velocity Boltzmann model for a nine-velocity lattice is considered. Compared to the conventional lattice Boltzmann (LB) schemes the collisions for the model are defined explicitly. Space and time discretization of the model is based on the collide and stream method; in addition, the regularization of the collision term is proposed. It is demonstrated that the regularized model can be represented as a two-relaxation-time LB model of a special type. The scheme is compared to the Onsager regularized (a specific filtered Galilean invariant model) and recursively regularized LB equations in terms of stability and dissipation properties, and linear stability analysis is performed. Several numerical experiments are carried out: double shear layer, lid-driven cavity flow, and propagation of acoustic and shear waves are considered for different grid resolutions, Mach and Reynolds numbers. It is shown that free parameters in the model corresponding to collision cross sections can be adjusted in such a way that the dissipation and stability properties can be controlled.

Linear and brute force stability of orthogonal moment multiple-relaxation-time lattice Boltzmann methods applied to homogeneous isotropic turbulence

Multiple-relaxation-time (MRT) lattice Boltzmann methods (LBM) based on orthogonal moments exhibit lattice Mach number dependent instabilities in diffusive scaling. The present work renders an explicit formulation of stability sets for orthogonal moment MRT LBM. The stability sets are defined via the spectral radius of linearized amplification matrices of the MRT collision operator with variable relaxation frequencies. Numerical investigations are carried out for the three-dimensional Taylor-Green vortex benchmark at Reynolds number 1600. Extensive brute force computations of specific relaxation frequency ranges for the full test case are opposed to the von Neumann stability set prediction. Based on that, we prove numerically that a scan over the full wave space, including scaled mean flow variations, is required to draw conclusions on the overall stability of LBM in turbulent flow simulations. Furthermore, the von Neumann results show that a grid dependence is hardly possible to include in the notion of linear stability for LBM. Lastly, via brute force stability investigations based on empirical data from a total number of 22 696 simulations, the existence of a deterministic influence of the grid resolution is deduced. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Cross-platform programming model for many-core lattice Boltzmann simulations

We present a novel, hardware-agnostic implementation strategy for lattice Boltzmann (LB) simulations, which yields massive performance on homogeneous and heterogeneous many-core platforms. Based solely on C++17 Parallel Algorithms, our approach does not rely on any language extensions, external libraries, vendor-specific code annotations, or pre-compilation steps. Thanks in particular to a recently proposed GPU back-end to C++17 Parallel Algorithms, it is shown that a single code can compile and reach state-of-the-art performance on both many-core CPU and GPU environments for the solution of a given non trivial fluid dynamics problem. The proposed strategy is tested with six different, commonly used implementation schemes to test the performance impact of memory access patterns on different platforms. Nine different LB collision models are included in the tests and exhibit good performance, demonstrating the versatility of our parallel approach. This work shows that it is less than ever necessary to draw a distinction between research and production software, as a concise and generic LB implementation yields performances comparable to those achievable in a hardware specific programming language. The results also highlight the gains of performance achieved by modern many-core CPUs and their apparent capability to narrow the gap with the traditionally massively faster GPU platforms. All code is made available to the community in form of the open-source project stlbm, which serves both as a stand-alone simulation software and as a collection of reusable patterns for the acceleration of pre-existing LB codes.