李 鵬 陳 銘 李明然

在水利工程中，因水流水頭高、流速大、紊動強、水流表面或過流邊界突變，造成大量的空氣摻入水流中，形成水氣兩相流，稱之為摻氣水流。隨著計算機技術(shù)的迅速發(fā)展，計算流體力學(xué)的數(shù)值模擬技術(shù)有了很大發(fā)展。本文以Fluent6.2軟件為基礎(chǔ)，采用標準紊流模型，對明渠水氣兩相流進行了數(shù)值模擬研究。

1. 模型建立及邊界條件

根據(jù)實驗工況，計算實例如圖1所示。明渠坡度為45°，槽長16.50m，槽寬0.20m，入口水深為8.0cm，流量為100L/s。分別對明渠陡槽的4個斷面（離入口的距離分別為6m、9m、13m和15m）的流速、紊動、摻氣濃度能進行計算。

圖1 明渠水氣兩相流計算模型

數(shù)值模擬網(wǎng)格用Fluent的Gambit前處理軟件生成，采用四邊形非均勻結(jié)構(gòu)化網(wǎng)格。水流入口的邊界分為空氣進口和水進口兩部分，其中，空氣進口邊界設(shè)定為壓力邊界條件，水進口邊界設(shè)定為速度入口邊界，其壓力均為標準大氣壓。由于不知道出庫水深，氣體和水流的邊界無法分開，只能作為單獨出口邊界處理，采用壓力出口邊界。壁面采用無滑移邊界，即U=0、V=0，近壁區(qū)，雷諾數(shù)較小，在粘性底層用壁函數(shù)處理。

2.控制方程及數(shù)值計算方法

2.1 控制方程

連續(xù)方程：

式中：下標a表示氣相，w表示水相，α 為容積分數(shù)，ρk為相密度，Uk為相速度，Γk為界面質(zhì)量傳輸率。

動量方程：

式中：p為壓強，g為重力加速度，μkeff為有效粘性系數(shù)，μkeff=μk+μt， 各相渦粘性系數(shù) μt=Cμtρk2/ε，Mk為相間作用力。

k方程：

其中：

ε方程：

2.2 數(shù)值計算方法

采用控制體積法對偏微分形式的控制方程組進行離散，得到代數(shù)形式的方程組，然后采用迭代法對方程求解。其中離散采用標準格式，壓力與速度的耦合采用SIMPLE方法計算。當計算各變量殘差小于1.0×10-4或殘差隨迭代次數(shù)的增加基本不變化時，可認為計算收斂。

3. 計算結(jié)果分析

3.1 摻氣濃度

圖2不同斷面摻氣濃度分布圖1-下游6m，2-下游9m，3-下游 13m，4-下游 15m

計算得到的4個斷面的摻氣濃度分布如圖2所示。從圖中可以看出，摻氣濃度沿程增加，并且在明渠底部增加幅度較小，而在上部增加幅度較大。每個斷面摻氣濃度變化特性是自渠底向水面逐漸增加。

3.2 流速分布

圖3 不同斷面流速分布圖1-下游6m，2-下游9m，3-下游13m，4-下游15m

從圖3中可以看出，計算的摻氣水流流速沿程逐漸增大，同一斷面流速自底面向水面不斷增大，隨后逐漸減小，與理論預(yù)測一致。流速自底面向水面方向不斷增大，其原因為邊壁的阻力作用，摻氣只能影響其變化梯度，在靠近水面的部分區(qū)域，流速沿垂向不斷減小，主要是受到摻氣的影響，摻氣使氣體邊界對水流速度的影響深入到了水流內(nèi)部。

3.3 紊動能

從圖4中看出水相紊動能在靠近壁面處大，為平均動能的90%左右，靠近水面處較小，且沿程有增大趨勢，且增加非常緩慢，這與理論分析結(jié)果相一致。從圖5可看出氣相紊動能在不同斷面上的分布情況。圖中可以看出，水相與氣相的紊動能分布規(guī)律相同，大小極為接近。

圖5 氣相紊動能在不同斷面的分布圖1-下游6m，2-下游9m，3-下游13m，4-下游15m

4. 結(jié)論

圖4 水相紊動能在不同斷面的分布圖1-下游 6m，2-下游9m，3-下游 13m，4-下游 15m

一是本文采用Fluent軟件中的紊流模型，利用VOF方法確定自由表面計算水氣兩相流是合理的。二是采用Fluent軟件計算的水相及氣相紊動能在靠近壁面處較大，隨水深增加逐漸減小。由于沒有紊動能分布的試驗資料，因此準確性還需進一步驗證。