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.     

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.     

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.     

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.     

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.     

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.     

Continuous and Discrete Adjoint Approach Based on Lattice Boltzmann Method in Aerodynamic Optimization Part I: Mathematical Derivation of Adjoint Lattice Boltzmann Equations

The significance of flow optimization utilizing the lattice Boltzmann (LB) method becomes obvious regarding its advantages as a novel flow field solution method compared to the other conventional computational fluid dynamics techniques. These unique characteristics of the LB method form the main idea of its application to optimization problems. In this research, for the first time, both continuous and discrete adjoint equations were extracted based on the LB method using a general procedure with low implementation cost. The proposed approach could be performed similarly for any optimization problem with the corresponding cost function and design variables vector. Moreover, this approach was not limited to flow fields and could be employed for steady as well as unsteady flows. Initially, the continuous and discrete adjoint LB equations and the cost function gradient vector were derived mathematically in detail using the continuous and discrete LB equations in space and time, respectively. Meanwhile, new adjoint concepts in lattice space were introduced. Finally, the analytical evaluation of the adjoint distribution functions and the cost function gradients was carried out.     

Simulation of Combustion Field with Lattice Boltzmann Method

Turbulent combustion is ubiquitously used in practical combustion devices. However, even chemically non-reacting turbulent flows are complex phenomena, and chemical reactions make the problem even more complicated. Due to the limitation of the computational costs, conventional numerical methods are impractical in carrying out direct 3D numerical simulations at high Reynolds numbers with detailed chemistry. Recently, the lattice Boltzmann method has emerged as an efficient alternative for numerical simulation of complex flows. Compared with conventional methods, the lattice Boltzmann scheme is simple and easy for parallel computing. In this study, we present a lattice Boltzmann model for simulation of combustion, which includes reaction, diffusion, and convection. We assume the chemical reaction does not affect the flow field. Flow, temperature, and concentration fields are decoupled and solved separately. As a preliminary simulation, we study the so-called counter-flow laminar flame. The particular flow geometry has two opposed uniform combustible jets which form a stagnation flow. The results are compared with those obtained from solving Navier–Stokes equations.     

Incorporation of smooth spherical bodies in the Lattice Boltzmann method

A method to simulate bodies suspended in a Lattice Boltzmann solvent is proposed. It is based on a generalized reaction force that enforces no-slip boundary conditions at the fluid–body interface as the limiting case of an iterative procedure. A smooth version of the Heaviside function allows to treat spherical particles of arbitrary size and produces smooth hydrodynamic forces as particles move in the continuum. Numerical tests demonstrate the accuracy of the method in reproducing the hydrodynamic field around a single particle and the fluid-mediated forces between pairs of particles. The drag force experienced by a particle moving in a straight channel and at various Reynolds numbers is studied as a non-trivial testcase.