基于格子boltzmann方法模拟热毛细对流

《基于格子boltzmann方法模拟热毛细对流》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、基于格子Boltzmann方法模拟热毛细对流重庆大学博士学位论文学生姓名：谢海琼指导教师：曾忠业：流体力学学科门类：工学教授专重庆大学航空航天学院二〇一五年九月SimulationoftheThermocapillaryConvectionwithLatticeBoltzmannMethodAThesisSubmittedtoChongqingUniversityinPartialFulfillmentoftheRequirementfortheDoctor’sDegreeofEngineeringByHaiqiongXieSupervisedbyProf.Ze

2、ngZhongSpecialty:FluidMechanicsCollegeofAerospaceEngineeringofChongqingUniversity,Chongqing,ChinaSeptember2015中文摘要摘要气液界面（自由液面）或者液液界面的流体力学问题广泛存在于自然界与工业生产实践中。流体的界面是一个相到另一个相之间很薄的过渡区域，分子间的相互作用力在这个过渡区域宏观上表现为表面张力，其大小一般是温度和浓度的函数。其中，流体界面上的非均匀温度引起的切向表面张力梯度会驱动流体流动，这种流动称为热毛细对流。在微重力或微尺度条件下，重力和浮力

3、对流体流动的影响极度地减弱，热毛细流可能成为主要的对流。在热毛细对流问题中，温度梯度引起的非均匀表面张力即驱动流体流动又影响界面形状，这两者的动态耦合使热毛细对流的数值模拟成为计算流体动力学长期以来的一个重要挑战。基于格子Boltzmann方法（LatticeBoltzmannMethod:LBM），本文发展了考虑界面动态变形的热毛细对流问题的计算方法，其中两种不相溶的流体流动采用多松弛不可压缩LBM方法模拟，界面的动态变形则采用相场法进行捕捉。通过自主开发的程序包，本文数值模拟研究了双液滴热毛细迁移与合并，以及双层环形液池内热毛细对流和界面变形耦合过程。本文的

4、主要工作与贡献如下：(1)首先简要地介绍了LBM方法，推导了格子Boltzmann方程恢复到宏观Navier-Stokes方程的过程，并建立了不可压缩多松弛LBM模型。随后，对经典流动问题进行数值模拟，验证了多松弛不可压缩LBM基础模型及其程序的正确性。(2)基于上述LBM基础模型，建立了两种不同的耦合模型用于模拟考虑界面动态变形的等温两相流问题，两种模型中均采用LBM方法求解两相流流场，而界面运动则分别通过FrontTracking方法和相场法捕捉。最后基于OpenMP并行技术对自主开发的三维问题程序包实现并行计算。(3)在LBM耦合相场法模型中进一步引入能量

5、方程研究热毛细对流问题。建立了轴对称坐标系下的LBM模型，两相流流场采用LBM进行模拟，相场方程和能量方程则采用有限差分法求解。在模型中根据相场方程得到的指示函数，通过引入界面力模型（ContinuumSurfaceForce：CSF）实现法向和切向表面张力的计算。(4)系统地研究了双液滴的热毛细迁移及合并过程。首次通过数值模拟揭示了Marangoni数和Capillary数对两液滴合并过程的影响。研究表明在随Marangoni数增大，尾随液滴更难追上前导液滴并与之合并。此外，随Capillary数增大前导液滴的变形更明显，同样尾随液滴亦更难追上前导液滴并与之合

6、并，同时液滴之I重庆大学博士学位论文间合并的持续时间也随Capillary数增大而增加。(5)系统地研究了多种因素对双层环形液池内热毛细对流和界面变形的影响。结果表明上层流体具有较大的密度或者粘性时，下层流体的流动强度被明显削弱，但液液界面的变形程度也随之增加。不同Capillary数对流动强度的影响较小但对界面变形的影响却非常明显。关键词：格子Boltzmann方法，两相流，热毛细对流，相场法，表面张力II英文摘要ABSTRACTThefluiddynamicsprobleminvolvinggas-liquidorliquid-liquidinterface

7、iswidelyobservedinnatureandindustrialproduction.Thefluidinterfaceisathintransitionregion,theinteractiveforcebetweenmoleculesbehavesmacroscopicallyasthesurfacetension.Ingeneral,thesurfacetensiondependsontheinterfacetemperatureorconcentration.Thetemperaturegradientontheinterfaceresults

8、inagradiento