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.

Theory, Analysis, and Applications of the Entropic Lattice Boltzmann Model for Compressible Flows.

The entropic lattice Boltzmann method for the simulation of compressible flows is studied in detail and new opportunities for extending operating range are explored. We address limitations on the maximum Mach number and temperature range allowed for a given lattice. Solutions to both these problems are presented by modifying the original lattices without increasing the number of discrete velocities and without altering the numerical algorithm. In order to increase the Mach number, we employ shifted lattices while the magnitude of lattice speeds is increased in order to extend the temperature range. Accuracy and efficiency of the shifted lattices are demonstrated with simulations of the supersonic flow field around a diamond-shaped and NACA0012 airfoil, the subsonic, transonic, and supersonic flow field around the Busemann biplane, and the interaction of vortices with a planar shock wave. For the lattices with extended temperature range, the model is validated with the simulation of the Richtmyer-Meshkov instability. We also discuss some key ideas of how to reduce the number of discrete speeds in three-dimensional simulations by pruning of the higher-order lattices, and introduce a new construction of the corresponding guided equilibrium by entropy minimization.

Wetting boundaries for a ternary high-density-ratio lattice Boltzmann method.

We extend a recently proposed ternary free-energy lattice Boltzmann model with high density contrast [Phys. Rev. Lett. 120, 234501 (2018)PRLTAO0031-900710.1103/PhysRevLett.120.234501] by incorporating wetting boundaries at solid walls. The approaches are based on forcing and geometric schemes, with implementations optimized for ternary (and, more generally, higher-order multicomponent) models. Advantages and disadvantages of each method are addressed by performing both static and dynamic tests, including the capillary filling dynamics of a liquid displacing the gas phase and the self-propelled motion of a train of drops. Furthermore, we measure dynamic angles and show that the slip length critically depends on the equilibrium value of the contact angles and whether it belongs to liquid-liquid or liquid-gas interfaces. These results validate the model capabilities of simulating complex ternary fluid dynamic problems near solid boundaries, for example, drop impact solid substrates covered by a lubricant layer.

Lattice Boltzmann model for compressible flows on standard lattices: Variable Prandtl number and adiabatic exponent.

A lattice Boltzmann model for compressible flows on standard lattices is developed and analyzed. A consistent two-population thermal lattice Boltzmann is used which allows a variable Prandtl number and a variable adiabatic exponent, and appropriate correction terms are introduced into the kinetic equations to compensate for deviations in the hydrodynamic limit. Using the concept of a shifted lattice, the model is extended to supersonic flows involving shock waves, and the shock-vortex interaction problem is simulated to show the accuracy of the proposed model. Numerical results demonstrate that the proposed model is a viable candidate for compressible flow simulations.

Derivation of regularized Grad's moment system from kinetic equations: modes, ghosts and non-Markov fluxes.

Derivation of the dynamic correction to Grad's moment system from kinetic equations (regularized Grad's 13 moment system, or R13) is revisited. The R13 distribution function is found as a superposition of eight modes. Three primary modes, known from the previous derivation (Karlin et al. 1998 Phys. Rev. E57, 1668-1672. (doi:10.1103/PhysRevE.57.1668)), are extended into the nonlinear parameter domain. Three essentially nonlinear modes are identified, and two ghost modes which do not contribute to the R13 fluxes are revealed. The eight-mode structure of the R13 distribution function implies partition of R13 fluxes into two types of contributions: dissipative fluxes (both linear and nonlinear) and nonlinear streamline convective fluxes. Physical interpretation of the latter non-dissipative and non-local in time effect is discussed. A non-perturbative R13-type solution is demonstrated for a simple Lorentz scattering kinetic model. The results of this study clarify the intrinsic structure of the R13 system.This article is part of the theme issue 'Hilbert's sixth problem'.

Water ring-bouncing on repellent singularities.

Texturing a flat superhydrophobic substrate with point-like superhydrophobic macrotextures of the same repellency makes impacting water droplets take off as rings, which leads to shorter bouncing times than on a flat substrate. We investigate the contact time reduction on such elementary macrotextures through experiment and simulations. We understand the observations by decomposing the impacting drop reshaped by the defect into sub-units (or blobs) whose size is fixed by the liquid ring width. We test the blob picture by looking at the reduction of contact time for off-centered impacts and for impacts in grooves that produce liquid ribbons where the blob size is fixed by the width of the channel.

Entropic multirelaxation lattice Boltzmann models for turbulent flows.

We present three-dimensional realizations of a class of lattice Boltzmann models introduced recently by the authors [I. V. Karlin, F. Bösch, and S. S. Chikatamarla, Phys. Rev. E 90, 031302(R) (2014)] and review the role of the entropic stabilizer. Both coarse- and fine-grid simulations are addressed for the Kida vortex flow benchmark. We show that the outstanding numerical stability and performance is independent of a particular choice of the moment representation for high-Reynolds-number flows. We report accurate results for low-order moments for homogeneous isotropic decaying turbulence and second-order grid convergence for most assessed statistical quantities. It is demonstrated that all the three-dimensional lattice Boltzmann realizations considered herein converge to the familiar lattice Bhatnagar-Gross-Krook model when the resolution is increased. Moreover, thanks to the dynamic nature of the entropic stabilizer, the present model features less compressibility effects and maintains correct energy and enstrophy dissipation. The explicit and efficient nature of the present lattice Boltzmann method renders it a promising candidate for both engineering and scientific purposes for highly turbulent flows.

Matrix lattice Boltzmann reloaded.

The lattice Boltzmann equation was introduced about 20 years ago as a new paradigm for computational fluid dynamics. In this paper, we revisit the main formulation of the lattice Boltzmann collision integral (matrix model) and introduce a new two-parametric family of collision operators, which permits us to combine enhanced stability and accuracy of matrix models with the outstanding simplicity of the most popular single-relaxation time schemes. The option of the revised lattice Boltzmann equation is demonstrated through numerical simulations of a three-dimensional lid-driven cavity.

Comment on "rectangular lattice Boltzmann method".

It is shown both analytically and numerically that the suggested lattice Boltzmann model on rectangular grids [J. G. Zhou, Phys. Rev. E 81, 026705 (2010)] leads to anisotropic dissipation of fluid momentum and thus does not recover the Navier-Stokes equations. Hence, it cannot be used for the simulation of hydrodynamics.

Adaptive simplification of complex multiscale systems.

A fully adaptive methodology is developed for reducing the complexity of large dissipative systems. This represents a significant step toward extracting essential physical knowledge from complex systems, by addressing the challenging problem of a minimal number of variables needed to exactly capture the system dynamics. Accurate reduced description is achieved, by construction of a hierarchy of slow invariant manifolds, with an embarrassingly simple implementation in any dimension. The method is validated with the autoignition of the hydrogen-air mixture where a reduction to a cascade of slow invariant manifolds is observed.

Quasiequilibrium lattice Boltzmann models with tunable bulk viscosity for enhancing stability.

Taking advantage of a closed-form generalized Maxwell distribution function [P. Asinari and I. V. Karlin, Phys. Rev. E 79, 036703 (2009)] and splitting the relaxation to the equilibrium in two steps, an entropic quasiequilibrium (EQE) kinetic model is proposed for the simulation of low Mach number flows, which enjoys both the H theorem and a free-tunable parameter for controlling the bulk viscosity in such a way as to enhance numerical stability in the incompressible flow limit. Moreover, the proposed model admits a simplification based on a proper expansion in the low Mach number limit (LQE model). The lattice Boltzmann implementation of both the EQE and LQE is as simple as that of the standard lattice Bhatnagar-Gross-Krook (LBGK) method, and practical details are reported. Extensive numerical testing with the lid driven cavity flow in two dimensions is presented in order to verify the enhancement of the stability region. The proposed models achieve the same accuracy as the LBGK method with much rougher meshes, leading to an effective computational speed-up of almost three times for EQE and of more than four times for the LQE. Three-dimensional extension of EQE and LQE is also discussed.

Generalized Maxwell state and H theorem for computing fluid flows using the lattice Boltzmann method.

Generalized Maxwell distribution function is derived analytically for the lattice Boltzmann (LB) method. All the previously introduced equilibria for LB are found as special cases of the generalized Maxwellian. The generalized Maxwellian is used to derive a different class of multiple relaxation-time LB models and prove the H theorem for them.