液滴撞击液膜形态演变过程模拟.docxVIP

中国工程热物理学会传热传质学学术会议论文编号：123180液滴撞击液膜形态演变过程模拟沈胜强，梁刚涛，郭亚丽(大连理工大学海洋能源利用与节能教育部重点实验室，辽宁大连 116024)(Tel Email:ylguo@dlut.edu.cn)摘要：采用CLSVOF方法对滴撞击平面液膜后的形态演变过程进行了二维数值模拟，描述了不同撞击速度下液滴撞击液膜后形态演变的过程，成功捕捉到了气泡夹带现象。分析了撞击速度对飞溅过程的影响和液膜厚度对撞击后形成冠状水花扩展半径的影响，结果显示，撞击速度越大，液滴与液膜撞击后出现飞溅的时间越早，扩展半径随液膜厚度的增大而减小。关键词：CLSVOF；液滴；撞击；液膜；形态演变0 前言液滴撞击固体表面液膜的现象广泛地存在于各种工业过程中,如冶金连铸喷雾冷却、燃油液滴冲击发动机汽缸壁、水平管降膜蒸发器中液滴撞击换热管表面等，液滴撞击液膜的状况是影响机械性能的重要因素，故有必要对液滴撞击液膜后的形态演变进行研究。该演变过程是一种典型的自由表面流动问题，研究该过程不仅能更好地了解自由表面流动问题的本质,而且对于指导实际生产具有非常积极的意义。对液滴撞击液膜的研究，多采用高速数码相机拍照的实验方法。郭加宏等[1]实验观测了液滴撞击物体表面液膜后产生水花、发生飞溅，特别是产生“钟形”水花等的流动现象。Rioboo[2]通过实验观察发现，液滴撞击液膜时随着撞击速度的变化会出现三种不同现象：液滴在液膜表面的沉积，撞击后形成冠状水花但不发生飞溅，形成水花并发生飞溅现象。Cossali[3]的实验结果表明液滴撞击液膜后形成的冠状水花上部外径小于下部直径，两者之间的差值随撞击速度的增大而增大。数值模拟虽然出现的比较晚，但随着计算机的发展而有长足进步。郭加宏等[4-6]采用单相格子 Boltzmann 方法分别模拟了两个液滴撞击同种薄膜的流动过程、垂直间隔的两液滴相继撞击液膜流动过程，以及二维液滴冲击流动液膜的含自由面的流动过程，分析了液滴的速度、两液滴的中心间距等因素的影响。严永华等[7]通过用格子Boltzmann方法（LBM）对液滴撞击液膜喷溅过程的模拟，得出随着We数和Re数的不断增大, 液滴撞击液膜后将产生铺展、喷溅以及喷溅并伴有小液滴脱落等不同现象,结果表明LBM 研究气液两相流动问题的可行性。Lee[8]采用Level Set方法对液滴撞击液膜的初期形态进行了二维数值模拟，发现冠状水花直径随液滴冲击速度和液膜厚度的增大而增大，但从整体角度讲，冲击速度和液膜厚度对冠状水花直径的影响不大。郭加宏[1]认为液滴冲击液膜时，液滴所具有的动能大小是产生水花和飞溅的主要原因，当液膜厚度增加，液滴冲击产生的飞溅现象减弱，液滴直径不同，但液滴冲击的We数和Re数相近时，液滴冲击薄液膜产生水花的流动现象相似。沈胜强等[9，10]模拟了液滴撞击固体壁面的行为，另外分析了液滴撞击液膜时撞击速度和液膜厚度的影响作用[11]。国内外学者对液滴撞击液膜现象的研究总体较少且主要集中于实验研究，同时由于实验条件的不同，结论也存在差异，尤其对于液滴撞击液膜后形态演变的微小细节特征难以悉数捕捉，因此有必要采用数值方法模拟液滴撞击液膜的形态演变，以揭示其内在机理。本文采用CLSVOF（Coupled Level Set and Volume Of Fluid）方法对液滴撞击液膜后的形态演变进行了二维数值模拟，捕捉到了数值模拟中不常见的气泡夹带现象，描述了不同撞击速度下液滴撞击液膜后形态演变的过程，分析了液膜厚度对冠状水花扩展半径的影响，揭示了扩展半径随时间变化的规律。1 物理数学模型液滴撞击液膜的形态演变过程是一个复杂的气液两相界面变形的两相流动过程，对该过程模拟的关键和难点问题是如何方便准确的捕捉两相流界面。CLSVOF方法是将VOF方法与Level Set方法有机结合，从而有效解决上述难题的一种方法，它利用VOF方法良好的守恒性和Level Set方法对界面曲率和法向量的精确计算，捕捉两相界面。主流场控制方程为：（1）（2）其中为界面曲率，可由下式计算：（3）定义为：（4）其中w为最小网格尺寸。图1 初始时刻模型如图1所示，假定撞击初始时刻液滴与液膜相切，液滴为球形，液滴直径d=2mm，液膜厚度h分别为0.6、1.0、1.4mm，液滴撞击速度v分别为2.75、3.43、4.32m·s-1。无量纲参数We数、Re数、Oh数、无量纲液膜厚度h*、无量纲半径r*和无量纲时间t*定义如下：，，，，，，（5）其中，、和分别为液体的密度、表面张力系数和动力粘度，t为时间。则三种液膜厚度对应的无量纲液膜厚度h*为0.3、0.5、0.7。液滴和液膜均采用水作为介质。采用轴对称计算，计算区域设为10mm×10mm，网格采用四边形均匀网