第九章 粘性不可压缩流体运动-1

1,第九章粘性不可压缩流体运动,2,第一节粘性不可压缩流体的运动特点,研究粘性不可压缩流体运动的特点：产生内磨擦及传热有关的能量耗损过程,主要解决的重点：流动过程阻力的产生机理、计算与控制；,考虑粘性作用的情况有：,(1)流体运动能量损耗的过程,(2)粘性力与惯性力同阶或较惯性力大得多的时候（如在边界层内）,3,粘性不可压缩均质流体运动方程组,连续性方程,运动方程,能量方程,本构方程,状态方程,4,粘性不可压缩均质流体运动方程组,连续性方程,NS方程,本构方程,涡旋运动方程,5,初始条件与边界条件,(1)初始条件：t=0时，流场中已知速度分布及压力分布,(2)边界条件：,静止固壁上：满足粘附条件,运动固壁上：满足,自由面上：满足,6,(1)有旋性：绝大部分粘性不可压缩流体运动都是有旋的,粘性流体运动的一般性质,(2)涡旋的扩散性：涡旋强的地方将向涡旋弱的地方输运涡量，直至涡量相等为止。,7,粘性流体运动的一般性质,(3)机械能的损耗性：由于粘性的存在，体力和面力所做的功只有一部分转化为动能，另一部分被粘性应力耗损变成了热能，单位体积内耗损的动能参耗损函数确定：,8,第二节粘性不可压缩流体运动方程组的求解途径,方程组的特点：二阶非线性偏微分方程组,连续性方程,NS方程,本构方程,解的存在和唯一性？,9,第二节粘性不可压缩流体运动方程组的求解,层流,湍流,准确解,粘性不可压缩流体运动,近似解,小Re数,大Re数,中Re数,统计理论,模式理论,混合长度理论,K-方程,RSM模型,10,9.2.1粘性不可压缩流体层流运动的准确解,粘性不可压缩流体在无限长柱形管道内的定常运动,已知：管截面上的形状及两个截面上的压力,求：管截面上速度分布、流量及管道中的阻力系数,11,连续性方程,一维流动,12,NS方程,13,NS方程,连续性方程,14,边界条件,静止固壁上：满足粘附条件,在截面a处，即x=0，满足：,在截面b处，即x=l，满足：,15,16,P为常数,在截面a处，即x=0，满足：,在截面b处，即x=l，满足：,17,压力沿轴向线性下降,泊松方程,18,二元二阶偏微分方程,(1)轴对称流动：圆心在原点的圆管中粘性流体运动,(2)平面运动：两个平行x-z坐标面的无限长平面间的粘性流体运动,一元二阶常微分方程,准确解,边界条件,直接积分,19,二元二阶偏微分方程,(1)轴对称流动：圆心在原点的圆管中粘性流体运动,结构轴对称,流动分布轴对称,20,(1)轴对称流动：圆心在原点的圆管中粘性流体运动,21,(1)轴对称流动：圆心在原点的圆管中粘性流体运动,在边壁,在中心,22,(1)轴对称流动：圆心在原点的圆管中粘性流体运动,粘性不可压缩流体运动轴对称圆管内定常层流,23,(a)速度分布剖面,24,(b)流量及平均速度,r,半径r处圆环的面积,25,(c)阻力系数,r=a时：,定义阻力系数：,26,(3)阻力系数,定义雷诺系数：,在轴对称圆管内定常层流状态下与实验值吻合,27,(2)两个平行板间的定常运动及库塔流,由上、下两个平行平板组成的二维渠道，粘性不可压缩流体在压差作用下在渠道内作定常流动，板间距离为2h,28,在上截面处，即y=h，满足：,在下截面处，即y=-h，满足：,29,二维泊素叶流动,30,纯剪切流动,在上截面处，即y=h，满足：,在下截面处，即y=0，满足：,U,纯剪切流动：上、下游没有压差，只有平板的拖动,31,库塔流：既有压差，又有平板的拖动,在上截面处，即y=h，满足：,在下截面处，即y=0，满足：,U,Pa,Pb,32,粘性不可压缩流体层流运动的准确解小结,NS方程,一元二阶常微分方程,33,9.2.2粘性不可压缩流体层流运动的近似解,普朗特边界层方程大Re数下层流运动近似解,惯性力远远大于粘性力,能否忽略粘性力的作用？,34,普朗特边界层方程大Re数下层流运动近似解,惯性力远远大于粘性力,如忽略粘性力的作用，简化N-S方程,与理想不可压缩流体运动方程相同,边界上不满足粘附条件：v=0,35,边界层(附面层)：当流体流过物体，或物体在流体中运动时，在物体表面和与之直接接触的薄层流体之间，由于粘性的存在，都会出现附着作用，而使这一层流体附着在物体表面，速度为零，与相邻的另一层流体之间便出现速度梯度。离开表面向外沿法线方向延伸，速度急剧增大，速度梯度则逐渐减小。速度梯度变化很大的那一层流体称为边界层或附面层。,边界层,外流区,36,整个绕流区,边界层,外流区,外流区：忽略粘性力的作用理想无旋,边界层：考虑粘性力的作用粘性有旋,边界线：与来流速度相差1%的流体质点连线,37,普朗特的观点：,外流区：粘性力远远小惯性力的作用，忽略粘性力的作用理想无旋（平面势流）,边界层：粘性力与惯性力同量级，考虑粘性力的作用粘性有旋，边界层厚度比特征长度L小得多，而且x方向速度分量沿法线方向的变化比切向大得多。（NS方程）,38,边界层内粘性不可压缩流体基本方程（二维）,39,边界层内粘性不可压缩流体基本方程（二维）,(1)通过量纲分析，发现比更高阶的无穷小，故可忽略此方向的压力变化，即：,即压力数值穿过边界层并不改变，同时忽略y方向的动量方程,(2)通过量纲分析，发现：,是更高阶的无穷小，可忽略,40,边界层内粘性不可压缩流体基本方程（二维）,41,边界条件,静止固壁上：满足粘附条件,在边界层边界y=处，满足：,U(x)是边界层外部边界上外流的速度分布,42,初始条件：,t=t0时刻，已知全部区域内的速度及压力分布,43,绕流区域内粘性不可压缩流体基本方程（二维）普朗特边界层方程,外部区域内理想流体的运动方程,44,绕流区域内粘性不可压缩流体基本方程（二维）普朗特边界层方程的核心思想,1提出了边界层的概念，合理的将整个绕流区划分为两个部分；,2在边界层内合理的将NS方程进行了简化；,3整个绕流区内压力的分布不受边界层分布的影响，与理想不可压缩流体无旋运动时相同,45,9.2.3半无穷长平板的层流边界层普朗特边界层方程的Blasius解,无限空间中一均匀气充以速度U沿板面方向定常地向一半无穷长且厚度为零的平板流来，求解在板面上边界层内的速度分布,1边界层以外的区域：速度场均匀定常且为常数U,46,2边界层内流体的运动方程,47,边界条件,静止固壁上：y=0,在边界层边界y=处，满足：,48,2边界层内流体的运动方程,引入流函数,二元二阶非线性偏微分方程组,连续性方程自动满足,49,2边界层内流体的运动方程,边界条件,静止固壁上：y=0,在边界层边界y=处，满足：,一元三阶非线性偏微分方程,50,2边界层内流体的运动方程,根据量纲分析，构造组合变量,使得：,51,2边界层内流体的运动方程,52,2边界层内流体的运动方程,边界条件,一元三阶非线性常微分方程,53,Blasiuse的解决方案：,普朗特边界层方程,一元三阶非线性常微分方程,二元二阶非线性偏微分方程组,一元三阶非线性偏微分方程组,54,Blasiuse的解层流边界层近似解,55,9.3边界层脱体现象及产生的条件,降压增速区,增压减速区,顺压区,逆压区,56,9.3边界层脱体现象及产生的条件,57,9.3边界层脱体现象及产生的条件,顺压区,逆压区,58,9.3边界层脱体现象及产生的条件,59,9.3边界层分离现象及产生的条件,结论1：沿流动方向存在逆压区是产生流动脱体现象的原因,60,粘性流体的绕流存在逆压区，但不产生流动分离现象,61,沿流动方向存在逆压区，但不产生流动分离现象,理想不可压缩流体无旋运动圆柱无环量绕流,62,9.3边界层分离现象及产生的条件,结论2：沿流动方向存在逆压区和壁附近的粘性滞止作用是产生流动脱体现象的必要条件,结论3：沿流动方向存在逆压区和壁附近的粘性滞止作用，并且逆压梯度足够大时，才产生流动脱体现象,63,