格子Boltzmann方法中的曲边界处理

研究了格子Boltzmann方法中实现曲边界条件的3种格式,对它们的精度和稳定性进行了分析和比较.通过二维Poiseuille流和等边三角域上空腔流的模拟,讨论了这3种格式的数值精度和稳定性.     

Boundary value problems for the steady Boltzmann equation

We discuss steady boundary value problems for the Boltzmann equation with inflow and diffusive boundary conditions in one, two, and three dimensions, with suitable truncations of the collision kernel. General existence and uniqueness results are obtained if the domain is sufficiently small. In one dimension, the existence of solutions on general intervals is obtained by abstract fixed-point theory.     

Boundary conditions for lattice Boltzmann simulations

A heuristic interpretation of no-slip boundary conditions for lattice Boltzmann and lattice gas simulations is developed. An improvement is suggested which consists of including the wall nodes in the collision operation.     

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     

Boundary Singularity of Moments for the Linearized Boltzmann Equation

We study the boundary singularity for stationary solutions of the linearized Boltzmann equation with hard-sphere potential. An asymptotic estimate for the gradient of the moments is established, which shows the logarithmic singularity near the boundary. Our formula holds for the solutions of the Milne and Kramers problems obtained by Bardos-Caflish-Nicolaenko (Commun. Pure Appl. Math. 49:323–452, 1986). Our theorem requires the Hölder continuity of the boundary data. In particular, it applies to the complete condensation problem for half space.     

A Characteristic Non-Reflecting Boundary Treatment in Lattice Boltzmann Method

In lattice Boltzmann methods, disturbances develop at the initial stages of the simulation, the decay characteristics depend mainly on boundary treatment methods; open boundary conditions such as equilibrium and bounce-back schemes potentially generate uncontrollable disturbances. Excessive disturbances originate from non-physical reflecting waves at boundaries. Characteristic boundary conditions utilizing the signs of waves at boundaries which suppress these reflecting waves, as well as their implementation in the lattice Boltzmann method, are introduced herein. The performance of our novel boundary treatment method to effectively suppress excessive disturbances is verified by three different numerical experiments.     

格子Boltzmann守恒型非平衡态的外推边界

针对非平衡态外推边界处理方法在某些条件下发生质量泄漏的问题,结合质量守恒定律对其进行修正,建立一种新的质量守恒型的边界处理格式.其基本思想是在边界处定义一个虚拟密度来修正平衡态分布函数,使得系统质量守恒.以重力作用下的密度场及磁流体模拟问题为例,利用该格式进行分析,讨论质量泄露问题的本质及各种因素的影响.通过对三种典型的流动:二维Poiseuille流、Couette流、Womersley流的模拟,验证该边界格式的计算精度和稳定性.     

Heat Conduction Beyond the Fourier Law

Technical Physics - The Fourier law correctly describes heat transport in most practical macroscopic problems. However, for heat transfer in rapid processes, heat transport on micro- and...     

Lees–Edwards Boundary Conditions for Lattice Boltzmann

Lees–Edwards boundary conditions (LEbc) for Molecular Dynamics simulations⁽¹⁾ are an extension of the well known periodic boundary conditions and allow the simulation of bulk systems in a simple shear flow. We show how the idea of LEbc can be implemented in isothermal lattice Boltzmann simulations and how LEbc can be used to overcome the problem of a maximum shear rate that is limited to less then 1/L _y (with L _y the transverse system size) in traditional lattice Boltzmann implementations of shear flow. The only previous Lattice Boltzmann implementation of LEbc⁽²⁾ requires a specific fourth order equilibrium distribution. In this paper we show how LEbc can be implemented with the usual quadratic equilibrium distributions.     

A Weak Formulation of the Boltzmann Equation Based on the Fourier Transform

In this article we present an alternative formulation of the spatially homogeneous Boltzmann equation. Rewriting the weak form of the equation with shifted test functions and using Fourier techniques, it turns out that the transformed problem contains only a three-fold integral. Explicit formulas for the transformed collision kernel are presented in the case of VHS models for hard and soft potentials. For isotropic Maxwellian molecules, a classical result by Bobylev is recovered, too.     

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

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

Boundary Slip and Surface Interaction： A Lattice Boltzmann Simulation

The factors affecting slip length in Couette geometry flows are analysed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactions. The main factors influencing the boundary slip are the strength of interactions between fluid-fluid and fluid-wall particles. Other factors, such as fluid viscosity, bulk pressure may also change the slip length. We find that boundary slip only occurs under a certain density （bulk pressure）. If the density is large enough, the slip length will tend to zero. In our simulations, a low density layer near the wall does not need to be postulated a priori but emerges naturally from the underlying non-ideal mesoscopic dynamics. It is the low density layer that induces the boundary slip. The results may be helpful to understand recent experimental observations on the slippage of micro flows.     

Nonlinear Boundary Layers of the Boltzmann Equation: I. Existence

 We study the half-space problem of the nonlinear Boltzmann equation, assigning the Dirichlet data for outgoing particles at the boundary and a Maxwellian as the far field. We will show that the solvability of the problem changes with the Mach number ℳ^∞ of the far Maxwellian. If ℳ^∞<−1, there exists a unique smooth solution connecting the Dirichlet data and the far Maxwellian for any Dirichlet data sufficiently close to the far Maxwellian. Otherwise, such a solution exists only for the Dirichlet data satisfying certain admissible conditions. The set of admissible Dirichlet data forms a smooth manifold of codimension 1 for the case −1<ℳ^∞<0, 4 for 0<ℳ^∞<1 and 5 for ℳ^∞>1, respectively. We also show that the same is true for the linearized problem at the far Maxwellian, and the manifold is, then, a hyperplane. The proof is essentially based on the macro-micro or hydrodynamics-kinetic decomposition of solutions combined with an artificial damping term and a spatially exponential decay weight. Received: 20 April 2002 / Accepted: 4 December 2002 Published online: 21 March 2003 Communicated by H.-T. Yau     

格子Boltzmann方法中三种流固耦合格式的对比

详细对比分析模拟流固两相流的格子Boltzmann方法中三类不同的流固耦合格式(半反弹方法、亚网格方法和曲边界方法)的特点,发现半反弹方法所获得的速度、角速度存在较大的扰动,而亚网格颗粒方法则能够获得平滑的发展曲线.     

Boundary conditions for Boltzmann equation and heat transfer between parallel plates

The methods developed for linear transport-relaxation equations are applied to establish boundary conditions for the linearized Boltzmann equation itself. Interfacial entropy production and reciprocity postulate play the decisive roles. Heat transfer between parallel plates as an example.     

A Half-space Problem for the Boltzmann Equation with Specular Reflection Boundary Condition

There are many open problems on the stability of nonlinear wave patterns to the Boltzmann equation even though the corresponding stability theory has been comparatively well-established for the gas dynamical systems. In this paper, we study the nonlinear stability of a rarefaction wave profile to the Boltzmann equation with the boundary effect imposed by specular reflection for both the hard sphere model and the hard potential model with angular cut-off. The analysis is based on the property of the solution and its derivatives which are either odd or even functions at the boundary coming from specular reflection, and the decomposition on both the solution and the Boltzmann equation introduced in [24, 26] for energy method.Research supported by the RGC Competitive Earmarked Research Grant, CityU 1142/01P.Research supported by the JSPS Research Fellowship for Foreign Researchers, the National Natural Science Foundation of China (10329101, 10431060), the National Key Program for Basic Research of China under grant 2002CCA03700, and the grant from the Chinese Academy of Sciences entitled Yin Jin Guo Wai Jie Chu Ren Cai Ji Jin.     

On Implementation of Boundary Conditions in the Application of Finite Volume Lattice Boltzmann Method

A new implementation of boundary condition based on the half-covolume and bounce-back rule for the non-equilibrium distribution function for the finite volume LBM is proposed here. The numerical simulation results for the expansion channel flow and driven cavity problem indicate that this method is workable for arbitrary meshes. In addition, the fourth order Runge–Kutta scheme is found to be a practical way in the LBM to accelerate the calculation speed.     

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model(DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation.It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model.Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model.To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified.     

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.     

On Stable Wall Boundary Conditions for the Hermite Discretization of the Linearised Boltzmann Equation

We define certain criteria, using the characteristic decomposition of the boundary conditions and energy estimates, which a set of stable boundary conditions for a linear initial boundary value problem, involving a symmetric hyperbolic system, must satisfy. We first use these stability criteria to show the instability of the Maxwell boundary conditions proposed by Grad (Commun Pure Appl Math 2(4):331–407, 1949). We then recognise a special block structure of the moment equations which arises due to the recursion relations and the orthogonality of the Hermite polynomials; the block structure will help us in formulating stable boundary conditions for an arbitrary order Hermite discretization of the Boltzmann equation. The formulation of stable boundary conditions relies upon an Onsager matrix which will be constructed such that the newly proposed boundary conditions stay close to the Maxwell boundary conditions at least in the lower order moments.