激波的格子玻尔兹曼方法模拟

格子玻尔兹曼方法已被广泛应用于模拟流体,但主要应用于模拟不可压流体。应用二维正方格阵上的3速9点模型模拟激波的传播。考察了激波关系式、激波的宽度和结构。模拟结果与由Chapman-Enskog展开导出的流体方程的解相符。这说明格子玻尔兹曼方法能够模拟可压缩流体。     

应用于非线性热传导方程的格子玻尔兹曼方法

将格子玻尔兹曼方法应用于非线性热传导方程的求解,详细推导一种新的Lattice Boltzmann模型,并给出新方法所对应的多尺度方案和宏观量形式.导热系数与温度之间满足多项式函数关系,计算中模拟了不同的参数情况,并与线性热传导方程的理论解进行比较.新的Lattice Boltzmann方法展现出极大的灵活性和普适性,具有很好的应用前景.     

直通道内两相流动的格子玻尔兹曼方法模拟

本文改进了多弛豫时间多组分伪势格子玻尔兹曼方法,使其可以适用于极大黏度比的情况.研究了二维直通道内的两相流动,分析了饱和度、黏度比、润湿性和毛细管数对相对渗透率和相渗曲线的影响.结果表明,大黏度相的相对渗透率随着黏度比的增大而增大,维持高流动能力的饱和度范围也随之扩大,小黏度相与之相反;黏性耦合效应对两相相对渗透率的影响存在极限,即不会导致小黏度相的相对渗透率低于最小值,也不会使大黏度相的相对渗透率高于最大值;非润湿条件会使大黏度相的流动能力增强,相对渗透率增加,并强化相渗曲线的非单调性;润湿条件会减弱大黏度相的流动能力,相对渗透率降低;毛细管数的增加会导致更明显的接触角滞后效应,使两相相对渗透率降低.     

超短激光打孔中快速相变的格子玻尔兹曼模拟

采用双重分布函数的格子玻尔兹曼模型，对单脉冲激光金属打孔过程中的快速相变传热进行研究.模型考虑了金属材料熔化后熔体的流动换热，并采用浸没移动边界方案对过程中的固液界面进行追踪.采用纯导热模型和考虑对流的换热模型计算，将结果和试验进行对比，结果表明：在激光打孔过程中熔体的流动对相变传热产生较大影响，采用考虑流动换热模型的结果与实验更接近.进而对熔化速度、熔化深度以及温度场的变化进行分析，并探讨不同激光工艺参数对相变过程的影响.模拟发现一个脉冲结束后，激光的脉宽越大，孔深越小，孔径越大，且最高温度较短脉冲激光越低.     

凝固前沿和气泡相互作用的大密度比格子玻尔兹曼方法模拟

建立了二维双组分两相流的大密度比格子玻尔兹曼方法 (lattice Boltzmann method, LBM)模型. 该模型基于改进的Shan-Chen伪势多相流LBM模型, 结合采用不同时间步长的方法, 实现密度比达800以上的气液两相流模拟. 为了对模型进行验证, 模拟了在不同气液相互作用系数和密度比条件下气泡内外压力差与其半径之间的关系, 其结果满足Laplace定律. 将所建立的大密度比LBM与介观尺度的元胞自动机(cellular automaton, CA)和有限差分法(FDM)相耦合, 用LBM模拟气液两相流, 用CA方法模拟固相生长, 用有限差分法模拟温度场, 采用LBM-CA-FDM耦合模型对定向凝固过程中凝固前沿的气泡与液-固界面之间的相互作用进行模拟研究. 结果表明, 绝热气泡的存在影响了温度场分布, 使得凝固前沿接近气泡时, 液-固界面凸起, 在不同的固相生长速度条件下, 出现凝固前沿淹没气泡或气泡脱离凝固前沿的不同情况, 模拟结果与实验结果符合良好. 关键词： 格子玻尔兹曼方法 元胞自动机 凝固 气泡     

偏晶合金液-液相分离的格子玻尔兹曼方法模拟

本文建立了一个模拟在弥散相液滴的扩散长大、碰撞凝并和Ostwald熟化等因素的作用下偏晶合金液-液相分离过程的二维格子玻尔兹曼方法 (lattice Boltzmann method, LBM) 模型.该模型结合了Shan-Chen的两相流模型和Qin的介观粒子相互作用势模型的优点,并在LB演化方程中引入了反映相变的源项.应用该模型模拟研究了偏晶合金液-液相分离过程中单液滴的生长、两液滴的合并和多液滴的生长规律.结果表明在两液相区中第二相单个液滴的生长是一个通过扩散从非平衡态到平衡态过渡的过程.两液滴合并 关键词： 偏晶合金 液-液相分离 格子玻尔兹曼方法     

纳米结构超疏水表面冷凝现象的三维格子玻尔兹曼方法模拟

采用三维多相流格子玻尔兹曼方法 (lattice Boltzmann method, LBM),对纳米结构超疏水表面液滴的冷凝行为进行模拟研究.通过Laplace定律和光滑表面的本征接触角理论对三维LBM模型进行定量验证.模拟分析了超疏水表面纳米阵列的几何尺寸和润湿性的局部不均匀性对冷凝液滴形核位置和最终润湿状态的影响规律.结果表明,较高的纳米阵列使液滴在纳米结构间隙的上部侧面和底部优先形核长大,通过采用上下不均匀的间隙可避免液滴在底部形核长大,而在上部侧面形核的冷凝液滴在生长过程中向上运动,其润湿状态由Wenzel态转变为Cassie态;较低的纳米阵列使液滴在纳米结构底部优先形核长大,液滴的最终润湿状态为Wenzel态;润湿性不均匀的纳米结构表面使液滴在阵列顶端亲水位置处优先形核长大,成为Cassie态.冷凝液滴在不同几何尺寸的纳米结构表面上的最终润湿状态的模拟结果与文献报道的实验结果符合良好.通过模拟还发现,冷凝液滴在生长过程中的运动行为与液滴统计平均作用力的变化有关.本文的LBM模拟再现了三维空间中液滴的形核、长大和润湿状态转变等物理现象.     

负压强机制的格子玻尔兹曼方法研究

采用格子玻尔兹曼方法(LBM)研究了气液相变动态演化过程,并揭示气液混合体系中负压强的产生机理.本文采用多相流模拟中的单-陈模型(Shan-Chen模型)研究相变问题,该模型通过一个伪势来表示不同相之间的相互作用,从而控制不同相的分离,结合粒子间的相互作用力可得到一个能够描述非理想气体的状态方程,通过研究此状态方程,确定了发生相变和产生负压强的临界条件.再结合LBM,对相变和负压强现象进行数值模拟,并在气液混合体系中对拉普拉斯定律进行了验证.从模拟结果中发现,当液滴与周围气体处于力平衡和热平衡时,液滴内外压强差与其半径之间的关系满足拉普拉斯定律;另外,在气液交界面处会产生负压强,为使得理论解释与数值模拟结果相符,对于此处负压强的起源问题,我们采用同样能够描述非理想气体的范德瓦尔斯方程结合分子动力学的方法,导出内压强的变化会导致负压强的出现,并通过解释内压强的产生原因,从而进一步了解了负压强产生的微观机制.     

非对称弯曲微流道中粒子惯性聚焦动态过程及流速调控机理研究

设计制作了一种具有非对称弯曲微流道结构的微流控芯片, 搭建实验平台定量表征聚苯乙烯粒子和血细胞沿流道的动态惯性聚焦过程, 并系统研究了流体流速和粒子尺寸对粒子聚焦特性的调控机理. 通过分析粒子荧光图谱和对应量化强度曲线, 将粒子沿流道长度的横向迁移过程分为形成聚焦和平衡位置调整两个阶段, 指出在整个聚焦过程中具有小曲率半径的流道结构起主导作用. 根据全流速段内粒子聚焦特性的演变, 重点分析潜在惯性升力和Dean 曳力的竞争机制, 提出了阐述粒子聚焦流速调控过程的三阶段模型. 进一步比较两种尺寸粒子聚焦位置和聚焦率随流速与流道长度的变化规律, 发现大粒子具有更好的聚焦效果和稳定性, 且两种粒子的相对位置可通过流速进行调整. 最后, 通过分析血细胞在非对称弯流道中的横向迁移特性, 验证了粒子惯性聚焦机理在复杂生物粒子操控方面的适用性. 上述结论为深入研究微流体环境下粒子的运动特性以及开发微流式细胞术等临床即时诊断器件提供了重要参考.     

基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究

针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法(immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS)的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.     

Multi-relaxation time lattice Boltzmann simulation of inertial secondary flow in a curved microchannel

The inertial secondary flow is particularly important for hydrodynamic focusing and particle manipulation in biomedical research.In this paper,the development of the inertial secondary flow structure in a curved microchannel was investigated by the multi relaxation time lattice Boltzmann equation model with a force term.The numerical results indicate that the viscous and inertial competition dominates the development of secondary flow structure development.The Reynolds number,Dean number,and the cross section aspect ratio influence significantly on the development of the secondary vortexes.Both the intensity of secondary flow and the distance between the normalized vortex centers are functions of Dean numbers but independent of channel curvature radius.In addition,the competition mechanism between the viscous and inertial effects were discussed by performing the particle focusing experiments.The present investigation provides an improved understanding of the development of inertial secondary flows in curved microchannels.     

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.     

Lattice Boltzmann simulation of transverse wave travelling in Maxwell viscoelastic fluid

A nine-velocity lattice Boltzmann method for Maxwell viscoelastic fluid is proposed. Travelling of transverse wavein Maxwell viscoelastic fluid is simulated. The instantaneous oscillating velocity, transverse shear speed and decay rateagree with theoretical results very well.     

Lattice Boltzmann simulation of fluid flows in two-dimensional channel with complex geometries

Boundary conditions (BCs) play an essential role in lattice Boltzmann (LB) simulations. This paper investigates several most commonly applied BCs by evaluating the relative L₂-norm errors of the LB simulations for two-dimensional (2-D) Poiseuille flow. It is found that the relative L₂-norm error resulting from FHML's BC is smaller than that from other BCs as a whole. Then, based on the FHML's BC, it formulates an LB model for simulating fluid flows in 2-D channel with complex geometries. Afterwards, the flows between two inclined plates, in a pulmonary blood vessel and in a blood vessel with local expansion region, are simulated. The numerical results are in good agreement with the analytical predictions and clearly show that the model is effective. It is expected that the model can be extended to simulate some real biologic flows, such as blood flows in arteries, vessels with stenosises, aneurysms and bifurcations, etc.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

椭圆柱体在牛顿流体中运动的格子Boltzmann方法模拟

用格子Boltzmann方法建立了椭圆柱体的二维动力学模型，利用所建立的模型，数值模拟了 牛顿流体中不同形状的椭圆柱体在相同初始条件下的运动和同一椭圆柱体在不同初始条件下 的运动，并通过比较相同条件下圆柱体的运动，讨论了椭圆柱体二维运动的特征，得到了一 些有意义的结果. 关键词： 格子Boltzmann方法 椭圆柱体 牛顿流体     

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.     

Lattice Boltzmann simulations of a dumbbell moving in a Poiseuille flow

In this paper, the lattice Boltzmann method is applied to simulate a dumbbell moving in a pressure-driven flow in a planar channel with the stress-integration method for the evaluation of hydrodynamic force acting on the cylinders. The simulation results show that the dumbbell also has the important feature of the Segr\'e--Silberberg effect like a particle in a Poiseuille flow. The dumbbell trajectories, orientations, the cylinders vertical velocities and angular velocities all reach their equilibrium values separately independent of their initial positions. It is also found that the dumbbell equilibrium positions depend on the flow Reynolds number, blockage ratio and elastic coefficient. This study is expected to be helpful to understand the dynamics of polymer solutions, polymer synthesis and reaction, etc.     

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.     

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.