基于Partical模型的固液兩相流貫流式水輪機(jī)流場數(shù)值模擬
發(fā)布時(shí)間:2021-03-24 22:53
為了研究燈泡貫流式水輪機(jī)在固液兩相流動(dòng)下過流部件的磨損程度,以某水電站實(shí)際運(yùn)行的水輪機(jī)為研究對象,基于Particle模型和非均相模型,對水輪機(jī)內(nèi)部流場進(jìn)行數(shù)值模擬。分析了體積分?jǐn)?shù)為1%,泥沙粒徑為0.01、0.05、0.10和0.25 mm時(shí),過流部件工作面和背面的固體顆粒相的滑移速度及固相體積分?jǐn)?shù)。結(jié)果表明:隨著顆粒粒徑的增大,槳葉上固體顆粒相的滑移速度減小、導(dǎo)葉上滑移速度增大;在漿葉進(jìn)水邊和輪緣,導(dǎo)葉的進(jìn)水邊和輪轂處固體顆粒相的體積分?jǐn)?shù)較大,磨損較為嚴(yán)重,與現(xiàn)場真機(jī)照片磨損位置相同,說明采用Particle模型和非均相模型能準(zhǔn)確地模擬固液兩相流水輪機(jī)內(nèi)的流動(dòng)規(guī)律。
【文章來源】:熱能動(dòng)力工程. 2020,35(02)北大核心CSCD
【文章頁數(shù)】:9 頁
【部分圖文】:
貫流式水輪機(jī)流體域計(jì)算模型
為提高計(jì)算精確度,進(jìn)口段、燈泡體和尾水管流體域?yàn)榻Y(jié)構(gòu)化網(wǎng)格,活動(dòng)導(dǎo)葉與轉(zhuǎn)輪流體域?yàn)榉墙Y(jié)構(gòu)化網(wǎng)格,經(jīng)網(wǎng)格無關(guān)性驗(yàn)證,如圖2所示,最終確定網(wǎng)格單元數(shù)為8 877 746,節(jié)點(diǎn)數(shù)為1 937 620,各部分網(wǎng)格數(shù)量如表1所示,其流體域網(wǎng)格如圖3所示。表1 水輪機(jī)各部件網(wǎng)格分布情況Tab.1 Grid layout of turbine components 名稱 網(wǎng)格數(shù)量 網(wǎng)格節(jié)點(diǎn) 網(wǎng)格類型 進(jìn)口段 65 588 69 260 六面體網(wǎng)格 燈泡體 30 240 34 188 六面體網(wǎng)格 活動(dòng)導(dǎo)葉 5 528 362 1 071 163 四面體網(wǎng)格 轉(zhuǎn)輪 3 091 201 594 621 四面體網(wǎng)格 尾水管 162 355 168 388 六面體網(wǎng)格
表1 水輪機(jī)各部件網(wǎng)格分布情況Tab.1 Grid layout of turbine components 名稱 網(wǎng)格數(shù)量 網(wǎng)格節(jié)點(diǎn) 網(wǎng)格類型 進(jìn)口段 65 588 69 260 六面體網(wǎng)格 燈泡體 30 240 34 188 六面體網(wǎng)格 活動(dòng)導(dǎo)葉 5 528 362 1 071 163 四面體網(wǎng)格 轉(zhuǎn)輪 3 091 201 594 621 四面體網(wǎng)格 尾水管 162 355 168 388 六面體網(wǎng)格在非定常計(jì)算中,泥沙密度設(shè)為2 650 kg/m3,為均勻圓球顆粒,多相流模型采用非均相多相流模型,液相使用RNG k-ε湍流模型,固相使用離散相零方程模型[14-15]。相間傳輸采用Particle模型,動(dòng)量交換考慮拽力,并使用Schiller Naumann函數(shù)模型,不考慮升力、虛擬質(zhì)量力及壁面潤滑力的影響。非定常計(jì)算以定常計(jì)算為初始值,湍流耗散力采用法夫爾平均拖拽力(Favre Averaged Drag Force),系數(shù)取0.9,湍流輸送模型使用佐藤增強(qiáng)粘度(Sato Enhanced Viscosity)。進(jìn)口選擇質(zhì)量流量條件,出口為壓力出口條件。連續(xù)相設(shè)為無滑移壁面邊界,固相設(shè)置為自由滑移壁面邊界。壁面接觸模型選取體積分?jǐn)?shù)模型,動(dòng)靜轉(zhuǎn)子交界面采用凍結(jié)轉(zhuǎn)子(Frozen Rotor)。設(shè)定物理時(shí)間步長為1/ω,即轉(zhuǎn)輪旋轉(zhuǎn)角速度的倒數(shù)。先進(jìn)行了定常計(jì)算,然后將此結(jié)果作為初始值,再進(jìn)行固液兩相流非定常計(jì)算,固液兩相計(jì)算的收斂精度為10-4。
【參考文獻(xiàn)】:
期刊論文
[1]泥沙直徑和濃度對混流式水輪機(jī)轉(zhuǎn)輪內(nèi)流場的影響[J]. 廖姣,賴喜德,張興. 水力發(fā)電學(xué)報(bào). 2017(05)
[2]抽黃工程用離心泵固液兩相流數(shù)值模擬[J]. 周世杰,廖偉麗,趙亞萍,阮輝,羅興锜. 水力發(fā)電學(xué)報(bào). 2017(05)
[3]基于固液兩相流的離心泵內(nèi)部流場數(shù)值分析[J]. 廖姣,賴喜德,廖功磊,張文明. 熱能動(dòng)力工程. 2017(05)
[4]高水頭混流式水輪機(jī)內(nèi)部固液兩相流數(shù)值分析[J]. 廖姣,張興,張文明. 人民長江. 2017(07)
[5]基于改進(jìn)歐拉算法的雙吸離心泵泥沙磨損特性研究[J]. 張自超,王福軍,陳鑫,廖翠林,徐洪泉,陸力. 農(nóng)業(yè)機(jī)械學(xué)報(bào). 2017(03)
[6]低濃度固液兩相流相間阻力修正模型研究[J]. 張自超,王福軍,陳鑫,肖若富,姚志峰. 農(nóng)業(yè)機(jī)械學(xué)報(bào). 2016(12)
[7]水輪機(jī)泥沙磨損兩相湍流場數(shù)值模擬[J]. 黃劍峰,張立翔,姚激,龍立焱. 排灌機(jī)械工程學(xué)報(bào). 2016(02)
[8]泥沙濃度及粒徑對水輪機(jī)轉(zhuǎn)輪內(nèi)部流動(dòng)影響的數(shù)值分析[J]. 張廣,魏顯著. 農(nóng)業(yè)工程學(xué)報(bào). 2014(23)
[9]后掠式葉片軸流泵固液兩相流數(shù)值模擬與優(yōu)化[J]. 施衛(wèi)東,邢津,張德勝,陳刻強(qiáng),程成. 農(nóng)業(yè)工程學(xué)報(bào). 2014(11)
[10]固液兩相流離心泵內(nèi)部流場數(shù)值模擬與磨損特性[J]. 汪家瓊,蔣萬明,孔繁余,宿向輝,陳浩. 農(nóng)業(yè)機(jī)械學(xué)報(bào). 2013(11)
博士論文
[1]燈泡貫流式水輪機(jī)協(xié)聯(lián)關(guān)系及性能研究[D]. 李正貴.蘭州理工大學(xué) 2014
本文編號:3098544
【文章來源】:熱能動(dòng)力工程. 2020,35(02)北大核心CSCD
【文章頁數(shù)】:9 頁
【部分圖文】:
貫流式水輪機(jī)流體域計(jì)算模型
為提高計(jì)算精確度,進(jìn)口段、燈泡體和尾水管流體域?yàn)榻Y(jié)構(gòu)化網(wǎng)格,活動(dòng)導(dǎo)葉與轉(zhuǎn)輪流體域?yàn)榉墙Y(jié)構(gòu)化網(wǎng)格,經(jīng)網(wǎng)格無關(guān)性驗(yàn)證,如圖2所示,最終確定網(wǎng)格單元數(shù)為8 877 746,節(jié)點(diǎn)數(shù)為1 937 620,各部分網(wǎng)格數(shù)量如表1所示,其流體域網(wǎng)格如圖3所示。表1 水輪機(jī)各部件網(wǎng)格分布情況Tab.1 Grid layout of turbine components 名稱 網(wǎng)格數(shù)量 網(wǎng)格節(jié)點(diǎn) 網(wǎng)格類型 進(jìn)口段 65 588 69 260 六面體網(wǎng)格 燈泡體 30 240 34 188 六面體網(wǎng)格 活動(dòng)導(dǎo)葉 5 528 362 1 071 163 四面體網(wǎng)格 轉(zhuǎn)輪 3 091 201 594 621 四面體網(wǎng)格 尾水管 162 355 168 388 六面體網(wǎng)格
表1 水輪機(jī)各部件網(wǎng)格分布情況Tab.1 Grid layout of turbine components 名稱 網(wǎng)格數(shù)量 網(wǎng)格節(jié)點(diǎn) 網(wǎng)格類型 進(jìn)口段 65 588 69 260 六面體網(wǎng)格 燈泡體 30 240 34 188 六面體網(wǎng)格 活動(dòng)導(dǎo)葉 5 528 362 1 071 163 四面體網(wǎng)格 轉(zhuǎn)輪 3 091 201 594 621 四面體網(wǎng)格 尾水管 162 355 168 388 六面體網(wǎng)格在非定常計(jì)算中,泥沙密度設(shè)為2 650 kg/m3,為均勻圓球顆粒,多相流模型采用非均相多相流模型,液相使用RNG k-ε湍流模型,固相使用離散相零方程模型[14-15]。相間傳輸采用Particle模型,動(dòng)量交換考慮拽力,并使用Schiller Naumann函數(shù)模型,不考慮升力、虛擬質(zhì)量力及壁面潤滑力的影響。非定常計(jì)算以定常計(jì)算為初始值,湍流耗散力采用法夫爾平均拖拽力(Favre Averaged Drag Force),系數(shù)取0.9,湍流輸送模型使用佐藤增強(qiáng)粘度(Sato Enhanced Viscosity)。進(jìn)口選擇質(zhì)量流量條件,出口為壓力出口條件。連續(xù)相設(shè)為無滑移壁面邊界,固相設(shè)置為自由滑移壁面邊界。壁面接觸模型選取體積分?jǐn)?shù)模型,動(dòng)靜轉(zhuǎn)子交界面采用凍結(jié)轉(zhuǎn)子(Frozen Rotor)。設(shè)定物理時(shí)間步長為1/ω,即轉(zhuǎn)輪旋轉(zhuǎn)角速度的倒數(shù)。先進(jìn)行了定常計(jì)算,然后將此結(jié)果作為初始值,再進(jìn)行固液兩相流非定常計(jì)算,固液兩相計(jì)算的收斂精度為10-4。
【參考文獻(xiàn)】:
期刊論文
[1]泥沙直徑和濃度對混流式水輪機(jī)轉(zhuǎn)輪內(nèi)流場的影響[J]. 廖姣,賴喜德,張興. 水力發(fā)電學(xué)報(bào). 2017(05)
[2]抽黃工程用離心泵固液兩相流數(shù)值模擬[J]. 周世杰,廖偉麗,趙亞萍,阮輝,羅興锜. 水力發(fā)電學(xué)報(bào). 2017(05)
[3]基于固液兩相流的離心泵內(nèi)部流場數(shù)值分析[J]. 廖姣,賴喜德,廖功磊,張文明. 熱能動(dòng)力工程. 2017(05)
[4]高水頭混流式水輪機(jī)內(nèi)部固液兩相流數(shù)值分析[J]. 廖姣,張興,張文明. 人民長江. 2017(07)
[5]基于改進(jìn)歐拉算法的雙吸離心泵泥沙磨損特性研究[J]. 張自超,王福軍,陳鑫,廖翠林,徐洪泉,陸力. 農(nóng)業(yè)機(jī)械學(xué)報(bào). 2017(03)
[6]低濃度固液兩相流相間阻力修正模型研究[J]. 張自超,王福軍,陳鑫,肖若富,姚志峰. 農(nóng)業(yè)機(jī)械學(xué)報(bào). 2016(12)
[7]水輪機(jī)泥沙磨損兩相湍流場數(shù)值模擬[J]. 黃劍峰,張立翔,姚激,龍立焱. 排灌機(jī)械工程學(xué)報(bào). 2016(02)
[8]泥沙濃度及粒徑對水輪機(jī)轉(zhuǎn)輪內(nèi)部流動(dòng)影響的數(shù)值分析[J]. 張廣,魏顯著. 農(nóng)業(yè)工程學(xué)報(bào). 2014(23)
[9]后掠式葉片軸流泵固液兩相流數(shù)值模擬與優(yōu)化[J]. 施衛(wèi)東,邢津,張德勝,陳刻強(qiáng),程成. 農(nóng)業(yè)工程學(xué)報(bào). 2014(11)
[10]固液兩相流離心泵內(nèi)部流場數(shù)值模擬與磨損特性[J]. 汪家瓊,蔣萬明,孔繁余,宿向輝,陳浩. 農(nóng)業(yè)機(jī)械學(xué)報(bào). 2013(11)
博士論文
[1]燈泡貫流式水輪機(jī)協(xié)聯(lián)關(guān)系及性能研究[D]. 李正貴.蘭州理工大學(xué) 2014
本文編號:3098544
本文鏈接:http://sikaile.net/kejilunwen/shuiwenshuili/3098544.html
最近更新
教材專著