王三武，夏德印，王華奧，鄒琳，金華勝

(武漢理工大學(xué)機電工程學(xué)院，湖北 武漢 430070)

氣動汽車是一種新型無污染綠色汽車。法國的MDI公司已經(jīng)研制出自行設(shè)計的以壓縮空氣為動力的氣動汽車樣機，Lin LIU，Xiaoli YU等[1]也探究了氣動汽車的可行性，均說明了氣動汽車的可行性及其優(yōu)勢。但不可否認(rèn)的是氣動汽車還有很多自身的缺點需要克服，如氣動汽車減壓過程氣體可用能損失過大是制約氣動汽車發(fā)展的重要因素之一。

很多學(xué)者對高壓氣動減壓裝置進行了設(shè)計、改進及能耗分析等研究。何為、吳玉庭等[2]采用有效能分析法對如何降低氣動汽車減壓過程的能耗等進行了分析，得到氣動汽車各種能效損失中減壓裝置引起的可用能損耗占有的比例最大的結(jié)論。張遠(yuǎn)深、宋有明等[3]對氫能源動力汽車的高壓氣動減壓閥進行了結(jié)構(gòu)優(yōu)化及特性分析，通過對減壓閥的各個主要部件進行分析，優(yōu)化了該閥的結(jié)構(gòu)，設(shè)計了一種動態(tài)性能良好的新型高壓氣動減壓閥。賈光政等[4]研制了應(yīng)用于氣動汽車的高壓氣體容積減壓裝置，介紹了該裝置的組成和工作原理并進行了實驗研究，證明了該高壓氣動容積減壓裝置具有良好的減壓、穩(wěn)壓效果。丁衛(wèi)華、劉昊等[5]應(yīng)用可用能分析方法,對氣動汽車動力系統(tǒng)主要環(huán)節(jié)的可用能消耗進行了分析，仿真計算結(jié)果表明，節(jié)流減壓可用能損耗約占各環(huán)節(jié)能耗值的50%。周杰、王建華[6]對氣體從高壓氣瓶到做功排出的整個過程進行了研究。訚耀保、陳潔萍等[7]分析了車載超高壓氣動減壓閥的工作機理、基本特性及其影響因素。訚耀保、沈力等[8]采用計算流體動力學(xué)方法并利用Gambit及Fluent軟件對某氫能源汽車兩級高壓氣動減壓閥進行了流場分析。

國內(nèi)外很多學(xué)者對氣動汽車減壓裝置做功過程中的能耗和流場進行了分析，但是針對減壓裝置能耗過大和流場穩(wěn)定性不高等方面的改進研究較少。本研究在前人研究的基礎(chǔ)上進行了更加全面深入的研究，對目前氣動汽車上比較常用的節(jié)流減壓裝置進行能耗分析，設(shè)計了新型容積式減壓裝置，對此減壓裝置的能量損耗和節(jié)流減壓作了對比，提出了二級容積減壓裝置，并對節(jié)流減壓和二級容積減壓裝置進行了流場分析。

1 節(jié)流減壓裝置能耗與流場分析

1.1 能耗分析

1.1.1可用能損耗數(shù)值計算

氣動汽車一般常用的減壓方式是利用減壓閥等節(jié)流元件使高壓氣體在流動過程中與節(jié)流口產(chǎn)生摩擦以消耗氣體能量來實現(xiàn)減壓，這種減壓方式稱為節(jié)流減壓。這種減壓方式會發(fā)生不可逆轉(zhuǎn)的壓力損失，在節(jié)流減壓過程中流體和外界環(huán)境交換的熱量很少，一般視為絕熱節(jié)流減壓。所以節(jié)流減壓過程中不可避免地會造成可用能損耗[9]。

經(jīng)過節(jié)流減壓的高壓氣體的可用能損耗為

(1)

式中：R為氣體常數(shù)；T0為環(huán)境溫度；p1為初始壓力；p2為終止壓力。

則任意狀態(tài)(p1，T1)的高壓氣體在(p0，T0)的環(huán)境狀態(tài)下，可以完全用來做功的總的可用能為

(2)

式中：T1為初始溫度；p0為環(huán)境壓力；k為氣體多變指數(shù)。

高壓氣體節(jié)流減壓過程可用能損失百分比為

(3)

1.1.2可用能損耗仿真分析

在上述數(shù)學(xué)模型的基礎(chǔ)上，設(shè)定環(huán)境溫度為293 K，壓力為0.101 MPa，在數(shù)學(xué)分析軟件Simulink中分別對節(jié)流前壓力p1為10 MPa，20 MPa，30 MPa，節(jié)流后壓力p2為1～3 MPa的高壓氣體節(jié)流減壓可用能損耗進行仿真[10]，結(jié)果見圖1和圖2。

圖1 節(jié)流減壓過程氣體可用能損失隨節(jié)流前后壓力的變化

圖2 節(jié)流減壓過程單位質(zhì)量氣體可用能損失百分比

從圖1中可以看出，對于不同氣源壓力的氣體，隨著節(jié)流后壓力的減小，單位質(zhì)量的氣體損失的可用能會增加，氣源壓力越高，節(jié)流前后氣體壓差越大，則可用能損失越嚴(yán)重。從圖2中可以看出，在節(jié)流后相同氣體壓力的情況下，節(jié)流前壓力越低，氣體的可用能損耗就越低。

由上述研究可知：當(dāng)氣動汽車節(jié)流前后壓力分別設(shè)置為30 MPa，3 MPa時，氣體可用能損耗依然高達40%，雖然可以通過降低氣源壓力值來減小氣體的可用能損耗，但是降低氣源壓力值的同時氣源可用能的總值也將降低，這將嚴(yán)重降低氣動汽車的續(xù)航能力，因此急需設(shè)計一種新的減壓系統(tǒng)來解決目前的難題。

1.2 流場分析

1.2.1網(wǎng)格劃分及參數(shù)選擇

對閥體流動區(qū)域的幾何結(jié)構(gòu)作了適當(dāng)?shù)娜A角和倒角處理后，在FLUENT前處理軟件ICEM中進行了網(wǎng)格劃分，整體網(wǎng)格采用四面體網(wǎng)格，為保證仿真結(jié)果的準(zhǔn)確性，對閥芯處進行了網(wǎng)格加密處理并設(shè)置了邊界層網(wǎng)格，網(wǎng)格數(shù)量大約4×106，網(wǎng)格劃分如圖3所示。

圖3 節(jié)流減壓閥網(wǎng)格劃分

為了使仿真保持較高的精度，計算模型選擇標(biāo)準(zhǔn)κ-ε模型，模擬仿真算法選擇SIMPLE算法，流場仿真過程采用有限體積法。在三維條件下對節(jié)流減壓閥內(nèi)流動狀況進行仿真分析，其中入口邊界條件設(shè)定為壓力入口(pressure-inlet)，入口壓力設(shè)定為30 MPa，出口邊界條件為壓力出口(pressure-outlet)，出口壓力設(shè)定為3 MPa，其余減壓閥外輪廓及閥芯閥桿面為壁面(wall)。

1.2.2仿真結(jié)果分析

圖4示出減壓閥內(nèi)部流場壓力分布。從圖4a可以看出，高壓腔氣體壓力基本穩(wěn)定在30 MPa，低壓腔上游的壓力穩(wěn)定在4～6 MPa，但中部流速降低，存在湍動能耗散，使得低壓腔出口處壓力基本穩(wěn)定在3 MPa；從圖4b可以看出，減壓口處壓力呈現(xiàn)階梯型的降低趨勢；從圖4c可以看出，減壓閥低壓腔從上游到下游基本呈現(xiàn)壓力降低的趨勢，閥芯上表面錐度位置受壓基本在3 MPa左右。由流場壓力分布圖可以看出：在此工況下，減壓閥低壓腔能輸出穩(wěn)定壓力值，且壓力分布趨勢基本滿足預(yù)期效果。

圖4 減壓閥內(nèi)流道壓力分布

圖5示出減壓閥內(nèi)部流場速度分布。從圖5a可以看出，高壓腔氣體流線基本能均勻充滿減壓閥內(nèi)流道，閥口處流線較稀疏，節(jié)流效果較明顯。從圖5b可以看出，整體速度流場分為三個階段，高壓腔內(nèi)氣體流速基本穩(wěn)定在30 m/s左右，為亞音速狀態(tài)，在節(jié)流口后2～4 mm處氣體流速達到537 m/s，為超音速狀態(tài)，此后沿著閥桿速度逐漸下降，出口處最大流速為260 m/s，為跨音速狀態(tài)。從圖5c可以看出，在閥芯上端面出現(xiàn)了漩渦，主要原因是由于閥芯錐面過渡為平面時氣流方向不改變引起的，但是該漩渦并不影響壓力的穩(wěn)定輸出。

圖5 減壓閥內(nèi)流場速度分布

綜合上述節(jié)流減壓流場分析可知，節(jié)流減壓過程較穩(wěn)定，壓力輸出穩(wěn)定，節(jié)流效果明顯。

2 容積減壓裝置能耗分析

2.1 裝置介紹

容積式減壓裝置在氣動汽車氣動系統(tǒng)中的作用如圖6所示，容積式減壓主要有定容積式和變?nèi)莘e式兩大類。本研究主要針對定容積式減壓裝置進行。

圖6 容積式減壓裝置在氣動系統(tǒng)中的應(yīng)用

2.2 數(shù)值計算

經(jīng)過容積減壓的高壓氣體的可用能損耗為

(4)

式中：cp為氣體比定壓熱容；T0為環(huán)境溫度；T1為初始溫度；k為氣體多變指數(shù)；Z為氣體壓縮因子；R為氣體常數(shù)；p1為初始壓力；p2為終態(tài)壓力。

任意狀態(tài)(p1，T1)的高壓氣體在(p0，T0)的環(huán)境狀態(tài)下，可以完全用來做功的總可用能為

(5)

式中：R為氣體常數(shù)；k為氣體多變指數(shù)；T0為環(huán)境溫度；T1為初始溫度；p0為環(huán)境壓力；p1為初始壓力。

高壓氣體容積減壓過程可用能的損耗比例為

(6)

2.3 仿真分析

在上述數(shù)學(xué)模型的基礎(chǔ)上，設(shè)定環(huán)境溫度為293 K，壓力為0.101 MPa，在數(shù)學(xué)分析軟件Simulink中分別對節(jié)流前壓力為10 MPa，20 MPa，30 MPa，節(jié)流后壓力為1～3 MPa的高壓氣體容積式減壓可用能損耗進行仿真，結(jié)果如圖7和圖8所示。

從圖7中可以看出，對于不同氣源壓力的氣體，隨著減壓后壓力的減小，單位質(zhì)量氣體損耗的可用能會增加，氣源壓力越高，減壓前后氣體壓差越大，則可用能損耗越嚴(yán)重。從圖8中可以看出，容積式減壓和節(jié)流式減壓相比單位氣體可用能損耗率減少了6%～22%。同一氣源壓力下，減壓后壓力越低，容積式減壓相對節(jié)流式減壓節(jié)能效果越好。

圖7 容積式減壓可用能損耗隨減壓前后壓力值的變化

3 兩級容積式減壓流動分析

3.1 減壓原理

圖9示出兩級容積減壓系統(tǒng)原理。每級減壓裝置分別由一個高頻電磁閥和一個減壓氣罐構(gòu)成，兩級減壓裝置串聯(lián)在一起，分別將各自氣體狀態(tài)反饋給控制器完成減壓控制。圖中p0為氣源壓力，q1，q2和q3分別為一、二級減壓氣罐入口和二級減壓出口氣體的質(zhì)量流量，T1和T2分別為兩級減壓氣罐內(nèi)氣體的溫度，V1和V2分別為兩級減壓氣罐的容積。兩級減壓裝置中的減壓過程均受控制器控制，各自完成同時充放氣和只有放氣兩個過程的交替，實現(xiàn)減壓。

圖9 兩級容積減壓系統(tǒng)原理

3.2 流場分析

容積式減壓利用高壓氣體膨脹減壓，采用流體仿真軟件FLUENT對兩級容積式減壓過程進行CFD仿真分析，目的在于探究減壓腔內(nèi)部氣體流動結(jié)構(gòu)分布，探明其減壓機理。

由于減壓過程屬于高壓超音速流動，流動過程設(shè)置為可壓縮流動，即選擇密度求解器，計算模型采用大渦模擬模型，仿真算法選擇SIMPLE算法，流場仿真過程采用有限體積法，采用結(jié)構(gòu)化網(wǎng)格并使其滿足y+條件，網(wǎng)格數(shù)量設(shè)定為5×106左右，流動過程求解質(zhì)量、動量及能量方程。兩級減壓過程邊界條件均設(shè)定入口邊界為壓力入口(pressure-inlet)，其余邊設(shè)定為壁面(wall)，一級減壓入口壓力20 MPa，初始壓力為大氣壓，當(dāng)瓶內(nèi)平均壓力達到8 MPa，入口關(guān)閉。二級減壓入口壓力為8 MPa，初始壓力為大氣壓，當(dāng)瓶內(nèi)平均壓力達到3 MPa時入口關(guān)閉。

氣瓶內(nèi)壓力分別達到8 MPa和3 MPa時速度分布如圖10所示。對比圖10a、圖10b發(fā)現(xiàn)，高壓緩沖較劇烈，因此噴射段越短，瓶內(nèi)上游速度分布越紊亂，形成碎渦越多。雖然低壓減壓相比高壓減壓入口噴射更長，但瓶內(nèi)中下游氣體流速更加均勻，速度梯度變化較小。

圖10 兩級減壓入口關(guān)閉前一刻速度云圖

氣瓶內(nèi)壓力分別達到8 MPa和3 MPa時的壓力分布如圖11所示。對比圖11a、圖11b發(fā)現(xiàn)，由于二級減壓入口壓力較低，其下游速度變化梯度較小，引起的湍流強度較小，因此二級減壓下游壓力變化梯度小。

通過對氣源壓力為20 MPa的兩級減壓過程進行流動仿真，發(fā)現(xiàn)第二級相對于第一級減壓容器內(nèi)氣體壓力場及速度場分布更均勻，減壓過程更加穩(wěn)定。

4 結(jié)論

a) 對節(jié)流減壓裝置在氣動汽車上的應(yīng)用作了能耗分析，發(fā)現(xiàn)節(jié)流式減壓的可用能損失高達27%～59%，且出口壓力越低，可用能損失越大，在30 MPa的氣源壓力和3 MPa的出口壓力最常用工況下，能量損耗高達40%，能耗偏大；

b) 節(jié)流減壓過程較穩(wěn)定，能穩(wěn)定輸出壓力，節(jié)流效果明顯；

c) 容積式減壓能耗損失相對于節(jié)流減壓減少了6%～22%，改進效果明顯；

d) 對于兩級容積減壓，第二級相對于第一級減壓容器內(nèi)氣體壓力場及速度場分布更均勻，減壓過程更加穩(wěn)定。