一般曲线坐标系平面二维水沙数学模型的研究与应用.pdfVIP

泥 沙 研 究 年 月 第 期 一般曲线坐标系平面二维水沙数学模型研究与应用 谢作涛 张小峰 袁晶 许全喜 杨芳丽 武汉大学水资源与水电工程科学国家重点实验室 湖北武汉 长江水利委员会水文局 湖北武汉 摘要 对笛卡尔坐标系下控制方程进行坐标变换 采用控制体积法与 算法对方程进行离散和求解 建 立了一般曲线坐标系平面二维全沙模型 模型控制方程离散时物理量布置采用同位网格 避免了采用交错网 格时所需的大量插值计算 提高了计算效率 模型较全面地考虑了非均匀悬移质及推移质运动 具有较好的 普遍适用性 通过长江沙市河段 年至 年实测水沙资料对模型进行了较系统的验证 计算结果表明 模型能够较合理反映天然河道水沙运动及河床变形规律 关键词 一般曲线坐标系 同位网格 控制体积法 算法 中图分类号 文献标识码 文章编号 前言 天然河道一般边界曲折 地形复杂 对于复杂河段的平面二维水沙运动数值模拟 多采用基于正交 曲线网格的正交曲线变换方法 但在工程实际运用中 有时很难形成完全正交的曲线网格 采用正交曲 线坐标模型 易引入近似误差 进而影响模型计算精度 与正交曲线变换相比 一般曲线变换不受计算 网格必须严格保证正交的限制 网格生成也较灵活 本文从笛卡尔坐标水沙控制方程出发 对控制方程 进行一般曲线变换 采用控制体积法与 算法对方程进行离散和求解 建立了一般曲线坐标系平 面二维水沙模型 并采用天然河道实测水沙资料对模型进行了较系统的检验 模型控制方程及坐标变换 笛卡尔坐标系控制方程 笛卡尔坐标系下沿水深平均的平面二维水沙运动控制方程可表示为 水流连续方程 方向水流运动方程 ? 方向水流运动方程 ? 泥沙连续性方程 河床变形方程 收稿日期 基金项目 国家重大基础研究计划 项目 国家自然科学基金 作者简介 谢作涛 男 湖北洪湖人 博士研究生 研究方向为河流海岸数值模拟 其中 为水深 和 为 和 方向的流速 ? 为水位 为 糙率 系数 为紊动粘性系数 为水体密度 ? 为含沙量 ? 为挟沙力 ? 为泥沙颗 粒沉速 ? 为泥沙干密度 ? 为恢复饱和系数 为泥沙扩散系数 分别对应分 组泥沙颗粒沉速 含沙量和挟沙力 分别为 和 方向的推移质输移率 其表达式为 一般曲线坐标系控制方程 一般曲线坐标下 任一物理量 的偏导数与笛卡尔坐标系偏导数对应的转化关系为 记曲线坐标中流速在 和 方向的分量为 则有 同时利用关系式 ? ? ? ? 再通过链式求导法 则 可得到一般曲线坐标的平面二维水沙控制方程 水流连续方程 方向水流运动方程 ? 方向水流运动方程 ? 泥沙连续性方程 河床变形方程 其中 表示偏导数 如 其它量意 义与笛卡尔坐标系相同 控制方程离散 目前求解不可压缩水流流动采用较多的是