程秀，莫健華，明志新，劉祖強(qiáng)

（1.華中科技大學(xué) 材料科學(xué)與工程學(xué)院，湖北 武漢 432700；2.湖北水利水電職業(yè)技術(shù)學(xué)院 機(jī)電工程系，湖北 武漢 432700）

機(jī)動(dòng)車加油管注油過(guò)程的Workbench Fluent仿真

程秀1,2，莫健華1，明志新2，劉祖強(qiáng)2

（1.華中科技大學(xué) 材料科學(xué)與工程學(xué)院，湖北 武漢 432700；2.湖北水利水電職業(yè)技術(shù)學(xué)院 機(jī)電工程系，湖北 武漢 432700）

隨著環(huán)保要求的日益嚴(yán)格，機(jī)動(dòng)車加油過(guò)程所產(chǎn)生的油氣排放污染控制越來(lái)越嚴(yán)格. 針對(duì)此問(wèn)題，選擇一體化軟件Workbench 對(duì)簡(jiǎn)化后的油管模型進(jìn)行網(wǎng)格劃分，再使用VOF多相流模型和RNGk-ε模型對(duì)油管的注油過(guò)程進(jìn)行模擬計(jì)算. 設(shè)置好邊界條件以及求解條件，建立仿真模型，對(duì)不同內(nèi)徑加油管在不同加油速度下的注油過(guò)程進(jìn)行兩相流模擬計(jì)算，討論其對(duì)液封形成的影響. 結(jié)果表明，加油速度越大，加油管徑越小，自液封越容易形成，但是速度大、管徑小容易導(dǎo)致反噴.

加油排放；計(jì)算流體力學(xué)；數(shù)值模擬；油管液封

在加油過(guò)程中，汽油蒸發(fā)的排放占整個(gè)機(jī)動(dòng)車油氣排放的20%，造成大量能源的損失和大氣的污染，因此必須重視[1-6]. 車載油氣回收系統(tǒng)（On-board Refueling Vapor Recovery，ORVR）在美國(guó)首先提出來(lái)，系統(tǒng)要求在加油管某處完全充滿汽油以形成自液封，使油氣不向大氣中排放而只能被液態(tài)汽油帶入油箱內(nèi)，以減少機(jī)動(dòng)車加油過(guò)程中油氣的排放[2,7].

機(jī)動(dòng)車注油過(guò)程瞬間萬(wàn)變，有三維的瞬態(tài)流動(dòng)、氣液相變流動(dòng)，還有自由表面紊流等. 計(jì)算流體力學(xué)CFD的快速發(fā)展，使得仿真機(jī)動(dòng)車加油過(guò)程變得簡(jiǎn)單可行. 美國(guó)的Stoneman教授[8]和加拿大溫莎大學(xué)[9]都利用CFD對(duì)加油管的加油過(guò)程進(jìn)行了模擬仿真，取得了顯著的成果，但是受限于早期計(jì)算機(jī)的水平，模擬結(jié)果跟實(shí)際加油過(guò)程有一定的差距，而且相關(guān)研究結(jié)果多是二維模擬. 中國(guó)在2008年才有人開展這方面的工作，國(guó)內(nèi)做得較好的是江蘇大學(xué)的何仁教授課題組以及北京石油化工學(xué)院的陳家慶教授課題組. 他們采用CFD數(shù)值模擬對(duì)機(jī)動(dòng)車加油過(guò)程中氣液兩相流動(dòng)特性進(jìn)行了分析，但他們使用的網(wǎng)格劃分軟件和仿真計(jì)算軟件是獨(dú)立的，從而影響了計(jì)算精度[2,10]. 本研究借助網(wǎng)格劃分和仿真計(jì)算一體化的軟件Workbench Fluent，對(duì)不同內(nèi)徑加油管在不同加油速度情況下的液封性能進(jìn)行數(shù)值模擬，分析研究了加油管內(nèi)徑和加油速度對(duì)液封性能的影響.

1 模型介紹

1.1 加油管的幾何模型

某型號(hào)的汽車加油管的幾何模型如圖1所示，4個(gè)拐角的存在有利于液封的形成. CFD數(shù)值模擬使用的模型在保證主要結(jié)構(gòu)特征的前提下需盡量簡(jiǎn)化并轉(zhuǎn)化成包裹的實(shí)體，對(duì)圖1的加油管進(jìn)行簡(jiǎn)化得到圖2所示的簡(jiǎn)化模型，保留4個(gè)拐角幾何結(jié)構(gòu)特征不變. 加油槍的噴嘴直徑是標(biāo)準(zhǔn)的，需要與加油管的管口相配套，根據(jù)市場(chǎng)使用情況選擇直徑為20 mm的噴嘴. 加油管的直徑選擇33 mm和25 mm兩種.

圖1 加油管的幾何模型

圖2 簡(jiǎn)化的加油管幾何模型

1.2 網(wǎng)格劃分和相關(guān)設(shè)置

采用Workbench Fluent 自帶的劃分網(wǎng)格軟件對(duì)加油管模型進(jìn)行網(wǎng)格劃分. 為減少計(jì)算時(shí)間和工作量只計(jì)算加油管部分，全部采用六面體網(wǎng)格以保證精度. 網(wǎng)格如圖3所示，不管是截面（圖3-a）還是表面（圖3-b)，網(wǎng)格規(guī)整，無(wú)明顯變形區(qū)域，說(shuō)明網(wǎng)格劃分合理.

圖3 油管網(wǎng)格模型

邊界條件設(shè)置：加油槍嘴的油出口設(shè)為速度入口（Velocity Inlet），加油管末端設(shè)為壓力出口（Pressure Outlet）；加油管的全部固體外壁面和加油槍噴嘴段的外壁作壁面(Wall) 處理. 速度入口處液態(tài)汽油的體積百分比設(shè)為1，湍流強(qiáng)度都取為5% ，水力直徑取各自的管道直徑.

空氣和汽油為計(jì)算材料. 空氣的密度為1.225 kg/m3，粘度為1.789 4×10-5kg· m-1·s-1；汽油的密度為720 kg/m3，粘度為5.4×10-4kg· m-1·s-1.

1.3 數(shù)學(xué)模型

汽油從加油管入口注入油管，并沿著加注管流向第一個(gè)拐角處，在重力及慣性作用下順流而下，依次流過(guò)第二至第四拐角，最后注入油箱. 加油過(guò)程復(fù)雜多變，涉及到非穩(wěn)態(tài)、多相流、湍流等眾多因素，不可能考慮所有的因素[11]. 因此假設(shè)：等溫條件下加油，氣液兩相之間無(wú)熱傳遞、傳質(zhì)現(xiàn)象，忽略浮力[2,7].力，從而導(dǎo)致壓力梯度.

1.3.1 多相流模型選擇

Workbench將求解器看作一個(gè)組件，F(xiàn)luent是其中的一個(gè)，它包括一系列的多相流模型，我們采用VOF（Volume of Fluid）多相流模型. 加油過(guò)程是氣液兩相湍流流動(dòng)，VOF模型能夠準(zhǔn)確計(jì)算層流、液體中有大氣泡流動(dòng)、水壩決堤的水流等復(fù)雜情況下流體的動(dòng)態(tài)特性，并能精確液-氣分界面的穩(wěn)態(tài)或瞬時(shí)分界面[12-13].

在VOF 模型中，氣相和液相可以看成是兩種互不相溶、不可壓縮的流體，其平均質(zhì)量守恒方程和平均動(dòng)量守恒方程[2,7]：

式中，ρ為密度，kg/m3；t為時(shí)間，s；x,u為速度矢量，m/s；P為壓力，Pa；μ為動(dòng)力學(xué)粘度，Pa· s；g為重力加速度，kg· m· s-1；fσ表面張力引起的體積力，N；i,j為坐標(biāo).

混合物的密度和動(dòng)力粘度可以通過(guò)平均體積分?jǐn)?shù)的方法來(lái)求[2]：

式中，α為體積分?jǐn)?shù).

氣液兩相之間的界面會(huì)產(chǎn)生不同的表面張力，而且是突變，如圖4所示. 假設(shè)液滴是一個(gè)球，半徑為r，球缺底半徑為1r.A、B點(diǎn)水平方向合力相互抵消，豎直方向的合力為必將產(chǎn)生一個(gè)附加壓

圖4 氣液界面模型

因此，壓力梯度必須與動(dòng)量方程中額外的體積力相等以維持平衡. 給氣液一個(gè)一定厚度的過(guò)渡區(qū)，該不連續(xù)的壓力突變可以表示為[2,7]：

式中，σ為表面張力，N.

那么在連續(xù)域內(nèi)的氣體體積分?jǐn)?shù)則由體積分?jǐn)?shù)的連續(xù)方程求解[2,7]：

氣液兩相界面處計(jì)算單元為0＜α＜1.

1.3.2 湍流模型的選擇

湍流模型對(duì)計(jì)算精度有顯著影響，因此必須選擇合適的湍流模型. Workbench Fluent求解器中提供的湍流模型包括標(biāo)準(zhǔn)k-ε模型以及其變型RNGk-ε模型、雷諾應(yīng)力模型（Reynolds Stress Model，RSM）和大渦模擬（Large Eddy Simulation，LES）等. 加油管因?yàn)榧佑退俣瓤斐尸F(xiàn)出強(qiáng)湍流狀態(tài)，并且拐角是急彎度，因此選擇RNGk-ε模型，其動(dòng)能方程和能量耗散方程分別如下[2,7,14]：

式中，k為湍動(dòng)能，m2/s2；Gk為平均速度梯度產(chǎn)生的湍動(dòng)能，m2/s2；Gb為浮力產(chǎn)生的湍動(dòng)能，m2/s2；ε為湍流耗散率；常量

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

加油時(shí)汽油首先以一定的速度從加油槍管口噴射而出，在重力及速度慣性的作用下沿著油管向下流動(dòng)，在加油管喉管處汽油與通風(fēng)管處的汽油蒸氣混合，并伴隨空氣的卷吸與汽油的相變，在拐角處發(fā)生碰撞沖擊，最終進(jìn)入油箱. 氣液兩相混合物在重力及速度慣性的作用下流入燃油箱中.

加油管的幾何結(jié)構(gòu)對(duì)加油過(guò)程有影響，但加油管機(jī)械結(jié)構(gòu)受制于車身，所以主要通過(guò)規(guī)范加油速度以及加油管的內(nèi)徑來(lái)實(shí)現(xiàn)自液封和抑制反噴. 文獻(xiàn)[7]提到：加油速度過(guò)小，填充時(shí)間長(zhǎng)，空氣有充分的時(shí)間卷入，液封不易形成；加油速度過(guò)大，又增加提前跳槍停止加油的可能性. 國(guó)標(biāo)中規(guī)定加油速度不能超過(guò)4.5 m/s，本研究選擇2～3 m/s. 已有的計(jì)算結(jié)果表明，汽油充滿整個(gè)油管的時(shí)間小于2 s.

2.1 加油管徑為33 mm的計(jì)算結(jié)果

實(shí)際加油過(guò)程中油管的流體力學(xué)特性與加油速度直接相關(guān). 33 mm的管徑相對(duì)來(lái)說(shuō)比較大，為此，選取加油速度為2 m/s、3 m/s. 加油管油帽段上部有一管口連通大氣，使出口壓力相對(duì)于大氣壓為0. 圖5是油管直徑為33 mm，兩種不同速度下第一拐角橫截面汽油占比云圖（紅色的代表汽油，藍(lán)色代表氣體）. 通過(guò)對(duì)比顏色，估算出汽油體積占比. 從圖5可以看出：速度為2 m/s的汽油體積占比約為80%，速度為3 m/s的汽油體積占比約為75%. 由于自重以及油從加油槍快速?zèng)_出，在第一拐角處發(fā)生激烈碰撞并產(chǎn)生旋流，使得湍流現(xiàn)象更加復(fù)雜，汽油在加油管中在第一拐角處出現(xiàn)的噴濺、碰撞、旋流等現(xiàn)象符合流體動(dòng)力學(xué)的規(guī)律，與機(jī)動(dòng)車加油過(guò)程的特性吻合.

圖5 d=33 mm ，不同速度下第一拐角油氣體積分?jǐn)?shù)截面圖

圖6是兩種加油速度下第一拐角處的速度云圖. 截面的節(jié)點(diǎn)數(shù)是一個(gè)定值，速度是矢量，從圖中可以看出有部分液體的速度因?yàn)榕鲎舶l(fā)生反彈，出現(xiàn)與宏觀速度方向相反的矢量. 粗略估算圖6-a的反向矢量個(gè)數(shù)小于圖6-b反向矢量的個(gè)數(shù)，圖6-b中的拐角速度矢量變向明顯增多，這說(shuō)明速度大，碰撞的幾率增大. 反向矢量多，反彈的幾率大，這必將增大發(fā)生反噴的幾率，從而造成提前“跳槍”.因此，必須嚴(yán)格控制加油槍的最大加油流量，不能一味追求快速，否則欲速則不達(dá).

圖6 d=33 mm ，不同速度下第一拐角的速度云圖

圖7是汽油體積分?jǐn)?shù)云圖，為了對(duì)比明顯，以藍(lán)色顯示體積分?jǐn)?shù)大于50%的氣體體積，完全斷開的空白區(qū)域視為汽油形成的液封. 由圖7可知：兩種速度都可以形成液封，但是速度大者加油管中的汽油體積率也大，形成液封的時(shí)間短、長(zhǎng)度大，液封效果好. 整個(gè)油管充滿油的時(shí)間小于2 s，比較短. 因此本文沒(méi)有探討液封形成需要的時(shí)間. 4個(gè)拐角的存在，使油管的流線復(fù)雜，降低了汽油在管內(nèi)的流速，有利于液封的形成.

圖7 d=33 mm ，不同速度下汽油體積分?jǐn)?shù)云圖（藍(lán)色代表氣體）

2.2 加油管徑為25 mm的計(jì)算結(jié)果

減小油管直徑也可以降低汽油在管內(nèi)的流速. 文獻(xiàn)[7]的結(jié)果顯示：直徑為25 mm的自液封效果較好，因此對(duì)這種規(guī)格的管徑也進(jìn)行仿真計(jì)算. 25 mm的管徑較小，采用加油速度為2 m/s和2.5 m/s.

圖8是管徑為25 mm的油管第一拐角處截面汽油體積占比云圖. 不同速度的汽油體積占比區(qū)分不明顯，但是速度大者的截面汽油分布不規(guī)則，說(shuō)明沖擊攪拌現(xiàn)象劇烈. 由于汽油在第一拐角經(jīng)過(guò)垂直加速后突然改變方向，在第二拐角產(chǎn)生強(qiáng)烈的能量耗散，速度迅速減小，沖擊力減小，有利于液封的形成. 但是管徑過(guò)小，加油速度不變，會(huì)導(dǎo)致在拐角處的劇烈沖擊碰撞，產(chǎn)生汽油回濺.

圖8 d=25 mm ，不同速度下第一拐角油氣體積分?jǐn)?shù)截面圖

從圖9中可以估算出兩種速度下反向矢量占比比較接近. 汽油從加油槍管口噴射流出，碰到加油管壁面后一部分碰撞反彈，一部分沿管壁下流. 管內(nèi)汽油呈強(qiáng)湍流流動(dòng)狀態(tài)，伴隨有碰撞、旋流等現(xiàn)象，速度大者，這個(gè)現(xiàn)象更加明顯，速度反向矢量顯著增加，但是在數(shù)據(jù)統(tǒng)計(jì)的時(shí)候因?yàn)榻y(tǒng)計(jì)方法的誤差，該差別不是很明顯.

圖9 d=25 mm ，不同速度下第一拐角的速度云圖

如圖10所示，油管直徑25 mm的加油管中第一拐角與第二拐角之間形成大段的液封區(qū)域，但是速度大者液封長(zhǎng)度明顯大于速度小者的. 這符合文獻(xiàn)[7]所述. 注油過(guò)程與管徑為33 mm管徑的注油過(guò)程類似，不同之處在于湍流、碰撞、旋流程度不同而已. 要區(qū)分速度對(duì)注油過(guò)程的影響，還需要進(jìn)一步的仿真計(jì)算. 液封形成原理及過(guò)程與前述相同.

圖10 d=25 mm ，不同速度下汽油體積分?jǐn)?shù)云圖（藍(lán)色代表氣體）

由圖7和圖10可知：自液封主要發(fā)生在第一至第二拐角的水平段，空白長(zhǎng)度越大，液封長(zhǎng)度越大，自液封效果越好；管徑相同，速度大者液封長(zhǎng)度大；加油速度相同，管徑小者液封效果稍好. 這個(gè)結(jié)果與文獻(xiàn)[7]的結(jié)果相符，同時(shí)說(shuō)明本文研究模型可靠.

3 結(jié)論

利用Workbench Fluent對(duì)兩種直徑的加油管在不同加油速度的流體力學(xué)特性進(jìn)行仿真計(jì)算，得到以下結(jié)論：

1）保持加油管的4個(gè)拐角幾何特征不變，增大加油速度和減小加油管直徑，有利于自液封的形成，但反噴傾向增大，會(huì)增加提前“跳槍”的可能性.

2）綜合考慮自液封和反噴的情況，相對(duì)較好的加油參數(shù)是加油管徑25 mm、加油速度2.5 m/s.

3）Workbench Fluent通過(guò)合理邊界條件的設(shè)置，可以對(duì)復(fù)雜的加油過(guò)程進(jìn)行有效的計(jì)算模擬，從而預(yù)測(cè)加油過(guò)程容易出現(xiàn)的問(wèn)題，在實(shí)際中加以改進(jìn)，可以有效縮短研發(fā)周期、降低成本.

本研究由于時(shí)間等方面的影響，考慮的參數(shù)不多，因此統(tǒng)計(jì)數(shù)據(jù)存在誤差，幾個(gè)結(jié)論只是定性結(jié)論. 因此在今后的工作中需要增加加油速度、管徑參數(shù)以及油管幾何結(jié)構(gòu)的考慮，以減少統(tǒng)計(jì)數(shù)據(jù)的誤差.

[1] LIU Huan, MAN Hanyang, TSCHANTZ M, et al. VOC from vehicular evaporation emissions∶ status and control strategy [J]. Environmental Science & Technology, 2015, 49(24)∶ 14424-14431.

[2] 汪智，何仁. 基于Fluent仿真的ORVR加油管液封設(shè)計(jì)[J]. 重慶理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2014, 28(2)∶16-21.

[3] YAMADA H, SATOSHI I, HIROSHI T. Refueling emissions from cars in Japan∶ compositions, temperature dependence and effect of vapor liquefied collection system [J]. International Journal of Social Economics, 2015,36(5)∶ 658-678.

[4] YANG Xiaofan，LIU Huan, CUI Hongyang, et al. Vehicular volatile organic compounds losses due to refueling and diurnal process in China∶ 2010-2050 [J]. Journal of Environmental Sciences, 2015, 33(7)∶ 88-96.

[5] 杜建波，方茂東，陸紅雨，等. 加油排放試驗(yàn)程序的設(shè)計(jì)和車載油氣回收系統(tǒng)開發(fā)的建議 [J]. 汽車工程，2014, 36(5)∶ 537-541.

[6] 朱玲，陳家慶，王耔凝. 車載加油油氣回收ORVR系統(tǒng)應(yīng)用進(jìn)展[J]. 油氣儲(chǔ)運(yùn)，2015, 34(5)∶ 469-476.

[7] 湯水清，陳家慶，劉美麗，等，ORVR系統(tǒng)加油管液封性能數(shù)值模擬[J]. 汽車工程學(xué)報(bào)，2016, 6(3)∶ 187-195.

[8] STONEMAN S. On the design of automotive fuel filler pipes [J]. Automotive Engineer, 1997, 22(1)∶ 32-36.

[9] ZHAO A G. A hybrid numerical study of automotive fuel tank filling [D]. Canada∶ University of Windsor, 2003.

[10] 陳家慶，張男，王金惠，等. 機(jī)動(dòng)車加油過(guò)程中氣液兩相流動(dòng)特性的CFD數(shù)值模擬[J]. 環(huán)境科學(xué)，2011,32(12)∶ 3710-3716.

[11] 任碧琪，常耀卿，劉明宇，等. 北京加油站加油速率與車載油氣回收系統(tǒng)的兼容性[J]. 環(huán)境工程，2015,33(S1)∶ 487-490.

[12] ISSA R I. Solution of the implicitly discretised fluid flow equations by operator-splitting [J]. Journal of Computational Physics, 1986, 62(1)∶ 40-65.

[13] GUEYFFIER D, LI Jie, NADIM A, et al. Volume-of-fluid interface tracking with smoothed surface stress methods for three-dimensional flows [J]. Journal of Computational Physics, 1999, 152(2)∶ 423-456.

[14] YAKHOT V, ORSZAG S A. Renormalization group analysis of turbulence I basic theory [J]. Journal of Scientific Computing, 1986, 1(1)∶ 3-51.

[責(zé)任編輯：熊玉濤]

An Oil Injection Process for Motor Vehicle Refueling Based on Workbench Fluent Simulation

CHENG Xiu1,2, MO Jian-hua1, MING Zhi-xin2, LIU Zu-qiang2
(1. School of Materials Science and Engineering, Huazhong University of Science and Technology Wuhan 432700, China; 2. Department of Mechanical and Electrical Engineering, Hubei Water Resource Technical College, Wuhan 432700, China)

With the increasingly stringent environmental requirements, control of pollution from exhaust of oil and gas in motor vehicle fueling process is increasingly stringent. To solve this problem,the integrated software Workbench is selected to generate grid mesh for the simplified tubing model and the VOF model and RNGk-εmodel is applied to simulate and calculate the fuel injection process of oil pipes. The boundary conditions and solution conditions are set and a simulation model is established to simulate and calculate the two-phase flow of the fuel injection process of oil pipes of different sizes at different rates, and its effect on liquid seal is discussed. The results show the greater the fueling rate and the smaller the fuel pipe diameter, the easier it is to form the liquid seal.But the greater rate of fueling and smaller diameter of the tube can easily lead to back spray.

oil and gas emissions from the fueling process; computational fluid dynamics; numerical simulation; filler pipe liquid seal

U473.9

A

1006-7302（2017）04-0032-08

2017-06-13

程秀（1976—），女，湖北廣水人，副教授，博士，主要從事計(jì)算仿真及機(jī)電一體化方面的工作.