范家驊，祁 偉，戴 清

（中國水利水電科學研究院，北京 100048）

1 異重流潛入及其工程問題

異重流的形成是兩種不同密度的流體相遇時，在一定條件下較重或較輕的流體從有壓水流或無壓水流向另一種不同密度的流體開始過渡到分層流（異重流）的相對運動過程。需要研究的問題是：從非分層流過渡到分層流需要遵循哪些條件。這種過渡現(xiàn)象即在潛入點下游形成底部異重流流動，或稱之為異重流渾水楔。

（1）水庫入口處。洪峰期間或非汛期挾沙水流進入水庫末端回水區(qū)，流至一定的水深處潛入庫底形成底部異重流，即異重流渾水楔。潛入點下游渾水楔長度隨洪峰歷時而變，洪峰歷時短時，渾水楔流不到壩址，全部淤在水庫內(nèi)；洪峰歷時長時，可通過壩底孔排出部分泥沙；而潛入點上游壅水區(qū)會形成三角洲泥沙淤積并導致河段水位抬高［1-2］。

（2）盲腸河段與河道的交匯處（船閘和引航道）。河道挾沙水流與盲腸河段內(nèi)靜止水體相遇因密度不同產(chǎn)生相對運動（即所謂交換水流），挾沙水流潛入底部并流進盲腸河段，同時上層清水則自盲腸河段流出進入河道。由于河道含沙水流這樣持續(xù)不斷地進入盲腸段而造成盲腸河段的累積淤積。

（3）河道交匯處。河道主流含沙量大于支流含沙量時，主流渾水潛入支流形成底部異重流并向支流上游方向運動，且上溯至一定距離，形成一段具有一定長度的渾水楔；而當支流含沙量大于主流含沙量時，支流渾水則潛入主流，也形成底部異重流和具有一定長度的渾水楔。工程設計或研究的內(nèi)容主要為楔內(nèi)泥沙沿程淤積和淤積處理方法。

（4）河口與海域的交匯處。海區(qū)鹽水與河道水流在河口處相遇，在一定條件下會形成向上游的鹽水楔運動；而鹽水中含沙量較大時（受海灘風浪掀沙影響），上溯水流與河道水流相遇，在一定條件下會形成含沙分層運動的渾水楔，渾水楔沿程泥沙淤積會抬高航道底高，需挖除淤積礙航泥沙，才能保持一定航深的航道通航。另外河道挾沙水流流出河口，因其含沙量較小，渾水密度小于海水密度，渾水會向海面擴散，形成上層異重流。

2 潛入現(xiàn)象研究回顧

異重流潛入現(xiàn)象的水槽實驗、水力理論分析和模型數(shù)值計算工作，國內(nèi)國外不少學者進行過研究，以下將分述和回顧各家研究成果；同時本文還將著重于渾水異重流在等寬水槽內(nèi)的現(xiàn)象，在筆者先前工作的基礎上做進一步的分析，與前人工作進行對比，以獲得對異重流潛入現(xiàn)象進一步的認識，并對潛入現(xiàn)象的水力分析中若干問題進行討論，為工程設計中水力計算提供參考［1-2］。

2.1 水庫壅水區(qū)渾水潛入點實驗 渾水異重流潛入點水槽實驗有范家驊［3］、蘆田和男［4］、曹如軒［5］、曹如軒與錢善琪［6］、俞維升［7］、姚鵬與王興奎［8-9］和焦恩澤［10］等。筆者分別在寬15 cm水槽和寬50 cm水槽內(nèi)進行過渾水潛入實驗，異重流潛入過程如圖1所示。

水槽實驗中一般觀測記錄潛入點的水深、含沙量、流速以及潛入點潛入后交界面線的沿程變化，并進行部分測次若干斷面的含沙量垂線分布和流速分布的測量。對潛入點測定數(shù)據(jù)進行量綱分析，可獲得潛入點斷面密度Froude數(shù)Fp的值。根據(jù)Schijf和Schonfeld［11］的異重流雙層漸變流方程，潛入點水深過渡到潛入點異重流水深之間會存在一個交界面線的拐點，該處密度Froude數(shù)等于1；而潛入點水深大于該處臨界密度Froude數(shù)中的水深，故可見潛入點處密度Froude數(shù)Fp小于1。資料分析得潛入點處密度Froude數(shù)：

圖1 異重流潛入圖

式中：up為異重流潛入點流速；hp為潛入點水深；g′=(?ρ/ρ)g。1961年筆者又在長50 m、寬0.5 m、高2 m水槽內(nèi)進行異重流潛入點實驗，進沙粒徑d90=0.0145～0.05mm，d50=0.0055～0.008 mm，觀測沿程渾水水深、含沙量和泥沙粒徑的變化，4次實驗后得Fp=0.58～0.66。

蘆田和男［4］進行渾水實驗，分析得：其中hp為潛入點水深；S為底部比降；q為單寬流量。曹如軒等曾對上式利用他們的實驗資料進行比較，其系數(shù)大于蘆田和男的0.365，一般在0.4～0.57之間，平均值為0.44。

曹如軒［5］利用3個水槽進行含沙量6.5～715 kg/m3的渾水潛入實驗，其中根據(jù)含沙量在30 kg/m3以下的資料得：而渾水含沙量大于100 kg/m3時Fp值較小，且隨著含沙量的增加Fp值減小。含沙量在100～360 kg/m3時，F(xiàn)p=0.4～0.2。曹如軒等對含沙量高且流體黏滯系數(shù)大的高含沙水流，利用Bingham體有關參數(shù)和流體阻力系數(shù)聯(lián)系來分析潛入點密度Froude數(shù)與含沙濃度的某種關系，分析得到符合實驗數(shù)據(jù)變化趨勢的結果。曹如軒與錢善琪［6］進行含粗沙高含沙量潛入點水槽實驗，得到與低含沙量潛入類似的Fp值結果。俞維升［7］用高岺土和石英沙挾沙水流以及鹽水進行潛入點水槽實驗，觀察到潛入點向下游移動至一定處趨于穩(wěn)定的現(xiàn)象，且平均Fp≈0.71；同時通過沿程流速測量與分析，計算得潛入后異重流流量沿程增加的現(xiàn)象。姚鵬與王興奎［8-9］在長而寬的水槽內(nèi)進行4種底坡的渾水異重流潛入點試驗，槽寬1.2 m、長63.8 m，觀察到潛入點Fp值因底坡加大而降低，變化范圍在0.67～0.57之間。焦恩澤［10］進行不同含沙量（12～479 kg/m3）的潛入點水槽實驗，其中大量試驗為高含沙量異重流實驗。同曹如軒實驗類似，F(xiàn)p值隨含沙量增加而變小，并利用Bingham體參數(shù)與潛入點建立一定關系。范家驊［12］考慮小浪底水庫的含沙量范圍可用曹如軒、焦恩澤等人的高含沙異重流潛入點數(shù)據(jù)，并采用量綱分析方法，得到含沙量在400 kg/m3以下范圍內(nèi)的Fp平均關系為Fp=0.9c0.1，c為含沙量。

近年來國內(nèi)關于水庫潛入點Fp的研究公開發(fā)表的論文近20余篇，主要為小浪底水庫高含沙量異重流的控制運用中進行異重流潛入點參數(shù)的研究，并對包括小浪底水庫及其它水庫的觀測資料、水槽實驗和模型試驗資料進行有關水沙因素與Fp的影響分析。如李濤和張俊華等［13］、李書霞等［14］、李濤和夏軍強等［15］的研究，可見在高含沙量異重流潛入點處，F(xiàn)p參數(shù)與含沙量值有一定影響。

各家實驗結果以及渾水模型、數(shù)值計算結果列于表1。

表1 渾水異重流潛入點水槽實驗、數(shù)值計算與水力分析

2.2 鹽水與冷熱水潛入點實驗 除了渾水實驗，尚有利用鹽水和冷熱水模擬渾水進入水庫壅水區(qū)潛入的水槽實驗，有Keulegan［17］、Singh與Shah［18］、Itakura與Kishi［19］、Fukuoka等［20］、Kan與Tamai［21］等人的鹽水實驗，以及Farrell與Stefan［22］、Akiyama與Stefan［23-24］等人的冷熱水實驗。各家實驗結果列于表2。

Keulegan［17］進行無潮河口鹽水水槽實驗，觀測進槽的鹽水異重流前鋒加速運動規(guī)律。令進口水深即為潛入點水深，以hp表示；異重流前鋒速度以ud表示。用資料分析可得：hp=2hd且0.57，故有

表2 鹽水和冷熱水潛入點水槽實驗

Singh與Shah［18］用鹽水在水槽內(nèi)進行了異重流潛入實驗。水槽長14 m，寬0.4 m，進槽單寬流量0.5～135 cm2/s，密度差為0.0005～0.013 g/cm3，異重流水深為1.5～17 cm，潛入點水深3～22.5 cm，水槽底坡S為0.005～0.02。他們使用量綱分析潛入點水深Re（Reynolds數(shù)）和底坡S諸因素的函數(shù)，并點繪潛入點水深與Re和S底坡的關系，表明在Re=600～11000以及S=0.0056～0.0215之間，此兩因素對潛入點水深無多大影響；在低流量和高密度的某些測次，水面有明顯的分界線，但未觀測到有向下游流動的異重流，這些測次的值小于2 cm，均未包括在分析之中。最后分析對較大值有：

以密度Froude數(shù)Fp可表示為：

Fukuoka等［20］在長6.9 m、高0.9 m、寬0.20 m水槽內(nèi)進行鹽水潛入點實驗。槽底比降為1/10的實驗7次、比降1/60的實驗7次，觀測潛入點交界面形狀、流速分布、鹽分分布等；并分析點繪潛入點水深與的關系；還加入石橋和巖崎在水庫中數(shù)次實驗數(shù)據(jù)，得平均值：Fu?kuoka等認為他們的實驗結果與他們在文中介紹Benjamin的理論分析結果接近，并計算出摻混系數(shù)E在0.046×10-2～1.8×10-2之間，還點繪實測和計算的沿程交界面高度的變化。

Kan與Tamai［21］在槽長6 m、高0.9 m、寬0.3 m、底坡為0.1～0.25的水槽內(nèi)進行鹽水異重流潛入實驗，進槽鹽水密度為1.00006～1.056 g/cm3，q=2.4～2.5 cm3/cm，得潛入點Fp平均值0.63，此值系26次數(shù)據(jù)0.45～0.92的平均值。

Farrell與Stefan［22］應用k-ε模型計算了潛入點水流現(xiàn)象，并進行了水槽冷熱水實驗。槽長40呎、寬16.5呎、底坡S=0.047，觀測潛入點處水深、流速和溫度差以及潛入?yún)^(qū)下游異重流水深、流速及溫度等。實驗共進行7次，實驗資料分析得：Fp=0.69。實驗還觀測了1次摻混系數(shù)：式中Qm為異重流流量，Q0為進槽流量。

Akiyama與 Stefan［23-24］在槽寬61 cm、深30.5 cm且具有一邊擴展角1°、3°和7°的水槽內(nèi)進行冷熱水異重流潛入點實驗，研究槽寬擴寬對潛入的影響，實測獲得Fp=0.56～0.89之間，平均值為0.68；實驗中還測得潛入點下游異重流水深hd和潛入點水深hp的比值在0.65～0.9之間；并實測摻混系數(shù)γ＝0～0.31，摻混系數(shù)的定義為Q0為進槽流量，Qd、Qp分別為異重流流量和潛入點流量。該文中點繪Fp和F0的關系，可看出在1°時，F(xiàn)p在0.6～0.7之間；3°時Fp=0.65～0.8之間；7°時Fp=0.6～0.8之間；總的范圍在0.6～0.8之間，可見F0對Fp的影響不大。關于此關系，Stefan與Johnson［25］在他們論文中補充了擴展角3°共計5個新的測點，點繪在一起，顯示Fp與F0和δ兩者無多大關系。

2.3 “交換水流”渾水及鹽水潛入實驗與水力分析 引航道口門渾水潛入形成渾水楔、船閘鹽水交換以及無潮河口入海處含鹽水流潛入河口段形成鹽水楔示意圖見圖2。

圖2 引航道河口門異重流潛入圖

這類船閘“交換水流”潛入實驗的研究重點在于：求取船閘閘門打開后鹽水進入儲有清水的船閘時潛入形成底層水流的鹽水和前鋒速度兩者的關系，可用密度Froude數(shù)的值表示，式中符號：前鋒速度uf（即異重流流速ud），河口水深H（即潛入點水深hp），清水與鹽水的密度差?ρ/ρ。

進行相關理論分析和鹽水實驗工作的有O'Brien 與 Cherno［26］、 Yih C S（易 家 訓）［27］、Keulegan［28］、Middleton［29］、Rozovsky［30］、Kersey和 Hsu［31］、 Lam［32］、Barr［33］、Rottman 與 Simp?son［34］。利用渾水進行引航道異重流潛入實驗的有范家驊［1］。此外陳國謙和李行偉［35］利用k－ε（BN6）模型進行了數(shù)值計算。實驗結果與水力分析列于表3。

至于各家水槽實驗，由于上下層水深不同以及交界的阻力系數(shù)的影響，F(xiàn)p值略大于或略小于0.25。

水庫壅水區(qū)和引航道口門這兩類異重流潛入，因為潛入處的幾何和水流條件不同，主要由于上層水流流速ua值在水庫中很小，分析時可忽略不計；而在引航道中ua接近ud。故兩者的潛入點密度Froude數(shù)Fp有差異。綜上所述：水庫中Fp在0.5～0.8之間，引航道中Fp在0.25左右。

表3 船閘交換水流鹽水、渾水實驗與水力分析

3 異重流潛入現(xiàn)象的理論分析成果

3.1 von Karman的理論分析成果 von Karman［36］應軍事部門的求詢，討論一種理想流體模型。圖3為較重液體ρ2在較輕液體ρ1中運動，考慮較輕液體在無限深度假設前提下，推導下層較重液體恒定運動的流動速度c1。

von Karman應用A點和B點的Bernoulli方程：

在以一定速度運動的框架之中水流為恒定，因為相對于這框架，水流為靜止：

因此有：

3.2 Benjamin的理論分析成果 Benjamin［37］認為von Karman利用能量方程并忽略滯流點前端頭部高密度水流受到與上層水流混合的能量損失是不能接受的，雖然Benjamin的推理得到同樣的c1結果。易家訓［38］的看法是von Karman雖然在這一點上是一個疏忽，但使用Bernoulli方程所取上下游兩點，如果在上游的點取的位置更上游一些，下游的點更下游一些，使用Bernoulli方程則無問題。

Benjamin分析兩平板之間的空穴流的空穴移動速度，分析中忽略黏性和表面張力，上游無窮遠處水深為d，流速為c1，下游無限遠處水深為h，流速為c2。概化圖見圖4。

沿表面的Bernoulli方程，O點表面的壓力為零：

無限遠處的壓力p0：

上下游遠處的斷面壓力分布為靜水壓力分布，動量守恒方程為：

圖3 von Karman異重流概化圖

圖4 Benjamin空穴流概化圖

連續(xù)條件為c1d=c2h，解得：

其次Benjamin考慮水流的能量損失：假定水流有水頭損失，流速c2在很遠的下游會變成均勻流，定義Δ為水頭損失，則有：

此式與力平衡所得的下式：

相等，可得：當h＜0.5d時，Δ為負值，故在實際中不能發(fā)生。利用以上各式，可得當h＞0.5d時，有：

3.3 Singh和Shah的理論分析成果 Singh和Shah［18］除了進行鹽水潛入點水槽實驗外，還用動量方程分析潛入現(xiàn)象，動量方程為：

式中：左邊代表動量的改變率，β1和β2為動量系數(shù)；v1和v2為兩斷面的平均流速；P1和P2為兩斷面的壓力；F1和F2為底部和交界面的剪力阻力；Ws為兩斷面間液體重量的分力。

Singh在分析中作出如下假定：水壓力為靜水壓力分布；每一點的水流方向平行于槽底；潛入長度L和底坡S均為小值，故第2斷面總水深等于第1斷面的水深，即潛入點水深；交界面上的阻力F2假定平行于槽底。概化圖見圖5。

為簡化分析，忽略F1、F2和Ws小值，可得：

對上式的分析采用三種簡化：第一種是忽略潛入點斷面與下游異重流斷面之間底部阻力和交界面阻力，并忽略動量流速分布系數(shù)，得：第二種考慮動量流速分布系數(shù)，計算得該系數(shù)β1和β2等于1.18，則：第三種考慮潛入點潛入?yún)^(qū)距離之內(nèi)的底部和交界面阻力的影響，Singh從均勻異重流中計算得平均總阻力系數(shù)f為底部阻力系數(shù)與交界面系數(shù)之和，其值為0.13，最后得：

Singh按三種關系式計算的潛入點水深值同實驗測得的潛入點水深值點繪在同一圖上對比，顯示三條平均線在水深較小部分計算值偏大，在水深較大時則偏差較小。

3.4 Savage和Brimberg的理論分析成果 Savage與Brimberg［39］采用類似Benjamin的方法分析潛入點的水力特性。應用Bernoulli原理，分析恒定二維兩種密度ρ1和ρ2的流體，斷面1處于潛入點上游遠處，流速為u1，水深為h1；斷面2處于潛入點的下游遠處，流速為u2，水深為h2。概化圖見圖6。設

O點為滯流點，沿自由水面的O點的上游應用Bernoulli原理；其次應用Bernoulli原理自O點沿交界面下游；并設在潛入點上游相當遠處和下游相當遠處水壓力為靜水壓力，運用動量平衡方程可導得：

圖5 Singh潛入點概化圖

圖6 Savage潛入點概化圖

令潛入點水深h1即為hp，可表示為：

Savage用以上推導的關系式同Singh的潛入點實驗水深比較，兩者較接近，但測點分散。實驗數(shù)據(jù)顯示除以外尚有其它參數(shù)的影響，如底坡和交界面剪應力。Savage考慮到第一部分的水力分析工作忽略了底部比降和剪應力，故利用Schijf與Schonfeld［11］的雙層漸變流一維運動方程來分析潛入點現(xiàn)象，分析中令水流為恒定及上層流速為零，對動力方程進行數(shù)值計算。

Savage為了定性描述潛入點下游的交界面形狀，利用底層漸變流方程，將潛入點處的底部作為坐標零點，垂向坐標z代表水深，縱向坐標為距離x，推求漸變流方程的解。這樣就類似于明渠漸變水流，可分為數(shù)種水面線類型。

從潛入點起向下游的交界面線在緩坡和急坡情況下可能出現(xiàn)的是接近M2和S2的曲線。數(shù)值計算的結果可建立Fp為f，α和S（底坡）的函數(shù)。對于緩坡有：

上式是在f=0.01～0.09，α=0.2～0.8，α=fi/f以及緩坡范圍內(nèi)的計算結果。

而在急流條件下，計算表明其性質(zhì)與緩坡情況類似。但有一點不同，當潛入點密度Froude數(shù)為Fp=0.8至0.845時，交界面曲線出現(xiàn)劇烈變化。Savage分析Singh的實驗資料，其潛入點Fp值小于0.8，在0.3～0.8之間，點繪的關系線，并以參數(shù)Fp=0.3～0.8作直線，測點絕大部分在Fp＝0.3～0.8線的范圍之內(nèi)。

3.5 Jain的理論分析成果 Jain［40］著文檢驗Savage的論文，提出一些不同意見。他認為用能量方程于潛入點附近的水流在物理上是不可能的。潛入點附近的能量耗損不可忽視，因此他利用動量方程和能量方程（考慮能量損失Δ）進行分析，概化圖見圖7。并得到以下關系式：

圖7 Jain潛入點概化圖

Jain指出Savage和Brimberg的分析，F(xiàn)1＝0.5為以上關系式當HL＝0和且β＝0.5時的特解。Jain對Δ值作出的假定，表示β與的相應關系。Jain計算出在不同β值時動量方程的解，并且指出對于某一已知β值存在一個F1的極大值（F1max），大于此值時不能發(fā)生潛入現(xiàn)象。

Jain的論文第二部分分析交界面曲線的計算，認為積分的方向是敏感的，計算應從控制斷面算起。他利用Schijf與Schonfeld［11］的雙層漸變流的一維方程進行計算，并假定其上層水流的平均流速為零。上層水面為水平，水深為a1，下層水深為a2，在潛入點處a1=0，如圖5所示圖形?？傻茫?/p>

上式僅包含兩個無量綱參數(shù)λ和ηc，可用不同的λ和ηc值和不同下游條件求得數(shù)值解，其積分是從下游控制斷面向上游方向計算。

對緩坡的計算，設一定的阻力系數(shù)α=0.5、λ=1 3和ηc＜1（即臨界水深小于正常水深，如ηc=0.5）以及若干下游條件向上游進行積分。所求得的對于大范圍的下游條件的交界面曲線在潛入點附近連接，且潛入點水深值與下游邊界條件無關。但當積分是從上游向下游進行計算，即使很小的差別，卻會造成不同的交界面曲線。

對陡坡（ηc＞1）計算兩種典型的交界面曲線ηc=1.3以及λ=1 3（α=0.5），交界面曲線S1和S3均為緩流。

對于潛入水深，無量綱潛入點水深ηp僅為ηc和λ（或α）的函數(shù)，有下列近似關系：

3.6 朱鵬程的理論分析成果 朱鵬程［41］首先對異重流潛入條件按以下概化圖形（見圖8）進行了理論分析，得：

得：

曹如軒等［5］用同樣方法，得上述相同公式。朱鵬程注意到上列動量方程導得的方程類似水躍前后的共軛水深，因此其間須經(jīng)過臨界水深。他認為要產(chǎn)生共軛水深，F(xiàn)p必須小于1且Fd必須大于1。

蘆田和男［4］對異重流潛入現(xiàn)象進行分析和渾水實驗，求解潛入點水深。解得：

3.7 Akiyama和Stefan的理論分析成果 Akiyama與Stefan［42-43］分析潛入現(xiàn)象時，首先利用動量公式并引進摻混系數(shù)γ，推導潛入點密度Froude數(shù)與潛入點下游異重流厚度等因子的關系。

圖9概化圖形中的CVⅠ動量方程假定可忽略混合區(qū)的底部阻力以及忽略比降的影響，得下式：

他們定義摻混系數(shù)γ（水庫清水混入底部異重流的水量）：

CVⅡ的動量方程：

圖8 朱鵬程異重流潛入點概化圖

圖9 Akiyama異重流潛入點概化圖

以上兩動量方程合并，得：

式（19）中Fd為潛入點下游的異重流密度Froude數(shù)，此式表明Fd和和γ的關系。

為了分析潛入點水沙因子與底部阻力、交界面阻力以及下游異重流的影響，Akiyama作出假定：潛入點下游異重流在急坡時Fd等于臨界值；在緩坡時Fd等于正常值Fn。他們利用Ellison與Turner［44］的具有摻混系數(shù)的異重流緩變方程，求解緩坡和急坡時的與摻混系數(shù)、阻力系數(shù)、底部比降以及異重流流速分布系數(shù)的關系式。由于所求得的關系式中阻力系數(shù)與摻混系數(shù)的各項關系式過長，為了同Singh等人實驗資料做比較，Akiyama假定在進口處無摻混現(xiàn)象，即γ=0，故得下式：

緩坡：

陡坡：

式中的異重流流速分布系數(shù)采用Ellison與Turner的實驗值。（Ellison與Turner用鹽水做實驗所得流速分布系數(shù)S1、S2值遠較Parker等［45］渾水異重流實驗測得值為?。?。Akiyama將此兩式與Singh的鹽水潛入實驗資料進行對比，實驗值與計算值符合。

3.8 Parker和Toniolo的理論分析成果 Parker和Toniolo［46］專門針對Akiyama與Stefan的論文［43］提出意見。他們指出：Akiyama文中在分析阻力系數(shù)和底部比降等影響時所假定潛入點下游的異重流運動在緩坡時的Fd值為正常值、在急坡時的Fd值為臨界值是錯誤的，也是沒有必要的。他們認為Akiya?ma的分析，用兩種CVⅠ和CVⅡ的動量守恒方程就足以闡明hd和hp的關系，而且能得到Fp和Fd兩者僅為摻混系數(shù)γ的關系式。

Parker分析CVⅠ的動量方程，導得以下兩式：

當γ=0時，式（22）可得與朱鵬程推導相同的式（16）。

從式（22）可計算k為Fp和γ的函數(shù)，如果hp、up、εγ（來水密度以及γ為已知，則hd、ud和εd（異重流密度可從式（22）、連續(xù)方程求得。

Parker根據(jù)CVⅡ的動量方程，導得：

將式（24）代入式（22）

Parker認為利用CVⅡ的動量方程，有了式（24）就可直接計算Fp，如僅有CVⅠ的動量方程，所求得的式（22）還須先知道Fp和γ（或Fd和γ）用于計算k值，而用CVⅡ動量方程求得的式（24）和式（25），F(xiàn)p和k兩者可從已知γ時求得。

4 異重流潛入現(xiàn)象的數(shù)學模型計算成果

方春明等［16］對水庫異重流潛入條件采用動量方程分析得：

根據(jù)兩種概化圖形，利用能量方程和動量方程進行理論分析，探討了異重流在水庫中潛入點密度Froude數(shù)小于1的水沙條件。方春明等還建立了二維水流和懸沙數(shù)學模型進行數(shù)值計算。為檢驗模型結果，選取曹如軒9組水槽潛入點實驗所測資料進行對比計算，計算得Fp=0.48～0.77，與實驗Fp=0.45～0.75比較，兩者接近。而且他們對影響水流和泥沙模型計算結果的幾個因素進行了潛入點流態(tài)和潛入點水深的對比計算；還計算了異重流潛入點的流場、清渾水交界面水深的縱向變化以及水面線的縱向變化，發(fā)現(xiàn)在壅水區(qū)水面線在潛入點附近出現(xiàn)倒比降。此計算結果與Savage提出的異重流潛入現(xiàn)象分析的概化圖形中下游斷面水位隆起高度Δ在定性上相符。

Farrell與Stafan［22，47］對水庫中異重流潛入現(xiàn)象用k～ε模型進行計算，模型包括水流動量方程和溫度方程，共進行12次紊流模型方程計算，得，即Fp=0.5。潛入點水深與Singh和Shah、Farrell和Stefan的實驗數(shù)據(jù)系數(shù)1.3相比大20%，F(xiàn)arrell和Stafan尚難解釋其差別的原因；其次這12次計算中，計算初期混合系數(shù)（用符號γ表示）在0.028～0.19之間，γ的變化可用下式表示，Q0為進槽流量，Qd為異重流流量。

Bournet等［48］的異重流潛入數(shù)值模型計算，應用水流動量方程、k-ε模型以及溫度輸移方程，計算流速場和溫度在等寬槽和水槽其中一邊擴展一個小的角度（1°～7°）兩種情況下的分布。在等寬槽的計算中引進摻混系數(shù)，其值在0.005～0.012之間，計算得潛入點水深此式hp計算值較Singh實驗值大些。

表4 水庫異重流潛入點理論分析與數(shù)值計算

表5 引航道口門異重流潛入點理論分析

以上各家理論分析和數(shù)值模型計算成果，匯總列于表4和表5。有關渾水異重流潛入實驗結果分析，將在本文第Ⅱ篇中給予介紹。

參 考 文 獻：

［1］范家驊.異重流與泥沙工程：實驗與設計［M］.北京：中國水利水電出版社，2011.

［2］范家驊，等.異重流運動的實驗研究［J］.水利學報，1959（5）：30-48.

［3］范家驊，等.異重流的研究和應用［M］.北京：水利電力出版社，1959.

［4］EGASHIRA S，ASHIDA K.Condition that suspended load plunges to form a gravity current［C］//19th Proc.of the Conference on Natural Disaster Science.Japan，1978.

［5］曹如軒，任曉楓，盧文新.高含沙異重流的形成與持續(xù)條件分析［J］.泥沙研究，1984（2）：1-10.

［6］曹如軒，錢善琪，郭崇，等.粗沙高含沙異重流的運動特性［J］.泥沙研究，1995（2）：64-73.

［7］俞維升.水庫沉滓運動特性之研究［D］.臺北：臺灣大學，1991.

［8］姚鵬，王興奎.異重流潛入規(guī)律研究［J］.水利學報，1996（8）：77-83.

［9］姚鵬.異重流運動的試驗研究［D］.北京：清華大學，1994.

［10］焦恩澤.黃河水庫泥沙［M］.鄭州：黃河水利出版社，2004.

［11］SCHIJF J B，SCHONFELD J C.Theoretical Considerations on the Motion of Salt and Fresh Water［C］//Proceed?ings of Minnesota International Hydraulics Convention.Minnesota：Minneapolis，1953：321-333.

［12］范家驊.關于水庫渾水潛入點判別數(shù)的確定方法［J］.泥沙研究，2008（1）：74-81.

［13］李濤，張俊華，馬懷寶，等.異重流潛入重力修正系數(shù)研究［J］.人民黃河，2012，34（7）：28-29.

［14］李書霞，夏軍強，張俊華，等.水庫渾水異重流潛入點判別條件［J］.水科學進展，2012，23（3）：363-368.

［15］李濤，夏軍強，張俊華，等.水庫異重流潛入點流速分布及其判別式改進［J］.工程科學與技術，2017，49（2）：62-68.

［16］方春明，韓其為，何明民.異重流潛入條件分析及立面二維數(shù)值模擬［J］.泥沙研究，1997（4）：68-75.

［17］KEULEGAN G H.Twelfth progress report on model laws for density currents：the motion of saline fronts in still water［R］.National Bureau of Standards：U.S.Dept.of Commerce，1958.

［18］SINGH B，SHAH C R.Plunging phenomenon of density currents in reservoirs［J］.La Houille Blanche，1971，26（1）：59-64.

［19］ITAKURA T，RISHI T.Study on the turbidity density current in a reservoir［C］//16th Symp.on Science of Natural Disaster.Japan，1979：233-234.

［20］FUKUOKA S，F(xiàn)UKUSHIMA Y，NAKAMURA K.Study on the plunge depth and interface form of density cur?rents in a two-dimensional reservoir［C］//土木學會論文報告集第302號 .1980：55-65.

［21］KAN K，TAMA N.On the plunging point and initial mixing of the inflow into reservoirs［C］//Proc.25th Japanese Conf.on Hydraulics.Japan，1981：631-636.

［22］FARRELL C J，STEFAN H G.Buoyancy induced plunging glow into reservoirs and coastal regions［R］.St.An?thony Falls Hydraulic Laboratory，Minnesota：University of Minnesota，1986.

［23］AKIYAMA J，STEFAN H.Gravity currents in lakes，reservoirs and coastal regions：Two-layer stratified flow analysis［R］.St.Anthony Falls Hydraulic Laboratory，Minnesota：University of Minnesota，1987.

［24］AKIYAMA J，STEFAN H.Onset of underflow in slightly diverging channels［J］.Journal of Hydraulic Engineer?ing-Asce，1987，113（7）：825-844.

［25］STEFAN H，JOHNSON T.Negative buoyant flow in diverging channel.Ⅲ.Onset of underflow［J］.Journal of Hydraulic Engineering，1989，115（4）：423-436.

［26］O'BRIEN M P，CHERNO J.Model law for motion of salt water through fresh［J］.Trans.ASCE，1934（99）：576-594.

［27］ YIH C S.A study of the characteristics of gravity waves at a liquid surface［D］.Iowa：State Univ.of Iowa，1947.

［28］KEULEGAN G H.An experimental study of the motion of saline water from locks into fresh water channels［C］//13th Progress Report on Model Laws of Density Current，National Bureau of Standards，1957.

［29］MEDDLETON G V.Experiments on density and turbidity currents.I.motion of the head［J］.Canadian J.of Earth Sciences，1966（3）：523-546.

［30］ROGOVSKY.International Symposium of Stratified Flows［C］.Novosibirsk，1972：440-458.

［31］KERSEY D G，HSü K.Energy relations of density current flows：an experimental investigation［J］.Sedimentolo?gy，1976（23）：761-789.

［32］LAM S T.Experimental investigation of lock-release gravity current［D］.Hong Kong：Univ.of Hong Kong，1995.

［33］BARR D I H.Densimetric exchange flow in rectangular channels.II：Some observations of the structure of lock exchange flow［J］.La Houille Blanche，1963（7）：757-766.

［34］ROTTMAN J W，SIMPSON J E.The initial development of gravity currents from fixed-volume release of heavy fluids［C］//IUTAM Symp.Delft，1983：347-359.

［35］陳國謙，李行偉，李植.開閘式湍動異重流［J］.中國科學（E輯），2002，32（6）：754-764.

［36］von KARMAN T.The engineer grapples with non-linear problems［J］.Bull.Am.Math.Soc.，1940（46）：615-683.

［37］BENJAMIN T B.Gravity current and related phenoma［J］.J.Fluid Mech.，1968（31）：209-248.

［38］YIH C S.Stratified Flows［J］.Annual Review of Fluid Mechanics，1969.

［39］SAVAGE S B，BRIMBERG J.Analysis of plunging phenomena in water reservoirs［J］.J.Hydraul.Res.，1975，13（2）：187-205.

［40］JAIN S C.Plunging phenomena in reservoirs［C］//Proc.Symp.On Surface Water Impoundments，Minnesota：Minneapolis，1981.

［41］朱鵬程.異重流的形成與衰減［J］.水利學報，1981（5）：52-59.

［42］AKIYAMA J，STEFAN H.Theory of plunging flow into a reservoir：Internal memoratum［R］.St.Anthony Falls Hydraulic Lab.，Univ.of Minnesota，1981.

［43］AKIYAMA J，STEFAN H.Plunging flow into a reservoir：theory［J］.J.Hydraul.Eng.，1984，110（4）：484 499.

［44］ELLISON T H，TURNER J S.Turbulent entrainment in stratified flows［J］.J.Fluid Mech.，1959（6）：423-428.

［45］PARKER G，GARCIA M，F(xiàn)UKUSHIMA Y，et al.Experiments on turbidity currents over an erodible bed［J］.J.Hydraul.Res.，1987，25（1）：123-147.

［46］PARKER G，TONIOLO H.Note on the analysis of plunging of density flows［J］.Journal of Hydraulic Engineer?ing，2007，133（6）：690-694.

［47］FARRELL C J，Stefan H.Mathematical modelling of plunging reservoir flows［J］.Journal of Hydraulic Re?search，1988，26（5）：525-537.

［48］BOURNET P，DARTUS D，TASSIN B，et al.Numerical investigation of plunging density current［J］.J.Hy?draul.Eng.，1999，125（6）：584-594.