基于和模式的大风作用下黄骅港波浪潮流泥沙的三维数值模拟.pdfVIP

泥 沙 研 究 年 月 第 期 基于 和 模式的大风作用下黄骅港 波浪 潮流 泥沙的三维数值模拟 赵 群 交通部天津水运工程科学研究所 天津 摘要 根据黄骅港外航道淤积问题的具体情况 本文采用三维细颗粒泥沙的 模式对黄骅港海域的 潮流 泥沙进行了模拟 数值计算结果与潮流为主要动力条件的现场实测结果进行了比较 二者吻合较好 在此基础上 本文采用 风浪模型和 对黄骅港骤淤进行了数值模拟 并与收集到的现场实测 资料进行比较 也得到较好的结果 数值模拟显示 波浪是造成黄骅港外航道泥沙淤积的主要动力因素 在以 波浪为主的动力条件下 含沙量在波浪破碎带达到最大 并且沿水深变化梯度较大 破碎带当地的泥沙跨过 航道落淤是黄骅港外航道淤积的主要原因 岸边泥沙随落潮流下泻不是黄骅港外航道大风骤淤的主要原因 关键词 黄骅港 三维潮流 细颗粒泥沙模型 中图分类号 文献标识码 文章编号 前言 黄骅港属于典型的粉沙质海岸 受大风影响淤积严重 其基本特征为大风骤淤 风 浪与泥沙淤积有 着密切的联系 波浪掀沙 潮流输沙 强浪弱流 波浪是影响港口外航道骤淤的主要因素 以波浪为主要动力的泥沙的运动和以潮流为主要动力的泥沙运动存在一定的区别 首先 由于 潮流的掀沙能力较弱 含沙量的垂线分布梯度不大 上下分布比较均匀 而在波浪作用下 尤其是崩破波 作用下 由于紊动能量基本集中在波谷以上 沿水深方向湍流的掺混不完全 含沙量受底面边界层的控 制 因此沿水深变化幅度较大 易形成底部高含沙水体 其次 由于含沙量的大小与波高成正比 与水深 成反比 因此 在潮流为主的动力条件下 含沙量在近岸浅水处较大 由岸边向深水处递减 而在波浪为 主的动力条件下 含沙量在破波带附近达到最大 考虑到以上情况 我们采用三维的细颗粒泥沙模型 和风浪模型 对黄骅港海域 的波浪 潮流 泥沙进行初步数学模拟 以期为进一步的航道淤积的研究奠定基础 三维潮流模型的控制方程及定解条件 基本控制方程 三维潮流模型以 模型为基础 该模型以速度 和压力 为待求变量 设 分别为 离岸方向和沿岸方向 自水底垂直向上 则雷诺平均的质量及动量守衡方程为 收稿日期 基金项目 国家自然科学基金 科技部项目 作者简介 赵群 女 浙江湖州人 博士 主要从事海岸动力学的研究工作 其中 为基准密度 为计算液体的密度 为垂直方向的湍流掺混系数 方程采用静压假设 上式中 大气压 为常数 为自由水面 的大小决定速度的垂直分布 由 湍 流模型来封闭 具体细节见文献 应当指出 对于三维水流模型 恰当的湍流模型是得到合理的流速 垂线分布的重要保证 实践证明 湍流模型适用范围广 精确度高 式 中 计算网格尺度以下的过程 分别由 和 来封闭 为 模型掺混系数 边界条件 在自由表面 运动学边界条件遵循 动力学边界条件要求边界层的切应力表达为速度梯度的函数 为自由表面摩擦力 在水底 底面剪应力为速度梯度函数 这里 为水底标高 为底面摩擦力 数值离散方法 模式采用有限差分及有限体积对控制方程及边界条件进行离散 为了既提高数值计算的精 度又减少对计算机内存的压力 在水平方向采用正交的曲线坐标