消能錐形閥過流狀態(tài)研究與優(yōu)化

秦武+李志鵬+喻哲欽+童成彪+巨偉

摘要：消能錐形閥是一種用于中、高水頭下的消能控流閥門，其內(nèi)部結(jié)構(gòu)設(shè)計的好壞直接影響著控流效果。為探究消能錐形閥在不同開度下，過流流量及內(nèi)部流場的變化情況。進行了模型試驗與CFD數(shù)值模擬計算，試驗中發(fā)現(xiàn)開度從100%關(guān)至85%時出現(xiàn)了流量反而上升、進口壓力下降的情況。利用CFD數(shù)值模擬的方法對內(nèi)部流場進行了可視化研究。研究結(jié)果表明，在全開狀態(tài)下錐形閥內(nèi)部存在以下問題：過強的局部高速高壓區(qū)域、低速流動死區(qū)范圍過大、渦流區(qū)過度影響高速過流部分。并對錐形閥進行了一定的結(jié)構(gòu)改進，優(yōu)化了其內(nèi)部流場情況。

關(guān)鍵詞：錐形閥；試驗；CFD；開度；過流能力；內(nèi)部結(jié)構(gòu)；優(yōu)化設(shè)計

中圖分類號：TV131.4文獻標(biāo)志碼：A文章編號：1672-1683（2017）01-0193-06

Abstract：Energy dissipating fixed cone valve is used under medium or high water head，and the design of its internal structure directly affects its flow control efficiency.To study the changes of the valve′s internal flow field and flow rate when the closed switch is opened to varied degrees，the writer conducted model tests and CFD numerical simulation computation，and found that when the switch was turned down from 100% open to 85% open，the flow rate increased and pressure around the inlet dropped.The writer also conducted visualization research on the valve′s internal flow field through CFD numerical simulation，and the results showed the following problems inside the valve when the switch was fully open：excessively high speed and high pressure in partial areas，extensive low-flow stagnation areas，excess impact of eddy area on the high-speed flow area.Finally，the writer optimized the valve′s structure and its internal flow field.

Key words：fixed cone valve；test；CFD；opening degree；flow capacity；internal structure；optimal design

泄水消能建筑是水利水電工程建設(shè)中的重要組成部分。在水利水電工程中，泄水消能建筑的建設(shè)成本達到了總建設(shè)成本的三成以上[1]， 所以優(yōu)化泄水消能設(shè)備使其在相同的條件下達到更好的使用效果具有重要的工程應(yīng)用價值。消能錐形閥一般應(yīng)用于較高水頭、大流量下的消能控流工作，其內(nèi)部流場較為復(fù)雜，通過內(nèi)部消能結(jié)構(gòu)促生高速流撞擊、漩渦、摩擦、擴散等現(xiàn)象，以此將過大的出流能量轉(zhuǎn)化消除[2-3]。魏文禮等人研究并提出了利用嵌套式結(jié)構(gòu)提高錐形閥的消能效果[4-5]。Dallin Stephens，Michael C.Johnson等人為提高錐形閥的性能效果，創(chuàng)新地在錐形閥內(nèi)部增加了金屬罩部件。 [6-7]。李穎，楊東升在對內(nèi)流式錐形閥的流場情況進行了仿真計算，并對結(jié)構(gòu)進行了優(yōu)化[8-9]。鄭淑娟等人對錐形閥內(nèi)部流場進行了預(yù)測得到了更精確的液動力計算公式[10-12]。Ji Kang Bo在應(yīng)用研究中驗證了流固耦合的多學(xué)科聯(lián)合優(yōu)化設(shè)計在錐形閥領(lǐng)域的可行性[13]。國外科研團隊則通過大量的試驗研究認為消能錐形閥可以對高水頭出流達到70%以上的消能效果[14-15]。蘭州理工大學(xué)閔為等人利用流固耦合的模擬方法得出合理選取半錐角度可以改善錐形閥的穩(wěn)定性、優(yōu)化出流條件[16-17]。Jalal M.Jalil等對錐形閥數(shù)值研究和試驗進行了聯(lián)合研究，對閥芯位置、渦流半徑及強度等對錐形閥整體性能影響進行了探討[18]。Kelly Roes則通過研究認為消能錐形閥內(nèi)部結(jié)構(gòu)尺寸與部件的選擇設(shè)計對消能錐形閥的整體消能效果起到關(guān)鍵作用[19]。

消能錐形閥除了消減上游過大水頭能量外還擔(dān)負著一定的流量調(diào)節(jié)作用。文本目的在于探究不同開度下消能錐形閥過流流量、進出口壓力及內(nèi)部流場的變化情況；并嘗試對其進行優(yōu)化改進。研究中，先進行不同開度下的錐形閥模型試驗。再利用CFD數(shù)值模擬的手段進一步對內(nèi)部流場進行分析，并根據(jù)分析結(jié)果對錐形閥進行一定的優(yōu)化改進。

1 錐形閥模型試驗

1.1 試驗?zāi)Ｐ团c裝置

本次試驗?zāi)Ｐ鸵阅炒笮烷y門企業(yè)提供的消能錐形閥為原型制作。根據(jù)試驗條件對原型進行一定比例的縮小處理。原型錐形閥進口直徑為2.5 m，縮小模型的幾何比例尺為30。本次搭建了完整的水力工程試驗臺進行試驗工作，使模型與原型滿足弗勞德數(shù)相似，其模型流速比例尺為5.48，流量比例尺為4 929.50。

圖1為消能錐形閥模型內(nèi)部結(jié)構(gòu)情況，為了更清楚地表示錐形閥開度定位情況，錐形閥內(nèi)部只繪制出錐形閥芯與啟閉套筒。該消能錐形閥通過啟閉套筒在水平方向的前后移動來改變開度的大小。圖1L中是啟閉套筒最大開啟行程，L1是當(dāng)前位置的開啟行程。當(dāng)套筒關(guān)閉至與錐形閥芯相貼合時，達到對過流的完全截止，即L1=0；同時，啟閉套筒也可開啟至最大行程開度，即L1=L。并將消能錐形閥的行程開度定義為K=L1/L×100%。

試驗?zāi)Ｐ筒捎酶邚姸鹊挠袡C玻璃材料制作，以更好的觀察錐形閥內(nèi)部的流動情況。在試驗裝置本上，次試驗采用高位供水水箱提供穩(wěn)定的恒壓進流。在錐形閥前后連接了測壓管測量其進出口壓力值的大小。并在試驗裝置出流水槽處設(shè)置了量水三角堰。模型試驗裝置見圖2。

壓力降低的反?，F(xiàn)象。故認為該錐形閥存在有全開狀態(tài)下過流能力不佳的問題，且該變化趨勢在實際應(yīng)用中將對壓力流量的調(diào)節(jié)產(chǎn)生不利影響。為了探究該現(xiàn)象發(fā)生的原因并優(yōu)化消能錐形閥的過流條件。在后續(xù)研究中采用CFD數(shù)值模擬的方法對不同開度下錐形閥內(nèi)部流場進行了相同工況條件下的數(shù)值模擬分析。

2 數(shù)值模擬分析

2.1 三維建模

本次采用Pro/E軟件進行錐形閥三維模型的建立。三維建模是為CFD流體計算提供幾何計算域，為了保證模擬與試驗的一致性，三維建模尺寸采用與試驗?zāi)Ｐ拖嗤壤哌M行縮小。本次研究的消能錐形閥，內(nèi)部部件包括：錐形閥芯、環(huán)形孔套、啟閉套筒。建模過程中先建立錐形閥整體流道的幾何模型。再分別建立啟閉套筒、閥芯、孔套等內(nèi)部結(jié)構(gòu)，通過模型切除的操作從整體流道中切除為非流通壁面區(qū)域的內(nèi)部部件，即得到所需要的三維計算域。

2.2 網(wǎng)格劃分與參數(shù)設(shè)置

網(wǎng)格劃分是CFD模擬計算中質(zhì)量要求最嚴格的一步，所建網(wǎng)格質(zhì)量的好壞將直接關(guān)系到是否能夠得到準確的收斂結(jié)果。本次研究采用ANSYS軟件組中的ICEM軟件進行網(wǎng)格劃分。由于錐形閥內(nèi)部結(jié)構(gòu)復(fù)雜，故采用幾何適應(yīng)性更好的非結(jié)構(gòu)化網(wǎng)格進行劃分。同時，為了更好的體現(xiàn)錐形閥內(nèi)部幾何結(jié)構(gòu)，計算網(wǎng)格尺寸設(shè)定為0.19。全開狀態(tài)下模型網(wǎng)格數(shù)量為1150萬，小開度下模型網(wǎng)格數(shù)量略微減少，40%開度下網(wǎng)格數(shù)量為1 000萬。經(jīng)網(wǎng)格檢查其質(zhì)量均在0.35以上，符合本次研究的計算精度要求。消能錐形閥三維內(nèi)部結(jié)構(gòu)及網(wǎng)格劃分情況見圖4。

完成網(wǎng)格劃分后，在CFX軟件中進行模擬計算及計算結(jié)果的后處理分析。CFX同屬于ANSYS軟件組，具有計算精度高、收斂效果好及計算速度快等優(yōu)點。模擬參數(shù)設(shè)置的關(guān)鍵在于確定合適的邊界條件進行數(shù)值計算。本次模擬采用收斂效果較好的進口速度與出口壓力來進行設(shè)置，內(nèi)部套筒、閥芯及孔套等結(jié)構(gòu)則設(shè)置為非流通的wall邊界。完成設(shè)置后進行迭代計算，均在500步以內(nèi)得到收斂結(jié)果。

2.3 模擬結(jié)果及誤差分析

本次模擬采用速度進口和壓力出口邊界設(shè)置，進口壓力是模擬關(guān)鍵的結(jié)果值，故可以通過對比分析數(shù)值模擬與模型試驗中進口壓力值，進行模擬計算精度的驗證。表2為模擬精度驗證結(jié)果，試驗與模型進口壓力對比見圖5。

從驗證結(jié)果可以看出，在較大開度下數(shù)值模擬具有很高的精度，70%以上開度誤差小于5%，在55%和40%的小開度下誤差上升到5%左右。從折線圖看，模擬隨開度改變壓力變化的趨勢與試驗一致。通過對CFD模擬與模型試驗方法、過程進行分析，得出可能造成誤差的原因有以下幾點：

首先，在小開度下由于流道對過流更強的阻礙作用，造成內(nèi)部流場的進一步復(fù)雜化，對數(shù)值模擬的準確計算造成了困難。其次，試驗?zāi)Ｐ驮诩庸ぶ圃焐想y免會與理想模型存在誤差。其三，試驗與模擬本身均存在難以避免的誤差。例如在試驗中實驗儀器與試驗人員觀察計量上存在不可避免的誤差；同時CFD數(shù)值模擬作為一種計算機仿真手段也難以對真實情況實現(xiàn)百分百準確的模擬。

綜上分析，由于本次研究的目的在于探究100%至85%開度時流量壓力的反常變化原因，且試驗與模擬在100%、85%、75%三個開度下均只有較小的誤差。同時，從整體來看模擬所得結(jié)果變化趨勢也與試驗一致，故認為該次模擬結(jié)果具有良好的研究可靠性。

3 分析與優(yōu)化建議

3.1 模擬結(jié)果分析

研究為探究流量壓力曲線出現(xiàn)反常變化的原因，對100%、85%、70%三個開度下內(nèi)部流場的CFD模擬結(jié)果進行對比分析。重點分析錐形閥內(nèi)部壓力和速度分布及流線變化情況，CFD模擬計算結(jié)果見圖6，其中182 cm水頭與133 cm水頭工況下的壓力、流線分布變化相似，故文中只給出182 cm水頭工況組模擬圖。

分析得出以下幾點結(jié)論。

（1）消能錐形閥各開度下內(nèi)部壓力分布趨勢基本一致。由于過流流道的突然改變使得在閥芯頭部及錐形形變處出現(xiàn)高壓區(qū)域。但是85%開度時錐形閥內(nèi)部最高壓力均小于100%及70%開度，認為錐形閥內(nèi)部結(jié)構(gòu)在85%時對過流流態(tài)的影響更小。

（2）消能錐形閥各開度下速度分布趨勢基本一致。流場在流道收縮處出現(xiàn)了局部的高速現(xiàn)象，流速的突變以及高速流的撞擊將影響整體流動狀態(tài)?？梢钥闯鲈?5%開度時錐形閥內(nèi)部最大速度均小于100%及70%開度，認為85%開度時實際的流通狀態(tài)更好。

（3）消能錐形閥大致流線變化情況相似。可以發(fā)現(xiàn)在閥芯后側(cè)和周側(cè)出現(xiàn)了幾處渦旋區(qū)域，且伴隨有流速極低的流動死區(qū)，這都將對流體的通過性能產(chǎn)生影響。并可看出100%與70%開度時，旋流與低速流動死區(qū)范圍更大，其對流道中部高速流動區(qū)域的通過性將產(chǎn)生更大阻礙作用。

3.2 錐形閥優(yōu)化建議

在錐形閥試驗與模擬研究中，均出現(xiàn)100%至85%開度下出現(xiàn)了反常的流量上升和閥前壓力下降現(xiàn)象。這在實際工程應(yīng)用中將對流量的控制調(diào)節(jié)產(chǎn)生不利影響。故研究團隊根據(jù)錐形閥設(shè)計經(jīng)驗及模擬分析結(jié)果對該類型消能錐形閥進行了一定的內(nèi)部結(jié)構(gòu)優(yōu)化，在環(huán)形孔套后部增加了錐形孔套部件以改善其內(nèi)部流場情況。

為驗證優(yōu)化后模型是否改善了在全開狀態(tài)下的內(nèi)部流場情況，對改進后模型進行了相同工況下的CFD數(shù)值模擬計算，并以之與優(yōu)化前模型模擬結(jié)果進行對比分析。模擬邊界條件采用與優(yōu)化前模型相同的進口速度與出口壓力進行設(shè)置。從模擬結(jié)果可以看出，優(yōu)化后模型有效的降低了局部最大壓力；同時減小了閥芯后側(cè)渦旋區(qū)域與流動死區(qū)的范圍，使影響較大的大面積渦旋區(qū)域轉(zhuǎn)化為面積和影響較小的多個小渦旋區(qū)域，降低了其對高速過流的阻礙作用。模擬對比分析結(jié)果見表3與圖7，模擬結(jié)果同樣以182cm水頭工況為例。

4 結(jié)論

本次為探究消能錐形閥在不同開度下過流流量及內(nèi)部流場的變化情況。進行了水力模型試驗和CFD數(shù)值模擬計算。試驗中發(fā)現(xiàn)開度從100%關(guān)至85%時出現(xiàn)過流流量上升和進口壓力下降的變化情況，這在實際工程應(yīng)用中將對過流調(diào)節(jié)產(chǎn)生不利影響。故利用CFD數(shù)值模擬的方法對不同開度下內(nèi)部流場進行了數(shù)值模擬分析，分析結(jié)果認為在全開開度下錐形閥內(nèi)部流場存在以下問題：過強的局部高壓高速區(qū)域、低速流動死區(qū)過大、渦流區(qū)過度影響高速過流部分。并對錐形閥進行了一定的結(jié)構(gòu)改進，降低了最高壓力、最大流速，減小了流動死區(qū)及渦流區(qū)對過流能力的影響，優(yōu)化了其內(nèi)部流場情況。

參考文獻：

參考文獻（References）：

[1] 高季章，董興林，劉繼廣.生態(tài)環(huán)境友好的消能技術(shù)-內(nèi)消能的研究與應(yīng)用[J].水利學(xué)報，2008（10）：1176-1182，1188。 （GAO Ji-zhang，DONG Xing-lin，LIU Ji-guang.Research and application of eco-friendly energy dissipation technology-internal energy dissipation[J].Journal of Hydraulic Engineering，2008（10）：1176-1182，1188.（in Chinese）） DOI：0559-9350（2008）10-1176-07

[2] 陳千文，易忠有.阿薩漢一級水電站消能閥選型設(shè)計[J].水力發(fā)電，2012（3）：18-20，57.（CHEN Qian-wen，YI Zhong-you.Selection of fixed cone valve for Asahan NO.1 Hydropower Station[J].Water Power 2012（3）：18-20，57（in Chinese））DOI：10.3969/j.issn.0559-9342.2012.03.006.

[3] 周琦，楊校禮，傅宗甫等.山區(qū)彎道水流特性及消能選型[J].水電能源科學(xué)，2013（1）：76-78，91.（ZHOU Qi，YANG Xiao-li，F(xiàn)U Zong-fu，et al.Study on Characteristics and Energy Dissipation of Mountainous Bend Flow[J].Water Resources and Power，2013（1）：76-78，91.（in Chinese））DOI：1000-7709（2013）01-007-04

[4] 魏文禮，李超.錐形閥附近水力特性的三維數(shù)值模擬研究[J].計算力學(xué)學(xué)報，2011（1）：108-112.（WEI Wen-li，LI Chao.3D numerical simulation of hydraulic characteristics around fixed-cone valve[J].Chinese Journal of Computational Mechanics，2011（1）：108-112.（in Chinese）） DOI：1007-4708（2011）01-0108-05

[5] 魏文禮，呂彬，劉玉玲.消能錐形閥結(jié)構(gòu)改進的三維數(shù)值模擬研究[J].應(yīng)用力學(xué)學(xué)報，2013（1）：109-114，151.（WEI Wen-li，LV Bin，LIU Yu-ling.Research on the improvement of structure for energy dissipation fixed-cone valve by 3D numerical simulation[J].Chinese Journal of Applied Mechanics，2013（1）：109-114，151.（in Chinese））DOI：1000- 4939（2013） 01-0109-06

[6] Dallin Stephens ，Michael C Johnson ，Zachary B Sharp.Design Considerations for Fixed-Cone Valve with Baffled Hood[J].Journal of hydraulic engineering.2012（2）：204-209.DOI：10.1061/（ASCE）HY.1943-7900.0000496

[7] Michael C.Johnson，M.ASCE，Rajesh Dham，M.ASCE.Innovative Energy-Dissipating Hood[J].Journal of Hydraulic Engineering ASCE，2006（8）：759-764.DOI：10.1061/（ASCE）0733-9429（2006）132：8（759）

[8] 李穎，楊東生.錐閥流場的數(shù)值模擬與結(jié)構(gòu)優(yōu)化[J].科技技術(shù)與工程，2011（21）：5180-5183.（LI Ying，YANG Dong-sheng.Numerical simulation and structure optimizing for the flow of poppet valve[J].Science Technology and Engineering，2011（21）：5180-5183.（in Chinese））DOI：10.3969/j.issn.1671-1815.2011.21.050.

[9] 楊東升.內(nèi)流式錐閥流場的數(shù)值模擬與結(jié)構(gòu)優(yōu)化[J].管道技術(shù)與設(shè)備，2012（3）：35-37.（YANG Dong-sheng.Numerical simulation and structure optimizing for the inward flow filed of poppet valve[J].Pipeline Technique and Equipment，2012（3）：35-37.（in Chinese）） DOI：10.3969/j.issn.1004-9614.2012.03.013.

[10] 鄭淑娟，權(quán)龍，陳青.閥芯運動過程液壓錐形閥流場的CFD計算與分析[J].農(nóng)業(yè)機械學(xué)報，2007（1）：168-172.（ZHENG Shu-juan，QUAN Long，CHEN Qing.Analysis and CFD simulation of the flow filed in a moving poppet valve[J].Transactions of the Chinese Society for Agricultural Machinery，2007（1）：168-172.（in Chinese））

[11] 鄭淑娟，權(quán)龍.錐臺型液壓錐形閥過流面積的CFD可視化計算[J].農(nóng)業(yè)機械學(xué)報，2014（4）：327-333.（ZHENG Shu-juan，QUAN Long.CFD visual analysis on cross section of truncated poppet valve[J].Transactions of the Chinese Society for Agricultural Machinery，2014（4）：327-333.（in Chinese））DOI：10.6041 /j.issn.1000-1298.2014.04.052

[12] 鄭淑娟，權(quán)龍，王文全.平底錐形閥內(nèi)流工況穩(wěn)態(tài)液動力的研究[J].液壓與氣動，2014（8）：55-59，108.（ZHENG Shu-juan，QUAN Long，WANG Wen-quan.Steady flow force of truncated poppet valve with convergent flow[J].Chinese Hydraulics & Pneumaticsy，2014（8）：55-59，108.（in Chinese））

[13] Ji Kang Bo.Multidisciplinary Design Optimization of a Hydraulic Poppet Valve Considering Fluid-Solid Coupling[J].Advanced Materials Research，2011，Vol.1377 （308），pp.559-562.

[14] Power-Bulk Solids Group.Cone valve discharger[J].Powder Bulk Solids.2003，21（1） 167-172.

[15] Bulk Solids Today Group.Cone Valve Silo Dischargers[J].Bulk Solids Today.2004，10（45）：112-117.

[16] 閔為，王崢嶸.不同閥座半錐角條件下的錐閥閥口流場仿真[J].蘭州理工大學(xué)學(xué)報，2012（6）：49-52.（MIN Wei，WANG Zheng-rong.Flow field simulation of valve port with different half conical angles of valve seat[J].Journal of Lanzhou University of Technology，2012（6）：49-52.（in Chinese））DOI：10.3969/j.issn.1673-5196.2012.06.011.

[17] 閔為，翼宏，王崢嶸，等.單級壓力調(diào)節(jié)閥的阻尼孔射流響應(yīng)特性研究[J].西安交通大學(xué)學(xué)報，2014（6）：80-85.（MIN Wei，YI Hong，WANG Zheng-rong.Orifice jet response characters in pressure regulating valve[J].Journal of Lanzhou University of Technology，2014（6）：80-85.（in Chinese））DOI：10.6752/xjtuxb201406014

[18] Jalal M Jalil，Sabah T Ahmed，YI Qinxue，al et.Experimental and numerical investigation of fluid flow of truncated conical poppet valve [J].International Journal of Fluid Power，2015，16 （1）：25-34

[19] Kelly Rose.Cone improvements：Advances in cone valve technology gives new levels of IBC containment[J].Solids & Bulk Handing.2009，35（7）：156-160