祝 龍， 周冬卉， 李 云， 宣國祥， 王曉剛

(南京水利科學(xué)研究院 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京 210029)

土壩潰決跌坎水流水動力特性數(shù)值模擬

祝 龍， 周冬卉， 李 云， 宣國祥， 王曉剛

(南京水利科學(xué)研究院 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京 210029)

黏性土壩漫頂潰決涉及多學(xué)科交叉，過程極其復(fù)雜，盡管國內(nèi)外大量物理模型試驗(yàn)成果表明其潰決多以“跌坎式”溯源沖蝕為主要特征，然而對該沖蝕發(fā)展形式下的水流-壩體微觀作用機(jī)制尚不清晰。水流作為漫頂潰決的沖刷主動力，對壩體潰決發(fā)展起著主導(dǎo)性作用，采用RNGk-ε紊流模型和VOF自由液面捕捉技術(shù)針對黏性土壩漫頂潰決代表性水流結(jié)構(gòu)——潰決跌坎水流開展了三維數(shù)值模擬研究，對跌坎水流的水流結(jié)構(gòu)、流態(tài)、水力特性指標(biāo)等進(jìn)行了細(xì)致分析，揭示了不同工況下壩體跌坎上的剪切應(yīng)力、流速分布規(guī)律，進(jìn)而從水動力學(xué)的角度對壩面進(jìn)行受力分析，初步推斷了黏性土壩漫頂潰決過程中各級跌坎的主要合并方式為“臺階水平面刷深下切”。研究成果為進(jìn)一步掌握黏性土壩漫頂潰決發(fā)展演變機(jī)理提供了理論基礎(chǔ)。

黏性土壩； 漫頂； 跌坎水流； 水動力特性； 數(shù)值模擬

截至2011年底，我國已建成各類水庫大壩98 002座[1]，是當(dāng)今世界擁有水庫大壩最多的國家。這些堤壩工程極大地促進(jìn)了我國社會經(jīng)濟(jì)發(fā)展，然而一旦發(fā)生潰壩事故，將造成不可估量的社會經(jīng)濟(jì)損失及慘重的人員傷亡。近年來，國內(nèi)外專家學(xué)者對潰壩問題給予了高度重視并開展了一系列科研工作[2-5]，如：美國1996年啟動的國家大壩安全計(jì)劃(NDSP)，歐盟自1998年起連續(xù)啟動的CADAM，IMPACT以及FLOODsite項(xiàng)目，芬蘭1999年啟動實(shí)施的RESCDAM項(xiàng)目等，并在試驗(yàn)的基礎(chǔ)上提出了“Headcut”潰壩溯源沖蝕發(fā)展機(jī)理。結(jié)合我國黏性土壩自身的特點(diǎn)，南京水利科學(xué)研究院開展了大量現(xiàn)場原型試驗(yàn)和室內(nèi)小尺度物理模型試驗(yàn)研究[6]，對潰壩機(jī)理進(jìn)行了有效發(fā)展和完善。相關(guān)成果表明，土壩漫頂時(shí)，在相當(dāng)長的時(shí)間內(nèi)并不會形成貫通壩體上下游的大潰口，而在壩體下游坡面的中上部易形成首個(gè)沖刷坑，造成初始沖刷并使得水流結(jié)構(gòu)發(fā)生相應(yīng)的改變，之后在其下部形成第二級沖坑，各沖坑以“跌坎剪剝式(Headcut)”沖蝕方式向上游逐步推進(jìn)。Powledge[7]和李云[8]等亦分別從理論分析和數(shù)值模擬的角度闡述了類似觀點(diǎn)，他們認(rèn)為潰壩發(fā)生時(shí)，初始沖坑可能發(fā)生在下游坡面的任意位置，尤其壩肩附近由于壩體體型突變而導(dǎo)致的漫頂水流挑流落點(diǎn)處。

由此可見，潰壩跌坎水流是黏性土壩漫頂潰決過程中的典型代表性水流型式之一，然而以往對于潰決水流的模擬大多假設(shè)壩體瞬間全部潰決或者是瞬間局部潰決，且研究內(nèi)容主要側(cè)重于土壩潰決峰值流量預(yù)測以及下游洪水演進(jìn)問題[9-12],對潰壩過程中壩體區(qū)域的水流特性缺乏系統(tǒng)的研究；此外，黏性土壩漫頂潰決涉及多學(xué)科交叉，過程極其復(fù)雜，考慮水土耦合從整體上進(jìn)行統(tǒng)一模擬具有極大的難度和不確定性。鑒于此，本文針對潰壩發(fā)展過程中的典型潰口形態(tài)，采用三維數(shù)值模擬技術(shù)對其水動力特性進(jìn)行研究，全面揭示水流沖蝕主動力，進(jìn)而從該角度對黏性土壩溯源蝕退機(jī)理的發(fā)展和完善提供理論基礎(chǔ)。

1 數(shù)學(xué)模型的建立

1.1 水流控制方程組

計(jì)算流體動力學(xué)無論具體采用何種形式，都建立在流體力學(xué)基本控制方程——連續(xù)性方程、動量方程和能量方程的基礎(chǔ)之上，必須遵循3個(gè)基本物理學(xué)原理：質(zhì)量守恒、牛頓第二定律以及能量守恒。在黏性土壩漫頂潰決的整個(gè)發(fā)生發(fā)展過程中，機(jī)械能與內(nèi)能的相互轉(zhuǎn)化可以忽略不計(jì)，因此在相關(guān)水力學(xué)問題研究中，主要滿足連續(xù)性方程和動量守恒方程即可。本文在建模時(shí)采用“FAVOR”技術(shù)，即在連續(xù)性方程和動量守恒方程中加入含有面積和體積分?jǐn)?shù)的參數(shù)，表示如下。

(1)

(2)

式中：ρ為流體密度；p為壓力；t為時(shí)間；u，v，w為速度矢量在x，y，z方向上的分量；Ax，Ay，Az表示x，y，z方向可以流動的面積分?jǐn)?shù)；VF為可流動的體積分?jǐn)?shù)；G為重力項(xiàng)；fx，fy，fz為3個(gè)方向的黏滯力。

1.2 RNG k-ε紊流模型

RNGk-ε模型是針對標(biāo)準(zhǔn)k-ε模型在模擬強(qiáng)旋流或帶有彎曲壁面的流動時(shí)出現(xiàn)失真現(xiàn)象而提出的，在模擬沖擊射流、分離流、二次流等復(fù)雜流動時(shí)具有較大的優(yōu)勢，考慮到潰壩跌坎沖射水流的強(qiáng)紊動特性，本文在數(shù)值計(jì)算中采用該模型對Navier-Stokes控制方程進(jìn)行封閉。同樣在模型中加入表征相應(yīng)體積、面積分?jǐn)?shù)的參數(shù)，則直角坐標(biāo)系下的湍動能kT和耗散率εT的輸運(yùn)方程為：

(3)

(4)

式中：kT為紊動能；PT為由速度梯度引起的紊動能產(chǎn)生項(xiàng)；GT為浮力引起的紊動能產(chǎn)生項(xiàng)；DT和Dε為紊動擴(kuò)散項(xiàng)；εT為紊動耗散率；C1s,C2s和C3s均為無量綱系數(shù)，在模型中，C1s取1.44，C2s取1.92，C3s取0.2。

速度梯度引起的紊動項(xiàng)PT由式計(jì)算，如下：

(5)

式中：Cs為參數(shù)，默認(rèn)取1；μ為動力黏滯系數(shù)。

浮力引起的紊動項(xiàng)GT由式計(jì)算，如下：

(6)

式中：p為壓力;Co為紊動參數(shù)。

擴(kuò)散項(xiàng)DT和Dε的表達(dá)式分別如下所示：

(7)

(8)

式中：υk和υε為紊動擴(kuò)散系數(shù)，根據(jù)運(yùn)動黏滯系數(shù)υT而定，其中υT可表示為：

(9)

式中：Cu默認(rèn)取0.085。計(jì)算過程中對εT進(jìn)行嚴(yán)格控制，避免造成與實(shí)際不符的巨大紊動能耗散。

1.3 數(shù)值計(jì)算方法

采用有限差分法對控制方程進(jìn)行離散，在空間上形成三維交錯(cuò)網(wǎng)格，將速度和面積分?jǐn)?shù)定義在網(wǎng)格邊界面中心，其他變量如壓強(qiáng)、密度、流體黏度等定義在控制體中心。對流項(xiàng)采用一階迎風(fēng)離散格式，黏性項(xiàng)采用標(biāo)準(zhǔn)的中心差分離散格式。將控制方程離散為代數(shù)方程之后，首先引入一個(gè)中間速度量，不考慮下一時(shí)刻壓力場對速度場產(chǎn)生的影響，引入當(dāng)前時(shí)刻的壓力修正值，通過求解動量方程獲得中間速度，再將由動量方程離散出來的相關(guān)關(guān)系代入連續(xù)方程中，生成含有壓力修正值的泊松方程。

GMRES算法計(jì)算精度高，收斂速度快且不易發(fā)散，因此本文選用該算法進(jìn)行數(shù)值求解計(jì)算。壁面采用無滑移邊界條件，近壁區(qū)域采用數(shù)學(xué)函數(shù)近似模擬，且由速度梯度產(chǎn)生的黏性剪切應(yīng)力通過加入面積分?jǐn)?shù)參數(shù)來考量，自由液面的追蹤采用VOF(Volume of Fluid)方法來實(shí)現(xiàn)。

1.4 計(jì)算區(qū)域及邊界條件設(shè)置

潰壩發(fā)展過程中典型潰口形態(tài)見圖1。數(shù)學(xué)模型針對圖1所示的黏性土壩典型潰決水流開展研究，基于對眾多試驗(yàn)現(xiàn)象的觀測，將其概化如圖2所示。

圖1 黏性土壩漫頂潰決現(xiàn)場原型試驗(yàn)(壩高9.7 m)Fig.1 Prototype test of clayey soil dam failure due to overtopping(dam height: 9.7 m)

圖2 潰壩跌坎水流結(jié)構(gòu)示意Fig.2 Schematic diagram of headcut flow structure of dam-break

圖3 數(shù)值模擬計(jì)算區(qū)域及網(wǎng)格劃分Fig.3 Numerical simulation region and mesh generation

為方便利用前期開展的物理模型試驗(yàn)[13]進(jìn)行驗(yàn)證，計(jì)算模型中取壩體高度0.2 m；上、下游原始坡度1∶2；在壩體下游坡面中上部設(shè)不同高度的垂直跌坎；壩前水庫模擬長度取0.5 m，以保證水位充分穩(wěn)定；壩后0.44 m處設(shè)置為出口，減小出口段對壩體水流的影響。模型計(jì)算區(qū)域及網(wǎng)格劃分如圖3所示，為了最大限度保證計(jì)算精度，整個(gè)計(jì)算域采用均一結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，網(wǎng)格長寬比1∶1，只在跌坎處等重點(diǎn)部位進(jìn)行局部加密，最大網(wǎng)格尺寸0.002 5 m，總體網(wǎng)格數(shù)量為(120～170) 萬。

數(shù)模計(jì)算時(shí)，上游入口采用靜水壓力邊界條件；下游出口則采用自由出流邊界條件，對出口水流不做過多限制；模型底部采用不可滑移壁面邊界條件；在沿壩軸線方向兩側(cè)設(shè)置為對稱邊界條件，即邊界上沒有不穩(wěn)定的特性和剪切存在。

1.5 模型驗(yàn)證

潰壩水流流經(jīng)跌坎，由于體型突變，形成沖射水流對下游壩面造成沖刷。前期通過物理模型對潰壩典型水流進(jìn)行了初步試驗(yàn)研究，將數(shù)模計(jì)算結(jié)果與模型試驗(yàn)結(jié)果進(jìn)行對比分析，如圖4～5所示。從圖4～5可見，數(shù)模計(jì)算結(jié)果(壩頂、坎后壩坡水面線以及壩頂水流流速)與物模成果基本一致，僅局部稍有差異，鑒于潰壩跌坎水流紊動強(qiáng)烈、結(jié)構(gòu)復(fù)雜，可認(rèn)為所建數(shù)學(xué)模型對潰壩跌坎水流進(jìn)行了準(zhǔn)確模擬，具備了可靠的計(jì)算精度。

圖4 水面線對比驗(yàn)證Fig.4 Contrast and verification of the water surface line

圖5 壩頂平均流速對比驗(yàn)證Fig.5 Contrast and verification of the mean flow velocity on dam crest

編號S1S2S3S4S5S6S7S8跌坎高度/cm4666681012漫頂水深/cm32345333

2 計(jì)算結(jié)果分析

跌坎沖射水流是黏性土壩潰決發(fā)展過程中的代表性水流結(jié)構(gòu)，對其水力特性開展研究可從根本上了解掌握壩體形態(tài)的發(fā)展演變過程，具有重大的科研價(jià)值。

2.1 數(shù)模計(jì)算工況

本文針對不同跌坎高度和不同漫頂水深工況開展全面三維數(shù)值模擬，對漫頂?shù)菜鞯臎_刷主動力進(jìn)行細(xì)致研究，從而為深入了解掌握黏性土壩潰決演變過程奠定基礎(chǔ)，具體的數(shù)模計(jì)算工況如表1所示。

2.2 跌坎整體水流結(jié)構(gòu)及流態(tài)分析

圖6為跌坎水流整體流場結(jié)構(gòu)，圖7為漫頂沖射水流的流線分布情況。

圖6 跌坎水流整體流場結(jié)構(gòu)(S3工況)Fig.6 Flow structure of headcut flow (S3 condition)

圖7 跌坎沖射水流流線分布(S3工況)Fig.7 Streamline distribution of headcut flow (S3 condition)

從圖6可見：①壩體上游庫水位基本不變，在到達(dá)壩頂前緣時(shí)開始有明顯的跌落現(xiàn)象，且漫頂水深在壩頂繼續(xù)降低；②在沖射水流落點(diǎn)上游側(cè)并非充滿水，而是將形成一定深度的坎下回流區(qū)，這與試驗(yàn)結(jié)果完全一致；③在回流區(qū)內(nèi)存在漩渦水流結(jié)構(gòu)，使得水域中心位置流速稍小，且在沖射主流與落點(diǎn)上游側(cè)水域之間存在較為明顯的高速水流與低速水流的剪切層；④坎后下游坡面水流水面線較為光滑，水流流速沿程迅速增大，然而在壩址處由于壩體體型突變，存在局部流速減小區(qū)域。

由圖7可見：①水流流線在庫區(qū)分布較為均勻，接近壩體時(shí)，受壩體的阻擋作用出現(xiàn)流線彎轉(zhuǎn)，并在壩頂處匯集；②跌坎落點(diǎn)上游側(cè)環(huán)狀流線結(jié)構(gòu)表明該處存在明顯的漩渦水流結(jié)構(gòu)。

2.3 跌坎沖射水流水面線分布

圖8(a)為同一跌坎高度下不同漫頂水深工況的水面線分布情況，從中可以看出：①水流流經(jīng)跌坎形成沖射水流直接沖擊跌坎水平面，并在落點(diǎn)上游側(cè)形成具有一定深度的坎下回流區(qū)；②對于不同漫頂水深，沖射水平距離(沖射水流落點(diǎn)距跌坎底部的距離)和回流區(qū)水深都隨漫頂水深的增大而增加；③在壩頂和坎后坡面，壩面水深沿程逐漸下降，并均隨漫頂水深的增加而增大。圖8(b)為不同坎高工況下的壩面水面線分布情況，可見：①壩頂水面線分布基本保持一致，說明跌坎高度對于壩頂水流形式影響不大；②沖射水流空中軌跡具有較大的相似性，落點(diǎn)位置隨著坎高的增加而逐漸增大；③不同跌坎高度下，坎后壩面水深沿程變化情況較為相似。

圖8 不同工況跌坎沖射水流水面線分布Fig.8 Water surface line distribution under different conditions

2.4 壩面剪切應(yīng)力分布

圖9 相同坎高(6 cm)不同漫頂水深壩面剪切應(yīng)力分布Fig.9 Shear stress distribution under different overtopping depths

圖10 相同漫頂水深(3 cm)不同跌坎高度下壩面剪切應(yīng)力分布Fig.10 Shear stress distribution under different headcut heights

2.4.1相同坎高不同漫頂水深下壩面剪切應(yīng)力分布 漫頂潰決水流產(chǎn)生的壁面剪切應(yīng)力是壩體產(chǎn)生沖坑及侵蝕破壞的直接原因，圖9繪制了相同坎高(6 cm)不同漫頂水深工況下壩面剪切應(yīng)力分布情況。從中可以看出：①剪切應(yīng)力在壩頂上游壩肩位置存在略減小趨勢，之后沿程迅速增大，并且在壩頂下緣的壩肩處存在剪切應(yīng)力極大值；②在跌坎水平面水流落點(diǎn)兩側(cè)分別出現(xiàn)了剪切應(yīng)力極大值，剪切應(yīng)力分布呈現(xiàn)“雙峰”特征，相比而言下游側(cè)的剪切應(yīng)力峰值較大。在沖射水流落點(diǎn)上游側(cè)，首先出現(xiàn)剪切應(yīng)力極大值，隨后逐漸減小，至跌坎底部時(shí)相應(yīng)剪切應(yīng)力很?。虎墼诘菜矫媾c下游壩坡連接部位由于體型變化再次出現(xiàn)剪切應(yīng)力的極大值，而若水流落點(diǎn)距離連接部位較近，則兩個(gè)剪切應(yīng)力極值有可能“合并”成為較大極值，如漫頂水深4和5 cm工況所示；④坎后坡面剪切應(yīng)力分布較為規(guī)則，沿程逐漸增大，而在壩趾處受壩體體型變化出現(xiàn)剪切應(yīng)力極小值；⑤盡管在不同壩體區(qū)域存在不同的剪切應(yīng)力分布形式，然而隨著漫頂水深的增大，壩面剪切應(yīng)力整體呈現(xiàn)增大的趨勢，且在壩坡中下部增幅逐漸加大。

2.4.2相同漫頂水深不同坎高壩面剪切應(yīng)力分布 本文亦對不同跌坎高度工況進(jìn)行了數(shù)值模擬，圖10為相同漫頂水深(3 cm)不同跌坎高度工況下的壩面剪切應(yīng)力分布。由圖10可以看出：①各工況壩頂剪切應(yīng)力分布基本相同，在壩頂上游壩肩位置出現(xiàn)小范圍減小趨勢，之后剪切應(yīng)力沿程迅速增大，并在壩頂下緣位置出現(xiàn)極大值；②不同坎高將導(dǎo)致不同沖射水流落點(diǎn)位置，沖射落點(diǎn)水平距離(落點(diǎn)距跌坎根部的水平距離)隨跌坎高度的增大而增大；③沖射水流在落點(diǎn)兩側(cè)同樣形成剪切應(yīng)力極大值，呈現(xiàn)“雙峰”型式，隨著跌坎高度的增大，水流落點(diǎn)兩側(cè)的剪切應(yīng)力極值也逐漸增大。然而在跌坎根部附近各工況剪切應(yīng)力分布基本一致，剪切應(yīng)力值均較小；④對于同一跌坎高度，坎后坡面剪切應(yīng)力依然沿程增大，然而對于不同坎高工況，坎后剪切應(yīng)力隨著坎高的增大整體呈減小趨勢。

2.4.3不同工況下跌坎垂直壁面剪切應(yīng)力分布 圖11(a)為相同跌坎高度不同漫頂水深工況下跌坎垂直面剪切應(yīng)力分布，由圖可見：①在跌坎垂直面上，剪切應(yīng)力呈“單峰”分布，在垂直面中部剪切應(yīng)力較大，跌坎底部及水域上部剪切應(yīng)力較??；②不同漫頂水深工況下，隨著漫頂水深的增大，垂直壁面上剪切應(yīng)力最大值逐漸增大，最大剪切應(yīng)力作用點(diǎn)位置也逐漸上移；③隨著漫頂水深的不斷增大，盡管垂直壁面剪切應(yīng)力整體呈現(xiàn)增大的趨勢，然而峰值變幅存在逐漸減小的趨勢。圖11(b)為相同漫頂水深不同跌坎高度工況下，跌坎垂直壁面剪切應(yīng)力的分布，由圖可見：①壩體跌坎垂直面上剪切應(yīng)力均呈“單峰”分布，自跌坎底部沿垂直方向剪切應(yīng)力先增大后減?。虎谠诘哺叨戎饾u增大的過程中，垂直壁面剪切應(yīng)力峰值存在先增大后減小的趨勢，這一規(guī)律表明：在跌坎逐漸加深的過程中存在一定的限值，當(dāng)坎高達(dá)到這一高度后，水流對于跌坎垂直面的作用力反而會越來越小。

跌坎垂直面上的剪切應(yīng)力分布直接關(guān)系到壩體縱向沖蝕發(fā)展過程，對比圖9，10和11可知，跌坎沖射水流在水平壁面上產(chǎn)生的最大水流剪切應(yīng)力比垂直面上剪切應(yīng)力要大很多，有的工況甚至可達(dá)幾十倍，由此可初步推斷黏性土壩漫頂潰決過程中各級跌坎主要以“臺階水平面刷深下切”的方式進(jìn)行合并。

圖11 不同工況下跌坎垂直壁面剪切應(yīng)力分布Fig.11 Shear stress distribution on vertical wall surface

2.5 跌坎底壁水平流速分布

跌坎水平面底壁流速分布直接反映了坎下回流區(qū)的水流結(jié)構(gòu)特性，圖12(a)為相同跌坎高度不同漫頂水深下跌坎水平面底壁附近流速分布，可見：①水流水平流速在落點(diǎn)處分別向兩側(cè)分流，并迅速增大至相應(yīng)最大流速?？傮w上看，反向水流流速小于主流方向水流流速；②反向水流流速迅速增大至最大流速，之后逐漸減小，至跌坎底部時(shí)，水流流速很小，整體上呈“單峰”型式。圖12(b)為相同漫頂水深不同跌坎高度下跌坎水平面底壁流速分布，從中可以看出：①沖射水流依然自水流落點(diǎn)處進(jìn)行分界，形成流向截然相反的兩股水流；②跌坎高度越大，沖射水流落點(diǎn)距跌坎底部距離越遠(yuǎn)，水流流速越大；③盡管不同工況下水流流速分布存在一定差別，然而在跌坎底部流速均較小，甚至接近于0，這與剪切應(yīng)力分布規(guī)律相類似，表明潰決水流對跌坎根部的直接沖刷作用相對較弱。

圖12 不同工況跌坎水平面底壁流速分布Fig.12 Velocity distribution of horizontal plane near bottom

3 結(jié) 語

漫頂水流作為壩體沖蝕的主導(dǎo)作用力，對于土壩潰決發(fā)生、發(fā)展過程至關(guān)重要，本文采用RNGk-ε紊流模型和VOF技術(shù)針對黏性土壩漫頂潰決過程中的代表性水流結(jié)構(gòu)——潰決跌坎水流開展了三維數(shù)值模擬研究，揭示了不同工況下壩體跌坎上的剪切應(yīng)力、流速分布規(guī)律，為進(jìn)一步發(fā)展完善黏性土壩漫頂潰決“跌坎式”溯源沖蝕機(jī)理提供了理論基礎(chǔ)。研究結(jié)果表明：

(1)跌坎沖射水流在落點(diǎn)上游側(cè)并非充滿水，而是形成一定深度的坎下回流區(qū)，回流區(qū)內(nèi)存在漩渦紊動水流結(jié)構(gòu)；相同坎高工況沖射水平距離和回流區(qū)深度均隨漫頂水深的增大而增加，不同坎高工況沖射水流結(jié)構(gòu)基本一致，具有較大相似性。

(2)壩面剪切應(yīng)力在壩頂上游壩肩位置存在略減小趨勢，之后沿程迅速增大，并且在壩頂下緣的壩肩處存在剪切應(yīng)力極大值；在跌坎水平面水流落點(diǎn)兩側(cè)會出現(xiàn)“雙峰”分布特征，且下游側(cè)峰值較大；在跌坎水平面與下游壩坡連接部位會出現(xiàn)剪切應(yīng)力的極大值，而坎后坡面剪切應(yīng)力沿程逐漸增大，在壩趾處則受體型變化影響出現(xiàn)極小值。

(3)跌坎垂直面上，剪切應(yīng)力呈“單峰”分布，在水深中部剪切應(yīng)力較大，底部及水域上部剪切應(yīng)力較??；隨漫頂水深的增大，剪切應(yīng)力峰值逐漸增大，最大剪切應(yīng)力作用點(diǎn)位置也逐漸上移；隨跌坎高度逐漸增大，剪切應(yīng)力峰值存在先增大后減小的趨勢。

(4)跌坎水平面底壁流速在落點(diǎn)處向兩側(cè)分流，并迅速增大至相應(yīng)最大流速。總體上看，反向水流流速呈“單峰”型式，小于主流方向水流流速，至跌坎底部時(shí)，水流流速很??；跌坎高度越大，沖射水流落點(diǎn)距跌坎底部距離越遠(yuǎn)，水平面底壁附近水流最大流速也越大。

(5)由數(shù)模結(jié)果分析可知，跌坎水平面上最大水流剪切應(yīng)力比垂直面上剪切應(yīng)力要大很多，甚至可達(dá)幾十倍，由此可初步推斷黏性土壩漫頂潰決過程中各級跌坎主要以“臺階水平面刷深下切”的方式進(jìn)行合并。

[1] 中華人民共和國水利部， 中華人民共和國國家統(tǒng)計(jì)局. 第一次全國水利普查公報(bào)[M]. 北京: 中國水利水電出版社， 2013. (The Ministry of Water Resources of the People’s Republic of China, National Bureau of Statistics of the People’s Republic of China. Bulletin of first national census for water[M]. Beijing: China Water Power Press, 2013. (in Chinese))

[2] 李云， 李君. 潰壩模型試驗(yàn)研究綜述[J]. 水科學(xué)進(jìn)展, 2009, 20(2): 304- 310. (LI Yun, LI Jun. Review of experimental study on dam-break[J]. Advances in Water Science, 2009, 20(2): 304- 310. (in Chinese))

[3] 朱勇輝， VISSER P J， VRIJLING J K， 等. 堤壩潰決試驗(yàn)研究[J]. 中國科學(xué): 技術(shù)科學(xué), 2011, 41(2): 150- 157. (ZHU Yonghui, VISSER P J, VRIJLING J K, et al. Experimental investigation on breaching of embankments[J]. Scientia Sinica (Technological), 2011, 41(2): 150- 157. (in Chinese))

[4] 郭軍. 歐美國家近期潰壩研究及發(fā)展動向[J]. 中國水利, 2005(4): 23- 26, 29. (GUO Jun. Recent studies and activities on dam breach in European countries and United States[J]. China Water Resources, 2005(4): 23- 26, 29. (in Chinese))

[5] HANSON G J, COOK K R, HUNT S L. Physical modeling of overtopping erosion and breach formation of cohesive embankments[J]. Transactions of the ASAE， 2005， 48(5)： 1783- 1794.

[6] ZHANG Jianyun, LI Yun, XUAN Guoxiang, et al. Overtopping breaching of cohesive homogeneous earth dam with different cohesive strength[J]. Science in China Series E: Technological Sciences, 2009, 52(10): 3024- 3029.

[7] POWLEDGE G R, RALSTON D C, MILLER P, et al. Mechanics of overflow erosion on embankments[J]. Journal of Hydraulic Engineering， 1989， 115(8)： 1056- 1075.

[8] 李云, 王曉剛, 劉火箭, 等. 土石壩漫頂過程水力特性分析[J]. 水動力學(xué)研究與進(jìn)展(A輯), 2012, 27(2): 147- 153. (LI Yun, WANG Xiaogang, LIU Huojian, et al. Overtopping hydraulic characteristics of earth-rock dam[J]. Chinese Journal of Hydrodynamics, 2012, 27(2)： 147- 153. (in Chinese))

[9] 魏文禮, 郭永濤, 王紀(jì)森.一維潰壩洪水波的高精度數(shù)值模擬[J]. 計(jì)算力學(xué)學(xué)報(bào), 2007, 24(3): 362- 364. (WEI Wenli, GUO Yongtao, WANG Jisen. One-dimensional dam-break numerical simulation based on ENO scheme[J]. Chinese Journal of Computational Mechanics, 2007, 24(3): 362- 364. (in Chinese))

[10] 黃金池， 何曉燕. 潰壩洪水的統(tǒng)一二維數(shù)學(xué)模型[J]. 水利學(xué)報(bào)， 2006， 37(2)： 222- 226. (HUANG Jinchi, HE Xiaoyan. Unified 2-D numerical model for simulating dam break wave propagation[J]. Journal of Hydraulic Engineering, 2006， 37(2)： 222- 226. (in Chinese))

[11] 張新華， 隆文非， 謝和平，等. 任意多邊形網(wǎng)格2D FVM模型及其在城市洪水淹沒中的應(yīng)用[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版)， 2007， 39(4)： 6- 11. (ZHANG Xinhua, LONG Wenfei, XIE Heping, et al. 2D FVM model based on an arbitrary polygonal mesh system for the simulation of flood inundation and its application in municipal regions[J]. Journal of Sichuan University(Engineering Science Edition), 2007, 39(4)： 6- 11. (in Chinese))

[12] 夏軍強(qiáng), 王光謙, LIN Binliang, 等. 復(fù)雜邊界及實(shí)際地形上潰壩洪水流動過程模擬[J]. 水科學(xué)進(jìn)展, 2010, 21(3): 289- 298. (XIA Junqiang, WANG Guangqian, LIN Binliang, et al. Two-dimensional modelling of dam-break floods over actual terrain with complex geometries using a finite volume method[J]. Advances in Water Science, 2010, 21(3)： 289- 298. (in Chinese))

[13] 祝龍, 李云, 宣國祥, 等. 潰壩典型水流水力特性試驗(yàn)研究[J]. 水資源與水工程學(xué)報(bào), 2015, 26(5): 132- 136. (ZHU Long, LI Yun, XUAN Guoxiang, et al. Experiment on hydraulic characteristics of typical dam break flow[J]. Journal of Water Resources and Water Engineering, 2015, 26(5): 132- 136. (in Chinese))

Numericalstudyofhydrodynamiccharacteristicsofdam-breakheadcutflow

ZHU Long， ZHOU Donghui， LI Yun， XUAN Guoxiang， WANG Xiaogang

(StateKeyLaboratoryofHydrology-WaterResourcesandHydraulicEngineering,NanjingHydraulicResearchInstitute,Nanjing210029,China)

The dam-break process of clay soil dam under overtopping condition, related to multi-disciplines, is very complex. Although a large number of domestic and foreign physical model test results show that the headcut erosion caused by overtopping is the main breach process of most clay soil dams, the micro mechanism between water flow and dam body in the form of erosion is still not clear. The water flow plays a leading role in the development of dam break as it’s the active force during overtopping failure. In the study, 3D numerical simulation research is carried out on the representative flow structure—the headcut flow, by using the RNGk-εturbulence model and the volume of fluid (VOF) method. The flow structure, flow pattern and hydraulic characteristics of the dam headcut are analyzed in detail. Also, the distributions of shear stress and flow velocity in different working conditions are revealed. By using the research results, the analysis of the forces on dam body can be carried out from the point of view of hydrodynamics. The main way of all-level headcuts merged is inferred initially. The research results can help to make a clear understanding of the dam-break process.

clay soil dam; overtopping; headcut flow; hydrodynamic characteristics; numerical simulation

TV871

A

1009-640X(2017)05-0001-08

10.16198/j.cnki.1009-640X.2017.05.001

祝龍， 周冬卉， 李云， 等. 土壩潰決跌坎水流水動力特性數(shù)值模擬[J]. 水利水運(yùn)工程學(xué)報(bào), 2017(5): 1-8. (ZHU Long, ZHOU Donghui, LI Yun, et al. Numerical study of hydrodynamic characteristics of dam-break headcut flow[J]. Hydro-Science and Engineering, 2017(5): 1-8. (in Chinese))

2016-11-21

國家自然科學(xué)基金委青年科學(xué)基金項(xiàng)目(51409173)；中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金(Y116009，Y114004)

祝 龍(1988—)， 男， 河北邢臺人， 工程師， 博士， 主要從事潰壩水力學(xué)、通航水力學(xué)研究。

E-mail： zhulong0101@163.com