三维水沙数值模拟中紊流随机理论应用初探.pdfVIP

泥 沙 研 究 年 月 第 期 三维水沙数值模拟中紊流随机理论应用初探 张新周 窦希萍 王惠民 许 慧 河海大学 江苏南京 南京水利科学研究院 江苏南京 摘要 在平面上采用一般曲线坐标 在垂向上采用 坐标变换 用模式分裂法求解自由面水位 采用窦国仁紊 流随机理论 构建紊流随机模型下的水流泥沙三维数学模型 采用有限体积法离散方程 并针对水工构筑物 附近水流运动特点 采用动水压力校正法求解耦合的压力和速度场 在泥沙起动流速里考虑了垂向流速的影 响 把计算结果与试验结果及基于 假设的紊流模型计算结果进行了对比 表明在同等计算条件下 基于紊流随机理论的紊流模型有更精确的模拟效果 能较好的反映三维水流泥沙的运动和输移情况 关键词 动水压力校正 紊流随机理论 模式分裂法 算法 假设 中图分类号 文献标识码 文章编号 近年来 有不少学者发表了关于三维水沙数学模型研究的文章 采用相对分层 模型计算了海湾和河网的三维水流 和 通过绝对分层的方法求解计算了河口的三维水流 窦希滨 应用紊流随机理论模型进行了桥墩局部冲刷的三维水流 泥沙数值模拟计算 叶坚 平面采 用一般曲线 垂向采用 坐标变换 建立了同位网格下动量差值的三维水流模型 很好的模拟了弯道和 分汊水流的三维运动 在国内 赵士清 提出了一个三维数学模型并应用于长江口潮流的计算 该模型 在垂线方向上将水体划分为若干层 从而将三维问题简化为求解一系列二维问题 既节省了存储量又节 省了工作量 李　来 窦希萍等 在正交曲线坐标系下 采用坐标变换 有限差分 交替隐式求解 动边界 处理等方法 建立了边界拟合网格的三维潮流数学模型 成功地应用在长江口潮流场的计算中 李国 斌 采用紊流双方程模型闭合雷诺方程组 用控制体积法离散三维偏微分方程 并由 计算程 式求解淹没丁坝流场不同深度位置的流速 紊动能和紊动粘滞系数的分布规律 夏云峰 在平面采用一 般曲线 垂向采用 坐标变换 建立了非交错网格变量布置的水深 流速迭代计算模式 并由水深校正 法求解水体自由表面和流速场耦合问题 陆永军等 根据紊流随机理论 导出了各向异性紊流的雷诺 应力的数值格式 并将精细壁函数应用于边壁处理 建立了三维紊流数学模型 并应用于三峡坝区冲淤 的模拟 童朝峰 使用平面一般曲线垂向 坐标变换 建立了交错网格变量布置的压力 流速校正计 算模式 并成功的应用于分汊口水流的计算 本文在已建立的二维水沙数学模型基础上 采用紊流随机理论建立紊流模型 考虑到水工构筑物附 近流速场和水位场发生剧烈变化 压力场不再按照静压假定分布 在求解中采用动水压力修正法求解流 速和动水压力的耦合方程 以得到水工构筑物附近更精确的水位场和流速场 在此基础上建立泥沙输移 的三维数学模型 对三维水流和泥沙运动进行细致研究 收稿日期 作者简介 张新周 男 河南博爱人 博士研究生 主要从事河流海岸动力学及泥沙运动等方面的研究 数学模型 水流方程 笛卡尔坐标系下的水流方程如下 连续方程 动量方程 式中 为时均流速在 方向上的分量 为质量力 为运动粘滞系数 为压强 为流体密度 为雷诺应力 紊流模型 由于方程组 和 式是不封闭的 就需要对雷诺应力项增加新的补充关系式 以使方程组封闭 从而发展出了许多紊流模型 着重介绍如下 基于 假设的双方程紊流模型 在 假设下 雷诺应力可写为 式中 为 符号 为紊动粘滞系数 可由紊动动能 和紊动耗散率 计算得出 式中 为经验常数 取值 基于紊流随机理论的紊流模型 在 假设下 对于不同方向的雷诺应力 紊动粘滞系数取值相同 事实上 在复杂的流动 中 紊动粘滞系数在各向的分布是不相同的