第6章 层流的解析解与近似解(1)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第6章 层流的解析解与近似解
粘性流动基本方程组的解析解有着它固有的数学困难,真正能做解析解的流动为数不多,而且都是比较简单的流动。 本章将介绍几种粘性流动的解析解,有助于我们开阔思路,认识多种实际流动的性质。 首先先介绍一下粘性流研究的意义和研究的特点以及粘性流动的基本方程组,接着介绍一些解析解。 在介绍解析解时先考虑常特性不可压缩流体,通过基本方程,解得流场的速度和温度分布,最后求出摩擦阻力系数和热交换系数。 为了认识可压缩流动的特性,介绍两种简单的可压缩流动的解析解。 另外本章只限于雷诺数不大的流动。
6.1 粘性流研究的意义
一切流体都具有粘性,但是人类最经常接触的流体,如水和空气其粘性都很小,要考虑粘性的影响就会使数学问题变得非常复杂;另外,对于这些粘性小的流体,忽略其粘性所得到的结果又能在一定程度上符合实际情况,因此,理想无粘性流体理论最先得到了发展,它比粘性流体理论要成熟得多。
应当指出,虽然理想流体理论取得了重大的成就,但在某些方面却有不可逾越的先天性缺陷。 例如,它不能预估管道流动的压力损失,也不能计算在流体中运动的物体所受到的阻力。 后一问题与著名的达朗伯疑题有关。 达朗伯对理想流体进行了严谨的研究后得出了如下结论:当任意形状的固体在静止的充满无限空间的无粘性流体中作匀速直线运动,它不承受沿运动方向的作用力,即物体所受阻力为零。 在他所做假设的前提下,这一结论的逻辑推理是完全正确的,但它却与实际完全不符,因为所有的物体在流动中运动时都受到阻力作用。 这从反面说明了考虑粘性的必要性。
例1 圆柱绕流
对于理想不可压缩流体,
()
22
214sin s p p p C U θρ∞∞
-=
=- 其中 p ∞——远前方静压,ρ——流体密度。
图6-1给出了上述理想流体的压力系数与实际测量值的比较。 图中的实验曲线对应于两个不同的Re 数。
图6-1 圆柱表面的压力分布,理想流体理论与实验测量数据的比较
由图6-1可见,在圆柱的前缘(0οθ=和360ο
)附近,理想流体的理论结果与实际符合较好。 但在后缘
(180οθ=)附近两者差别则相当大。 对于理想流体,圆柱前后的流动是完全对称的,所以理论阻力为零。 但是实测的压力分布前后不对称,圆柱后部的实测压力系数低与前部对应点处的值,使圆柱受到向后作用的力,即
压差阻力。 另外,实际流体也引起表面摩擦阻力。 理想流体理论不能计算出这些阻力,这是它与实际流动情况的重要差别。
图6-2真实流绕圆柱的流动
由图6-1还可看出,理想流体结果与亚临界雷诺数流动的差别较大,与超临界雷诺数流动的差别较小。 实际上流体在圆柱体后部处于减速增压流动阶段,由于粘性耗散,使边界层内底层流体动能不断消耗,无力克服迎面高压。 这股流体将在该处与固体壁面脱离,这种现象称为边界层分离。 流体分离后,静压不易再有较大的回升,并在其后形成宽的尾迹,见图6-2。 在图6-1中实际流体在圆柱体后缘呈现出的低压区就是这样产生的。 分离点的位置以及尾迹流的宽度和特性取决于雷诺数Re 的数值。 亚临界雷诺数通常对应于层流流动,流体易于分离,而超临界雷诺数通常对应于湍流流动,流体有较强的承受逆压力梯度的能力,不易分离。 这就是图6-1中不同的雷诺数有不同的压力分布曲线的原因。
图 6-3 圆柱的阻力系数随雷诺数的变化
图6-3表示无量纲阻力系数()
2/D C F U R ρ∞=与雷诺数Re /UD ν=的关系曲线,其中F 为单位长度圆柱所受到的阻力,D 为圆柱直径。 由图可见,亚临界雷诺数时,阻力系数很大,随着雷诺数增加,阻力系数下降,在5
Re 510=⨯附近,阻力系数急剧降低,这对应于由层流边界层转变为湍流边界层。 阻力系数的这种变化与图6-1中压力系数分布随雷诺数的变化是一致的。
例2 二维机翼绕流
二维机翼是指沿展向无限长,且翼型不变的机翼。
圆柱绕流是非线性体的典型例子,机翼绕流则是流线型体的典型例子。
图6-4给出了儒科夫斯基翼型表面的压力分布。 这是在理想流体与实际测量有相同的升力条件下进行的比较。 由图可见,这里的理想流动的结果比圆柱绕流的情况好得多。 几乎沿翼型的整个表面理想流体的结果都与实验符合,只是在翼型的尾部的上表面有较大的差别,这也是沿流动升压使边界层分离的结果。 图6-5给出了儒科夫斯基翼型的升力系数和阻力系数随攻角的变化。 由图可见,攻角在10C -到10C 的范围内,理想流体导出的升力系数与实验符合得很好,这时没有发生严重的分离。 至于阻力的计算,则和圆柱绕流的情况一样,理想流体理论不能得出有用的结果。
图6-4儒科夫斯基翼型表面的压力分布在
流体理想与实际测量有相同的升力条件下 图6-5儒科夫斯基翼型的升力系数 理论值与实测值的比较 和阻力系数随攻角的变化
从上面两个例子可见,理想流体理论虽在某些方面(如圆柱体前缘附近的压力分布,翼型的压力分布和升力等)能得出与实际情况大体符合的结果,但不能用这种理论来预估阻力,它也不能处理不同雷诺数引起的差别以及分离等问题,而在许多工程技术问题中人们是很关心这些问题的。 因此需要研究有粘性的实际流体的运动和力的作用关系,即粘性流体的运动学和动力学。
6.2 粘性流体研究的特点(以不可压粘性流μ不变为例)
6.2.1 粘性流体有旋(只要壁面相对流场运动就是有旋运动)
理想流体运动一般为无旋运动,但也可作有旋运动。 根据亥姆霍兹定理,质量力有势的正压理想流体的涡量和环量具有守恒性,如果初始时刻或入口截面上运动是无旋的,则整个流场都是无旋的,反之则都有旋。 均匀流绕物体流动或物体在静止介质中运动时,从理想流动的观点来看,全流场都是无旋流动。 理想流体的有旋运动出现在质量力无势的斜压流体中,这类运动在气象学中会碰到。 与此相反,粘性流体运动除个别情况外,都是有旋运动,而且涡量和环量没有守恒性,在流动过程中,涡量不断生成,传输和衰减。 粘性运动的有旋性可通过实验观察到,也可从基本方程出发,从数学上得到证明。 下面从不可压缩流体的N-S 方程出发,用反证法来证明有旋性。
根据矢量分析和不可压缩流体的连续方程,可得
()()∆=∇∇⋅-∇⨯∇⨯=-∇⨯Ωv v v
()2
2
v ⋅∇=∇-⨯Ωv v v