原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
1 基本理论和方法
1.1 模型的基本方程
1.1.1 水流模型方程
水流连续方程：
(1) 水流运动方程： (2) (3) 其中：x、y为空间坐标，t为时间,h为水深，z为水位z=z0+h,z0为床面高程，u,v为流速在x、y方向上的分量，f为柯氏力系数f=2ωsinφ,ω为地球自转速度，φ为纬度，Cs为谢才系数，曼宁公式Cs=h1/6/n,n为糙率，g为重力加速度，λ为涡动粘滞系数，，拉普拉斯算子。 (4) ??? 河床冲淤演变过程： (5) D=Φ1Cω* (6) 其中：Φ1为校正值。
1.2.2 泥沙起悬量P的确定
??? 泥沙起悬量P与饱和挟沙状态下的含沙量C*和泥沙沉速ω*有关，可表示为：
P=Φ2C*ω* (7) ??? C*通过Lane—Kalinske公式确定为：
C*=5.55I0[0.5(u*T/ω)exp(-(ω*/u*T)2)]1.61 (8) 式中u*T为潮流与波浪合成摩阻流速，
u2*T=(u2+v2)/C2s+fwu20/4 (9) u0=πH/[Tsh(2πh/L)]，波浪底部最大轨迹流速 ??? ab=H/2[sh(2πh/L)]-1，近底波浪质点振幅，r*=120d ??? ,泥沙的平均沉速，ω1)与最小(d2)颗粒的沉速；I0为床面活动泥沙含量。 ??? 这样，将床面层附近泥沙起悬量和沉降量区分开来处理的方法，具有明确的物理意义，可以清楚地反映出悬移泥沙与床面层内运动泥沙及床面活动泥沙发生相互交换过程中矛盾的两个方面。在波浪掀沙和潮流输沙过程中，平衡只是相对的、暂时的，而不平衡则是绝对的、长期的。泥沙起悬量作为一个起悬函数(Pickup Function)的形式出现，有助于反映泥沙在床面层附近的交换机理，同时采用Lane-Kalinske公式，其优越之处在该公式在建立时，充分考虑了底部水流湍流脉动特性，反映了一定的力学机理，从目前流体力学研究来看，所不完善的是尚未区分大尺度拟序结构和小尺度随机脉动的作用，这方面我们正进行深入研究，但作为工程应用，该公式已经过了许多实验室和现场天然实测资料的验证，证明是可靠的[1,2]。总之，这样选取，有助于逐步弄清泥沙在床面层附近的交换机理，也有助于从力学机理上反映掀沙和输沙过程的不平衡特性。
1.3 模型方程的数值求解方法
1.3.1 水流方程的数值解
??? (1)杂交有限元方法
??? 杂交有限元方法是三角形单元和等参四边形单元相交错结合的有限元方法，这种方法对于具有复杂边界和长宽尺寸悬殊的水域是非常合适的[3,4]。
??? 设单元中的近似函数取：
(10) 其中Φi(x,y)为基函数，Z、U、V取相同的基函数；如果单元为三角形，Φi(x,y)为线性三角元基函数；如果单元是等参四边形，则Φi(x,y)为四结点等参基函数。
??? 将单元插值近似函数(10)代入方程(1)、(2)、(3)的Galerkin积分式中，可得杂交有限元方程：
(11) (12) (13) ??? 式中A、B1、C1、D1、G、B2、C2、D2、B4、C4分别为系数矩阵，其元素分别为： Aij=(ΦiΦj),，，，，，，，，，
(a)三角形单元
(b)等参四边形单元
(c)三角形-等参四边形单元
图1 AUn+1=(A-ΔtB1-ΔtG)Un-(C1·Δt-ΔtfA)Vn-ΔtD1Zn (14) AVn+1=(-ΔtfA-ΔtB2)Un-(C2·Δt-A+ΔtG)Vn-ΔtD2Zn (15) 然后通过离散式(11)，利用Un+1、Vn+1和Zn，求得Zn+1值： AZn+1=AZn-(ΔtB3+ΔtC3)Zn-(ΔtB4Un+1+ΔtC4Vn+1) (16) (17) 其中Φi(x,y)为基函数；如果单元为三角形，Φi(x,y)为线性三角元基函数；如果单元是等参四边形，则Φi(x,y)为四结点等参基函数。
??? 采用Garlerkin方法，方程(4)可写为：
(18) ??? 然后在计算节点Pi的相关单元中寻找其迎风单元，并将dcdt在迎风单元上离散，设(x1,y1),(xiβ,yiβ),(xiγ,xiγ)为节点Pi迎风单元ei的三个顶点坐标，Δi为迎风单元ei的面积。
??? 其中： ，这样，由式(18)可得： (19) ??? Pik(k=1,2,…,ni)为Pi的相关节点，其中：aiik=aiiβ·δikiβ+aiiγ·δikiγ，这里，E代表Pi节点的相关单元点数。 ??? (2)海床演变的数值模拟 ??? 利用向前差分格式离散床面随时间的变化项，然后代入方程(6)，整理得： (20)
2.1模型域的确定及剖分 河口海域模型域：以工农兵闸及以
文档评论（0）