会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

水利工程论文-河势贴体河道平面二维正交网格生成方法的研究及应用.doc水利工程论文-河势贴体河道平面二维正交网格生成方法的研究及应用.doc -- 2 元

宽屏显示 收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文河势贴体河道平面二维正交网格生成方法的研究及应用摘要基于河势概念和Hermite三次插值函数,提出了河势贴体河道平面二维正交四边形网格的生成方法与边界层坐标系下水深平均流体力学控制方程和SIMPLER算法联合使用,建立了河道平面二维数学模型进行了葛洲坝枢纽至磨盘溪河段二维网格生成及水流泥沙数学模型的实例研究讨论了该网格生成方法的主要特征。关键词正交网络数学模型河势Hermite函数边界层SIMPLER算法1河势贴体网格河道平面二维数学模型网格生成方法研究中两个关键问题是1网格与河道拟合的贴体问题2二维网格、控制方程和数值方法三者之间的匹配问题。目前,常用的河道二维正交网格生成方法是边界拟合坐标系方法,即河道Thompson法1,它主要是通过物理平面天然河道平面与变换平面数模计算平面之间Poisson方程边值问题数值解实现二维正交网格的生成。式中x,y为物理平面网格坐标ξ,η为变换平面网格坐标D为物理平面区域。河道Thompson法存在的主要问题有1二维网络主要考虑与河道岸线的拟合,多数情况下与河道的河势或主流线之间的拟合不理想2复杂洲、滩及岸线河道边界、非恒定流动岸和数值求解引起的动边界等情况下,网格与河道岸线之间的拟合同样会出现实际偏离3控制方程变换引起的数值求解困难和数值诱发耗散问题4对于常见的宽、窄相间的河道平面形态,二维网格不均匀间距可能导致的数值计算精度问题5与河道一维断面数学模型的嵌套和联合存在接口困难等等。为了尽可能地避免河道Thompson法的上述问题和困难,本文提出了河势贴体河道平面二维正交网格生成方法。采用这一网格生成方法及相应的河道二维数学模型,不仅可以进行天然河道大多数边界条件下的水深平均平面二维水流、泥沙及温排和污排等计算分析及应用研究,而且还为长河段河道二维及一、二维嵌套数学模型的研究和应用提供了新思路。2网格生成方法张瑞瑾教授2认为河道水流的流态或河势具有很广的含义,一切标志河道水流总体倾向的现象,都被纳入这一概念之中。将河流动力学中基本和核心的河势概念,引入到河道二维网格生成方法及其河道二维数学模型的研究之中,具有更加坚实的理论依据和物理基础。河势贴体网格生成方法的总体思想为放弃传统的网格生成方法中严格要求网格与河道岸线相拟合的思路,采用Hermite三次插值函数3,生成河道沿程纵向与河势或主流线相拟合的河势拟合线曲线,并使得河宽方向横向网格线直线与网格控制断面相吻合,从而构造出平面二维正交四边形网格。由此生成的二维网格,一方面避开了河道复杂和变动的洲、滩及河岸岸线,另一方面体现了控制河道水流运动的河势概念。2.1Hermite插值函数Hermite插值函数不仅要求插值节点上的函数值相等,而且还要求节点处一阶甚至高阶导数相等。两个插值节点情况下的Hermite三次插值函数可表述如下3需要满足的节点条件为函数表达式为系数表达式为式中H3x,H3′x为Hermite三次插值函数及其导数α,β为插值系数x,y,m分别为插值节点的坐标、函数和一阶导数。2.2网格生成步骤河势贴体河道平面二维正交四边形网格生成方法包括如下三个步骤。2.2.1确定网格控制断面和节点选取研究河段的进出口断面、河段内水位/水文站点或测流断面、河势控制断面以及需要重点研究河段的控制断面等作为网格控制断面在所选取的网格控制断面上确定网格控制节点,这些节点可以任意选取在控制断面的左右岸、深泓点、主流点、中心点等处,所选择的网格控制节点即为数学上Hermite三次插值函数的计算节点。2.2.2生成河势拟合曲线利用上述Hermite三次插值函数,可以生成一条既通过网格控制节点,又垂直于网格控制断面的河道纵向网格控制曲线。通过多次调整网格控制断面和节点,使得所生成的网格控制曲线与研究河段的河势或主流线相拟合,将最终生成的河道纵向网格控制曲线确定为河势拟合曲线。这一步是河势贴体网格生成方法的关键和核心。2.2.3构造平面二维网格选取和调整纵向和横向网格间距,构造由平行于河势拟合曲线的曲线簇和垂直于河势拟合曲线的直线簇包括网格控制断面的河势贴体平面二维正交四边形网格。3河道二维数学模型河势贴体河道平面二维正交四边形网格生成方法与边界层坐标系下水深平均流体力学控制方程4以及合适的数值方法如SIMPLER算法5之间联合使用,可以建立河道二维数学模型。本文给出了河道二维水流泥沙数学模型的控制方程和数值方法。3.1二维水流泥沙控制方程采用边界层坐标系下简化的河道二维浅水控制方程和泥沙对流扩散方程4,6式中x,y,t为边界层坐标系下平面二维坐标和时间坐标R为x轴沿程曲率半径U,V为流速分量h,Z为水深和水位g,n,vt为重力加速度常数、河道糙率和紊动粘性系数S为含沙量S为水流挟沙力εs为泥沙扩散系数α为泥沙恢复饱和系数ω为泥沙沉速。上述控制方程中二阶导数项进行了适当的简化和合并当R→∞时,上述控制方程可转换为直角坐标系下水深平均平面二维水流泥沙控制方程有关定解条件和参、系数函数或方程如水流挟沙力公式等与直角坐标系下二维水沙数模基本一致7,8。3.2数值求解方法图1葛洲坝枢纽至磨盘溪河段河势图1997年9月河道地形图2葛洲坝枢纽至磨盘溪河段河道二维水沙数模河势贴体正交四边形网格图3葛洲坝枢纽至磨盘溪河段二维水沙数模计算的流速等值线二维水流控制方程的数值离散和迭代求解基于SIMPLER5算法,泥沙对流扩散方程的数值离散和迭代求解基于有限控制容积法5,离散方程迭代求解方法具体包括以三对角追赶法TDNA法为核心的逐行法和高斯塞德尔点迭代法,并配合块修正和欠松驰修正技术等。计算过程中采用了动边界模拟技术,具体处理措施包括每次迭代根据二维网格节点的计算水深值,均要判断和区分水域和陆域节点对于陆域节点采用边界隔墙法5处理,并让陆域节点始终保持一个较小的富余水深等。4应用实例河势贴体网格生成的主要目的是为河道平面二维数学模型的研究提供二维计算网格及离散节点。作为应用实例,本文给出了葛洲坝枢纽至磨盘溪近坝河段图1二维网格生成的具体过程及二维水流泥沙数学模型研究的部分成果。该项研究的主要目的是分析论证葛洲坝下游壅水工程措施如胭脂坝左汊布设潜坝的壅水效果及其对防洪和航运的影响,而采用壅水工程措施的主要目的是为了解决三峡工程运用初期枯水位时葛洲坝下游引航道通航水深不足问题。依据上述河势贴体网格生成方法及步骤,首先选取B0坝轴线,Y34庙咀,Y37宜昌,Y39,Y41宝塔河,Y44,Y46,Y49和Y50磨盘溪共9个河段河势控制或水文/水位断面作为网格控制断面,选取9个断面的中点作为网格控制节点然后,由Hermite三次插值函数生成了本河段的河势拟合曲线AB图1最后,通过确定纵向x,横向y的网格间距,构成本河段河道平面二维计算网格图2。最终生成的二维计算网格节点数为13741,河势方向x方向网格间距为100~200m,断面方向y方向网格间距为50m。采用生成的河势贴体河道平面二维正交网格,利用上述的河道二维水流泥沙数学模型,即可进行葛洲坝枢纽至磨盘溪近坝河段的二维水流泥沙数模计算分析研究。计算程序为自编HELIU11程序,有关研究内容及成果请详见长江科学院九五三峡工程泥沙问题研究子题分报告葛洲坝枢纽下游近坝段整治二维水流泥沙数学模型研究。本文仅给出了宜昌流量13500m3/s时该河段二维数模计算的流速等值线图3和1993年11月1日至1997年8月31日冲淤验证计算中的冲刷部位等值线图4。5讨论本文提出的河势贴体河道平面二维正交网格生成方法具有如下几个主要特点1引入河流动力学的河势概念,凸现了河道二维数学模型研究中的主要矛盾,避开了河道岸线拟合这一复杂但相对次要的矛盾2生成的二维正交网格与边界层坐标系下河道水深平均流体力学控制方程及相应合适的数值方法如SIMPLER算法三者之间匹配较好3网格保留了河道一维断面形式,为长河段河道二维以及河道一、二维嵌套数学模型的嵌套和联合提供了便利的二维网格和一、二维接口条件4生成的二维网格具有河宽方向上分布均匀,河势方向上可调整性大,操作上简便易行等优点5该网格生成方法不适用于强弯及鹅头型分汊等平面形态的少数河段。参考文献1吴江航,韩庆书.计算流体力学的理论、方法及应用M.北京科学出版社,1988.2张瑞瑾.河道水流运动的基本特性A.张瑞瑾论文集C.北京中国水利水电出版社,1996.3李庆扬,王能超,易大义.数值分析M.武汉华中工学院出版社,1982.4吴望一.流体力学M.北京北京大学出版社,1982.5美S.V.帕坦卡.传热与流体流动的数值计算M.北京科学出版社,1980.6DONGYAOHUA.ThegeneratingmethodofriverregimefittedorrivermainstreamfittedorthogonalgridanditsapplicationsA.LeeJ.H.,JayawardenaA.W.&WangZ.Y..Proceedingsofthe2ndInternationalSymposiumonEnvironmentalHydraulicsC.NetherlandsA.A.BalkemaRotterdamBrookfield,1999.7董耀华.河道温度及浓度场的平面二维数模研究J.长江科学院院报,1997,22226.8董耀华.空腔回流区水沙特性的计算分析J.泥沙研究,1999,23439.
编号:201312122211015984    大小:28.28KB    格式:DOC    上传时间:2013-12-12
  【编辑】
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:10次
奋斗不息上传于2013-12-12

官方联系方式

客服手机:13961746681   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关资源

相关资源

相关搜索

专业文献   建筑水利   精品文档   水利工程  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5