原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
长江汉口河段平面二维水沙数学模型初步研究.doc
------------《水资源研究》第24卷 第3期(总第88期)2003年9月--------------------
长江汉口河段平面二维水沙数学模型初步研究
张 革 联
（长江水利委员会 长江中游水文水资源勘测局，湖北 武汉 430012）
摘 要： 在长江汉口河段二维泥沙数学模型的建立与求解过程中，采用了曲线四边形9点等参单元法，任意剖分计算区域，并构造出形函数将水流连续方程、运动方程、泥沙扩散方程和河床变形方程离散，根据来水来沙过程并利用边界条件，模拟了较长时间的河床变形。
关键词：形函数； 曲线四边形； 坐标变换； 边界条件
1 基本方程
水流连续方程：
(1)
水流运动方程：
(2)
泥沙扩散方程：
(3)
河床变形方程：
(4)
水流挟沙率公式：
(5)
合流速公式：
(6)
式中h为水深；u、v为横向及纵向水流速度；Z为水位；ZB为河床高
程；n为满宁二维糙率；νt为涡流系数；εS为泥沙扩散系数；α为泥沙综合系数；ωi为分组沉速；ρ′为泥沙干容重。
2 基本方程的求解
2.1 单元插值形函数
采用曲线四边形9点等参单元，任意剖分计算区域，并构造出形函数为完全双二次插值函数：
(7)
故任意函数可表示为i为主元结点编号。
基本单元可用图1、图2表示，经坐标变换后，即可得到对应的矩形单元，变换过程的Jacobian矩阵为：
(8)
2.2 基本方程的离散
根据微分法则，对于任意函数的一阶偏导数有：
(9)
于是
(10)
式中J为坐标变换的Jacobian阵；J-1为其逆阵。
令
｜J｜为Jacobian行列式，那么，(10)式可写为：
(11)
再求二阶偏导数：
求得微商后，即可将原方程离散。
2.3 Neumann边界条件的处理
(1) 对于流速计算，以边界条件
以及作为Neumann条件进行离散，其中Z为水位；为穿过边界结点的法线，原Neumann条件即可表示为：
（13）
式中α为法线方向角。
(2) 对含沙量计算同样有第2类边界条件，即：
(14)
3 数学模型的实际运用
3.1 计算河段范围及河型特点
本文计算河段为汉江河口区域下太古断面以上约10 km范围，属汇流型微弯河型。
3.2 计算网格的划分及边界条件处理
计算前将整个区域划分成40×100个网格，计算流场时以下游出口断面水位过程Z=Z(t)作为出口Dirichlet边界条件，进口断面以流量过程Q=Q(t)为进口边界条件。
3.3 模型参数
3.3.1 阻力系数
二维阻力公式为：
(15)
式中n为糙率；J为比降，以“0”为脚标的为垂线均值；f(η)函数是通过冲积平原的大量实测资料整理的一种相对河面函数[1]。
3.3.2 水流挟沙力
垂线单宽挟沙率公式[1]为：
(16)
3.3.3 挟沙力、河床质级配[1]
在冲淤平衡下，考虑垂线含沙量Rouse公式所建立的挟沙率级配与河底之河床质级配之间关系为：
(17)
(18)
(19)
式中ΔPi、ΔPbi分别为第i粒径组的挟沙率级配及河床质级配；ωi为第i粒径组沉速 ；为摩阻流速；k为卡门常数；σv为水流垂向脉动强度系数；N为泥沙分组粗径数。
3.3.4 悬移质中床沙质粒径范围的确定[1]
关于悬移质中床沙质粒径范围的确定一般用悬浮指标：
(20)
式中为摩阻流速；β为系数。
床沙质对应之最小及最大沉速为：
(21)
式中ωmin、ωmax用沉速公式计算，鉴于冲积平原河道的河床质与悬移质中床沙质粒径范围基本相同，可考虑忽略推移质影响，实行较稳定的几组粒径。
3.3.5 涡流紊动系数 νt的确定
采用νt与水力因素之间的下列关系来计算。其中c为与本河段有关的经验系数。
3.4 模型的验证及展望
通过运用1986年6月至1993年9月的实际来水来沙过程及初始地形资料为起始计算，对河床地形进行为时7 a的二维河床细部模拟，计算结果基本满意，流速计算效果见表1。（略）
今后此模型有待进一步通过更多的实际工程问题得到完善与提高，特别对于网格划分，粒径范围的确定及阻力系数和紊动系数的确定有待进一步地准确，适合当地河道情况。
参考文献：（略）
文档评论（0）