yyyy-M-流体动力学基础流体动力学基础1/108第三章 流体动力学基础yyyy-M-流体动力学基础流体动力学基础2/108本章主要内容本章主要内容描述流体运动的方法描述流体运动的方法流场的若干概念流场的若干概念质点导数与系统导数质点导数与系统导数流体运动的基本物理定律及基本方程流体运动的基本物理定律及基本方程平行直线流断面上的压强关系式平行直线流断面上的压强关系式定常流动中的机械能关系定常流动中的机械能关系运动流体与固体壁面间的作用力运动流体与固体壁面间的作用力层流与湍流层流与湍流yyyy-M-流体动力学基础流体动力学基础3/108物质体方法和场方法物质体方法和场方法是从不同是从不同观点出发描述同一流体运动规观点出发描述同一流体运动规律的两种方法，律的两种方法，它们之间必然它们之间必然存在着相互的联系和某种对应存在着相互的联系和某种对应关系。关系。第三节第三节 质点导数与系统导数质点导数与系统导数yyyy-M-流体动力学基础流体动力学基础4/108t=1t=1时刻的场时刻的场yyyy-M-流体动力学基础流体动力学基础5/108t=2t=2时刻的场时刻的场yyyy-M-流体动力学基础流体动力学基础6/108t=2t=2时刻的场时刻的场yyyy-M-流体动力学基础流体动力学基础7/108对于在某一时刻（如对于在某一时刻（如t 时刻），时刻）， 占据着流场中某位置（如占据着流场中某位置（如 r 位位置）的流体质点的运动表现，置）的流体质点的运动表现，两种方法两种方法所描述的是同一个客所描述的是同一个客观存在观存在。第三节第三节 质点导数与系统导数质点导数与系统导数yyyy-M-流体动力学基础流体动力学基础8/108在经历在经历dt 时间后，原来占据时间后，原来占据r 位置的那个流体质点运动到达位置的那个流体质点运动到达r+dr 的新位置。的新位置。Euler法需要质点的坐标。法需要质点的坐标。Lagrange法不需要质点的坐标法不需要质点的坐标第三节第三节 质点导数与系统导数质点导数与系统导数yyyy-M-流体动力学基础流体动力学基础9/108牛顿第二定律的数学表达中，需要牛顿第二定律的数学表达中，需要用到用到同一个流体质点同一个流体质点（或同一个流（或同一个流体系统）的运动物理量对时间的变体系统）的运动物理量对时间的变化率，如何化率，如何追踪到我们研究的质点追踪到我们研究的质点的坐标的坐标是场方法能够应用牛顿第二是场方法能够应用牛顿第二定律的前提。定律的前提。yyyy-M-流体动力学基础流体动力学基础10/108一、质点导数一、质点导数yyyy-M-流体动力学基础流体动力学基础11/108质点导数质点导数流体质点具有的物理量对时间流体质点具有的物理量对时间的变化率，也称为质点的的变化率，也称为质点的随体导数。随体导数。任一物理量任一物理量 的导数的导数 yyyy-M-流体动力学基础流体动力学基础12/108注意到注意到 yyyy-M-流体动力学基础流体动力学基础13/108yyyy-M-流体动力学基础流体动力学基础14/108物理量物理量 的质点导数表达式的质点导数表达式 局部导数局部导数当地导数当地导数位变导数位变导数迁移导数迁移导数由物理量的非定常性产生的由物理量的非定常性产生的由场的不均匀性引起的由场的不均匀性引起的物理量物理量 的质点导数表达式的质点导数表达式 注意：注意：物理量物理量 可以可以是是标标量、矢量或量、矢量或张张量。量。 yyyy-M-流体动力学基础流体动力学基础15/108yyyy-M-流体动力学基础流体动力学基础16/108yyyy-M-流体动力学基础流体动力学基础17/108二、系统导数二、系统导数yyyy-M-流体动力学基础流体动力学基础18/108系统导数系统导数流体系统流体系统所具有的物理量所具有的物理量随时随时间的变化率。间的变化率。系统系统是一团流体质点的集合；在运动过是一团流体质点的集合；在运动过程中，系统始终包含着确定的流体质点，程中，系统始终包含着确定的流体质点，有确定的质量，而这一团流体的表面常常有确定的质量，而这一团流体的表面常常会不断地变形。会不断地变形。 yyyy-M-流体动力学基础流体动力学基础19/108控制体控制体流场中某一确定的空间区域，这个区流场中某一确定的空间区域，这个区域的域的周界称为控制面。周界称为控制面。虚线为系统虚线为系统 实线为控制体实线为控制体yyyy-M-流体动力学基础流体动力学基础20/108设设 是定是定义义在体在体积积内内单单位体位体积积所具有所具有的物理量，的物理量，yyyy-M-流体动力学基础流体动力学基础21/108则则体体积积V内内该该物理量的物理量的总总量量为为 。yyyy-M-流体动力学基础流体动力学基础22/108系系统统占有的空占有的空间间在不断改在不断改变变它的大小和形状它的大小和形状 也随也随时间时间而而变变。yyyy-M-流体动力学基础流体动力学基础23/108此此积积分分对时间对时间的的导导数数 即即为为系系统导统导数。数。yyyy-M-流体动力学基础流体动力学基础24/108根据导数的定义，有根据导数的定义，有（a）yyyy-M-流体动力学基础流体动力学基础25/108（b）yyyy-M-流体动力学基础流体动力学基础26/108（c）式（式（b b）右端第一项可表示为）右端第一项可表示为（b）yyyy-M-流体动力学基础流体动力学基础27/108式（式（b b）中第二项积分中的前一项可写为）中第二项积分中的前一项可写为（d）（b）yyyy-M-流体动力学基础流体动力学基础28/108式（式（b）中第二项积分中的后一项可表示为）中第二项积分中的后一项可表示为（e）流入流入（b）yyyy-M-流体动力学基础流体动力学基础29/108将式（将式（d）和式（）和式（e）合并，并注意到）合并，并注意到 ， 是整个控制体的面积，于是式（是整个控制体的面积，于是式（b）右端第）右端第二项成为二项成为 （f）yyyy-M-流体动力学基础流体动力学基础30/108将式（将式（c）和式（）和式（f）代入（）代入（b）中，得）中，得这就是流体系统内物理量对时间的随体导数公这就是流体系统内物理量对时间的随体导数公式式, ,也称为输运定理（或输运公式）也称为输运定理（或输运公式） yyyy-M-流体动力学基础流体动力学基础31/108上式即是上式即是按拉格朗日法求系按拉格朗日法求系统统内物理量的内物理量的时间变时间变化化率率转换为转换为按欧拉法按欧拉法计计算的公式。算的公式。系系统统内部物理量的内部物理量的时间变时间变化率由两部分化率由两部分组组成：成：控制体内控制体内 对对时间时间的的变变化率；化率；单单位位时间时间内内经过经过控制面控制面 的的净净通量。通量。yyyy-M-流体动力学基础流体动力学基础32/108第四第四节节 流体运流体运动动的基本的基本物理定律及基本方程物理定律及基本方程yyyy-M-流体动力学基础流体动力学基础33/108一、连续性方程（质量守恒）一、连续性方程（质量守恒）yyyy-M-流体动力学基础流体动力学基础34/108分析分析该微元六面体内该微元六面体内 质量质量的变化。的变化。速度分量为速度分量为u、v、wyyyy-M-流体动力学基础流体动力学基础35/108由表面由表面ABCD流出的流出的质质量是量是 yyyy-M-流体动力学基础流体动力学基础36/108在在dt 时间时间内在内在x方向方向净净流入的流入的质质量量为为 同理，同理，在在 y 方向方向和和 z 方向方向上上 dt 时间时间通通过过相相应应表面表面净净流入的流入的质质量分量分别为别为yyyy-M-流体动力学基础流体动力学基础37/108dt 时间时间通通过该过该微元六面体微元六面体所有方向所有方向净净流入的流入的质质量量为为 yyyy-M-流体动力学基础流体动力学基础38/108在在dt 时间内，因控制体内密度变化引起时间内，因控制体内密度变化引起的的总质量的增加量总质量的增加量为为 yyyy-M-流体动力学基础流体动力学基础39/108根据质量守恒定律得根据质量守恒定律得或或微分形式的连续性方程微分形式的连续性方程意义：表达了任何存在的流体流动必须满足的质量守恒条件意义：表达了任何存在的流体流动必须满足的质量守恒条件yyyy-M-流体动力学基础流体动力学基础40/108几种特殊情形下的几种特殊情形下的简简化形式化形式 （1 1）定常流动，非稳态项为零，只剩：）定常流动，非稳态项为零，只剩： （2）不可压缩流体，密度随着时间和空间均不变化）不可压缩流体，密度随着时间和空间均不变化yyyy-M-流体动力学基础流体动力学基础41/108在有限通道中作定常流在有限通道中作定常流动动的情况，的情况，质质量守恒可用量守恒可用简单简单的形的形式表达，而不必采用微分方程的形式。式表达，而不必采用微分方程的形式。流体做定常流流体做定常流动时动时，流流场场内的各物理量均不随内的各物理量均不随时间时间改改变变，则则1 断面断面上上单单位位时间时间内流入的流体内流入的流体质质量量应应等于等于2 断面断面流出流出的流体的流体质质量，量， 即即 沿流道的上下游两个过流断面沿流道的上下游两个过流断面若是不可压缩的均质流体，则连续性方程为：若是不可压缩的均质流体，则连续性方程为：yyyy-M-流体动力学基础流体动力学基础42/108二、运动微分方程二、运动微分方程（牛顿第二定律）（牛顿第二定律）yyyy-M-流体动力学基础流体动力学基础43/108流体的运流体的运动与其所受外力之与其所受外力之间的关系，的关系，应应遵循遵循动动量守恒定量守恒定律律牛牛顿顿第二定律。牛第二定律。牛顿顿第二定律的数学表达式即第二定律的数学表达式即为为运运动动方程。方程。1、理想流体的运、理想流体的运动动微分方程微分方程理想流体理想流体：流体的粘性被忽略不计流体。：流体的粘性被忽略不计流体。在运动的理想流体中，取出一个微元平行六面体的微团，在运动的理想流体中，取出一个微元平行六面体的微团，它的各边长度分别为它的各边长度分别为dx、dy和和dz，如图所示；，如图所示；理想流体运动时不产生内摩擦力，所以作用在流体微团上理想流体运动时不产生内摩擦力，所以作用在流体微团上的外力只有质量力和压强。的外力只有质量力和压强。微分方程的微分方程的导导出出yyyy-M-流体动力学基础流体动力学基础44/108微元六面体内理想流体微团微元六面体内理想流体微团的受力（的受力（y方向）方向） 假设假设六面体形心的坐标六面体形心的坐标为为x、y、z，速度、压强分别为速度、压强分别为u、v、w、p，单位质量力为单位质量力为 f . 先分析先分析 y 方向的运动方向的运动 左右两个平面中心点处的压左右两个平面中心点处的压强分别为强分别为及及作用在流体微团上的质量力在作用在流体微团上的质量力在 y 方向的分量方向的分量 沿沿y方向牛顿第二定律的表达式为方向牛顿第二定律的表达式为 yyyy-M-流体动力学基础流体动力学基础45/108将上式各将上式各项项除以流体微除以流体微团团的的质质量量 得到得到 y 方向的运动微分方程方向的运动微分方程类似地，在类似地，在 x 和和 z 两个方向上，分别有两个方向上，分别有理想流体的运动微分方程理想流体的运动微分方程yyyy-M-流体动力学基础流体动力学基础46/108在一般情况下，作用在流体上的质量力在一般情况下，作用在流体上的质量力fx、 fy 和和 fz 是已知是已知的，对理想不可压缩流体其密度的，对理想不可压缩流体其密度为一常数。在这种情况下，为一常数。在这种情况下，方程组中有四个未知数方程组中有四个未知数u、v、w和和p，而方程仅有三个，为，而方程仅有三个，为此需加上不可压缩流体的连续性方程，这样方程组封闭，此需加上不可压缩流体的连续性方程，这样方程组封闭，从理论上提供了求解的可能性。从理论上提供了求解的可能性。 方程组的封闭性问题方程组的封闭性问题yyyy-M-流体动力学基础流体动力学基础47/1082、实际流体的运动微分方程、实际流体的运动微分方程实际流体运动时，作用于微元六面体各个面上的应力不仅实际流体运动时，作用于微元六面体各个面上的应力不仅有法向应力，有法向应力，还有由于粘性产生的切向应力还有由于粘性产生的切向应力。yyyy-M-流体动力学基础流体动力学基础48/108根据牛顿第二定律，写出沿根据牛顿第二定律，写出沿y y方向的运动微分方程方向的运动微分方程 yyyy-M-流体动力学基础流体动力学基础49/108化简后，得化简后，得 y 方向的运动微分方程方向的运动微分方程同理可得同理可得以应力形式表示的粘性流体的运动微分方程以应力形式表示的粘性流体的运动微分方程yyyy-M-流体动力学基础流体动力学基础50/108方程组方程组不封闭不封闭九个应力三个速度分量均未知，12个未知数三个运动方程+1个连续性方程，4个方程不封闭。需要确定粘性流体中各应力与速度变形速率的关系yyyy-M-流体动力学基础流体动力学基础51/108对对于一般的流体，于一般的流体，各各应应力与速度力与速度变变形速率形速率间间的关系的关系为为： yyyy-M-流体动力学基础流体动力学基础52/108对对于一般的流体，于一般的流体，各各应应力与速度力与速度变变形速率形速率间间的关系的关系为为： yyyy-M-流体动力学基础流体动力学基础53/108将将应应力关系式，代入其力关系式，代入其应应力表示的运力表示的运动动微分微分方程中，方程中，化化简简得得不可不可压缩压缩粘性流体的运粘性流体的运动动微分方程，又称微分方程，又称纳维纳维-斯托克斯托克斯，斯，简称简称N-S方程方程yyyy-M-流体动力学基础流体动力学基础54/108方程组方程组封闭封闭在解决实际问题时，除上述方程外，还应包括连续性方程。四个方程四个未知数。在数学上却遇到很大困难，主要是由于纳维-斯托克斯方程是二阶的非线性偏微分方程，只有在少数情况下才可以求得方程的解析解。