高粘度流體截止閥開度優(yōu)化模擬研究

胡艷凱， 何備林

(1.陜西國(guó)防工業(yè)職業(yè)技術(shù)學(xué)院機(jī)械工程學(xué)院， 西安710300;2.西安理工大學(xué)機(jī)械與精密儀器工程學(xué)院， 西安710048)

引言

截止閥作為最重要的截?cái)囝愰y門之一，在航空航天領(lǐng)域發(fā)揮著重要作用[1]。得益于科學(xué)技術(shù)的發(fā)展，學(xué)者們采用實(shí)驗(yàn)和計(jì)算流體力學(xué)等方法對(duì)各型截止閥的開度因素進(jìn)行了大量的模擬研究，成果豐碩。例如：羅向龍等[2]針對(duì)某型空氣源熱泵開展了研究，得出了在相同的充灌量下除濕量和單位輸入功率除濕量隨電磁閥開度變化的曲線。王志華等[3]對(duì)R134a空氣源熱泵熱水器樣機(jī)進(jìn)行了實(shí)驗(yàn)測(cè)試，在最佳充灌量下研究了膨脹閥開度對(duì)系統(tǒng)性能的影響，結(jié)果表明增大膨脹閥的開度可以提高系統(tǒng)的COP。楊永忠等[4]基于響應(yīng)曲面法研究柴油機(jī)性能，結(jié)果表明隨著EGR閥開度的增大，柴油機(jī)動(dòng)力性下降，經(jīng)濟(jì)性惡化；當(dāng)VNT與EGR耦合時(shí)，在最大轉(zhuǎn)矩工況下和EGR閥關(guān)閉及開度較小時(shí)，隨著VNT噴嘴環(huán)開度的減小，轉(zhuǎn)矩呈現(xiàn)增大的趨勢(shì)。李智等[5]采用數(shù)值模擬方法研究粒子在有風(fēng)閥的通風(fēng)管道中的沉積規(guī)律，分析風(fēng)閥開度對(duì)不同粒徑粒子沉積效果的影響，結(jié)果表明風(fēng)閥開度越小，粒子在風(fēng)管中的沉積量越大。Qian J Y等[6]對(duì)PCGV穩(wěn)態(tài)閥芯位移進(jìn)行了實(shí)驗(yàn)和仿真比較研究，分析了彈簧剛度對(duì)流量特性和開閉響應(yīng)時(shí)間的影響，并提出了彈簧選擇的修正方程，研究有利于PCGVS或類似彈簧閥門的進(jìn)一步設(shè)計(jì)工作。文獻(xiàn)[7-15]在不同的應(yīng)用領(lǐng)域，針對(duì)各種型號(hào)的電子膨脹閥的開度進(jìn)行多因素分析，揭示了眾多流場(chǎng)規(guī)律。

但是對(duì)適用于高粘度流體的截止閥開度的研究卻并未見公開報(bào)道，所以本文基于計(jì)算流體力學(xué)的動(dòng)網(wǎng)格技術(shù)對(duì)高粘度流體截止閥的動(dòng)態(tài)關(guān)閉過程進(jìn)行模擬分析，對(duì)截止閥的設(shè)計(jì)和應(yīng)用具有重要的指導(dǎo)意義。

1 截止閥開度初步設(shè)計(jì)

高粘度流體閥門研究的目的是解決普通閥門不適用于高粘度液體的問題。一方面：由于高粘度流體本身粘度非常高，其粘性引起的摩擦阻力遠(yuǎn)大于普通流體，而市面上的閥門結(jié)構(gòu)都比較復(fù)雜，流道對(duì)流體流動(dòng)造成非常大的形狀阻力，但結(jié)構(gòu)簡(jiǎn)單的被動(dòng)單向閥又要求流體損失動(dòng)能克服密封裝置做功才能打開閥口，這個(gè)過程給本就不易流動(dòng)的高粘度流體帶來額外的流動(dòng)阻力，造成不必要的水頭損失；另一方面：市面上的閥門基本都是根據(jù)常用國(guó)標(biāo)尺寸系列或市場(chǎng)熱銷類型為標(biāo)準(zhǔn)來生產(chǎn)的，所以很難采購(gòu)到非標(biāo)的能夠應(yīng)用于本研究項(xiàng)目的合適截止閥。

考慮到上述原因，根據(jù)閥門設(shè)計(jì)手冊(cè)和閥門設(shè)計(jì)使用經(jīng)驗(yàn)，高粘度流體閥門設(shè)計(jì)采用常閉式截止閥的閥門形式，選擇推拉式電磁鐵來作為驅(qū)動(dòng)器。

為了減小高粘度流體的過流阻力，閥門設(shè)計(jì)為主動(dòng)開啟式，以盡量減小流體在閥門處的水頭損失；閥門設(shè)計(jì)著眼于能實(shí)現(xiàn)對(duì)管道的快速通斷控制而不需要調(diào)節(jié)流速流量，并且結(jié)構(gòu)足夠簡(jiǎn)單，且過流阻力小，以及維修保養(yǎng)簡(jiǎn)單等特點(diǎn)，因此選用截止閥為最優(yōu)方案。結(jié)構(gòu)簡(jiǎn)單帶來的另外一個(gè)優(yōu)點(diǎn)就是加工難度低，制造周期短。

由于高粘度流體流動(dòng)阻力大，所以將閥芯初步設(shè)計(jì)為球形，其圓輪廓有利于較小流束的阻力。文獻(xiàn)[16]提出：當(dāng)截止閥小開度工作時(shí)，流經(jīng)閥門內(nèi)壁與閥桿間的速度通常高達(dá)數(shù)十米每秒，導(dǎo)致閥門密封面或閥桿遭到損壞。同時(shí)，為獲得較小的過流阻力，需盡可能設(shè)計(jì)較大的閥門完全開度，使流束在密封線附近有充足的過流空間。綜合考慮后，將最大開度初步設(shè)計(jì)為5 mm。

截止閥結(jié)構(gòu)如圖1 所示，圖中1為HCNE型推拉式電磁鐵，2為復(fù)位彈簧，3為限位卡扣，4為螺釘，5為閥蓋，6為閥座，7為閥芯，8為動(dòng)密封壓蓋，9為O型密封圈，10為密封圈壓環(huán)，11為閥芯氟膠密封圈，其內(nèi)流場(chǎng)腔體基本尺寸見表1。

圖1 截止閥結(jié)構(gòu)圖

表1 截止閥內(nèi)流場(chǎng)腔體基本尺寸

2 動(dòng)網(wǎng)格設(shè)置

2.1 動(dòng)網(wǎng)格理論

有限元網(wǎng)格模型是將虛擬物理模型進(jìn)行有機(jī)的細(xì)化，分成有限個(gè)微小單位元，每個(gè)單位元通過共邊或者共點(diǎn)的方式相互連接，構(gòu)成整體模型。這些單位元既可以是二維的面單元，也可以是三維的體單元。求解方式是在邊界條件已知的條件下，邊界處單位元可通過插值運(yùn)算的方式求解出其他未知節(jié)點(diǎn)的運(yùn)動(dòng)參量，然后傳遞給共邊或共節(jié)點(diǎn)的相鄰單位元，依次類推可以得到由有限個(gè)運(yùn)動(dòng)參量已知的空間點(diǎn)組成的整個(gè)流場(chǎng)的運(yùn)動(dòng)場(chǎng)情況。但是當(dāng)某個(gè)模型的某個(gè)部位為運(yùn)動(dòng)的，單位元就會(huì)變形，當(dāng)單位元拓?fù)渎首銐虼髸r(shí)，通過插值計(jì)算得出的運(yùn)動(dòng)參量精度就會(huì)非常差，導(dǎo)致模型分析結(jié)果可信度低。動(dòng)網(wǎng)格技術(shù)由此應(yīng)運(yùn)而生，用戶通過該技術(shù)可合理設(shè)置網(wǎng)格拓?fù)渎室灾付ňW(wǎng)格運(yùn)動(dòng)范圍。當(dāng)模型運(yùn)動(dòng)時(shí)，如果網(wǎng)格變形超出這個(gè)范圍，那么這個(gè)局部網(wǎng)格就會(huì)重新被強(qiáng)制調(diào)整，有可能是進(jìn)行局部光滑處理，也有可能是局部增加或較少，抑或是局部重新劃分。經(jīng)過調(diào)整后的網(wǎng)格會(huì)重新獲得很高的質(zhì)量，可計(jì)算出足夠高精度的運(yùn)動(dòng)參量值。

有限元中動(dòng)網(wǎng)格技術(shù)的具體實(shí)施是通過合理指定模型的運(yùn)動(dòng)件和指定運(yùn)動(dòng)方式相配合來完成的。指定運(yùn)動(dòng)方式可以通過Fluent軟件自帶的邊界函數(shù)和調(diào)用外部自定義程序(UDF)相結(jié)合來實(shí)現(xiàn)；由于Fluent軟件自帶的邊界函數(shù)操作不便，函數(shù)類型有限，所以本文采用調(diào)用外部自定義程序的方式來控制運(yùn)動(dòng)件的運(yùn)動(dòng)方式。

由于Fluent軟件計(jì)算原則嚴(yán)格遵循連續(xù)性方程，所以模擬計(jì)算的動(dòng)網(wǎng)格不可以運(yùn)動(dòng)到5 mm，否則會(huì)出現(xiàn)動(dòng)網(wǎng)格和靜網(wǎng)格的直接接觸，造成邊界類型出現(xiàn)邏輯錯(cuò)誤，導(dǎo)致網(wǎng)格重構(gòu)失敗。故設(shè)置模擬最小開度為0.01 mm，即閥芯從完全打開的位置向密封圈移動(dòng)4.9 mm。

2.2 閥芯運(yùn)動(dòng)程序設(shè)計(jì)

由于Fluent軟件的底層程序由C語言構(gòu)成，所以對(duì)于外部自定義程序也需要使用C語言來進(jìn)行編程。UDF程序的主要部分命令流如下所示：

DEFINE_CG_MOTION(valve, dt, vel, omega, time, dtime)%調(diào)用宏命令

{

real NV_VEC (A);%指定閥芯動(dòng)網(wǎng)格面

FILE *f_frc;

{

if(loc_prev>=4.9)%循環(huán)語句判斷條件

v_prev= -1;%指定速度值

else

return;

}

loc_prev=loc_prev+v_prev*dtime;%閥芯路徑更新

vel[2] =v_prev;%閥芯速度賦值

}

3 截止閥關(guān)閉過程模擬分析

3.1 時(shí)間步長(zhǎng)敏感性分析

對(duì)于瞬態(tài)分析而言，時(shí)間步長(zhǎng)的設(shè)置是一個(gè)矛盾體，即時(shí)間步長(zhǎng)與計(jì)算步數(shù)成反比。設(shè)置足夠小的時(shí)間步長(zhǎng)有利于反映流場(chǎng)中非常細(xì)微-快速的變化，但同時(shí)使計(jì)算時(shí)間呈正比增加，導(dǎo)致效率成倍下降；反之，設(shè)置較大的時(shí)間步長(zhǎng)，就有可能忽略了微小的流場(chǎng)渦旋作用，使流場(chǎng)的模擬準(zhǔn)確度大打折扣，導(dǎo)致計(jì)算結(jié)果精度下降。因此，對(duì)時(shí)間步長(zhǎng)敏感性的檢驗(yàn)分析變得非常必要。

參考相關(guān)文獻(xiàn)設(shè)置經(jīng)驗(yàn)，本節(jié)分別設(shè)置時(shí)間步長(zhǎng)為0.01 s、0.001 s、0.0001 s、0.000 01 s和0.000 001 s，控制其他變量不變來進(jìn)行模擬對(duì)比分析。不同時(shí)間步長(zhǎng)的截止閥壓降如圖2 所示。

圖2 不同時(shí)間步長(zhǎng)的截止閥壓降

從圖2 可知，迭代時(shí)間步長(zhǎng)從0.01 s減小到0.0001 s時(shí)，壓降折線的斜率比較大，波動(dòng)比較明顯；當(dāng)時(shí)間步長(zhǎng)小于1e-5后，壓降變化平緩。其原因是：當(dāng)時(shí)間步長(zhǎng)比較大時(shí)，液體從進(jìn)口流入層流遇到形狀阻礙后，由于局部單位元迭代計(jì)算的細(xì)分量不夠小，造成計(jì)算值描述的湍流還未經(jīng)充分發(fā)展就被傳遞給下游的單位元，導(dǎo)致模擬計(jì)算失真。流層之間的湍動(dòng)不夠強(qiáng)烈，局部流層混亂程度小使機(jī)械能損失小，導(dǎo)制截止閥在這種步長(zhǎng)下的壓降小。但隨著時(shí)間步長(zhǎng)的減小，迭代計(jì)算捕捉到更小尺度的漩渦和局部二次流，此時(shí)模擬對(duì)于湍流的描述相對(duì)前者而言有一定程度改善，造成壓降增大；當(dāng)時(shí)間步長(zhǎng)足夠細(xì)小的時(shí)候，迭代計(jì)算能夠?qū)⒔^大部分的湍流都捕捉到；隨著時(shí)間步長(zhǎng)繼續(xù)減小，剩下的還沒有被描述的局部漩渦也越來越少；同時(shí)，過小的漩渦對(duì)整體流場(chǎng)的影響可以忽略不計(jì)，所以此時(shí)壓降變化并不明顯。

另一方面，在模擬計(jì)算時(shí)，隨著時(shí)間步長(zhǎng)的減小，整個(gè)模擬計(jì)算時(shí)間和計(jì)算機(jī)的物理計(jì)算內(nèi)存占用率都成倍增加。所以考慮到對(duì)計(jì)算精度和計(jì)算效率的綜合影響，本文選取時(shí)間步長(zhǎng)為0.000 01 s。

3.2 模擬分析邊界條件設(shè)置

高粘度流體屬于微可壓流體，故采用壓力求解器，并激活瞬態(tài)分析模式；模型設(shè)置為RNG k-e；設(shè)置流體粘度為6 Pa·s，密度為1160 kg/m3；

為避免恒定流量進(jìn)口邊界條件掩蓋了流量變化的反映，選擇入口條件為pressure-inlet，其中大小設(shè)置為500 000 Pa，標(biāo)量設(shè)置Turbuluent Intensity為0.05、Hydraulic Diameter為0.01 m；出口壓力設(shè)置為30 000 Pa，標(biāo)量同進(jìn)口。

在define中導(dǎo)入U(xiǎn)DF程序，并激活動(dòng)網(wǎng)格中smoothing、layering和remeshing，設(shè)置spring constant favtor為0.4，minimum length scale為0.2 mm，maximum length scale為1 mm，maximum cell skewness為0.7；動(dòng)網(wǎng)格域設(shè)置為閥芯面，type設(shè)置為rigid body，在meshing options里設(shè)置cell height為0；設(shè)置時(shí)間步長(zhǎng)為0.000 01 s，每50步記錄一次參數(shù)。

3.3 截止閥開度模擬分析

為充分了解截止閥在關(guān)閉過程中，閥芯不同高度對(duì)閥門流量以及過流性能的影響，定義閥芯以1 m/s的速度勻速關(guān)閉進(jìn)行模擬，每間隔0.5 mm輸出一個(gè)瞬態(tài)相對(duì)壓強(qiáng)圖。閥門不同開度下的模擬壓強(qiáng)如圖3 所示。

圖3 閥門不同開度下的模擬壓強(qiáng)圖

從圖3 可知，隨著閥門開度的降低，瞬態(tài)模擬壓強(qiáng)圖的極值呈現(xiàn)近似指數(shù)式地增長(zhǎng)。正是由于極差非常大，所以出圖時(shí)沒有將左側(cè)的刻度標(biāo)尺進(jìn)行統(tǒng)一規(guī)劃。圖3 (a)至圖3 (d)的壓場(chǎng)分布比較規(guī)則，壓強(qiáng)由進(jìn)口沿流線向出口遞減，橫向梯度均勻，沒有出現(xiàn)局部渦流，說明該閥門開度沒有對(duì)高粘度流體的過流狀態(tài)造成大的影響。圖3 (e)至圖3 (j)中均出現(xiàn)空白區(qū)域，即壓力場(chǎng)中出現(xiàn)沿流線方向相反的壓強(qiáng)梯度，從圖中閥芯的左上角可以看出很高的局部壓強(qiáng)，甚至高于進(jìn)口壓強(qiáng)，達(dá)到了刻度標(biāo)尺上限之上的高度，且能看見越來越明顯的不規(guī)則漩渦。

為進(jìn)一步說明閥芯上方區(qū)域的湍動(dòng)情況，下面給出內(nèi)流場(chǎng)速度矢量圖，如圖4 所示。

圖4 內(nèi)流場(chǎng)速度矢量圖

從圖4 (a)～圖4 (d)可知，閥芯位置在1.5 mm和0.5 mm時(shí)，閥芯側(cè)方的速度都出現(xiàn)了明顯的逆流線方向矢量，且后者更加明顯；同時(shí)，逆流矢量長(zhǎng)度大于周圍順流線速度矢量；閥芯位置在0.5 mm時(shí)局部湍動(dòng)速度，大于閥芯位置為1.5 mm時(shí)的速度，該變化趨勢(shì)與壓強(qiáng)的變化趨勢(shì)(圖3 )一致，由此說明閥芯兩側(cè)的局部湍動(dòng)異常劇烈。

產(chǎn)生上述結(jié)果的原因是過流面積大幅減小，進(jìn)口的壓強(qiáng)累積到閥口處產(chǎn)生局部射流，引起局部壓場(chǎng)重新分布，即形狀阻力導(dǎo)致邊界層分離引發(fā)強(qiáng)烈的局部湍流。文獻(xiàn)[17]也出現(xiàn)了與本文類似的現(xiàn)象：閥芯出口側(cè)和閥體出口段下游位置出現(xiàn)了低壓區(qū)域，閥芯開啟后形成的流動(dòng)間隙是壓降產(chǎn)生的主要原因。閥芯在閥桿的帶動(dòng)下運(yùn)動(dòng)，在結(jié)構(gòu)力學(xué)上類似于懸臂梁結(jié)構(gòu)，局部湍動(dòng)會(huì)引發(fā)閥芯的喘振，這也是造成噪聲和氣蝕的重要原因之一，同時(shí)也造成懸臂梁結(jié)構(gòu)左右受力不均勻，長(zhǎng)期工作的情況下非常不利于閥桿上方的動(dòng)密封，并嚴(yán)重影響閥芯的壽命和密封效果。

局部高強(qiáng)度湍動(dòng)發(fā)生在閥芯上方，而高速射流發(fā)生在密封線附近，這個(gè)現(xiàn)象非常符合流體動(dòng)力學(xué)原理，即：湍動(dòng)的發(fā)生是由于剪切壓降引起的，而跟該區(qū)域流速大小無關(guān)。

閥芯運(yùn)動(dòng)到4.9 mm位置時(shí)，截止閥內(nèi)部流場(chǎng)的瞬態(tài)波動(dòng)壓強(qiáng)高達(dá)2.12 MPa(圖3 (j))。這個(gè)瞬態(tài)高壓會(huì)對(duì)閥芯和閥體施加一個(gè)沖擊力，嚴(yán)重時(shí)會(huì)損壞系統(tǒng)結(jié)構(gòu)造成經(jīng)濟(jì)損失和員工人身安全傷害，所以建議在工程中高速高壓的工況下，將截止閥與安全閥或泄壓閥配合使用，以削弱流體瞬態(tài)沖擊的影響來保證系統(tǒng)安全性。

瞬態(tài)模擬輸出的閥芯運(yùn)動(dòng)面相對(duì)閥門進(jìn)口邊界的阻力系數(shù)曲線如圖5 所示。

圖5 瞬態(tài)模擬阻力系數(shù)曲線

從圖5 可知，從瞬態(tài)流動(dòng)時(shí)間0.0025 s起，阻力系數(shù)開始出現(xiàn)脈動(dòng)現(xiàn)象，且越來越劇烈，說明流場(chǎng)出現(xiàn)明顯的湍流波動(dòng)，說明這幾種開度對(duì)高粘度流體的過流造成很大的影響，且隨著開度的降低，阻礙作用越明顯。

進(jìn)口段中間壓強(qiáng)明顯比進(jìn)口邊界和閥內(nèi)腔體都高一個(gè)色級(jí)(圖3 (i))，說明此時(shí)已經(jīng)產(chǎn)生“水錘”作用；由于閥門開度很低，由進(jìn)口流入的高粘度液體流到閥芯處時(shí)受到嚴(yán)重阻礙而無法順利向下游繼續(xù)流動(dòng)，但因慣性作用流體并不能立即停止下來，源源不斷的來流就會(huì)向前撞擊和擠壓，給閥門的阻礙結(jié)構(gòu)一個(gè)液動(dòng)錘擊力；這個(gè)力并沒有被閥門吸收，因此動(dòng)能聚集成壓能后只能向其他壓強(qiáng)較低的方向傾泄，即反流向方向；壓能形成一股壓力波沿反流向方向傳遞，圖3 (i)中正好處于這個(gè)反流向傳遞狀態(tài)，圖5 中也反映出了這一變化，即圖中的最大波峰。由于圖3 (i)中捕捉的是閥門0.5 mm開度的時(shí)刻，對(duì)應(yīng)瞬態(tài)流動(dòng)的0.0045 s；而圖5 中最大波峰出現(xiàn)在瞬態(tài)流動(dòng)時(shí)間的0.006 s處，出現(xiàn)時(shí)間差的原因有二：(1)由于流體慣性引起的壓能聚集在0.0045 s之前已經(jīng)開始，同相的反向壓力波傳遞早在異相的閥結(jié)構(gòu)發(fā)生泊松耦合變化之前就已發(fā)生；(2)流體液動(dòng)壓力引起閥門固體結(jié)構(gòu)行程流固耦合作用的過程也需要傳遞時(shí)間，并不是絕對(duì)地同時(shí)發(fā)生，ANSYS軟件的迭代計(jì)算方式在本質(zhì)上很好地體現(xiàn)了這一微小時(shí)間差。

圖6 密封圈截面質(zhì)量流量率圖

密封圈處截面的質(zhì)量流量率監(jiān)控圖如圖6 所示。

從圖6 可知，閥門開度剛開始下降時(shí)，對(duì)高粘度流體的質(zhì)量流量率基本沒有影響，但是當(dāng)開度下降到三分之二附近的位置后，閥門的質(zhì)量流量出現(xiàn)非常明顯的上升，這與圖3 反映的變化趨勢(shì)高度一致。

出現(xiàn)這種現(xiàn)象的原因是由于閥芯在關(guān)閉過程中，其輪廓下降到密封吼線以下時(shí)使密封線處的有效過流面積減小，但是進(jìn)口處壓強(qiáng)保持恒定，所以閥芯處的壓能聚集使減小的有效過流面積內(nèi)的過流速度明顯提高，導(dǎo)致質(zhì)量流量率上升。

閥門不同開度下的壓強(qiáng)圖如圖7 所示。

圖7 閥門不同開度下的壓降圖

從圖7 可知，在0到0.0015 s的時(shí)間步內(nèi)，壓降曲線比較平緩，增長(zhǎng)幅度很小；但在0.0015 s以后，每段0.005 s間隔內(nèi)的曲線斜率都比前段的斜率要高，說明流阻從0.0015 s以后開始越來越大。

綜上所述，截止閥的最大開度優(yōu)化為3.5 mm時(shí)更加合適。一方面，3.5 mm的開度對(duì)于高粘度流體的過流基本不造成阻礙作用，并保留一定的安全余量；另一方面，減小閥門最大開度就是減小驅(qū)動(dòng)件電磁鐵的最大行程。最大開度為3.5 mm時(shí)，在保證足夠驅(qū)動(dòng)力的前提下可大幅減小電磁鐵的尺寸和耗用功率，科學(xué)合理地降低了制造成本，并能獲得更小巧的外形尺寸。

3.4 閥芯形狀優(yōu)化設(shè)計(jì)分析

通過模擬分析發(fā)現(xiàn)閥芯在趨于關(guān)閉時(shí)，閥芯上方出現(xiàn)了非常強(qiáng)烈的湍流，影響截止閥的工作穩(wěn)定性。通過在球形閥芯外輪廓左右設(shè)置兩個(gè)監(jiān)測(cè)點(diǎn)，提取動(dòng)壓強(qiáng)，并繪制變化曲線，如圖8 所示。

圖8 閥芯左右側(cè)監(jiān)測(cè)點(diǎn)的動(dòng)壓曲線

從圖8 可知，隨著流場(chǎng)時(shí)間的增加，左右兩點(diǎn)的動(dòng)壓差距急劇擴(kuò)大，極差可達(dá)40 MPa。這會(huì)給閥芯一個(gè)很大的橫向壓力使閥芯偏離豎直方向而跑偏，引起截止閥關(guān)閉效果不佳等問題。

故對(duì)閥芯的形狀進(jìn)行優(yōu)化設(shè)計(jì)以改善這種局部湍流狀況，閥芯形狀優(yōu)化前后對(duì)比如圖9 所示。

圖9 閥芯優(yōu)化前后對(duì)比圖

為了削弱圓柱繞流的形狀阻力影響，將閥芯上部設(shè)計(jì)成流線型肋板的形式，使原來產(chǎn)生局部對(duì)稱湍流的空間被閥芯實(shí)體代替，同時(shí)，將截止閥的最大開度設(shè)計(jì)為3.5 mm，由此得到優(yōu)化后的閥芯形狀如圖9 (b)所示。其他條件不變，對(duì)優(yōu)化后的閥芯進(jìn)行模擬，每0.5 mm輸出一次壓強(qiáng)，如圖1 0所示。

圖1 0閥芯優(yōu)化后不同開度模擬壓強(qiáng)圖

從圖1 0可知，優(yōu)化后的閥芯上方?jīng)]有出現(xiàn)局部湍流，說明閥芯優(yōu)化效果顯著。但是圖1 0(f)中的壓強(qiáng)明顯低于圖1 0(e)和圖1 0(g)，說明水錘作用依然發(fā)生。

圖1 1閥芯優(yōu)化后瞬態(tài)模擬阻力系數(shù)曲線

閥芯優(yōu)化后瞬態(tài)模擬阻力系數(shù)曲線如圖1 1所示。

從圖1 1可知，優(yōu)化后的閥芯在關(guān)閉后期出現(xiàn)壓力脈動(dòng)現(xiàn)象，但數(shù)值等級(jí)得到了大幅減小。雖然閥芯附近的局部湍流問題得到改善且閥芯升力系數(shù)大幅減小，但脈動(dòng)壓力影響仍然存在，所以后續(xù)還需要繼續(xù)從流固耦合的角度進(jìn)行進(jìn)一步的研究。

4 結(jié)論

(1)通過設(shè)置不同的時(shí)間步長(zhǎng)分別模擬分析了瞬態(tài)流場(chǎng)，以壓降的穩(wěn)定性為判斷依據(jù)得出最佳的時(shí)間步長(zhǎng)。

(2)通過模擬分析揭示了由截止閥關(guān)閉引發(fā)的水錘效應(yīng)暫態(tài)失穩(wěn)現(xiàn)象，說明過流液產(chǎn)生的交變液動(dòng)沖擊力不利于截止閥的的結(jié)構(gòu)穩(wěn)定性。

(3)以流阻為判斷依據(jù)，通過分析閥芯運(yùn)動(dòng)過程中的壓降和流量等流場(chǎng)參量的變化，得到了截止閥芯最大開度的最佳值。

(4)通過在截止閥芯上方增加流線形筋的方式來削弱局部湍流的影響，并通過模擬得出優(yōu)化后的閥芯能有效地改善內(nèi)流場(chǎng)的湍流強(qiáng)度。