缓冲击波的计算与分析.pdfVIP

缓冲击波的计算与分析 1 2 1 娄妍 ，王如云 ，班长英 1 河海大学理学院（210098） 2 河海大学交通与海洋工程学院（210098） E-mail：louyan@ 摘 要：本文在缓冲击波理论 Von Karman 积分的基础上，通过分析给出了缓冲击波附体的 偏转角条件，指出了缓冲击波理论关系式的适用范围，并分析给出最大边墙偏转角的理论数 值。同时本文在缓冲击波微分式的基础上分析给出一种缓冲击波数值计算方法，并给出计算 波后物理量的计算公式，该方法计算简便，适于由计算机编程实现，较好的解决了 Von Karman 积分式不易应用以及简化式适用范围不广的不足。 关键词：急流，缓冲击波，偏转角条件 1 前 言 水工明槽急流中，当边墙向水流外部偏转时，由于水流失去边墙的依托，水面出现跌落， 产生扰动，形成缓冲击波。早在 20 世纪 30 年代 Von Karman 首先推导得到缓冲击波积分式[1,2] , 该积分式目前被各教科书普遍采用，但是由于该积分式较复杂，直接求解较困难，故常利用 辅助曲线图进行求解[3] ，但图形数据误差较大，且不便于应用。为解决 Karman 积分不易求 解的困难，90 年代刘韩生等人对缓冲击波计算公式进行了研究，得到简化积分式[4] 。简化 积分式虽然计算简便，易于实际应用，但是在适用范围和计算精度方面都存在一定问题。为 此本文分析给出一种缓冲击波的数值计算方法，并给出计算波后物理量的计算公式。该方法 计算简便易于编程实现。 我们知道缓冲击波理论是在水流附体 （水流附着边墙）情况下推导出来的，但当边墙向 外偏转角度大到一定程度时，在偏转点后面会形成水流脱离边墙状况[1] ，在那里会出现回流 漩涡，空蚀等复杂情况，从而致使缓冲击波基本理论在这些区域不成立。而迄今为止，对于 偏转角的限制只是在工程设计中给定一些大致的范围[5] ，还没有人给出保证缓冲击波基本 理论成立的最大边墙偏转角度。为此，本文在基本理论关系式的基础上，通过分析给出了缓 冲击波附体的偏转角条件，指出了最大边墙偏转角的理论数值。 2 缓冲击波关系式的简单分析 - 1 - 图 1 缓冲击波示意 缓冲击波是指边墙向水流外侧转折产生的水流扰动现象。图 1 所示为一明渠扩散段，河 道偏转起始点为A，在A点之后河道向水流外侧偏转，偏转角度记为θ，由于向外侧偏转，为 便于理解计算，偏转角取负值，即为−θ，称为边壁偏转角。假设来流为均匀急流，流速记 2 v 为v ，水深为h ，来流佛汝德数为Fr 1 ，且Fr 1 。受河道偏转影响，产生缓冲击 1 1 1 1 gh1 波。缓冲击波可以看作是由连续分布的干扰线所组成。干扰线形成的波高总和为∆h 。根据 再利用比能不变假定可以得到缓冲击波 两条流线之间线段的动量变化方程及经过一些变换， 微分式： dh 2(H 0 −h) h （1）