空气动力学：第4章粘性流体动力学基础.ppt

EXIT4.1、流体的粘性及其对流动的影响、流体的粘性及其对流动的影响4.2、雷诺实验、层流与湍流、雷诺实验、层流与湍流4.3、粘性流体的应力状态、粘性流体的应力状态4.4、广义牛顿内摩擦定理（本构关系）、广义牛顿内摩擦定理（本构关系）4.5、粘性流体运动方程、粘性流体运动方程——Navier-Stokes方程方程4.6、流动相似及相似准则、流动相似及相似准则第第第第4 4章章章章 粘性流体动力学基础粘性流体动力学基础粘性流体动力学基础粘性流体动力学基础EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响1、流体的粘滞性、流体的粘滞性 在静止状态下，流体不能承受剪力但是在运在静止状态下，流体不能承受剪力但是在运动状态下，流体可以承受剪力，而且对于不同种流动状态下，流体可以承受剪力，而且对于不同种流体所承受剪力大小是不同的体所承受剪力大小是不同的 流体的粘滞性是指，流体在运动状态下抵抗剪流体的粘滞性是指，流体在运动状态下抵抗剪切变形能力切变形能力 流体的剪切变形是指流体质点之间出现相对运流体的剪切变形是指流体质点之间出现相对运动。

因此流体的粘滞性是指抵抗流体质点之间的相动因此流体的粘滞性是指抵抗流体质点之间的相对运动能力对运动能力EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 流体抵抗剪切变形能力，可通过流层之间的剪流体抵抗剪切变形能力，可通过流层之间的剪切力表现出来这个剪切力称为内摩擦力）流切力表现出来这个剪切力称为内摩擦力）流体在流动过程中，必然要克服内摩擦力做功，因此体在流动过程中，必然要克服内摩擦力做功，因此流体粘滞性是流体发生机械能损失的根源流体粘滞性是流体发生机械能损失的根源牛顿的内摩擦定律（牛顿的内摩擦定律（Newton，，1686年）年） F=µAU/hFhUEXIT  =µ/ ——流体的运动粘性系数流体的运动粘性系数量纲量纲[  ]=L2/T ，单位，单位m2/s 水：水： 1.139 10-6，空气：，空气：1.461 10-5 4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流层之间的内摩擦力与接触面上的压力无关。

流层之间的内摩擦力与接触面上的压力无关 设设 表示单位面积上的内摩擦力（粘性切应力），表示单位面积上的内摩擦力（粘性切应力），则则 µ——流体的动力粘性系数量纲流体的动力粘性系数量纲[µ]=M/L/T，单位，单位N·s/m2=Pa·s 水：1.139 10-3，空气：，空气：1.7894 10-5 EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响一般流层速度分布不是直线，如下图所示一般流层速度分布不是直线，如下图所示 u yyF=µAdu/dy  =µdu/dydu/dy —— 表示单位高度流层的速度增量，称为流表示单位高度流层的速度增量，称为流速梯度EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 速度梯度速度梯度du/dy物理上也表示流体质点剪切变形物理上也表示流体质点剪切变形速度，如图所示速度，如图所示 d  =dudt/dy d /dt=du/dy流体切应力与速度梯度的一般关系为流体切应力与速度梯度的一般关系为d dydudtu+duuEXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响1 —   = 0+µdu/dy2 —   =µ（（du/dy））0.53 —   =µdu/dy 4 —   =µ(du/dy)2 5 — 理想流体理想流体1---bingham流体，泥浆、血浆、牙膏等流体，泥浆、血浆、牙膏等2---伪塑性流体，尼龙、橡胶、油漆等伪塑性流体，尼龙、橡胶、油漆等3---牛顿流体，水、空气、汽油、酒精等牛顿流体，水、空气、汽油、酒精等4---胀塑性流体，生面团、浓淀粉糊等胀塑性流体，生面团、浓淀粉糊等5---理想流体，无粘流体。

理想流体，无粘流体 du/dy123450EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 自然界中流体都是有粘性的，因此粘性对流体自然界中流体都是有粘性的，因此粘性对流体运动的影响是普遍存在的但对于具体的流动问题，运动的影响是普遍存在的但对于具体的流动问题，粘性所起的作用并不一定相同特别是象水和空气粘性所起的作用并不一定相同特别是象水和空气这样的小粘性流体，对于某些问题忽略粘性的作用这样的小粘性流体，对于某些问题忽略粘性的作用可得到满意的结果因此为了简化起见，提出了理可得到满意的结果因此为了简化起见，提出了理想流体的概念和理论想流体的概念和理论2、粘性流体运动特点、粘性流体运动特点 以下用若干流动事例说明粘性流动与无粘流动以下用若干流动事例说明粘性流动与无粘流动的差别EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响（（1）绕过平板的均直流动）绕过平板的均直流动 当理想流体绕过平板（无厚度）时，平板对流当理想流体绕过平板（无厚度）时，平板对流动不产生任何影响，在平板表面，允许流体质点滑动不产生任何影响，在平板表面，允许流体质点滑过平板，但不允许穿透平板（通常称作为不穿透条过平板，但不允许穿透平板（通常称作为不穿透条件）。

平板对流动无阻滞作用，平板阻力为零平板对流动无阻滞作用，平板阻力为零EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 但如果是粘性流体，情况就不同了由于存在粘但如果是粘性流体，情况就不同了由于存在粘性，紧贴平板表面的流体质点粘附在平板上，与平板性，紧贴平板表面的流体质点粘附在平板上，与平板表面不存在相对运动（既不允许穿透，也不允许滑动）表面不存在相对运动（既不允许穿透，也不允许滑动），这就是说，在边界面上流体质点必须满足不穿透条，这就是说，在边界面上流体质点必须满足不穿透条件和不滑移条件件和不滑移条件 随着离开平板距离的增大，流体速度由壁面处的随着离开平板距离的增大，流体速度由壁面处的零值迅速增大到来流的速度这样在平板近区存在着零值迅速增大到来流的速度这样在平板近区存在着速度梯度很大的流动，因此流层之间的粘性切应力就速度梯度很大的流动，因此流层之间的粘性切应力就不能忽略，对流动起控制作用这个区称为边界层不能忽略，对流动起控制作用这个区称为边界层平板对流动起阻滞作用，平板的阻力不为零。

即平板对流动起阻滞作用，平板的阻力不为零即EXIT沿平板的边界层实验演示沿平板的边界层实验演示4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响EXIT无滑移实验演示无滑移实验演示与物面的粘附条件（无滑移条件）是粘性与物面的粘附条件（无滑移条件）是粘性流体运动有别与理想流体运动的主要标志流体运动有别与理想流体运动的主要标志4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响EXIT4.14.1流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响（（2）圆柱绕流）圆柱绕流 理想流体绕流圆柱时，在圆柱上存在前驻点理想流体绕流圆柱时，在圆柱上存在前驻点A，，后驻点后驻点D，最大速度点，最大速度点B、、C中心流线在前驻点分中心流线在前驻点分叉，后驻点汇合叉，后驻点汇合EXIT4.14.1流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 根据根据Bernoulli方程，在流体质点绕过圆柱的流方程，在流体质点绕过圆柱的流动过程中，在动过程中，在A—B（（C）区，流体质点在）区，流体质点在A点流速点流速为零，压强最大，以后质点的压强沿程减小，流速为零，压强最大，以后质点的压强沿程减小，流速沿程增大，到达沿程增大，到达B（（C）点流速最大，压强最小。

该）点流速最大，压强最小该区属于增速减压区，是顺压梯度区区属于增速减压区，是顺压梯度区EXIT4.14.1流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 在在B（（C））—D区，流体质点的压强沿程增大，区，流体质点的压强沿程增大，流速沿程减小，到达流速沿程减小，到达D点压强最大，流速为零点压强最大，流速为零该该区属于减速增压区，是逆压梯度区区属于减速增压区，是逆压梯度区EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 在流体质点绕过圆柱的过程中，只有动能、压在流体质点绕过圆柱的过程中，只有动能、压能的相互转换，而无机械能的损失在圆柱面上压能的相互转换，而无机械能的损失在圆柱面上压强分布对称，无阻力存在著名的达朗贝尔佯谬强分布对称，无阻力存在著名的达朗贝尔佯谬或达朗贝尔疑题）或达朗贝尔疑题）EXIT 对于粘性流体的绕流，与理想流体绕流存在很对于粘性流体的绕流，与理想流体绕流存在很大的差别由于流体与固壁表面的粘附作用，在物大的差别。

由于流体与固壁表面的粘附作用，在物面近区将产生边界层面近区将产生边界层4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响受流体粘性的阻滞作用，流体质点在由受流体粘性的阻滞作用，流体质点在由A点到点到B点的点的流程中，将消耗部分动能用之克服摩擦阻力做功，流程中，将消耗部分动能用之克服摩擦阻力做功，以至使其无法满足由以至使其无法满足由B点到点到D点压力升高的要求，压点压力升高的要求，压强分布与理想流体绕流不同强分布与理想流体绕流不同EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响随着雷诺数随着雷诺数 的不同，绕流压强分布也不同的不同，绕流压强分布也不同 但它们与理想流动的有一个共同区别在于粘性流动中但它们与理想流动的有一个共同区别在于粘性流动中圆柱体的背风面圆柱体的背风面Cp为负值。

为负值下图给出了不同雷诺数时圆柱表面的压强系数下图给出了不同雷诺数时圆柱表面的压强系数Cp分布EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响在下游点在下游点θθ=0=0°°处处Cp=+1.0 理想流体理想流体Cp=-1.1 层流，层流，Cp=-0.7 湍流，湍流，Cp=-0.1 湍流，湍流，EXIT 背风面的压强小于迎风面背风面的压强小于迎风面的压强，压强分布不对称的压强，压强分布不对称使圆柱体受到流体给它的使圆柱体受到流体给它的阻力，阻力的方向指向下阻力，阻力的方向指向下游4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 进一步的研究表明，圆柱进一步的研究表明，圆柱表面的粘性流体流动在靠近尾部表面的粘性流体流动在靠近尾部的某个区域将产生边界层分离，的某个区域将产生边界层分离，从而形成圆柱体下游一个由旋涡从而形成圆柱体下游一个由旋涡组成的尾流区。

尾流区的流动性组成的尾流区尾流区的流动性质随质随Re数的变化有所不同数的变化有所不同EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响 绕球流动的流场显示绕球流动的流场显示EXIT无逆压梯度的平板边界层沿流向会变厚但不分离，无逆压梯度的平板边界层沿流向会变厚但不分离，有摩擦阻力（平板边界层流动实验演示：）有摩擦阻力（平板边界层流动实验演示：）4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响EXIT有逆压梯度时，边界层变厚并可能分离有逆压梯度时，边界层变厚并可能分离（低速扩压段中的边界层与分离实验演示：）（低速扩压段中的边界层与分离实验演示：）4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响EXIT圆柱绕流有较大的逆压梯度，必然分离圆柱绕流有较大的逆压梯度，必然分离（圆柱绕流分离细节实验演示：）（圆柱绕流分离细节实验演示：）4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响EXIT4.1 4.1 流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响流体的粘性及其对流动的影响综合上述所述，结论如下：综合上述所述，结论如下： （（1）粘性摩擦切应力与物面的粘附条件（无滑移）粘性摩擦切应力与物面的粘附条件（无滑移条件）是粘性流体运动有别与理想流体运动的主要条件）是粘性流体运动有别与理想流体运动的主要标志。

标志2）粘性的存在是产生阻力的主要原因粘性的存在是产生阻力的主要原因 （（3）边界层的分离必要条件是，流体的粘性和逆）边界层的分离必要条件是，流体的粘性和逆压梯度 （（4）粘性对于研究阻力、边界层及其分离、旋涡）粘性对于研究阻力、边界层及其分离、旋涡的扩散等问题起主导作用，不能忽略的扩散等问题起主导作用，不能忽略EXIT 雷诺（雷诺（Osborne Reynolds，，1842~1921，英国工程师兼物理学，英国工程师兼物理学家，维多利亚大学教授）最早详细家，维多利亚大学教授）最早详细研究了管道中粘性流体的流动状态研究了管道中粘性流体的流动状态及其影响因素及其影响因素4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流 1880年，年，用管用管径径2.54cm、长度、长度1.372m玻璃管进玻璃管进行了著名的流态转捩试验，并于行了著名的流态转捩试验，并于1883年在一篇论文年在一篇论文中明确指出了管中水流存在层流和湍流（紊流）两中明确指出了管中水流存在层流和湍流（紊流）两种流态EXIT4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流 层流是指流体质点不互相混杂，流体质点作有条层流是指流体质点不互相混杂，流体质点作有条不紊的有序的运动。

不紊的有序的运动 紊流是一种杂乱无章、互相混掺，不规则的随机紊流是一种杂乱无章、互相混掺，不规则的随机运动EXIT4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流EXIT 实验发现，不同的流态对于流动的摩擦阻力、实验发现，不同的流态对于流动的摩擦阻力、压力损失、速度分布等影响很大压力损失、速度分布等影响很大4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流流态从层流到湍流的过渡称为转捩流态从层流到湍流的过渡称为转捩 实验表明流态的转捩不是单单取决于某一个流动实验表明流态的转捩不是单单取决于某一个流动参数参数V ,μ等，而是取决于无量纲的组合量等，而是取决于无量纲的组合量 Re这就是管道粘性运动的相似参数，雷诺数这就是管道粘性运动的相似参数，雷诺数EXIT 在非管道流动中，也存在这两种不同的流态：在非管道流动中，也存在这两种不同的流态：层流与湍流，从层流到湍流的转捩也与雷诺数大小层流与湍流，从层流到湍流的转捩也与雷诺数大小有关 雷诺数之所以对流态起着重要作用，从而对粘雷诺数之所以对流态起着重要作用，从而对粘性流体运动的其他特性起着重要作用，在于雷诺数性流体运动的其他特性起着重要作用，在于雷诺数具有很强的物理意义。

具有很强的物理意义4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流 雷诺数的物理意义：雷诺数代表作用在流体微雷诺数的物理意义：雷诺数代表作用在流体微团上的惯性力与粘性力之比团上的惯性力与粘性力之比EXIT惯性力正比于质量乘加速度惯性力正比于质量乘加速度其中，质量正比于密度乘尺度的其中，质量正比于密度乘尺度的 3 次方次方: ～～ρρL L3 3加速度正比于速度加速度正比于速度 2 2 次方和尺度－次方和尺度－1 1次方：～次方：～V V2 2L L-1-1从而惯性力正比于：从而惯性力正比于： ～～ρV2L2 粘性力正比于剪应力乘面积粘性力正比于剪应力乘面积其中，剪应力正比于粘性系数乘速度其中，剪应力正比于粘性系数乘速度1次方和尺度的次方和尺度的－－1次方次方: ～～μVL-1 -1 ，， 面积正比于尺度面积正比于尺度2 2次方：～次方：～L L2 2从而粘性力正比于：从而粘性力正比于： ～～μVL因此惯性力与粘性力之比正比于：因此惯性力与粘性力之比正比于：ρVL/μ , 此即雷此即雷诺数4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流EXIT4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流 惯性力的作用是促使质点失稳，扰动放大；粘性惯性力的作用是促使质点失稳，扰动放大；粘性力的作用是对质点起约束作用的，是遏制扰动的。

力的作用是对质点起约束作用的，是遏制扰动的 流动为紊流时，流动为紊流时，Re数较大，此时质点惯性力大数较大，此时质点惯性力大于粘性力，流动失去稳定于粘性力，流动失去稳定 流动为层流时，流动为层流时，Re数较小，此时质点粘性力大数较小，此时质点粘性力大于惯性力，流动稳定，层次分明于惯性力，流动稳定，层次分明EXIT4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流 在管道紊流中，由于流体质点的随机脉动，致在管道紊流中，由于流体质点的随机脉动，致使流层之间的动量发生不断的交换，快层流体速度使流层之间的动量发生不断的交换，快层流体速度减慢，慢层流体速度增大，造成时均流速分布更加减慢，慢层流体速度增大，造成时均流速分布更加均匀EXIT 层流与层流与 湍流的对比湍流的对比管道雷诺实验管道雷诺实验平板边界层：上：平板边界层：上： 湍流湍流 下：层流下：层流4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流EXIT 层流与湍流的区别层流与湍流的区别 层流层流 湍流湍流4.2 4.2 雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流雷诺实验、层流与湍流1. 外观外观色线规则，流动分色线规则，流动分层，外表光滑层，外表光滑流动紊乱不规则，外表流动紊乱不规则，外表粗糙粗糙2. Re较小较小较大较大3. 质量与质量与动量交换动量交换层间只限于分子层间只限于分子间的较小的扩散间的较小的扩散宏观微团纵、横向大宏观微团纵、横向大的质量、动量交换的质量、动量交换4. 速度分速度分布布较尖瘦的抛物线分较尖瘦的抛物线分布，壁面附近速度布，壁面附近速度和梯度都相对较小和梯度都相对较小较饱满的对数分布，壁面较饱满的对数分布，壁面附近速度和梯度相对较大附近速度和梯度相对较大5. 损失损失随随Re增加而降低增加而降低随随Re增加转捩时损失增加增加转捩时损失增加6. 剪应力剪应力牛顿应力及雷诺应力牛顿应力及雷诺应力牛顿应力牛顿应力EXIT4.3 4.3 粘性流体的应力状态粘性流体的应力状态粘性流体的应力状态粘性流体的应力状态1、理想流体和粘性流体作用面受力差别、理想流体和粘性流体作用面受力差别 流体处于静止状态，只能承受压力，几乎不能承受拉力流体处于静止状态，只能承受压力，几乎不能承受拉力和剪力，不具有抵抗剪切变形的能力。

理想流体在运动状态和剪力，不具有抵抗剪切变形的能力理想流体在运动状态下流体质点之间可以存在相对运动，但不具有抵抗剪切变形下流体质点之间可以存在相对运动，但不具有抵抗剪切变形的能力因此，作用于流体内部任意面上的力只有正（法）的能力因此，作用于流体内部任意面上的力只有正（法）向力，无切向力向力，无切向力 粘性流体在运动状态下，流体质点之间可以存在相对运粘性流体在运动状态下，流体质点之间可以存在相对运动，流体具有抵抗剪切变形的能力因此，作用于流体内部动，流体具有抵抗剪切变形的能力因此，作用于流体内部任意面上力既有正向力，也有切向力任意面上力既有正向力，也有切向力EXIT 在粘性流体运动中，由于存在切向力，过任意一点单位在粘性流体运动中，由于存在切向力，过任意一点单位面积上的表面力就不一定垂直于作用面，且各个方向的大小面积上的表面力就不一定垂直于作用面，且各个方向的大小也不一定相等因此，作用于任意方向微元面积上合应力可也不一定相等因此，作用于任意方向微元面积上合应力可分解为法向应力和切向应力分解为法向应力和切向应力4.3 粘性流体的应力状态粘性流体的应力状态 如果作用面的法线方向与坐标轴重合，则合应力可分解如果作用面的法线方向与坐标轴重合，则合应力可分解为三个分量，其中垂直于作用面的为法应力，另外两个与作为三个分量，其中垂直于作用面的为法应力，另外两个与作用面相切为切应力，分别平行于另外两个坐标轴，为切应力用面相切为切应力，分别平行于另外两个坐标轴，为切应力在坐标轴向的投影分量。

在坐标轴向的投影分量2、粘性流体中的应力状态、粘性流体中的应力状态 EXIT4.3 粘性流体的应力状态粘性流体的应力状态 由此可见，用两个下标可把各个应力分量的作用面方位由此可见，用两个下标可把各个应力分量的作用面方位和投影方向表示清楚其中第一个下标表示作用面的法线方和投影方向表示清楚其中第一个下标表示作用面的法线方向，第二个下标表示应力分量的投影方向向，第二个下标表示应力分量的投影方向如，对于如，对于x面的合应力可表示为面的合应力可表示为y面的合应力表达式为面的合应力表达式为z面的合应力表达式为面的合应力表达式为EXIT 如果在同一点上给定三个相互垂直坐标面上的应力，那么如果在同一点上给定三个相互垂直坐标面上的应力，那么过该点任意方向作用面上的应力可通过坐标变换唯一确定过该点任意方向作用面上的应力可通过坐标变换唯一确定上面这个结论可利用对微元六面体的动量矩定理得到证明上面这个结论可利用对微元六面体的动量矩定理得到证明4.3 粘性流体的应力状态粘性流体的应力状态 我们把三个坐标面上的九个应力分量称为该点的应力状我们把三个坐标面上的九个应力分量称为该点的应力状态，由这九个应力分量组成的矩阵称为应力矩阵（或应力张态，由这九个应力分量组成的矩阵称为应力矩阵（或应力张量）。

量）这九个应力分量并不全部独力，其中的六个切向应力是两两这九个应力分量并不全部独力，其中的六个切向应力是两两相等的，所以独力的一共是三个法向的，三个切向的相等的，所以独力的一共是三个法向的，三个切向的EXIT（（1）在理想流体中，不存在切应力，三个法向应力相等，）在理想流体中，不存在切应力，三个法向应力相等，等于该点压强的负值即等于该点压强的负值即（（2）在粘性流体中，任意一点的任何三个相互垂直面上的）在粘性流体中，任意一点的任何三个相互垂直面上的法向应力之和为一个不变量，并定义此不变量的平均值为法向应力之和为一个不变量，并定义此不变量的平均值为该点的平均压强的负值即该点的平均压强的负值即（（3）在粘性流体中，任意面上的切应力一般不为零在粘性流体中，任意面上的切应力一般不为零4.3 粘性流体的应力状态粘性流体的应力状态EXIT4.4 4.4 广义牛顿内摩擦定理（本构关系）广义牛顿内摩擦定理（本构关系）广义牛顿内摩擦定理（本构关系）广义牛顿内摩擦定理（本构关系）Stokes（（1845年）根据牛顿内摩擦定理的启发，（粘性流体作直年）根据牛顿内摩擦定理的启发，（粘性流体作直线层状流动时，流层之间的切应力与速度梯度成正比）并在做了线层状流动时，流层之间的切应力与速度梯度成正比）并在做了一些合理的假设之后，一些合理的假设之后， 将牛顿内摩擦定律进行推广，提出广义将牛顿内摩擦定律进行推广，提出广义牛顿内摩擦定理牛顿内摩擦定理——应力应变率关系（或称本构关系）：应力应变率关系（或称本构关系）： 这个关系将六个应力与微团的变形率直接联系（线性关系）。

满这个关系将六个应力与微团的变形率直接联系（线性关系）满足上述关系的流体称为牛顿流体足上述关系的流体称为牛顿流体 EXIT对于不可压缩流体，上述应力应变率关系可化简为：对于不可压缩流体，上述应力应变率关系可化简为：4.4 广义牛顿内摩擦定理（本构关系）广义牛顿内摩擦定理（本构关系）EXIT4.5 4.5 粘性流体运动方程粘性流体运动方程粘性流体运动方程粘性流体运动方程————NavierNavier-Stokes-Stokes方程方程方程方程1、流体运动的基本方程、流体运动的基本方程 利用牛顿第二定理推导以应力形式表示的流体运动微分利用牛顿第二定理推导以应力形式表示的流体运动微分方程向推导欧拉方程一样，在流场中取一个微元六面体方程向推导欧拉方程一样，在流场中取一个微元六面体进行分析，以进行分析，以x方向为例，建立运动方程方向为例，建立运动方程作用在作用在ABCD和和A’B’C’D’两个侧面的两个侧面的法向力差是：法向力差是：作用在作用在ABB’A’和和CDC’D’两个侧面的两个侧面的切向力差是：切向力差是：EXIT作用在作用在ADA’D’和和BCB’C’两个侧面的切向力差是：两个侧面的切向力差是：仍然设单位质量彻体力分量为：仍然设单位质量彻体力分量为：fx , fy , fz ,按照牛顿第按照牛顿第二定律：二定律：是欧拉法表示的加速度或速度的物质导数。

是欧拉法表示的加速度或速度的物质导数4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT或：或：同理：同理：将反映粘性应力与变形率关系的广义牛顿内摩擦定理代入将反映粘性应力与变形率关系的广义牛顿内摩擦定理代入上式右端，即得到粘性流动的运动方程上式右端，即得到粘性流动的运动方程 N－－S 方程：方程：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT其中其中 是拉普拉斯算子：是拉普拉斯算子：可见，对于理想流右端的粘性项为零，方程化为欧拉方程可见，对于理想流右端的粘性项为零，方程化为欧拉方程4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程写成向量形式写成向量形式 为为EXIT当当不可压时，根据连续方程：不可压时，根据连续方程：则不可压粘流的则不可压粘流的 N－－S方程写为：方程写为：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT用用 三个方向的单位向量三个方向的单位向量 i 、、j、、k 分别乘上三式并相加，分别乘上三式并相加，可得不可压粘流可得不可压粘流 N-S方程比较简捷的向量形式：方程比较简捷的向量形式：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程其中其中为速度分量为速度分量为哈密顿算子为哈密顿算子为拉普拉斯算子为拉普拉斯算子为单位质量的质量力为单位质量的质量力EXIT 与第二章一样，这个方程中速度的随体导数可以加以分与第二章一样，这个方程中速度的随体导数可以加以分解，把涡量分离出来，写成格罗米柯形式的方程也称为解，把涡量分离出来，写成格罗米柯形式的方程也称为Lamb型方程。

这样有利于研究流体的有旋性：型方程这样有利于研究流体的有旋性：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程2、、 Bernoulli积分积分 与与Bernoulli积分理想流体运动方程类似，积分积分理想流体运动方程类似，积分N-S方方程假定程假定:（（1）不可压缩粘性流体；（）不可压缩粘性流体；（2）定常流动；（）定常流动；（3））质量力有势；（质量力有势；（4）沿流线积分沿流线积分 沿流线积分沿流线积分N-S方程，可推导出粘性流体的能量方程方程，可推导出粘性流体的能量方程与理想流体能量不同的是，方程中多了一项因粘性引起的与理想流体能量不同的是，方程中多了一项因粘性引起的损失项，表示流体质点克服粘性应力所消耗的能量损失项，表示流体质点克服粘性应力所消耗的能量EXIT4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程 在粘性不可压缩定常流动中，任取一条流线，在流线在粘性不可压缩定常流动中，任取一条流线，在流线上某处取一微段上某处取一微段ds，该处所对应的流速为，该处所对应的流速为 对不可压对不可压N-S方程的三个分量分别乘方程的三个分量分别乘 dx、、dy、、dz后相加得：后相加得：EXIT定常情况下，上式左端即：定常情况下，上式左端即：注意到沿流线有流线方程：注意到沿流线有流线方程：前式右端第一项可由流线方程进行变量替换，得：前式右端第一项可由流线方程进行变量替换，得：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT彻体力有势彻体力有势ΠΠ，因此有：，因此有：不可压缩流动，有：不可压缩流动，有：粘性项写成为粘性项写成为第二、三项同理可替换为：第二、三项同理可替换为：三项之和为：三项之和为：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT 与理想流体能量微分方程相比，在上式中多了一项与粘与理想流体能量微分方程相比，在上式中多了一项与粘性有关的项，物理上表示单位质量流体质点克服粘性应力所性有关的项，物理上表示单位质量流体质点克服粘性应力所做的功，代表机械能的损失，不可能再被流体质点机械运动做的功，代表机械能的损失，不可能再被流体质点机械运动所利用。

故称其为单位质量流体的机械能损失或能量损失故称其为单位质量流体的机械能损失或能量损失 对于质量力只有重力的情况，方程的形式变为对于质量力只有重力的情况，方程的形式变为 方程两边同除以方程两边同除以g，，得到得到表示单位重量流体总机械能量沿流线的变化表示单位重量流体总机械能量沿流线的变化从而不可压从而不可压N-S方程，在定常沿流线条件下可写为：方程，在定常沿流线条件下可写为：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT 如果令如果令 能量方程变为能量方程变为 单位重量流体所具有的机械能为单位重量流体所具有的机械能为 ；单位重量流体粘；单位重量流体粘 性力所做的功为性力所做的功为 沿着同一条流线积分，得到：沿着同一条流线积分，得到：4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT 4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程 上式说明，在粘性流体中，沿同一条流线上单上式说明，在粘性流体中，沿同一条流线上单位重量流体的所具有的机械能总是沿程减小的，不位重量流体的所具有的机械能总是沿程减小的，不能保持守恒（理想流体时，总机械能是保持守恒的，能保持守恒（理想流体时，总机械能是保持守恒的，无机械能损失），减小的部分代表流体质点克服粘无机械能损失），减小的部分代表流体质点克服粘性应力做功所消耗的机械能量。

粘性流体的性应力做功所消耗的机械能量粘性流体的Bernoulli积分方程说明，粘性流体在流动中，无论积分方程说明，粘性流体在流动中，无论势能、压能和动能如何转化，但总机械能是沿程减势能、压能和动能如何转化，但总机械能是沿程减小的，总是从机械能高的地方流向机械能低的地方小的，总是从机械能高的地方流向机械能低的地方EXITN－－S方程为非线性偏微分方程，它的求解一般需要借助计方程为非线性偏微分方程，它的求解一般需要借助计算机用数值方法求解而在一些简单的粘流问题上，算机用数值方法求解而在一些简单的粘流问题上，N－－S方程也有解析解方程也有解析解例：求解二维平行壁之间的不可压粘性流动，二壁固定例：求解二维平行壁之间的不可压粘性流动，二壁固定2bxy解解: 设流动定常，彻体力可略设流动定常，彻体力可略二维不可压二维不可压 N N––S S 方程写为：方程写为：3. N-S方程的解析解举例方程的解析解举例*4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT由于由于 ，第二个方程化为，第二个方程化为：：即在流动横截面压强不变。

又第一个方程化为即在流动横截面压强不变又第一个方程化为：：对对 y 积分，注意到积分，注意到 不是不是 y 的函数，对的函数，对 y 积分时当常数看积分时当常数看4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT由边界条件定常数由边界条件定常数 C1 和和 C2 ：：y=±b 处，处，u=0，定得，定得 C1＝＝0，，C2＝－＝－b2/2，于是：，于是：即即 u 在在y 向作抛物线分布中心点流速为：向作抛物线分布中心点流速为：表明沿表明沿x轴轴 是个负值，即压强是逐步下降的一段长是个负值，即压强是逐步下降的一段长度度 L 上的压降是：上的压降是：这个压降是用于克服壁面摩擦阻力的这个压降是用于克服壁面摩擦阻力的4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT管道平均流速为：管道平均流速为：壁面摩擦应力为：壁面摩擦应力为：一段长一段长 L 的壁面上摩擦应力是：的壁面上摩擦应力是：两侧壁面上的总摩擦力是两侧壁面上的总摩擦力是这个力刚好等于压降乘以通道面积，说明流动的损失完全消耗这个力刚好等于压降乘以通道面积，说明流动的损失完全消耗在克服壁面摩擦上了。

在克服壁面摩擦上了4.5 粘性流体运动方程粘性流体运动方程——Navier-Stokes方程方程EXIT4.6 流动相似及相似准则流动相似及相似准则* 有了有了N－－S方程，我们可以看一看考虑粘性的作用在内，方程，我们可以看一看考虑粘性的作用在内，流动要符合什么条件才能相似要流动相似，毫无疑问首流动要符合什么条件才能相似要流动相似，毫无疑问首先要求对流动起扰动作用的物体形状要相似，即尺寸成比先要求对流动起扰动作用的物体形状要相似，即尺寸成比例假设有两个流场，一个记为例假设有两个流场，一个记为1，另一个记为，另一个记为2，两个流，两个流动的各项参数之间的比例关系为：动的各项参数之间的比例关系为：其中的其中的 r 是比例系数，下标代表了是比例系数，下标代表了相应的参数相应的参数12EXIT 将上述比例关系代入第二个流场的将上述比例关系代入第二个流场的x方向方向N－－S方程：方程：可得：可得：4.6 流动相似及相似准则流动相似及相似准则*EXIT方程的各不同类型的项都出来一个各比例常数所组成的数，方程的各不同类型的项都出来一个各比例常数所组成的数，如果这些数相等，那么流场如果这些数相等，那么流场 2 的方程就变成流场的方程就变成流场 1 的方程的方程了。

方程一样，无量纲的边界条件也一样，流动就相似了方程一样，无量纲的边界条件也一样，流动就相似了所以流场相似要求：所以流场相似要求：这个式子包含四个独力的等式关系这个式子包含四个独力的等式关系第一、二两部分相等得：第一、二两部分相等得：即：即：4.6 流动相似及相似准则流动相似及相似准则*EXIT第二、三两部分相等得：第二、三两部分相等得：第二、四两部分相等得：第二、四两部分相等得：或：或：第二、五两部分相等得：第二、五两部分相等得：4.6 流动相似及相似准则流动相似及相似准则*EXIT上面四个等式规定了四个相似参数，或称相似准则这些相上面四个等式规定了四个相似参数，或称相似准则这些相似参数都是无量纲的数似参数都是无量纲的数 称为斯托罗哈尔数，是非定常运动的相似准则称为斯托罗哈尔数，是非定常运动的相似准则称为傅劳德数，是水面船只运动须遵循的相似准则称为傅劳德数，是水面船只运动须遵循的相似准则称为马赫数，是高速可压缩流动须遵循的相似准则称为马赫数，是高速可压缩流动须遵循的相似准则称为雷诺数，是粘性流动须遵循的相似准则称为雷诺数，是粘性流动须遵循的相似准则。

4.6 流动相似及相似准则流动相似及相似准则*EXIT当两个流场的相似参数完全相等时，两个流场相似当两个流场的相似参数完全相等时，两个流场相似事实上当考虑到其他一些力学和热力学关系时还会有相应的事实上当考虑到其他一些力学和热力学关系时还会有相应的相似参数要求满足然而要想同时满足这些相似准则几乎是相似参数要求满足然而要想同时满足这些相似准则几乎是不可能的在实用中只能就具体问题抓住主要矛盾，满足一不可能的在实用中只能就具体问题抓住主要矛盾，满足一个（顶多两个）相似准则个（顶多两个）相似准则雷诺数与粘性流动的流态有密切关系，从而影响到跟粘性有雷诺数与粘性流动的流态有密切关系，从而影响到跟粘性有关的物理量如摩擦阻力等，因此雷诺数是研究粘性流体运动关的物理量如摩擦阻力等，因此雷诺数是研究粘性流体运动时须尽量满足的相似准则时须尽量满足的相似准则即使是雷诺数一个相似准则要准确满足也是不容易的，例如即使是雷诺数一个相似准则要准确满足也是不容易的，例如在风洞中作飞机模型实验，模型尺寸远小于实际飞机，要满在风洞中作飞机模型实验，模型尺寸远小于实际飞机，要满足雷诺数相似需要加大实验速度足雷诺数相似需要加大实验速度 v，增大，增大ρρ 或减小或减小 μ。

一来增大一来增大 v 是受是受动力所限的，二来力所限的，二来 v 过高高时M数也不数也不对了4.6 流动相似及相似准则流动相似及相似准则*EXIT本章基本要求本章基本要求1.了解流体的粘性及其对流动的影响了解流体的粘性及其对流动的影响2.了解雷诺实验、掌握雷诺数的定义与意义、层流与湍了解雷诺实验、掌握雷诺数的定义与意义、层流与湍流的流的 特征与区别特征与区别3.了解粘性流体的应力状态了解粘性流体的应力状态4.了解广义牛顿内摩擦定理（本构关系）了解广义牛顿内摩擦定理（本构关系）5. 了解粘性流体运动方程了解粘性流体运动方程---N-S方程，掌握方程，掌握N-S各项所各项所代表的意义代表的意义。