挑流水舌空中形态及含水浓度计算方法..docVIP

挑流水舌空中形态与含水浓度的数值解法 柳海涛，孙双科，王晓松，夏庆福 (中国水利水电科学研究院水力学研究所北京100038) 摘要：本文基于自由水面掺气理论，建立了计算掺气水舌的空中形态、流速和含水浓度的数学模型，并运用东江水电站的原型观测资料对模型进行验证。结果表明，该模型能够有效地解决复杂挑坎体型下的掺气水舌问题，计算成果也可用于下游消能和雾化问题的研究。 关键词：掺气水舌；数学模型；原型观测；含水浓度 The numerical solution of the surface figure and water concentration of aerated jets Liu haitao, Sun shuangke, Wang xiaosong, Xia qingfu (China Institute of Water Resources and Hydropower Research, Beijing, 100038) Abstract: In this paper, a numerical model was developed for calculating the surface figure, velocity and water concentration distribution of aerated jet in atmosphere, based on the theory of free surface air entrainment. The good agreement between the calculated results and the prototype observation data from Dong-jiang project proved that, the model is very effective to solve the problems of aerated jets from the complex-type flip buckets. The result is also useful to the study of the downstream energy dissipation and jet flow atomization. Key words: aerated jet; numerical model; prototype observation; water concentration 挑流消能工自1933年被首次应用于西班牙里科拜约重力拱坝的溢洪道出口以来，已有70多年的发展历史。由于其结构简单，不需要在下游修建河床防护工程，在国内外得到了广泛的应用。近20年来，随着我国水利水电建设技术的不断发展，出现了一批高水头大流量的挑流泄洪消能工。同以往的中小工程不同，为提高空中消能率和减轻下游冲刷，这些大型消能工采用了各种扩散性较好的异型挑坎，使其在泄洪时形成的高速挑射水流，能够充分掺气、裂散，甚至翻卷。这些都明显改变了水舌的空中及入水形态，并且引起的大范围雾化降雨，对当地的生产和生活也产生了不利影响。因此，对于高坝挑流水舌的研究，目前已不再局限于对坝面掺气和水舌挑距的初步估算，还需要对水舌断面形态与含水浓度的沿程变化进行深入研究，为下游冲刷及雾化的分析评估提供参考。在这方面，已有不少最新的成果[1~6]。本文在前人研究的基础上，将通过对水舌表面掺气机理进行理论分析，建立计算水舌表面特性沿程变化的数值模型，该模型除了考虑重力、空气阻力、浮力等因素的影响外，还考虑水舌掺气引起的体积膨胀。通过计算，可得到掺气水舌的空中轮廓、流速以及含水浓度的沿程变化，并采用相关的实测资料对模型参数进行率定和验证。 1. 挑流水舌运动微分方程 水量守恒方程 在水舌沿流线方向取一微元微元的侧边界无水体穿过，： (1) 其中，为水的密度为微元的时均流速；为断面体积含水浓度A为微元的横断面面积，为微元外湿周，为厚度。由于为常数，有： 水气连续方程 平均运动速度的函数： (3) 其中，右端第一项为气泡通过水流紊动获得的运动速度，系数与水舌脱离固壁前的水力条件有关，可视为常数；右端第二项表示气泡用以克服表面张力、水的粘滞性以及内外压差等所需的最小流速。 考虑到实际工程中，由于水舌充分紊动，上式中右端项略去不计，故又可以简化为： ]，上式中的系数称为卷吸系数，为时均流速。 由此可知，单位时间，单位面积上从外表面进入微元的气泡质量为，为微元外边界上的含气浓度，可以认为，为外边界速度；内边界上流出的气体质量为，分别为微元内边界处的含气浓度和速度，。 微元横向与纵向边界