城市污水是一種具有非牛頓特性的冪率流體，在其中含有大量的污物，包括大小尺度泥沙類(lèi)和懸浮纖維類(lèi)，流動(dòng)特性較水有顯著差異。考慮到污水中包括的固相污物特殊性，并且為了分析和研究方便，將城市污水視為均質(zhì)流體，采用單相非牛頓流體模型進(jìn)行分析研究，將固相污物的影響歸結(jié)到流體黏度的影響之中。由于V型調(diào)節(jié)球閥閥芯的V型切口（見(jiàn)圖1）和球閥密封面之間具有剪切作用，可以切斷纖維狀的流體，密封性能好，因此可用于城市污水流動(dòng)特性研究。該型調(diào)節(jié)球閥還具有流體阻力小、流通能力大的特點(diǎn)。調(diào)節(jié)閥內(nèi)部結(jié)構(gòu)復(fù)雜，而且還在不斷的變化更新，隨著計(jì)算機(jī)技術(shù)和數(shù)值計(jì)算的發(fā)展，運(yùn)用CFD技術(shù)對(duì)流體流動(dòng)進(jìn)行模擬仿真分析已經(jīng)成為一種普遍應(yīng)用的工程分析手段，通過(guò)CFD通用軟件fluent對(duì)某一型號(hào)的V型調(diào)節(jié)球閥內(nèi)部流體介質(zhì)的流動(dòng)特性進(jìn)行模擬仿真，分析不同性質(zhì)流體在調(diào)節(jié)閥內(nèi)部的流阻情況。

圖1 V型調(diào)節(jié)球閥的閥芯結(jié)構(gòu)

    1 流體計(jì)算模型

    在研究單相介質(zhì)流體流動(dòng)特性時(shí)，運(yùn)動(dòng)規(guī)律遵循三個(gè)基本守恒定律，即質(zhì)量守恒定律、動(dòng)量守恒定律及能量守恒定律，該課題對(duì)流體進(jìn)行數(shù)值模擬研究時(shí)采用適合工程問(wèn)題的k－ε湍流模型，對(duì)流體研究的控制方程如下：

    連續(xù)性方程為

         （1）

    動(dòng)量方程為

        （2）

    湍動(dòng)能k的輸運(yùn)方程為

        （3）

    湍流耗散率ε的輸運(yùn)方程為

        （4）

    湍動(dòng)能和湍流耗散率生成項(xiàng)Gk和Gε的計(jì)算公式為

        （5）

    其中：υt為渦黏性系數(shù)．在用進(jìn)行數(shù)值計(jì)算時(shí)，k－ε模型中經(jīng)驗(yàn)常數(shù)的取值通常情況為Cμ＝0.09，c1＝1.44，c2＝1.92，σk＝1.0，σε＝1.3。

    非牛頓冪率流體本構(gòu)方程為

       （6）

    其中：τ為剪切應(yīng)力（Pa）；K為稠度系數(shù)（N•sⁿ/m²）；γ為應(yīng)變速度（s^-1）；n為流動(dòng)指數(shù)。其中當(dāng)0＜n＜1時(shí)為剪切細(xì)化流；當(dāng)n＞1時(shí)為剪切稠化流。

    2 數(shù)值模擬邊界條件設(shè)置

    通過(guò)三維建模軟件solidworks建立V型調(diào)節(jié)球閥在不同開(kāi)度下的流道模型（模型參數(shù)由某公司提供），將建立的流道實(shí)體模型導(dǎo)入fluent前處理軟件Gambit中進(jìn)行數(shù)值計(jì)算前處理，建立流道的一半作為實(shí)體計(jì)算模型。為了提高計(jì)算精度，劃分網(wǎng)格時(shí)采用兼容性較強(qiáng)的非結(jié)構(gòu)化四面體網(wǎng)格，劃分后的網(wǎng)格數(shù)大約為26萬(wàn)。當(dāng)閥芯開(kāi)度為50℃時(shí)流道網(wǎng)格模型如圖2所示，進(jìn)口設(shè)置為壓力入口，出口設(shè)置為壓力出口。

圖2 開(kāi)度為50℃時(shí)流道網(wǎng)格模型

    3 數(shù)值求解

    設(shè)置進(jìn)出口壓差分別為0.2MPa、0.5MPa、1MPa、1.5MPa時(shí)，分別對(duì)牛頓流體（水）和非牛頓冪律流體（污水）剪切細(xì)化流和剪切稠化流通過(guò)調(diào)節(jié)閥不同開(kāi)度下的流體流動(dòng)進(jìn)行數(shù)值模擬計(jì)算，對(duì)非牛頓流體進(jìn)行數(shù)值模擬計(jì)算時(shí)，通過(guò)以下步驟調(diào)出隱藏的非牛頓流體計(jì)算模型：define＞models＞vis－cous＞turbulence－expert＞turb－non－newtonian＞Y，湍流模型采用標(biāo)準(zhǔn)k－ε模型。當(dāng)壓差為0.5MPa時(shí)壓力分布云圖分別如圖3、圖4、圖5所示。

    從圖3、圖4、圖5中可以看出，由于在節(jié)流處產(chǎn)生渦漩流動(dòng)，從而引起流動(dòng)速度變化，流體流動(dòng)速度變化又引起該處壓力變化，因此在調(diào)節(jié)閥閥芯和流道節(jié)流處壓力變化*明顯，對(duì)比入口處壓力可以看出在節(jié)流處壓力值有所減小，主要是因?yàn)樵诠?jié)流處流體流動(dòng)速度突然增大從而引起壓力突然減小，不同性質(zhì)流體流經(jīng)調(diào)節(jié)閥時(shí)壓力變化比較明顯部位都集中在節(jié)流處，在調(diào)節(jié)閥入口和出口端流體流動(dòng)平穩(wěn)，壓力值沒(méi)有明顯變化。

圖3 開(kāi)度50℃壓差為0.5MPa的剪切細(xì)化流壓力云圖

圖4 開(kāi)度50℃壓差為0.5MPa的剪切稠化流壓力云圖

圖5 開(kāi)度50℃壓差為0.5MPa時(shí)水的壓力云圖

    4 流阻分析

    閥門(mén)的流阻系數(shù)ξ取決于閥門(mén)產(chǎn)品的尺寸、結(jié)構(gòu)以及內(nèi)腔形狀，其內(nèi)部每一個(gè)元件都可以看作產(chǎn)生阻力的元件系統(tǒng)，當(dāng)有流體通過(guò)閥門(mén)時(shí)，整個(gè)流動(dòng)系統(tǒng)中的流阻系數(shù)還和流體本身的特性有關(guān)．流體在閥門(mén)中的流阻系數(shù)計(jì)算公式為

         （7）

    其中：Δp為閥門(mén)前后的壓力損失（Pa）；ξ為閥門(mén)流阻系數(shù)；ρ為流體密度（kg/m³）；u為流體在管道內(nèi)的平均流速（m/s）。

    通過(guò)數(shù)值模擬得到壓差分別為0.2MPa、0.5MPa、1MPa、1.5MPa時(shí)，調(diào)節(jié)閥不同開(kāi)度下的流阻系數(shù)值，見(jiàn)圖6、表1、表2、表3。

表1 不同壓差下水流阻系數(shù)ξ

表2 不同壓差下污水（n＝0.1666）流阻系數(shù)ξ

表3 不同壓差下污水（n＝1.666）流阻系數(shù)ξ

    從圖6中我們可以看出，調(diào)節(jié)閥流阻系數(shù)隨開(kāi)度增大呈現(xiàn)出減小趨勢(shì)，當(dāng)調(diào)節(jié)閥為全開(kāi)狀態(tài)時(shí)流阻系數(shù)*小。閥門(mén)開(kāi)度較小時(shí)，由于閥芯、閥座所造成的流通面積小，形成局部阻力，并且在節(jié)流處會(huì)形成渦旋流動(dòng)，因此在開(kāi)度較小時(shí)流阻系數(shù)值較大。

圖6 壓差為0.5MPa時(shí)不同開(kāi)度下剪切細(xì)化流、

剪切稠化流和水流阻系數(shù)變化曲線

    對(duì)比表1、表2和表3的數(shù)據(jù)，可以看出對(duì)于不同流體，閥門(mén)流阻系數(shù)值不同，當(dāng)具有牛頓流體性質(zhì)水通過(guò)閥門(mén)時(shí)，流阻系數(shù)較非牛頓性質(zhì)污水小。由于牛頓流體的剪切應(yīng)力和應(yīng)變速度呈線性關(guān)系，黏度隨流體流動(dòng)不發(fā)生變化；而非牛頓流體剪切應(yīng)力隨應(yīng)變速度變化呈現(xiàn)非線性變化。黏度變化與剪切應(yīng)力有關(guān)，受污水所含固體污物顆粒影響，流動(dòng)過(guò)程中會(huì)產(chǎn)生較大黏性作用，黏度變化引起流阻改變。因此，具有非牛頓特性污水較牛頓流體特性水流阻系數(shù)值有明顯不同。

    對(duì)比表2和表3可以看出，在非牛頓流體流動(dòng)指數(shù)變化情況下，當(dāng)開(kāi)度較大時(shí)，由于流體流動(dòng)所受剪切應(yīng)力變化較小，因此稠化流流阻系數(shù)較細(xì)化流流阻系數(shù)值變化不明顯；在小開(kāi)度時(shí)，由于過(guò)流面積小引起剪切應(yīng)力變化大，黏度變化比較明顯，因此剪切稠化流流阻系數(shù)較剪切細(xì)化流變化明顯，其流阻系數(shù)值較剪切細(xì)化流流阻系數(shù)值大。

    5 結(jié)論

    通過(guò)模擬分析不同開(kāi)度和不同壓差條件下不同性質(zhì)流體通過(guò)V型調(diào)節(jié)球閥時(shí)的流動(dòng)特性可以看出：

    （1）牛頓流體特性水通過(guò)V型調(diào)節(jié)球閥時(shí)流阻系數(shù)比非牛頓特性污水通過(guò)調(diào)節(jié)閥時(shí)的流阻系數(shù)小。對(duì)比剪切細(xì)化流和剪切稠化流可以看出在小開(kāi)度時(shí)，稠化流通過(guò)調(diào)節(jié)閥時(shí)流阻系數(shù)較細(xì)化流大，在開(kāi)度較大時(shí)，流阻系數(shù)隨流動(dòng)指數(shù)變化不明顯；

    （2）通過(guò)數(shù)值模擬結(jié)果可以看出，當(dāng)調(diào)節(jié)閥開(kāi)度增大時(shí)，流體通過(guò)閥門(mén)時(shí)的流阻系數(shù)逐漸減小，在開(kāi)度*大時(shí)流阻系數(shù)值*小；

    （3）在同一開(kāi)度下，當(dāng)壓差變化時(shí)，不同性質(zhì)流體通過(guò)調(diào)節(jié)閥時(shí)流阻系數(shù)值變化不明顯。