Nonexistence of H Theorem for Some Lattice Boltzmann Models

In this paper, we provide a set of sufficient conditions under which a lattice Boltzmann model does not admit an H theorem. By verifying the conditions, we prove that a number of existing lattice Boltzmann models does not admit an H theorem. These models include D2Q6, D2Q9 and D3Q15 athermal models, and D2Q16 and D3Q40 thermal (energy-conserving) models. The proof does not require the equilibria to be polynomials.     

Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer-Meshkov Instability

The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) model proposed in [Europhys. Lett. 90 (2010) 54003]. We discuss the reason why the Gram-Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 (2004) 035701(R)] is only valid in subsonic flows.The von Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer-Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.     

Gas Flow in Microchannels – A Lattice Boltzmann Method Approach

Gas flow in microchannels can often encounter tangential slip motion at the solid surface even under creeping flow conditions. To simulate low speed gas flows with Knudsen numbers extending into the transition regime, alternative methods to both the Navier–Stokes and direct simulation Monte Carlo approaches are needed that balance computational efficiency and simulation accuracy. The lattice Boltzmann method offers an approach that is particularly suitable for mesoscopic simulation where details of the molecular motion are not required. In this paper, the lattice Boltzmann method has been applied to gas flows with finite Knudsen number and the tangential momentum accommodation coefficient has been implemented to describe the gas-surface interactions. For fully-developed channel flows, the results of the present method are in excellent agreement with the analytical slip-flow solution of the Navier–Stokes equations, which are valid for Knudsen numbers less than 0.1. The present paper demonstrates that the lattice Boltzmann approach is a promising alternative simulation tool for the design of microfluidic devices.     

Aeroacoustic Simulations Using Compressible Lattice Boltzmann Method

This paper presents a lattice Boltzmann (LB) method based study aimed at numerical simulation of aeroacoustic phenomenon in flows around a symmetric obstacle. To simulate the compressible flow accurately, a potential energy double-distribution-function (DDF) lattice Boltzmann method is used over the entire computational domain from the near to far fields. The buffer zone and absorbing boundary condition is employed to eliminate the non-physical reflecting. Through the direct numerical simulation, the flow around a circular cylinder at $Re$=150, $M$=0.2 and the flow around a NACA0012 airfoil at $Re$=10000, $M$=0.8, $α$=$0^◦$ are investigated. The generation and propagation of the sound produced by the vortex shedding are reappeared clearly. The obtained results increase our understanding of the characteristic features of the aeroacoustic sound.     

Factorization symmetry in the lattice Boltzmann method

A non-perturbative algebraic theory of the lattice Boltzmann method is developed based on the symmetry of a product. It involves three steps: (i) Derivation of admissible lattices in one spatial dimension through a matching condition which imposes restricted extension of higher-order Gaussian moments, (ii) A special quasi-equilibrium distribution function found analytically in closed form on the product-lattice in two and three spatial dimensions, and which proves the factorization of quasi-equilibrium moments, and (iii) An algebraic method of pruning based on a one-into-one relation between groups of discrete velocities and moments. Two routes of constructing lattice Boltzmann equilibria are distinguished. The present theory includes previously known limiting and special cases of lattices, and enables automated derivation of lattice Boltzmann models from two-dimensional tables, by finding the roots of one polynomial and solving a few linear systems.     

Thermodynamic Foundations of Kinetic Theory and Lattice Boltzmann Models for Multiphase Flows

This paper demonstrates that thermodynamically consistent lattice Boltzmann models for single-component multiphase flows can be derived from a kinetic equation using both Enskog's theory for dense fluids and mean-field theory for long-range molecular interaction. The lattice Boltzmann models derived this way satisfy the correct mass, momentum, and energy conservation equations. All the thermodynamic variables in these LBM models are consistent. The strengths and weaknesses of previous lattice Boltzmann multiphase models are analyzed.     

Application of Lattice Boltzmann Method to sensitivity analysis via complex differentiation

Sensitivity analysis of flows computed by Lattice Boltzmann Method via Complex Differentiation is proposed. The theoretical work is illustrated and the proposed method assessed considering the D2Q9 scheme, along with the differentiation of the Lattice Boltzmann Equation. Boundary condition implementation is also detailed. Some examples illustrate the capability of the proposed method.     

On the Rate of Entropy Production for the Boltzmann Equation

We show that there exists a wide class of distribution functions (with moments of any order as close to their equilibrium values as we like) which can provide an abnormally low rate of entropy production. The result is valid for the Boltzmann equation with any cross section (|V|, ) satisfying a mild restriction. The functions are constructed in an explicit form and we discuss some applications of our results.     

Lattice Boltzmann simulation of solid particles suspended in fluid

The lattice Boltzmann method, an alternative approach to solving a fluid flow system, is used to analyze the dynamics of particles suspended in fluid. The interaction rule between the fluid and the suspended particles is developed for real suspensions where the particle boundaries are treated as no-slip impermeable surfaces. This method correctly and accurately determines the dynamics of single particles and multi-particles suspended in the fluid. With this method, computational time scales linearly with the number of suspensions,N, a significant advantage over other computational techniques which solve the continuum mechanics equations, where the computational time scales asN ³. Also, this method solves the full momentum equations, including the inertia terms, and therefore is not limited to low particle Reynolds number.     

Exact solutions of the nonlinear Boltzmann equation

A review is given of research activities since 1976 on the nonlinear Boltzmann equation and related equations of Boltzmann type, in which several rediscoveries have been made and several conjectures have been disproved. Subjects are (i) the BKW solution of the Boltzmann equation for Maxwell molecules, first discovered by Krupp in 1967, and the Krook-Wu conjecture concerning the universal significance of the BKW solution for the large(v, t) behavior of the velocity distribution functionf (v, t); (ii) moment equations and polynomial expansions off (v, t); (iii) model Boltzmann equation for a spatially uniform system of very hard particles, that can be solved in closed form for general initial conditions; (iv) for Maxwell and non-Maxwell-type molecules there exist solutions of the nonlinear Boltzmann equation with algebraic decrease at ; connections with nonuniqueness and violation of conservation laws; (v) conjectured super-H-theorem and the BKW solution; (vi) exactly soluble one-dimensional Boltzmann equation with spatial dependence.Reference due to C. Cercignani.     

格子Boltzmann方法中三维运动边界的统一模型

格子Boltzmann方法(LBM)中边界条件的处理很复杂,在现有的边界条件处理方法中,动力学格式能够精确满足宏观边界条件,但由于要解一个不定方程,必须引入附加假设确保方程非奇异.作为动力学格式和反弹格式的一种扩展,提出一种处理三维任意速度运动边界的统一模型,其中人口速度和固体壁面速度是该模型的特殊情形.给出用于三维15速度的表达式.为了检验该模型,模拟对角顶盖驱动三维空腔流,并将结果与有限差分法计算的结果进行比较,说明所提出的统一模型是合理可行的.     

Lattice Boltzmann method for the generalized Kuramoto-Sivashinsky equation

In this paper, a lattice Boltzmann model with an amending function is proposed for the generalized Kuramoto-Sivashinsky equation that has the form u_t+uu_x+αu_xx+βu_xxx+γu_xxxx=0. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. It is found that the numerical results agree well with the analytical solutions.     

Natural Convection in a Concentric Annulus: A Lattice Boltzmann Method Study with Boundary Condition-Enforced Immersed Boundary Method

In this paper, a boundary condition-enforced IBM is introduced into the LBM in order to satisfy the non-slip and temperature boundary conditions, and natural convections in a concentric isothermal annulus between a square outer cylinder and a circular inner cylinder are simulated. The obtained results show that the boundary condition-enforced method gives a better solution for the flow field and the complicated physics of the natural convections in the selected case is correctly captured. The calculated average Nusselt numbers agree well with the previous studies.     

辐射传输的格子Boltzmann模拟

辐射动力学理论是描述辐射传输是重要手段,基于此,本文建立了辐射能和辐射动量守恒方程,并基于Chapman-Enskog多尺度展开方法实现了从辐射传输Boltzmann方程到宏观方程的推导,进而建立了适于一维辐射传输的2分量格子Boltzmann模型。数值结果与精确解吻合较好,表明本文提出的LBM方法具有很好的准确性和稳...     

H-theorem for a linear kinetic theory

A strongH-theorem is proved for the approximate linear kinetic theory of Bawzdziewicz and Cichocki, obtained by truncating a transformed hierarchy of evolution equations. For an ith truncation we define an entropy functional that is strictly increasing in time, unless the ith reduced distribution function depends on position coordinates only. It also follows that the only stationary solution of the linear kinetic theory is the equilibrium solution. In addition, we show that the usual symmetry properties of equilibrium time correlation functions are preserved by the approximate kinetic theory under consideration.On leave of absence from Institute of Physics, Szczecin University, Wielkopolska 15, Szczecin, Poland.     

H-theorem for the (modified) nonlinear Enskog equation

We construct anH-function suitable for a system of dense hard spheres satisfying the (modified) nonlinear Enskog equation and we show that _t H 0. The equality sign holds only when the system has reached absolute equilibrium, in which caseS=– k_B H becomes the exact equilibrium entropy of the hard-sphere fluid.     

Lattice Boltzmann simulations of micron-scale drop impact on dry surfaces

A lattice Boltzmann equation (LBE) method for incompressible binary fluids is proposed to model the contact line dynamics on partially wetting surfaces. Intermolecular interactions between a wall and fluids are represented by the inclusion of the cubic wall energy in the expression of the total free energy. The proposed boundary conditions eliminate the parasitic currents in the vicinity of the contact line. The LBE method is applied to micron-scale drop impact on dry surfaces, which is commonly encountered in drop-on-demand inkjet applications. For comparison with the existing experimental results [H. Dong, W.W. Carr, D.G. Bucknall, J.F. Morris, Temporally-resolved inkjet drop impaction on surfaces, AIChE J. 53 (2007) 2606–2617], computations are performed in the range of equilibrium contact angles from 31° to 107° for a fixed density ratio of 842, viscosity ratio of 51, Ohnesorge number (Oh) of 0.015, and two Weber numbers (We) of 13 and 103.     

一个三维多松弛时间全速域格子Boltzmann模型

构建一个既适用于低速不可压流体又适用于高速可压缩流体的三维自由参数多松弛时间格子Boltzmann模型.模型中,根据SO（3）群的不可约表述基函数构造转化矩阵,根据恢复可压Navier-Stokes方程的需要选取非守恒矩平衡值.通过von Neumann稳定性分析模型参数对数值稳定性的影响,并给出建议选择范围.模型经过基准问题的验证,模拟结果与解析解及其它数值结果符合较好.