周蘇芬,易子靖,閆旭峰,楊青遠(yuǎn),2,王協(xié)康

(1.四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,四川成都 610065; 2.香港城市大學(xué)土木與建筑工程系,香港 310027)

山區(qū)寬窄相間河道平面二維水流數(shù)值模擬

周蘇芬1,易子靖1,閆旭峰1,楊青遠(yuǎn)1,2,王協(xié)康1

(1.四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,四川成都 610065; 2.香港城市大學(xué)土木與建筑工程系,香港 310027)

針對山區(qū)河道受地質(zhì)條件及河床演化的影響,其平面形態(tài)常呈現(xiàn)寬窄相間的漸變形態(tài)特征,采用貼體坐標(biāo)法剖分正交貼體曲線網(wǎng)格求解域,用有限體積法離散控制方程,以SIMPLER計(jì)算程式求解耦合方程,建立山區(qū)寬窄相間河道平面二維水流數(shù)學(xué)模型。將模型計(jì)算值與寬窄相間水槽水流試驗(yàn)實(shí)測值及水動力學(xué)SMS軟件計(jì)算結(jié)果進(jìn)行比較,結(jié)果表明:該模型能夠較好地模擬寬窄相間河道平面二維水流運(yùn)動特性,寬窄相間河道河寬漸變將引起水深和流速的變化,在河道最窄處水深降到最低,流速最大;在河道最寬處水深達(dá)到峰值,流速最小。

寬窄相間河道;平面二維數(shù)值模擬;貼體坐標(biāo);有限體積法;SIMPLER算法

山區(qū)河流在峽谷地區(qū)受地質(zhì)條件約束,經(jīng)過極長時(shí)間調(diào)整常形成寬窄相間形態(tài),如岷江金馬河段最寬處為1010m,而其他地方河寬多在300～600 m之間。寬窄相間的平面形態(tài)顯著影響河道水沙運(yùn)動及其演化特性。Lucy等[1]基于野外調(diào)查認(rèn)為河道寬度變化是影響河床演化的關(guān)鍵因素之一。Singha等[2]基于試驗(yàn)指出,河道形態(tài)變化將重新調(diào)整水流紊動能的分布特性。Armellini等[3]基于數(shù)值計(jì)算得出,受岸灘形態(tài)的影響,水流結(jié)構(gòu)出現(xiàn)極為明顯的分離區(qū)。Paiement-Paradis等[4]通過試驗(yàn)表明河道寬窄變化引起的水流減速或加速運(yùn)動對泥沙推移質(zhì)輸移有重要影響。閆旭峰等[5]基于室內(nèi)試驗(yàn)指出漸變河道形態(tài)顯著影響寬窄相間區(qū)域的局部水頭損失。

目前,平面二維水流數(shù)學(xué)模型得到了較成熟的發(fā)展,為利用數(shù)值模型研究寬窄相間河道的水流特性奠定了較好的基礎(chǔ)。早期有荷蘭德爾夫特水力研究所的RIVMOR 2D模型,美國愛荷華大學(xué)水利所的MOBED2模型,德國卡爾斯魯厄大學(xué)的FAST2D模型及美國陸軍工程兵水利工程實(shí)驗(yàn)室和楊百翰大學(xué)等合作開發(fā)的水動力學(xué)SMS模型及丹麥水利所的MIKE21模型。Arruda等[6-8]將IB法(Immersed boundary method)應(yīng)用于平面二維水流數(shù)值模擬計(jì)算中,并取得了較好的結(jié)果。李大鳴等[9]考慮上、下游洪水邊界的控制條件及糙率隨河道流量變化的調(diào)整,提出了河道洪水演進(jìn)的二維水流數(shù)學(xué)模型的計(jì)算方法。張細(xì)兵等[10]采用質(zhì)量集中的簡化處理和預(yù)估校正的時(shí)間推算法,較好地解決了平面二維水流泥沙數(shù)值求解的速度問題。張大偉等[11]通過有限體積法建立了反映大壩瞬時(shí)潰決水流特點(diǎn)的數(shù)學(xué)模型,并以Malpasset水庫對模型進(jìn)行驗(yàn)證。葉如意等[12]利用平面二維水動力學(xué)模型對曹娥江至寧波的引水工程進(jìn)行模擬,同時(shí)對模型進(jìn)行驗(yàn)證和優(yōu)化。方神光等[13]結(jié)合純隱格式的混合有限分析法及干濕判別標(biāo)準(zhǔn)與凍結(jié)法建立了二維水流數(shù)學(xué)模型并將其應(yīng)用于實(shí)際工程中,取得了較準(zhǔn)確的計(jì)算結(jié)果。王君超[14]利用有限體積法建立了平面二維水沙數(shù)值模型,計(jì)算結(jié)果能較好地反映天然河床的變形規(guī)律。為了更深入地探討寬窄相間河道的水流特性,本文在前人研究的基礎(chǔ)上,以寬窄相間水流運(yùn)動為研究對象,通過室內(nèi)水槽試驗(yàn)及數(shù)值模擬,進(jìn)一步研究平面二維水流模型以模擬山區(qū)河道水流的運(yùn)動特性。

1 正交貼體曲線網(wǎng)格的生成

采用橢圓形方程法生成二維正交貼體曲線網(wǎng)格,為避免產(chǎn)生波狀速度場,動量方程的離散采用交錯(cuò)網(wǎng)格體系,即將流速布置在控制體積的交界面上,與水深、水位等錯(cuò)開半個(gè)網(wǎng)格。

物理平面上的不規(guī)則區(qū)域與計(jì)算平面上的規(guī)則區(qū)域之間的坐標(biāo)變換滿足Possion方程:

式中:x、y為笛卡兒坐標(biāo)系下求解域網(wǎng)格節(jié)點(diǎn)坐標(biāo);ξ、η為貼體坐標(biāo)系下求解域網(wǎng)格節(jié)點(diǎn)坐標(biāo);P、Q為調(diào)節(jié)因子,即網(wǎng)格密度控制函數(shù),其作用是調(diào)節(jié)網(wǎng)格線位置的變化,使網(wǎng)格的疏密根據(jù)需要分布并可使變換后的網(wǎng)格曲線更好地?cái)M合邊界,保證正交,P、Q采用Thompson等[15]建議的函數(shù)表達(dá)式。

通過求解Possion方程生成試驗(yàn)漸變河段網(wǎng)格,見圖1。

2 山區(qū)河道平面二維數(shù)值模型

2.1 控制方程

針對山區(qū)河道漸變特點(diǎn),采用貼體坐標(biāo)法對生成的二維正交貼體曲線網(wǎng)格求解域進(jìn)行剖分,利用有限體積法離散控制方程,以SIMPLER計(jì)算程式求解耦合方程。

將水流連續(xù)方程和動量方程寫成統(tǒng)一的偏微分方程:

圖1 試驗(yàn)漸變河段網(wǎng)格

式中:Cξ、Cη為拉梅系數(shù);H為水深;Φ為通用變量;t為時(shí)間;u、v分別為x、y方向的流速;Γ為擴(kuò)散系數(shù);S為源項(xiàng)。

對S進(jìn)行負(fù)坡線性化處理,即S=SC+ΦOSO(SC為源項(xiàng)的常數(shù)部分;SO為系數(shù);ΦO為中心點(diǎn)的變量)。根據(jù)有限體積法,將通用偏微分方程(式(3))在控制體積內(nèi)積分,得到離散方程如下[16]:

式中:a為系數(shù);Φ為節(jié)點(diǎn)處變量;b為相鄰節(jié)點(diǎn)變量與源項(xiàng)產(chǎn)生的常數(shù);下標(biāo)O、E、W、N、S表示節(jié)點(diǎn)位置,O表示中心點(diǎn),E、W、N、S分別表示東、西、北、南節(jié)點(diǎn)。

離散格式采用的是由Patankar[17]提出的冪函數(shù)格式,則表示如下:

式中:Pa為派克里特?cái)?shù),為對流強(qiáng)度與擴(kuò)散強(qiáng)度之比。

計(jì)算初始流速場為零,初始水位場根據(jù)已知進(jìn)出口斷面初始水位內(nèi)插出求解域中各點(diǎn)的初始水位,斷面不考慮橫比降。計(jì)算河段進(jìn)口邊界給定流量,下游出口邊界給定各點(diǎn)水位過程線。

2.2 模型的建立

基于Visual C++語言建立平面二維水流數(shù)值模型,模型涉及的主要變量和參數(shù)有各點(diǎn)水流流速u、v,床底高程zb,水位z,初始流速u0、v0,初始水位z0,各網(wǎng)格點(diǎn)x、y、z坐標(biāo),正交曲線坐標(biāo)系中的拉梅系數(shù)Cξ、Cη及河床的糙率n。

3 模型的驗(yàn)證與計(jì)算

通過模型計(jì)算、試驗(yàn)測量以及SMS軟件[5]計(jì)算結(jié)果的比較來驗(yàn)證模型的可靠性,其中,SMS軟件計(jì)算數(shù)據(jù)采用文獻(xiàn)[5]中的計(jì)算成果。

3.1 SMS軟件

SMS軟件包括RMA2、TABS-MD、AD-CIRC、CGWAVE、STWAVE、HIHEL等計(jì)算模塊。RMA2是一個(gè)沿水深積分后平均的二維有限元水動力數(shù)值模型,可以用于計(jì)算具有自由表面緩流的水面高程、水深、平面流速的二維流場,適用于恒定流和非恒定流。該模塊根據(jù)Galerkin的殘值權(quán)重理論的有限元法及Newton-Raphson迭代法求解,再根據(jù)Navier-Stokes方程推求水深,其中綜合糙率用曼寧公式計(jì)算。RMA2基于靜水條件假定運(yùn)行,并忽略重力方向的平均加速度,它僅局限于水平面上的二維數(shù)值模擬。

3.2 試驗(yàn)工況

寬窄相間河道模型水槽位于四川大學(xué)山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室模型實(shí)驗(yàn)場A區(qū)的漸變水槽自循環(huán)系統(tǒng)中(圖2)。水槽內(nèi)壁采用水泥光滑抹面,糙率為0.013。水槽全長28 m,坡度為0.2％。水槽進(jìn)口斷面處為磚砌花墻,能起到平緩水流的作用;花墻距離漸變起始點(diǎn)5.5m;最后一個(gè)最窄斷面至水槽末端7.5m;水槽最寬為1.4 m,最窄為0.6 m;水槽邊墻高0.5 m,邊線為正弦曲線。水流方向如圖2所示。試驗(yàn)采用美國SonTec公司的16M聲學(xué)多普勒流速儀(ADV)測量三維瞬時(shí)流速,整個(gè)水槽共設(shè)37個(gè)測量斷面(圖3),在河寬擴(kuò)展處水流比較平穩(wěn),斷面間距較大,在縮窄處測量斷面較密集,每個(gè)斷面隨槽寬的漸變設(shè)置不同數(shù)量的測量垂線。試驗(yàn)有兩種工況:工況1流量為0.112 m3/s,下游為自由出流;工況2流量為0.100 m3/s,下游為自由出流。

圖2 漸變水槽自循環(huán)系統(tǒng)(單位:m)

3.3 模型驗(yàn)證

將工況1實(shí)測斷面中軸線上的平均流速和水深測量結(jié)果分別與本文模型計(jì)算結(jié)果及SMS軟件計(jì)算結(jié)果進(jìn)行比較,見圖4和圖5(圖中ˉu為平均流速)。從圖4和圖5可以看出在進(jìn)出口段,本文模型計(jì)算結(jié)果比SMS軟件計(jì)算結(jié)果與實(shí)測結(jié)果吻合得更好。總體來說,平均流速和平均水深的模型計(jì)算值與實(shí)測值及SMS軟件計(jì)算值較為一致。由此可知,該模型對漸變水流特性的模擬效果較好。

圖3 試驗(yàn)測量斷面分布

圖4 實(shí)測斷面中軸線沿程平均流速比較

圖5 實(shí)測斷面中軸線沿程水深比較

圖6 縮窄斷面速度場比較

圖6 為典型縮窄斷面速度場計(jì)算值與實(shí)測值的比較。由于所測水槽水流對稱,故僅測量水槽一半寬度的流速。造成圖中測量點(diǎn)與計(jì)算點(diǎn)坐標(biāo)之間差異的主要原因是計(jì)算網(wǎng)格的疏密程度與測量斷面疏密程度不一致,如果計(jì)算網(wǎng)格足夠密,將能找到與測量斷面一致的計(jì)算斷面。從圖中可以看出,計(jì)算速度場與實(shí)測速度場的分布規(guī)律及大小基本相同,即流速順著邊界方向變化而變化,在河寬縮窄段速度場出現(xiàn)收縮,流速增大;在河寬漸擴(kuò)段場出現(xiàn)擴(kuò)散,流速減小,這些現(xiàn)象進(jìn)一步驗(yàn)證了模型的可靠性。

3.4 模型計(jì)算

利用數(shù)值模型模擬工況2的水流特性,計(jì)算結(jié)果與實(shí)測值的比較如圖7～9所示。從圖7可知,流速和水深會隨著河寬變化而發(fā)生相應(yīng)的變化,在最窄斷面處水深最小,流速最大。從圖8可以看出,工況2最窄斷面處Fr最大,并且其值大于1.0,隨著河道逐漸展寬,Fr逐漸減小,其值小于1.0,可見在河寬縮窄段水流處于急流狀態(tài),在河寬漸闊段水流處于緩流狀態(tài)。由此可以說明,河寬的改變會影響水流的流態(tài)。從流速等值線(圖9)也可以看出,在最窄斷面處流速最大,但在河段最寬斷面靠近邊壁附近對稱地出現(xiàn)了相對于原斷面較低流速的區(qū)域,此時(shí)形成了回流區(qū)。

圖7 斷面中軸線沿程平均流速與水深

圖8 Fr沿程變化曲線

圖9 流速等值線(單位:m/s)

4 結(jié) 語

針對山區(qū)河道寬窄相間的特點(diǎn),采用貼體坐標(biāo)法將河道不規(guī)則的邊界區(qū)域轉(zhuǎn)化成規(guī)則的計(jì)算區(qū)域,并利用SIMPLER計(jì)算程式求解耦合方程,基于Visual C++語言建立了山區(qū)寬窄相間河道平面二維水流數(shù)值模型。模型計(jì)算值與實(shí)測值及SMS軟件計(jì)算值的比較表明:①寬窄相間河道河寬漸變將引起水深和流速的變化,在河道最窄處水深降到最低,流速最大;在河道最寬處水深達(dá)到峰值,流速最小。此外,水流運(yùn)動方向隨著邊界的收縮和擴(kuò)散將出現(xiàn)相應(yīng)的偏移。②本文建立的平面二維水流模型對于模擬寬窄相間河道水流特性具有較高的可靠性。③該模型較好地揭示了寬窄相間河道水流流態(tài)的變化特征及其回流區(qū)特性。

[1]LUCY C,QUINE T A,NICHOLAS A.An experimental investigation of autogenic behavior during alluvial fan evolution[J].Geomorphology,2010,115(3/4):278-285.

[2]SINGHA A,BALACHANDAR R.Structure of wake of a sharp-edged bluff body in a shallow channel flow[J]. Journal of Fluids and Structures,2011,27(2):233-249.

[3]ARMELLINI A,CASARSAL,GIANNATTASIOP. Separated flow structures around a cylindrical obstacle in a narrow channel[J].Experimental Thermal and Fluid Science,2009,33(4):604-619.

[4]PAIEMENT-PARADIS G,MARQUIS G,ROY A.Effects of turbulence on the transport of individual particles as bedload inagravel-bedriver[J].EarthSurface Processes and Landforms,2011,36(1):107-116.

[5]閆旭鋒,易子靖,劉同宦,等.漸變河道水流結(jié)構(gòu)及局部水頭損失特性分析[J].長江科學(xué)院院報(bào),2011,28 (9):1-5.(YAN Xufeng,YI Zijing,LIU Tonghuan,et al.Flow structure and characteristics of local head loss in transition channel[J].Journal of Yangtze River Scientific Research Institute,2011,28(9):1-5.(in Chinese))

[6]ARRUDA J M,LIMA E-SILVA A L F,NETO A S,et al. Numerical simulation of channel flows using the immersed boundary method[C]//IDELSOHN S R,SONZOGNI V E,CARDONAA.South-AmericanCongresson Computational Mechanics.Santa Fe-Paraná,Argentina: [s.n.],2002.

[7]LIMAE-SILVA A F E,NETO A S,DAMASCENO J J R. Numerical simulation of two-dimensional flows over a circular cylinder using the immersed boundary method [J].Journal of Computational Physics,2003,189(2): 351-370.

[8]LIMAE-SILVA A F E,SILVAAR R D,NETO A S. Numerical simulation of two-dimensional complex flows around bluff bodies using the immersed boundary[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering,2007,29(4):379-387.

[9]李大鳴,陳虹,李世森.河道洪水演進(jìn)的二維水流數(shù)學(xué)模型[J].天津大學(xué)學(xué)報(bào),1998,31(4):439-446.(LI Daming,CHEN Hong,LI Shishen.A 2-D numerical model of propelling flood in the river[J].Journal of Tianjin University,1998,31(4):439-446.(in Chinese))

[10]張細(xì)兵,殷瑞蘭.平面二維水流泥沙數(shù)值模擬[J].水科學(xué)進(jìn)展,2002,13(6):665-669.(ZHANG Xibing,YIN Ruilan.Planar 2-Dflowandsedimentmathematical modeling[J].Advances in Water Science,2002,13(6): 665-669.(in Chinese))

[11]張大偉,程曉陶,黃金池.大壩瞬時(shí)潰決水流數(shù)值模擬:以Malpasset水庫為例[J].水利水電科技進(jìn)展,2009, 29(5):1-4.(ZHANG Dawei,CHENG Xiaotao,Huang Jinchi.Numerical flow modeling of instant dam failure:a case study of the Malpasset Reservoir[J].Advances in Science and Technology of Water Resources,2009,29 (5):1-4.(in Chinese))

[12]葉如意,王攀.曹娥江至寧波引水工程平面二維水流數(shù)值模擬[J].水資源保護(hù),2010,26(5):24-28.(YE Ruyi,WANG Pan.Two-dimensional numerical simulation of flow for water diversion project from Cao'e River to Ningbo[J].Water Resources Protetion,2010,26(5):24-28.(in Chinese))

[13]方神光,黃勝偉.純隱格式的混合有限分析法在廣雅橋水流數(shù)值模擬中的應(yīng)用[J].水利水電科技進(jìn)展,2010, 30(6):53-57.(FANG Shenguang,HUANG Shengwei. Application of implicit schemehybridfiniteanalytic method in numerical simulation of flow effect by Guangya Bridge[J].Advances in Science and Technology of Water Resources,2010,30(6):53-57.(in Chinese))

[14]王君超.平面二維河流數(shù)值模擬方法及其應(yīng)用[J].中國科教創(chuàng)新導(dǎo)刊,2010(11):93-94.(WANG Junchao.A two-dimensional numericalsimulationmethodandits application in the river[J].China Education Innovation Herald,2010(11):93-94.(in Chinese))

[15]THOMPSON J F,WARSI Z U A,WAYNE C M. Numerical grid generation:foundation and applications [M].New York:Elsevier Science Publishing Co.,1985.

[16]陳月華.干支流交匯河段水流特性計(jì)算研究[D].南京:南京水利科學(xué)研究院,2007.

[17]PATANKAR SV.Acalculationprocedurefortwodimensionalellipticsituations[J].NumericalHeat Transfer,1981,4(4):409-425.

Two-dimensional numerical simulation of flows in wide and narrow alternated channels in mountainous areas

// ZHOU Sufen1,YI Zijing1,YAN Xufeng1,YANG Qingyuan1,2,WANG Xiekang1(1.State Key Laboratory of Hydraulics and Mountain River Engineering,Sichuan University,Chengdu610065,China;2.Department of Civil and Architectural Engineering,City University of Hong Kong,Hong Kong310027,China)

The plane configuration of river channels in mountainous areas often exhibits gradual characteristics of wide and narrow alternated variations under the effect of geological conditions and riverbed evolution.A two-dimensional mathematical model for wide and narrow alternated river channels in mountainous areas was established by dividing computational domains using the body-fitted coordinates,discretizing governing equations using the finite volume method and solving equation by the SIMPLER algorithm.The model simulations were compared with the observations of flume experiment and the calculated ones by use of the hydrodynamic surface water modeling system(SMS).The results show that the proposed model can satisfactorily depict the two-dimensional flow motion characteristics of wide and narrow alternated river channels.The gradual wide variations of river channels will lead to change of their water depth and velocity. The flow velocity at the narrowest section reaches is the largest and the depth is the smallest.The flow velocity at the widest section is the smallest and the depth is the largest.

wide and narrow alternated channel;two-dimensional numerical simulation;body-fitted coordinates;finite volume method;SIMPLER algorithm

10.3880/j.issn.10067647.2013.01.005

TV131.3

A

10067647(2013)01002205

2012-03-06 編輯:駱超)

國家自然科學(xué)基金(41171016);四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室開放基金(1008)

周蘇芬(1989—),女,江西新余人,碩士研究生,主要從事水力學(xué)及河流動力學(xué)研究。E-mail:512595913@qq.com

王協(xié)康(1970—),男,江西吉安人,研究員,主要從事水力學(xué)及河流動力學(xué)研究。E-mail:wangxiekang@scu.edu.cn