曹义刚 付萌萌 杨喜昶 李登峰 王晓霞

(郑州大学物理学院,郑州 450001)

Kelvin-Helmholtz 不稳定性(KHI)增长的动力学分析是一个活跃的研究领域.本文解析研究了流体在横截面不同的直管道中流动时,热传导对KHI 的影响.结果表明:管道中上下流体的界面相对切向速度会随着波数的增加先增加后减小,并且小的界面热传导系数导致相对切向速度随波数的减小更多,不同于横截面相同的直管道结果.另外,热传导会提高KHI 的增长率,与横截面相同的直管道一致.研究结果可以为实际管道中流体不稳定性的分析以及管道的通风设计和供暖等工程研究提供一定参考.

1 引言

流体不稳定性是流体力学的一个基本问题,通常发生在不同流体或者同种流体不同密度的界面.常见的界面不稳定性有Reyleigh-Taylor 不稳定性(RTI)[1,2]、Richtmyer-Meshkov 不稳定性(RMI)[3,4]和Kelvin-Helmholz 不稳定性(KHI)[5,6].RTI 发生在密度梯度和加速度方向相反的流体界面,当冲击波通过两流体的界面时会产生RMI.RTI 和RMI均垂直于界面方向,而KHI 则是由平行剪切作用于界面时引起的.上述不稳定性在天体物理[7]、等离子体[8]、超流体[9]、磁流体[10]以及惯性约束聚变(ICF)[11−13]等领域中发挥着关键作用.

理解任何一种流体不稳定性不仅对该不稳定性本身很重要,而且对其他流体不稳定性的理解也很重要.事实上,KHI 不仅会加剧非线性RTI 和RMI 的发展,是界面蘑菇状结构演变的关键,而且对由RTI 和RMI 产生的流动向湍流的过渡起着至关重要的作用[14,15].

KHI 线性不稳定性增长的动力学分析是一个活跃的研究领域.在过去的几十年里,人们广泛研究了各种因素对KHI 线性演化的影响.这些因素包括压缩[16]、黏性[17]、表面张力[18]、热传导[19,20]、倾斜界面[21]等.热传导现象是由物质中大量的分子热运动互相撞击,从而使能量从物体的高温部分传至低温部分,或由高温物体传给低温物体的过程.热传导效应在KHI 中不可忽视.之前人们[19]在横截面不变的直管道中通过将热传导项加入到速度势中,求解管道中上下流体界面的伯努利方程,从而得到了界面热传导系数与上下流体的界面相对切向速度和KHI 增长率的关系.人们发现:增加界面热传导系数会抑制上下流体的界面相对切向速度,但会提高KHI 的增长率[19,20].实际上,管道的横截面和管道中的流体速度均会发生改变[22].管道横截面或流体速度的突然变化均会造成流体界面处的能量损失,从而引起水头损失[22−24].水头损失一般分为沿程水头损失和局域水头损失,而由管道横截面突变产生的水头损失为局部水头损失[22].另外,流体流动中,由于管道的横截面发生变化,流体内部各粒子的动能和势能相互转换,流体内部的热传导(这里主要是热导率)也会受到影响.通过计算热传导对KHI 的影响,可以有效地设计日常通风供暖等流体管道.

目前,人们对横截面变化的直管道中(或者水头损失存在下)热传导对流体不稳定的影响研究甚少.本文将解析研究热传导对横截面不同的直管道中KHI 的影响,并与横截面相同的直管道结果相比较.论文安排如下:第2 节给出模型;第3 节对结果进行分析和讨论;第4 节给出结论.

2 模型

假设一个水平放置的直管道,如图1 所示.管道具有两种不同的截面A1和A2,其中A1为窄管道的横截面,A2为宽管道的横截面,这两部分的管径分别是D1和D2.管道中有两种不同的流体,密度小的流体在上,密度大的流体在下,两种流体之间存在一个界面(y=0).假设上下流体密度、黏度、平均流速和温度分别为ρu和ρl,µu和µl,Uu和Ul,Tu和Tl.此外,还假定流体不可压缩且两种流体之间没有旋转和摩擦,即两种流体存在势流以及流体之间无耗散.在流体流动过程中,水头损失主要由管道横截面的突变引起.

图1 两种不同横截面的直管道示意图Fig.1.Schematic diagram of a straight pipe with two different cross sections.

上下流体的界面满足:

式中函数F(x,y,t)中;t为时间;η为解界面扰动振幅,表示为

其中η0为常数,c.c.代表复共轭.这里假定管道壁面为光滑壁面.对于上下均匀流体,小的界面扰动振幅满足速度势[19,25]:

式中,k,ω为扰动的波数和频率;α′为界面热传导系数,定义为界面处的热流强度与单位长度的潜热之比[19];hu(hl)为流体界面到上(下)壁面的距离.

因为流体是无旋的,所以速度势(3)式和(4)式满足拉普拉斯方程[26]:

流体界面的控制方程为

流体界面的动力学方程为

式中g为重力加速度,γ为界面张力系数.上下流体之间的正压力为Pu−Pl,可以表达为[20]

局域水头损失可表示为[22]

式中水头损失系数ζ与管道两部分的截面比A1/A2有关[27],即ζ=(1−A1/A2)2;uu,x和ul,x分别是上、下流体在产生水头损失时的瞬时速度.uu,x和ul,x可由速度势(3)式和(4)式的负梯度得到

将(3)式、(4)式和(9)式—(13)式代入(8)式,可以得到

通过(14)式,可得到与ω有关的色散方程:

对照(15)式与(14)式,可得到各项系数为

令ω=ωR+iωI,可得

联立方程(22)和方程(23),得

再将ωR代入到方程(22)中,就得到关于增长率ωI的四次方程:

这里

KHI 是由于扰动的增长率大于零时引起的.如果增长率小于零,扰动将随时间快速衰减.然而,当增长率大于零时,扰动随时间呈指数增长[18].因此只有当增长率为零时界面稳定.KHI 的增长率由扰动频率的虚部直接表征.令ωI=0 (即稳定界面),并将其代入(25)式,于是有

将(16)式—(20)式代入(31)式,可以得到不同传热系数下,上下流体界面相对切向速度与波数之间的关系,见(A1)式,其中V=Uu−Ul是上下流体的相对切向速度.

在考虑KHI 的增长率时,由于有黏和无黏情况下得到的结果几乎是一样的[21],因此令µu=µl=0,并且将(16)式—(20)式代入(25)式中,就可以得到不同传热强度下,KHI 的增长率随波数的变化关系,见(A2)式.

3 结果分析与讨论

为了定量分析热传导效应,取扰动振幅η=0.03 m,管径D1=0.05 m,其他参数取自参考文献[25]:γ=0.06 N/m,ρu=1.2 kg/m3,ρl=1000 kg/m3,µu=0.00018 Pa·s,µl=0.01 Pa·s,hu=0.025 m,hl=0.025 m.界面热传导系数α′=α×103kg/(m3·s)的取值参照文献[19],这里α是无量纲界面热传导系数.此外,由文献[22]可以看出,在水头损失系数等于0 和0.2 时,能够明显得到管道横截面的变化对KHI 的影响效果.因此,图2 给出了ζ=0,0.2 时,不同界面热传导系数下上下流体的界面相对切向速度与波数的关系曲线.从图2可以发现:上下流体的界面相对切向速度V随波数k的增加而增加,当波数达到一定值的时候,界面相对切向速度达到最大,尔后随着波数的继续增加相对切向速度反而会降低.再者,我们发现,当ζ=0时,即不存在水头损失(或管道横截面不发生变化)时,相对切向速度随着波数增加,最终趋于一致.而当ζ=0.2 且相对切向速度开始随波数增大而减小时,小的界面热传导系数会比大的界面热传导系数对相对切向速度的抑制更大.当ζ=0.2时,界面相对切向速度要低于ζ=0 时的相对切向速度.这是由于当流体流经横截面变化的管道时能量会损失,从而导致较低的相对切向速度.

图2 上下流体的界面相对切向速度与波数的关系Fig.2.Wavenumber dependence of the relative tangential velocity at the interface between upper and lower fluids.

图3 给出了ζ=0,0.2 时,不同界面热传导系数下KHI 的增长率随波数的关系曲线,其中参数取自参考文献[19]:γ=0.06 N/m,ρu=1.2 kg/m3,ρl=0.001 kg/m3,hu=0.025 m,hl=0.025 m. 我们发现,当ζ=0 时,增长率随波数的增长先增大后减小并最终趋于不变.这是因为界面动能随波数增加而增大,导致增长率增大,当波数接近临界值kc时,KHI 被破坏,增长率随波数的增加而急剧减小,趋近于零[19,22].当ζ=0.2 时,即存在水头损失(或管道横截面发生变化)时,KHI 的增长率要大于不存在水头损失时的增长率,这是由于在管道横截面变化时,能量损耗的一部分转化为内能,促进界面热传导对KHI 增长率的影响.比较不同热传导系数的曲线,可以明白这一点.

4 结论

本文建立了一个横截面变化的直管道模型,并分析了热传导对在该管道中流动的上下两种流体KHI 的影响.结果表明,在管道横截面变化的情况下增加界面热传导系数会抑制界面相对切向速度,与文献[19,20]中报道的横截面不变的直管道情况一致;与之不同的是,当波数到达一定值,相对切向速度开始随着波数的增加而减小,并且小的界面热传导系数在波数较大的情况下对相对切向速度的抑制效果更大.另外,增加界面热传导系数会促进KHI 的增长率,这与文献[19]中报道的横截面不变的直管道的情况一致,但是相对于横截面不变的直管道,热传导对KHI 增长率的促进效果更为明显.

值得指出的是,本文没有考虑不同横截面处流体的密度、黏度以及平均流速的不同对结果的影响.考虑上述不同可能会对结果造成较大影响.另外,尽管本文的结论是由圆形横截面的管道得出,但是也适用于非圆形横截面的直管道.本文的结果对实际管道中流体不稳定性的研究以及管道的通风设计和供暖等工程领域具有一定的参考价值.

附录A 上下流体界面相对切向速度和KHI 的增长率与波数的关系

上下流体界面相对切向速度与波数之间的关系满足以下方程:

KHI 的增长率随波数的变化关系满足: