Adjoint lattice Boltzmann equation for parameter identification

The lattice Boltzmann equation is briefly introduced using moments to clearly separate the propagation and collision steps in the dynamics. In order to identify unknown parameters we introduce a cost function and adapt control theory to the lattice Boltzmann equation to get expressions for the derivatives of the cost function vs. parameters. This leads to an equivalent of the adjoint method with the definition of an adjoint lattice Boltzmann equation.To verify the general expressions for the derivatives, we consider two elementary situations: a linearized Poiseuille flow to show that the method can be used to optimize parameters, and a nonlinear situation in which a transverse shear wave is advected by a mean uniform flow. We indicate in the conclusion how the method can be used for more realistic situations.     

Multigrid solution of the steady-state lattice Boltzmann equation

Efficient solution strategies for the steady-state lattice Boltzmann equation are investigated. Stable iterative methods for the linearized lattice Boltzmann equation are formulated based on the linearization of the lattice Boltzmann time-stepping procedure. These are applied as relaxation methods within a linear multigrid scheme, which itself is used to drive a Newton solver for the full non-linear problem. Although the linear multigrid strategy provides rapid convergence, the cost of a linear residual evaluation is found to be substantially higher than the cost of evaluating the non-linear residual directly. Therefore, a non-linear multigrid approach is adopted, which makes use of the non-linear LBE time-stepping scheme on each grid level. Rapid convergence to steady-state is achieved by the non-linear algorithm, resulting in one or more orders of magnitude increase in solution efficiency over the LBE time-integration approach. Grid-independent convergence rates are demonstrated, although degradation with increasing Reynolds number is observed, as in the case of the original LBE time-stepping scheme. The multigrid solver is implemented in a modular fashion by calling an existing LBE time-stepping routine, and delivers the identical steady-state solution as the original LBE time-stepping approach.     

Bluff body flow simulation using lattice Boltzmann equation with multiple relaxation time

Numerical simulations using multiple-relaxation-time lattice Boltzmann model (MRT-LBM) are carried out for a long slender rigid circular cylinder in a cross flow to examine three-dimensional wake effect on the flow-induced forces. A mesh refinement technique is applied in the MRT-LBM calculation. The aim is to assess the validity and efficiency of the MRT-LBM model in three-dimensional calculation. In order to simulate the practical situation correctly, wall boundary conditions are specified at both ends of the cylinder. The aspect ratio of the slender cylinder is 16. The calculation is compared with results obtained from a finite volume method (FVM) and a lattice BGK model [Bhatnagar PL, Gross EP, Krook M. A model for collision processes in gases. 1. Small amplitude processes in charged and neutral one-component systems. Phys Rev 1954;94:511–25] with refined grid. Good agreement is obtained. It is found that the MRT-LBM is more efficient and faster in three-dimensional calculations.     

Progress in the understanding of the fluctuating lattice Boltzmann equation

We give a brief account of the development of methods to include thermal fluctuations into lattice Boltzmann algorithms. Emphasis is put on our recent work [B. Dünweg, U.D. Schiller, A.J.C. Ladd, Phys. Rev. E 76 (2007) 036704] which provides a clear understanding in terms of statistical mechanics.     

A lattice Boltzmann model for the Korteweg-de Vries equation with two conservation laws

A lattice Boltzmann model for the Korteweg-de Vries (KdV) equation is presented by using the higher-order moment method. In contrast to the previous lattice Boltzmann model to the KdV equation, our method has higher-order accuracy. Two key steps in the development of this model are the addition of a momentum conservation condition, and the construction of a correlation between the first conservation law and the second conservation law. The numerical example shows the higher-order moment method can be used to raise the truncation error of the lattice Boltzmann scheme.     

基于格子Boltzmann方法的一维Burgers方程的数值模拟

基于格子Boltzmann方法(LBM)的一维Burgers方程的数值解法,已有2-bit和4-bit模型。文中通过选择合适的离散速度模型构造出恰当的平衡态分布函数; 然后, 利用单松弛的格子Bhatnagar-Gross-Krook模型、Chapman-Enskog展开和多尺度技术, 提出了用于求解一维Burgers方程的3-bit的格子Boltzmann模型,即D1Q3模型,并进行了数值实验。实验结果表明,该方法的数值解与解析解吻合的程度很好,且误差比现有文献中的误差更小,从而验证了格子Boltzamnn模型的有效性。     

Simulation of floating bodies with the lattice Boltzmann method

This paper is devoted to the simulation of floating rigid bodies in free surface flows. For that, a lattice Boltzmann based model for liquid–gas–solid flows is presented. The approach is built upon previous work for the simulation of liquid–solid particle suspensions on the one hand, and on an interface-capturing technique for liquid–gas free surface flows on the other. The incompressible liquid flow is approximated by a lattice Boltzmann scheme, while the dynamics of the compressible gas are neglected. We show how the particle model and the interface capturing technique can be combined by a novel set of dynamic cell conversion rules. We also evaluate the behaviour of the free surface–particle interaction in simulations. One test case is the rotational stability of non-spherical rigid bodies floating on a plane water surface–a classical hydrostatic problem known from naval architecture. We show the consistency of our method in this kind of flows and obtain convergence towards the ideal solution for the heeling stability of a floating box.     

Multi-level lattice Boltzmann model on square lattice for compressible flows

A compressible lattice Boltzmann model is established on a square lattice. The model allows large variations in the mean velocity by introducing a large particle-velocity set. To maintain tractability, the support set of the equilibrium distribution is chosen to include only four directions and three particle-velocity levels in which the third level is introduced to improve the stability of the model. This simple structure of the equilibrium distribution makes the model efficient for the simulation of flows over a wide range of Mach numbers and gives it the capability of capturing shock jumps. Unlike the standard lattice Boltzmann model, the formulation eliminated the fourth-order velocity tensors, which were the source of concerns over the homogeneity of square lattices. A modified collision invariant eliminates the second-order discretization error of the fluctuation velocity in the macroscopic conservation equation from which the Navier–Stokes equation and energy equation are recovered. The model is suitable for both viscous and inviscid compressible flows with or without shocks. Two-dimensional shock-wave propagations and boundary layer flows were successfully simulated. The model can be easily extended to three-dimensional cubic lattices.     

The fractional Boltzmann transport equation

The fractional Boltzmann transport equation is derived making use of the fractional Hamilton’s equations based on the fractional actionlike variational approach. By simply defining a distribution function and inspecting its time derivative, many important results in statistical physics can be derived.     

Simulations of the dynamo effect with the lattice Boltzmann method

Simulations of the dynamo effect require the simultaneous integration of the Navier–Stokes equation and of the induction equation of electrodynamics. We present a hybrid method in which the Navier–Stokes equation is solved with a lattice Boltzmann method and the induction equation is treated with a spectral method.     

A lattice Boltzmann equation method for real fluids with the equation of state known in tabular form only in regions of liquid and vapor phases

We propose several simple interpolations of the isotherms for real fluids in the region below the binodal curve, where data concerning the equation of state is absent, especially in the thermodynamically prohibited region. All interpolations satisfy the boundary conditions at the points on the binodal curve. The Maxwell rule is also fulfilled. As an example, we construct several isotherms for real water. The data for the isotherms of water, in the liquid and vapor states, is given in tabular form. All smooth interpolations of the isotherms show similar hydrodynamic behavior of two-phase systems in LBE simulations. The reduced specific volumes of the liquid and vapor phases and the reduced pressure on the binodal curve obtained in the LBE simulations for the different interpolations agree well with the experimental data for real EOS of water. The surface tension depends on the form of the interpolation of the isotherm under the binodal curve. Hence, the value of surface tension can be varied in some range by changing the interpolation curve. Actually, our variant of the LBE method allows one to obtain the values of the liquid and vapor densities at the interface corresponding to the saturation curve of real fluids with high accuracy. At low temperatures, the large values of the liquid-to-vapor density ratio can be obtained, in accordance with the EOS of real fluids.     

Adjoint-based fluid flow control and optimisation with lattice Boltzmann methods

A lattice Boltzmann (LB) framework to solve fluid flow control and optimisation problems numerically is presented. Problems are formulated on a mesoscopic basis. In a side condition, the dynamics of a Newtonian fluid is described by a family of simplified Boltzmann-like equations, namely BGK–Boltzmann equations, which are linked to an incompressible Navier–Stokes equation. It is proposed to solve the non-linear optimisation problem by a line search algorithm. The needed derivatives are obtained by deriving the adjoint equations, referred to as adjoint BGK–Boltzmann equations. The primal equations are discretised by standard lattice Boltzmann methods (LBM) while for the adjoint equations a novel discretisation strategy is introduced. The approach follows the main ideas behind LBM and is therefore referred to as adjoint lattice Boltzmann methods (ALBM). The corresponding algorithm retains most of the basic features of LB algorithms. In particular, it enables a highly-efficient parallel implementation and thus solving large-scale fluid flow control and optimisation problems. The overall solution strategy, the derivation of a prototype adjoint BGK–Boltzmann equation, the novel ALBM and its parallel realisation as well as its validation are discussed in detail in this article. Numerical and performance results are presented for a series of steady-state distributed control problems with up to approximately 1.6 million unknown control parameters obtained on a high performance computer with up to 256 processing units.     

Optimizing lattice Boltzmann simulations for unsteady flows

We present detailed analysis of a lattice Boltzmann approach to model time-dependent Newtonian flows. The aim of this study is to find optimized simulation parameters for a desired accuracy with minimal computational time. Simulation parameters for fixed Reynolds and Womersley numbers are studied. We investigate influences from the Mach number and different boundary conditions on the accuracy and performance of the method and suggest ways to enhance the convergence behavior.     

多重网格格子Boltzmann方法的并行算法

针对复杂流动数值模拟中的格子Boltzmann方法存在计算网格量大、收敛速度慢的缺点,提出了基于三维几何边界的多重笛卡儿网格并行生成算法,并基于该网格生成方法提出了多重网格并行格子Boltzmann方法（LBM）。该方法结合不同尺度网格间的耦合计算,有效减少了计算网格量,提高了收敛速度;而且测试结果也表明该并行算法具有良好的可扩展性。     

Implementing lattice Boltzmann computation on graphics hardware

The Lattice Boltzmann Model (LBM) is a physically-based approach that simulates the microscopic movement of fluid particles by simple, identical, and local rules. We accelerate the computation of the LBM on general-purpose graphics hardware, by grouping particle packets into 2D textures and mapping the Boltzmann equations completely to the rasterization and frame buffer operations. We apply stitching and packing to further improve the performance. In addition, we propose techniques, namely range scaling and range separation, that systematically transform variables into the range required by the graphics hardware and thus prevent overflow. Our approach can be extended to acceleration of the computation of any cellular automata model.     

Some results on energy-conserving lattice Boltzmann models

We consider the problem of “energy conserving” lattice Boltzmann models. A major difficulty observed in previous studies is the coupling between the viscous and thermal waves even at moderate wave numbers. We propose a theoretical framework based on the knowledge of the partial equivalent equations of the lattice Boltzmann scheme at several orders of precision. With the help of linearized models (inviscid and dissipative advective acoustics and classical acoustics), we suggest natural sets of relations for the parameters of lattice Boltzmann schemes. The application is proposed for three two-dimensional schemes. Numerical test cases for simple linear and nonlinear waves establish that the main difficulty in the previous contributions can now be overcome.     

Optimal relaxation collisions for lattice Boltzmann methods

The optimal relaxation time of about 0.8090 has been proposed to balance the efficiency, stability, and accuracy at a given lattice size of numerical simulations with lattice Boltzmann methods. The optimal lattice size for a desired Reynolds number can be refined by reducing the Mach number for incompressible flows. The functioned polylogarithm polynomials are defined and used to express the lattice Boltzmann equations at different time scales and to analyze the impact of relaxation times and lattice sizes on truncation errors. Smaller truncation errors can be achieved when relaxation times are greater than 0.5 and less than 1.0. The steady-state lid-driven cavity flow was chosen to validate the code of lattice Boltzmann procedures. The applications of the optimal relaxation parameters numerically balance the stability, efficiency, and accuracy through Hartmann flow. The optimal relaxation time can also be used to select the initial lattice size for the channel flow over a square cylinder with a given Mach number.     

On a superconvergent lattice Boltzmann boundary scheme

In a seminal paper [20], Ginzburg and Adler (1994) analyzed the bounce-back boundary conditions for the lattice Boltzmann scheme and showed that it could be made exact to second order for the Poiseuille flow if some expressions depending upon the parameters of the method were satisfied, thus defining so-called “magic parameters”. Using the Taylor expansion method that one of us developed, we analyze a series of simple situations (1D and 2D) for diffusion and for linear fluid problems using bounce-back and “anti bounce-back” numerical boundary conditions. The result is that “magic parameters” depend upon the detailed choice of the moments and of their equilibrium values. They may also depend upon the way the flow is driven.     

Consistent initial conditions for lattice Boltzmann simulations

We propose a simple and effective iterative procedure to generate consistent initial conditions for the lattice Boltzmann equation (LBE) for incompressible flows with a given initial velocity field u₀. Using the Chapman-Enskog analysis we show that not only the proposed procedure effectively solves the Poisson equation for the pressure field p₀ corresponding to u₀, it also generates at the same time the initial values for the nonequilibrium distribution functions {f_α} in a consistent manner. This procedure is validated for the decaying Taylor–Green vortex flow in two dimensions and is shown to be particularly effective when using the generalized LBE with multiple relaxation times.