劉文君^*，崔磊
(清華大學(xué)環(huán)境科學(xué)與工程系，北京，100084)

　　摘要: 本文介紹了應(yīng)用計(jì)算流體力學(xué)（Computational Fluid Dynamics，CFD）優(yōu)化給水廠清水池設(shè)計(jì)的原理和應(yīng)用方法，并以863示范工程為物理模型，通過(guò)建立中試裝置的CFD模型，得到清水池?cái)?shù)學(xué)模型的液齡分布函數(shù)和累計(jì)液齡分布函數(shù)。研究比較了CFD計(jì)算和實(shí)際中試模型示蹤試驗(yàn)的結(jié)果，表明二者的誤差在可接受的范圍之內(nèi)，因此證明應(yīng)用CFD方法優(yōu)化清水池設(shè)計(jì)是可行的。
　　關(guān)鍵詞：CFD模型；清水池；水力效率；停留時(shí)間

The Principles and Pilot Tests of Applying Computational Fluid Dynamics (CFD) Method in Optimizing the Design of Clearwell
LIU Wen-jun﹡, CUI Lei

Abstract: This paper introduces the principles and application procedures of the Computational Fluid Dynamics (CFD) method in optimizing the design of clearwell. In this study, a CFD modeling of the pilot-scale clearwell prototyped an 863 demo-project was established, which described the detailed information of the local characteristic of the flow domain which could not be visualized in trace tests. The residence time distribution (RTD) curves and cumulative residence time distribution (CRTD) curves of the pilot-scale clearwell were obtained, which was compared with those from the pilot-scale trace-tests. The errors between the results from CFD calculation and trace-tests were within the acceptable range, which demonstrated that CFD method can be used in optimizing the clearwell design.

Key words:CFD model; clearwell; hydraulic efficiency; residence time

　　目前國(guó)內(nèi)已開(kāi)始研究改善清水池的流體力學(xué)特性以提高消毒效果和控制消毒副產(chǎn)物的產(chǎn)生^[1,2]。CFD(Computational Fluid Dynamics)就是利用計(jì)算機(jī)求解流體流動(dòng)的各種守恒偏微分方程組的技術(shù)，其中涉及流體力學(xué)尤其是湍流力學(xué)、計(jì)算方法以及計(jì)算機(jī)圖形處理等技術(shù)。流體的流動(dòng)可用Navier－Stokes方程描述，因此CFD中最基本的數(shù)學(xué)物理模型即為粘性流體的Navier－Stokes方程。
　　目前在清水池流體力學(xué)特性研究中普遍通過(guò)示蹤劑試驗(yàn)求得停留時(shí)間函數(shù)模型(RTD)來(lái)表達(dá)清水池的水力效率^[1-3]。但通過(guò)示蹤試驗(yàn)測(cè)定清水池時(shí)間分布函數(shù)費(fèi)時(shí)費(fèi)力，而且有時(shí)可能不能實(shí)現(xiàn)。隨著電腦運(yùn)行能力和計(jì)算流體力學(xué)方法的不斷進(jìn)步，對(duì)流場(chǎng)和傳質(zhì)的高精度數(shù)值模擬已成為可能。通過(guò)CFD模擬計(jì)算可以很容易、精確且直觀地得出清水池內(nèi)水流流態(tài)和速度分布，揭示出很多物理模型無(wú)法表達(dá)的信息，如局部流速、短路流形態(tài)、單個(gè)顆粒停留時(shí)間和運(yùn)行軌跡等，優(yōu)化清水池設(shè)計(jì)。本研究以國(guó)家863示范工程深圳水務(wù)集團(tuán)筆架山水廠＃2清水池為模擬對(duì)象，采用CFD方法計(jì)算清水池在不同結(jié)構(gòu)配置條件下的水力學(xué)特征，并與中試模型試驗(yàn)結(jié)果進(jìn)行比較，在證明CFD方法可行的前提下，用CFD方法優(yōu)化給水廠清水池設(shè)計(jì)。本文主要介紹CFD方法原理和CFD計(jì)算結(jié)果與中試結(jié)果的比較，CFD方法的應(yīng)用結(jié)果將另文介紹。

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

1.1　 數(shù)學(xué)模型的整體框架
　　建立模型的步驟分為數(shù)學(xué)描述、數(shù)值求解、邊界和初始條件設(shè)置三部份。數(shù)學(xué)模型建立后便是數(shù)學(xué)模型的解析，包括迭代求解、數(shù)據(jù)采集與后處理、檢驗(yàn)?zāi)M準(zhǔn)確度三部分，流程圖如圖1。
　　圖1概括了建立數(shù)學(xué)模型的基本過(guò)程，整個(gè)過(guò)程分為建立數(shù)學(xué)模型和解析數(shù)學(xué)模型兩部分，而建立數(shù)學(xué)模型又分為幾何模型和物理模型兩方面，其中物理模型又包括連續(xù)相流體和分散相顆粒兩個(gè)子模型。

圖1 數(shù)學(xué)模型整體框架
Fig.1 　Mathematical model frameworks

1.2　 兩相流數(shù)學(xué)模型基本概念
1.2.1　 連續(xù)相（水相）紊流模型
　　由于一般流體處于紊流狀態(tài)，因此本模型采用紊流RNG k－ε 模型。它最重要的兩個(gè)方程就是關(guān)于兩個(gè)最重要參數(shù)－動(dòng)能k和能量耗散ε的方程。

　　

　　其中G_k是由層流速度梯度而產(chǎn)生的湍流動(dòng)能，G_b是由浮力而產(chǎn)生的湍流動(dòng)能，Y_M是由于在可壓縮湍流中過(guò)渡的擴(kuò)散產(chǎn)生的波動(dòng)，C₁，C₂，C₃是常量，α_k和α_ε是k方程和ε方程的湍流Prandtl數(shù)，S_k和S_e由用戶定義^[5]。
1.2.2分散相（示蹤劑）水力輸送模型
　　本文對(duì)顆粒相采用歐拉－拉格朗日分散相模型，對(duì)顆粒的運(yùn)行軌跡和停留時(shí)間進(jìn)行模擬，其最重要的方程是顆粒在流場(chǎng)中的受力方程^[4]：

　　　　 　　　　　(3)

　　其中u為連續(xù)相流速，u_p為顆粒相流速，ρ和ρ_p分別為流體和顆粒密度。F_d(u－u_p)表示了單位質(zhì)量顆粒所受的曳力。G_x(ρ_p－ρ)/ρ_p表示單位質(zhì)量顆粒所受的重力。F_x為顆粒受到的其他力。
1.3 　數(shù)學(xué)模型涉及的相關(guān)公式
　　CFD軟件根據(jù)以Navier－Stokes方程為基礎(chǔ)的流體流動(dòng)和能量平衡關(guān)系進(jìn)行計(jì)算求解，使用者可以根據(jù)自己的模型選擇特定的平衡關(guān)系包含到自己的模型中。實(shí)際上，根據(jù)模型合理選擇特定的平衡關(guān)系是一種合理的簡(jiǎn)化。
　　Fluent軟件中，這些平衡關(guān)系來(lái)源于Navier- Stokes方程中關(guān)于質(zhì)量、能量和動(dòng)量的守恒，還有如熱量傳遞、物質(zhì)混合或反應(yīng)、k－ε紊流模型等^[4]。
1.3.1　 質(zhì)量守恒方程
　　對(duì)于多組分系統(tǒng)，質(zhì)量守恒方程描述如下：

　　　　　 　　　　　(4)

　　其中m_i為系統(tǒng)中各組分的局部質(zhì)量分?jǐn)?shù)，J_i‘,i為組分i‘由于濃度梯度引起的擴(kuò)散通量，R_i‘為由于化學(xué)反應(yīng)而產(chǎn)生或損耗的質(zhì)量速率，S_i‘為其他源項(xiàng)的質(zhì)量速率。
1.3.2　 動(dòng)量守恒方程
　　在慣性（非加速）坐標(biāo)系中i方向上的動(dòng)量守恒方程為：

　　　　　　　(5)

　　其中p是靜壓，τ_ij是應(yīng)力張量，ρg_i和F_i分別為i方向上的重力體積力和外部體積力（如離散相相互作用產(chǎn)生的升力）。F_i包含了其他的模型相關(guān)源項(xiàng)，如多孔介質(zhì)和自定義源項(xiàng)。
1.4　 建立網(wǎng)格
　　如上所述，描述流場(chǎng)的一系列偏微分方程通過(guò)有限體積方法將流場(chǎng)劃分為眾多小區(qū)域后，再由求解器進(jìn)行計(jì)算求解^[5]。
　　網(wǎng)格的密度和質(zhì)量可以由用戶自己決定，并且影響最終的求解過(guò)程。過(guò)于粗糙的網(wǎng)格，往往只給出流場(chǎng)一些簡(jiǎn)單的輪廓，而無(wú)法給出詳細(xì)且精確的結(jié)果，甚至無(wú)法揭示出流場(chǎng)的本質(zhì)特征。
　　當(dāng)設(shè)置好模型幾何圖形的邊界條件后，就可以用上面提到的平衡方程和邊界條件對(duì)所有網(wǎng)格進(jìn)行迭代計(jì)算，從而得到特定問(wèn)題的收斂解。
1.5　 網(wǎng)格精度
　　由于網(wǎng)格決定著模型的準(zhǔn)確度，所以生成合理的網(wǎng)格是建立數(shù)學(xué)模型的一個(gè)關(guān)鍵步驟。生成網(wǎng)格不但要考慮其細(xì)密程度來(lái)提高模擬過(guò)程的細(xì)節(jié)準(zhǔn)確程度，又要避免由于網(wǎng)格過(guò)于復(fù)雜帶來(lái)的漫長(zhǎng)計(jì)算時(shí)間。一般來(lái)說(shuō)，需要找到一個(gè)最優(yōu)化的網(wǎng)格密度來(lái)滿足模型的計(jì)算精度和計(jì)算時(shí)間，再提高精度也不會(huì)更多地反映出需要的流場(chǎng)細(xì)節(jié)，而減少精度又會(huì)掩蓋住流場(chǎng)的本質(zhì)特征。
1.6 　網(wǎng)格結(jié)構(gòu)
　　基于有限體積方法，對(duì)于那些流場(chǎng)特性變化不明顯的區(qū)域，網(wǎng)格精度可以予以降低，而對(duì)于變化顯著的地方如邊界層、拐角處等地方，網(wǎng)格需要提高精度，減小有限體積，從而提高反映該區(qū)域流場(chǎng)特性變化的詳細(xì)信息的準(zhǔn)確性。圖2為基于示范工程＃2清水池的實(shí)際物理中試模型建立的幾何模型網(wǎng)格圖。
1.7　 物理原型
　　CFD模擬的中試清水池幾何尺寸和水力特性如表2。

圖2 　七擋板清水池?cái)?shù)學(xué)模型網(wǎng)格圖
Fig.2 Nets of the mathematical model of seven-baffle clearwell

　　從圖2中可以看出邊界層和拐角處的網(wǎng)格設(shè)置較細(xì)密。表1給出該網(wǎng)格的基本幾何信息。

表1 七擋板清水池?cái)?shù)學(xué)模型網(wǎng)格幾何信息
Table 1 Geometrical parameters of the mathematical model of seven-baffle clearwell

網(wǎng)格信息 設(shè)置條件 網(wǎng)格信息 設(shè)置條件 入口管徑/m 0.05 網(wǎng)格/個(gè) 30611 入口管長(zhǎng)/m 0.3 節(jié)點(diǎn)/個(gè) 15907 池長(zhǎng)/m 3 最小體積/m³ 4.618e-05 池寬/m 1 最大體積/m³ 6.817e-05 擋板長(zhǎng)/m 2.5 總體積/m³ 3.03

表2 物理模型幾何尺寸和水力特性
Table 2 Geometrical and hydraulic parameters of the pilot tests

　 長(zhǎng)度/m 寬度/m 高度/m 管徑/mm HRT/s 流量/m³/s 進(jìn)口流速/m/s 自由液面高度/m 中試模型 3 1 0.6 50 2497.5 4.44×10^-4 0.2266 0.37

　　無(wú)論對(duì)于物理模型還是數(shù)學(xué)模型，確定流場(chǎng)的T₁₀/T值^[1]和RTD曲線都需要進(jìn)行示蹤試驗(yàn)。示蹤試驗(yàn)方法見(jiàn)文獻(xiàn)^[2]。
1.8　 設(shè)置邊界和初始條件
　　對(duì)于數(shù)值模擬的求解過(guò)程，正確設(shè)置邊界和初始條件是使模擬得到收斂，特別是得到穩(wěn)態(tài)解的重要保障。
1.8.1　 進(jìn)水口流速
　　確定進(jìn)口處的流速則需要根據(jù)質(zhì)量守恒的原則，即單位時(shí)間內(nèi)三維物理模型進(jìn)入流場(chǎng)的質(zhì)量/體積等于單位時(shí)間二維數(shù)學(xué)模型進(jìn)入流場(chǎng)的質(zhì)量/體積，則有以下公式：

　　　　　　　　　(6)

　　其中，U為平均流速；
　　H為二維入水槽高度；
　　D為二維入水槽寬度。
　　物理模型入口管徑為0.05m，流量為1.6m³/h，水面高度0.37m。
　　可以算得二維模型的入口處進(jìn)水平均流速為

　　

1.8.2　 確定紊流模型參數(shù)
　　入口處紊流動(dòng)能k定義為：

　　k=aU²　　　　　　

　　單位質(zhì)量的能量耗散率ε定義為：

　　　　或　ε=1.5k^1.5/Γ_H　　　

　　其中ρ為流體密度，μ為流體粘度，C_μ為一常數(shù)0.09，R_μ為紊流和流體粘度之比，Г_H為水力半徑。
　　流場(chǎng)邊界設(shè)置見(jiàn)表3。

表3 數(shù)學(xué)模型邊界條件設(shè)置
Table 3 Boundary conditions of the mathematical model

變量 進(jìn)口 出口 墻面（外墻和擋板） 類型 入口速度 出口壓強(qiáng) 無(wú)滑壁面 速度U /m/s 0.0148 　 0 紊流動(dòng)能k 6.58e-07 0 0 能量耗散ε 2.814e-11 0 0 壓強(qiáng)P /　Pa 　 101325 　 對(duì)顆粒作用 約束 逃逸 約束

1.9　 求解過(guò)程
　　計(jì)算過(guò)程可以不將顆粒相對(duì)于流體（連續(xù)相）的影響考慮在內(nèi)，也可以考慮顆粒對(duì)流體的作用效果。這兩種算法分別稱為非耦合計(jì)算和耦合計(jì)算。由于考慮到顆粒相濃度很小，和連續(xù)相相比，顆粒對(duì)流體的作用可以忽略，所以本文的模型采取兩相非耦合模型進(jìn)行計(jì)算。
1.10　 數(shù)據(jù)采集和后處理
　　模型的兩相經(jīng)過(guò)求解計(jì)算后，可以得到每個(gè)控制體積內(nèi)兩相的各項(xiàng)屬性，包括流體的流速、壓強(qiáng)、顆粒的流速以及顆粒的軌跡和在流場(chǎng)中的停留時(shí)間。這里，對(duì)物理原型的采樣方法適當(dāng)簡(jiǎn)化，即以相同的時(shí)間步長(zhǎng)對(duì)顆粒出口情況進(jìn)行統(tǒng)計(jì)，來(lái)得到其時(shí)間停留密度曲線和累計(jì)液齡分布函數(shù)曲線。T₁₀則由F( t )上最接近F( t )=0.1的兩個(gè)點(diǎn)線性插值求得。

2　 模擬結(jié)果和準(zhǔn)確度驗(yàn)證

2.1　 流場(chǎng)的模擬
　　模擬的流場(chǎng)速度等值線圖見(jiàn)圖3。
　　如圖3所示，通過(guò)流場(chǎng)出入口的流體流速變化明顯。流場(chǎng)內(nèi)部，特別是擋板間的長(zhǎng)流道部分，流體均勻流特性較好，拐角部分流速變慢，越遠(yuǎn)離擋板，速度越小。擋板的設(shè)置可以增強(qiáng)流體的均勻性，但同時(shí)也會(huì)造成短路流，拐角部分矢量圖的局部放大圖如圖4。

圖3 數(shù)學(xué)模型速度等值線圖
Fig.3 Velocity field of mathematical model

圖4 數(shù)學(xué)模型（基準(zhǔn)）拐角部分矢量圖
Fig.4 Vector field of the mathematical model in turning area

　　從圖中可以看出，在靠近外墻的部分，造成了很多與流體速度逆向的短路流。短路流的存在會(huì)影響示蹤劑在流場(chǎng)內(nèi)的停留時(shí)間，進(jìn)入短路流的顆粒需要很長(zhǎng)時(shí)間才能擺脫，重新回到主流體的行進(jìn)過(guò)程當(dāng)中。
2.2　 示蹤劑的模擬
　　對(duì)于單個(gè)粒子來(lái)說(shuō)，進(jìn)入流場(chǎng)后，隨著流體流動(dòng)，直至從出口出逃逸。典型單個(gè)粒子在流場(chǎng)中受到短路流影響的軌跡見(jiàn)圖5。
　　顆粒進(jìn)入擋板拐角處外側(cè)和拐角后流體內(nèi)側(cè)速度很小的短路流后，增加了其停留時(shí)間，可見(jiàn)由于擋板設(shè)置造成的短路流對(duì)于流場(chǎng)的水力運(yùn)輸效率產(chǎn)生了負(fù)面的影響。表4給出了物理模型和數(shù)學(xué)模型的一些水力效率評(píng)價(jià)參數(shù)比較。圖6、圖7為CFD模擬計(jì)算結(jié)果與中試示蹤試驗(yàn)結(jié)果比較。可以看出，CFD模擬計(jì)算結(jié)果與試驗(yàn)實(shí)際測(cè)出的液齡分布函數(shù)和累接液齡分布函數(shù)基本一致。結(jié)果可以接受，但二者得出的T₁₀/T值誤差為26.7%，還需進(jìn)一步提高CFD的精度。

圖5 數(shù)學(xué)模型（基準(zhǔn)）單顆粒軌跡線圖
Fig.5 Single particle trace of the mathematical model

表4 數(shù)學(xué)模型模擬效果比較
Table 4 The comparison of the mathematical model and pilot tests

　 T/s T₁₀/T T₅₀/T T₉₀/T 物理模型 2497.5 0.635 0.855 1.067 數(shù)學(xué)模型 2497.5 0.805 0.906 1.225 百分比誤差/％ 0 26.7 6 14.8

圖6 數(shù)學(xué)、物理模型累計(jì)液齡分布函數(shù)比較
Fig.6 The comparison of the cumulative residence time distribution curves from CFD calculation and pilot tests

圖7 數(shù)學(xué)、物理模型液齡分布函數(shù)比較
Fig.7 The comparison of the cumulative residence time distribution curves from CFD calculation and pilot tests

3　 結(jié)論

　　本文介紹了采用計(jì)算流體力學(xué)（CFD）方法模擬計(jì)算清水池內(nèi)流體特性的原理、方法和應(yīng)用。通過(guò)比較中試示蹤試驗(yàn)和CFD計(jì)算結(jié)果證明CFD模擬計(jì)算的結(jié)果精度在可接受的范圍，可以應(yīng)用于清水池設(shè)計(jì)的優(yōu)化。

參考文獻(xiàn)：

[1] 劉文君. 清水池改進(jìn)原理與應(yīng)用 給水排水，2004, 30(5): 10~12.
[2] 金俊偉，劉文君. 影響清水池t₁₀/T值的因素試驗(yàn)研究 給水排水，2004, 30(12): 36~39.
[3] Gil F. Crozes, James P. Hagstrom, Mark M. Clark. Improving Clearwell Design for CT Compliance. AWWA Research Fundation and AWWA, 1997.
[4] C.N. Hass. Disinfection under Dynamic Conditions: Modification of Hom’s Model for Decay. Environ. Sci. Technol.1994, 28(6):1367~1369.
[5] Fluent Inc. Fluent 6.0 Software Documentation. 2001.
[6] Fluent Inc. Gambit 2.1.0 Software Documentation. 2001.

---

資助項(xiàng)目：國(guó)家863項(xiàng)目“南方地區(qū)飲用水安全保障技術(shù)（2002AA601120）”
作者簡(jiǎn)介：劉文君(1968～)，男，博士，副教授，研究方向：水質(zhì)科學(xué)與工程。
*通訊作者： E-mail:wjliu@tsinghua.edu.cn