VOF耦合Level-set方法模拟液滴演化及静电场下锥射流的形成讲义.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：1206020 VOF 耦合 Level-set 方法模拟液滴演化及 静电场下锥射流的形成 王军锋 * 左子文 霍元平 （江苏大学能源与动力工程学院，镇江 212013） (Tel:0511 Email:wangjunfeng@) 摘要：液滴的演化作为雾化的最基本形式，研究其动力学特性为深入探究静电雾化机理提供了重要的 理论依据。静电雾化的重要形式之一是锥射流，研究其形成过程及电动力学特性对于探讨分析多股射 流等其他复杂多相流具有重要意义。利用 Level-set 耦合 VOF 方法，在假设流体不可压缩并做平流运 动的情况下，对单个液滴生成过程以及静电场下锥射流的形成进行了数值模拟。模拟结果与实验结果 较好吻合，验证了 Level-set 耦合 VOF 方法模拟两相流问题的有效性。同时分析了液滴生成过程及锥 射流状态下的速度场等分布，为模拟多场耦合等复杂条件下的液滴二次雾化奠定了基础。 关键词： Level-set 方法；VOF 方法；液滴形成；静电场；锥射流 0 前言 喷雾广泛应用于日常生活、医学、工农业生产、消防以及航空航天等新型领域。静 电雾化与常规雾化相比，荷电雾化具有显著的优点，如液滴粒径更小，沉积率高，可控 性强等，因此受到人们越来越多的关注。近些年来，人们提出了多种数值方法来模拟研 究常规或非常规环境下液滴破碎、气泡生长等涉及气液、液液的两相流问题，比如 Front Tracking 方法[1]、Marker Particle 方法[2]、VOF 方法[3]，以及 1988 年提出的 Level set 方 法[4]。尤其是 VOF 方法与 Level set 方法易于处理界面的拓扑变化而得到广泛使用。 VOF 方法基于固定欧拉网格计算和追踪每个单元特定相的体积率。其优点在于能够 很好的保证物理量的守恒。但是却很难准确计算多相流交界面的法向和曲率。并且，VOF 函数在交界面是不连续的，这增加了它向高维推广的难度。Level-set 方法利用定义的符 号距离函数捕捉追踪两相流交界面。Level-set 函数始终平滑连续，易于处理复杂界面变 形及其拓扑结构改变[5,6,7]。但是该方法在计算时有物理量损失，不是一种守恒方法。 本文利用 VOF 方法确定交界面单元、Level-set 方法捕捉交界面的耦合方法来模拟 液滴的演化及静电场下锥射流的形成，并与实验进行对比，验证了 Level-set 耦合 VOF 方法对于模拟液滴脱落破碎气液两相流问题的有效性，为后续的多场耦合等复杂环境下 液滴破碎的模拟研究打好基础。 1 Level-set 耦合 VOF 方法 对于两相流场， ),( txvα 满足 @#$%I1507010002ZX00*^% 0= ? ? + ? ? + ? ? y v x u t ααα (1) 网格上流体体积函数 Cij 为 ),( tx rα 在网格上的积分： dVtx V C ijlij ij ∫Δ= ),( 1 vα (2) C=1 为满网格，C=0 为空网格，0C1 为半网格。流体体积函数 C 满足： 0= ? ? + ? ? + ? ? y Cv x Cu t C (3) 这个方程就是 VOF 方程。 Level-set 函数定义为符号距离函数?，构造零等值面 ),( tx? 为两相流交界面 )(tΓ ， 即： }0),(:{)( =Ω∈=Γ txxt ? 。 在两相流场中： ? ? ? ? ? Ω∈ Γ∈ Ω∈+ = 2 1 ,- )0( ,0 , ),( xd x xd tx? (4) d 表示 x 到交界面 )(tΓ 的距离， 1Ω 和 2Ω 分别表示第一种介质和第二种介质所处区域。 计算表面张力过程中需要法线的方向及曲率由下式估算得到： 0= ? ? = ? ? ?nv (5) 0= ? ? ??= ? ? ?κ (6) Level-set 函数可改写为类似 VOF 方程的函数： 0)( =??+ ? ? ?? u t v (7) 其中，uv为速度场。动量方程为： [ ] guu puu t u T vvv vv v ρ??σκδμ ρρ +???+??? +??=??+ ? ? )()( )()( (8) @#$%I1507010002ZX00*^% 其中，如果 )5.1(, h= αα? ，（h 为网格间距），则 α απ??δ 2 )/cos(1)( += ，否则 0)( =?δ （σ 为表面张力系数）。 在计算方程(7)的过程中，由于误差原因，不能维持 1=?? ，得到的 ),( txφ 可能不 再是 t 时刻点 x 到界面的距离，从而导致迭代过程中误差的累加。所以引入几重新初始 化，首先通过 VOF 模型