双扩散自然对流的格子Boltzmann模拟

建立了一个模拟双扩散自然对流系统的格子Boltzmann模型,用此模型对受温度和浓度梯度驱动的方腔流动进行了模拟,研究了浓度Rayleigh数对传热传质的影响,并与其它文献的结果进行了比较.     

辐射传输的格子Boltzmann模拟

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

蒸汽剪切驱动液滴行为的格子Boltzmann模拟

滴状冷凝过程中,存在蒸汽流动对液滴的吹扫作用,液滴在蒸汽剪切作用下克服壁面黏附变形和运动,液滴运动速度越大,冷凝传热性能越高。但是液滴在蒸汽作用下变形和运动的细节还不清晰,蒸汽速度对液滴变形和运动的影响机理还不明确。本文采用自由能格子Boltzmann方法研究了在不同蒸汽速度剪切作用下,液滴在具有不同润湿性固体表面上的变形和运动过程,分析了蒸汽速度和接触角对液滴变形和运动的影响机制,结果显示随着蒸汽速度的增加,液滴变形越大,液滴在固体表面的运动速度越大,停留时间越短,有利于液滴的移除和表面更新,相同蒸汽速度的作用下,液滴在接触角大的固体表面上变形和运动速度越大,也有利于液滴的移除和表面更新。从而定性或半定量地揭示了蒸汽速度影响蒸汽滴状冷凝传热的物理机制。     

基于格子Boltzmann方法的二维气泡群熟化过程模拟

Ostwald熟化(ripening)是指局部热力学平衡状态下,颗粒/液滴/气泡系统为了减小界面能而自发地进行颗粒群尺度分布调整的过程,具有重要研究价值.针对目前数值模拟研究不充分的现状,本文采用格子Boltzmann方法,对相变速率主控的二维蒸气泡系统演化开展了数值模拟研究.模拟结果与本文推导的二维气泡群演化标度律符合较好,证实了格子Boltzmann方法对复杂相变-物质输运过程捕捉的准确性.研究同时表明,蒸气泡系统演化过程中物质输运为液相压力不平衡所驱动,并且在小气泡“溃灭”过程中水动力学作用会影响气泡群半径分布函数的局部细节;气-液状态方程参数对熟化过程的影响效果分析显示,气液两相比内能差是驱动相变的核心要素,此差异越大相变速率越快,该结论进一步诠释了化学势驱动熟化过程的物理图像.     

二维空腔黏性流的格子Boltzmann方法模拟

用13速六方格子BhatnagarGrossKrook(缩写为BGK)模型模拟二维空腔黏性流.给出了上边界流体作匀速运动时,具有不同雷诺数的空腔黏性流的流场速度分布情况,模拟了在雷诺数Re=3000时,流场中的涡旋形成过程及流场稳定后,腔内密度、压力和温度的分布情况 关键词：     

微通道中气体流动的格子Boltzmann数值模拟

采用可压缩格子Boltzmann模型及非平衡外推边界条件,数值模拟微通道中的气体在滑移区域(Kn≤0.1)内的流动,计算结果包括出口速度剖面、通道中心压力分布以及质量流率等,与理论结果及其他实验结果符合得很好.还模拟了180°弯曲通道中的气体流动.结果表明,滑移速度的存在抑制了边界层的分离,因此在弯曲处不存在漩涡.计算结果还表明,弯道的存在显著影响了气体的质量流率.     

液滴撞击液膜过程的格子Boltzmann方法模拟

本文采用格子Boltzmann方法对液滴撞击液膜过程进行了研究, 主要考察了雷诺数(Re)、韦伯数(We)、相对液膜厚度 (h) 以及表面张力 (σ) 等物理参数对界面运动过程的影响. 首先, 随着Re数和We数的增加, 可以明显观察到液滴撞击液膜过程中形成的皇冠状水花以及卷吸现象; 当Re数较大时, 液体会发生飞溅, 由液体飞溅形成的小液滴则会继续下落, 并与液膜再次发生碰撞. 其次, 当相对液膜厚度较小时, 液滴撞击液膜并最终导致液膜断裂; 然而随着相对液膜厚度的增大, 尽管撞击过程溅起的液体会越来越多, 但是液膜并不会发生断裂. 再次, 随着表面张力的增大, 界面变形阻力增大, 撞击过程中产生的界面形变也逐渐减弱. 最后还发现皇冠(由液滴溅起形成)半径r 随时间满足r/(2R) ≈ α√Ut/(2R), 这一结果与已有结论是一致的.     

室内可吸入颗粒物运动特性的格子 Boltzmann方法模拟

为探究室内可吸入颗粒物的运动特性,使用格子Boltzmann方法,在颗粒物运动概率模型的基础上考虑布朗力对颗粒物的作用,利用改良的LB-CA（Lattice Boltzmann-cellular automata）模型,模拟了Re数分别为400、1 000和2 000时粒径为0.01 μm、0.1 μm和1 μm的颗粒物在上送上回和上送侧回两种回风形式中的运动特性.结果显示：颗粒物的空间分布范围随着Re数的增大而增大,小粒径的颗粒物受到气流湍动和扩散作用的影响更明显;颗粒物的均方位移（mean square displacement,MSD）与Re数、颗粒物粒径的大小成反比,而同样的Re数下,颗粒物在上送侧回的回风形式中其MSD较大.总体上看上送侧回的回风形式具有较低的悬浮颗粒数和更高的室内空气品质.     

底部局部加热多孔介质自然对流传热的格子Boltzmann模拟

运用格子Boltzmann方法研究了底部局部加热多孔介质方腔的自然对流传热.方腔的上壁面为低温热源,下壁面为局部高温热源,左右壁面为绝热条件.重点分析了高温热源位置a及尺寸b对多孔介质方腔自然对流传热性能的影响,提出了平均Nusselt数Nu和位置a及尺寸b的拟合关系式.研究结果表明:高温热源位置及尺寸对多孔介质方腔内自然对流传热性质的影响很大,且存在最佳高温热源位置(a=4/16)和尺寸(b=0.75),以达到最强的对流换热强度(Nu_(max)≈10.35)和最大的对流换热量(Q_(max)≈5.69).     

用格子Boltzmann方法数值模拟三维空化现象

空化是一种微观、瞬时、随机、多相的复杂现象,其过程中所产生的极端条件以及伴随的一系列空化效应,将对液流系统产生破坏性和建设性两方面的作用.采用基于Shan-Chen模型的单组分多相流格子Boltzmann方法对水体中的三维空化现象进行了数值模拟,研究了低压下水体中气核半径与空化现象的相互关系,成功再现了低压下水体中微小气核发展成气泡的过程,并进一步研究了水体依次流经低压区、高压区时空化产生、发展、溃灭的全过程.数值模拟结果和理论预测结果符合良好. 关键词： 单组分多相流 格子Boltzmann方法 三维空化     

三维格子Boltzmann方法研究上升双气泡的运动特性

应用格子Boltzmann三维模型，对双气泡在静水中的运动进行数值研究.采用八点差分和十八点差分格式分别求解一阶▽φ和二阶▽²φ可以有效避免气液密度比过大造成的数值不稳定问题.结果表明：当两个相同直径的气泡在上升时，位置靠上的气泡形状变化像单气泡上升一样，而位置靠下的气泡会受到前一个气泡尾迹的影响，并有很明显的形状变化.当两个气泡直径不同时，不管初始位置如何，大气泡总会对小气泡造成强烈的影响.     

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

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

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

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

格子Boltzmann方法模拟双液滴撞击液膜的流动过程

采用单相自由面格子Boltzmann方法中流场计算,自由面条件和表面张力处理方法,研究双液滴撞击液膜的流动过程,实现对垂直间隔的两液滴相继撞击液膜流动过程的模拟.数值结果表明,液滴相继冲击液膜时,液滴间上下距离对液膜状态有很大的影响,液滴间距不同,液膜形状会出现很大差异.数值结果与实验结论定性一致.     

一个三维多松弛时间全速域格子Boltzmann模型

构建一个既适用于低速不可压流体又适用于高速可压缩流体的三维自由参数多松弛时间格子Boltzmann模型.模型中,根据SO（3）群的不可约表述基函数构造转化矩阵,根据恢复可压Navier-Stokes方程的需要选取非守恒矩平衡值.通过von Neumann稳定性分析模型参数对数值稳定性的影响,并给出建议选择范围.模型经过基准问题的验证,模拟结果与解析解及其它数值结果符合较好.     

Lattice Boltzmann Simulation of 3D Nematic Liquid Crystal near Phase Transition

Phase transition between nematic and isotropic liquid crystal is a very weak first order phase transition.We avoid to use the normal Landau-de Gennes‘s free energy that reduces a strong first order transition,and set up a data base of free energy calculated by means of Tao-Sheng-Lin‘s extended molecular field theory that can explain the experiments of the equilibrium properties of nematic liquid crystal very well.Then we use the free energy method of lattice Boltzmann developed by Oxford group to study the phase decomposition,pattern formation in the flow of the liquid crystal near transition temperature.     

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.     

纳米SiO2材料声子热输运的格子玻尔兹曼法模拟

根据德拜物理模型,采用格子玻尔兹曼方法（LBM）对纳米二氧化硅薄膜内的声子热输运特性进行了模拟分析,得到了薄膜内的温度响应特性;在此基础上,分析了其法向有效导热系数。计算结果表明,当努森数大于0．01时,薄膜边界处出现温度跳跃,呈现出明显的微纳米尺度传热特性;当薄膜厚度小于20nm时,减小厚度可使其有效导热系数迅速降低...