王 冠 ，鄧建飛 ，寇琳媛 ，金鵬飛 ，朱學(xué)軍 ，馬玉山

（1.寧夏大學(xué) 機(jī)械工程學(xué)院，銀川 750000；2.寧夏大學(xué) 寧夏智能裝備 CAE 重點實驗室，銀川 750021；3.吳忠儀表有限責(zé)任公司，寧夏吳忠 751100）

0 引言

隨著節(jié)能減排工程的發(fā)展，高壓減壓閥在工程領(lǐng)域得到了廣泛的應(yīng)用［1-3］。其中多級減壓調(diào)節(jié)閥以流量控制穩(wěn)定，耐高壓等優(yōu)點廣泛地應(yīng)用在各種工業(yè)系統(tǒng)中［4-5］。多級降壓節(jié)流元件是多級降壓調(diào)節(jié)閥關(guān)鍵控壓結(jié)構(gòu)，其流道窄而壓力大，對其進(jìn)行結(jié)構(gòu)改進(jìn)以適應(yīng)工程領(lǐng)域嚴(yán)苛的工作條件、延長使用壽命對工程領(lǐng)域有重要意義。

目前，降壓性能一直是調(diào)節(jié)閥研究的熱點，圍繞實際工程問題，眾多學(xué)者對調(diào)節(jié)閥降壓結(jié)構(gòu)開展了深入的研究工作。侯聰偉等［6］對高壓減壓閥進(jìn)行了參數(shù)化分析，分析了相對角度、孔板厚度、孔板數(shù)量和板孔直徑對多級降壓調(diào)節(jié)閥的控制性能的影響，結(jié)果表明相對角設(shè)為時，蒸汽流經(jīng)多孔帶冠閥芯，可以獲得最大的減壓壓力，紊流度最低。高佳男等［7］針對氣體低壓中小口徑管道流量調(diào)節(jié)，設(shè)計了一種新的軸流式流量調(diào)節(jié)閥，同時采用計算流體力學(xué)（CFD）方法對所設(shè)計的調(diào)節(jié)閥內(nèi)部流動和流量調(diào)節(jié)特性進(jìn)行了仿真分析，得到流量特性曲線，結(jié)果表明，所設(shè)計的調(diào)節(jié)閥流量調(diào)節(jié)特性與設(shè)計曲線相符。黃凱等［8］基于標(biāo)準(zhǔn)K-ε湍流模型對四組不同閥芯結(jié)構(gòu)（平底、小弧、大弧、波浪形）調(diào)節(jié)閥流阻性能進(jìn)行研究，發(fā)現(xiàn)大弧形與波浪形閥芯結(jié)構(gòu)調(diào)節(jié)閥的截阻性能較平底閥芯調(diào)節(jié)閥的流阻系數(shù)值減小9.71%，同時閥芯處于全開度時，平底閥芯截阻性能具有一定的優(yōu)勢。姚世衛(wèi)等［9］利用FLUENT流體計算仿真軟件對不同開度時閥口的流量特性進(jìn)行了數(shù)值模擬，并通過試驗對閥口的流量特性進(jìn)行了驗證。

眾多學(xué)者通過計算流體力學(xué)的方法對閥門的內(nèi)部流場進(jìn)行了數(shù)值模擬分析研究。宋學(xué)官等［10］建立了直接操作安全泄壓閥（SRV）精確的CFD模型，模擬了系統(tǒng)從閥門開啟到關(guān)閉的全過程，得到了可壓縮流體通過SRV的詳細(xì)圖像，包括閥座區(qū)域的小尺度流動特征。錢金元等［11］研究了新型高壓減壓閥（HPRV）及其連接管道內(nèi)多級多孔板的馬赫數(shù)，分析了可逆等熵過程的馬赫數(shù)，提出了HPRV多級多孔板的設(shè)計方法，并進(jìn)行了數(shù)值模擬。NIKO HERAKOVI? 等［12］基于數(shù)值模擬和實驗分析，提出改變閥室和閥芯的幾何形狀，可以顯著降低流體力的軸向分量，從而降低必要的驅(qū)動力。YOUNG JOON AN等［13］利用CFD方法對LNG船用系統(tǒng)抗汽蝕控制閥進(jìn)行了三維不可壓縮湍流流動的數(shù)值分析，研究和分析了具有復(fù)雜流場的調(diào)節(jié)閥的流動特性，通過CFD分析顯示出改進(jìn)的流型，減少了汽蝕。

通過對以上文獻(xiàn)進(jìn)行分析可以知道降壓結(jié)構(gòu)對閥門控制性能起決定性作用，然而具體針對多級減壓閥降壓結(jié)構(gòu)部分進(jìn)行結(jié)構(gòu)優(yōu)化和分析的文獻(xiàn)較少。本文通過數(shù)值模擬，分析串式多級降壓調(diào)節(jié)閥閥芯的流量特性，以及研究不同降壓閥芯參數(shù)（轉(zhuǎn)彎倒角和流道傾角）對內(nèi)部流場的影響。為改善節(jié)流元件的節(jié)流效果提供理論依據(jù)，對其他串式多級降壓調(diào)節(jié)閥設(shè)計也有一定的借鑒意義。

1 數(shù)值計算方法

1.1 控制方程與湍流模型

流體控制方程與湍流模型對內(nèi)部流場數(shù)值模擬分析有重要影響。在分析多級降壓閥內(nèi)部流動時，首先要求解質(zhì)量守恒方程和動量守恒方程，此外，在進(jìn)行流阻特性分析時應(yīng)考慮流量系數(shù)以及流阻系數(shù)方程。

連續(xù)性方程：

動量方程：

流量系數(shù)：

流阻系數(shù)：

式中 v ——運動黏度；

Q —— 流量，m3/h，Q=Avin；

A——管道橫截面積，m2；

vin——管道進(jìn)口流速，m/s；

ρ ——對應(yīng)的介質(zhì)密度，kg/m3；

ρ0—— 15 ℃時水的密度，kg/m3，其中ρ0=999 kg/m3；

ΔP —— 閥前后壓差，Pa；

va——閥門平均流速，m/s。

通過多級降壓調(diào)節(jié)閥閥內(nèi)介質(zhì)的流動狀態(tài)為湍流流動，因為雷諾數(shù)數(shù)值大于106。RNG K-ε湍流模型與標(biāo)準(zhǔn)k-ε模型相似，考慮了渦流對湍流的影響，提高了渦流計算精度，同時RNG K-ε湍流模型被證明能更準(zhǔn)確地描述高壓降壓閥內(nèi)部的復(fù)雜流動［10-24］，因此本次數(shù)值模擬采用RNG K-ε湍流模型，其對應(yīng)的湍動能K和湍流耗散率ε方程如下：

其中Rε為：

有效湍流黏度 μeあ為：

式中 Gk—— 由平均速度梯度引起的湍動能生成項；

其中經(jīng)驗常數(shù)αk= αε=1.393，Cε=1.68；經(jīng)驗常數(shù)修正值 C1ε=1.42；Cμ=0.084 5；η =SK/ε，η0=4.38，β=0.012。

1.2 幾何結(jié)構(gòu)與計算模型

圖1（a）示出串式多級降壓調(diào)節(jié)閥的結(jié)構(gòu)，主要部件包括閥體、閥芯、閥芯套、閥桿等。其中流體流向為左進(jìn)右出，流體介質(zhì)從閥芯底部的凹口流入閥門降壓結(jié)構(gòu)，從閥芯套上頂端4個圓孔流出，完成對流體的降壓，主要降壓區(qū)域為閥芯與閥芯套組成的降壓結(jié)構(gòu)，如圖1（b）所示，降壓結(jié)構(gòu)由3個凹口以及出口組成，每個凹口由四級降壓結(jié)構(gòu)組成，包括轉(zhuǎn)彎和流道傾角部分，如圖1（c）所示。

圖1 串式多級降壓閥結(jié)構(gòu)及流道示意Fig.1 Schematic diagram of structure and flow channel of the tandem multistage step-down valve

由于本次數(shù)值模擬只對閥芯降壓結(jié)構(gòu)部分進(jìn)行參數(shù)化分析，故只對閥芯流道部分建模，即閥芯與閥芯套組成的降壓結(jié)構(gòu)部分，其結(jié)構(gòu)如圖1（b）所示。圖1（c）示出多級降壓調(diào)節(jié)閥的每一級凹口流道，對圖1（c）所示的降壓結(jié)構(gòu)建立不同參數(shù)的降壓結(jié)構(gòu)流道模型，參數(shù)包括轉(zhuǎn)彎倒角L和流道傾角θ。

表1為多級降壓調(diào)節(jié)閥閥芯結(jié)構(gòu)參數(shù)，建立不同參數(shù)的流道模型，分析參數(shù)對流體流動的影響，并建立不同開度的流道模型，以分析其流阻特性。參數(shù)分析中流道倒角為 2，4，6，8，10 mm，傾角為36°，48°，60°，72°。在進(jìn)行傾角分析時，倒角均設(shè)置為0 mm，進(jìn)行倒角分析時，傾角均設(shè)置為48°。

表1 流道的參數(shù)分析Tab.1 Analysis of flow channel parameters

1.3 邊界條件和網(wǎng)格無關(guān)性

閥門中的流體介質(zhì)是渣油加氫裝置中的熱高分油。邊界條件采用壓力進(jìn)口和壓力出口。進(jìn)口壓力為18.7 MPa，出口壓力為2.9 MPa，流體介質(zhì)密度為 680.9 kg/m3，動力黏度為 1.4×10-4kg/ms。由于結(jié)構(gòu)的對稱性，用實際流場的一半進(jìn)行計算，對稱中心平面作為對稱面，其他表面均設(shè)置為光滑無滑移壁面。

網(wǎng)格劃分時，網(wǎng)格質(zhì)量和網(wǎng)格大小均是影響計算精度的重要條件。采用ANSYS mesh劃分網(wǎng)格，因為降壓結(jié)構(gòu)部分的復(fù)雜性，所以采用自動網(wǎng)格劃分的方式進(jìn)行網(wǎng)格劃分。在數(shù)值模擬之前，需要對網(wǎng)格進(jìn)行獨立性檢驗。通過改變單元尺寸的大小，調(diào)整計算區(qū)域的網(wǎng)格數(shù)目，得到合適的網(wǎng)格劃分。選擇原結(jié)構(gòu)進(jìn)行網(wǎng)格獨立性檢驗，如圖2所示，在網(wǎng)格尺寸為1.5 mm和1 mm時，出口的質(zhì)量流量相對誤差為0.05%，故認(rèn)為可以進(jìn)行數(shù)值模擬。由于考慮到降壓結(jié)構(gòu)部分的復(fù)雜性，所以本次數(shù)值模擬采用網(wǎng)格大小為1 mm進(jìn)行計算。

圖2 網(wǎng)格無關(guān)性驗證Fig.2 Grid independence verification

2 結(jié)果與分析

2.1 流量特性

對不同開度下閥芯的流道結(jié)構(gòu)，分析其流動特性?；诓煌_度的穩(wěn)態(tài)仿真結(jié)果，分別做出流量系數(shù)和流阻系數(shù)的特性曲線，如圖3，4所示。從圖3，4中可以看出，隨著開度的變化，流量系數(shù)和流阻系數(shù)有著相反的變化趨勢，流量系數(shù)隨著開度的增大而增大，流阻系數(shù)隨著開度的增大而減小，這一變化也可以從理論公式推導(dǎo)得出，將流量系數(shù)和流阻系數(shù)公式聯(lián)立得：

圖3 流量系數(shù)特性曲線Fig.3 Characteristic curve of flow coefficient

圖4 流阻系數(shù)特性曲線Fig.4 Characteristic curve of flow resistance coefficient

圖5 流量系數(shù)與流阻系數(shù)關(guān)系曲線Fig.5 Relation curve of flow coefficient and flow resistance coefficient

從圖3，4還可以看出，多級降壓調(diào)節(jié)閥在中小開度下有較好的線性流量特性，隨著開度的增大，在開度達(dá)到80%之后，開度繼續(xù)增大，其流動阻力變化不大，因此對應(yīng)的流量系數(shù)變化較小。對于流阻系數(shù)，開度在40%之前，因其閥芯和閥芯套的結(jié)構(gòu)特點，隨著開度的增大，流阻系數(shù)迅速減小，在40%開度以上時，流阻系數(shù)變化不明顯。

從單級流道示意圖中可以看出流道的降壓部分主要由轉(zhuǎn)彎和拐角組成，為了在高壓差等惡劣工況下獲得更好的減壓過程，研究不同結(jié)構(gòu)參數(shù)對其內(nèi)部流動特性的影響。主要研究節(jié)流元件的結(jié)構(gòu)參數(shù)，有流道傾角和流道倒角兩部分。

2.2 流道傾角對流體流動的影響

流道傾角會改變流體在流道中的流動狀態(tài)。因此，針對流道傾角對流體在多級降壓結(jié)構(gòu)中內(nèi)部流動的影響開展研究。

圖6示出不同流道傾角對稱平面內(nèi)的壓力分布。當(dāng)流體進(jìn)入節(jié)流結(jié)構(gòu)時，其流動截面迅速收縮。此時流體速度增加，壓力降低。根據(jù)能量守恒定律，流體在節(jié)流結(jié)構(gòu)中流動時，由于彎道近壁面會形成渦流，從而造成能量的消耗，壓力無法回到原來的值。從圖中可以看出，當(dāng)流道傾角小于48°時，近壁面的渦流區(qū)明顯增大；當(dāng)流道傾角大于48°時，隨著流道傾角的增大，近壁面的渦流區(qū)明顯減小。

圖6 不同流道傾角對稱面壓力分布Fig.6 Pressure distribution on symmetrical surfaces with different flow channel inclination angles

圖7示出不同角度X方向上的壓力分布。從壓力分布圖中可以看出，節(jié)流結(jié)構(gòu)對流體的整體壓降達(dá)到逐級降壓的效果。流道傾角從36°變化到72°時，節(jié)流結(jié)構(gòu)內(nèi)的流動阻力不同，壓降也發(fā)生變化。其中，當(dāng)x=0.08 m時，流道傾角越大，第一凹口的壓降越大，且此后每級凹口都表現(xiàn)出相同的變化趨勢，即流道傾角越大，每一級凹口分配的壓降越大，在流道傾角為72°時，每一凹口的壓降達(dá)到最大，在凹口處完成主要降壓，從而使得在出口處壓降減小，進(jìn)而能減小在出口處產(chǎn)生空化的概率。同時，由圖可知，隨著流道傾角的增加，在每一凹口處，壓力恢復(fù)逐漸增大，當(dāng)流道傾角為48°時每一凹口的壓力恢復(fù)沒有另外3個角度時劇烈，且壓降比較平穩(wěn)。

圖7 不同角度X方向壓力分布示意Fig.7 Schematic diagram of pressure distribution in X direction at different angles

在不同流道傾角的情況下，多級降壓閥節(jié)流結(jié)構(gòu)內(nèi)部速度、湍流動能和湍流耗散率的最大差值見表2。從表中可以看出，傾角在48°時，速度、湍流動能和湍流耗散率的最大差值均為最小值，由以上可以確定的是，傾角在48°時為最佳的流道傾角值。

表2 不同流道傾角的速度和湍流參數(shù)最大差值Tab.2 Maximum difference values of velocity and turbulence parameters with different channel inclination angles

2.3 流道倒角對流體流動的影響

倒角的目的是將直線過渡變?yōu)槠交^渡。與直線過渡相比，平滑過渡可以使拐角處的流體流動更加均勻，湍流強(qiáng)度更小，因此倒角可以幫助改變多級降壓調(diào)節(jié)閥內(nèi)部的流體流動狀態(tài)。倒角參數(shù)分析為從0 mm增加到10 mm，以2 mm為間隔，分別分析其對多級降壓調(diào)節(jié)閥內(nèi)流體流動的影響。

圖8示出不同倒角長度對稱面單級凹口速度分布。由圖可知，隨著倒角的增大，流道近壁面的渦旋區(qū)域在逐漸減小，但是卻在流道倒角處形成了高速射流，其最大速度分別為100.11，124.31，124.01，129.04，121.56，127.17 m/s。隨著倒角的增大，流道倒角處射流區(qū)域在逐漸增大，這是因為隨著倒角的增大，流體進(jìn)入倒角部分的緩沖區(qū)域減少，在倒角部分造成的壓力損失減小，因此導(dǎo)致射流區(qū)域逐漸增大。

圖8 不同倒角長度對稱面單級凹口速度分布Fig.8 Velocity distribution of single notch on symmetrical planes of different chamfered lengths

圖9示出了不同倒角出口處的平均速度，從圖中可以看出，隨著倒角的逐漸增大，出口處的平均速度在逐漸減小，且近似呈線性關(guān)系，經(jīng)線性擬合得到：

圖9 不同倒角出口處的平均速度曲線Fig.9 Average velocity curves at different chamfer exits

式中 v ——出口的平均速度，m/s；

L ——倒角的長度，mm。

從擬合公式可知：隨著倒角的增大，出口處的平均速度在呈線性關(guān)系的減小。

圖10示出了不同倒角長度下多級降壓調(diào)節(jié)閥閥芯內(nèi)部湍流動能和湍流耗散率的最大差值。

圖10 湍流參數(shù)最大差值曲線Fig.10 Maximum difference curve of turbulence parameters

從圖中可以看出，在流道倒角為6 mm時閥芯內(nèi)部湍流動能和湍流耗散率的最大差值均為最小值，分別為 28.86，831 511.1 m2/s2，其中倒角設(shè)置為6 mm時，相對于倒角為0 mm時，平均湍動能降低73%，極大程度的降低了流體介質(zhì)的湍流度。倒角在8 mm和10 mm時湍流動能和湍流耗散率突然增大是因為倒角增大，導(dǎo)致流體在倒角處速度的急劇增加，形成高速射流，使得流體湍流程度變得非常劇烈。圖11示出不同倒角時湍動能平均值與均方差對比，其中均方差用來衡量一組數(shù)據(jù)的離散程度的統(tǒng)計量。由圖11可知，在橫坐標(biāo)方向上，不同倒角的湍動能平均值差異不是很大，在縱坐標(biāo)方向上，可以看出，均方差與平均值的差異，即數(shù)據(jù)的偏離程度，也代表了湍動能的偏離程度，其差值越大表明湍流強(qiáng)度越強(qiáng)，也即流體流動狀態(tài)更不穩(wěn)定。從圖10，11曲線變化可知，倒角對流體內(nèi)部流動有很大影響，這是因為倒角可以幫助改變多級降壓結(jié)構(gòu)內(nèi)部的流體流動狀態(tài)，適當(dāng)?shù)牡菇强梢允沟昧黧w流動時更加平緩。

圖11 不同倒角湍動能平均值與均方差對比Fig.11 Comparison of mean and variance of turbulent kinetic energy with different chamfer angles

圖12示出不同凹口及出口壓比曲線。從圖中可以看出，隨著流道倒角的增加，凹口處壓降在逐漸增大，出口處壓降在逐漸減小，且流道內(nèi)的壓降在同一凹口時，倒角越大，壓降越大，說明倒角對流道壓降分布有顯著的影響。在第一凹口時，倒角從下到上的降壓順序為 0，2，4，6，8，10 mm分布，分別為21%，25%，25%，26%，26%，28%，以上可知，在第一級凹口時，壓降均沒有超過30%，同理可知，在第2和3級凹口時，均有倒角越大，壓降越大，且壓降均分布在30%附近，所以在出口處，倒角越大，壓降越小，其從下到上的排列順序為 10，8，6，2，4和 0 mm，壓比分別為 10%，13%，15%，18%，19%和20%，出口處壓降2 mm比4 mm大是因為在第3級凹口時4 mm的壓降小于2 mm壓降。當(dāng)在出口處壓降很大時，會在出口處形成高速射流和發(fā)生汽蝕現(xiàn)象，從而沖擊壁面，造成節(jié)流結(jié)構(gòu)的破壞，而適當(dāng)增加倒角則可以降低出口處分配的壓降，從而減小發(fā)生汽蝕的概率。由以上可以得出結(jié)論存在一個最佳倒角值為6 mm，可以將內(nèi)部湍流動能和湍流耗散率降到最低，同時更合理的分配壓降，減小出口處發(fā)生汽蝕的概率。

圖12 不同凹口及出口壓比曲線Fig.12 Curves of pressure ratios at different notches and outlets

3 結(jié)論

（1）多級降壓調(diào)節(jié)閥整體流量系數(shù)隨著開度的增大而增大。在80%開度以上時，流阻變化很小導(dǎo)致流量系數(shù)變化不大。開度在30%之前，流阻很大，隨著開度的增大，流阻系數(shù)迅速下降，并在中、大開度下保持較低的過流損失，由此指導(dǎo)閥門的常開開度為30%～80%之間。

（2）對于流道傾角對流體流動的影響，存在一個最佳值48°，使得流體流動時湍流動能和湍流耗散率減小到最小值，同時能減小出口處的壓力波動。

（3）渦流的產(chǎn)生總是伴隨著能量的消耗，為了降低湍流程度和湍流耗散率，最好在閥芯中倒角位置處加工1個6 mm倒角，降低出口處分配的壓降，進(jìn)而減小出口處發(fā)生汽蝕的概率。