热对流条件下双颗粒沉降的直接数值模拟

在任意拉格朗日-欧拉(ALE)算法模拟等温惰性颗粒两相流的基础上,增加对能量方程的联立求解,对热对流条件下双颗粒的沉降进行了直接数值模拟.结果表明:双颗粒在等温流体、热流体和冷流体中运动形式是不同的,造成不同的主要因素是热对流引起了颗粒不同的动态尾迹,在热流体中颗粒尾部形成了涡的脱落,在冷流体中颗粒尾部形成了羽流.     

溶解与热对流对固体颗粒运动影响的直接数值模拟

对牛顿流体内溶解与热对流对单颗粒在垂直管道中的沉降运动进行了直接数值模拟.流体运动由守恒方程计算,密度和黏性的变化考虑流场温度变化的影响,通过积分黏性应力和压力获得颗粒的受力跟踪颗粒运动,溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定.通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合.分别模拟了颗粒在等温流体、热流体、冷流体及颗粒溶解四种情况下的沉降过程.结果表明,在一定雷诺数内,热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒表面形态的变化引起了颗粒的横向摆动,并使颗粒沉降速度发生了变化. 关键词： 溶解 热对流 颗粒两相流 直接数值模拟     

溶解椭圆颗粒沉降的介观尺度数值模拟

在任意拉格朗日欧拉(ALE)算法模拟有热对流影响的颗粒两相流动的直接数值模拟基础上, 通过建立颗粒溶解速度和颗粒表面热流密度的关系, 对溶解的椭圆颗粒在垂直管道内牛顿流体中的沉降进行了直接数值模拟. 计算结果表明:与等温惰性椭圆颗粒沉降相比, 流体的对流运动、颗粒质量以及形状的变化等因素使溶解的椭圆颗粒在不同初始角度沉降时, 颗粒沉降动态尾迹、颗粒受力、颗粒沉降速度等都发了较大变化.     

热对流条件下固液密度比对颗粒沉降运动影响的直接数值模拟

本文应用任意拉格朗日-欧拉(ALE)算法对固液两相流流场中考虑热对流的非等温颗粒在竖直通道中的沉降运动进行了数值模拟.在牛顿流体中通过积分黏性应力和压力获得颗粒的受力跟踪颗粒运动,使用有限元方法数值求解流场的N-S方程和能量方程,模型不需经验假设.通过模拟来研究颗粒沉降的运动规律和热对流下固液密度比对固液两相流的影响作用.结果表明随着固液密度比的增加,颗粒经历了稳定沉降、周期性摆动,不规则摆动等过程;热对流使颗粒的摆动幅度和沉降速度发生变化;热对流对颗粒的影响作用随着固液密度比的增加而减小.     

热动力学格子Boltzmann模型

提出了一类热动力学格子Boltzmann模型,其分布函数中含有内能项与高阶速度项。以三个HPP模型为基础,建立了一个三迭加HPP格子Boltzmann模型,并给出其局部平衡分布函数。设计了热动力学现象模拟中的外力处理方法和温度边界条件处理方法,用该模型成功地模拟了Benard热对流现象。     

用Lattice Boltzmann方法模拟1/4圆腔内的定常层流运动

分析格子Boltzmann方法中二阶精度的曲线边界处理方法.应用格子Boltzmann方法及其边界处理方法模拟1/4圆腔内的定常层流运动,引入流线图和等涡线图分析流场随Re数的变化.并且发现当Re数在10~100区间内变化时,随着Re数的增大,顺时针旋转流场的涡心位置偏离x轴的角度逐渐减小,而逆时针旋转流场的涡心位置偏离x轴的角度却越来越大.     

矩形腔内纳米流体自然对流的两相格子Boltzmann模拟

在考虑纳米颗粒和基液之间相互作用力的基础上,建立了两相格子Boltzmann模型,通过对比,结果表明本文建立的两相格子Boltzmann模型与相关文献吻合较好.本文利用该模型研究了矩形腔内纳米流体的流动和换热,讨论了瑞利数、纳米颗粒组分对换热的影响,并分析了纳米颗粒在流场中的分布.     

滴形颗粒在通道中沉降的格子Boltzmann-虚拟区域方法模拟

采用格子Boltzmann-虚拟区域方法对滴形颗粒在垂直通道中的沉降过程进行直接数值模拟.通过格子Boltzmann方法求解N-S方程,流体与固体之间的相互作用通过虚拟区域方法描述.研究雷诺数在10^-2到100范围内颗粒形状因子对其摩擦系数和阻力系数的影响.为便于比较,给出了圆形颗粒沉降的结果.结果发现,当雷诺数小于1的时候,颗粒的摩擦系数始终保持常数,而当雷诺数大于1时摩擦系数随雷诺数的增大而增大.此外,当雷诺数小于约30时圆形颗粒的摩擦系数和阻力系数均小于滴形颗粒,而当雷诺数大于30时情况正好相反.颗粒周围的压力分布证明了这一结论.     

微尺度稀薄效应下颗粒沉降的格子Boltzmann方法模拟

基于多松弛格子Boltzmann模型,对竖直细长微通道内颗粒自由沉降过程进行模拟,分析气体稀薄效应、初始位置以及颗粒间相互作用对微颗粒沉降特性的影响.研究表明：随Knudsen数增大,微通道内气体稀薄效应增强,颗粒表面气体滑移速度增大,气相流体有效粘度减小,颗粒相同运动状态下受到气体阻力相应减小,颗粒沉降平衡速度明显增大;不同初始位置颗粒沉降过程存在明显差异,初始位置偏离中心线颗粒将发生水平方向位移且呈振荡趋势,最终稳定于中心线平衡位置;在微尺度双颗粒沉降DKT现象过程中,气体稀薄效应影响颗粒运动特性,后颗粒跟随过程明显增长.     

双颗粒在溶解条件下沉降的多相流动特性

对牛顿流体内溶解双颗粒在垂直管道中的沉降运动进行了直接数值模拟. 流体运动由守恒方程计算, 密度和黏性的变化考虑流场温度变化的影响, 通过积分粘性应力和压力获得颗粒的受力跟踪颗粒运动, 溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定. 通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合. 对双颗粒在等温流体无溶解条件和非等温流体溶解条件下的沉降过程进行了计算. 结果表明, 在一定雷诺数内, 热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒质量、颗粒表面形态的变化引起了颗粒的横向摆动, 并使颗粒沉降速度发生了变化.     

Lattice Boltzmann Simulations of Particle-Particle Interaction in Steady Poiseuille Flow

The moving behaviour of two- and three-particles in a pressure-driven flow is studied by the lattice Boltzmann simulation in two dimensions. The time-dependent values, including particles＇ radial positions, translational velocities, angular velocities, and the x-directional distance between the particles are analysed extensively. The effect of flow Reynolds number on particle motion is also investigated numerically. The simulation results show that the leading particle equilibrium position is closer to the channel centre while the trailing particle equilibrium position is closer to the channel wall. If Reynolds number Re is less than 85.30, the larger flow Reynolds number results in the smaller x-directional equilibrium distance, otherwise the x-directional distance increases almost linearly with the increase of time and the particles separate finally. The simulation results are helpful to understand the particle-particle interaction in suspensions with swarms of particles.     

Provide a suitable range to include the thermal creeping effect on slip velocity and temperature jump of an air flow in a nanochannel by lattice Boltzmann method

The thermal creeping effect on slip velocity of air forced convection through a nanochannel is studied for the first time by using a lattice Boltzmann method. The nanochannel side walls are kept hot while the cold inlet air streams along them. The computations are presented for the wide range of Reynolds number, Knudsen number and Eckert number while slip velocity and temperature jump effects are involved. Moreover appropriate validations are performed versus previous works concerned the micro–nanoflows.The achieved results are shown as the velocity and temperature profiles at different cross sections, streamlines and isotherms and also the values of slip velocity and temperature jump along the nanochannel walls. The ability of the lattice Boltzmann method to simulate the thermal creeping effects on hydrodynamic and thermal domains of flow is shown at this study; so that its effects should be involved at lower values of Eckert number and higher values of Reynolds number especially at entrance region where the most temperature gradient exists.     

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 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.     

Numerical simulation of particle motion in an ultrasound field using the lattice Boltzmann model

In this paper we investigate the motion of small particles suspended in a fluid through which an ultrasound field is propagating. The application of the lattice Boltzmann model to this problem is considered using a two dimensional model. Particles in an ultrasound field are observed to move with a mean particle motion. Further, the time-averaged force on a fixed cylinder is computed and found to be in good agreement with a theoretical expression for the radiation force. Simulations are performed with a single particle, although the approach can equally be applied for a larger number of particles.     

Numerical simulation of particle motion in an ultrasound field using the lattice Boltzmann model

In this paper we investigate the motion of small particles suspended in a fluid through which an ultrasound field is propagating. The application of the lattice Boltzmann model to this problem is considered using a two dimensional model. Particles in an ultrasound field are observed to move with a mean particle motion. Further, the time-averaged force on a fixed cylinder is computed and found to be in good agreement with a theoretical expression for the radiation force. Simulations are performed with a single particle, although the approach can equally be applied for a larger number of particles.     

Influence of wall motion on particle sedimentation using hybrid LB-IBM scheme

We integrate the lattice Boltzmann method (LBM) and immersed boundary method (IBM) to capture the coupling between a rigid boundary surface and the hydrodynamic response of an enclosed particle laden fluid. We focus on a rigid box filled with a Newtonian fluid where the drag force based on the slip velocity at the wall and settling particles induces the interaction. We impose an external harmonic oscillation on the system boundary and found interesting results in the sedimentation behavior. Our results reveal that the sedimentation and particle locations are sensitive to the boundary walls oscillation amplitude and the subsequent changes on the enclosed flow field. Two different particle distribution analyses were performed and showed the presence of an agglomerate structure of particles. Despite the increase in the amplitude of wall motion, the turbulence level of the flow field and distribution of particles are found to be less in quantity compared to the stationary walls. The integrated LBM-IBM methodology promised the prospect of an efficient and accurate dynamic coupling between a non-compliant bounding surface and flow field in a wide-range of systems. Understanding the dynamics of the fluid-filled box can be particularly important in a simulation of particle deposition within biological systems and other engineering applications.     

基于总能形式的耦合的双分布函数热晶格玻尔兹曼数值方法

双分布函数热晶格玻尔兹曼数值方法在微尺度热流动系统中得到广泛的应用. 本文基于晶格玻尔兹曼平衡分布函数低阶Hermite展开式, 创新性地提出了包含黏性热耗散和压缩功的耦合的双分布函数热晶格玻尔兹曼数值方法, 将能量场内温度的变化以动量源的形式引入晶格波尔兹曼动量演化方程, 实现了能量场与动量场之间的耦合. 研究了考虑黏性热耗散和压缩功的和不考虑的两种热自然对流模型, 重点分析了不同瑞利数和普朗特数下流场内的流动情况以及温度、速度和平均努赛尔数的变化趋势. 本文实验结果与文献结果一致, 验证了本文数值方法的可行性和准确性. 研究结果表明: 随着瑞利数和普朗特数的增大, 方腔内对流传热作用逐渐增强, 边界处形成明显的边界层; 考虑黏性热耗散和压缩功的模型对流作用相对增强, 黏性热耗散和压缩功对自然对流的影响在微尺度流动过程中不能忽略.