可冲刷坡面滚波数值模拟研究.docVIP

可冲刷坡面滚波数值模拟研究 李侃禹1,曹志先1[基金项目: 欧盟第六框架计划资助项目037024; 国家自然科学基金重点项 作者简介: 李侃禹1983-, 女, 研究生; 武汉, 武汉大学水资源与水电工程科学国家重点实验室430072 *联系人/通讯作者, e-mail: zxcao@]*,刘青泉2 1 武汉大学水资源与水电工程科学国家重点实验室,湖北 武汉 430072;2 中国科学院力学研究所,北京100080) 摘要: 本文应用一维水沙耦合数学模型研究了可冲刷坡面上滚波的水动力学特性。模型的基本控制方程采用完整的、基于守恒定律的一维浅水动力学方程,运用能够捕捉激波和泥沙运动不连续性的WAF TVD二阶数值格式离散控制方程。通过复演定床滚波的运动特点,对模型进行了验证。对可冲刷坡面及有下渗情况下的滚波进行了数值模拟,结果表明床面形态对滚波水动力学特性的影响显著,下渗仅影响波幅的大小。 关键词: 滚波; 可冲刷坡面; 数值模拟; 下渗; 弗汝德数 1 引言 滚波是明渠水流自由面失稳后发育的一系列波动现象,它对坡面土壤侵蚀有着重要的影响。在坡面比降较陡时,坡面水流基本上以急流、紊流状态顺坡下泻,坡面流边壁状况的不规则变化或雨滴打击扰动等因素对水流所产生的任何微小扰动都会导致坡面水流以一种特殊的流动形式?滚波?沿坡向下运动,较大的波追及较小的波聚叠而成更大的波并继续向前传播。早在1934年Cronish[1]就推测出了滚波,Thomas[2]最早分析了大振幅的滚波,Dressler[3]于1949年得到了滚波演化方程的间断性解析解,认为谢才公式中的无量纲阻力系数需要达到床面比降的1/4才能形成滚波。Brock[4]于1969年所做的水槽试验成为这方面的经典工作。Brock[4]发现,由均匀流失稳后形成的滚波不是周期性的,波形、波速在不断发生变化;滚波的演化由3个阶段组成:①初始的小振幅阶段,波周期基本不变;②过渡阶段,波周期开始增加和出现波的聚合过程;③最后阶段,此时发生显著的聚合现象,可以是小波聚合大波,也可以是大波吞并小波;同时,波周期也随着波高的增加而增加。最近,Liu等[5]、Balmforth和Mandre[6]都对滚波的结构进行了深入的理论分析,Zanuttigh和Lamberti[7]利用基于HLL算子的WAF TVD数值方法求解一维浅水动力学方程,较好地模拟了Brock[4]试验中滚波演化的3个阶段。我们注意到:现有关于滚波的理论、实验和数值研究局限于定床条件。然而,滚波可能发生在较易冲刷的坡面上如干旱和半干旱地区暴雨产流产沙过程,水流可能挟带一定数量的泥沙,且在运动过程中伴随着下渗。但迄今尚未见有关可冲刷坡面上滚波的研究报导。本文应用一维水沙耦合数学模型Cao等 [8] 对可冲刷坡面及有下渗情况下的滚波形成与演化进行数值研究。 2 数学模型 2.1基本控制方程 本文在Cao等 [8]的一维浅水水沙耦合数学模型基本控制方程的基础上,考虑了下渗及表面张力。控制方程可以整理成如下守恒形式: 1 ,2 3 4 式中U为守恒向量;F为通量向量;S为源项向量;t为时间;x为空间坐标;h为单元水深;u为断面平均流速;z为床面高程;c为断面平均含沙量;g为重力加速度;为阻力坡度;p为床面泥沙孔隙率;E,D为水体底部与床面边界上的泥沙上扬通量和沉降通量;为水沙混合密度;为床面泥沙饱和湿密度;,为水和泥沙的密度;为表面张力[9];为表面张力系数,其单位为N/m;R为曲率半径;为下渗率。 2.2 封闭模式 应用谢才系数C计算阻力坡度: 5 泥沙沉降通量按Richardson and Zaki[10]公式计算:6 式中,为单颗粒泥沙在清水中的沉降速度;d为泥沙颗粒粒径;为近床体积含沙量,可以根据平均含沙量计算,即,;指数,,为清水运动粘性系数,。 对于水深较小及床面比降较陡的坡面流,泥沙上扬通量的定量化研究远未充分。这里,选择Sander[11]公式,泥沙上扬通量可按下式计算: 7 式中,。由于目前对坡面流挟沙能力的研究还不成熟,故本文对上扬通量进行了不同系数的修正,分别取值为0.005和0.02。 2.3 数值格式 本文运用够捕捉激波和泥沙运动不连续性的WAF TVD二阶数值格式,求解完整的、基于守恒定律的一维浅水水动力学方程。 应用算子分裂法[12]离散方程1: 8 式中为时间步长;为空间步长;i为空间节点号;k为时间节点号;为数值通量,详尽算法参考Cao等[8]。 床面变形方程4离散为: 9 3 模型验证以下验证本文的一维水沙耦合数学模型能否准确、完整地模拟滚波特性。 3.1水流失稳条件验证 计算条件为:矩形水槽长1000 m,渠道的谢才系数为C5