A two-dimensional lattice Boltzmann model for uniform channel flows

In this paper a novel two-dimensional lattice Boltzmann model (LBM) is developed for uniform channel flows. The axial velocity is solved from a momentum diffusion equation over the cross-sectional plane. An extrapolation boundary condition is also introduced to enhance the no-slip boundary in the momentum equation. This boundary treatment can also be applied to LBM simulations of other diffusion processes. The algorithm and boundary treatment are validated by simulations of steady Poiseuille and pulsatile Womersley flows in circular pipes. The numerical convergence and accuracy are comparable to those of existing models. Moreover, comparison with general three-dimensional lattice Boltzmann simulations demonstrates the advantages of our two-dimensional model, including lower computational resource requirements (memory and time), easier boundary treatment for arbitrary cross-sectional shapes, and no velocity constraint. These features are attractive for practical applications with uniform channel flows.     

Multi-relaxation-time lattice Boltzmann model for axisymmetric flows

In this paper, a lattice-Boltzmann equation (LBE) with multi relaxation times (MRT) is presented for axisymmetric flows. The model is an extension of a recent model with single-relaxation-time [Guo et al., Phys. Rev. E 79, 046708 (2009)], which was developed based on the axisymmetric Boltzmann equation. Due to the use of the MRT collision model, the present model can achieve better numerical stability. The model is validated by some numerical tests including the Hagen-Poiseuille flow, the pulsatile Womersley flow, and the external flow over a sphere. Numerical results are in excellent agreement with analytical solutions or other available data, and the improvement in numerical stability is also confirmed.     

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.     

Stable lattice Boltzmann schemes with a dual entropy approach for monodimensional nonlinear waves

We follow the mathematical framework proposed by Bouchut (2003) [22] and present in this contribution a dual entropy approach for determining equilibrium states of a lattice Boltzmann scheme. This method is expressed in terms of the dual of the mathematical entropy relative to the underlying conservation law. It appears as a good mathematical framework for establishing a “H-theorem” for the system of equations with discrete velocities. The dual entropy approach is used with D1Q3 lattice Boltzmann schemes for the Burgers equation. It conducts to the explicitation of three different equilibrium distributions of particles and induces naturally a nonlinear stability condition. Satisfactory numerical results for strong nonlinear shocks and rarefactions are presented. We prove also that the dual entropy approach can be applied with a D1Q3 lattice Boltzmann scheme for systems of linear and nonlinear acoustics and we present a numerical result with strong nonlinear waves for nonlinear acoustics. We establish also a negative result: with the present framework, the dual entropy approach cannot be used for the shallow water equations.     

An evaluation of lattice Boltzmann schemes for porous medium flow simulation

We quantitatively evaluate the capability and accuracy of the lattice Boltzmann equation (LBE) for modeling flow through porous media. In particular, we conduct a comparative study of the LBE models with the multiple-relaxation-time (MRT) and the Bhatnagar–Gross–Krook (BGK) single-relaxation-time (SRT) collision operators. We also investigate several fluid–solid boundary conditions including: (1) the standard bounce-back (SBB) scheme, (2) the linearly interpolated bounce-back (LIBB) scheme, (3) the quadratically interpolated bounce-back (QIBB) scheme, and (4) the multi-reflection (MR) scheme. Three-dimensional flow through two porous media—a body-centered cubic (BCC) array of spheres and a random-sized sphere-pack—are examined in this study. For flow past a BCC array of spheres, we validate the linear LBE model by comparing its results with the nonlinear LBE model. We investigate systematically the viscosity-dependence of the computed permeability, the discretization error, and effects due to the choice of relaxation parameters with the MRT and BGK schemes. Our results show unequivocally that the MRT–LBE model is superior to the BGK–LBE model, and interpolation significantly improves the accuracy of the fluid–solid boundary conditions.     

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

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

An improved entropic lattice Boltzmann model for parallel computation

In this paper, we suggest two kinds of approximation methods based on Taylor series expansion which can solve the non-linear equation in entropic lattice Boltzmann model without using any iteration methods such as Newton–Raphson method. The advantage of our methods is to be able to avoid the load imbalance in parallel computation which occurs due to the differences of iteration number on each calculation grid. In this study, ELBM simulations using our methods were compared with those using Newton–Raphson method for the channel flow past a square cylinder in Re = 1000 and the validity of the results and computational effort were investigated. As a result, it was found that the solutions obtained by our methods are qualitatively and quantitatively reasonable and CPU time is shorter than those obtained by Newton–Raphson method.     

A high-performance lattice Boltzmann implementation to model flow in porous media

We examine the problem of simulating single and multiphase flow in porous medium systems at the pore scale using the lattice Boltzmann (LB) method. The LB method is a powerful approach, but one which is also computationally demanding; the resolution needed to resolve fundamental phenomena at the pore scale leads to very large lattice sizes, and hence substantial computational and memory requirements that necessitate the use of massively parallel computing approaches. Common LB implementations for simulating flow in porous media store the full lattice, making parallelization straightforward but wasteful. We investigate a two-stage implementation consisting of a sparse domain decomposition stage and a simulation stage that avoids the need to store and operate on lattice points located within a solid phase. A set of five domain decomposition approaches are investigated for single and multiphase flow through both homogeneous and heterogeneous porous medium systems on different parallel computing platforms. An orthogonal recursive bisection method yields the best performance of the methods investigated, showing near linear scaling and substantially less storage and computational time than the traditional approach.     

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.     

A thermal lattice Boltzmann model with diffuse scattering boundary condition for micro thermal flows

A lattice Boltzmann model for simulating isothermal micro flows has been proposed by us recently [Niu XD, Chew YT, Shu C. A lattice Boltzmann BGK model for simulation of micro flows. Europhys Lett 2004;67(4):600]. In this paper, we extend the model to simulate the micro thermal flows. In particular, the thermal lattice Boltzmann equation (TLBE) [He X, Chen S, Doolen GD. A novel thermal model for the lattice Boltzmann method in incompressible limit. J Comput Phys 1998;146:282] is used with modification of the relaxation times linking to the Knudsen number. The diffuse scattering boundary condition (DSBC) derived in our early model is extended to consider temperature jump at wall boundaries. Simple theoretical analyses of the DSBC are presented and the results are found to be consistent with the conventional velocity slip and temperature jump boundary conditions. Numerical validations are carried out by simulating two-dimensional thermal Couette flows and developing thermal flows in a microchannel, and the obtained results are found to be in good agreement with those given from the direct simulation Monte Carlo (DSMC), the molecular dynamics (MD) approaches and the Maxwell theoretical prediction.     

Towards simulating urban canyon circulations with a 2D lattice Boltzmann model

The Lattice Boltzmann (LB) method is a novel fluid modelling technique developed from cellular automata. Instead of numerically solving the continuum Navier–Stokes equations, it simulates the interactions of mesoscopic particle populations p_α using discrete speeds and positions to obtain the macroscopic velocity, density and temperature fields. Localised at neighbouring grid nodes, the method handles complex geometries and multiple fluids more easily than traditional continuum CFD methods.Rothman and Zaleski (Lattice-Gas Cellular Automata: Simple Models of Complex Hydrodynamics (1997) Cambridge University Press, Cambridge) discuss LB method theory and development in more detail.To demonstrate the power of the technique, a 2D LB model is first used to perform urban canyon configuration studies at Reynolds number Re=100 for Height to Width (H/W) ratios from 0.125 to 2. Then, thermal lid driven cavity simulations for Re=100 and Rayleigh number Ra=2000 are performed for different locations of a relatively hot wall. The simulated flow fields appear qualitatively consistent with physical flows observed in wind tunnel and field studies, and indicate that LB methods generate results comparable to traditional CFD methods for the selected flow situations.     

基于元胞自动机各向异性扩散模型的图像分割算法

介绍了一种基于元胞自动机各向异性扩散模型的图像分割算法.在扩散模型基础上,引入跨膜介质,构建各向异性扩散模型;然后模拟热扩散方程的活动轮廓模,建立热量场,通过 LBADM 模型对图像进行分割获得目标分割边缘.实验结果表明,该算法模型能够取得闭合的分割曲线,同样能够很好的处理图像拓扑结构的变化.通过与水平集方法和窄带水平集方法进行计算速度实验对比,表明该算法能够大大减少了分割的计算量.     

A lattice Boltzmann study of gas flows in a long micro-channel

In this paper, the pressure-driven flow in a long micro-channel is studied via a lattice Boltzmann equation (LBE) method. With the inclusion of the gas–wall collision effects, the LBE is able to capture the flow behaviors in the transition regime. The numerical results are compared with available data of other methods. Furthermore, the effects of rarefaction and compressibility on the deviation of the pressure distribution from the linear one are also investigated.     

Isotropic color gradient for simulating very high-density ratios with a two-phase flow lattice Boltzmann model

This study presents the integration of isotropic color gradient discretization into a lattice Boltzmann Rothman–Keller (RK) model designed for two-phase flow simulation. The proposed model removes one limitation of the RK model, which concerns the handling of O(1000) large density ratios between the fluids for a wide range of parameters. Taylor’s series expansions are used to characterize the difference between an isotropic gradient discretization and the commonly used anisotropic gradient. The proposed color gradient discretization can reduce, by one order of magnitude, the spurious current problem that affects the interface between the phases. A set of numerical tests is conducted to show that a rotationally invariant discretization enables widening of the parameter range for the surface tension. Surface tensions from O(10⁻²) to O(10⁻⁸), depending on the density ratio, are accurately simulated. An extreme density ratio of O(10,000) is successfully tested for a steady bubble with an error of 0.5% for Laplace’s law across a sharp interface, with a thickness of about 5–6 lattice units.     

A viscosity counteracting approach in the lattice Boltzmann BGK model for low viscosity flow: Preliminary verification

Due to numerical instability, the lattice Boltzmann model (LBM) with the Bhatnagar–Gross–Krook (BGK) collision operator has some limitations in the simulation of low viscosity flows. In this paper, we propose a viscosity counteracting approach for simulating a moderate viscosity flow. An extra negative viscosity term is introduced to counteract part of the moderate viscosity by using the lattice Boltzmann equation with a source term. The counteracting viscosity term is treated as a non-uniform unsteady source. The stability is enhanced; thus small viscosity flows can be simulated. Model verification consists of benchmark cases such as those of Poiseuille flow, Couette flow, waterhammer waves, Taylor–Green vortex flow, and lid-driven cavity flow. The flow patterns, error characteristics, and representative parameters are carefully analyzed. It is shown that this approach can simulate flows with lower viscosities than may be simulated using the normal LBGK model; the second-order accuracy of the LBGK model is definitely retained, although a little dissipation is added. These preliminary studies prove the effectiveness and accuracy of the model. Sophisticated analysis and further verification of the stability mechanism will be done in the near future.     

Multi-dimensional texts in a one-dimensional medium

This paper discusses one of the tools which may be used for representing texts in machine-readable form, i.e. encoding systems or markup languages. This discussion is at the same time a report on current tendencies in the field. An attempt is made at reconstructing some of the main conceptions of text lying behind these tendencies. It is argued that, although the conceptions of texts and text structures inherent in these tendencies seem to be misguided, text encoding is nevertheless a fruitful approach to the study of texts. Finally, some conclusions are drawn concerning the relevance of this discussion to themes in text linguistics.Claus Huitfeldt studied philosophy at the University of Trondheim, writing his dissertation on the nature of transcendental arguments. He then worked for several years at the Norwegian Computing Centre for the Humanities, at the Norwegian Wittgenstein Project, and as Research Fellow in philosophy, before becoming Director of the Wittgenstein Archives at the University of Bergen. He has published a number of papers on text encoding.     

The computation of strain rate tensor in multiple-relaxation-time lattice Boltzmann model

The multiple-relaxation-time (MRT) lattice Boltzmann (LB) model is an important class of LB model with lots of advantages over the traditional single-relaxation-time (SRT) LB model. Generally, the computation of strain rate tensor is crucial for the MRT-LB simulations of some complex flows. At present, only two formulae are available to compute the strain rate tensor in the MRT LB model. One is to compute the strain rate tensor using the non-equilibrium parts of macroscopic moments (Yu formula). The other is to compute the strain rate tensor using the non-equilibrium parts of density distribution functions (Chai formula). The mathematical expressions of these two formulae are so different that we do not know which formula to choose for computing the strain rate tensor in the MRT LB model. To overcome this problem, this paper presents a theoretical study of the relationship between Chai and Yu formulae. The results show that the Yu formula can be deduced from the Chai formula, although they have their own advantages and disadvantages. In particular, the Yu formula is computationally more efficient, while the Chai formula is applicable to more lattice patterns of the MRT LB models. Furthermore, the derivation of the Yu formula in a particular lattice pattern from the Chai formula is more convenient than that proposed by Yu et al.     

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.     

A multiphase lattice Boltzmann study of buoyancy-induced mixing in a tilted channel

We study the buoyancy-induced interpenetration of two immiscible fluids in a tilted channel by a two-phase lattice Boltzmann method using a non-ideal gas equation of state well-suited for two incompressible fluids. The method is simple, elegant and easily parallelizable. After first validating the code for simulating Rayleigh–Taylor instabilities in a unstably-stratified flow, we applied the code to simulate the buoyancy-induced mixing in a tilted channel at various Atwood numbers, Reynolds numbers, tilt angles, and surface tension parameters. The effects of these parameters are studied in terms of the flow structures, front velocities, and velocity profiles. For one set of parameters, comparisons have also been made with results of a finite volume method. The present results are seen to agree well with those of a finite volume method in the interior of the flow; however near the boundary there is some discrepancy.