，，，

(1.大連海事大學(xué) 輪機(jī)工程學(xué)院，遼寧 大連 116026；2．中國(guó)船舶工業(yè)集團(tuán)公司 船舶系統(tǒng)工程部，北京 100036)

傳統(tǒng)船舶空氣動(dòng)力學(xué)研究大多采用風(fēng)洞實(shí)測(cè)的方法確定流場(chǎng)速度和壓力分布。該方法造價(jià)高，耗資大并且在風(fēng)洞試驗(yàn)中存在著動(dòng)力相似和幾何相似的影響，還要考慮風(fēng)洞邊界條件的影響和湍流、風(fēng)速、風(fēng)向、雷諾數(shù)等。實(shí)驗(yàn)結(jié)果要進(jìn)行換算，還存在測(cè)量點(diǎn)有限、采集數(shù)據(jù)的測(cè)量誤差等問(wèn)題，對(duì)于紊亂的流動(dòng)區(qū)域難以獲得準(zhǔn)確的結(jié)果[1-2]。利用數(shù)值模擬的方法研究船舶內(nèi)外流場(chǎng)的流動(dòng)狀況，對(duì)流場(chǎng)參數(shù)進(jìn)行分析，并與傳統(tǒng)的研究方法結(jié)合，可以有效地改善船舶性能、節(jié)約研究資金，提高研究效率。

目前，國(guó)內(nèi)開(kāi)展的進(jìn)氣室內(nèi)流場(chǎng)研究中，針對(duì)進(jìn)氣濾清裝置和蝸殼內(nèi)的流場(chǎng)結(jié)構(gòu)及其優(yōu)化的研究相對(duì)較多[3-4]，但是考慮進(jìn)氣室整體結(jié)構(gòu)以及消音器、濾清器等氣動(dòng)阻力影響的研究很少。同時(shí)，航行中的船體由于其外流場(chǎng)所處的環(huán)境的變化、甲板上的各種設(shè)施使空氣流場(chǎng)中包含了流體分離、脫落渦、循環(huán)等諸多復(fù)雜現(xiàn)象[5-6]。這些流動(dòng)的復(fù)雜性直接導(dǎo)致了進(jìn)氣系統(tǒng)內(nèi)部流動(dòng)的不確定性，其最大的影響就是動(dòng)力系統(tǒng)進(jìn)氣流量的變化。在某些條件下會(huì)引起進(jìn)氣流量不足致使動(dòng)力系統(tǒng)無(wú)法正常工作。進(jìn)氣系統(tǒng)進(jìn)氣流量的變化通過(guò)濾清器前后壓差的變化得以體現(xiàn)，因此往往要在濾清器前后安裝壓力探測(cè)裝置，當(dāng)壓差大于一定值就會(huì)報(bào)警，以便及時(shí)采取補(bǔ)救措施，增加進(jìn)氣流量。但是在實(shí)際工作中，濾清器前后壓差還沒(méi)有達(dá)到臨界值，報(bào)警器就會(huì)報(bào)警，給實(shí)際工作帶來(lái)不必要的麻煩。因此研究濾清器前后壓差隨進(jìn)口流動(dòng)參數(shù)的變化規(guī)律，并探求壓差變化與測(cè)點(diǎn)位置間的關(guān)系是十分必要的。

本文用定常數(shù)值模擬的方法研究額定進(jìn)氣條件下進(jìn)氣系統(tǒng)的總體性能；在文獻(xiàn)[6]的基礎(chǔ)上，以周期函數(shù)的形式模擬外流場(chǎng)的變化，對(duì)進(jìn)氣室內(nèi)部流場(chǎng)進(jìn)行非定常研究。重點(diǎn)研究進(jìn)口氣流條件變化對(duì)濾清器前后壓差的影響，為壓差探測(cè)裝置的安裝位置以及敏感度的設(shè)置提出建議。

1 計(jì)算模型和邊界條件

1.1 計(jì)算模型

本文研究的某船舶進(jìn)氣系統(tǒng)由前進(jìn)氣室、后進(jìn)氣室、消音器、豎井和進(jìn)氣蝸殼等結(jié)構(gòu)組成。計(jì)算模型由Gambit生成，見(jiàn)圖1。其中百葉窗、濾清器和豎井處采用結(jié)構(gòu)化網(wǎng)格，而對(duì)于結(jié)構(gòu)較復(fù)雜的消音器以及前后進(jìn)氣室則采用非結(jié)構(gòu)化網(wǎng)格。計(jì)算域整體網(wǎng)格數(shù)約為159.6萬(wàn)。

1.2 邊界條件

1.2.1 定常計(jì)算邊界條件

氣體流經(jīng)進(jìn)氣系統(tǒng)的過(guò)程中，受百葉窗、濾清器等結(jié)構(gòu)的阻力影響很大，造成較大的壓力損失。因此，為了盡可能真實(shí)地反映進(jìn)氣系統(tǒng)的性能，這些部件的阻力是不能被忽略的。但是在數(shù)值模擬中，真實(shí)地再現(xiàn)這些部件的具體結(jié)構(gòu)，比如慣性級(jí)的結(jié)構(gòu)，勢(shì)必會(huì)使計(jì)算網(wǎng)格大量增加，增加計(jì)算難度和周期。因此，本研究采用Fluent軟件包內(nèi)的Fan邊界條件來(lái)直接給出這些結(jié)構(gòu)后的壓力損失，從而使計(jì)算簡(jiǎn)化。

圖1 進(jìn)氣系統(tǒng)模型

Fan邊界是Fluent軟件包中的一個(gè)模型。該模型假設(shè)在流場(chǎng)中存在一個(gè)無(wú)限小的風(fēng)扇(見(jiàn)圖2)，對(duì)其后氣流產(chǎn)生作用，使其后氣流存在一定的壓升。這樣一種壓力的不連續(xù)可以用與速度相關(guān)的方程表示，該方程可以是常量，也可以是多項(xiàng)式、分段函數(shù)或者自定義的某種函數(shù)。本文通過(guò)給定一個(gè)常量來(lái)表示相應(yīng)的結(jié)構(gòu)阻力。

圖2 Fan模型作用示意

具體邊界條件設(shè)置見(jiàn)圖3。百葉窗進(jìn)口給定質(zhì)量進(jìn)口，進(jìn)口流量為燃?xì)廨啓C(jī)額定工作流量(本文中所涉及的流量、壓力均以與額定流量的比值形式表示)。進(jìn)氣蝸殼出口給定壓力出口。百葉窗出口、濾清器出口以及進(jìn)氣蝸殼進(jìn)口給定Fan邊界條件。消音器出口和豎井進(jìn)口交界按混和平面法處理。其它固體邊界作壁面處理。給定濾清器壁面粗糙度系數(shù)1，消音器隔板給定粗糙度系數(shù)0.5，粗糙度高度0.001 m。計(jì)算采用二階迎風(fēng)格式，湍流模型選取標(biāo)準(zhǔn)的k-ε模型，其中k取0.8，湍流耗散率ε取0.8。

圖3 邊界條件示意

1.2.2 非定常計(jì)算進(jìn)口條件

由文獻(xiàn)[6]知，風(fēng)速恒定時(shí)，進(jìn)氣室百葉窗前的總壓變化與船體和風(fēng)向的夾角近似成正弦關(guān)系。本文在研究進(jìn)氣室內(nèi)部流場(chǎng)隨外界條件變化時(shí)做如下簡(jiǎn)化：假定航向、航速和風(fēng)速恒定，而風(fēng)向以一定頻率變化。在這種條件下求解內(nèi)部流場(chǎng)參數(shù)隨進(jìn)口總壓的變化情況，假設(shè)進(jìn)口總壓波動(dòng)按正弦形式變化，即

p*=101 450+Asin(ωt)

(1)

式中:ω——風(fēng)向變化的角速度；

t——時(shí)間步長(zhǎng)；

A——總壓變化的幅值。

研究中假定船的航向和航速不變，風(fēng)向改變可以使進(jìn)口總壓分布涵蓋船在不同風(fēng)向時(shí)的大多數(shù)進(jìn)口總壓值，而風(fēng)向變化的角速度越大也就意味著風(fēng)向改變得越快?？紤]到本文研究的重點(diǎn)是濾清器前后壓差隨進(jìn)口條件變化的情況，而實(shí)際中濾清器前后允許的最大壓差為250 Pa，因而確定進(jìn)口總壓波動(dòng)曲線中A為250 Pa，同時(shí)，給定為0.497 1 rad/s，也就是風(fēng)向轉(zhuǎn)動(dòng)1周的周期為12.63 s，則進(jìn)口壓力曲線見(jiàn)圖4。

圖4 百葉窗進(jìn)口總壓波動(dòng)曲線

非定常計(jì)算采用雙時(shí)間步法，對(duì)出口流量進(jìn)行監(jiān)測(cè)，當(dāng)流量呈現(xiàn)規(guī)則的周期性變化時(shí)認(rèn)為計(jì)算進(jìn)入收斂階段，開(kāi)始對(duì)計(jì)算結(jié)果進(jìn)行保存。

2 計(jì)算結(jié)果分析

2.1 定常計(jì)算結(jié)果分析

為了能夠定量分析各位置的壓力變化，選定如圖5所示截面。各截面與對(duì)應(yīng)的百葉窗出口、濾清器進(jìn)出口以及進(jìn)氣蝸殼進(jìn)口均相距5 mm。各截面位置的總壓與給定值的比較見(jiàn)表1。可以看出，百葉窗和濾清器前后的壓力損失與給定值很接近，只有蝸殼前后的壓力損失與給定值差別較大。這是因?yàn)槲仛ぬ帤饬靼l(fā)生90°折轉(zhuǎn)，流場(chǎng)的不均勻性必然會(huì)增加，從而增加流動(dòng)損失?？梢酝茢啵谖仛じ浇鲃?dòng)非常復(fù)雜，損失也很大。從其它位置壓力損失的計(jì)算結(jié)果來(lái)看，文中采用的Fan邊界條件基本可以準(zhǔn)確地模擬壓力損失。

圖5 描述阻力損失的截面

表1 各裝置壓力損失

由百葉窗進(jìn)口至進(jìn)氣蝸殼出口的氣流壓差曲線見(jiàn)圖6。

圖6 壓差曲線

從圖6可見(jiàn)，壓差的變化有三次階越。第一次壓力損失最大，即濾清器慣性級(jí)的影響；第二次損失相對(duì)最小，是后進(jìn)氣室氣流轉(zhuǎn)向進(jìn)入消音器所引起的；第三次壓力損失主要是在蝸殼部分產(chǎn)生的，除了空氣阻力引起的損失外，絕大部分損失是由于氣流變向，流場(chǎng)不均勻引起的。

2.2 進(jìn)氣系統(tǒng)非定常研究

為防止進(jìn)氣室濾清器前后壓差過(guò)大，特在濾清器前后布置總壓測(cè)點(diǎn)，當(dāng)測(cè)點(diǎn)壓差大于臨界值時(shí)，報(bào)警器就會(huì)報(bào)警，從而開(kāi)啟應(yīng)急裝置增加濾清器后總壓。為了系統(tǒng)研究壓差與測(cè)點(diǎn)位置間的關(guān)系，選取如圖7所示的G1～G7點(diǎn)作為總壓差的監(jiān)測(cè)點(diǎn)。其中，G2點(diǎn)為船體實(shí)際工作中壓差探測(cè)器的布置位置，G5點(diǎn)在G2點(diǎn)正上方，G1、G3、G4、G6分別位于濾清器窗口的四角處，G7位于小濾清器窗口底邊的中央位置。濾清器后部測(cè)點(diǎn)與前測(cè)點(diǎn)位置一一對(duì)應(yīng)。計(jì)算中求得各個(gè)時(shí)刻的濾清器前后測(cè)點(diǎn)的總壓值，然后將各點(diǎn)對(duì)應(yīng)時(shí)刻的總壓值相減得到各測(cè)點(diǎn)的壓差隨時(shí)間變化曲線。

圖7 濾清器前后壓差測(cè)點(diǎn)布置

為方便表述，將濾清器前總壓用“front”表示，濾清器后總壓用“back”表示。不同時(shí)刻各測(cè)點(diǎn)濾清器前后的總壓見(jiàn)圖8。顯然，測(cè)點(diǎn)前后的總壓分布與百葉窗進(jìn)口相比存在一定的相位偏移，但測(cè)點(diǎn)前的相位偏移很小，而測(cè)點(diǎn)后偏移較大。

從測(cè)點(diǎn)前后總壓分布曲線可以看出，測(cè)點(diǎn)位置不同，測(cè)點(diǎn)前后總壓的分布存在很大差異。在靠近進(jìn)氣室的底部各測(cè)點(diǎn)處，濾清器前總壓普遍高于濾清器后的總壓，尤以G3測(cè)點(diǎn)處最為明顯；而中間高度處的測(cè)點(diǎn)(G4～G7)總壓分布則剛好相反，濾清器后的總壓要高一些?？梢猿醪酱_定，在G3測(cè)點(diǎn)處獲得壓差全部為負(fù)值，也就是總壓有損失；而在G7測(cè)點(diǎn)處獲得總壓則全部為正值，說(shuō)明此處總壓升高。另外，除總壓損失或總壓升高會(huì)產(chǎn)生壓差外，測(cè)點(diǎn)前后總壓變化的不同步同樣是造成壓差的原因。

圖8 燃燒空氣濾清器前后總壓分布曲線

3 結(jié)論

1)Fan邊界條件基本可以準(zhǔn)確地模擬各關(guān)鍵部件前后的壓力損失，加上研究中對(duì)進(jìn)氣系統(tǒng)的完整建模，從而使流場(chǎng)計(jì)算結(jié)果更接近真實(shí)情況。氣流流經(jīng)進(jìn)氣系統(tǒng)有三次明顯的壓力損失，第一次壓力損失最大，即濾清器慣性級(jí)的影響，第二次和第三次壓力損失主要由氣流轉(zhuǎn)向引起。

2)濾清器前后產(chǎn)生壓差的原因有兩個(gè):①氣體經(jīng)過(guò)濾清器后壓力損失；②氣體經(jīng)過(guò)濾清器后壓力損失很小，但是相位不同造成壓差。不同測(cè)點(diǎn)處得到的壓差曲線差別很大，存在完全正壓差或完全負(fù)壓差的情況

3)基于測(cè)量濾清器前后壓差的目的，測(cè)點(diǎn)不應(yīng)放置在G7附近。G2位置正負(fù)壓差變化較平均，是比較理想的觀測(cè)點(diǎn)。

[1] 姚育成,李萬(wàn)平,李良軍.高雷諾數(shù)情況下鈍體繞流的數(shù)值模擬[J].華中科技大學(xué)學(xué)報(bào)，2003,31(2):106-108.

[2] XIENAN,DODWORTH K,VASSALOS D. Computation of free surface turbulent flow around a wigley hull [J]. Journal of Ship Mechanics，2001,5(6):1-8.

[3] 張?zhí)煲?王傳溥.船用燃?xì)廨啓C(jī)進(jìn)氣過(guò)濾系統(tǒng)的計(jì)算流體力學(xué)(CFD)研究[J].哈爾濱工程大學(xué)學(xué)報(bào),2000,21(1):20-24.

[4] 陳 康,黃德波.CFD技術(shù)在三體船阻力性能研究中的應(yīng)用[J].哈爾濱工程大學(xué)學(xué)報(bào)，2006,27(3):362-366.

[5] 吳裕平.CFD在艦船甲板流場(chǎng)計(jì)算與特性研究中的應(yīng)用[J].直升機(jī)技術(shù),2006,147(3):1-4.

[6] 于 洋,孫 鵬,鐘兢軍.風(fēng)向?qū)?cè)向進(jìn)氣裝置外流場(chǎng)影響的數(shù)值研究[J].工程熱物理學(xué)報(bào)，2009,30(3):389-392.