Lattice Boltzmann Solid Particles in a Lattice Boltzmann Fluid

We define a lattice Boltzmann model of solid, deformable suspensions immersed in a fluid itself described in terms of the lattice Boltzmann method. We discuss the rules governing the internal dynamics of the solid object as well as the rules specifying the interaction between solid and fluid particle. We perform a numerical drag experiment to validate the model. Finally we consider the case of a population of flexible chains in suspension in a shear stress flow and study the influence on the velocity profile.     

Analysis of Lattice Boltzmann Initialization Routines

LB simulations can be affected by the arising of initial layers due to an inconsistent initialization of the discrete LB populations. We present some previously proposed initialization routines built to overcome that problem; using the asymptotic expansion technique, we show how their features can be analyzed and, in some cases, how accuracy and computational efficiency can be improved     

Towards a Renormalized Lattice Boltzmann Equation for Fluid Turbulence

A coarse-grained Lattice Boltzmann equation is examined in which the effects of unresolved (subgrid) scales are formally incorporated within a renormalized relaxation time of the collision operator. Actual values of the renormalized relaxation time are analyzed for the practical case of high-Reynolds flows past slant bodies (airfoils).     

Lattice Boltzmann simulations of flow past a cylindrical obstacle

We present lattice Boltzmann simulations of flow past a cylindrical obstacle. Our study is based on the Lévy walk model of turbulence in a lattice Boltzmann model. We discuss pressure around the cylinder with laminar and turbulent incident flows, as well as the dependence of the von Karman street on the analog of integral scale in our model.     

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.     

辐射传输的格子Boltzmann模拟

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

微通道内气体流动的格子Boltzmann方法研究

采用格子Boltzmann方法模拟了微通道在滑移区内不同Knudsen数下的微气体Poiseuille流,分析了微气体流动的速度分布以及流量与压降的关系,并给出了相对滑移长度和Poiseuille数随Knudsen数的变化特性。研究结果表明,微气体Poiseuille流的速度轮廓呈抛物线分布,但是边界速度大于0,出现速...     

非等温流向振荡单圆柱绕流的格子Boltzmann模拟

本文采用格子Boltzmann方法模拟对非等温横向振荡圆柱绕流(Re=100)进行模拟,详细考察振荡频率(F=f_o/f_s=0.2～8.0)、振幅(A/D=0.1～0.4,A为圆柱流向最大位移,D为圆柱直径)对绕流换热的影响。结果表明:低频率、低振幅流向振荡运动将使圆柱换热增强,高频率、高振幅振荡使圆柱换热效果减弱;圆柱表面平均N_u脉动均方根值随着振荡频率、振幅的增大而增大;在锁定区域内,平均N_u及其脉动均方根值均出现极大值。     

Lattice Boltzmann Study of a Vortex Ring Impacting Spheroidal Particles

Interaction of vortex rings with solid is an important research topic of hydrodynamic. In this study, a multiple-relaxation time (MRT) lattice Boltzmann method (LBM) is used to investigate the flow of a vortex ring impacting spheroidal particles. The MRT-LBM is validated through the cases of vortex ring impacting a flat wall. The vortex evolution due to particle size, the aspect ratio of a prolate particle, as well as Reynolds $(Re)$ number are discussed in detail. When the vortex ring impacting a stationary sphere, the primary and secondary vortex rings wrap around each other, which is different from the situation of the vortex ring impacting a plate. For the vortex ring impacting with a prolate spheroid, the secondary vortex ring stretches mainly along the long axis of the ellipsoid particle. However, it is found that after the vortex wrapping stage, the primary vortex recovers along the short axis of the particle faster than that in the long axis, i.e., the primary vortex ring stretches mainly along the short axis of the particle. That has never been addressed in the literature.     

Lattice Boltzmann computational fluid dynamics in three dimensions

The recent development of the lattice gas method and its extension to the lattice Boltzmann method have provided new computational schemes for fluid dynamics. Both methods are fully paralleled and can easily model many different physical problems, including flows with complicated boundary conditions. In this paper, basic principles of a lattice Boltzmann computational method are described and applied to several three-dimensional benchmark problems. In most previous lattice gas and lattice Boltzmann methods, a face-centered-hyper-cubic lattice in four-dimensional space was used to obtain an isotropic stress tensor. To conserve computer memory, we develop a model which requires 14 moving directions instead of the usual 24 directions. Lattice Boltzmann models, describing two-phase fluid flows and magnetohydrodynamics, can be developed based on this simpler 14-directional lattice. Comparisons between three-dimensional spectral code results and results using our method are given for simple periodic geometries. An important property of the lattice Boltzmann method is that simulations for flow in simple and complex geometries have the same speed and efficiency, while all other methods, including the spectral method, are unable to model complicated geometries efficiently.     

Lattice Boltzmann Simulation of Sedimentation of a Single Elastic Dumbbell in a Newtonian Fluid

Based on the lattice Boltzmann method (LBM), the sedimentations of a single elastic dumbbell in a Newtonian fluid under different initial positions and orientations, and also that of the elastic dumbbells with different free lengths of the spring under the same initial conditions have been simulated. All of the numerical results show that the final orientations of the elastic dumbbells are in the same horizontal direction, and the final positions of their centroids are all on the centerline of the tube no matter what the initial positions and orientations of the elastic dumbbell or the free lengths of the spring are. When the elastic dumbbell finally falls down vertically, the two circular cylinders of the elastic dumbbell rotate around their own symmetry-axis respectively, and their angular velocities are equal but opposite to each other. For the sedimentations of the elastic dumbbells with different free lengths of the spring, the shorter of the free length is, the faster the final angular velocity and vertical velocity of the circular cylinder will be.     

A Lattice Boltzmann Model for Anisotropic Crystal Growth from Melt

We coupled the lattice Boltzmann method with enhanced collisions for hydrodynamics with a model for the anisotropic liquid/solid phase transition. The model is based on a simple reaction model. As a test we have performed calculations for dendritic growth of a crystal into an undercooled melt.     

Lattice Boltzmann method with the cell-population equilibrium

The central problem of the lattice Boltzmann method （LBM） is to construct a discrete equilibrium. In this paper, a multi-speed 1D cell-model of Boltzmann equation is proposed, in which the cell-population equilibrium, a direct non- negative approximation to the continuous Maxwellian distribution, plays an important part. By applying the explicit one-order Chapman-Enskog distribution, the model reduces the transportation and collision, two basic evolution steps in LBM, to the transportation of the non-equilibrium distribution. Furthermore, 1D dam-break problem is performed and the numerical results agree well with the analytic solutions.     

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.     

From Reactive Boltzmann Equations to Reaction–Diffusion Systems

We consider the reactive Boltzmann equations for a mixture of different species of molecules, including a fixed background. We propose a scaling in which the collisions involving this background are predominant, while the inelastic (reactive) binary collisions are very rare. We show that, at the formal level, the solutions of the Boltzmann equations converge toward the solutions of a reaction-diffusion system. The coefficients of this system can be expressed in terms of the cross sections of the Boltzmann kernels. We discuss various possible physical settings (gases having internal energy, presence of a boundary, etc.), and present one rigorous mathematical proof in a simplified situation (for which the existence of strong solutions to the Boltzmann equation is known).     

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.     

Asymptotic Analysis of Lattice Boltzmann Boundary Conditions

In this article, we use a general method for the analysis of finite difference schemes to investigate lattice Boltzmann algorithms for Navier–Stokes problems with Dirichlet boundary conditions. Several link based boundary conditions for commonly used lattice Boltzmann BGK models are considered. With our method, the accuracy of the algorithms can be exactly predicted. Moreover, the analytical results can be used to construct new algorithms which is demonstrated with a corrected bounce back rule that requires only local evaluations but still yields second order accuracy for the velocity. The analysis is applicable to general geometries and instationary flows     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

格子Boltzmann方法对Rayleigh-Benard流的模拟与非线性分析

本文利用格子Boltzmann方法对不同Ra下的Rayleigh-Benard自然对流进行了数值模拟,得到的Nu与前人的结果一致。数值结果表明,Ra=2000时,格子Boltzmann方法得到了三个静态分岔解,其中一个为非稳定解;在Ra=5×10~5时格子Boltzmann方法得到了类似于周期的振荡解。本文采用插值格子Boltzmann方法,在Ra=2000,取不同插值节点进行计算时,得到了一种非稳态解。