馮建剛,錢向棟,張 睿

(河海大學(xué)水利水電學(xué)院,江蘇 南京 210098)

城市輸水泵站受城市土地限制和周圍建筑物的影響,往往很難按水力條件良好的要求布置,前池易出現(xiàn)漩渦、回流等不良流態(tài),水泵進(jìn)水條件惡化,影響泵站的安全可靠運(yùn)行[1],因此對(duì)此類泵站的前池流態(tài)進(jìn)行研究十分必要。近年來隨著計(jì)算流體力學(xué)和計(jì)算機(jī)技術(shù)的快速發(fā)展,數(shù)值模擬被越來越多的學(xué)者所認(rèn)可,成為實(shí)際工程中研究水體流動(dòng)特性的重要手段之一[2-6]。田艷等[7]采用RNGk-ε湍流模型對(duì)城市排水泵站進(jìn)行全流場(chǎng)模擬,揭示了城市排水泵站前池內(nèi)存在的復(fù)雜水流流態(tài)。對(duì)于這類布置較為局促的城市輸水泵站,為改善前池流態(tài),學(xué)者們對(duì)整流措施進(jìn)行了探究。魯俊[8]對(duì)泵站前池水力特性進(jìn)行數(shù)值模擬,發(fā)現(xiàn)其結(jié)果與試驗(yàn)值相符,并提出了整流工程措施。馮建剛等[9]在上海市長(zhǎng)江引水三期工程水源泵站的整體水力模型試驗(yàn)中,對(duì)底坎和八字形導(dǎo)流墩措施進(jìn)行了比較分析。黃繼宏等[10]通過分析不同導(dǎo)流墩參數(shù)對(duì)泵站前池水流流態(tài)的影響,明確了優(yōu)化導(dǎo)流墩布置的相關(guān)參數(shù)。羅燦等[11]對(duì)泵站正向進(jìn)水前池底坎的整流機(jī)理進(jìn)行了數(shù)值模擬。然而，已有研究主要集中于單類整流措施尺寸參數(shù)對(duì)前池流態(tài)的影響,對(duì)組合式整流措施應(yīng)用于城市輸水泵站的研究尚不多見。本文結(jié)合某底部涵管進(jìn)水形式的城市輸水泵站,擬組合導(dǎo)流墩和底坎,探究出適用于此類泵站的整流措施。通過對(duì)泵站前池流態(tài)進(jìn)行數(shù)值模擬研究,分析前池存在的不良流態(tài)及其成因,在對(duì)比分析的基礎(chǔ)上提出改善前池流態(tài)的工程措施,并對(duì)數(shù)值模擬結(jié)果進(jìn)行試驗(yàn)驗(yàn)證。

1 泵站模型的建立

1.1 泵站進(jìn)水建筑物概況

泵站前池及進(jìn)水池結(jié)構(gòu)形式如圖1所示。該泵站由進(jìn)水涵管引水進(jìn)入前池。前池呈平面擴(kuò)散形,擴(kuò)散角為20°,池底高程為-10.00 m。泵站設(shè)有3臺(tái)立式斜流泵機(jī)組,編號(hào)分別為1號(hào)、2號(hào)、3號(hào),單泵設(shè)計(jì)流量為6.21 m3/s。每臺(tái)水泵均單獨(dú)設(shè)進(jìn)水池,以隔墩分隔,池底與前池等高。

圖1 前池及進(jìn)水池結(jié)構(gòu)形式(單位:m)

1.2 數(shù)值模型建立

泵站前池流態(tài)處于復(fù)雜的紊流狀態(tài),雷諾數(shù)通常大于104,為不可壓流動(dòng)。采用RNGk-ε模型和SIMPLE算法求解流速場(chǎng)分布。應(yīng)用雷諾時(shí)均方程,得到定常流動(dòng)在笛卡兒坐標(biāo)系下三維不可壓縮流體的連續(xù)性方程和動(dòng)量方程[12]:

(1)

(2)

式中：ui、uj(i、j=1,2,3)分別為x、y、z方向的速度分量;xi、xj(i、j=1,2,3)分別為空間坐標(biāo)x、y、z;gi為沿i方向的質(zhì)量力;ρ為流體密度;p為流體微元體上的壓力;ν為水的運(yùn)動(dòng)黏性系數(shù);νt為紊動(dòng)黏性系數(shù)。

在笛卡兒坐標(biāo)系下,采用Pro/E建立實(shí)體模型,利用ANSYS ICEM進(jìn)行非結(jié)構(gòu)化網(wǎng)格的劃分,網(wǎng)格數(shù)400萬左右。

進(jìn)口邊界設(shè)置為總壓,出口邊界采用質(zhì)量流條件。由于前池內(nèi)水面波動(dòng)不大,故自由水面采用剛蓋假定。壁面采用無滑移邊界條件,當(dāng)采用基于渦黏假設(shè)的湍流模型來封閉雷諾平均N-S(RANS)方程時(shí),需要對(duì)近壁面流動(dòng)采用特殊的處理方式。與傳統(tǒng)的標(biāo)準(zhǔn)壁面函數(shù)相比,可伸縮壁面函數(shù)對(duì)Y+(第1層網(wǎng)格質(zhì)心到壁面的無量綱距離)的要求并不十分嚴(yán)格,故計(jì)算過程中采用可伸縮壁面函數(shù),將壁面處的物理量與湍流核心區(qū)相應(yīng)物理量聯(lián)系起來。

1.3 物理模型設(shè)計(jì)

物理模型按重力相似準(zhǔn)則設(shè)計(jì),綜合考慮模型水流在阻力平方區(qū)的要求以及模型泵的選擇,確定模型線性比尺λl=10。試驗(yàn)系統(tǒng)采用自循環(huán)開敞式,由取水管、前池、進(jìn)水池、模型泵機(jī)組、出水管路、回水渠、控制及量測(cè)系統(tǒng)等組成。在前池、進(jìn)水池典型斷面位置布置量測(cè)水位的測(cè)壓管,測(cè)壓管連接測(cè)針桶,用水位測(cè)針讀數(shù),水面讀數(shù)誤差為±0.1 mm;在每臺(tái)水泵出水管上安裝超聲波流量計(jì)量測(cè)水泵流量,流量計(jì)精度為0.5級(jí),量測(cè)相對(duì)誤差為±0.5%,同時(shí)每臺(tái)泵分別設(shè)置量水堰槽對(duì)單泵流量進(jìn)行比測(cè),在回水渠設(shè)置量水堰量測(cè)泵站總流量;流速采用多普勒三維剖面點(diǎn)式流速儀(ADV)量測(cè),ADV的精度為±1.0%,流速的相對(duì)誤差為±1.0%,將進(jìn)水池進(jìn)口斷面A作為一典型過流斷面,進(jìn)行流速分布的量測(cè)。流速分布量測(cè)點(diǎn)的布置見圖2,其中水平方向上布置9條施測(cè)垂線,垂線方向上從池底到水面分5層布置。

圖2 流速測(cè)點(diǎn)布置示意圖(單位:mm)

2 泵站前池流態(tài)數(shù)值模擬

2.1 前池典型截面的選取

平面上,選取泵喇叭口入口截面(對(duì)應(yīng)高程-8.90 m)作為典型截面;立面上,選取3號(hào)泵管直徑所在截面作為典型截面。典型截面示意圖見圖3。

圖3 前池典型截面示意圖

2.2 前池進(jìn)水流態(tài)分析

結(jié)合輸水泵站運(yùn)行特點(diǎn),將設(shè)計(jì)高水位的運(yùn)行工況作為研究對(duì)象。該工況為3泵全開,前池水位為1.70 m。以典型截面的速度矢量圖反映無整流措施時(shí)的前池流態(tài)，如圖4所示。

圖4 無整流措施時(shí)的前池流態(tài)

根據(jù)模擬結(jié)果,在平面上,主流居于中間,前池兩側(cè)產(chǎn)生脫壁回流且回流區(qū)的大小與擴(kuò)散段長(zhǎng)度幾乎相等;進(jìn)水池流量分配不均,兩側(cè)進(jìn)水池進(jìn)流較少。在立面上,前池表面有大尺度的回流現(xiàn)象,進(jìn)水池底部流速很大,同時(shí)中上層有多個(gè)漩渦存在。造成前池不良流態(tài)的主要原因是有壓涵管進(jìn)水,主流集中且貼于前池底部,流速較大,雖然該泵站前池?cái)U(kuò)散角較小,擴(kuò)散段較長(zhǎng),水流在前池仍得不到充分的調(diào)整。

3 泵站前池流態(tài)整流措施

3.1 八字形導(dǎo)流墩整流措施

3.1.1 布置形式

導(dǎo)流墩通過其導(dǎo)流作用,減小前池平面擴(kuò)散角,均化水流,能夠有效消除脫壁回流、偏流等不良流態(tài)。該泵站進(jìn)水居中,前池兩側(cè)存在大尺度回流,進(jìn)水池流量分配不均,故導(dǎo)流墩整流措施需在能有效調(diào)整流量分配的同時(shí)消除回流現(xiàn)象。這里參照馮建剛等[13]的研究成果,將八字形導(dǎo)流墩布置在涵管出口處,墩頭與出口距離D=2 m,墩長(zhǎng)L=8 m,兩墩的夾角α=24°,后墩頭正對(duì)第2排結(jié)構(gòu)柱的中心,墩子高度與第2層橫梁底部平齊，如圖5所示。

圖5 八字形導(dǎo)流墩整流措施布置形式(單位：m)

3.1.2 前池進(jìn)水流態(tài)分析

從圖6可知,采用八字形導(dǎo)流墩整流措施后,涵管出口密集居中的主流受導(dǎo)流墩作用,向兩邊分流,并在擴(kuò)散段不斷地向兩側(cè)擴(kuò)散,前池兩側(cè)有了來流的補(bǔ)充,使得回流區(qū)近乎消失,回流現(xiàn)象得到明顯抑制;經(jīng)過較長(zhǎng)的擴(kuò)散段擴(kuò)散后,流量分配得到了一定程度的調(diào)整,3個(gè)進(jìn)水池流量大致相同,但各進(jìn)水池流速分布均勻性較差,主要表現(xiàn)為各流道均存在偏流現(xiàn)象。在立面上,流態(tài)沒有得到明顯的改善,前池表層大尺度回流仍存在,進(jìn)水池底部流速仍較大。

圖6 采用八字形導(dǎo)流墩整流措施后的前池流態(tài)

圖7 所選斷面和各截面示意圖(單位：m)

圖8 采用八字形導(dǎo)流墩整流措施整流前后各進(jìn)水流道所選斷面軸向速度云圖

3.1.3 泵進(jìn)水條件分析

在泵進(jìn)水管前取一斷面,如圖7(a)所示,斷面與泵管軸心在x方向的距離為3 m,分析該斷面軸向速度云圖,直觀地比較泵進(jìn)水狀況的優(yōu)劣;同時(shí),在該斷面上取3個(gè)等高同寬的截面Ⅰ、Ⅱ和Ⅲ,如圖7(b)所示,截面的高度為1.8 m,所在區(qū)域正對(duì)泵喇叭口進(jìn)水區(qū)域,分別計(jì)算各個(gè)截面的法向流速分布均勻度Vnu，計(jì)算公式[14]為

圖8是前池水位為1.70 m時(shí)整流前后各進(jìn)水流道所選斷面軸向速度云圖。整流前,由于主流居中偏下,泵進(jìn)水的最大流速值出現(xiàn)在中間機(jī)組,相對(duì)而言,兩邊機(jī)組進(jìn)水流速要小得多;采用八字形導(dǎo)流墩整流措施后,中間機(jī)組進(jìn)水流速明顯變小,3臺(tái)泵進(jìn)水平均流速大致相等。但各進(jìn)水流道進(jìn)水流速分布不均,偏于某一邊。比較表1中整流前后各截面法向流速分布均勻度,發(fā)現(xiàn)采用八字形導(dǎo)流墩整流措施后,3個(gè)截面的法向流速分布均勻度均下降了20%左右,可見泵喇叭口進(jìn)水流速分布均勻性較整流前變差。

表1 八字形導(dǎo)流墩方案整流前后各截面法向流速分布均勻度

3.2 組合式整流措施

3.2.1 布置形式

在有壓涵管出口布置八字形導(dǎo)流墩,能有效地將密集的主流向前池兩側(cè)分流,消除回流區(qū),且3個(gè)進(jìn)水池流量分配大致相同,但每個(gè)進(jìn)水流道均存在偏流現(xiàn)象,泵喇叭口進(jìn)水流速分布均勻性變差。單一的八字形導(dǎo)流墩措施并不能有效解決城市輸水泵站有壓涵管進(jìn)流存在的不良流態(tài)問題,需另加整流措施。底坎的整流機(jī)理是人為地制造坎后立面漩滾,破壞平面回流,使坎后流態(tài)得到重新調(diào)整,整流效果明顯,同時(shí)底坎結(jié)構(gòu)簡(jiǎn)單,施工簡(jiǎn)便[15-17]。

組合式整流措施是在涵管出口處布置八字形導(dǎo)流墩的基礎(chǔ)上,在前池適當(dāng)位置同時(shí)布置底坎。根據(jù)前池典型截面速度矢量分布和進(jìn)水流道所選斷面軸向速度分布制定出多個(gè)初步布置方案:將底坎布置在位置一、位置二或位置三(每個(gè)位置的間距S=2.5 m),底坎的坎高H設(shè)定為1 m、1.2 m或1.4 m,底坎寬度B設(shè)為0.7 m,如圖9所示。經(jīng)過數(shù)值計(jì)算結(jié)果的分析比較,發(fā)現(xiàn)將底坎布置在位置一時(shí),由于其離八字形導(dǎo)流墩較近,導(dǎo)流墩的擴(kuò)散作用被弱化,主流變得居中;布置在位置三時(shí),坎后漩滾區(qū)域延伸至進(jìn)水流道前,流態(tài)來不及得到充分調(diào)整?？哺逪=1 m時(shí),整流效果較H=1.2 m時(shí)差;坎高H=1.4 m時(shí),挑流現(xiàn)象較為明顯,造成表面流速較大。因此最終確定底坎的布置方案為底坎布置在位置二,寬度B=0.7 m,坎高H=1.2 m。

3.2.2 前池進(jìn)水流態(tài)分析

采取組合式整流措施后,由于底坎的作用,使得主流在坎后發(fā)生立面漩滾,流速得到再分布。從圖10可知,在平面上,坎后有一回流區(qū),這是坎后立面漩滾造成的,過了漩滾區(qū)后,流速分布較均勻;進(jìn)水池流量分配大致相同,且每個(gè)進(jìn)水流道流速分布較均勻,不存在明顯的偏流現(xiàn)象。在立面上,坎前主流收縮,水流偏向自由表面[18],原方案(導(dǎo)流墩方案)前池表面的大尺度回流被擾亂,回流現(xiàn)象得到一定程度的抑制;進(jìn)水池底部流速明顯減小,泵喇叭口前流速分布趨于均勻。

圖10 采用組合式整流措施后前池流態(tài)

3.2.3 泵進(jìn)水條件分析

圖11是無整流措施和采用組合式整流措施后各進(jìn)水流道所選斷面軸向速度云圖。比較后發(fā)現(xiàn),采用組合式整流措施后,3個(gè)進(jìn)水流道流量分配大致相同,斷面的流速最大值明顯減小。從表2可知,采用組合式整流措施后,截面Ⅰ和截面Ⅲ的法向流速分布均勻度提高明顯,截面Ⅱ的法向流速分布均勻度也有略微提高,總體上每臺(tái)泵進(jìn)水流速分布均勻性均較高,且2臺(tái)邊機(jī)組的進(jìn)水條件有了較大改善。

圖11 采用組合式整流措施整流前后各進(jìn)水流道所選斷面軸向速度云圖

有無整流措施Vnu/%截面Ⅰ截面Ⅱ截面Ⅲ無整流措施68．8481．8863．72組合式整流措施84．7084．2688．59

4 試驗(yàn)驗(yàn)證

圖12 物理模型前池流態(tài)

圖13 進(jìn)水流道典型斷面流速計(jì)算值與測(cè)量值

模 型各層測(cè)點(diǎn)平均流速/(m·s-1)1層2層3層4層5層數(shù)值模型0．120．160．140．130．08物理模型0．160．190．120．100．07

為驗(yàn)證數(shù)模計(jì)算結(jié)果的可靠性,進(jìn)行了泵站物理模型試驗(yàn)。通過示蹤粒子來觀測(cè)前池流態(tài),并對(duì)典型過流斷面流速分布進(jìn)行了測(cè)量。選取采用組合式整流措施后3泵全開的工況,前池流態(tài)如圖12所示，進(jìn)水池典型過流斷面x方向(由管涵向進(jìn)水池方向)流速計(jì)算值與測(cè)量值見圖13。對(duì)比圖13(a)(b)發(fā)現(xiàn),計(jì)算值與測(cè)量值總體變化范圍均在0.04～0.23 m/s之間,但計(jì)算值每一層不同立面間流速變化幅度稍大。將典型過流斷面每一層的平均流速計(jì)算值和測(cè)量值統(tǒng)計(jì)入表3,發(fā)現(xiàn)平均流速的最大值均出現(xiàn)在第2層,分別為0.16 m/s和0.19 m/s;最小值均出現(xiàn)在第5層,分別為0.09 m/s和0.07 m/s?？傮w上,典型過流斷面每一層的平均流速計(jì)算值和測(cè)量值大小較相近,變化趨勢(shì)呈現(xiàn)一致性,由此可見,上述數(shù)值計(jì)算方法可行,整流措施優(yōu)化研究結(jié)果可信。

5 結(jié) 論

a. 城市輸水泵站布局較為局促,往往很難按水力條件良好的要求布置,進(jìn)流在前池?cái)U(kuò)散段難以得到充分?jǐn)U散和調(diào)整,前池兩側(cè)存在大尺度回流現(xiàn)象,進(jìn)水池進(jìn)流流量分配不均,水泵進(jìn)水條件較為惡劣。

b. 前池內(nèi)增設(shè)八字形導(dǎo)流墩和底坎的組合式整流措施,能顯著改善前池流態(tài),有效消除前池兩側(cè)回流,坎后流態(tài)重新調(diào)整,各進(jìn)水池流量分配大致相同,流速分布均勻且平穩(wěn),水泵進(jìn)水條件良好。

c. 流速的數(shù)值模擬計(jì)算值與物理模型測(cè)量值大小相近,變化趨勢(shì)一致,表明采用數(shù)值計(jì)算方法研究泵站前池整流是可行的,數(shù)值計(jì)算結(jié)果可信。

[1] 徐輝,張林.側(cè)向進(jìn)水泵站前池整流技術(shù)研究綜述[J].水利水電科技進(jìn)展,2008,28(6):84-88.(XU Hui,ZHANG Lin.Review of flow pattern modifying technology for pump sump with side-inlet[J].Advances in Science and Technology of Water Resources,2008,28(6):84-88.(in Chinese))

[2] 羅燦,錢均,劉超,等.非對(duì)稱式閘站結(jié)合式泵站前池導(dǎo)流墩整流模擬及試驗(yàn)驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(7):100-108.(LUO Can, QIAN Jun, LIU Chao, et al. Numerical simulation and test verification on diversion pier rectifying flow in forebay of pumping station for asymmetric combined sluice-pump station project[J].Transactions of the Chinese Society of Agricultural Engineering,2015,31(7):100-108.(in Chinese))

[3] 成立,祁衛(wèi)軍,羅燦,等.Y形導(dǎo)流墩幾何參數(shù)對(duì)泵站前池整流效果的影響[J].水利水電科技進(jìn)展,2014,34(1):68-72.(CHENG Li,QI Weijun,LUO Can, et al. Effect of geometric parameters of Y-shaped diversion piers on flow pattern in forebay of pumping station[J].Advances in Science and Technology of Water Resources,2014,34(1):68-72.(in Chinese))

[4] 戎貴文,袁岳,肖柏青,等.局部斷面收縮的潰壩水流三維數(shù)值模擬[J].水利水電科技進(jìn)展,2015,35(6):20-25.(RONG Guiwen,YUAN Yue,XIAO Baiqing, et al.Three-dimensional numerical simulation of dam-break water flow with local cross-section reduction[J].Advances in Science and Technology of Water Resources,2015,35(6):20-25.(in Chinese))

[5] 周大慶,郭優(yōu),姜德政.基于槳葉調(diào)節(jié)的軸流轉(zhuǎn)槳式水輪機(jī)模型飛逸數(shù)值模擬[J].水利水電科技進(jìn)展,2016,36(4):13-19.(ZHOU Daqing,GUO You,JIANG Dezheng.Numerical simulation of runaway transients of kaplan turbine model based on blade regulation[J].Advances in Science and Technology of Water Resources,2016,36(4):13-19.(in Chinese))

[6] 向小華,吳曉玲,牛帥,等.基于顯式有限體積法的一維河網(wǎng)模型[J].水利水電科技進(jìn)展,2015,35(4):6-9.(XIANG Xiaohua,WU Xiaoling,NIU Shuai, et al.Construction of one-dimensional river network model based on explicit finite volume method[J].Advances in Science and Technology of Water Resources,2015,35(4):6-9.(in Chinese))

[7] 田艷,楊平,魯俊.城市排水泵站前池流態(tài)的數(shù)值模擬研究[J].水泵技術(shù),2008(4):24-26.(TIAN Yan, YANG Ping, LU Jun. Numerical simulation of flow patterns in forebays of the urban drainage pumping stations[J].Pump Technology,2008(4):24-26.(in Chinese))

[8] 魯俊.城市排水泵站前池整流技術(shù)CFD研究[D].南京:河海大學(xué),2006.

[9] 馮建剛,李杰.大型城市水源泵站前池流態(tài)及改善措施試驗(yàn)[J].水利水電科技進(jìn)展,2010,30(2):70-74.(FENG Jiangang,LI Jie.Flow patterns in forebays of large urban water source pumping stations and countermeasures[J].Advances in Science and Technology of Water Resources,2010,30(2):70-74.(in Chinese))

[10] 黃繼宏,王曉升,張曉毅.城市排水泵站前池導(dǎo)流墩整流措施的數(shù)值模擬[J].江蘇農(nóng)業(yè)科學(xué),2013,41(5):360-362.(HUANG Ji Hong, WANG Xiao Shen, ZHANG Xiao Yi. Numerical simulation of rectification measures of diversion pier in the forebay of urban drainage pumping stations[J].Jiangsu Agricultural Sciences,2013,41(5):360-362.(in Chinese))

[11] 羅燦,成立,劉超.泵站正向進(jìn)水前池底坎整流機(jī)理數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2014,32(5):393-398.(LUO Can, CHENG Li, LIU Chao. Numerical simulation of mechanism for sill rectifying flow in pumping station intake[J].Journal of Drainage and Irrigation Machinery Engineering,2014,32(5):393-398.(in Chinese))

[12] 王福軍.計(jì)算流體動(dòng)力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2007:114-126.

[13] 馮建剛,王曉升,佟宏偉,等.大型城市供水泵站前池流態(tài)改善措施研究[J].給水排水,2010,36(11):51-54.(FENG Jiangang,WANG Xiaosheng,TONG Hongwei, et al. Improving of the flow patterns in forebays of the large urban water supply pumping stations[J].Water & Wastewater Engineering,2010,36(11):51-54.(in Chinese))

[14] 李坦,靳世平,黃素逸,等.流場(chǎng)速度分布均勻性評(píng)價(jià)指標(biāo)比較與應(yīng)用研究[J].熱力發(fā)電,2013,42(11):60-63.(LI Tan,JIN Shiping,HUANG Suyi, et al. Evaluation indices of flow velocity distribution uniformity: comparison and application[J].Thermal Power Generation,2013,42(11):60-63.(in Chinese))

[15] 周大慶,劉躍飛.基于VOF模型的軸流泵機(jī)組起動(dòng)過程數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2016,34(4):307-312.(ZHOU Daqing,LIU Yuefei.Numerical simulation of axial pump unit startup process using VOF model[J].Journal of Drainage and Irrigation Machinery Engineering,2016,34(4):307-312.(in Chinese))

[16] 張琳,施衛(wèi)東,張德勝,等.基于大渦模擬的混流泵不穩(wěn)定流動(dòng)數(shù)值研究[J].排灌機(jī)械工程學(xué)報(bào),2017,35(4):303-308.(ZHANG Lin，SHI Weidong，ZHANG Desheng，et al. Analysis of instability flow on mixed-flow pump based on LES[J].Journal of Drainage and Irrigation Machinery Engineering,2017,35(4):303-308.(in Chinese))

[17] 成立,劉超,周濟(jì)人,等.泵站前池底壩整流數(shù)值模擬研究[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,29(3):42-45.(CHENG Li, LIU Chao, ZHOU Jiren, et al.Numerical simulation of sill flows in the forebay of pumping station[J].Journal of Hohai University(Natural Sciences),2001,29(3):42-45.(in Chinese))

[18] 馮旭松.泵站前池底坎整流及坎后流動(dòng)分析[D].揚(yáng)州:揚(yáng)州大學(xué),1996.