- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
泥 沙 研 究
年 月 第 期
一般曲线坐标系平面二维水沙数学模型研究与应用
谢作涛 张小峰 袁晶 许全喜 杨芳丽
武汉大学水资源与水电工程科学国家重点实验室 湖北武汉 长江水利委员会水文局 湖北武汉
摘要 对笛卡尔坐标系下控制方程进行坐标变换 采用控制体积法与 算法对方程进行离散和求解 建
立了一般曲线坐标系平面二维全沙模型 模型控制方程离散时物理量布置采用同位网格 避免了采用交错网
格时所需的大量插值计算 提高了计算效率 模型较全面地考虑了非均匀悬移质及推移质运动 具有较好的
普遍适用性 通过长江沙市河段 年至 年实测水沙资料对模型进行了较系统的验证 计算结果表明
模型能够较合理反映天然河道水沙运动及河床变形规律
关键词 一般曲线坐标系 同位网格 控制体积法 算法
中图分类号 文献标识码 文章编号
前言
天然河道一般边界曲折 地形复杂 对于复杂河段的平面二维水沙运动数值模拟 多采用基于正交
曲线网格的正交曲线变换方法 但在工程实际运用中 有时很难形成完全正交的曲线网格 采用正交曲
线坐标模型 易引入近似误差 进而影响模型计算精度 与正交曲线变换相比 一般曲线变换不受计算
网格必须严格保证正交的限制 网格生成也较灵活 本文从笛卡尔坐标水沙控制方程出发 对控制方程
进行一般曲线变换 采用控制体积法与 算法对方程进行离散和求解 建立了一般曲线坐标系平
面二维水沙模型 并采用天然河道实测水沙资料对模型进行了较系统的检验
模型控制方程及坐标变换
笛卡尔坐标系控制方程
笛卡尔坐标系下沿水深平均的平面二维水沙运动控制方程可表示为
水流连续方程
方向水流运动方程 ?
方向水流运动方程 ?
泥沙连续性方程
河床变形方程
收稿日期
基金项目 国家重大基础研究计划 项目 国家自然科学基金
作者简介 谢作涛 男 湖北洪湖人 博士研究生 研究方向为河流海岸数值模拟
其中 为水深 和 为 和 方向的流速 ? 为水位 为 糙率
系数 为紊动粘性系数 为水体密度 ? 为含沙量 ? 为挟沙力 ? 为泥沙颗
粒沉速 ? 为泥沙干密度 ? 为恢复饱和系数 为泥沙扩散系数 分别对应分
组泥沙颗粒沉速 含沙量和挟沙力 分别为 和 方向的推移质输移率 其表达式为
一般曲线坐标系控制方程
一般曲线坐标下 任一物理量 的偏导数与笛卡尔坐标系偏导数对应的转化关系为
记曲线坐标中流速在 和 方向的分量为 则有
同时利用关系式 ? ? ? ? 再通过链式求导法
则 可得到一般曲线坐标的平面二维水沙控制方程
水流连续方程
方向水流运动方程
?
方向水流运动方程
?
泥沙连续性方程
河床变形方程
其中 表示偏导数 如 其它量意
义与笛卡尔坐标系相同
控制方程离散
目前求解不可压缩水流流动采用较多的是
文档评论(0)