河口水流,波浪,潮流,泥沙,河床变形二维数学模型.pdfVIP

水 利 学 报 SHUILI XUEBAO 第5期 河口水流、波浪、潮流、泥沙、河床变形二维数学模型 - Uvf 马福喜 李文新 了 华北水利水电学院、郑州．450045) (北京大学，北京，1011871 摘 要 根据河口 波浪掀沙，潮流输沙”的泥沙运动机制，引入 “波浪幅射应力”反应波浪对承体和泥沙的 驱动力 考虑波浪对底部泥沙作用力、波浪对水流攘沙能力的影响及风场对水流和泥沙运动的影响 问时考虑 泥沙粒径大小、承体含盐度、含沙量及水流紊动程度等因素对絮凝沉速的影响 用ADI算法求解沿水探平均方 程，破开算子法解悬移质不平衡输沙方程及河床变形方程，引入通度概念处理岸滩不规则边界，建立河盯水流、 波浪、潮流、泥沙、河床变形二维数学模型 通过实侧资料及卫星摇感图像对所建数学模型进行捡骑，表明它 可以较好地预测河口泥沙的运动 关键词 波浪幅射应力，不规则边界 通度，·_· 一 翌坚 f可叫建v 中圉号 TV148 “波浪掀沙，潮流输沙”是波流共同作用下的河VI泥沙运动机制 波流共同作用使悬沙浓度场町 比无浪情况下增加40％～100％ ．因此，在河Ii：I地区，波浪和潮流对泥沙运动的控制作用具有同等 重要的意义．但波浪周期远小于潮夕周期，并且与潮流数学模型的计算时间步长相比也小的多．凼 此，本文在建立波流共同作用下的潮流、泥沙数学模型时，不考虑波浪的瞬时变化，只考虑波浪对潮 流和泥沙的综合效应 致于在建模过程中所要用到的波要素在空间的分布 (如波高在计算区域中的分 布等)，则由三维波浪传播的数学模型l2 计算结果给出．在研究波浪对潮流运动的综合效应时．引入 “波浪幅射应力”来反应波浪对水体的驱动力 在研究波浪对泥沙运动的综合效应时，除考虑 “波浪 幅射应力”对泥沙运动的影响外，还考虑波浪对底部泥沙作用力及波浪对水流挟沙能力的影响．此 外，在河VI地区，水面风应力也是影响潮流和沿岸流的重要因素 本文考虑了恒定风场对水流和泥沙 运动的影响．最后，一般情况下河VI地区的悬移质泥沙的颗粒粒径都较小 (0．01ram)， 沙絮凝 对沉速有影响．本文考虑了泥沙颗粒粒径大小、水体含盐度、含沙浓度及水流紊动程度等因素对絮凝 沉速的影响 对流动区域不规则固体边界，为避免通常的用台阶型折线近似曲线不规则边界，本文引 入通度概念_]J，用一段弦来逼近它 在计算时，用ADI算法求解沿水深平均的水流控制方程，及破 开算子法解悬移质不平衡输沙方程和河床变形方程．由此建立水流、波浪、潮流、泥沙、河床变形二 维数学模型 1 水流控制方程 具有不规则边界的三维水流连续方程为 ： v ·(Ou)=0 其中：0是流场区域的通度，它表明流场区域能够被流体通过的程度，在流体内部0=1．在流体外 部，即固体边界中0=0．这里采用通常的静压假定，将上述连续方程沿水深进行平均，得水深平均 连续方程： 研 所喜磊 ；器 ¨月 ’日收到1系华北水利水电学院重点学科、重点实验室基金资助项目1同时藏珠江水利委 会水刮科学 一 39 — + + ：0， (1) a 0 式中 对应于离散单元的面通度，为阿格中能够被流体通过的面积 (网格面积减去网格中固体或障 碍物的面积)与整个网格面积之比，定义在网格中