第九章 绕流与射流重点阐述不可压缩粘性流体绕流二维和回转物体绕流现象及其绕流阻力的计算,分析工业生产中常遇到的紊流射流问题。
§9-1 绕流阻力与阻力系数当粘性流体绕流物体时,物体总是受到压力和摩擦力的作用。
作用在整个物体一表面上的压力和摩擦力的合力F 可分解为两个分力,即绕流物体的未受干扰时来流速度∞V 方向上的分力F D ,及垂直来流速度∞V 方向上的分力F L 。
对于在静止流体中运动的物体来讲,由于F D 与物体运动方向相反,是阻碍物体运动的力,故称之为绕流阻力;F L 称为绕流升力。
于是D L F F F +=绕流阻力和升力二者都包含摩擦力和压力两个分量,因此,物体所受摩擦力和压力的大小及二者的变化是分析绕流阻力的基础。
一、绕流阻力一般分析物体壁面所受摩擦阻力是粘性直接作用的结果,所受压力又称压差阻力,是粘性间接作用的结果,当粘生流体绕流物体时,边界层分离是引起压差阻力的主要原因。
下面以圆柱绕流为例来说明绕流阻力的变化规律。
在绕流未分离的情况下,由理想流体所确定的物面上的压强分布如图6-12所示,在第六章的第四节详细地讨论过这个解,物体所受压力阻力为零。
在绕流圆柱体发生严重分离的情况下,由于柱体后部背流面存在分离区,此时主流区的边界处在分离区的外缘,柱面上的压强分布不同于未分离时的压强分布,从分离点开始,柱体后部受到的流体压强大约等于分离处的压强,而不能恢复到理想流体绕圆柱体流动时应有的压强数值,从而产生对圆柱体的压差阻力。
图9-1(b)所示是有边界层分离的圆柱面上的无因次压强分布,实验曲线见图6-12中的II 、III 曲线。
对于摩擦阻力,其形成过程比较清楚。
实验表时,象机翼、船只和其它一些流线型物DF体都有较大的摩擦阻力。
钝体如圆柱、球、桥墩和汽车等都有较大的甚至压倒优势的压差阻力。
由于压差阻力的大小与物体的形状有很大关系,因此,压差阻力又称为形状阻力。
二、阻力系数虽然绕流物体阻力的形成过程从物理观点看完全清楚,但要想从理论上通过面积分求解一个任意形状物体的阻力是十分困难的,目前都是由实验测得,工程上习惯借助无因次阻力系数来确定总阻力的大小,即A V C F 2DD 21∞=ρ (1) A V F C 2D D 21∞=ρ (2)式中A 为物体的投影面积,当物体主要受压差阻力时,采用物体垂直于来流速度方向的投影面积,即迎流面积。
物体的阻力系数的大小,主要取决于雷诺数Re 的大小和物体的形状,也与物体在流场中的方位密切相关。
由相似定律知道,对于不同的不可压缩流体中的几何相似体,如果Re相等且在流场中的方位相同,则它们的阻力系数相等。
因此,在不可压缩流体中,对于与来流方向具有相同方位的几何相似体,如果Re 相等且在流场中的方位相同,则它们的阻力系数相等。
因此,在不可压缩流体中,对于与来流方向具有相同方位的几何相似体,其阻力系数C D 只是Re(vlV ∞=Re ,式中l 为特征尺寸)的函数,即 ()Re D f C = (3)图9-2和图9-3给出了无限长圆柱体和球体阻力系数与Re 的实验关系曲线。
由图可知,在不同的Re 下,流动现象的差异和阻力系数的大小是明显的。
下面仍以圆柱为例,具体分析随着Re 的变化绕流现象的变化过程及阻力系数的大小。
(1) 在Re ≤1的范围内,流动如图9-4(a),边界层没有分离,其特点为圆柱表面上下、前后流动对称且呈层流流态。
流动阻力来源于柱面摩擦阻力的合力,C D 与Re 成反比,如图中直线部分。
(2) 在3~5<Re<30~40范围内,流动如图9-4(b),其特点为在圆柱背流面发生边界层分离且形成对称驻涡区。
绕流阻力由摩擦阻力和压差阻力两部分组成,它们具有同等的重要性。
(3) 30~40<Re<60~90范围内,随着Re 的增大,圆柱背流面分离区逐渐变宽,对称涡(a)Re ≤1;(b)3~5<Re<30~40;(c)30~40<Re<60~90;(d)60~90<Re<1.5×105;(e)Re>1.5×105区出现摆动,如图9-4(c)所示。
此范围内,压差阻力在总阻力中占的比例逐渐增大,虽然摩擦阻力和压差阻力有同等重要地位,但压差阻力已逐渐占主要部分。
(4) 在60~90<Re<1.5×105范围内,圆柱背流面旋涡交替脱落,形成两排向下游运动的涡列,称之为卡门涡街,如图9-4(d)所示。
随着Re的增大,边界层分离点最终移向迎流面,该流动情况称为绕流物体的亚临界状态,如图9-4(e)所示。
绕流阻力仍由摩擦阻力和压差阻力组成,但压差阻力已明显地占绝大部分。
(5) 在Re>1.5×105条件下,随着的Re的增大,分离点前的边界层由层流转变为紊流,紊流边界层的强烈混合效应使得分离点向后移,尾涡区变窄,绕流得以改善,如图9-4(f)所示。
该流动情况称为绕流的超临界状态,阻力系数亦随之下降,即从Re=1.5×105~5×105,C D从1.2急剧下降到0.3,如图9-2所示,由于物体阻力以压差阻力为主,故绕流总阻力是下降的。
应指出,层流边界层可以人为地转变为紊流边界层,即亚临界状态可人为地提前转变为超临界状态。
转变的办法是扰动来流,增加上游的紊流度,普朗特曾用下面的实验证实了这一现象,他在紧靠圆球上层流边界层分离点的稍前面套上一圈细金属丝,人工地把层流边界层转变为紊流边界层,则Re在小于3×105的亚临界时,阻力就显著下降,此时分离点从原来的圆球前驻点后约80°处向后移到约110°~120°。
通过上述分析,可归纳如下:为了减小绕流阻力,应设法避免边界层分离。
在不发生边界层分离的情况下,边界层应尽可能保持层流。
如果已发生分离,则应在发生分离处稍前设法使边界层流转变为紊流,以使分离点后移。
由于发生边界层分离的条件是沿流向压强增大和动能不足,因此,如果在压强增大处采取减压措施或使流体增加流速,均能消除边界层分离而降低阻力,这种措施称为边界层控制,这里不再介绍,可能参阅有关文献。
§9-1 卡门涡街不可压缩粘性流体绕流圆柱体,在60~90<Re<15×105的绕流条件下,圆柱背流面将产生卡门涡街(Karman vortex street),本节的主要内容是分析卡门涡街的脱落频率以及由卡门涡街引起的声学共振现象。
一、卡门涡街的脱落频率1911年,匈牙利科学家卡门在德国专门研究了这种圆柱背后涡流的运动规律。
研究表明,圆柱背后生成的旋涡是交替脱落、在尾流中交叉排列,旋转方向相反的有规则的两行涡,故称为卡门涡街或卡门涡列。
卡门涡街的稳定条件为:两行旋涡间距与同行相邻两旋涡间距之比h /l =0.281。
如图9-5所示,此时的Re ≈150。
近年来,许多学者对卡门涡街进行了研究,认为当流体绕流单根圆柱体时,在250<Re<1.5×105范围内,柱后交替脱落的旋涡频率为⎟⎠⎞⎜⎝⎛−=∞Re 7.191198.0dV f (1)式中f ——涡街脱落频率;d ——圆柱直径;∞V ——来流速度;上式中的f 、∞V 、d 可组成一无因次系数,用St 表示,即∞=V fdSt (2) 式中的St 称作斯特罗哈数,是用该方面的早期研究者斯特罗哈(V ·Strouhal)名字命名的,根据罗斯柯(A ·Roshko)1954年的实验结果,在Re=800~1.5×105范围内,St 数近似等于常数,即St ≈0.21。
如图9-6所示。
绕流圆柱体形成的卡门涡街尾流,当Re 约大于300时为紊流尾流,卡门涡街不断消失在紊流尾流中,但当Re ≈3.3×105达到超临界时,有规则的旋涡脱落不再存在,这种情况一直要持续到Re=3.5×106,大于这个数值时,又会形成卡门涡街,这时的St 数为0.27。
根据卡门涡街的上述性质,可制成卡门涡街流量计,即垂直插入管道一根圆柱验测杆,管内流体流经验测杆时,在其下游产生卡门涡街,测得了涡街的脱落频率,便可求出管内流速,进而确定管内流量,测定卡门涡街脱落频率的方法有热敏电阻丝法,超声波束法等。
二、卡门涡街作用力对卡门涡街的研究尚发现,由于柱面上的涡以一定的频率交替脱落,因而柱面压力和切向力也以相同的频率发生规则的变化,产生交变的横向力,横向力的方向总是指向刚脱落旋涡的那一侧面。
这是由于刚脱落完旋涡绕流得以改善的柱面与准备脱落涡的柱面压力分布不同(见图6-12),前者无因次压强系数明显小于后者,尤其在小于亚临界雷诺数时,垂直于来流方向的横向变力比来流方向的横向交变力要大得多,后者可不考虑。
因此这个垂直来流方向的横向交变力可按下式计算t A V C F k k ωρsin 22∞= (3)式中F k ——卡门涡街作用力;C k ——卡门作用力系数,当Re=102~107时,C k =1; A ——圆柱投影面积;ω——圆频率,由卡门涡街脱落频率决定,ω=2πf ; t ——时间。
如果这个施加在圆柱(管)上的横向交变力与圆柱(管)的固有频率相等,就会引起柱(管)的共振,使其振幅越来越大直至破坏。
1940年华盛顿州的塔可马吊桥被风吹毁就是这个原因。
三、卡门涡街声振在日常生活中,常听到风吹输电线嘘嘘发响的呜叫声,这种呜响也是由于卡门涡街的交替脱落引起空气中压强脉动所造成的声波。
在工程设备中(例如管式空气预热器),空气横向绕流管束,卡门涡街的交替脱落会引起管箱中气柱的振动。
特别是当旋涡脱落频率与管箱中的声学驻波振动频率相等时,便会发生声学共振现象,产生严重的噪声,并使器壁在脉动压力作用下弯曲变形,甚至振裂。
最严重的情况是气室的声学驻波振动频率、管束的固有频率、卡门涡街的脱落频率三者相合时,将造成设备的严重破坏。
通常,消除声学共振破坏的措施是提高设备气室的声学驻波频率,也就是顺着流体流动的方向加若干块隔板,将设备气室的横向尺寸分隔成若干段,提高其声学共振频率,使之与卡门涡街的声振频率错开。
这种简单的方法实践证明是行之有效的,但具体做时要通过试验及必要的计算来解决。
上述讨论的着眼点是流体绕流单根圆柱(管)的情况,实际工程设备中都是排列有序的管排,这种情况下的卡门涡街脱落过程将更复杂,与管子排列方式密切相关。
参见曲线图9-7,可根据管束的排列方式查取St数,然后计算卡门涡街的脱落频率。
§9-3 球形物体的自由悬浮速度下面研究圆球在静止流体中的运动情况。
圆球在静止流体中由静止开始自由下落,由于重力作用而逐渐加速,同时阻力也随之加大,当圆球的重量等于作用在圆球上流体的浮力与阻力之和时,即满足:W =F B +F D圆球在流体中以等速V f 自由降落(类似跳降落伞,主要是压差阻力),临界速度V f 又称为圆球的自由降落速度。