第6章 流体流动微分方程ppt课件.ppt

《第6章 流体流动微分方程ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

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

2、这一基本方法推广应用于三维情况，建立一般条件下的流体运动微分方程。就其目的而言，积分方程反映流动过程中流体总质量、总动量和总能量的变化，而本章要建立的流动微分方程，目的在于流场分布的详细信息，以揭示宏观流动现象的内在规律。第六章流体流动微分方程６.１　连续性方程６.２　以应力表示的运动方程６.３　粘性流体运动微分方程６.４　流体流动微分方程的应用6.1.1直角坐标系中的连续性方程连续性方程反映流动过程遵循质量守恒这一事实。对于流场中的微元体，质量守恒原理可以类似于控制体仿照式（4-9）表述为（6-1）为了获得（6-1

3、）的数学表达式—连续性方程，不妨对图６-１所示的微元体进行分析。该微元体取自流场中的任意点Ａ，微元体在x、y、z方向的边长分别为dx、dy、dz，其六个面两两相互平行且分别垂直于x、y、z。流体在Ａ点的密度为ρ，速度为ν，其x、y、z方向的分量分别为。一般而言，速度ν和密度ρ均为坐标x、y、z和时间ｔ的函数。已知，流体穿越某一表面时的质量流量等于质量通量与表面积的乘积，而质量通量则为流体密度与流体在该表面上的法向速度的乘积，因此考察微元体上的输入与输出，首先要确定微元面上的法向速度。如图６-１所示,对于在流场中任意点

4、Ａ所取的微元体，因为与Ａ点相邻的三个微元面上的流体或流动参数反映的是Ａ点的参数，所以在这三个微元面上，流体密度均为ρ（Ａ点密度），且每一个面上流体的三个速度分量都为（Ａ点速度）。其中，对于dydx微元面，因其与ｘ轴垂直，该微元面上的三个速度分量中，是法向速度，其产生的法向通量为，而另为两个速度分量则平行于dydz平面，与质量输出输入无关（故图中dydz微元面上的质量通量）。同理，在垂直于y、z方向的微元面dxdz和dxdy上，法向速度分别为，质量通量分别为，如图6-1所示。按速度与坐标方向一致为正的约定，、都是输入

5、通量。于是将这三个通量分别乘以相应的面积dydz、dxdz、dxdy后相加，可得输入微元体的质量流量为相应地，当流体从与Ａ点不相邻的、分别垂直于ｘ、ｙ、ｚ方向的三个微元面上流出时，由于分别经过dx、dy、dz的距离后其输出时的质量通量将发生变化，如图６-１所示所以输出微元体的质量流量为由上述两项可得两一方面，对于图６－１所示的微元体，其瞬间质量为ρdxdydz，所以（６-２）（6-3）连续性方程将式（6-2）和式（6-3）代入微元体质量守恒文字表达式（6-1）可得直角坐标系中的连续性方程为或以矢量简洁表示为（6-4a

6、）（6-4b）其中，是质量通量的ρν的散度，是矢量算子。由于导出方程（6-4）的过程中没有对流体和流动状态作任何假设，故该方程对层流和湍流、牛顿流体和非牛顿流体均适用。将方程（6-4a）展开并引用第2章中的随体导数（质点导数）概念，可将连续性方程表示为另一种形式其中，是速度矢量ν的散度；是密度ρ随体导数，按第2章中随体导数的定义有或（6-5a）（6-5b）不可压缩流体的连续性方程对于不可压缩流体，因密度ρ=const，所以连续性方程简化为或（6-6a）（6-6b）在物理意义上，速度的散度表示为单位体积的流量在单位时间

7、内的体积增量，通常称为体变形率．对于不可压缩流体，不管其体积形状如何变化，其体积的大小不会变，故体变形率为零，即。也正是这一特点，对于不可压缩流体，无论是为稳态流动还是非稳态流动，其连续性方程都是一样的。不可压缩流体的连续性方程不仅形式简单，而且应用广泛，因为工程实际中除了经常遇到不可压缩流体外，不少可压缩流体的流动亦可常密度流体处理。由连续性方程(6-6)可知，对于不可压缩流体沿ｘ方向的一维流动，，其连续性方程就是。这正是第５章中分析不可压缩流体一维流动时曾经用到的条件。6.1.2柱坐标和球坐标系中的连续性方程在工

8、程实际中，除了直角坐标系外，出于描述的方便还经常采用柱坐标（如圆管流动问题）和球坐标（如球体绕流问题）。在此不加推导地写出这两种坐标系中的连续性方程，以供使用。图6－2对于以ｒ为径向坐标、θ为周向坐标、ｚ为轴向坐标的柱坐标体系见图6－2（a），其连续性方程为（6-7）其中，分别为r、θ、z坐标方向的速度分量。特别地，对于不可压缩流体，柱坐标系下