剪切流动条件下液滴变形的数值模拟研究.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：056074 剪切流动条件下液滴变形的数值模拟 林长志郭烈锦“ 西安交通大学多相流国家熏点实验室，西安，710049 TEL：029- Email：￡!li鲤!丛：!j地：￡血：￡! I蔓：奉文通过数值方法研究剪切流动条件下悬浮液滴变形的动力学问题。采用扩散界面和连续表面 张力模型模拟液滴界面瞬形状及界面的破碎。流体问的界面采用质量浓度直观表达。瞬态、小可压缩 end 研究了稳态液滴形状及液滴破碎的Evolvingpinching机理。计算表明奉计算方法刚以很好的捕 获液辅变形的动力学特征并与试验观察一致。 关t宇：液滴，界坷，变形，破碎： 1、前言 液滴动力学现象广泛存在予大量的工业应用中，如石油化工行业中的萃取、分离和 输运过程。液滴的生成增加了接触界面的面积从而提高了相之间的质量和热量的传输。 准确掌握界面的运动变化机理对于工业控制和优化设备性能具有重要意义。但是由丁．界 砸运动的1}线性及密度、粘度、表面张力等物理量的影响，对界面的运动变化难以准确 的进行理论和试验研究，因此数值模拟是研究液滴界面动力学现象的有力工具【1，2】。 界面变形的动力学问题是一个完全的三维问题，准确的界面描述是数值模拟的笑键。 界面的数值模拟方法主要有界面追踪和界面扑获。界面追踪对界面上的标志点进行追 踪，包括边界积方法【3】、前沿追踪法【4】等。边界积分法的优点是可以将计算的维 数降低一维．适用于无粘性或Stokes流动条件下小变形及界面模拟，但是数值不稳足 性敏感，当界面出现夹断或聚并时不能准确：I：作。前沿追踪法通过设置标志点描述界面 的运动，具有较高的准确性，但是当界面变形较大时。标志点的增加或删减会产生数值 的不稳定。界面扑获法通过一个标照函数自动的扑获界面，包括水平集法【5】、流体 基本思想是认为界面具有非零的厚度，通过引入一个具有物理意义参数(如质量浓度) 来表征界面，从而实现界面的追踪。扩散界面法是模拟界面动力学的问题的有力C具， 等。本文通过采用扩散界面方法数值研究悬浮液滴的变形问题，引入质量浓度追踪界面 的运动变化。 2、物理问墨 2．i问题描述 本文研究悬浮丁不相溶、不可压缩体相液体中孤立液滴在剪切流动条件1．-的变形问 407 题，如图1所示，两种流体不可压缩，互不相溶。液清与体相流体的粘度分别为锄、‰ 密度相同(岛=户k)a液滴的初始半径为玎。计算域在J，，方向为周期边界条件，上 下壁面为无滑移边界条件。，为施加的剪切率，Er为界面张力。 图i．剪切条件下液滴变形动力学筒图 2．2数学模型 液滴的运动变化过程是完全的三维动力学问题，需要直接数值解三维控胄0方程。本 ll 文采用扩散界面方法，通过Cahn—Hi 态流动条件卜液滴与体相流体的运动，表面张力的计算采用continuumsurfaceforce V·“=0 (1) (2) V。[节(c)(V卅V计7)]+只 u,+u-Vu=-Vp十丙1 q+”吼=去V’∞(c)耻) (3) ∥=厂(c)一CAc (4) 只是界面张力只一嚣阢隔j|VcIVc- 2．3数值计算方法 格算法．解Navier—Stokes方程的近似投影法和计算表面张力的CSF算法。 (1)非线性多重网格算法 f{{均匀矩形交错网