水電站近壩區(qū)水沙輸移二維數(shù)值模擬
發(fā)布時(shí)間:2021-11-09 14:37
以薩拉康水利工程為例,采用平面二維水沙模型對(duì)近壩區(qū)域的水沙運(yùn)動(dòng)進(jìn)行計(jì)算研究,分析了敞泄排沙條件下上引航道、泄洪閘、廠房等區(qū)域的泥沙淤積分布。結(jié)果表明,薩拉康水利樞紐工程運(yùn)行30年后,泄洪閘前的泥沙淤積不會(huì)影響閘門安全運(yùn)行;上引航道口門區(qū)主要受主流回流的影響,產(chǎn)生回流淤積;廠房段的攔沙坎與沖沙底孔可有效保護(hù)廠房進(jìn)水口。
【文章來源】:水電能源科學(xué). 2020,38(05)北大核心
【文章頁(yè)數(shù)】:5 頁(yè)
【部分圖文】:
圖6X2的沖淤對(duì)比Fig.6Comparisonoferosionandsedimentation
;④干濕邊界的處理。(2)泥沙模型求解。采用TVD格式的有限體積法來求解泥沙輸運(yùn)方程,以保證守恒與提高計(jì)算精度。由于不采用ELM方法求解對(duì)流項(xiàng),因此計(jì)算的時(shí)間步長(zhǎng)受到了Courant數(shù)穩(wěn)定條件的限制,將水流模型的時(shí)間步長(zhǎng)Δt=tn+1-tn劃分為若干個(gè)Δt′,每個(gè)分步Δt′大小由整個(gè)流場(chǎng)對(duì)流作用的強(qiáng)弱決定,則方程第m分步到第m+1分步采用有限體積法離散。有限體積法中網(wǎng)格i與相鄰網(wǎng)格j的關(guān)系見圖1。圖1有限體積法中網(wǎng)格i與相鄰網(wǎng)格j的關(guān)系Fig.1RelationshipbetweengridiandadjacentgridjinFVM此外,當(dāng)水流通過彎道時(shí),受離心慣性力與重力的共同作用,產(chǎn)生徑向和豎向流速,與主流疊加在一起形成螺旋流。螺旋流引起的橫向輸沙影響著河相的形成,需采用橫向流速分布公式對(duì)泥沙的輸移方向進(jìn)行修正。實(shí)際中常采用羅索夫斯基的基于縱向流速沿垂線分布的對(duì)數(shù)規(guī)律推出的橫向流速垂線分布公式[7],本文采用簡(jiǎn)化后的羅索夫斯基的橫向流速垂線分布公式:vr=ακ2U2槡+V2hr0.4-0.1n槡gκh1/6()(5)式中,vr?yàn)閺澋罊M向流速;U2槡+V2為垂向平均流速,U、V分別為沿x、y方向的流速;κ為卡門常數(shù);r?yàn)榇咕所在點(diǎn)的曲率半徑;h為水深;α為橫向流速垂向修正系數(shù),0.4≤α≤1,其值的選取與橫向流速對(duì)應(yīng)的相對(duì)水深有關(guān),且值越大,越靠近河床底部。"#"計(jì)算?
?船閘等建筑物組成。水庫(kù)正常蓄水位220m,庫(kù)容8.63×108m3,水庫(kù)長(zhǎng)度99.2km,電站總裝機(jī)容量660MW。"#"#"河段地形和網(wǎng)格劃分由于水電樞紐工程近壩區(qū)域水流結(jié)構(gòu)復(fù)雜,故采用非結(jié)構(gòu)化三角形網(wǎng)格劃分計(jì)算區(qū)域,并對(duì)引航道、泄洪閘、廠房段局部加密,河段單元網(wǎng)格共計(jì)31159個(gè),單元節(jié)點(diǎn)共16121個(gè)。利用庫(kù)區(qū)的1∶5000地形資料對(duì)網(wǎng)格進(jìn)行插值,生成的二維水沙數(shù)學(xué)模型模擬的地形及樞紐布置示意圖見圖2。圖2壩區(qū)二維水沙數(shù)學(xué)模型模擬河段地形Fig.2Theriverterrainsimulatedbytwo-dimensionalflowandsedimentmodelindamarea"#"#$邊界條件二維模型的進(jìn)口水沙邊界以距壩前3.5km的斷面作為入口邊界條件,以正常蓄水位及排沙運(yùn)行方式作為出口水位控制條件,壩址設(shè)計(jì)洪水及下游水位見表1。表1薩拉康壩址設(shè)計(jì)洪水、下游水位Tab.1DesignfloodanddownstreamwaterlevelatSarakangdamsite頻率設(shè)計(jì)值/%流量/(m3·s-1)水位/m225000215.85522900215.002019400213.3233.317800212.555016300211.72"#$模型驗(yàn)證"#$#!水流模型驗(yàn)證為驗(yàn)證平面二維淺水方程的適用性,選取上游距壩42m處斷面的水位與壩前斷面,比較不同洪水頻率下壩前水位壅高情況,見表2。表2不同洪水頻率下壩前水位壅高對(duì)比Tab.2C?
【參考文獻(xiàn)】:
期刊論文
[1]水庫(kù)壩區(qū)沖刷漏斗的形成機(jī)理[J]. 胡德超,池龍哲,楊瓊,王敏. 浙江大學(xué)學(xué)報(bào)(工學(xué)版). 2015(02)
[2]考慮彎道環(huán)流影響的平面二維水沙數(shù)學(xué)模型(Ⅱ)——模型的驗(yàn)證[J]. 易雨君,王兆印,張尚弘. 水力發(fā)電學(xué)報(bào). 2010(01)
[3]考慮彎道環(huán)流影響的平面二維水沙數(shù)學(xué)模型(Ⅰ)——模型的建立[J]. 易雨君,王兆印,張尚弘. 水力發(fā)電學(xué)報(bào). 2010(01)
[4]計(jì)算彎道環(huán)流的通用公式[J]. 宋志堯. 泥沙研究. 2003(04)
[5]通航建筑物引航道泥沙問題及其處理途徑[J]. 談廣鳴. 水電能源科學(xué). 1998(04)
[6]平面二維水沙數(shù)學(xué)模型[J]. 吳偉明,胡春燕,楊國(guó)錄. 水利學(xué)報(bào). 1995(10)
碩士論文
[1]基于MIKE(21、3)的河南天池抽水蓄能電站下水庫(kù)庫(kù)區(qū)水沙數(shù)值模擬研究[D]. 馬騰飛.華北水利水電大學(xué) 2016
本文編號(hào):3485550
【文章來源】:水電能源科學(xué). 2020,38(05)北大核心
【文章頁(yè)數(shù)】:5 頁(yè)
【部分圖文】:
圖6X2的沖淤對(duì)比Fig.6Comparisonoferosionandsedimentation
;④干濕邊界的處理。(2)泥沙模型求解。采用TVD格式的有限體積法來求解泥沙輸運(yùn)方程,以保證守恒與提高計(jì)算精度。由于不采用ELM方法求解對(duì)流項(xiàng),因此計(jì)算的時(shí)間步長(zhǎng)受到了Courant數(shù)穩(wěn)定條件的限制,將水流模型的時(shí)間步長(zhǎng)Δt=tn+1-tn劃分為若干個(gè)Δt′,每個(gè)分步Δt′大小由整個(gè)流場(chǎng)對(duì)流作用的強(qiáng)弱決定,則方程第m分步到第m+1分步采用有限體積法離散。有限體積法中網(wǎng)格i與相鄰網(wǎng)格j的關(guān)系見圖1。圖1有限體積法中網(wǎng)格i與相鄰網(wǎng)格j的關(guān)系Fig.1RelationshipbetweengridiandadjacentgridjinFVM此外,當(dāng)水流通過彎道時(shí),受離心慣性力與重力的共同作用,產(chǎn)生徑向和豎向流速,與主流疊加在一起形成螺旋流。螺旋流引起的橫向輸沙影響著河相的形成,需采用橫向流速分布公式對(duì)泥沙的輸移方向進(jìn)行修正。實(shí)際中常采用羅索夫斯基的基于縱向流速沿垂線分布的對(duì)數(shù)規(guī)律推出的橫向流速垂線分布公式[7],本文采用簡(jiǎn)化后的羅索夫斯基的橫向流速垂線分布公式:vr=ακ2U2槡+V2hr0.4-0.1n槡gκh1/6()(5)式中,vr?yàn)閺澋罊M向流速;U2槡+V2為垂向平均流速,U、V分別為沿x、y方向的流速;κ為卡門常數(shù);r?yàn)榇咕所在點(diǎn)的曲率半徑;h為水深;α為橫向流速垂向修正系數(shù),0.4≤α≤1,其值的選取與橫向流速對(duì)應(yīng)的相對(duì)水深有關(guān),且值越大,越靠近河床底部。"#"計(jì)算?
?船閘等建筑物組成。水庫(kù)正常蓄水位220m,庫(kù)容8.63×108m3,水庫(kù)長(zhǎng)度99.2km,電站總裝機(jī)容量660MW。"#"#"河段地形和網(wǎng)格劃分由于水電樞紐工程近壩區(qū)域水流結(jié)構(gòu)復(fù)雜,故采用非結(jié)構(gòu)化三角形網(wǎng)格劃分計(jì)算區(qū)域,并對(duì)引航道、泄洪閘、廠房段局部加密,河段單元網(wǎng)格共計(jì)31159個(gè),單元節(jié)點(diǎn)共16121個(gè)。利用庫(kù)區(qū)的1∶5000地形資料對(duì)網(wǎng)格進(jìn)行插值,生成的二維水沙數(shù)學(xué)模型模擬的地形及樞紐布置示意圖見圖2。圖2壩區(qū)二維水沙數(shù)學(xué)模型模擬河段地形Fig.2Theriverterrainsimulatedbytwo-dimensionalflowandsedimentmodelindamarea"#"#$邊界條件二維模型的進(jìn)口水沙邊界以距壩前3.5km的斷面作為入口邊界條件,以正常蓄水位及排沙運(yùn)行方式作為出口水位控制條件,壩址設(shè)計(jì)洪水及下游水位見表1。表1薩拉康壩址設(shè)計(jì)洪水、下游水位Tab.1DesignfloodanddownstreamwaterlevelatSarakangdamsite頻率設(shè)計(jì)值/%流量/(m3·s-1)水位/m225000215.85522900215.002019400213.3233.317800212.555016300211.72"#$模型驗(yàn)證"#$#!水流模型驗(yàn)證為驗(yàn)證平面二維淺水方程的適用性,選取上游距壩42m處斷面的水位與壩前斷面,比較不同洪水頻率下壩前水位壅高情況,見表2。表2不同洪水頻率下壩前水位壅高對(duì)比Tab.2C?
【參考文獻(xiàn)】:
期刊論文
[1]水庫(kù)壩區(qū)沖刷漏斗的形成機(jī)理[J]. 胡德超,池龍哲,楊瓊,王敏. 浙江大學(xué)學(xué)報(bào)(工學(xué)版). 2015(02)
[2]考慮彎道環(huán)流影響的平面二維水沙數(shù)學(xué)模型(Ⅱ)——模型的驗(yàn)證[J]. 易雨君,王兆印,張尚弘. 水力發(fā)電學(xué)報(bào). 2010(01)
[3]考慮彎道環(huán)流影響的平面二維水沙數(shù)學(xué)模型(Ⅰ)——模型的建立[J]. 易雨君,王兆印,張尚弘. 水力發(fā)電學(xué)報(bào). 2010(01)
[4]計(jì)算彎道環(huán)流的通用公式[J]. 宋志堯. 泥沙研究. 2003(04)
[5]通航建筑物引航道泥沙問題及其處理途徑[J]. 談廣鳴. 水電能源科學(xué). 1998(04)
[6]平面二維水沙數(shù)學(xué)模型[J]. 吳偉明,胡春燕,楊國(guó)錄. 水利學(xué)報(bào). 1995(10)
碩士論文
[1]基于MIKE(21、3)的河南天池抽水蓄能電站下水庫(kù)庫(kù)區(qū)水沙數(shù)值模擬研究[D]. 馬騰飛.華北水利水電大學(xué) 2016
本文編號(hào):3485550
本文鏈接:http://sikaile.net/kejilunwen/shuiwenshuili/3485550.html
最近更新
教材專著