辽宁石油化工大学工程流体力学第八章.pptVIP

第二节 特征线法 第三节 管路水击计算 第四节 变水头泄流和排空及充满时间的计算 第八章 一维不稳定流动 第一节 一维不稳定流动基本方程 第二节 特征线法 第三节 管路水击计算 第四节 变水头泄流和排空及充满时间的计算 第一节 一维不稳定流动基本方程 一、连续性方程 如图8-1所示 流入流出控制体的净增质量 ： 内控制体中的质量变化为 ： 流入控制体的净增质量应等于控制体内质量随时间的变化 ： 得到连续性方程 图8-1 连续性方程控制体 或： 上式即为一维不稳定流动的连续性方程 确定密度的相对变化率： 压缩系数 ： 体积模量 ： 由微分定理和 可求得： 考虑压力波传递速度a 上式又可写成： 二、动量方程 如图8-2所示 压力P密度 可得： 由动量定理 得 设为水平管道则a＝0，有 图8-2 动量方程隔离体 由于 可得： 由于有： 又称波动方程 波动方程可写为： 特征线法是偏微分方程的一种数值解法，它将偏微分方程变换成特殊的全微分方程，然后对全微分方程积分，得到便于数值处理的有限差分方程。 一、 单变量的一阶拟线性偏微分方程 设方程为 为了把方程变为常微分方程，令特征方向 则方程可变为： （ 特征曲线 ） 沿此方向 的全微分为 是式 称为特征关系式或相容方程 可得： 二、两个变量的一阶拟线性偏微分方程 设方程组为 且： 沿特征方向可求得： 代入方程组可得： 偏导数的系数行列式 令特征方向为 则得到： 即一阶拟线性双曲型偏微分方程组将存在两族特征线方向 由此我们得四个常微分程 及 三、波动方程的特征线解法 将波动方程 与上一节的方程组比较得到： 所以 得到： 将系数H ，D ，E ，F 及G代入相容方程得： 因为 ＝ 所以波动方程经特征线解法变换的四个常微分方程： 一、管道水击现象 水击现象：在有压管路系统中，由于阀门突然关闭或开启（或其它原因）使管内流速发生突然变化，从而引起管内压力急剧交替升降的现象称为水击现象。 主要原因：液体的惯性和压缩性。 在不考虑能量损失情况下，阀门处的压力 随时间变化的波形如图8-5所示，并将周而复 始地重复下去。 实际上，在压力波的传播过程中有能量损失， 因此压力波将不断地衰减 。 水击波传播周期为： 图8-5 阀门处压力变化过程 水击现象的发生，对管路系统十分有害，因此必须防止它的发生成设法减弱它的影响，具体可采用以下措施： (1) 延长关闭阀门时间，使 ，这样可以避免直接水击的发生； (2) 缩短受水击影响的管长来降低水击压力； (3) 减小阀门关闭前的管道中流速以减小水击压力； (4) 在管路上安装气囊，以吸收压能，减小水击压力； (5) 水击压力与水击波传播速度有关。减小水击波速度就能减小水击压力，当液体已确定的情况下，为了减小水击压力，尽量选管径大、管壁薄而又富于弹性的管道。 另外还有一些减弱水击的方法，在实际中具体采用什么方法要根据具体情况而定。 二、水击压力计算 1 阀门的瞬时关闭 如图8-8所示 整理得 式中C= 称为水击波传播速度。 2 阀门瞬时部分关闭 3 阀门逐渐关闭 图8-8 近阀处液层 三、水击波的传播速度 对于一维管道流动，时间为 ，压力波传播距离为 ，同时，由于管材具有弹性，所以A为 ， 为 ，在 段内液体质量增量为： 由虎克定律 和薄壁筒拉力公式 可得到： 令 ，它相当于液体内的音速。则得： 水击波的传播速度计算公式 一、变水头泄流及排空所需时间的计算 讨论等截面积的容器中的流体泄流及排空问题 如图8-9所示 泄流流量： （ ） 设在微小时间段 内，液面下降了 的高度， 令容器的断面面积为 则有： 由于水头为 图8-9 不稳定泄流 则有：