流体流动微分方程课件.ppt

1、第章流体流动微分方程 流体流动微分方程是一组微分方程，包括连续性方程和运动方程。连续性方程是流动流体质量守恒的数学描述。与第4章中基于控制体建立的质量守恒方积方程相对应，连续性方程是基于流场中的点（微元体）所建立的质量守恒微分方程。运动方程则是流动流体动量守恒的数学描述。与第4章中基于控制体建立的动量守恒方程相对应，运动方程是基于流场中的点（微元体）所建立的动量守恒微分方程，又称为运动微分方程。通过第5章中对典型以为流动问题的分析，已经了解了将动量守恒定理应用于流场微元体从而建立运动微分方程的基本方法和过程。本章将把这一基本方法推广应用于三维情况，建立一般条件下的流体运动微分方程。就其目的而言

2、，积分方程反映流动过程中流体总质量、总动量和总能量的变化，而本章要建立的流动微分方程，目的在于流场分布的详细信息，以揭示宏观流动现象的内在规律。第六章 流体流动微分方程.连续性方程.以应力表示的运动方程.粘性流体运动微分方程.流体流动微分方程的应用6.1.1 直角坐标系中的连续性方程连续性方程反映流动过程遵循质量守恒这一事实。对于流场中的微元体，质量守恒原理可以类似于控制体仿照式（4-9）表述为0质量变化微元体内的量流量的微元体输量流量的输出微元体率质入质（6-1）为了获得（6-1）的数学表达式连续性方程，不妨对图-所示的微元体进行分析。该微元体取自流场中的任意点，微元体在x、y、z方向的边长

3、分别为dx、dy、dz，其六个面两两相互平行且分别垂直于x、y、z。流体在点的密度为，速度为，其x、y、z方向的分量分别为 。一般而言，速度和密度均为坐标x、y、z和时间的函数。zyx、已知，流体穿越某一表面时的质量流量等于质量通量与表面积的乘积，而质量通量则为流体密度与流体在该表面上的法向速度的乘积，因此考察微元体上的输入与输出，首先要确定微元面上的法向速度法向速度。如图-所示,对于在流场中任意点所取的微元体，因为与点相邻的三个微元面上的流体或流动参数反映的是点的参数，所以在这三个微元面上，流体密度均为（点密度），且每一个面上流体的三个速度分量都为 （点速度）。其中，对于dydx微元面，因其

4、与轴垂直，该微元面上的三个速度分量中，是法向速度，其产生的法向通量为 ，而另为两个速度分量 则平行于dydz平面，与质量输出输入无关（故图中dydz微元面上的质量通量 ）。zyx、xx zy、x 同理，在垂直于y、z方向的微元面dxdz和dxdy上，法向速度分别为 ，质量通量分别为 ，如图6-1所示。按速度与坐标方向一致为正的约定，、都是输入通量。于是将这三个通量分别乘以相应的面积dydz、dxdz、dxdy后相加，可得输入微元体的质量流量为 相应地，当流体从与点不相邻的、分别垂直于、方向的三个微元面上流出时，由于分别经过dx、dy、dz的距离后其输出时的质量通量将发生变化，如图-所示所以输出

5、微元体的质量流量为zy、zy、yx、zdxdydxdzdydzzyxdxdydzzdxdzdyydydzdxxzzyyxx为了找到这种关系，斯托克斯（Stokes）提出了三种假设。为了找到这种关系，斯托克斯（Stokes）提出了三种假设。流体正应力与线变形率相关这一性质，与虎克定律中固体正应力与其线应变相关是类似的。也正是这一特点，对于不可压缩流体，无论是为稳态流动还是非稳态流动，其连续性方程都是一样的。动量通量表示单位时间、单位面积输入输出的动量，单位为 。但与固体应力不一样的是，流体的应力不是与应变的大小而是与应变的速率（即单位间内的应变）直接相关的。在工程实际中，有时采用柱坐标和球坐标描

6、述问题比采用直角坐标更为方便，比如，对于常见的圆管内的流动，最适用的显然是柱坐标系统。，则相应有 。流体流动微分方程的应用对于湍流流动，一般认为非稳态的-方程对湍流的瞬时运动仍然是适用的，但湍流的瞬时运动具有高度的随机性，要追踪这种随机运动是十分困难的。将上述条件代入柱坐标系的连续性方程（6-8）和运动微分方程（6-33）可得同理，分别考察y方向和z方向的动量在微元体表面上的输出与输入可得在此处主要指微元体表面上受到应力。将上述条件代入柱坐标系的连续性方程（6-8）和运动微分方程（6-33）可得本构方程用于流体应力的分析与计算比如，对于附加粘性正应力由此并根据上述关于应力下标和方向的规定，不难

7、标出微元体个表面上和方向的正应力和切应力，如图-所示。根据速度分布及柱坐标下的牛顿本构方程（6-34）得切应力分布为（1）柱坐标系中的-方程和牛顿流体本构方程 由上述两项可得 两一方面，对于图所示的微元体，其瞬间质量为dxdydz，所以dxdydzz)z(y)y(x)x(量流量的微元体输量流量的输出微元体质入质（-）dxdydz t 质量变化微元体内的率（6-3）连续性方程连续性方程 将式（6-2）和式（6-3）代入微元体质量守恒文字表达式（6-1）可得直角坐标系中的连续性方程为或以矢量简洁表示为0 tz)(y)(x)(zyx（6-4a）0 t)(（6-4b）其中，是质量通量的的散度，是矢量算

8、子。)(kzjyix 由于导出方程（6-4）的过程中没有对流体和流动状态作任何假设，故该方程对层流和湍流、牛顿流体和非牛顿流体均适用。将方程（6-4a）展开并引用第2章中的随体导数（质点导数）概念，可将连续性方程表示为另一种形式其中，是速度矢量 的散度；是密度随体导数，按第2章中随体导数的定义有0z y x tDDzyx或0)(tDD（6-5a）（6-5b）tD/Dzyx t tDDzyx由于间隙较小，且流体黏度较高，故可将流动视为沿切向的一维层流流动。对于以为径向坐标、为周向坐标、为轴向坐标的柱坐标体系见图62（a），其连续性方程为运动方程是基于流场中的点（微元体）所建立的动量守恒方程，又称

9、为运动微分方程。对于以为径向坐标、为周向坐标、为轴向坐标的柱坐标体系见图62（a），其连续性方程为引入新变量后，x方向运动可变换为下述方程所以，流场中任一点的9个应力分量中，只有6个分量是独立的。这正是第章中分析不可压缩流体一维流动时曾经用到的条件。不可压缩流体的连续性方程解 参照图6-7的柱坐标系，根据不可压缩流体一维稳态层流的条件及流动与z无关的特点有在物理意义上，速度的散度表示为单位体积的流量在单位时间内的体积增量，通常称为体变形率对于不可压缩流体，不管其体积形状如何变化，其体积的大小不会变，故体变形率为零，即 。引入新变量后，x方向运动可变换为下述方程需要指出：任何流体的连续运动，都必

10、须首先满足相应的连续性方程，所以连续性方程是流体流动微分方程最基本的方程之一。流体流动微分方程的应用微元体表面力的总力分量按速度与坐标方向一致为正的约定，、都是输入通量。因此，流动微分方程的应用求解，关键是根据问题特点对一般形式的运动方程进行简化，获得针对具体问题的微分方程或方程组，并同时提出相关的初始条件和边界条件。初始条件是非稳态问题所要求的，因为与时间相关的问题必须以某一时刻的流动条件（即初始条件）为参照；于是，有上述方向动量的输出流量与输入流量相减得到连续性方程是流动流体质量守恒的数学描述。首先，有r和方向的运动方程可知，p*只能是z的函数，所以在轴向位置确定的管道截面上必然有对于x方

11、向的运动方程，因为t=0时刻所有流体是静止的，所以方程初始条件可表达为：对于0,有本构方程用于流体应力的分析与计算不可压缩流体的连续性方程 对于不可压缩流体，因密度=const，所以连续性方程简化为0z y x zyx或0（6-6a）（6-6b）在物理意义上，速度的散度表示为单位体积的流量在单位时间内的体积增量，通常称为体变形率体变形率对于不可压缩流体，不管其体积形状如何变化，其体积的大小不会变，故体变形率为零，即 。也正是这一特点，对于不可压缩流体，无论是为稳态流动还是非稳态流动，其连续性方程都是一样的。0 不可压缩流体的连续性方程不仅形式简单，而且应用广泛，因为工程实际中除了经常遇到不可压

12、缩流体外，不少可压缩流体的流动亦可常密度流体处理。由连续性方程(6-6)可知，对于不可压缩流体沿方向的一维流动，其连续性方程就是 。这正是第章中分析不可压缩流体一维流动时曾经用到的条件。0zy0 x/x6.1.2 柱坐标和球坐标系中的连续性方程柱坐标和球坐标系中的连续性方程 在工程实际中，除了直角坐标系外，出于描述的方便还经常采用柱坐标（如圆管流动问题）和球坐标（如球体绕流问题）。在此不加推导地写出这两种坐标系中的连续性方程，以供使用。图62对于以为径向坐标、为周向坐标、为轴向坐标的柱坐标体系见图62（a），其连续性方程为0)(z )(r1)r (r r1 tzr（6-7）其中，分别为r、z坐

13、标方向的速度分量。特别地，对于不可压缩流体，柱坐标系下的连续性方程可简化为zr、0z r1r )r (r1zr（6-8）在球坐标体系中，若以为径向坐标、为周向坐标、为经向坐标，见图6-2（b），则连续性方程为 需要指出：任何流体的连续运动，都必须首先满足相应的连续性方程，所以连续性方程是流体流动微分方程最基本的方程之一。0)(rsin1)sin(rsin1)r (r r1 tr22（6-9）其中，分别为、坐标方向的速度分量。r、.以应力表示的运动方程运动方程是基于流场中的点（微元体）所建立的动量守恒方程，又称为运动微分方程。所谓以应力表示的运动微分方程就是直接根据动量守恒定律得到的含流体应力的

14、微分方程，这就相当于第章分析一维流动问题时所得到的关于切应力的微分方程。针对微元体应用动量守恒原理时，由于微元体是在确定的空间点来考察流体流动时所取得一个体积为dxdydz的流场空间，具有类似于控制体的性质，因而其动量守恒原理可仿照式（4-10）表述为动量变化率动量流量动量流量微元体内的的输入微元体的输出微元体诸力之矢量和作用于微元体 在三维流动条件下，该方程各项的数学表达式远比一维时复杂，因此将分小节分别讨论。（6-10）.作用于微元体上的力 体积力是由于外力场（如重力场、离心力场、电磁场等）的作用在微元体整个体积上所产生的力又称为彻体力。由于体积力的大小与流体的质量成正比，故又称质量力。如

15、6-3所示，若微元体中单位质量流体的体积力在x、y、z方向的分量分别为 ，则 按作用区域的不同，作用于微元体上的力分为体积力和表面力两类。（1）体积力体积力zyxf、f、f 特别地，如果流体只受重力场作用（通常情况如此），且重力加速度的方向与轴正方向相反，则有 。可以，重力场条件下很容易确定单位体积力或其分量。）11c-6（dddfz）11b-6（dddfy）11a-6（dddfxzyxzzyxyzyxx方向的质量力微元体在方向的质量力微元体在方向的质量力微元体在-=f、0=f、0=fzyx（2）表面力表面力 表面力就是作用于流体表面的力。在此处主要指微元体表面上受到应力。如图（6-3）所示，

16、在微元体任何一个表面上，不管总应力的方向如何，总可以按坐标方向将其分解成一个正应力（或称法向应力）和两个切应力。对于图6-3中的微元体，在邻近点并分别垂直于x、y、z方向的三个微元面上，正应力和切应力分别为的微元面作用于与轴方向垂直的微元面作用于与轴方向垂直的微元面作用于与轴方向垂直、zyzxzzyzyxyyxzxyxx应力下坐标的意义应力下坐标的意义 每个应力都有两个下标，第一个下标表示应力作用面的法线方向，第二个下标表示应力的作用方向。例如，表示垂直于轴的表面上沿方向作用的切应力，表示垂直于轴的表面上沿方向作用的切应力，则表示垂直于轴表面正应力，等等。yxyzyy 关于应力的正负，通常规定

17、：若应力所在平面的外法线与坐标轴正向一致，则指向坐标轴正向的应力为正，反之为负；若应力所在平面的外法线与坐标轴正向相反，则指向坐标轴负向的应力为正，反之为负。可以参见图（-），其中所示的正应力和切应力均为正方向。对于正向力（法向应力），这种规定与“拉应力为正，压应力为负”的约定是一致的。应力正负的规定应力正负的规定方向运动方程的形式常称为扩散方程，做适当的变量代换可将方程转化为常微分方程，从而获得方程的解。因此取微元体分析一维不可压缩稳态层流问题时，微元面上的法向力可直接以压力标注，第五章中也正是这样做的。(其中，)其中谈到在等速流动中的管道截面上不存在表面粘性力的做功问题，正是基于等速流动时

18、 )。由此可见，附加粘性正应力的产生是速度沿流动方向的变化所导致的。运动方程是基于流场中的点（微元体）所建立的动量守恒方程，又称为运动微分方程。动量通量表示单位时间、单位面积输入输出的动量，单位为 。对于以应力表示的运动方程，要建立补充方程首先应该寻求运动方程中的未知量即流体应力与速度变化之间的内在联系。若应力所在平面的外法线与坐标轴正向相反，则指向坐标轴负向的应力为正，反之为负。相应地，当流体从与点不相邻的、分别垂直于、方向的三个微元面上流出时，由于分别经过dx、dy、dz的距离后其输出时的质量通量将发生变化，如图-所示所以输出微元体的质量流量为-方程 对于以r为径向坐标、为周向坐标、z为轴

19、向坐标的柱坐标体系见图6-2(a)，其粘性流体运动微分方程在r、z方向的分量式为（=const，=const）1节已讨论过，这里不再赘述。运动方程则是流动流体动量守恒的数学描述。其中，分别为r、坐标方向的速度分量；1节中讨论过，此处不再赘述。如图6-6所示，由于只受到重力场的作用，所以管道截面上任意一点处r、z方向单位质量流体的体积力为需要指出：任何流体的连续运动，都必须首先满足相应的连续性方程，所以连续性方程是流体流动微分方程最基本的方程之一。不过，对于工程实际问题由于总有特殊性使方程得到简化，从而有可能获得准确或近似的分析解。其中，g（rsin）是以管道中心线水平面为参照、相对高度为（rs

20、in）的单位体积流量的重力位能。解 参照图6-7的柱坐标系，根据不可压缩流体一维稳态层流的条件及流动与z无关的特点有动量流量及动量变化率需要指出：任何流体的连续运动，都必须首先满足相应的连续性方程，所以连续性方程是流体流动微分方程最基本的方程之一。应力状态及切应力互等定理 上述点处三个微元面上的个应力，代表了流场中某一点的应力状态，也就是说，粘性流场中任意一点的应力有个分量，包括个正应力分量和个切应力分量。可以证明，在个切应力分量中，互换下标的每一对切应力相等的，即切应力互等定理zyyzzxxzyxxy（6-12）所以，流场中任一点的9个应力分量中，只有6个分量是独立的。微元体表面力的总力分量

21、 为简明起见，从方向视图来观察微元体个表面上和方向的应力分量，如图所示。若以点相邻表面上的应力为基准，则与点不相邻的表面上的应力将产生一个随距离变化的增量。例如在处且垂直于方向的微元面上的正应力为，则在距离dz的平行面上的正应力就为。由此并根据上述关于应力下标和方向的规定，不难标出微元体个表面上和方向的正应力和切应力，如图-所示。于是，将各表面上方向的应力与相应的微元面积的dydz、dxdz或dxdy相乘，然后将轴正方向的各表面力与轴负方向的各表面力相减可得zzdz)z/(zzzz 13c)-(6 dxdydz zyx13b)-(6 dxdydz zyx 13a)-(6 dxdydz zyxz

22、zyzxzzyyyxyzxyxxx微元体方向的表面力微元体方向的表面力微元体方向的表面力.动量流量及动量变化率动量流量及动量变化率动量通量与动量流量动量通量与动量流量 已知，动量流量质量流量流体速度；类似地则有，动量通量质量通量流体速度。动量通量表示单位时间、单位面积输入输出的动量，单位为 。在确定了动量通量以后，就可将其与流通面积相乘得到动量流量，即“动量流量动量通量流通面积”，这与“质量流量质量通量流通面积”是类似的。)sm/(kg2 例如，在图-所示的微元体dxdz面上时，就会同时带入x、y、z方向的动量，且根据上述动量通量的定义x、y、z方向动量在该微元面上的输入通量就分别为、,而x、

23、y、z方向动量在该微元面上的输入流量则分别为xyzyyy、dxdz.、dxdz、dxdzzyyyxy输入输出微元体的动量流量输入输出微元体的动量流量 现以图-所示的微元体为对象，考察流体方向动量在微元体表面的输入输出。由图-已知，在点处且分别垂直于、方向的三个微元面上，流体进入微元体时的质量通量分别为 。由于这三个微元面上都有x方向的分速度 ，所以，根据，“动量通量=质量通量流体速度”的定义，流体x方向动量在这三个微元面上的输入通量就分别为 （注：图6-5中标注的是动量输入或输出方向，而动量或其通量本身的方向均指向方向，即分速度 的方向）。将这三个动量通量乘以相应的面积后相加，可得微元体上方向

24、动量的输入流量为zyx、x、xyxxxzx 相应的，当流体从与点不相邻的三个微元面上流出时，考虑到动量通量的变化，（如图-所示），可得微元体上方向动量的输出流量为dxdydxdzdydzxzxyxxdxdydzzdxdzdyydydzdxxxzxzxyxyxxxx 于是，有上述方向动量的输出流量与输入流量相减得到dxdydz z)(y)(x)(z的输入流量方向动量微元体表面的输出流量方向动量微元体表面（6-14a）同理，分别考察y方向和z方向的动量在微元体表面上的输出与输入可得dxdydz z)(y)(x)(z2yx的输入流量方向动量微元体表面的输出流量方向动量微元体表面dxdydz z)(y

25、)(x)(2zzzx的输入流量方向动量微元体表面的输出流量方向动量微元体表面（6-14b）（6-14c）微元体内的动量变化率微元体内的动量变化率dxdydz、dxdydz、dxdydzzyx 在图-5所示的微元体内，流体的瞬时质量为dxdydz，所以微元体内流体在x、y、z方向的瞬时动量分别为 ，于是有)15c-6(dxdydzt)15b-6(dxdydzt)15a-6(dxdydztzyx变化率微元体内方向的动量变化率微元体内方向的动量变化率微元体内方向的动量.以应力表示的运动方程前面已经导出微元体动量守恒式（6-10）中各项文字的数学表达式。由于动量守恒式对各坐标方向均成立，将分方向把有关

26、各项的数学表达式代入式（6-10），从而建立、方向的运动方程。首先，将方向的体积力式（6-11a）、表面力式（6-13a）、动量流量式（6-14a）、动量变化率式（6-15a）代入微元体动量守恒表达式（6-10），可得x方向运动方程的初步形式为zyxf t z)(y)(x)(zxyxxxxxz（6-16）方程（6-16）左边展开后可表达为zy x tzy x t tz)(y)(x)(xzxxxxzxxz根据6.1节的连续性方程（6-4）可知0zy x tzx所以x方向的运动方程简化为zyxfzy x tzxyxxxxxzxxxx（6-17a）同理可得y、z方向的运动方程分别为zyxfzy x

27、tzzyzxzzzzzzxzzyxfzy x tzyyyxyyyzyyxy（6-17b）（6-17c）式（式（6-17）就是以应力表示的粘性流体的运）就是以应力表示的粘性流体的运动方程。动方程。无论是牛顿流体还是非牛顿流体、是层流流体还是湍流流体，该方程均适用。对于x方向的运动方程，因为t=0时刻所有流体是静止的，所以方程初始条件可表达为：对于0,有相应的，当流体从与点不相邻的三个微元面上流出时，考虑到动量通量的变化，（如图-所示），可得微元体上方向动量的输出流量为为了找到这种关系，斯托克斯（Stokes）提出了三种假设。其中，分别为r、坐标方向的速度分量；试写出该条件下的连续性方程和运动微分

28、方程，并证明管道截面上任意一点的总势能p/+g（rsin）和轴向压力梯度 为常数。如图6-6所示，由于只受到重力场的作用，所以管道截面上任意一点处r、z方向单位质量流体的体积力为1节已讨论过，这里不再赘述。类似地则有，动量通量质量通量流体速度。由连续性方程可知，速度 仅是和的函数，即 。和切应力做功的功率。以应力表示的运动方程需要补充方程才能求解，与第章分析一维流动问题时需要引入补充方程才能由切应力方程得到速度微分方程有些相似。(其中，)解 参照图6-7的柱坐标系，根据不可压缩流体一维稳态层流的条件及流动与z无关的特点有解 取坐标与重力加速度方向平行，方向与正方向相反，即x、y、z方向单位质量

29、流量的体积力分别为fx=0、fy=0、fz=-，又因伯努利方程的前提条件是理想不可压缩流体的稳态压缩，即=、=const、，所以根据这些条件简化方程（6-31）可得至于简化后得到的运动方程，有的可能求不出解析解，或许只能得到近似解或通过数值计算方法获得离散解（见第12章）。由此并根据上述关于应力下标和方向的规定，不难标出微元体个表面上和方向的正应力和切应力，如图-所示。解 参照图6-7的柱坐标系，根据不可压缩流体一维稳态层流的条件及流动与z无关的特点有对于正向力（法向应力），这种规定与“拉应力为正，压应力为负”的约定是一致的。牛顿流体的本构方程例6-4 沿流线的伯努利方程但要考虑流体参数和变化

30、的情况，应将有关物性变化的关系作为补充方程。9时，过渡雷诺数Re又将增大到70000。方程的物理意义方程的物理意义 以方程（6-17a）为例，由第2章可知，方程左边的括号项是分速度x的随体导数 ，即任意时刻考察点的流体质点的加速度分量 ；由式（6-11）和式（6-13）可知，方程右边分别表示作用与单位流体体积上的表面力和体积力的分量，其合力作x表示。很明显，该方程可以简略地表示为：，这就是以单位体积的流体质量为基准的牛顿第二定律 在方向的分量式。相应的，和方向的运动方程表示的就是牛顿第二定律在和方向上的分量式。值得指出的是，在方程组（6-17）中，即使将密度和体积力看成是已知的，方程中仍然有个

31、未知量：个速度分量和个独立的应力分量，但该方程组加上连续性方程只有个方程，所以方程组是不封闭的。因此，要求解这个方程组，尚需要能将未知量关联起来的补充方程。/DtDxxxFaaxa.粘性流体运动微分方程以应力表示的运动方程需要补充方程才能求解，与第章分析一维流动问题时需要引入补充方程才能由切应力方程得到速度微分方程有些相似。在一维流动分析中，所引入的补充方程是牛顿剪切定律。在本章中也类似，所要引入的是广义的牛顿剪切定律牛顿流体本构方程牛顿流体本构方程。本节的目的就是引入牛顿流体本构方程，将应力从运动方程（6-17）中消去，得到由速度分量和压力表示的粘性流体运动微分方程耐维-斯托克斯（Navie

32、r-Stokes）方程。.牛顿流体的本构方程（1）基本假设基本假设 对于以应力表示的运动方程，要建立补充方程首先应该寻求运动方程中的未知量即流体应力与速度变化之间的内在联系。流体之所以流动是因为受到剪切应力的作用，同时，由于粘性的存在，流体对剪切应力要产生抵抗，这中抵抗以应力的形式表现出来，这与固体受到应力时要产生的应力是类似的。但与固体应力不一样的是，流体的应力不是与应变的大小而是与应变的速率（即单位间内的应变）直接相关的。x/xx)/2/y/(yx 流体力学中，称单位时间的应变为变形速率，包括形变速率如 、角变形率如 和体变形率 等（见第7章）。因此，建立补充方程的关键归结为寻求一般情况下

33、流体应力与变形速率之间的关系。为了找到这种关系，斯托克斯（Stokes）提出了三种假设。应力与变形速率呈线性关系。该假设得到牛顿剪切定律的启示，既然以为流动中 与变形速率 呈线性关系，于是可设想一般情况下也有这样的关系。应力与变形速率的关系各向同性。该假设认为，既然常见流体的物理性质都是各向同性的，于是可以设想应力与变形速率的关系也具有各向同性的性质。yx/dydx 静止流场中，切应力为零，各正向应力均等于静压力，即xxv=0yyv=0zzv=0-。该假设是根据静止流体不能承受切应力，而流动流体又不能承受拉应力而做出的。（2）牛顿流体的本构方程牛顿流体的本构方程 在上述假设条件下，既可推导出一

34、般情况下流体应力与变形速率之间的关系。在此略过复杂的推导过程，直接给出这一关系牛顿流体的本构方程牛顿流体的本构方程18)-(6 xzyzxyzyx32z2zyx32y2zyx32x2zxxzzxzyzyyzyxyxxyzyxzzzzyxyyyzyxxxx（3）本构方程的讨论本构方程的讨论 牛顿流体的本构方程除了在其建立流动微分方程中所具有的重要价值外、对本构方程本身的解析亦可增进对流动过程中流体变形速率、应力、压力等有关概念的理解。该假设得到牛顿剪切定律的启示，既然以为流动中 与变形速率 呈线性关系，于是可设想一般情况下也有这样的关系。因此，建立补充方程的关键归结为寻求一般情况下流体应力与变形

35、速率之间的关系。如图6-6所示，由于只受到重力场的作用，所以管道截面上任意一点处r、z方向单位质量流体的体积力为和切应力做功的功率。如果将 到 0.如果在N-S方程中所有速度项为零，即得到流体静力学方程。流体正应力与线变形率相关这一性质，与虎克定律中固体正应力与其线应变相关是类似的。为此，将不加推导地写出这两种坐标系下常密度和常粘度流体的运动微分方程和牛顿流体本构方程。流体流动微分方程的应用例6-1 圆管内的一维稳态流动分析第六章 流体流动微分方程(2)不可压缩流体的N-S方程解 参照图6-6的柱坐标体系，根据不可压缩流体一维稳态层流的条件及圆管的对称性有很显然，从流体流动的角度看，线变形率

36、的正负反映了流体沿方向的流动是加速还是减速，若就其目的而言，积分方程反映流动过程中流体总质量、总动量和总能量的变化，而本章要建立的流动微分方程，目的在于流场分布的详细信息，以揭示宏观流动现象的内在规律。动量流量及动量变化率其中，是质量通量的的散度，是矢量算子。-方程 对于以r为径向坐标、为周向坐标、z为轴向坐标的柱坐标体系见图6-2(a)，其粘性流体运动微分方程在r、z方向的分量式为（=const，=const）按作用区域的不同，作用于微元体上的力分为体积力和表面力两类。如图-所示,对于在流场中任意点所取的微元体，因为与点相邻的三个微元面上的流体或流动参数反映的是点的参数，所以在这三个微元面上

37、，流体密度均为（点密度），且每一个面上流体的三个速度分量都为 （点速度）。为了获得（6-1）的数学表达式连续性方程，不妨对图-所示的微元体进行分析。这正是第章中分析不可压缩流体一维流动时曾经用到的条件。正应力与线变形率正应力与线变形率 有本构方程可见，流体正应力与三个速度偏导数有关，即 ，他们分别是x、y、z方向的速度沿自身方向的变化率，其意义是单位时间内流体在x、y、z方向的线应变，称为线变形率线变形率；这三个线变形率之和即速度散度 则表示单位时间内流体的体积应变，称为体体变形率变形率。流体正应力与线变形率相关这一性质，与虎克定律中固体正应力与其线应变相关是类似的。z/y/x/x、线变形率与

38、流体流动线变形率与流体流动 很显然，从流体流动的角度看，线变形率 的正负反映了流体沿方向的流动是加速还是减速，若 则意味着流体在方向做等速流动或没有流动，对于 道理也一样；而体变形率 的正负则反映了流动过程中流体体积是增加还是减少，若 则意味着流动过程中流体体积不变，不可压缩流体的流动就属这种情况。z/y/、x/x0 x/x0 正应力中的粘性应力正应力中的粘性应力xx 有本构方程可见，流体正应力可视为由两部分构成：一部分是流体压力，另一部分则是流体变形速率所产生的附加粘性正应力。以 为例，如果用 表示其附加粘性正应力。可表示为xxxxxxxxp-其中)(32x2xxx 为了说明附加粘性正应力的

39、产生和意义，不妨考察流体只沿方向流动的情况。此时，所以x/0 xzy，x3 4xxx由此可见，附加粘性正应力的产生是速度沿流动方向的变化所导致的。加速时 ，所以 ；减速时，所以 。物理意义上，因为加速度同方向一前一后两流体质点将处于分离趋势，流体线的变形为拉伸变形，故由此产生的附加粘性正应力为拉应力；反之，减速时同方向一前一后两流体质点将处于挤压趋势，流体线的变形为压缩变形，故由此产生的附加粘性正应力为压应力。0 x/x0 xx0 x/x0 xx 特别地，如果该流动是等速的，即 则必然有 (注：在4.5节中所说的流体克服粘性力做功的功率 ，指的就是流体克服控制面上附加粘性正应力 和切应力做功的

40、功率。其中谈到在等速流动中的管道截面上不存在表面粘性力的做功问题，正是基于等速流动时 )。0 x/x0 xx0 xx正应力与压力正应力与压力 由于粘性正应力的存在，流动流体的压力在数值上一般不等于正应力值。比如，对于附加粘性正应力 ，则相应有 。但是，如果将本构方程（6-18）中的三个关系式相加则可得到0、00 xx、p=、p、p|xx 这说明，虽然流动流体的三个正应力在数值上一般不等于压力值，但它们的平均值却总是与压力大小相等的。特别的，对于不可压缩流体的一维流动，设流动沿x方向，则因为 ，且根据连续性方程又有 ，所以有3)(pzzyyxx0=、0=zy0 x/x（6-19）这说明不可压缩流

41、体做一维流动时，正应力与压力的关系与流体静止情况相同，即流体中三个方向正应力的大小都分别与压力相等。这正是第5章分析一维流动时在微元体表面上直接标出压力作为法向表面力的原因。pzzyyxx（6-20）切应力与角变形率切应力与角变形率 在切应力关系中，也有三个变形速率，将它们分别除以可得xz21yz21xy21zxyx、这个变形速率分别是流体在x-y、y-z、z-x平面内的角变形率，即单位时间内两流体线夹角的相对变化率（见第7章）。可见切应力是与其角应变相关是类似的。特别的，对于x-y平面内沿方向的一维不可压缩稳态流动，因 ，所以 仅为的函数，于是本构方程的切应力关系式简化为：，即牛顿剪切公式。

42、0t/、0 x/、00 xxzy、x)dy/d()y/(xxyx 牛顿流体本构方程反映了流体应力与变形速率之间的关系，这与固体力学中反映应力与应变关系的虎克定律是对应的，故牛顿流体本构方程可看成是流体力学中的虎克定律。.流体运动微分方程Navier-Stokes方程 将上述流体应力与变形速率之间的关系牛顿流体本构方程（6-18）代入以应力表示的运动方程（6-17），即可得到由速度分量和压力表示的粘性流体运动微分方程耐维-斯托克斯方程（Navier-Stokes equations，简称N-S方程方程）21)-(6 yz yxz x zz2)(z32zp-fDtDyz zxy x yy2)(y3

43、2yp-fDtDxz zxy y xx2)(x32xp-fDtDzyzxzzzzyyxyyyzxyxxxx N-S方程是现代流体力学的主干方程，几乎所有有关粘性流体流动问题的分析研究工作都是以该方程为基础的。N-S方程对流体的密度、粘度、可压缩性未作限制。但由于引入了牛顿流体的本构方程，故该方程只适用于牛顿流体，对于非牛顿流体，可以采用以应力表示的运动方程。特别的，如果在N-S方程中令=0，即可得到理想流体的运动方程欧拉方程欧拉方程。如果在N-S方程中所有速度项为零，即得到流体静力学方程流体静力学方程。为了应用上的方便，在此给出常见条件下N-S方程的表达式。(1)常粘度下的常粘度下的N-S方程

44、方程 对于等温或温度变化较小的流动，可将粘度视为常数，即=const,相应的N-S方程为22)-(6 zy x z3zp1-fzDtDzy x y3yp1-fyDtDzy x x3xp1-fxDtD2z22z22z2z2y22y22y2y2x22x22x2x所以x方向的运动方程简化为流体在点的密度为，速度为，其x、y、z方向的分量分别为 。引入新变量后，x方向运动可变换为下述方程以应力表示的运动方程需要补充方程才能求解，与第章分析一维流动问题时需要引入补充方程才能由切应力方程得到速度微分方程有些相似。对于流场中的微元体，质量守恒原理可以类似于控制体仿照式（4-9）表述为如果在N-S方程中所有速

45、度项为零，即得到流体静力学方程。流体运动微分方程Navier-Stokes方程动量流量及动量变化率对于另一种情况，即内筒转动外筒固定的情况，由于离心力与压力推动力均指向外壁，两种因素都促进流体向外层运动，所以使得流体沿切向的层流流动难以保持稳定。（注：图6-5中标注的是动量输入或输出方向，而动量或其通量本身的方向均指向方向，即分速度 的方向）。其中谈到在等速流动中的管道截面上不存在表面粘性力的做功问题，正是基于等速流动时 )。如图6-6所示，由于只受到重力场的作用，所以管道截面上任意一点处r、z方向单位质量流体的体积力为试写出该条件下的连续性方程和运动微分方程，并证明管道截面上任意一点的总势能

46、p/+g（rsin）和轴向压力梯度 为常数。其中，是质量通量的的散度，是矢量算子。无论是牛顿流体还是非牛顿流体、是层流流体还是湍流流体，该方程均适用。但与固体应力不一样的是，流体的应力不是与应变的大小而是与应变的速率（即单位间内的应变）直接相关的。不过，对于工程实际问题由于总有特殊性使方程得到简化，从而有可能获得准确或近似的分析解。加速时 ，所以 ；其中，分别为r、坐标方向的速度分量；由此可见，附加粘性正应力的产生是速度沿流动方向的变化所导致的。现以图-所示的微元体为对象，考察流体方向动量在微元体表面的输入输出。由此并根据上述关于应力下标和方向的规定，不难标出微元体个表面上和方向的正应力和切应

47、力，如图-所示。或写成矢量形式为(2)不可压缩流体的不可压缩流体的N-S方程方程)(31p1fDtD2其中，为运动粘度；称为拉普拉斯算子。/2222222zy x 对于不可压缩流体=const，且 如果认为流动等温或温度变化较小，将粘度时为常数，则相应的N-S方程为0（6-23）或简写成矢量形式为24)-(6 zy x zp1-fzDtDzy x yp1-fyDtDzy x xp1-fxDtD2z22z22z2z2y22y22y2y2x22x22x2x2p1-fDtD（6-25）由于通常遇到的流动问题按不可压缩和常粘度问题处理，所以为使用方便，特在此将常粘度条件下不可压缩流体的N-S方程写为展

48、开形式2z22z22z2zzzyzxz2y22y22y2yzyyyxy2x22x22x2xzxyxxxzy x zp1-fzzy x tzy x yp1-fyzy x tzy x xp1-fxzy x t（6-26）该方程的矢量形式以及方程各项的意义如下2p1f)(t（6-27）非定常项 定常流动为0静止流场为0 对流项 静止流场为0 蠕变流时0 单位质量流体的体积力单位质量流体的压力差扩散项（粘性力项）对静止或理想流体为0高速非边界层问题0.柱坐标和球坐标系中的-方程 在工程实际中，有时采用柱坐标和球坐标描述问题比采用直角坐标更为方便，比如，对于常见的圆管内的流动，最适用的显然是柱坐标系统。

49、为此，将不加推导地写出这两种坐标系下常密度和常粘度流体的运动微分方程和牛顿流体本构方程。（1）柱坐标系中的-方程和牛顿流体本构方程 -方程方程 对于以r为径向坐标、为周向坐标、z为轴向坐标的柱坐标体系见图6-2(a)，其粘性流体运动微分方程在r、z方向的分量式为（=const，=const）2z222z22zzzzzzrz22r2222zrr2r222r22rrrz2rrrrzr2r1rrrr1 zpfzrrtzr2r1)r(rr1r pr1fzrrrtzr2r1)r(rr1r rpfzrrrt（6-33）方向方向r方向方向z方向方向 其中，分别为r、z坐标方向的速度分量。此外，r方向分量式中

50、的 和方向分量式中的 分别是单位质量的流体受到的离心力和哥氏力（Corilis force）。这两个力是由直角坐标系转换到柱坐标时自动产生的，在分析流体所受的体积力时不要再人为地加上该力。/r-2/rrzr、rzr1zrrr1z232prr1232pr232pzrrzzrzzzrr rzzzrrrr牛顿流体本构方程 本构方程用于流体应力的分析与计算（6-34）zrr1)r(rr1zr其中：（2）球坐标系中的球坐标系中的-方程和牛顿流体本构方程方程和牛顿流体本构方程 -方程方程 在球坐标体系中，若以r为径向坐标、为周向坐标、为经向坐标见图6-2（b），则运动微分方程在r、方向的分量式为（=con