Compressibility in lattice Boltzmann on standard stencils: effects of deviation from reference temperature

Entropic equilibrium for the lattice Boltzmann method: Hydrodynamics and numerical properties

The entropic lattice Boltzmann framework proposed the construction of the equilibrium by taking into consideration minimization of a discrete entropy functional. The effect of this entropic equilibrium on properties of the resulting solver has been the topic of discussions in the literature. Here we present a rigorous analysis of the hydrodynamics and numerics of the entropic equilibrium. We demonstrate that the entropic equilibrium features unconditional linear stability, in contrast to the conventional polynomial equilibrium. We reveal the mechanisms through which unconditional linear stability is maintained, most notable of which are adaptive propagation velocity of normal modes and the positive-definite nature of the dissipation rates of hydrodynamic eigenmodes. We further present a simple local correction to considerably reduce the deviations in the effective bulk viscosity.

Thermal Lattice Boltzmann Flux Solver for Natural Convection of Nanofluid in a Square Enclosure

In the present study, mathematical modeling was performed to simulate natural convection of a nanofluid in a square enclosure using the thermal lattice Boltzmann flux solver (TLBFS). Firstly, natural convection in a square enclosure, filled with pure fluid (air and water), was investigated to validate the accuracy and performance of the method. Then, influences of the Rayleigh number, of nanoparticle volume fraction on streamlines, isotherms and average Nusselt number were studied. The numerical results illustrated that heat transfer was enhanced with the augmentation of Rayleigh number and nanoparticle volume fraction. There was a linear relationship between the average Nusselt number and solid volume fraction. and there was an exponential relationship between the average Nusselt number and Ra. In view of the Cartesian grid used by the immersed boundary method and lattice model, the immersed boundary method was chosen to treat the no-slip boundary condition of the flow field, and the Dirichlet boundary condition of the temperature field, to facilitate natural convection around a bluff body in a square enclosure. The presented numerical algorithm and code implementation were validated by means of numerical examples of natural convection between a concentric circular cylinder and a square enclosure at different aspect ratios. Numerical simulations were conducted for natural convection around a cylinder and square in an enclosure. The results illustrated that nanoparticles enhance heat transfer in higher Rayleigh number, and the heat transfer of the inner cylinder is stronger than that of the square at the same perimeter.

Simulation of the FDA nozzle benchmark: A lattice Boltzmann study

BACKGROUND AND OBJECTIVE

Contrary to flows in small intracranial vessels, many blood flow configurations such as those found in aortic vessels and aneurysms involve larger Reynolds numbers and, therefore, transitional or turbulent conditions. Dealing with such systems require both robust and efficient numerical methods.

METHODS

We assess here the performance of a lattice Boltzmann solver with full Hermite expansion of the equilibrium and central Hermite moments collision operator at higher Reynolds numbers, especially for under-resolved simulations. To that end the food and drug administration's benchmark nozzle is considered at three different Reynolds numbers covering all regimes: (1) laminar at a Reynolds number of 500, (2) transitional at a Reynolds number of 3500, and (3) low-level turbulence at a Reynolds number of 6500.

RESULTS

The lattice Boltzmann results are compared with previously published inter-laboratory experimental data obtained by particle image velocimetry. Our results show good agreement with the experimental measurements throughout the nozzle, demonstrating the good performance of the solver even in under-resolved simulations.

CONCLUSION

In this manner, fast but sufficiently accurate numerical predictions can be achieved for flow configurations of practical interest regarding medical applications.

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.

Extended Lattice Boltzmann Model

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.

Linear stability and isotropy properties of athermal regularized lattice Boltzmann methods

The present work proposes a general methodology to study stability and isotropy properties of lattice Boltzmann (LB) schemes. As a first investigation, such a methodology is applied to better understand these properties in the context of regularized approaches. To this extent, linear stability analyses of two-dimensional models are proposed: the standard Bhatnagar-Gross-Krook collision model, the original precollision regularization, and the recursive regularized model, where off-equilibrium distributions are partially computed thanks to a recursive formula. A systematic identification of the physical content carried by each LB mode is done by analyzing the eigenvectors of the linear systems. Stability results are then numerically confirmed by performing simulations of shear and acoustic waves. This work allows drawing fair conclusions on the stability properties of each model. In particular, among the aforementioned models, recursive regularization turns out to be the most stable one for the D2Q9 lattice, especially in the zero-viscosity limit. Two major properties shared by every regularized model are highlighted: (1) a mode filtering property and (2) an incorrect, and broadly anisotropic, dissipation rate of the modes carrying physical waves in under-resolved conditions. The first property is the main source of increased stability, especially for the recursive regularization. It is a direct consequence of the reconstruction of off-equilibrium populations before each collision process, decreasing the rank of the system of discrete equations. The second property seems to be related to numerical errors directly induced by the equilibration of high-order moments. In such a case, this property is likely to occur with any collision model that follows such a stabilization methodology.

Efficient supersonic flow simulations using lattice Boltzmann methods based on numerical equilibria

A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on a numerical equilibrium that has been extensively used by discrete velocity methods since the late 1990s. Here, it is extended to reproduce an arbitrary number of moments of the Maxwell-Boltzmann distribution. These extensions to the standard 5-constraint (mass, momentum and energy) approach lead to the correct simulation of thermal, compressible flows with only 39 discrete velocities in 3D. The stability of this BGK-LBM is reinforced by relying on Knudsen-number-dependent relaxation times that are computed analytically. Hence, high Reynolds-number, supersonic flows can be simulated in an efficient and elegant manner. While the 1D Riemann problem shows the ability of the proposed approach to handle discontinuities in the zero-viscosity limit, the simulation of the supersonic flow past a NACA0012 aerofoil confirms the excellent behaviour of this model in a low-viscosity and supersonic regime. The flow past a sphere is further simulated to investigate the 3D behaviour of our model in the low-viscosity supersonic regime. The proposed model is shown to be substantially more efficient than the previous 5-moment D3Q343 DDF-LBM for both CPU and GPU architectures. It then opens up a whole new world of compressible flow applications that can be realistically tackled with a purely LB approach. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

Impact of collision models on the physical properties and the stability of lattice Boltzmann methods

The lattice Boltzmann method (LBM) is known to suffer from stability issues when the collision model relies on the BGK approximation, especially in the zero viscosity limit and for non-vanishing Mach numbers. To tackle this problem, two kinds of solutions were proposed in the literature. They consist in changing either the numerical discretization (finite-volume, finite-difference, spectral-element, etc.) of the discrete velocity Boltzmann equation (DVBE), or the collision model. In this work, the latter solution is investigated in detail. More precisely, we propose a comprehensive comparison of (static relaxation time based) collision models, in terms of stability, and with preliminary results on their accuracy, for the simulation of isothermal high-Reynolds number flows in the (weakly) compressible regime. It starts by investigating the possible impact of collision models on the macroscopic behaviour of stream-and-collide based D2Q9-LBMs, which clarifies the exact physical properties of collision models on LBMs. It is followed by extensive linear and numerical stability analyses, supplemented with an accuracy study based on the transport of vortical structures over long distances. In order to draw conclusions as generally as possible, the most common moment spaces (raw, central, Hermite, central Hermite and cumulant), as well as regularized approaches, are considered for the comparative studies. LBMs based on dynamic collision mechanisms (entropic collision, subgrid-scale models, explicit filtering, etc.) are also briefly discussed. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.