流体力学(热能)第5章 不可压缩流体动力学基础.

1、第七章第七章 不可压缩流体动力学基础不可压缩流体动力学基础本章讨论三元流动，主要内容是有关流体运动的基本概念和基本原理，本章讨论三元流动，主要内容是有关流体运动的基本概念和基本原理，以及不可压缩流体流动的基本方程。以及不可压缩流体流动的基本方程。积分形式的基本方程用于解决控制面上的流动参数问题。积分形式的基本方程用于解决控制面上的流动参数问题。微分方程可用于解决流微分方程可用于解决流 动参数在流场中的分布问题。动参数在流场中的分布问题。 平移运动 旋转运动 线变形、 变形运动 角变形dxABDCEFMuyuxdy7-1 流体微团运动的分析流体微团运动的分析一、运动形式一、运动形式1、流体微团、

2、流体微团：指体积微小，随流体一起运动的一团流体物质。与流体质点不：指体积微小，随流体一起运动的一团流体物质。与流体质点不同，虽体积微小，但包含无数个流体质点。各质点间存在着相对位置的变化。同，虽体积微小，但包含无数个流体质点。各质点间存在着相对位置的变化。2、基本运动形式、基本运动形式二、运动分析二、运动分析 以二元流动的情况为例，研究几种以二元流动的情况为例，研究几种 基本运动形式的速度表达式。基本运动形式的速度表达式。 如图，方形流动微团如图，方形流动微团xuMABCD2dxxuuxxyu2dxxuuyy2dyyuuxx2dyyuuyy2dxxuuxx2dxxuuyy2dyyuuxx2dy

3、yuuyyxuyuzu各侧边中点各侧边中点A、B、C、D的流速分量分别为的流速分量分别为1、平移运动速度、平移运动速度2、线变形速度、线变形速度xudxdtdxdtxuxxxX方向线变形速度方向线变形速度dxxuxA、C两点速度差值：两点速度差值：差值为正，发生伸长变形差值为正，发生伸长变形。xuxxyuyyzuzzABDCEFM)(21)(21)(21zuyuxuzuyuxuyzxzxyxyz222zyxzyxkji大小xzy3、旋转角速度、旋转角速度 逆时针为正逆时针为正 对角线对角线EMF的旋转角速度定义为的旋转角速度定义为 整个流体微团在整个流体微团在oxy平面上的旋转角速度。平面上的

4、旋转角速度。方向：右手定则方向：右手定则)(21)(21yuxuyuxuxuxuxyxyyzyz)(21)(21)(21yuxuxuzuzuyuxyzzxyyzx4、角变形速度、角变形速度z直角边直角边AMC（或（或BMD）与对角线）与对角线EMF的夹角的变形速度定义的夹角的变形速度定义为流体微团的角变形速度，记为为流体微团的角变形速度，记为 ，表示在，表示在xoy平面上的角变形速度。平面上的角变形速度。三元流动：三元流动：的下标表示发生角变形的所在平面的法线方向。的下标表示发生角变形的所在平面的法线方向。三、亥姆霍兹速度分解定理三、亥姆霍兹速度分解定理 （了解）（了解）设流体微团内某点设流体

5、微团内某点M0(x,y,z),速度为速度为 、 、 ，则邻边则邻边M0的另一点的另一点M (x+dx,y+dy,z+dz)的速度为的速度为 展开展开 .，变换整理得，变换整理得 0 xu0yu0zuxxxduuu0yyyduuu0zzzduuu0 xdudzdydxdzdyuuyzxyzxx0dxdzdydxdzuuzxyzxyy0dydxdzdydxuuxyzxyzz0平移旋转线变形角变形0yu0 xu0zu0MMkjizyx2zuyuyzxxuzuzxyyuxuxyztzyx,2涡量场xzy7-2 有旋流动有旋流动0zyxxyz无旋流动：无旋流动： （详见第（详见第8章）章）有旋流动有旋流

6、动 ： 、 、 至少有一个不等于零。至少有一个不等于零。 涡量：涡量： 涡线，涡线方程，涡量连续性方程涡线，涡线方程，涡量连续性方程 涡通量：斯托克斯方程，汤姆逊定理涡通量：斯托克斯方程，汤姆逊定理1、涡量、涡量u0zyxzyx0)(u ：哈米尔顿算子，矢性微分算子zyxdzdydx3、涡线，、涡线， 涡线方程涡线方程 涡线：表示某一瞬时流体质点旋转角速度向量方向的曲线。涡线：表示某一瞬时流体质点旋转角速度向量方向的曲线。2、涡量连续性微分方程、涡量连续性微分方程二、涡通量二、涡通量 涡量在涡量在 投影为投影为 ，则，则 为涡通量。为涡通量。nndAAn1、定义：、定义：A1A2Anndxdy

7、dzdxdydzdAAdJAzAyAxAnA对有旋转流动，在同一瞬间，通过同一涡管的各截面的涡通量相等。对有旋转流动，在同一瞬间，通过同一涡管的各截面的涡通量相等。dAdAAnAn21dzudyudxusduzysxsdAdAdAdAdxdyyuxudzdxxuzudydzzuyudzudyudxuAnzzyyxAxAxyzxyzzysx2、涡通量的计算、涡通量的计算 （1）速度环量：流速沿封闭曲线）速度环量：流速沿封闭曲线s的积分。的积分。s正向为逆时针方向。正向为逆时针方向。斯托克斯公式：斯托克斯公式：AsJ0dtdn s为流场中任意封闭曲线为流场中任意封闭曲线 A是是S所围成的曲面所围成

8、的曲面 是曲面是曲面A的外法线单位向量。的外法线单位向量。结论：沿任意封闭曲线结论：沿任意封闭曲线S的速度环量等于通过以该曲线的速度环量等于通过以该曲线为边界的曲面为边界的曲面A的涡通量。的涡通量。斯托克斯定理斯托克斯定理（2）汤姆逊定理）汤姆逊定理 （了解）（了解） 在理想流体的涡量场中，如果质量力具有单值的势函数，那么沿在理想流体的涡量场中，如果质量力具有单值的势函数，那么沿由流体质点所组成的封闭曲线的速度环量不随时间而变。即：由流体质点所组成的封闭曲线的速度环量不随时间而变。即：推论：质量力具有单值势函数的理想流体的流动，如果在某一推论：质量力具有单值势函数的理想流体的流动，如果在某一时

9、刻是有旋流，那么此前、此后也是有旋流。如果为无旋流，时刻是有旋流，那么此前、此后也是有旋流。如果为无旋流，那么此前、此后也是无旋流。那么此前、此后也是无旋流。自学：自学：P186例例7-20zuyuxuzyxdxdydzdtxudydzdtdxxuudydzdtdxxuuxxxxx)2()2(7-3 不可压缩流体连续性微分方程不可压缩流体连续性微分方程净流体体积净流体体积=流出流入流出流入2、分析推导：微元分析法、分析推导：微元分析法 依据质量守恒定律，取微小平行六面体，中心依据质量守恒定律，取微小平行六面体，中心M（x,y,z)， ux、uy、uz， x方向，方向，dt时间时间0zururu

10、ruzrr（柱面坐标形式）（柱面坐标形式）1、方程：、方程：同理同理y方向：方向：dxdydzdtyuyz方向：方向：dxdydzdtzuz根据不可压缩连续性条件，根据不可压缩连续性条件，dt时间内，时间内，x、y、z方向方向 流出流出- -流入流入=000)(zuyuxudxdydzdtzuyuxuzyxzyx0zurururuzrrcossinsincossincosuuuuuuzzryrxryrxcossinsincosyxyxruuuuuu柱面坐标系：柱面坐标系：柱面坐标与直角坐标换算关系：柱面坐标与直角坐标换算关系：190 图图7-73、应用、应用（1）方程对恒定流、非恒流都适用，是

11、判断流动连续性的条件。）方程对恒定流、非恒流都适用，是判断流动连续性的条件。（2）是以后运动微分方程求解的一个条件。）是以后运动微分方程求解的一个条件。例例1：p1907-6例例2：判断下列流场是否满足不可压缩流体的连续性方程：判断下列流场是否满足不可压缩流体的连续性方程（1） （2）解：解： 满足。满足。22)2(xyyuxxyuyx2sin212cos2rurrur02cos412cos212cos222rrzurururuzrr（3）此式给出了流体通过某固定点时，流体的三个速度分量之间的关系。）此式给出了流体通过某固定点时，流体的三个速度分量之间的关系。表明对不可压缩流体，单位时间内流入

12、与流出某空间点的流体体积之差为表明对不可压缩流体，单位时间内流入与流出某空间点的流体体积之差为零，即体积（质量）守恒。零，即体积（质量）守恒。xzxyxxp,zzzyzxp,yzyyyxp,3个压应力6个切应力xzxxpxyxyz7-4 粘性流体的运动微分方程粘性流体的运动微分方程 纳维纳维斯托克斯方程斯托克斯方程一、粘性流体的内应力一、粘性流体的内应力粘性流体运动时，所受表面应力包括法向应力和切应力。粘性流体运动时，所受表面应力包括法向应力和切应力。流场内任一点的应力可表示为流场内任一点的应力可表示为二、以应力表示的运动微分方程（根据牛顿第二定律）二、以应力表示的运动微分方程（根据牛顿第二定

13、律）dtduzyxpXxzxyxxx11dtduxzypYyxyzyyy11dtduyxzpZzyzxzzz11（7-4-1）讨论讨论：加上连续性方程，：加上连续性方程，4个方程，个方程，12个未知量，无法求解。个未知量，无法求解。需找其它关系式，这些其它关系即是应力与变形速度的关系。需找其它关系式，这些其它关系即是应力与变形速度的关系。dtdxuyudtdyxz21、切应力与角变形速度的关系、切应力与角变形速度的关系因此三元流动的牛顿内摩擦定律可以写成如下形式：因此三元流动的牛顿内摩擦定律可以写成如下形式：由牛顿内摩擦定律：由牛顿内摩擦定律：xoy平面上平面上:三、应力与变形速度的关系三、应

14、力与变形速度的关系)()()(zuyuzuxuxuyuyzyzzyxzxzzxyxyxxy式（式（7-5-1）式（式（7-5-1） 即广义牛顿摩擦定律，使式（即广义牛顿摩擦定律，使式（7-4-1）中的）中的12个未知数消去个未知数消去6个。个。2、法向应力和线变形速度的关系、法向应力和线变形速度的关系 在流体微团的法线方向上的线变形速度，使法向应力（压应力）的大在流体微团的法线方向上的线变形速度，使法向应力（压应力）的大小与理想流体相比有所改变，产生附加压应力。小与理想流体相比有所改变，产生附加压应力。可以证明，对于不可压缩流体，附加法应力与线变形速度的关系可以证明，对于不可压缩流体，附加法应

15、力与线变形速度的关系:zuyuxuzzzyyyxxx222)(31zzyyxxppppzuppyuppxuppztzzytyyxtxx222（1）(7-5-3)（2）平均压应力）平均压应力定义点压强定义点压强：(7-5-4)式中，式中， 、 、 表示法向应力，表示法向应力， 表示压强，表示压强， 表示理想流体压强。表示理想流体压强。xxpyypzzpptp（3）)(32)(31zuyuxupppppzyxtzzyyxx(7-5-5)代入（代入（7-5-4）（4）)(322)(322)(322zuyuxuzuppzuyuxuyuppzuyuxuxuppzyxzzzzyxyyyzyxxxx（7-5

16、-6）tzzyyxxpppp0zuyuxuzyx0),(zyxuuzyuu0zuyuxuzyxtzzyyxxpppp a）对于理想流体）对于理想流体b）不可压缩流体）不可压缩流体c ） 对于均匀流对于均匀流 流速沿流线是常数流速沿流线是常数d ）对于粘性流体方程）对于粘性流体方程 、 、 三个未知数变为一个三个未知数变为一个 ，原则上方程已可求解，原则上方程已可求解了。了。 xxpyypzzpp（7-5-6）式讨论：）式讨论：tpp 四、四、 N-S方程方程把（把（7-5-1）式和（）式和（7-5-6）式代入（）式代入（7-4-1）式，消去应力）式，消去应力对不可压缩流体有对不可压缩流体有 代

17、入得代入得展开展开 当地加速度当地加速度 位移加速度位移加速度0zuyuxuzyxdtduzuyuxuzpZdtduzuyuxuypYdtduzuyuxuxpXzzzzyyyyxxxx)(1)(1)(1222222222222222222zuuyuuxuutudtduxzxyxxxx（7-6-1）zuuyuuxuutuzuyuxuzpZzuuyuuxuutuzuyuxuypYzuuyuuxuutuzuyuxuxpXzzzyzxzzzzyzyyyxyyyyxzxyxxxxxx)(1)(1)(1222222222222222222（7-6-3）1、N-S方程方程柱坐标系：柱坐标系：P1977-6-

18、4式式2、讨论：、讨论：（1）与连续性方程联立，）与连续性方程联立，4个未知数、个未知数、4个方程，原则上可求解速个方程，原则上可求解速度分量和压强。度分量和压强。 N-S方程是不可压缩流体最普遍的运动微分方程。方程是不可压缩流体最普遍的运动微分方程。（2）二阶非线性非齐次的偏微分方程组，难以求出精确解。只对）二阶非线性非齐次的偏微分方程组，难以求出精确解。只对一些较简单的情况求出精确解。大多数问题是借助计算机技术求一些较简单的情况求出精确解。大多数问题是借助计算机技术求出近似解。出近似解。 P198例例7-7 自学自学例：已知流速场例：已知流速场 试求试求t=0.5时空间点（时空间点（2,5

19、,3）处的流体质点的加速度。）处的流体质点的加速度。 解：解：0,zyxuxztuyztu5 .1922txzyzzuuyuuxuutuaxzxyxxxx25.1722tyzxzzuuyuuxuutuayzyyyxyy0zuuyuuxuutuazzzyzxzz04.26222zyxaaaa0zuuyuuxuutuzpZzuuyuuxuutuypYzuuyuuxuutuxpXzzzyzxzyzyyyxyxzxyxxx1110zyxuuu010101zpZypYxpX7-5 理想流体运动微分方程及其积分理想流体运动微分方程及其积分分析分析N-S方程方程1、理想流体、理想流体（7-7-1）2、静止时

20、、静止时得流体平衡微分方程（得流体平衡微分方程（2-7-1a）3、变换（、变换（7-7-1）式，在方程中第一式的加速度项加）式，在方程中第一式的加速度项加 之后，之后，整理得：整理得： 同理变换第二式、第三式，得：同理变换第二式、第三式，得： （7-7-2）恒定流：恒定流： ，并设质量力有势函数，并设质量力有势函数W，则，则zuuyuuxzxy,)(2)2(12yzzyxuuuxtuxpX)(2)2(1)(2)2(1)(2)2(1222xyyxzzxxzyyzzyxuuuztuzpZuuuytuypYuuuxtuxpX0tututuzyxzyzyuuupWx2)2(2 分别乘以分别乘以dx、d

21、y、dz，相加得，相加得若若 ， ，则，则 ，理想恒定流能量方程。，理想恒定流能量方程。xzxzuuupWy2)2(2yxyxuuupWz2)2(2dzuudyuudxuuupWdyxyxxzxzzyzy222)2(2zyxzyxuuudzdydxupWd2)2(2gzW0)2(2upWdconst22gupz所以所以 是导出理想恒定流能量方程的条件。是导出理想恒定流能量方程的条件。行列式等于零，则任一行全等于零或任两行成比例。行列式等于零，则任一行全等于零或任两行成比例。讨论：讨论：（1） 流体静止。流体静止。（2） 流线方程流线方程同一条流线上同一条流线上 （3） 无旋流无旋流 无旋流空间

22、各点，处处满足能量方程。无旋流空间各点，处处满足能量方程。（4） 涡线方程涡线方程 涡线上满足理想流体能量方程。涡线上满足理想流体能量方程。0zyxzyxuuudzdydx0, 0, 0zyxuuuzyxudzudyudxcgupz22元流能量方程0, 0, 0zyxzyxdzdydx（5） 螺旋流动螺旋流动 （涡线与流线相重合）（涡线与流线相重合）螺旋流动中，全部流动均满足理想流体能量方程。螺旋流动中，全部流动均满足理想流体能量方程。 kuuuzzyyxx7-6 流体运动的初始条件和边界条件流体运动的初始条件和边界条件解二阶偏微分方程，需确定方程的定解条件（初始条件和边界条件）。解二阶偏微分

23、方程，需确定方程的定解条件（初始条件和边界条件）。目前，计算流体力学已广泛应用于解决工程中的流动问题，如何正确目前，计算流体力学已广泛应用于解决工程中的流动问题，如何正确合理给出初始条件和边界条件尤为重要。具体条件依赖于具体的流动。合理给出初始条件和边界条件尤为重要。具体条件依赖于具体的流动。以粘性不可压缩流体流动为例以粘性不可压缩流体流动为例1、初始条件、初始条件zyxutzyxuuxxx,00若恒定流动，不必给出。若恒定流动，不必给出。2、边界条件、边界条件边界包括固体壁面，两种流体介质的分界面，管道的出入口等。边界包括固体壁面，两种流体介质的分界面，管道的出入口等。（1）固体壁面静止）固体壁面静止0),(fzyxuuu固壁无滑移条件固壁无滑移条件（2）不同液体的分界面）不同液体的分界面 ， 两侧液体的速度、压强保持连续