淹没丁坝冲刷坑内三维水流运动数值分析.pdfVIP

淹没丁坝冲刷坑内三维水流运动数值分析 周宜林许芝国唐洪武雷燕 (河海大学水利水电工程学院，南京210098) 摘要：淹没丁坝冲刷坑内三维水流泥沙运动是一种比较复杂的流动现象，而且复杂的水流现象与丁坝局部 冲刷密切相关。对冲刷坑内三维水流运动进行数值模拟，是丁坝设计、防护措施研究的重要手段之一。本 文应用已建立并经过检验的大涡模拟数学模型，对淹没丁坝冲刷坑内三维水流运动进行了研究。研究表明， 数学模型能够较好地重现丁坝冲刷坑内大尺度涡体的运动规律，对认识丁坝冲刷坑内水沙运动规律、指导 工程设计具有应用价值。 关键词：丁坝：大涡模拟；冲刷坑：局部冲刷 1前言 沿河(两岸)设置丁坝是一种传统的和比较有效的治河方法。丁坝主要用于护岸、导流、增加枯 水河道水深以及维护航道稳定等。丁坝在国内外河流中有较多的具体应用。尤其在黄河中游，为了保 护河岸和控制河势，沿河两岸布置了很多组丁坝群；在长江口几十公里的河段内，设置了丁坝，束水 攻沙，增大航运水深。在我国南方河流，珠江和赣江水系等，利用丁坝对航道进行整治。在美国密西 西比河等河流，采用丁坝对河道进行整治。在欧洲和日本的许多中、小河流上均设有丁坝。 丁坝设置后，丁坝附近水流运动特性将发生变化，由此必将引起丁坝附近的局部冲刷和淤积，影 响主流位置和流向、河道演变等。丁坝附近水流具有三维特征，副流明显存在，丁坝附近的冲刷和淤 积与水流的这种特征有关。因此，研究丁坝冲刷坑内三维水流运动，不仅可以搞清丁坝附近水流特征， 而且为河床演变分析(数值模拟)打下基础，还可以为丁坝设计和施工提供依据。 丁坝冲刷坑内水流运动是一种比较复杂的水流现象，一是流动现象复杂，存在各种不同尺度的 涡体；二是边界条件复杂，存在自由水面和不规则的固体边界。目前，对丁坝冲刷坑内水流的数值模 拟研究成果很少，本文即采用经过验证的大涡模拟数学模型，对丁坝冲刷坑内水流特性进行研究。 2数学模型 2．1基本方程 丁坝附近水流运动数学模型采用大涡模拟方程。大涡模拟是将物理量分解为大尺度和小尺度物理 大尺度涡体的运动由大涡模拟方程来描述和计算，小尺度涡体以小尺度应力的形式对大涡模拟方程做 贡献，并通过大尺度可解速度场建立的数学模型，封闭小尺度应力。大涡模拟方程，包括质量连续方 程和动量方程。 连续方程： 堕：o 孤2 动量方程： 鲁+掣强一詈丧嘶¨坶。】 (2) 作者简介：周宜林(1970一)：男，博士，河海大学水利水电工程学院副教授，主要从事泥沙运动基本理论、河流模拟 以及水环境模拟等方面研究。E-mall‘：zhouyilin@hhu．edu．cnTel：025ax：025 式中：S拥=吉0uz／孤。+au。／孤，)；U为滤波后 得到的大尺度速度；g为重力加速度： y为清水运 动粘滞系数；y，为涡粘性系数；Z=1，2，3；(单位 质量)水体压强P=p／p，P为压强，P为水的密 度。 涡粘性系数采用Smagorinsky模型…： V，=(Ga)lsl (3) 式中：C。为系数，对于明渠流，一般取C．s=0．1； 图1丁坝附近河床地形 lSl=42S胁Sfm；A=3VAxAyAz，Ax、缈、 分别为x、Y、Z方向的网格尺寸。 自由水面的确定是三维数学模型研究中一个 重要内容，自由水面边界条件包括运动边界条件 和动力学边界条件。自由水面位置的描述采用如 下的自由水面运动方程： ah：Ws_UW=——一1，4——sU丝一1，s丝 (4)(4) I—竺—o‰·翌卜—竺)一叫 氆dx 两 式中：h为自由水面高程；“5、1，5、W5分别为 自由水面纵向、横向、垂向速度。 (a)平面分区 自由水面处的压强等于大气压强，即： 只=乞 (5) 上式中：只为自由水面处压强，只为大气压强。 2．2数值方法 (1)大涡模拟方程： 本文采用有限差分计算方法‘21，计算网