工程流体力学第七章理想不可压缩流体的有旋流动和无旋流动

1、本章内容：本章内容：在许多工程实际问题中，流动参数在许多工程实际问题中，流动参数不仅在流动方向上发生变化，而且在垂直于流不仅在流动方向上发生变化，而且在垂直于流动方向的横截面上也要发生变化。要研究此类动方向的横截面上也要发生变化。要研究此类问题，就要用问题，就要用多维流的分析方法多维流的分析方法。本章主要讨。本章主要讨论理想流体多维流动的基本规律，为解决工程论理想流体多维流动的基本规律，为解决工程实际中类似的问题提供理论依据，也为进一步实际中类似的问题提供理论依据，也为进一步研究粘性流体多维流动奠定必要的基础。研究粘性流体多维流动奠定必要的基础。首先推导在笛卡儿坐标系中微分形式的连续性方程首先

2、推导在笛卡儿坐标系中微分形式的连续性方程 0dAvdVtCVCSn(3-22)zv在在 方向上，单位时间通过方向上，单位时间通过EFGHEFGH面流入的流体质量为面流入的流体质量为 xdydzdxvxvxx2单位时间通过单位时间通过ABCDABCD面流出的流体质量面流出的流体质量 dydzdxvxvxx2 则在则在 方向单位时间内通过微元体表面的净通量为（方向单位时间内通过微元体表面的净通量为（b b）- -（a a），），即即 dxdydzvxx（c1c1）x（a a）（b b）同理可得同理可得 和和 方向单位时间通过微元体表面的净通量分别为方向单位时间通过微元体表面的净通量分别为yzdxd

3、ydzvyydxdydzvzz（c2c2） （c3c3） 因此，单位时间流过微元体控制面的总净通量为因此，单位时间流过微元体控制面的总净通量为dxdydzvzvyvxdAvzyxnCS（c c） 微元六面体内由于密度随时间的变化而引起的质量的变化率为：微元六面体内由于密度随时间的变化而引起的质量的变化率为： dxdydztdxdydztdVtCVCV 将式（将式（c c），（），（d d）代入式（代入式（3-223-22），同除以），同除以 ，则可得到微分形，则可得到微分形式的连续性方程式的连续性方程dxdydz0tvzvyvxzyx0)(tv或或 连续性方程表示了单位时间内控制体内流体质量的

4、增量等于流体在控连续性方程表示了单位时间内控制体内流体质量的增量等于流体在控制体表面上的净通量。它适用于理想流体和粘性流体、定常流动和非定常制体表面上的净通量。它适用于理想流体和粘性流体、定常流动和非定常流动。流动。 在定常流动中，由于在定常流动中，由于0t0zyxvzvyvx对于不可压缩流体（对于不可压缩流体（ = =常数）常数）0zvyvxvzyx0 v或或在其它正交坐标系中流场中任一点的连续性方程和柱坐标系中的表示式为在其它正交坐标系中流场中任一点的连续性方程和柱坐标系中的表示式为 ： 0)()(1)(1zrvzvrvrrrt对于不可压缩流体对于不可压缩流体 01rvzvvrrvrzr式

5、中式中 为极径；为极径； 为极角为极角r球坐标系中的表示式为球坐标系中的表示式为: :0sin1)sin(sin1)(122vrvrrrvrtr0cot2sin11rvrvvrvrrvrr式中式中 为径矩；为径矩； 为纬度；为纬度； 为径度。为径度。r【例【例7-17-1】0zvyvxvzyx044zvyxzyxzvz44 ),()(4yxfzyxvz0),(yxfzyxvz)(4vxyyxvx22zyvy220z0zvzzvxy0z0zv 流体与刚体的主要不同在于它具有流动性，极易变形。流体与刚体的主要不同在于它具有流动性，极易变形。因此，流体微团在运动过程中不但象刚体那样可以有移动和转因此

6、，流体微团在运动过程中不但象刚体那样可以有移动和转动，而且还会发生变形运动。一般情况下，流体微团的运动可动，而且还会发生变形运动。一般情况下，流体微团的运动可以分解为移动，转动和变形运动。以分解为移动，转动和变形运动。 dxDAdyyBCx（1 1）平移运动：如图所示，）平移运动：如图所示，矩形矩形ABCDABCD各角点具有相同的速度各角点具有相同的速度 。导致矩形导致矩形ABCDABCD平移平移x x = = dtdt, , y y = = dtdt, , 其其ABCDABCD的形状不变。的形状不变。（2 2）线变形运动：如图所示，线变形运动取决于速度分量在它所在方向）线变形运动：如图所示，

7、线变形运动取决于速度分量在它所在方向上的变化率（即线变形速率上的变化率（即线变形速率 和和 ），），导致矩形导致矩形ABCDABCD的变形量：的变形量：yxvv ,xvyvxvxyvydtdxxvxx22dtdyyvyy22（3 3）角变形运动和旋转运动：如图）角变形运动和旋转运动：如图7-47-4（c c）、（）、（d d）所示，当所示，当 dtxvdxdtdxxvyy)2(2tandtyvdydtdyyvxx)2(2tanxvyvyx当当矩形矩形ABCDABCD只发生角变形运动，如图所示。只发生角变形运动，如图所示。 xvyvyx当当矩形矩形ABCDABCD只发生旋转运动，形状不变只发生旋

8、转运动，形状不变在一般情况下在一般情况下 xvyvyx的同时，还会发生角变形运动。如图的同时，还会发生角变形运动。如图(e)(e)所示。所示。亦就是亦就是矩形矩形ABCDABCD在发生旋转运动在发生旋转运动在角变形运动和旋转运动同时发生的情况下，将会有以下关系式：在角变形运动和旋转运动同时发生的情况下，将会有以下关系式： 2121于是沿于是沿z z轴流体微团的旋转角速度分量：轴流体微团的旋转角速度分量： yvxvdtdtdtxyz2121同理，沿同理，沿x x，y y轴流体微团的旋转角速度分量分别为轴流体微团的旋转角速度分量分别为: : zvyvyzx21xvzvzxy21流体微团的旋转角速度

9、定义为流体微团的旋转角速度定义为: : Vkjizyx21其中，流体微团的旋转角速度分量及模量为：其中，流体微团的旋转角速度分量及模量为： yvxvxvzvzvyvxyzzxyyzx212121222zyx 流体微团沿流体微团沿z z轴的角变形速度分量：轴的角变形速度分量： yvxvdtdtdtxyz2121同理，可有流体微团角变形速度之半分量及其模量为：同理，可有流体微团角变形速度之半分量及其模量为： yvxvxvzvzvyvxyzzxyyzx212121222zyx 前面在流体微团的分析中，已给出前面在流体微团的分析中，已给出E E点的速度点的速度为为 : :dzzvdyyvdxxvvvd

10、zzvdyyvdxxvvvdzzvdyyvdxxvvvzzzzzEyyyyyExxxxxE如果在上式的第一式右端加入两组等于零的项：如果在上式的第一式右端加入两组等于零的项： dyxvdyxvyy2121dzxvdzxvzz2121其值不变。经过简单组合，可将该式写成其值不变。经过简单组合，可将该式写成 ：dzxvzvdyyvxvdzxvzvdyyvxvdxxvvvzxxyzxxyxxxE)(21)(21)(21)(21同理，有：同理，有： dyzvyvdxxvzvdyzvyvdxxvzvzzvvvdxyvxvdzzvyvdxyvxvdzzvyvyyvvvyzzxyzzxzzzExyyzxyy

11、zyyyE)(21)(21)(21)(21)(21)(21)(21)(21dzzvdyyvdxxvvvxxxxxE将旋转角速度分量，角变形速度之半代入以上三式，便可将式写将旋转角速度分量，角变形速度之半代入以上三式，便可将式写成成 ：)22()22(2)22()22(2)22()22(dxdydydxdzzvvvdzdxdxdzdyyvvvdydzdzdyxdxxvvvyxxyzzzExzzxyyyEzyyzxxxE 上式表明：各速度分量的第一项是平移速度分量，第二、上式表明：各速度分量的第一项是平移速度分量，第二、三、四项分别是由线变形运动、角变形运动和旋转运动所引起三、四项分别是由线变形运

12、动、角变形运动和旋转运动所引起的线速度分量。此关系也称为海姆霍兹的线速度分量。此关系也称为海姆霍兹( (HelmholtzHelmholtz) )速度分解速度分解定理，该定理可简述为：定理，该定理可简述为： 在某流场在某流场O O点邻近的任意点点邻近的任意点A A上的速度可以分成三个部分：上的速度可以分成三个部分：分别为分别为与与O O点相同的平移速度点相同的平移速度（平移运动）；（平移运动）；绕绕O O点转动在点转动在A A点点引起的速度引起的速度（旋转运动）；（旋转运动）；由于变形（包括线变形和角变形）由于变形（包括线变形和角变形）在在A A点引起的速度点引起的速度（变形运动）。（变形运动

13、）。 根据流体微团在流动中是否旋转，可将流体的流动分为两根据流体微团在流动中是否旋转，可将流体的流动分为两类：有旋流动和无旋流动。类：有旋流动和无旋流动。 数学条件：数学条件： 当当 021V021V当当 无旋流动无旋流动 有旋流动有旋流动 通常以通常以 是否等于零作为判别流动是否有旋或无旋是否等于零作为判别流动是否有旋或无旋的判别条件。的判别条件。 V 在笛卡儿坐标系中：在笛卡儿坐标系中： kyvxvjxvzvizvyvVxyzxyz即当流场速度同时满足：即当流场速度同时满足： zvyvyzxvzvzxyvxvxy时流动无旋。时流动无旋。 需要指出的是，需要指出的是，有旋流动和无旋流动仅由流

14、体微团本身是否有旋流动和无旋流动仅由流体微团本身是否发生旋转来决定，而与流体微团本身的运动轨迹无关。发生旋转来决定，而与流体微团本身的运动轨迹无关。 如图（如图（a a），），流体微团的运动为旋转的圆周运动，其微团自流体微团的运动为旋转的圆周运动，其微团自身不旋转，流场为无旋流动；图（身不旋转，流场为无旋流动；图（b b）流体微团的运动尽管为流体微团的运动尽管为直线运动，但流体微团在运动过程中自身在旋转，所以，该流直线运动，但流体微团在运动过程中自身在旋转，所以，该流动为有旋流动。动为有旋流动。（a） （b） 【例例7-27-2】 某一流动速度场为某一流动速度场为 ， ，其中，其中 是是不为零

15、的常数，流线是平行于不为零的常数，流线是平行于 轴的直线。试判别该流动是有轴的直线。试判别该流动是有旋流动还是无旋流动。旋流动还是无旋流动。 【解】【解】 由于由于 021zvyvyzx02121ayvxvxyz所以该流动是有旋运动。所以该流动是有旋运动。 ayvx0zyvvaxx021xvzvzxy 理想流体运动微分方程式是研究流体运动学的重要理论基理想流体运动微分方程式是研究流体运动学的重要理论基础。可以用础。可以用牛顿第二定律牛顿第二定律加以加以推导推导。 fzyxfff,zyxA (x,y,z)2dxxpp2dxxppf 微元体在质量力和表面力的作用下产生的加速度微元体在质量力和表面力

16、的作用下产生的加速度 ，根，根据牛顿第二定律据牛顿第二定律 ：adtdvmFxxdtdvdxdydzdydzdxxppdydzdxxppdxdydzfxx)2()2(两端同除以微元体的质量两端同除以微元体的质量 , ,并整理有：并整理有： dxdydzdtdvzpfdtdvypfdtdvxpfzzyyxx111写成矢量式：写成矢量式： dtvdpfgrad1将加速度的表达式代入有：将加速度的表达式代入有： zvvyvvxvvtvzpfzvvyvvyvvtvypfzvvyvvxvvtvxpfzzzyzxzzyzyyyxyyxzxyxxxx111其矢量式为其矢量式为 ：vvtvpf)(1grad

17、上式为理想流体运动微分方程式上式为理想流体运动微分方程式，物理上表示了作用在物理上表示了作用在单位质量流体上的质量力、表面力和惯性力相平衡。该式推单位质量流体上的质量力、表面力和惯性力相平衡。该式推导过程中对流体的压缩性没加限制，导过程中对流体的压缩性没加限制，故可适用于理想的可压故可适用于理想的可压流体和不可压缩流体，适用于有旋流动和无旋流动。流体和不可压缩流体，适用于有旋流动和无旋流动。 将上式作恒等变形，便可以直接由运动微分方程判定流动是将上式作恒等变形，便可以直接由运动微分方程判定流动是有旋还是无旋流动，在上式的第一式右端同时有旋还是无旋流动，在上式的第一式右端同时 xvvyy 1xp

18、fxvzvvxvyvvxvvxvvxvvtvxzxzyxyzzyyxxx据旋转角速据旋转角速度分量公式度分量公式 122 122 122222zpfvvvztvypfvvvytvxpfvvvxtvzxyyxzyzxxzyxyzzyxzvvyvvxvvtvxpfxzxyxxxx1xvvzz如果流体是在有势的质量力作用下，流场是正压性的，则：如果流体是在有势的质量力作用下，流场是正压性的，则： xfxyfyzfz此时存在一压强函数此时存在一压强函数: : pPFd将压强函数对坐标的偏导数有将压强函数对坐标的偏导数有: : xpxPF 1ypyPF1zpzPF1将上述关系代入兰姆运动微分方程，得将上

19、述关系代入兰姆运动微分方程，得: : 22 22 22222xyyxzFzxxzyFyzzyxFvvtvPvzvvtvPvyvvtvPvx写成矢量形关系式写成矢量形关系式 vv222tPvF二、欧拉积分二、欧拉积分 当理想正压性流体在有势的质量力作用下作定常无旋流当理想正压性流体在有势的质量力作用下作定常无旋流动时，上式右端为零。若在流场中任取一有向微元线段动时，上式右端为零。若在流场中任取一有向微元线段dldl，其在三个坐标轴的投影分别为其在三个坐标轴的投影分别为dxdx、dydy、dzdz，将它们分别依次将它们分别依次乘上式并相加，得乘上式并相加，得: : 0d2d2d2222zPvzyP

20、vyxPvxFFF02d2FPv 上式为欧拉积分的结果，表明理想正压性流体在有势的质上式为欧拉积分的结果，表明理想正压性流体在有势的质量力作用下作定常无旋流动时，单位质量流体的总机械能在流量力作用下作定常无旋流动时，单位质量流体的总机械能在流场中保持不变。场中保持不变。 积分积分 CPvF22三、伯努利积分三、伯努利积分 当理想正压性流体在有势的质量力作用下作定常有旋流当理想正压性流体在有势的质量力作用下作定常有旋流动时，式右端第一项等于零。由流线的特性知，此时流线与动时，式右端第一项等于零。由流线的特性知，此时流线与迹线重合，在流场中沿流线取一有向微元线段迹线重合，在流场中沿流线取一有向微元

21、线段dldl，其在三个，其在三个坐标轴上的投影分别为坐标轴上的投影分别为dxdx= =v vx xdtdt，dydy= =v vy ydtdt，d dz z= =v vz zdtdt，将它们将它们的左、右端分别依次乘式的左、右端，相加有的左、右端分别依次乘式的左、右端，相加有 02d2FPv积分有积分有 CPvF22 该积分为伯努里积分该积分为伯努里积分。表明理想正压性流体在有势的质表明理想正压性流体在有势的质量力作用下作定常有旋流动时，单位质量流体的总机械能沿量力作用下作定常有旋流动时，单位质量流体的总机械能沿流线保持不变。通常沿不同流线积分常数值有所不同。流线保持不变。通常沿不同流线积分常

22、数值有所不同。 本节本节主要讲述主要讲述理想流体有旋运动的理论基础，重点是速理想流体有旋运动的理论基础，重点是速度环量及其表征环量和旋涡强度间关系的度环量及其表征环量和旋涡强度间关系的斯托克斯定理斯托克斯定理。 涡量用来描述流体微团的旋转运动。涡量的定义为：涡量用来描述流体微团的旋转运动。涡量的定义为： V2涡量是点的坐标和时间的函数。它在直角坐标系中的投影为涡量是点的坐标和时间的函数。它在直角坐标系中的投影为 ：zvyvyzxxvzvzxyyvxvxyz 在流场的全部或部分存在角速度的场，称为在流场的全部或部分存在角速度的场，称为涡量场涡量场。如同在如同在速度场中引入了流线、流管、流束和流量

23、一样。在涡量场中同样速度场中引入了流线、流管、流束和流量一样。在涡量场中同样也引入涡线、涡管、涡束和旋涡强度的概念。也引入涡线、涡管、涡束和旋涡强度的概念。 1 1涡线涡线: :涡线是在给定瞬时和涡量矢量相切的曲线。如图所示。涡线是在给定瞬时和涡量矢量相切的曲线。如图所示。 根据涡通量矢量与涡线相切的条件，涡线的微分方程为根据涡通量矢量与涡线相切的条件，涡线的微分方程为: : ),(),(),(tzyxdztzyxdytzyxdxzyx2 2涡管、涡束：涡管、涡束：在涡量场中任取一不是涡线的封闭曲线，在同一在涡量场中任取一不是涡线的封闭曲线，在同一时刻过该曲线每一点的涡线形成的管状曲面称作涡管

24、。如图时刻过该曲线每一点的涡线形成的管状曲面称作涡管。如图7-87-8所所示。示。截面无限小的涡管称为微元涡管。涡管中充满着的作旋转运截面无限小的涡管称为微元涡管。涡管中充满着的作旋转运动的流体称为涡束，微元涡管中的涡束称为微元涡束或涡丝。动的流体称为涡束，微元涡管中的涡束称为微元涡束或涡丝。3 3旋旋涡强度（涡通量）涡强度（涡通量） 在涡量场中取一微元面积在涡量场中取一微元面积dAdA，见图见图，其上流体微团的涡其上流体微团的涡通量为通量为 ， 为为dAdA的外法线方向，定义的外法线方向，定义 2ndAdAnAddJn2)cos(2为任意微元面积为任意微元面积dAdA上的旋上的旋涡强度，也称

25、涡通量。涡强度，也称涡通量。 任意面积任意面积A A上的旋上的旋涡强度为：涡强度为： dAdAJnAA2 如果面积如果面积A A是涡束的某一横截面积，就称为涡束旋涡强度，是涡束的某一横截面积，就称为涡束旋涡强度，它也是旋转角速度矢量的通量。旋涡强度不仅取决于，而且取它也是旋转角速度矢量的通量。旋涡强度不仅取决于，而且取决于面积决于面积A A。 1 1速度环量速度环量: :在流场的某封闭周线上，如图在流场的某封闭周线上，如图7-97-9（b b），），流体速度流体速度矢量沿周线的线积分，定义为速度环量，用符号矢量沿周线的线积分，定义为速度环量，用符号 表示，即表示，即: : )(dzvdyvdx

26、vldvzyx 速度环量是一代数量，它的正负与速度的方向和线积分的绕速度环量是一代数量，它的正负与速度的方向和线积分的绕行方向有关。行方向有关。对非定常流动，速度环量是一个瞬时的概念，应根对非定常流动，速度环量是一个瞬时的概念，应根据同一瞬时曲线上各点的速度计算，积分时为参变量。据同一瞬时曲线上各点的速度计算，积分时为参变量。 2 2斯托克斯（斯托克斯（StokesStokes）定理定理: :在涡量场中，沿任意封闭周线的在涡量场中，沿任意封闭周线的速度环量等于通过该周线所包围曲面面积的旋涡强度，即速度环量等于通过该周线所包围曲面面积的旋涡强度，即:JdAAdldvnAA2 这一定理将旋涡强度与

27、速度环量联系起来，给出了通过速这一定理将旋涡强度与速度环量联系起来，给出了通过速度环量计算旋涡强度的方法。度环量计算旋涡强度的方法。 【例例7-37-3】已知二维流场的速度分布为】已知二维流场的速度分布为 ， ，试试求绕圆求绕圆 的速度环量。的速度环量。 yvx3xvy4222Ryx【解】【解】 此题用极坐标求解比较方便，坐标变换为此题用极坐标求解比较方便，坐标变换为: : cosrx sinry 速度变换为速度变换为 sincosyxrvvv,sincosxyvvv22sin3cos4rrv2022)sin3cos4(rdrr dr)sin3cos4(20222 2202227cos6rdr

28、r【例【例7-47-4】 一二维涡量场，在一圆心在坐标原点、半径一二维涡量场，在一圆心在坐标原点、半径 的圆区域内，流体的涡通量的圆区域内，流体的涡通量 。若流体微团在。若流体微团在半径半径 处的速度分量处的速度分量 为常数，它的值是多少？为常数，它的值是多少？ mr1 . 0smJ/4 . 02rv【解】由斯托克斯定理得【解】由斯托克斯定理得 ：Jrvrdv202smrJv/21 . 024 . 021 1汤姆孙（汤姆孙（ThomsonThomson）定理定理 理想正压性流体在有势的质量力作用下，沿任何封闭流体周线理想正压性流体在有势的质量力作用下，沿任何封闭流体周线的速度环量不随时间变化，

29、即：的速度环量不随时间变化，即：0dtdJ 假设在理想不可压缩的重力流体中，有一象刚体一样以等角假设在理想不可压缩的重力流体中，有一象刚体一样以等角速度速度 绕自身轴旋转的无限长铅垂直涡束，其涡通量为绕自身轴旋转的无限长铅垂直涡束，其涡通量为J J。涡束涡束周围的流体在涡束的诱导下绕涡束轴等速圆周运动，由斯托克斯周围的流体在涡束的诱导下绕涡束轴等速圆周运动，由斯托克斯定理知，定理知， 。由于直线涡束无限长，该问题可作一个平面问题。由于直线涡束无限长，该问题可作一个平面问题研究。可以证明涡束内的流动为有旋流动，称为涡核区，其半径研究。可以证明涡束内的流动为有旋流动，称为涡核区，其半径为为 ；涡束

30、外的流动区域为无旋流动，称为环流区。；涡束外的流动区域为无旋流动，称为环流区。 br在环流区内，速度分布为在环流区内，速度分布为: : 0rvrvv2brr 在环流区内，压强分布由伯努里方程式导出。环流区内半径为在环流区内，压强分布由伯努里方程式导出。环流区内半径为 的点和无穷远处的伯努里方程的点和无穷远处的伯努里方程: : rpvp22式中的式中的 即为即为 ， 为无穷远处的压强。将为无穷远处的压强。将 代入上式得代入上式得: : vvpv222282rpvpp 由上式可知，在涡束外部的势流区内，随着环流半径的减小，由上式可知，在涡束外部的势流区内，随着环流半径的减小，流速上升而压强降低；在

31、涡束边缘上，流速达该区的最高值，而流速上升而压强降低；在涡束边缘上，流速达该区的最高值，而压强则是该区的最低值，即压强则是该区的最低值，即 bbrv2222282bbbrpvpp涡束内部的速度分布为涡束内部的速度分布为: : 0rvrvv)(brr 由于涡束内部为有旋流动，伯努利积分常数随流线变化，故由于涡束内部为有旋流动，伯努利积分常数随流线变化，故其压强分布可由欧拉运动微分方程导出。对于平面定常流动，欧其压强分布可由欧拉运动微分方程导出。对于平面定常流动，欧拉运动微分方程为拉运动微分方程为: : ypyvvxvvxpyvvxvvyyyxxyxx11将涡核内任意点的速度投影到直角坐标上，则有

32、，代入上式得将涡核内任意点的速度投影到直角坐标上，则有，代入上式得: : xpx12ypy12将将 和和 分别乘以以上二式，相加后得分别乘以以上二式，相加后得: : dxdy)(1)(2dyypdxxpydyxdx或或 )2(222yxddp积分得积分得: : CvCrCyxp2222222121)(21在与环流区交界处，在与环流区交界处， ，代入上式，得积分，代入上式，得积分常数：常数： bbbbrvvpprr,222bbbvpvpC得涡核区的压强分布为得涡核区的压强分布为 ：2222222121bbrrpvvpp 由上式可知涡管中心的压强最低，其大小为由上式可知涡管中心的压强最低，其大小为

33、 ，涡核区边缘至涡核中心的压强差为涡核区边缘至涡核中心的压强差为2bcvppbbcbppvpp221 由以上讨论可知，涡核区和环流区的压强差相等，其数值由以上讨论可知，涡核区和环流区的压强差相等，其数值均为均为 。涡核区的压强比环流区的的低。在涡束内部。涡核区的压强比环流区的的低。在涡束内部, ,半径愈半径愈小小, ,压强愈低，沿径向存在较大的压强梯度，所以产生向涡核中压强愈低，沿径向存在较大的压强梯度，所以产生向涡核中心的抽吸作用，涡旋越强，抽吸作用越大。自然界中的龙卷风心的抽吸作用，涡旋越强，抽吸作用越大。自然界中的龙卷风和深水旋涡就具有这种流动特征，具有很大的破坏力。在工程和深水旋涡就具

34、有这种流动特征，具有很大的破坏力。在工程实际中有许多利用涡流流动特性装置，如锅炉中的旋风燃烧室、实际中有许多利用涡流流动特性装置，如锅炉中的旋风燃烧室、离心式除尘器、离心式超声波发生器、离心式泵和风机、离心离心式除尘器、离心式超声波发生器、离心式泵和风机、离心式分选机等。式分选机等。 2bv21 无旋条件是速度有势的充要条件。无旋必然有势，有势无旋条件是速度有势的充要条件。无旋必然有势，有势必须无旋。所以无旋流场又称为必须无旋。所以无旋流场又称为有势流场有势流场。速度势的存在与。速度势的存在与流体是否可压缩、流动是否定常无关。流体是否可压缩、流动是否定常无关。 在柱坐标系中在柱坐标系中rvrr

35、v1zvztzr,势函数的特性势函数的特性（1 1）有势流动中沿）有势流动中沿ABAB曲线的速度线积分等于终点曲线的速度线积分等于终点B B和起点和起点A A的速度势之差的速度势之差 ABBAzyBAxABdzzdyydxxdzvdyvdxv)()(上式表明，由于速度势是单值的，则该线积分与积分路径无关。上式表明，由于速度势是单值的，则该线积分与积分路径无关。 （2 2）速度沿封闭周线积分时，周线上的速度环量等于零速度沿封闭周线积分时，周线上的速度环量等于零 。0)ddzvdyvdxvzyx 式中式中 为拉普拉斯算子。满足拉普拉斯方程的函数为为拉普拉斯算子。满足拉普拉斯方程的函数为调和函数，故

36、速度势是调和函数，故速度势是调和函数调和函数。 2（2 2）势函数是调和函数。）势函数是调和函数。02222222zyx在柱坐标系中在柱坐标系中 在笛卡儿坐标系中在笛卡儿坐标系中 01122222222zrrrx平面不可压缩流体平面不可压缩流体 1、等流函数线为流线等流函数线为流线当当 常数时，常数时，0dyvdxvdyydxxdxyxyvvyxdxdy 即：即：在笛卡儿坐标系中，平面、不可压缩流体的连续性方程可写成在笛卡儿坐标系中，平面、不可压缩流体的连续性方程可写成0yvxvVyx若定义某一个函数若定义某一个函数),(yx，令，令yvxxvy2、流体通过两流线间单位高度的体积流量等于两条流线的流体通过两流线间单位高度的体积流量等于两条流线的流函数之差流函数之差。 不论是理想流体还是粘性流体，不论是有旋的还是无旋的不论是理想流体还是粘性流体，不论是有旋的还是无旋的流动，只要是不可压缩（或定常可压缩）流体的平面（或轴对流动，只要是不可