计算域及流动条件:
管道长度为5m,直径为0.5m,该管道为低水深,即 hL/D = 0.14。流体注入速度为0.2m/s,空气注入速度为20m/s,湍流状态雷诺数 ReG = 1.6×105。
图1 管道网格划分
VirtualFlow 计算模拟得到的界面边界层空间演化结果如图2所示,可以看到剪切界面上气液两相存在显著的速度差异。
(a) 涡度
(b) 速度
(c)涡粘性的等高线
图2 界面边界层的空间演化
从图2可以清楚地看到界面上边界层的演变,可以认为本模拟采用的参数取值已经使流动获得充分发展。计算结果说明VirtualFlow能较好地预测由剪切力引起的气液界面发生了碎裂。
图3 三维涡量等势线
图4 正负涡量等势图(横截面图)
图2至图4表明在剪切面处存在小振幅的扰动,并最终导致液滴分离,这样的现象也可从以下计算结果中看到,如图5。
图5 环状流中出现雾滴夹带现象
图2-图4可以观察到界面上存在小振幅扰动,并最终导致液滴脱离。这种不稳定性(Kelvin-Helmholtz型)在x/D =1.4(分隔板下游55厘米)开始出现,并且第一次的液滴分离仅出现流动方向1.5m处,即对应x/D=3。这些液滴随后被空气夹带流动,可以很直观的从图4观察到该现象。
图6 液滴形成及环状流界面形变(横截面图)
图6为VirtualFlow软件计算得到的错流截面图,展示了液滴形成及环状流的界面形变。等势线对应于流动方向的速度,并且随着接近界面而衰减,液体速度由于太小而未被明确显示。