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

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

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-剖面二维非恒定悬移质泥沙扩散方程的数值方法.doc

  • 资源星级:
  • 资源大小:74.34KB   全文页数:7页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-剖面二维非恒定悬移质泥沙扩散方程的数值方法.doc

水利工程论文剖面二维非恒定悬移质泥沙扩散方程的数值方法摘要通过讨论剖面二维非恒定泥沙扩散方程的数值方法,建立了一种用于求解含沙量分布沿程变化的差分格式ZC格式并通过一个具体的数值例子说明了计算的方法步骤。关键词扩散方程差分格式精度稳定性1引言数学模拟方法正在成为研究河流泥沙问题的重要手段。目前,一维数学模型发展较成熟,已广泛应用于模拟长河段的长期变形,但它只能给出河段平均冲淤深度的沿程变化,如需了解短河段的河床变形细节,则要采用二维以至三维数学模型。不论是一维数学模型还是平面二维维数学模型,都不能反映含沙量沿垂线的分布状况,并忽略了含沙量沿垂线分布对垂线平均含沙量变化过程的影响。要解决这类问题,必须建立剖面二维数学模型。这种模型主要通过解剖面二维泥沙扩散方程来研究悬移质泥沙沿水深的分布及含沙量的变化过程,对水电站进口和其它引水工程的引水口高程的确定都能提供较好的数值模拟。泥沙扩散方程实际上是一个变系数的二阶线性偏微分方程,这样的方程在各种复杂边界条件下求解是极为困难的。求扩散方程的解析解在数学上存在着难以克服的困难,往往只能通过对方程的简化,才能得到一些简单边界条件下的解析解,在这方面,A.A.Kalinske、野满隆治、W.E.Dobbins、俞维强、张启舜、韦直林2等都做了有益的尝试求扩散方程的数值解曾经因为缺乏高效率的计算工具而难以实现,直到60年代后,随着计算机的广泛应用,在各种复杂边界条件下求扩散方程的数值解不但成为可能,而且得到迅速的发展,在这方面,曹志先、崔侠4等做了大量工作,取得了很多成果。数值方法相对于解析方法在求解偏微分方程上有着明显的优势,即简单灵活、计算方便快捷,但要寻找一种精度高、稳定性好、计算方便的差分格式也并非易事。本文拟在前人研究的基础上着重讨论剖面二维泥沙扩散方程的数值解问题,希望能提供一种精度高、稳定性好、计算方便的数值解。2基本方程剖面二维泥沙扩散方程的形式为1式中x,y为水流方向和铅直方向的维轴u,v分别为沿水流方向和铅直方向的时均流速εsx,εsy分别为水流方向和铅直方向的泥沙扩散系数ω,S分别为泥沙静水沉速和含沙量。对于式1的求解,研究者一般会对它进行不同程度的简化,为此引入以下假定中的一种或几种A.非恒定流可以概化为梯级式恒定流,即B.在一个时段内,认为泥沙运动可以概化为处于恒定状态,即C.在二维流动中,纵向扩散系数与方程其他项相比,可以忽略不计,即认为方程右端第一项可以忽略D.认为悬移质泥沙粒径均一,即ωconstE.认为水流为二维均匀流,即v0。为简单起见,我们讨论的范围限于水流条件为二维非恒定均匀流,悬移质泥沙粒径均匀,为此引入假定C、D、E。这时,泥沙扩散方程为2目前,对εs的变化规律研究得不很充分,一般假定εsβεm3其中εm为动量传递系数,β为修正值。由勃兰特尔掺长理论可得εmκuy/1y/h4式中κ为卡门常数,u为摩阻流速。对于u,我们取卡曼勃兰特尔对数流速分布公式umaxu/u1/κlnh/y5令Wωβκu/h12y/h,2可变形为63差分方程3.1网格的剖分tτ,X方向的xh1,Y方向的空间步长为Δyh2,这样形成如下网格Dh{xj,yl,tn|xjjh1,yllh2,tnnτ|}其中j0,...,Nl0,...,Mn≥03.2构造差分格式通过对流方程和扩散方程的差分格式的构造,我们可以得到对流扩散方程的差分格式。由于隐式格式稳定性好,考虑CrankNicholson型隐式格式。为此,引入差分算子记号δ2ySnj,lSnj,l12Snj,lSnj,l1LxSnj,lSnj1,lSnj1,lLySnj,lSnj,l1Snj,l1为了看得更清楚,暂且取h1h2h.对式6离散,则CN格式为Sn1j,lSnj,l/τu/4hLxSnj,lSn1j,lεs/2h2δ2ySn1j,lSnj,lW/4hLySn1j,lSnj,l7CN格式的精度是二阶的,绝对稳定。但对于二维问题,由7导出的方程组,其系数矩阵不是三对角矩阵,不能用追赶法求解。因此,考虑构造交替方向的隐式格式命名为ZC格式Sn1/2j,lSnj,l/τ/2u/2hLxSnj,lεs/h2δ2ySn1/2j,lW/2hLySn1/2j,l8Sn1j,lSn1/2j,l/τ/2u/2hLxSn1j,lεs/h2δ2ySn1/2j,lW/2hLySn1/2j,l9可以看出,计算Sn1j,l是由两步组成的,每一步仅是一个方向的隐式,故用两次追赶法即可。3.3精度分析现在,我们考虑ZC格式的精度。先设法消去过渡值Sn1/2j,l,为此,将8和9两式相加,可得Sn1j,lSnj,lτu/2hLxSn1j,lSnj,lτεs/2h2δ2ySn1/2j,lτW/4hLxSn1/2j,l10将8和9两式相减,可得Sn1/2j,lSn1j,lSnj,l/2τεs/4h2δ2ySnj,lSn1j,lτW/8hLySnj,lSn1j,l11把式11代入10,变形整理,可得Sn1j,lSnj,l1τuεs/8h3Lxδ2yτ2uW/16h2LxLySnj,lSn1j,lτuεs/2h2δ2yτW/4hLyτu/4hLx12设Sx,y,t是12的精确解,并假定Sx,y,t关于t三次连续可微,关于x,y四次连续可微,那么利用Taylor级数展开可得Sxj,yl,tn1Sxj,yl,tn/τ1τ2εs/8h3Lxδ2yτ2uW/16h2LxLySxj,yl,tn1Sxj,yl,tn/ττεs/2h2δ2yτW/4hLyτu/4hLxOτ2h213由此可见,ZC格式具有二阶精度。3.4稳定性分析现在,我们来讨论ZC格式的稳定性。为此,把式12变形整理得1τεs/2h2δ2yτW/4hLy1τu/4hLxSn1j,l1τεs/2h2δy2τW/4hLy1τu/4hLxSnj,l14由式14可得出过渡因子为Gτ,k12τεs/h2sin2k2h/2iτW/2hsink2h1iτu/2hsink1h/12τεs/h2sin2k2h/2iτW/2hsink2h1iτu/2hsink1h15令a2τεs/h2sin2k2h/2,bτW/2hsink2h,cτu/2hsink1h,则|G|2|1a2b22bi/1a2b2|2|1c22ci/1c2|211a2b224b2/1a2b2214a24a4b24a3/1a2b2216显然,对于任意的τ,h,|Gτ,k|21,ZC格式是绝对稳定的。4数值计算4.1边界条件我们考虑初边值问题。1初始条件用Rouse公式给出含沙量沿垂线分布Sx,y,0Sa5a/hazhy/yz17式中zω/ku为悬浮指标,Sa为近底含沙量,h为水深,一般取a0.01~0.05h。2水面条件yh183底部边界条件ya19式中Sa为近底挟沙力,即输沙平衡时的近底售沙量Sa.4近底含沙量计算近底含沙量在求解泥沙扩散方程时具有边界条件性质,它选取的正确与否,意味着所给边界条件是否正确。实际工程中一般缺乏实测资料,近底含沙量不易测定。这里,我们利用水流挟沙力和含沙量沿垂线分布公式来反求近底含沙量4。已知断面平均挟沙力为Sku3/ghωm20假定SaαS21输沙平衡时,含沙量沿垂线分布用Rouse公式17表示,用17表达挟沙力的垂线分布,然后沿垂线积分得断面平均挟沙力为22将22与21比较,可得234.2计算步骤为方便计算,将式8和9式变形整理,并对X,Y方向取不同的空间步长τεs/2h22τW/4h2Sn1/2j,l1τεs/h221Sn1/2j,lτεs/2h22τW/4h2Sn1/2j,l1τu/4h1LxSnj,lSnj,l24τu/4h1Snj1,lSnj,lτu/4h1Snj1,lτεs/2h22ζ2ySn1/2j,lτW/4h2LySn1/2j,lSn1/2j,l25在一个时间层第n层内,计算分两步进行第一步,对式24用追赶法求第n1/2层的过渡值。令C1τu/4h1,C21,C3τu/4h1,E1τεs/2h22,E2τW/4h2D1τεs/2h22τW/4h2,D2τεs/h221,D3τεs/2h22τW/4h2l1时,D1Sn1/2j,0D2Sn1/2j,1D3Sn1/2j,2C3LxSnj,1Snj,1l2时,D1Sn1/2j,0D2Sn1/2j,1D3Sn1/2j,2C3LxSnj,1Snj,1lM时,D1Sn1/2j,M1D2Sn1/2j,MD3Sn1/2j,M1C3LxSnj,MSnj,M令HlSnj,lC3LxSnj,l1lMH1Snj,1C3LxSnj,1D1Sn1/2j,0HMSnj,MC3LxSnj,MD3Sn1/2j,M1其中Sn1/2j,0和Sn1/2j,M1由边界条件给出,则用矩阵形式表示为(26)第二步,再对25式用追赶法求第n1层的值令FjE1δ2ySn1/2j,lE2LySn1/2j,lSn1/2j,l1jNF1E1δ2ySn1/21,lE2LySn1/21,lSn1/21,lC1Sn10,lFNE1δ2ySn1/2N,lE2LySn1/2N,lSn1/2N,lC3Sn1N1,l其中Sn10,l和Sn1N1,l由边界条件给出,则同理可得矩阵方程27这样,按此步骤一层层地计算。4.3数值模拟合理性分析受所掌握的实测资料的限制,目前尚无法对本文提出的算法与含沙量沿垂线分布的实测值进行对比。我们用库里阿雷克沉沙池5的实测资料作了垂线平均值沿程变化的比较。该沉沙池的主要数据为池深h1.53m平均流速u0.12m/s泥沙沉速ω0.0176cm/s悬浮指标Z0.01。计算时取卡门常数κ0.4,a0.05h。表1给出了计算值和实测值,结果表明,计算值和实测值比较符合。表1断面平均含沙量验证单位kg/m3Verificationsofcrosssectionalaveragesedimentconcentrations距离m020040060080010001200计算值3.0122.5732.0711.6391.2090.8490.539

注意事项

本文(水利工程论文-剖面二维非恒定悬移质泥沙扩散方程的数值方法.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5