清華水力第五章讲义.doc

第五章 有旋流动和有势流动 本章首先从运动学的角度对有旋流动的流场特性作进一步的讨论和分析，然后从动力学的角度介绍在质量力有势，流体为理想正压流体的条件下，有关涡通量的保持性定理。论述势流理论的基本内容，引出不可压流体平面流动的流函数概念，重点讨论不可压流体平面无旋流动的速度势函数与流函数的关系以及求解势流问题的奇点叠加方法。 §5—1 有旋流动 有旋流动与有势流动的判别就在于流速场的旋度是否为零。对于有旋流动，将流速场的旋度称为涡量，它是流体微团旋转角速度矢量的两倍，. 正象流线是流速场的矢量线一样，定义涡线是涡量场的矢量线。涡线的微分方程为： ， 即 . 对应于流速场中流管、流量的概念可以建立涡量场中的涡管、涡通量概念。涡管的涡通量又称为涡管强度I . 斯托克斯公式 表明了速度环量与涡通量之间的关系。其中为曲面A的法向，A周界L的正向与成右手系。 由于 ，所以涡量场是无源（管形）场。这表明在同一时刻，穿过同一涡管的各断面的涡通量都是相同的。换句话说，同一时刻，一根涡管对应一个涡管强度。这是个纯运动学范畴的定理。表明涡管不能在流体中产生与消失，要么成环形，要么两端位于流场的自由面或固体边界。 封闭流体线上的速度环量对于时间的变化率等于此封闭流体线上的加速度环量。记 ，则 . 注意这里的L是由确定的流体质点组成的一条封闭线，是一个系统，在流动中会改变位置和形状。相应的速度环量对时间的变化率指的是全导数。 现简要地证明上面的结论，为使符号意义更明确，把对空间微分符号记为，符号d专门用于对时间的微分。 最后一个等号是因为速度为单值函数。 上述运动学规律若加上下面的动力学条件，即质量力有势，流体为理想不可压缩流体，那么根据欧拉方程易知加速度是有势的，封闭流体线上的加速度环量必为零，速度环量将不随时间变化，. 称为开尔文定理。 由开尔文定理推知：某时刻组成涡管的流体质点将永远组成涡管，且其涡管强度在流动中保持不变。称为亥姆霍兹定理，其成立的条件应为开尔文定理成立所需的条件。 开尔文定理说明，若质量力有势，流体为理想不可压缩流体，那么涡通量不会产生，初始时刻为无旋的流动将永远保持为无旋而有旋流动的涡通量则有保持性，既不会消失，也不会扩散。开尔文定理也反过来说明了之所以在实际流体的运动中会有旋涡的产生、发展和消失，以及涡量在流场中的扩散现象，粘性的存在应该是最重要的因素。 §5—2 有势流动 由开尔文定理可知，理想正压流体从静止或无旋状态开始的流动将保持为有势流动。所以有势流动往往是以理想流体为前提条件的。 一. 无旋流动的速度势函数 无旋流动 ,根据斯托克斯公式，沿流场中任意一条封闭曲线L的速度环量为零，即 . 自起点M0至终点M的速度环量与路径无关，在起点固定的条件下，是终点位置的函数，所以可定义 为速度势函数. 不难证明：，，. 这表明 ， 是全微分。 根据定义，速度势函数可相差一任意常数而不会影响对流场的描述。极坐标中速度势函数的微分为 . 在确认流动的速度场无旋，即 的前提下，可由速度场根据定义 求取势函数，选择由几段坐标线组成的路径以使计算简单，千万不可误认为做三个独立的不定积分。也可以用寻找全微分 的办法求速度势函数。 不可压流体无旋流动的速度势函数满足拉普拉斯方程：. 满足拉普拉斯方程的函数称为调和函数。 二. 不可压流体平面流动的流函数 不可压流体平面流动的连续方程为 ，它可改写为 ，这恰好说明矢量场 是无旋场，必有相应的势函数，我们就定义其势函数为原流速场的流函数： . 将平面上一段有向微元弧长 顺时针转900，方向为之法向，大小为，可记为. 根 据流函数定义知 ，这说明流函数的微分为穿过微元弧长的流量，所以把 称为流函数。 表示穿过M0至M连线的流量，它与连线路径无关，在起点M0确定的情况下是终点M的坐标的函数。这对于不可压流体的平面流动是容易理解的，而三维流动就得不到这样的结论。根据定义确定流函数时选取不同的起点M0，流函数将相差一个常数，但同样不会影响对流场的描述。 （常数）是不可压流体平面流动的流线方程。两点流函数的差表示穿过两点间任意连线的流量，如图中所示，若 ，表示有流量自M1M2连线左侧流进右侧，由此可确定流动方向。 画出穿过微元弧长的流量示意图，可以帮助记忆流函数定义。如在直角坐标系中，容易写出，所以. 而在极坐标中，则有，所以 . 如果不可压流体平面流动是无旋的，那么 ，说明流函数满足拉普拉斯方程，是调和函数。 流函数的概念本与流动是否无旋无关，在这里引出，是为了下面建立不可压流体平面无旋流动速度势函数与流函数关系的需要。 三. 不可压流体平面无旋流动的速度势函数与流函数的关系 不可压流体平面无旋流动既有速度势函数又有流函数，它们都满足拉普拉斯方程，都是调和函数