固液界面边界滑移对液滴振荡行为影响的VOF模拟.docx

该【固液界面边界滑移对液滴振荡行为影响的VOF模拟 】是由【科技星球】上传分享，文档一共【11】页，该文档可以免费在线阅读，需要了解更多关于【固液界面边界滑移对液滴振荡行为影响的VOF模拟 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。固液界面边界滑移对液滴振荡行为影响的VOF模拟
 
 
尹宗军，张椿英，苏 蓉，徐 辉，陈倩楠
(安徽信息工程学院 机械工程学院，安徽 芜湖 241100)
0 引言
流体力学表明，，而滑移长度λ是一个假设的边界内距离，在该距离处，[4]综述了液体流动的边界滑移及其相关问题的发展过程，，相关文献表明粗糙壁面、微通道、外部气体环境、高压高剪切以及驱动压强差对边界滑移长度都具有较大影响[5-15].赵士林[16][17]采用耗
散粒子动力学的方法，[18][19]研究了微通道疏水表面的滑移现象，，即疏水性越强，壁面滑移越明显，[20]，研究了五种状态的Cassie态表面上表面形貌、气液接触比、流道高度与驱动压强差对边界滑移长度的影响.
近年来，由于实验研究成本过高，、CLSVOF(Coupled Level Set and Volume of Fluid)、相场动力学等方法分别被运用到液滴撞击领域中，，模拟手段对改变撞击参数更为简单，且能给出液滴内部的动力学信息，例如压力、速度等，，已有的研究多数是基于静态网格，也就是网格化后，网格不再变化，，就需要细化整个计算域，，在物理量变化剧烈的地方自动细化网格，，当下迫切需要解决的技术问题是：(1)如何使用模拟计算取代实验仍然能获得精确的结果；(2)如何在模拟计算中用最少的网格精确获得界面形变，，利用动态网格自适应手段减少不必要的计算，且在相界面处自动加密，，本文基于VOF界面追踪和自
适应动态网格技术探究了液滴撞击壁面的振荡行为，进一步揭示固液界面边界滑移对液滴振荡特性的影响.
1 控制方程与参数设置
控制方程
在VOF模型中，不同的流体组分共用一套相体积分数方程，通过引进相体积分数这一变量，，假设二者之间为不可压缩、互不融合的层流运动，定义液相液滴为主相(ρL和μL用于表示液滴的密度和黏度，相分数设为c=1)，气相空气为第二相(ρV和μV分别表示空气的物理特性，相分数设为c=0)，而c介于0和1之间时代表两相界面区域；即可将两相流动等效为一种混合物(密度ρ和黏度μ).该混合物的密度ρ和黏度μ可以用体积分数变量c分别表示为ρ=cρL+(1-c)ρV，μ=cμL+(1-c)，流动控制方程包括连续性方程，带表面张力的动量方程以及VOF相分数方程[21]：
▽·u=0,
(1)
(2)
(3)
式中，u≡(u,v)为速度矢量，p为压强，D定义为Dij≡(∂iuj+∂jui)/2，g为重力加速度，δ为迪拉克算子，κ表示界面的平均曲率，n表示从液相流出的界面的单位法线.
离散方程
▽·un=0,
(4)
(5)
(6)
上述离散方程可进一步采用古典时间分解投影算法简化为[22]：
▽·un+1=0,
(7)
(8)
(9)
(10)
式中，ρi为密度，ui、μi、Di以及σi为对应时间节点i(i可取n、n+1/2和n+1)的值，u*(7)式和(10)式，我们有以下Poisson方程：
调查结果显示，有7%的学生参加了3项以上的课题研究，28%的学生参与了1～2项的课题研究。硕士研究生参与课题研究对其科研水平的提升有很大的帮助，经常参与科研课题，能够激发学生的学习积极性和科研创新能力。
(11)
GERRIS软件是求解描述流体流动的偏微分方程的自由软件程序，：(1)求解含时不可压缩变密度Euler、Navier-Stokes方程；(2)自适应网格细化,分辨率根据流体的特征动态调整；(3)时间空间二阶离散精度；(4)可以灵活的附加源项；(5)支持MPI并行、动态负载平衡、并行脱机可视化；(6).(11)式使用GERRIS中的基于四元/，(8)式可以重新改写为
(12)
这是一个Helmholtz型方程，可以使用GERRIS中的多级Poisson解算器求解.
边界条件
(13)
滑移边界条件和无滑移边界条件，这两种边界条件的法向流速都为零(即v=0)，区别在于这两种边界切向速度u是否为零；即无滑移条件就是指在靠近
壁面处流体的流动速度为零u=0，
u=u0+λ∂u/∂y，
(14)
计算域与初始条件
首先在GERRIS软件中定义理想的计算域尺寸，液滴撞击壁面过程可以简化为一个Lx×Ly=2L×L二维计算域，L=8 mm，如图1(a)、左侧以及上侧边界指定为压力出口边界，(液滴，直径D0)，液滴距离初始壁滴距离h=D0，以便液滴周围的气体流动能够在液滴撞击表面之前以物理方式发展，即保证将插入动量中的计算表面张力项始终具有最高的精确度.
液滴在铺展过程中液滴的形态是时变的，这就要求网格分辨率沿液滴/气体界面变化，因此动态自适应笛卡尔网格细化技术所用网格细化的判据可以表示为
‖c‖Δx/max‖u‖>δ,
(15)
式中，‖c‖是VOF体积分数变量c的值，Δx为子网格的尺寸，max‖u‖为流体穿过局部单元的最高速度，δ为阈值参量().若(15)式满足，，若满足下式
‖c‖Δx/max‖u‖
(16)
GERRIS软件自动将网格粗化一级，，若给定加密等级为10是指对构成计算的子正方形格子进行210=1024次切分，也就是说如果正方形格子边长为1 mm， ，加密等级越高，网格越密，计算精度越高，；其余区域加密等级选择6，故相应的计算单元尺寸为L/210和L/26(×10-×10-1mm).初始的网格使用GERRIS软件生成.
初始时，压强p均设置为0，(b)所示，应用的数值网格的放大区域显示了移动的液滴以及网格的高分辨率，(c)给出了液滴撞击过程中(i)碰撞、(ii)铺展以及(iii)反弹三种运动特征，表明网格能自适应地跟随液滴.
2 结果与分析
液滴撞击光滑固体表面的运动特征
我们将水滴(D0= mm，V0= m/s)与光滑固体超疏水表面(静态接触角θE=157°)的VOF模拟与文献[23]，水和空气的密度和黏度分别设定为ρL=×102kg·m-3，ρV= kg·m-3；μL=×10-3N·s·m-3，μV=×10-5N·s·m-3，表面张力设置为γLV=×10-2N·m-.
图1 液滴的滑移边界条件及运动特征
图2 液滴撞击光滑固体表面的运动演化
图2(a)和(b)所示分别为VOF模拟与文献[23]的液滴撞击光滑固体表面的运动演化图，图2(c)显示了铺展因子β的时间变化曲线(定义为β=L/D0，L为当前的润湿长度).当液滴接触壁面时，= ms为水滴刚接触到固体表面时刻，此时β= ms后，液滴底部出现褶皱型，由于疏水表面的黏附力影响，，在t= ms时，液滴底部达到最大润湿面积，最大铺展因子βmax== ms时，液滴呈现扁平状，，由于固体表面的阻挡，液滴的法向动量发生变化，液滴周向铺展以扩大表面积，直到动能完全转化为表面能量；这种不平衡状态是暂时的，在表面张力的作用下，液滴已经周向铺展的流体会向中心回缩， ms内，这种反向惯性力迫使湿滴面积减小，= ms时，(c)进一步表明液滴撞击光滑固体表面过程中模拟和实验的液滴铺展因子β变化具有一致性.
固液界面边界滑移对液滴振荡的影响
为了使我们的结论更具普遍性，我们使用了两组(V0= m/s,θE=70°和V0= m/s,θE=110°)，在VOF计算中我们使用了五种不同的边界条件：λ= mm、 mm、 mm、 ，选择直径D0= mm的水滴作为仿真模拟对象.
图3 不同滑移边界条件下铺展因子β的时间演化
图3显示了在不同滑移边界条件下，当V0= m/s,θE=70°(亲水表面，近毛细管铺展)和V0= m/s,θE=110°(疏水表面，近反弹)时，(i)，这有利于描述液滴的最大铺展因子βmax；图3右侧用于表示100 ms内液滴的振荡行为(ii)，这有利于描述液滴振荡行为的基本物理量，，无论是近毛细扩散还是近反弹，，(a)-(i)和图3(b)-(i)均表明，在第一个振荡周期下，