空气动力学-第4章资料.ppt

本章基本要求 了解流体的粘性及其对流动的影响 了解雷诺实验、掌握雷诺数的定义与意义、层流与湍流的 特征与区别 了解粘性流体的应力状态 了解广义牛顿内摩擦定理（本构关系） 5. 了解粘性流体运动方程---N-S方程，掌握N-S方程各项所代表的意义 日本设计的机械蜻蜓 俄罗斯设计的机械蜻蜓 美国设计的 机械苍蝇 微型飞行器的飞行雷诺数只有几百到几万的量级（102～104） 空气中的悬浮尘埃其运动雷诺数则更低甚至可以小于1 需要再次强调：雷诺数代表惯性力与粘性力之比只是宏观量级上的比例关系，根据雷诺数的大小可以判断流动中何种因素占主导作用，但绝不能认为 Re=1 表示流动的惯性力与粘性力刚好相等。 雷诺数对流动影响的进一步举例 斯托克斯阻力定律： 高雷诺数时物体受到的流动阻力正比于: ρV2 L2 低雷诺数时物体受到的流动阻力正比于: μV L 1、理想流体和粘性流体作用面受力差别 4.3、粘性流体的应力状态 静止或理想流体内部任意面上只有法向力，无切向力 粘性流体内部任意面上力既有正向力，也有切向力 在粘性流体运动中，过任意一点任意方向单位面积上的表面力不一定垂直于作用面，可分解为法向应力和切向应力 如果作用面的法线方向与坐标轴重合，则合应力可分解为三个分量，分别为法应力分量和切应力分量 4.3、粘性流体的应力状态 2、粘性流体中的应力状态 从而三个面的合应力可表示为 x面 : y面: z面: 4.3、粘性流体的应力状态 由此可见，用两个下标可把各个应力分量的作用面方位和投影方向表示清楚。其中第一个下标表示作用面的法线方向，第二个下标表示应力分量的投影方向。 如果在同一点上给定三个相互垂直坐标面上的应力，那么过该点任意方向作用面上的应力可通过坐标变换唯一确定。 上述九个应力分量可写为： 这九个应力分量并不全部独力，其中的六个切向应力是两两相等的，所以独立的一共是三个法向的，三个切向的。 这个结论可利用对微元六面体的动量矩定理得到证明，思路是：一对剪应力对微元产生的力矩将与彻体力力矩和微元质量的动量矩平衡，而后二者都正比于微元的体积乘以微距离，是一个高阶小量可略去，从而得到这一对剪应力相等。 4.3、粘性流体的应力状态 注：有的教材将法向应力记为: 关于应力的几个要点： （1）在理想流体及静止流体中不存在切应力，三个法向应力相等（各向同性），等于该点压强的负值。即： （2）在粘性运动流体中，任意一点的任何三个相互垂直面上的法向应力之和为一个不变量，并定义此不变量的平均值为该点的平均压强的负值。即： （3）在粘性运动流体中，任意面上的切应力一般不为零。 4.3、粘性流体的应力状态 4.4、广义牛顿内摩擦定理（本构关系） Stokes（1845年）根据牛顿内摩擦定理的启发（粘性流体作直线层状流动时，层间切应力与速度梯度成正比），在一些合理的假设下将牛顿内摩擦定律进行推广，提出广义牛顿内摩擦定理----应力应变率关系（本构关系）： 这个关系将六个应力与微团的变形率直接联系（线性关系）。满足上述关系的流体称为牛顿流体。 对于不可压缩流体，上述应力应变率关系可化简为： 4.4、广义牛顿内摩擦定理（本构关系） 本构关系满足： 4.5、粘性流体运动方程---Navier-Stokes方程 1、流体运动的基本方程 利用牛顿第二定理推导以应力形式表示的流体运动微分方程。像推导欧拉方程一样，在流场中取一个微元六面体进行分析，以x方向为例，建立运动方程。现在由于是粘性流体，作用在中心P点处不仅有法向应力，而且还有切向应力，控制面上的应力可用中心点处应力泰勒召开表示。 作用在ABCD和A’B’C’D’两个侧面的法向力差是： 作用在ABB’A’和CDC’D’两个侧面的切向力差是： 作用在ADA’D’和BCB’C’两个侧面的切向力差是： 仍然设单位质量彻体力分量为：fx , fy , fz , 按照牛顿第二定律： 是欧拉法表示的加速度或速度的物质导数。 4.5、粘性流体运动方程---Navier-Stokes方程 或： 同理： 将反映粘性应力与应变率关系的广义牛顿内摩擦定理代入上式右端，即得到粘性流动的运动方程 N－S 方程： （纳维Navier, C. L. M. H. 1785-1836, 法国力学家、工程师；斯托克斯Stokes, G. G. 1819-1903, 英国力学家、数学家） 4.5、粘性流体运动方程---Navier-Stokes方程 其中 是拉普拉斯算子： 可见，对于理想流右端的粘性项为零，方程化为欧拉方程。 4.5、粘性流体运动方程---Navier-Stokes方程 当不可压时，