资源描述：
《一维明渠非恒定流的数值计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、地下水蒸发影响相当n一,)的公式计算间距,一般不致造成较大误差,如对其计算结果稍作加大修正,则效果更好.丫认件叶玉叶瑟石图Zb与1/”关系参考文献,.tl〕张友义考虑地下水不同蒸发与埋深关系指数的田间排水沟间距计算水利水电科学研究,、.,.院科学研究论文集第10集(灌溉排水)水利电力出版社1982年、,.,汇2]瞿兴业张友义考虑蒸发影响和脱盐要求的田间排水沟(管)间距计算水利学报x981年第5期.、、、,〔3〕瞿兴业张友义苏锦星后同德地下水蒸发与埋深关系指数n=3的非稳定渗流排水计算.水利学报,1983年第9期。、,。〔4〕张友义苏敏珍地下水蒸发与埋深关系指数n二4的非稳定渗流排水计算
2、水利水电科学,、、。,研究院科学研究论文集第25集(水资源灌溉与排水水利史)水利电力出版社1986年。,.,.〔5」沙金煊地下排水系统的近似计算水利学报1985年第3期、、,.,198.上6〕李庆扬王能超易大义数值分析(第四章)华中工学院出版社2年一维明渠非恒定流的数值计算汪定扬郑一凡明道华(长江水利水电科学研究院),,,1980年提出插值元法后先后用此法计算了一维行波一维热传导二维圆往绕流,三维圆球绕流等问题,在单元布置较稀条件下,计算结果的前二位数与解析解相.,词〔”1984年开始计算非线性偏微分方程问题如葛洲坝三江引航道一维非恒定,,,,t’3J三狭库流区滑波涌浪以及三峡初设二级
3、船闸中间渠道非恒定流等同题计算结果都_l三模型试验结果基本一致,而用加权差分法达不到此精度.能与原型观测值和水,“,1“,“。最近又计算了三维重力坝溢流问题(有自ilI面)二维明渠908027。及360弯,,.〔`5」取得了进一步的成果道非恒定流以及二维滑波涌浪问题本文主要介绍用摘位.元法建立的一维明渠非恒定流数学模型的计算格式,方法和部分成果,供大家参考、一插值元法(1)的基本公式,,对于一维明渠非恒定流来说插值元法(l)的基本单元形式如图l所示其中口,N断面为出,l.斯面为进口共N+个断面,`时’旬轴’}、.它、犷挤,跳八j卜协`拼.、、困1单元格式荟二尤一熟,,二t一幻·中:令任
4、一函数在i点附近的插值式..一;)气óǔ一亡一卜一S亡一,)2,承舀普)杏刀(刀一刀一功户刀(刀一刀)`;甲一s(P:;尹十,2尹+刀:,刀:花瓦二瓦又丁卞一普)普(勺一刀(一刀)12一:`.(刀一刀)(刀一刃(彗套)(普一舀)止己`·飞3;.1:」2刀刃普普一1,2,3,4(为单元上局部编号)(2,,根据式(l)2)可求得在夕时刻i点函数(苦一。,一。)的偏导数为9尹9尹、,一1,一,一C;俨盖+C二价二十C:尹{二C生中{+C几笋{+C;尹;百牙一币歹式,c、一一俞熬万不而尝加骊不’C:一一一丁祥打不炭扮面’C,一一嗜资斗备斌午,.:一刁`,1``其中舀x套=一刁x由于时间步长一般
5、取等间距,故在j时刻主点函数对时间的偏导数为一2一’9尹旦竺}(3尹:+沪{一4尹{)·a0一2』tt一沙刀1曰.~一、一维明渠非恒定流的计算根据式(3),(4)可求得一维圣维南方程离散化后的计算格式.一维明渠缓变:非恒定流圣维南方程9尸刁O一一,,凡又-不云一一O+口乙CI凡9“忑之。拜.抚9祝9(}环~+-一~一)!只~~二下厂+万二十g嘴甘一于子丁于I一兀U奋夕X一夕公乙夕毛七五卜,设渠道断,:求解时面为梯形则须补充下述三个关系式丫,+,h~(一B+户4oF)/Zm(7)、`,R=F(/B+21(l/+m)(8)·,’牲,,33,K=一/CR=/R(9):;F;抚;;n式中F过水
6、断面积Q一以流量其中平均流速h水深糙——.——;;R;率m边坡系数水力半径B底宽——,—:将渠道分为N+l个断面下面介绍进出口及中间各断面的计算公式口,(一)进(0断面)的计算前述式(4)在刁t和刁x较小条件下才有很好精,玄,刁,度因它是点的导数表示式但一般工程实际问题中沂和x都较大此时可用一般,,,:差分格式代之若进口条件是流量过程Q(t)连续方程以下述形式为好即一`刁云+(Fu刁xl=o10(F右一F丢)/{{一Q仕))/()u,名三Q(亡)/F名(1一a)或。{一Q(幻/F`一ag(川一孔扫沂/刁x.(b)na一0.51.0是考虑。l断面间水面坡降对川的影响,当渠道中流速大于Q(
7、一一t)/F石时,需按此考虑.川可根据连续方程,又可根据运动方程求,即:1。。。,连续方程瓦一(一B+了B乙十4饥F石)/2饥.,一`一`.。了一,。。当m<015时h~h十(F孟一F石)/(B+2仍h孟)一’··::,uu运动方程h~川+刁x(否一孟)/g(毋)+(u{u;一u石u孟)/2纽..~al“1,al一0川h+l(一)h6一10(12),.。,一1,,al,,计算时若取则格式是最稳定的但波幅却偏小<06时波幅增大易引起计算过程的不稳
