Improving the Stability of the Multiple-Relaxation-Time Lattice Boltzmann Method by a Viscosity Counteracting Approach

Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method (LBM). The multiple-relaxation-time (MRT) model of the LBM can improve the accuracy and stability, but is still subject to numerical instability when simulating flows with large single-grid Reynolds number (Reynolds number/grid number). The viscosity counteracting approach proposed recently is a method of enhancing the stability of the LBM. However, its effectiveness was only verified in the single-relaxation-time model of the LBM (SRT-LBM). This paper aims to propose the viscosity counteracting approach for the multiple-relaxation-time model (MRT-LBM) and analyze its numerical characteristics. The verification is conducted by simulating some benchmark cases: the two-dimensional (2D) lid-driven cavity flow, Poiseuille flow, Taylor-Green vortex flow and Couette flow, and three-dimensional (3D) rectangular jet. Qualitative and Quantitative comparisons show that the viscosity counteracting approach for the MRT-LBM has better accuracy and stability than that for the SRT-LBM.     

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

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

Extended Lattice Boltzmann Model

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.     

基于格子Boltzmann大密度比模型的多孔介质内气泡动力学行为数值模拟

基于格子Boltzmann两相流大密度模型,研究气泡穿过多孔介质的动力学行为.研究发现:当孔隙率较大时,气泡只变形不破裂,能完整地通过多孔介质;而孔隙率较小时,气泡变形更加剧烈且发生破裂,穿过多孔介质所需的时间更长.另外,当障碍物表面接触角(θ)较小时,气泡均能完整地通过多孔介质,随着接触角的增大,气泡开始发生破裂,且...     

Lattice Boltzmann Simulation for Complex Flow in a Solar Wall

In this letter,we present a lattice Boltzmann simulation for complex flow in a solar wall system which includes porous media flow and heat transfer,specifically for solar energy utilization through an unglazed transpired solar air collector(UTC).Besides the lattice Boltzmann equation(LBE) for time evolution of particle distribution function for fluid field,we introduce an analogy,LBE for time evolution of distribution function for temperature.Both temperature fields of fluid(air) and solid(porous media) are modeled.We study the effects of fan velocity,solar radiation intensity,porosity,etc.on the thermal performance of the UTC.In general,our simulation results are in good agreement with what in literature.With the current system setting,both fan velocity and solar radiation intensity have significant effect on the thermal performance of the UTC.However,it is shown that the porosity has negligible effect on the heat collector indicating the current system setting might not be realistic.Further examinations of thermal performance in different UTC systems are ongoing.The results are expected to present in near future.     

近壁面气泡运动特性的数值计算

采用体积加速度模型确定水下爆炸气泡运动的初始条件,基于MSC.DYTRAN有限元软件开发了定义流场初始条件和边界条件的子程序,对气泡在水平刚性壁面附近的运动特性进行了数值模拟,通过对比发现数值计算结果与实验结果具有较好的一致性,证明了初始条件定义、子程序开发和有限元模型建立的正确性和数值计算的准确性。以此为基础,研究了水深、泡心与刚性底面之间的距离对气泡动态特性和射流速度的影响,通过数值计算得到了一些有规律的曲线。计算模型、方法及结果对相关的工程研究和计算具有一定的参考价值。     

三维双气泡融合的格子Boltzmann模拟

针对基于自由能模型的格子Boltzmann方法,推导D3Q15格子模型对应的平衡态分布函数;采用该模型模拟双气泡的融合过程.结果表明,气泡的融合不仅与它们之间的初始间距有关,还与表面张力有关.表面张力越大,气泡融合的临界距离也越大.此外,研究气泡融合速度与初始间距、表面张力及粘性系数的关系.     

过热液体中蒸汽泡上升过程的格子Boltzmann三维数值模拟

应用格子Boltzmann相变模型,在三维空间研究蒸汽泡在过热液体中生长、上升和变形等动力学行为.为研究传热传质对蒸汽泡运动的影响,对比模拟相同条件下气泡在等温环境中上升的物理过程.结果表明：蒸汽泡在过热液体中上升发生的变形程度较小,意味着相变对蒸汽泡的影响和表面张力一样使汽泡保持初始的形状.蒸汽泡在过热液体中的上升速度较小,说明随着汽泡生长拖拽力的影响比浮力大.蒸汽泡生长率在初始阶段达到最大值,随后会趋于一个恒定的值.随着汽泡体积增大和上升速度的增加,其对流场的扰动也越来越剧烈.蒸汽泡生长和上升引起的对流运动对温度场的演化造成很大的影响.     

A Lattice Boltzmann Model for Two-Dimensional Magnetohydrodynamics

A lattice Boltzmann model （LBM） has been developed for simulating magnetohydrodynamics （MHD） along the line of Dellar [J. Comput. Phys. 179 （2002）95]. In this model the magnetic field is presented by a vector valued magnetic distribution function which obeys a vector Boltzmann equation. The truncated error of the equilibrium distribution in the present model is up to order O（u^4） in velocity u rather than the usual 0（u^3）. For verification, the model is applied to solve the shock tube problem and the main features of the flow predicted by the model are found to compare well with the corresponding results obtained with high-order semi-discrete schemes [J. Comput. Phys. 201 （2004） 261]. The numerical experiments have also shown that the present LBM model with the equilibrium distribution truncated at O（u^4） performs much better in terms of numerical stability than those truncated at O（u^3）.     

General Propagation Lattice Boltzmann Model for the Boussinesq Equation

A general propagation lattice Boltzmann model is used to solve Boussinesq equations. Different local equilibrium distribution functions are selected, and the macroscopic equation is recovered with second order accuracy by means of the Chapman–Enskog multi-scale analysis and the Taylor expansion technique. To verify the effectiveness of the present model, some Boussinesq equations with initial boundary value problems are simulated. It is shown that our model can remain stable and accurate, which is an effective algorithm worthy of promotion and application.     

Diffusion Lattice Boltzmann Scheme on a Orthorhombic Lattice

We present a diffusion lattice Boltzmann (DLB) scheme which is derived from first principles. As opposed to the traditional lattice BGK schemes the DLB is valid for orthorhombic lattices and it has two eigenvalues of the collision operator. It is shown that the diffusion coefficient depends only on one eigenvalue of the collision operator. Hence, the DLB scheme can be optimized with means of the additional eigenvalue of the collision operator and with different lattice spacing along the principal axes. The properties of the DLB scheme concerning consistency, stability, and accuracy are studied with eigenmode analysis. This analysis shows that the DLB scheme is consistent with diffusion for a wide range of diffusion coefficients, it has unconditional stability, and that it has third-order accuracy. Furthermore, it is shown that accuracy is improved by setting the additional eigenvalue to zero and by densifying the lattice spacing along the direction of the density gradient.     

Energy-Conserving Lattice Boltzmann Thermal Model in Two Dimensions

A discrete velocity model is presented for lattice Boltzmann thermal fluid dynamics. This model is implemented and tested in two dimensions with a finite difference scheme. Comparison with analytical solutions shows an excellent agreement even for wide temperature differences. An alternative approximate approach is then presented for traditional lattice transport schemes     

Lattice Boltzmann Model for the Incompressible Navier–Stokes Equation

In this paper a lattice Boltzmann (LB) model to simulate incompressible flow is developed. The main idea is to explicitly eliminate the terms of o(M ²), where M is the Mach number, due to the density fluctuation in the existing LB models. In the proposed incompressible LB model, the pressure p instead of the mass density ρ is the independent dynamic variable. The incompressible Navier–Stokes equations are derived from the incompressible LB model via Chapman–Enskog procedure. Numerical results of simulations of the plane Poiseuille flow driven either by pressure gradient or a fixed velocity profile at entrance as well as of the 2D Womersley flow are presented. The numerical results are found to be in excellent agreement with theory.     

考虑双向耦合效应的格子Boltzmann气固两相湍流模型

在流体粒子概率密度函数输运方程中考虑颗粒对流体的反作用力,发展了考虑双向耦合效应的LB气固两相流模型,引入Smagorinsky亚格子模型模拟高雷诺数气相流场.对经典后台阶气固两相流动进行模拟,气相和颗粒相速度分布与实验结果进行比较,发现考虑双向耦合效应的LB气固两相流模型结果明显优于单向耦合结果.进一步研究不同惯性颗粒在流场中的弥散特性,小颗粒(St~O(0.1))对流体的跟随性较好,在流场中分布较为均匀;而St~O(1)的颗粒难被流场涡卷吸进入涡内,呈现倾向性弥散现象;大颗粒(St~O(10))由于自身惯性进入流场涡,在流场中分布较为均匀.     

Coolant Wetting Simulation on Simplified Stator Coil Model by the Phase-Field Lattice Boltzmann Method

Stator coils of automobiles in operation generate heat and are cooled by coolant poured from above. The flow characteristic of the coolant depends on the coil structure, flow condition, solid–fluid interaction, and fluid property, which has not been clarified due to its complexities. Since straight coils are aligned and layered with an angle at the coolant-touchdown region, the coil structure is simplified to a horizontal square rod array referring to an actual coil size. To obtain the flow and wetting characteristics, two-phase fluid flow simulations are conducted by using the phase-field lattice Boltzmann method. First, the flow onto the single-layered rod array is discussed. The wetting area is affected both by the rod gap and the wettability, which is normalized by the gap and the averaged boundary layer thickness. Then, the flow onto the multi-layered rod arrays is investigated with different rod gaps. The top layer wetting becomes longitudinal due to the reduction of the flow advection by the second layer. The wetting area jumps up at the second layer and increases proportionally to the below layers. These become remarkable at the narrow rod gap case, and finally, the dimensionless wetting area is discussed at each layer.     

激励介质中非线性化学波的格子Boltzmann模型

构造了用于非线性化学波的格子Boltzmann模型.通过设置无对流速度场,使用多重尺度和Chapman Enskog展开,得到了平衡态分布函数的各向同性解.算例考虑了用划痕起搏,在ε²尺度上给出了Selkov系统的模拟结果,再现了远离热力学平衡态的螺旋波结构的经典结果,并与传统数值方法及实验结果进行了比较.     

Deformation and Motion of a Red Blood Cell in a Shear Flow Simulated by a Lattice Boltzmann Model

A lattice Boltzmann model is presented to simulate the deformation and motions of a red blood cell (RBC) in a shear flow. The curvatures of the membrane of a static RBC with different chemical potential drops calculated by our model agree with those computed by a shooting method very well. Our simulation results show that in a shear flow, a biconcave RBC becomes highly flattened and undergoes tank-treading motion. With intrinsically parallel dynamics, this lattice Boltzmann method is expected to find wide applications to both single and multi-vesicles suspension as well as complex open membranes in various fluid flows for a wide range of Reynolds numbers.     

液滴撞击不同湿润性固壁面上静止液滴的格子Boltzmann研究

采用单组分多相的伪势格子Boltzmann方法,在大小液滴粒径比为1.5的情况下,对大液滴竖直撞击壁面上静止小液滴的过程进行模拟,研究亲水与超疏水壁面上大液滴竖直碰撞小液滴的过程,得到液滴铺展因子和相对高度随时间的变化.结果表明:增大We数会使液滴的铺展因子增大,铺展直径变大,相对高度减小;并且随着We数的增加,在超疏...     

Simulating Solid-Liquid Phase-Change Heat Transfer in Metal Foams via a Cascaded Lattice Boltzmann Model

A cascaded lattice Boltzmann (CLB) model is constructed for simulating heat transfer in metal-foam-based solid-liquid phase change materials (PCMs). The present model captures the phase interface implicitly via the enthalpy methodology, and to avoid iterations in simulations, the CLB equation of the PCM employs the enthalpy as the basic evolution variable through modifying the cascaded collision process. Numerical results demonstrate the effectiveness and practicability of the CLB model for investigating heat transfer in solid-liquid PCMs with metal foams. The effects of the inertial coefficient, permeability and porosity on the melting process are investigated. The results indicate that the empirical correlations of inertial coefficient and permeability based on packed beds overestimate the melting rate at high porosities. Moreover, the porosity has significant impact on phase-change processes. The melting rate increases as the porosity of the metal foam decreases since heat conduction through high thermal conductive metal foam dominates the total heat transfer.     

A Novel Lattice Boltzmann Model For Reactive Flows with Fast Chemistry

A novel lattice Boltzmann model, in which we take the ratio of temperature difference in the temperature field to the environment one to be more than one order of magnitude than before, is developed to simulate two-dimensional reactive flows with fast chemistry. Different from the hybrid scheme for reactive flows [Comput. Phys. Commun. 129 （2000）267], this scheme is strictly in a pure lattice Doltzmann style （i.e., we solve the flow, temperature, and concentration fields using the lattice Boltzmann method only）. Different from the recent non-coupled lattice Boltzmann scheme lint. J. Mod. Phys. B 17（2003）197], the fluid density in our model is coupled directly with the temperature. Excellent agreement between the present results and other numerical data shows that this scheme is an efficient numerical method for practical reactive flows with fast chemistry.