廖志濤　汪怡平　李　浩

(武漢理工大學(xué)汽車學(xué)院1)　武漢　430070)　(一汽大眾汽車有限公司2)　長春　130011)

基于CFD的某微型客車背部塵土污染機(jī)理研究

廖志濤1)汪怡平1)李浩2)

(武漢理工大學(xué)汽車學(xué)院1)武漢430070)(一汽大眾汽車有限公司2)長春130011)

摘要：針對汽車車身塵土污染現(xiàn)象，介紹汽車表面污染的研究意義及現(xiàn)狀.以某微型客車為研究對象，利用數(shù)值模擬方法研究車身周圍的流場結(jié)構(gòu)，分析引起背部塵土污染的因素并對背部塵土污染結(jié)果進(jìn)行預(yù)測.采用Realizable k-ε湍流模型，以及DPM(離散相)模型，探討塵土顆粒在流場中的受力情況，在朗格朗日坐標(biāo)系下對顆粒運(yùn)動(dòng)微分方程積分得到其運(yùn)動(dòng)軌跡，從而揭示汽車表面污染的形成機(jī)理.針對微型客車背部污染形成機(jī)理，提出幾種減少背部污染的措施.

關(guān)鍵詞：表面污染；數(shù)值模擬；湍流；離散相；運(yùn)動(dòng)軌跡

0引言

隨著道路條件的變化，以及人們對汽車乘坐舒適性要求的提高，汽車表面污染問題已經(jīng)成為空氣動(dòng)力學(xué)的重要方面.汽車在使用過程中難免會接觸泥水和塵土，如果汽車車身周圍的流場結(jié)構(gòu)不合理，隨車身周圍氣流運(yùn)動(dòng)的泥土或者顆粒會附著在汽車表面，不僅影響美觀，還會影響駕駛員視野，從而威脅行車安全[1].與此同時(shí)，汽車電子化的發(fā)展趨勢使汽車裝備越來越多的攝像設(shè)備，這些設(shè)備表面的清潔程度會影響設(shè)備的工作性能，以及乘客的舒適性[2].保持車身表面光潔，會降低洗車的頻率，從而減少汽車的使用成本.因此，探究汽車行駛時(shí)周圍的流場分布，分析車身表面的壓力、速度分布，研究塵土粒子的運(yùn)動(dòng)軌跡，具有重要的意義.

隨著汽車空氣動(dòng)力學(xué)的發(fā)展，汽車表面塵土污染已經(jīng)成為越來越活躍的研究熱點(diǎn).2011年捷豹路虎公司的Adrian Gaylard用格子玻爾茲曼方法對轎車和SUV表面污染進(jìn)行了數(shù)值模擬，揭示引起車身側(cè)壁及背部污染因素.2013年美國EXA公司的Gaylard等[3]針對某SUV車型背部塵土污染問題進(jìn)行風(fēng)洞試驗(yàn)，在后處理過程中用紫外線染色劑顯示背部污染物厚度，并輔以道路試驗(yàn)驗(yàn)證了數(shù)值模擬方法對結(jié)果預(yù)測的可靠性.2014年上汽通用五菱公司的Yang等[4]針對對五菱榮光底部結(jié)構(gòu)進(jìn)行空氣動(dòng)力學(xué)分析，對擋泥板優(yōu)化設(shè)計(jì)和增加車身底部發(fā)動(dòng)機(jī)蓋，減小整車氣動(dòng)阻力和發(fā)動(dòng)機(jī)遭受泥土污染的程度.2012和2014年，金益鋒等[5-6]分別對某微型客車和SUV的車身污染及降低污染的措施進(jìn)行了充分的探討.

本文在前人研究的基礎(chǔ)上建立一種有效的數(shù)值模擬方法，對某微型客車背部污染機(jī)理進(jìn)行分析，為后續(xù)的相關(guān)研究提供一定的理論指導(dǎo).

1汽車表面塵土污染概述

汽車在道路上行駛時(shí)，表面污染的來源有3種，見圖1，根據(jù)污染來源的不同可以分為：雨水污染、外部污染和自身污染[7].本文研究的重點(diǎn)集中在汽車行駛過程中自身造成的污染.汽車的表面塵土污染程度與車身外形，以及道路條件有關(guān)，通過合理的外形設(shè)計(jì)可以減少汽車在行駛過程中的塵土污染程度.目前的研究方法主要有試驗(yàn)研究、理論分析，以及數(shù)值模擬等.

圖1　造成汽車表面污染的3種來源

塵埃在車身粘附的過程主要分3個(gè)階段：塵土擴(kuò)散、傳播和顆粒分離.車輛在道路上行駛時(shí)，車身周圍的流場涉及湍流與多相流，并且這種復(fù)雜流動(dòng)現(xiàn)象包含了氣體、液體，以及固體懸浮顆粒等的參與，同時(shí)各相流動(dòng)之間可能還存在耦合作用.要建立一套精確模擬車身周圍流場的理論非常困難，因此針對具體問題進(jìn)行適度的簡化對問題本身具有非常重要的意義.本文研究微型客車背部塵土污染，做出如下假設(shè):(1)考慮汽車在干燥路面條件下的車身表面污染問題，將車身周圍的流場簡化為氣相與顆粒相混合的兩相流動(dòng)；(2)流場湍流流動(dòng)是充分發(fā)展的，氣相流場是無溫差且恒定的；(3)顆粒相的體積分?jǐn)?shù)小于10%，不考慮顆粒相之間的相互作用.

微型客車受車身設(shè)計(jì)及路面條件影響，其車身側(cè)壁，以及背部非常容易遭受塵土的污染，其中背部是污染較為嚴(yán)重的區(qū)域.背部污染的來源主要有：(1)旋轉(zhuǎn)的后車輪將地面上的塵土卷起，隨著車身周圍氣流運(yùn)動(dòng)繞過C柱流向汽車背部形成污染；(2)由于車輪輪腔內(nèi)外存在壓力差，輪腔內(nèi)的液滴或塵土顆粒受壓差力作用進(jìn)入到車身周圍流場中，隨氣流進(jìn)入車輪尾部區(qū)域形成污染；(3)進(jìn)入底盤與地面間隙的氣流攜帶著顆粒由于受空間和底盤部件的影響形成阻滯，從車身兩側(cè)或者車輛尾部溢出并向上流動(dòng)，并附著在車身表面.

2數(shù)值模擬理論

2.1基本方程及湍流模型方程

汽車車速遠(yuǎn)低于聲速，因此將汽車周圍的流場看作三維不可壓縮粘性等溫流場，研究方法一般是采用雷諾時(shí)均的N-S方程.研究汽車表面塵土污染問題時(shí)，車身周圍的流場除了滿足連續(xù)性方程和N-S方程外，還必須滿足湍流模型方程.Realizablek-ε湍流模型考慮了旋轉(zhuǎn)和曲率對流動(dòng)的影響，能有效模擬旋轉(zhuǎn)剪切流、含有射流的自由流以及邊界層流動(dòng)等[8].本文采用Realizablek-ε湍流模型來分析汽車表面塵土污染.

2.2離散相動(dòng)力學(xué)模型

在車身周圍夾雜著顆粒的氣流中，離散相的體積分?jǐn)?shù)小于10%，因此采用拉格朗日方法對顆粒的軌跡進(jìn)行追蹤.汽車塵土污染研究中，顆粒的運(yùn)動(dòng)軌跡對塵土污染的程度至關(guān)重要，在拉氏坐標(biāo)系下對顆粒作用力微分方程進(jìn)行積分可以求解顆粒的軌跡.顆粒的作用力平衡方程(顆粒慣性=作用在顆粒上的各種力)在笛卡爾坐標(biāo)系下的形式(X方向)為

(1)

式中：u為流體相的速度，m·s-1；up為顆粒相速度，m·s-1；ρ為流體密度，kg·m-3；ρp為顆粒密度，kg·m-3,FD(u-up)為顆粒的單位質(zhì)量曳力，

(2)

在車身近壁面該力主要為由速度梯度引起的升力[10]：

(3)

其中:υ1/2為流體運(yùn)動(dòng)粘度；dp為顆粒直徑；up′為顆粒下一點(diǎn)的速度；dij，dlk，dkl為流體變形速度張量[11]；K=2.894.本文研究汽車在干燥路面條件下行駛時(shí)的塵土污染問題，對塵土顆粒只考慮重力、氣流阻力和流場速度梯度引起的升力.

3數(shù)值仿真計(jì)算

3.1計(jì)算模型的建立

本文以某微型客車為研究對象，在不影響計(jì)算精度的前提下，對三維數(shù)模進(jìn)行適當(dāng)?shù)暮喕?，忽略門把手、后視鏡、雨刮器等附件，并對底盤進(jìn)行平整化處理.為使來流均勻以及尾流充分發(fā)展，在確定計(jì)算域時(shí)使來流口距車身3倍車長，出流口距車尾7倍車長，左右各3倍車寬，計(jì)算域高度為5倍車高.仿真求解的計(jì)算域見圖2.

圖2　仿真求解的計(jì)算域

3.2網(wǎng)格劃分

利用ANSYS ICEM CFD采用四面體和三棱柱的混合網(wǎng)格形式對計(jì)算區(qū)域進(jìn)行網(wǎng)格劃分，在車身附近利用密度盒進(jìn)行加密處理，此處參數(shù)變化梯度較大；在參數(shù)變化梯度小的地方采用較稀疏的網(wǎng)格，節(jié)約計(jì)算資源，兼顧精度與計(jì)算效率；同時(shí)為滿足壁面函數(shù)的要求，在車

身表面劃分與車身表面平行的三棱柱網(wǎng)格來模擬附面層，網(wǎng)格劃分完成后網(wǎng)格數(shù)量達(dá)500多萬.

3.3邊界條件設(shè)置以及求解

本論文采用基于壓力的瞬態(tài)求解器，湍流模型為Realizablek-ε湍流模型，標(biāo)準(zhǔn)壁面函數(shù)，二階迎風(fēng)格式對控制方程進(jìn)行離散化，壓力-速度耦合方式為SIMPLE算法.邊界條件設(shè)置見表1.

表1　邊界條件參數(shù)設(shè)置

3.4離散相模型的設(shè)置

離散相模型參數(shù)包括離散相噴射口、初始條件、邊界條件和介質(zhì)屬性等，這些參數(shù)的正確設(shè)置對汽車表面污染模擬的結(jié)果至關(guān)重要.根據(jù)造成微型客車背部污染的3方面的來源，本文進(jìn)行數(shù)值計(jì)算時(shí)主要設(shè)置3個(gè)方面的離散相噴射口，輪胎噴射口、輪腔噴射口和底盤尾部離散相噴射口等.考慮到輪胎前1/4胎面攜帶的顆粒較少，并且攜帶的少量顆粒速度方向朝向地面，因此輪胎噴射口設(shè)置為3/4胎面.離散相噴射口見圖3.

圖3　離散相噴射口

離散相模型的參數(shù)設(shè)置是在set injection properties中完成的，對射流源的射入類型、顆粒的直徑分布、射入位置以及速度等進(jìn)行相關(guān)設(shè)置，仿真過程考慮離散相與連續(xù)相之間的耦合作用，采用面射流源的噴射方式，車輪噴射源、輪腔噴射源和底盤噴射源的質(zhì)量流率分別為0.01，0.008和0.005 kg/s，具體參數(shù)設(shè)置見表2.

表2　離散相模型參數(shù)設(shè)置

離散相模型的邊界條件有捕捉、逃逸、反彈和通過等[12]，對計(jì)算域入口和出口設(shè)置為逃逸，其余壁面設(shè)置為反彈，顆粒直徑分布服從Rosin-rammler分布，最小粒子直徑為3×10-5m，最大粒子直徑5×10-4m，平均粒子直徑為1.5×10-4m.

4汽車尾部污染分析

4.1尾部區(qū)域速度分布

后車輪及輪腔速度流線見圖4，從圖中可以清晰地看出輪胎和輪腔引起的氣流經(jīng)過汽車側(cè)壁后部繞過制動(dòng)燈進(jìn)入車輛尾部流場區(qū)域，可以看到在車輛側(cè)壁和尾部的氣流分離區(qū)速度較高.因此，從后輪及輪腔噴射面射出的顆粒一方面加重了汽車側(cè)壁后部區(qū)域的污染，另外一方面會隨著氣流進(jìn)入尾部流場對汽車背部造成污染.

圖4　后輪及輪腔噴射面速度流線圖

底盤噴射面引起的速度流線見圖5，氣流經(jīng)過底盤末端繞到尾部流場區(qū)域，攜帶在氣流中的污染顆粒會隨著氣流沉積在車輛背部區(qū)域.為了更好地展示車身周圍的流場結(jié)構(gòu)，選取縱向?qū)ΨQ面A-A和Z=0.5 m的截面，圖6是縱向?qū)ΨQ面內(nèi)的速度矢量圖，可以看到車輛尾部存在一個(gè)較大的氣流低速區(qū)域，區(qū)域內(nèi)存在兩個(gè)向內(nèi)地對稱旋轉(zhuǎn)的渦，攜帶著污染顆粒的氣流隨著下部的渦上卷，附著在汽車背部，進(jìn)而形成表面污染.圖7是Z=0.5 m的截面內(nèi)速度矢量圖，從圖中可以看出尾部中央速度比兩側(cè)低，速度較低的區(qū)域內(nèi)存在方向相反的渦，污染顆粒從兩側(cè)隨氣流向中央運(yùn)動(dòng).

圖5　底盤噴射面速度流線圖

圖6　縱向?qū)ΨQ面速度矢量圖

圖7　Z=0.5 m截面內(nèi)速度矢量圖

4.2離散相顆粒的運(yùn)動(dòng)軌跡

圖8為某一時(shí)刻從后車輪和輪腔射出的顆粒的軌跡圖，從圖中可以看出顆粒在車輛側(cè)壁尾部濃度較高，并且部分顆粒隨著時(shí)間的推移會隨著氣流向車輛尾部運(yùn)動(dòng).為了更清晰地展示汽車背部污染的動(dòng)態(tài)過程，圖9截取了仿真過程中的4個(gè)時(shí)刻從底盤噴射口射出的粒子的軌跡圖，可以看出從底盤噴射口射出的粒子是造成汽車背部污染的主要原因，隨著仿真時(shí)間的不斷推進(jìn)，一方面更多的粒子從噴射面射出，另外一方面已經(jīng)存在流場中的顆粒隨著氣流進(jìn)入尾流區(qū)，并不斷上卷，在t=0.09 s時(shí)后保險(xiǎn)杠區(qū)域已經(jīng)存在了明顯污染，并且污染區(qū)域有向背部上方擴(kuò)散的趨勢.

圖8　車輪及車腔顆粒軌跡追蹤圖

圖9　底盤噴射面顆粒在4個(gè)不同時(shí)刻運(yùn)動(dòng)軌跡追蹤圖

4.3縱向?qū)ΨQ面內(nèi)的壓強(qiáng)分布

圖10為車輛尾部縱向?qū)ΨQ面的壓力分布云圖，由圖中可以看出大部分的區(qū)域的壓強(qiáng)都比標(biāo)準(zhǔn)大氣壓小，并且壓力從較遠(yuǎn)的位置到尾部不斷減小，氣流從車身頂部和底部匯集到尾部區(qū)域并形成低壓的渦區(qū)，而且渦區(qū)的旋轉(zhuǎn)使低壓有向整個(gè)車身背部表面擴(kuò)散的趨勢，正是這種壓力差及擴(kuò)散作用使攜帶在氣流中的塵埃沉積在背部車身表面，從而形成污染.

圖10　縱向?qū)ΨQ面內(nèi)的壓力分布云圖

5結(jié)論

1) 由底盤噴射源射出的顆粒是造成汽車背部污染的主要原因，由于車身尾部存在渦流區(qū)，粒子隨氣流在此運(yùn)動(dòng)會發(fā)生翻卷并在壓差和擴(kuò)散的作用下附著在背部車身表面.

2) 從后車輪和輪腔射流源射出的顆粒一方面對汽車側(cè)壁后部造成污染，另外一方面部分污染顆粒會隨著氣流進(jìn)入尾流區(qū)，加重背部區(qū)域的污染.

3) 針對汽車背部污染的來源，一方面可以采取措施以減少進(jìn)入地面間隙的氣流，如安裝擾流器和前擋板等裝置；另外一方面可以將車頂?shù)纳线吔堑箞A或在車頂上邊角安裝導(dǎo)流片，以使從車頂流抑制底部流向上翻起并對后窗上的塵土進(jìn)行沖刷.

參 考 文 獻(xiàn)

[1]谷正氣.汽車空氣動(dòng)力學(xué)[M].北京：人民交通出版社,2005.

[2]JILESEN J, GAYLARD A P, DUNCAN B, et al. Simulation of rear and body side vehicle soiling by road sprays using transient particle tracking[R]. SAE Technical Paper,2013.

[3]GAYLARD A P, DUNCAN B. Simulation of rear glass and body side vehicle soiling by road sprays[J]. SAE International Journal of Passenger Cars-Mechanical Systems,2011,4(1):184-196.

[4]YANG W, ZHOU X, PENG J, et al. Optimization of the underbody layout of a small van for better aerodynamics using digital simulation[R]. SAE Technical Paper,2014.

[5]金益鋒.基于CFD的某微型客車減阻與塵土污染研究[D].長沙：湖南大學(xué),2012.

[6]廖磊.車輪濺水及其對車身表面污染的仿真研究[D].長春：吉林大學(xué),2014.

[7]HAGEMEIER T, HARTMANN M, THéVENIN D. Practice of vehicle soiling investigations: A review [J]. International Journal of Multiphase Flow,2011,37(8):860-875.

[8]張師帥.計(jì)算流體動(dòng)力學(xué)及其應(yīng)用:CFD軟件的原理與應(yīng)用[M].武漢：華中科技大學(xué)出版社,2011.

[9]朱紅鈞,林元華.Fluent 12流體分析及工程仿真[M].北京：清華大學(xué)出版社,2011.

[10]谷正氣,王振,張勇，等.某轎車車身塵土污染數(shù)值仿真與控制[J].長安大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,35(5):145-152.

[11]羅昔聯(lián),顧兆林.基于DPM模型的街谷內(nèi)顆粒物擴(kuò)散特性研究[J].中國科學(xué)院研究生院學(xué)報(bào),2007,24(5):578-583.

[12]張萬濤,余宏明.正交試驗(yàn)設(shè)計(jì)方法在庫岸滑坡敏感性分析中的應(yīng)用[J].安全與環(huán)境工程,2009,16(5):13-16.

Investigation of Rear Surface Contamination Mechanism of a Micro-bus Based on CFD

LIAO Zhitao1)WANG Yiping1)LI Hao2)

(AutomotiveEngineeringCollegeofWuhanUniversityofTechnology,Wuhan430070,China)1)(FAW-VolkswagenAutomotiveCo.Ltd,Changchun130011,China)2)

Abstract：Taking a micro-bus as the study target, the numerical simulation method is utilized to investigate the flow structure around car body. The factors that cause rear surface contamination are observed and the simulation results are predicted. A combination of Realizable k-ε turbulence model and DPM is adopted and the particle force condition is considered. The trajectory of the particles is displayed with integral on particle motion differential equations on the Lagrange coordinate to reveal the formation mechanism of the surface contamination. Finally, several countermeasures are proposed to control the rear surface dirt deposition based on the mechanism of contamination.

Key words：surface contamination; numerical simulation; turbulence; DPM; trajectory

收稿日期：2016-03-12

中圖法分類號：U461.1

doi:10.3963/j.issn.2095-3844.2016.03.035

廖志濤(1991- )：男，碩士生，主要研究領(lǐng)域?yàn)槠嚳諝鈩?dòng)力學(xué)