強降雨條件下車輛空氣動力學特性分析

殷碩 于夢閣 陳煥明 王亞倫

摘要：? 為研究強降雨環境下的汽車空氣動力學特性問題，本文以標準MIRA階梯背汽車模型為研究對象，建立汽車外流場空氣動力學數值計算模型。同時，利用ANSYS Fluent軟件，研究在不同降雨強度、不同偏航角對車身氣動力、速度特性、壓力特性的影響，并追蹤了離散相雨滴的運動軌跡及其與車身的接觸，分析了車身表面附近區域雨滴質量濃度分布和變化。研究結果表明，在無降雨環境下，隨著車速的增加，車身表面正壓區和負壓區覆蓋范圍均明顯增大，部分雨滴在車身表面發生明顯的反彈現象，車速對氣動阻力影響較大，而對氣動阻力系數沒有實質性影響;相同車速條件下，降雨強度越大，氣動阻力（系數）增幅越大，并且與降雨強度近似呈線性關系;在一定降雨強度下，隨著偏航角的增大，整車氣動阻力（系數）和氣動側力（系數）相比無雨時增大，左側面表面壓力逐漸增大，汽車后側窗邊緣區域壓力減小，車頂高流速區域覆蓋范圍增大，并向車身右側延伸，車身右側低流速區域覆蓋范圍變大;降雨使得汽車左側迎側風區域的表面壓力增大，尾箱蓋區域的壓力分布發生改變，驗證了該模型的準確性。該研究為強降雨環境下汽車氣動研究和應用提供了參考。

關鍵詞：? 數值計算; 車輛氣動特性; 車速; 降雨強度; 偏航角

中圖分類號： U461.1? 文獻標識碼： A

收稿日期： 20210508; 修回日期： 20210620

基金項目： 山東省高等學校科技計劃資助項目（J18KA048）

作者簡介：? 殷碩（1997），男，碩士研究生，主要研究方向為車輛空氣動力學。

通信作者：? 于夢閣（1985），女，博士，副教授，主要研究方向為車輛空氣動力學和系統動力學。Email： yumengge0627@163.com

汽車空氣動力學分析與汽車的安全性、動力性、燃油經濟性以及美觀性關系密切，行駛過程中汽車所處的外部環境較為復雜，在強風雨條件下，汽車氣動性及燃油經濟性變差，對汽車行駛安全造成影響。世界衛生組織指出每年道路交通安全事故造成約上百萬人意外死亡[1]，中國道路交通事故統計年報統計數據顯示，一部分交通事故與降雨等惡劣環境有關[23]，因此風雨環境下的道路交通安全問題應加以重視。國內外學者對降雨條件下汽車空氣動力學領域進行了相關研究。李長城等人[4]分析了降雨強度大小對平均車速以及車速概率學分布規律的影響，并采用非線性回歸方法建立降雨條件下的車速預測模型，獲得了較高的擬合度;馮偉等人[5]采用數值模擬的方法研究了轎車在降雨、側風等環境下的流場特性，結果表明降雨使阻力和縱傾力矩增加，升力和橫擺力矩增加，可能會影響轎車行駛安全;王保成[6]研究了降水對高速公路車輛的影響因素;馬筱櫟[7]利用HLM模型分析了降水對山區高速公路不同路段的影響，發現平曲線半徑增大會對該路段運行車速起到積極作用，而坡度增大會對運行車速產生消極影響。另外在車輛空氣動力學其他相關領域也有相關研究成果。唐文淵[8]利用ANSYS Fluent對車身外部造型進行了分析，并依據汽車外流場性質，對汽車外形進行了優化;劉旺等人[9]利用STARCCM+研究了尾翼對汽車空氣動力學性能的影響，結果顯示，加裝尾翼可顯著提高車輛的空氣動力學性能;王毅剛等人[10]進行了實車風洞實驗，研究了汽車外流場特性和噪聲分布特性;方剛等人[11]推導出了基于橢圓形方程的車輛外流場網格計算生成方法，為車輛流場分析提供了重要參考。以上研究對降雨環境下車輛外流場涉及較少，尤其是短期內大降雨強度情況、偏航角及強風雨等作用于汽車具體機理的相關研究。基于此，本文建立汽車外流場數值計算模型，利用ANSYS Fluent軟件，分析了不同降雨強度和不同偏航角對車身氣動力、速度特性、壓力特性的影響，追蹤了離散相雨滴的運動軌跡及其與車身的接觸，分析了車身表面附近區域雨滴質量濃度的分布和變化。該研究為強降雨環境下的汽車氣動分析提供了理論依據。

1 控制方程

1.1 連續相模型

在實際工程應用中，基于雷諾方程的湍流模型廣泛應用，汽車外流場馬赫數較低，其流動為定常不可壓縮流動，因此應用雷諾時均NavierStokes方程對流場進行求解，其控制方程如下[12]：

連續性方程為

Uixi=0（1）

動量方程為

xiUiUj=-1ρPxi+v2Uixjxj-xiu′iu′j（2）

式中，Ui為i方向流體運動速度;Uj為j方向流體運動速度;xi為i方向的坐標，xj為j方向的坐標，i，j=1，2，3，表示坐標分量[12];ρ為空氣密度;v為粘性系數;P為壓強;-u′iu′j為雷諾應力。

風雨環境下汽車外部流場的氣液兩相流動屬于湍流，本文湍流模型選用SST kω模型，該模型作為一種雙方程模型，在近壁面區域，采用kω模型對流場進行求解;在遠壁面區域，采用kε模型對流場進行求解[13]，而SST kω模型結合二者優勢，取長補短，具有適應性強、計算精度高、計算效率高等優點[14]，在實際工程中的應用越來越廣泛。SST kω湍流模型控制方程為[15]

kt+uikxi=1ρpk-βkω+1ρxiμi+σkμikxi（3）

ωt+uiωxi=1ρp∞-βω2+1ρxiμi+σ∞μiωxi+21-G1σ∞21ωkxiωxi（4）

式中，k為湍動能;t為時間;ui為流動速度;xi為x軸方向速度分量;ρ為空氣密度;pk為生成項;ω為比耗散率;β為封閉系數;μi為湍流粘度;σk為常數項;G1為混合函數;p∞為無窮遠處空氣壓強;σ∞和σ∞2為常數項[15]。

1.2 離散相模型

降雨環境下，雨滴在空氣中所占比例較小，屬于離散相，因此利用可形變部件模型（discrete phase model，DPM）進行數值計算。DPM模型的力平衡方程為[16]

dupdt=FDu-up+gxρp-ρρp（5）

式中，u為連續相空氣速度;up為離散相雨滴速度;ρp為離散相雨滴密度;gx為x方向重力加速度分量;FDu-up為單位質量雨滴阻力[16]。氣動阻力為[16]

FD=3CDu-up4dp（6）

式中，dp為雨滴直徑;CD為空氣阻力系數[16]。

雨滴粒徑為[17]

D=AIpn-1n1n（7）

式中，常數A=1.3;I為降雨強度常數，p=0.232;常數n=2.25[17]。

雨滴末速度為[18]

V=9.581-exp-D1.771.147（8）

2 數值計算模型

2.1 幾何模型及計算區域

幾何模型選用標準階梯背式MIRA汽車模型，該模型在汽車外流場分析領域得到廣泛應用，其幾何模型尺寸如圖1所示。模型車長L為4.165 m，車寬W為1.625 m，車高H為1.421 m，正面迎風面積為1.858 m2，側面迎風面積為3.827 m2。流場的計算區域如圖2所示。計算域長度為15 L（15倍車長），寬度為10 W（10倍車寬），高度為6 H（6倍車高），在車頂上方2 m處設置降雨平面，經驗證該高度對連續相空氣和離散相降雨均無干擾，在無偏航角工況下，降雨平面尺寸為30 m×10 m，在有偏航角工況下，降雨平面左側寬度有所增加，尺寸為30 m×13 m。

設在連續相空氣邊界條件，在無偏航角工況下，車身為壁面邊界，前端面為速度入口，后端面為壓力出口，底面為壁面邊界，左右兩側面均為對稱邊界，頂面為對稱邊界;有偏航角工況下，前端面和左側面為速度入口，后端面和右側面為壓力出口，頂面為對稱邊界。離散相降雨的邊界條件車身設為wallfilm邊界，應用stantonrutland碰撞飛濺模型，速度入口和壓力出口均設為escape逃逸邊界，底面設置為trap捕集邊界。

2.2 網格劃分及驗證

對計算域進行六面體網格劃分，計算網格劃分如圖3所示。對于無偏航角工況，計算域體網格尺寸設置為256 mm，并設置4個體網格加密區域，其網格尺寸由外到內分別為128，64，32，16 mm，車身模型主體面網格尺寸設置為16 mm，輪胎面網格尺寸設置為8 mm，在車身外表面劃分20層邊界層網格，并且近壁面第1層邊界層網格滿足y+<1，網格數為1 287萬。在有偏航角工況下，體網格加密區域寬度增加，其余網格參數與無偏航工況網格一致，網格數為1 420萬。

由于車速在60～100 km/h范圍內處于雷諾數自模區，即雷諾數只與車身外形有關，且不隨車速變化而變化，因此以80 km/h工況下的流場流動情況為代表。在車速為80 km/h且無偏航角的條件下，流場穩定后，測得的整車阻力系數為0.304，參考同濟大學TJ2風洞實驗數據為0.296[19]，相對誤差為2.70%，參考斯圖加特大學IVK風洞實驗數據為0.305[20]，相對誤差為0.33%，滿足工程應用要求。在車速為80 km/h且不同偏航角工況下，將數值模擬結果與同濟大學TJ2風洞實驗數據進行對比[19]，氣動力系數對比如圖4所示。由圖4可以看出，阻力系數和側力系數均符合工程應用要求。

3 數值計算結果及分析

3.1 汽車氣動特性分析

無降雨環境下，行駛的汽車車身表面壓力分布和外流場速度分布如圖5所示。由圖5可以看出，車頭、引擎蓋尾端和前擋風玻璃下端區域為正壓區，車身表面壓力較大，車速為60 km/h時，壓力值不小于100 Pa，車頭邊緣和前后擋風玻璃邊緣區域為負壓區，車身表面壓力較小，壓力值在-300 Pa左右，車頂和側玻璃區域整體為負壓區;隨著車速增加，正壓區和負壓區覆蓋范圍均增大，車頂和側面玻璃區域壓力略減小，車頭邊緣、車頂和車底附近空氣流速較高。車速為100 km/h時，流速超過30 m/s，車頭前端和車尾附近空氣流速較低;隨著車速增加，高流速區域空氣流速增大，覆蓋范圍變大，低流速區域變化不明顯。

降雨環境下，有、無偏航角工況下雨滴路徑變化如圖6所示。由圖6可以看出，降雨可以覆蓋整個車身，雨滴路徑分布較均勻，說明計算域和降雨平面的取值合理;雨滴在重力和風場的共同作用下運動，雨滴下落與車身外表面接觸，部分雨滴附著在車身表面，并在風場和重力作用下流動發生明顯的反彈，其中車頭、引擎蓋、前擋風玻璃、車身左側等迎風區域反彈現象較明顯，反彈的雨滴在風場的作用下繼續運動。

3.2 強降雨對氣動特性的影響

為研究強降雨對汽車空氣動力學特性的影響，在車速為60，80，100 km/h，降雨強度為100，300，500 mm/h工況下，進行汽車多相流數值計算，不同降雨強度下的氣動參數變化曲線如圖7所示。由圖7可以看出，相同降雨強度下，隨著車速的提高，汽車氣動阻力逐漸增加，最大達180%，而氣動阻力系數基本維持在同一水平，說明車速對氣動阻力影響較大，而對氣動阻力系數沒有實質性影響;車速相同的情況下，隨著降雨強度的增大，汽車氣動阻力和氣動阻力系數相對無降雨工況增加，但增加幅度相對較小，在60，80，100 km/h車速下，氣動阻力系數分別增加0.894%～1.987%，0.526%～1.500%，0.525%～1.535%，說明降雨強度越大，氣動阻力系數增幅越大，且氣動阻力系數與降雨強度近似呈線性關系。

車身表面壓力是與汽車設計相關的重要參數。當車速為60 km/h，降雨強度分別為0，100，300，500 mm/h工況下，不同降雨強度對車身表面壓力分布的影響如圖8所示。由圖8可以看出，隨著降雨強度的增加，車頂和尾箱蓋區域壓力增大，負壓區域面積減少，由于車頭、引擎蓋、前擋風玻璃區域為迎風面，其壓力受風場影響較大，故降雨對該區域表面壓力的影響不明顯。

為了直觀呈現雨滴在車身附近的分布和運動情況，在車速為60 km/h，降雨強度分別為100，300，500 mm/h的條件下，不同降雨強度下車身表面附近流場雨滴質量濃度分布如圖9所示。由圖9可以看出，當降雨強度為100 mm/h時，前擋風玻璃區域液相分布較明顯，雨滴質量濃度約為0.16 kg/m3，其余空間內的雨滴質量濃度相對較小;當降雨強度達到300 mm/h時，車頭、引擎蓋區域開始出現較明顯的液相分布，前擋風玻璃、車頂、尾箱蓋區域出現液滴飛濺現象;當降雨強度為500 mm/h時，液相分布更加明顯。隨著降雨強度的增加，汽車車身附近區域的雨滴質量濃度增大，最大雨滴質量濃度出現在前擋風玻璃附近區域。

3.3 偏航角對氣動特性的影響

為研究偏航角對汽車空氣動力學特性的影響，對車速為80 km/h，降雨強度為100 mm/h，偏航角為3°，9°，15°工況下的汽車空氣動力學數值進行計算，不同偏航角下汽車氣動參數變化曲線如圖10所示。由圖10可以看出，隨著偏航角的增大，氣動阻力系數和氣動側力系數隨之增大;在降雨強度為100 mm/h下，整車氣動阻力系數和氣動側力系數在3°，9°，15°偏航角工況下與無降雨時相比，增加幅度分別為2.57%，2.37%，1.70%，4.19%，3.65%和1.42%。研究結果表明，在無降雨條件下，小偏航角工況氣動力系數的基值較小，從而其數值增加幅度相對較大。

當車速為80 km/h，降雨強度為100 mm/h，不同偏航角下車身壓力分布側視圖和流場速度分布正視圖如圖11所示。由圖11可以看出，隨著偏航角增大，車身左側所受側風增大，左側面表面壓力逐漸增大，汽車后側窗邊緣區域氣體流速增大，該區域表面壓力減小;隨著偏航角增大，汽車左側所受側風增大，車頂高流速區域覆蓋范圍增大，并向車身右側延伸，車身右側低流速區域覆蓋范圍變大。

當車速為80 km/h，偏航角為15°時，有無降雨工況下車身表面壓力分布云圖如圖12所示。由圖12可以看出，降雨使汽車左側迎側風區域的表面壓力增大，尾箱蓋區域的壓力分布發生改變，其余區域壓力變化較不明顯。

4 結束語

本文采用ANSYS Fluent軟件，分析了MIRA階梯背汽車在不同降雨強度、不同偏航角下對車身氣動力、速度特性、壓力特性的影響，研究了車身表面附近區域雨滴質量濃度分布和變化。研究結果表明，在無降雨環境下，車頭、引擎蓋尾端和前擋風玻璃下端區域為正壓區，車頭邊緣和前后擋風玻璃邊緣區域為負壓區;隨著車速的增加，車身表面正壓區和負壓區覆蓋范圍均明顯增大;車頭邊緣、車頂和車底附近空氣流速較高，車頭前端和車尾附近空氣流速較低;隨著車速增加，高流速區域空氣流速增大，覆蓋范圍變大。相同降雨強度條件下，隨著車速的增加，汽車氣動阻力逐漸增加;車速對氣動阻力影響較大，而對氣動阻力系數沒有實質性影響;相同車速條件下，隨著降雨強度的增加，汽車氣動阻力和氣動阻力系數相對無降雨工況增加，降雨強度越大，氣動阻力（系數）增幅越大，且氣動阻力（系數）與降雨強度近似呈線性關系;隨著降雨強度的增加，汽車車身附近區域的雨滴濃度增大，最大雨滴質量濃度出現在前擋風玻璃附近區域。隨著偏航角的增大，氣動阻力（系數）和氣動側力（系數）隨之增大，左側車身表面壓力逐漸增大，車頂高流速區域覆蓋范圍增大，并向車身右側延伸，車身右側低流速區域覆蓋范圍變大;降雨使車身左側迎側風區域的表面壓力增大，尾箱蓋區域的壓力分布發生改變。該研究為汽車動力學設計和整車研發提供了參考。

參考文獻：

[1] 王博宇， 李杰偉. 中國交通事故的統計分析及對策[J]. 當代經濟， 2015（20）： 116119.

[2] 李大軍. 汽車運輸的安全問題分析[J]. 時代汽車， 2020（22）： 191192.

[3] 公安部交通管理局. 中華人民共和國道路交通事故統計年報（2010年度）[R]. 北京： 公安部交通管理局， 2011.

[4] 李長城， 劉小明， 榮建. 降雨條件下高速公路車輛行駛速度特性[J]. 北京工業大學學報， 2015， 41（3）： 412418.

[5] 馮偉， 杜禮明， 倪守隆. 風雨環境下轎車行駛中的氣動特性研究[J]. 機電工程， 2015， 32（4）： 474478.

[6] 王保成. 降雨對高速公路行車安全的影響[J]. 中國水運（下半月）， 2015， 15（7）： 229230.

[7] 馬筱櫟， 樊博. 降雨對山區高速公路運行車速的影響研究[J]. 公路與汽運， 2020（5）： 3336.

[8] 唐文淵， 朱泓澄， 單鴻煊， 等. 汽車車身氣動造型數值仿真與優化設計[J]. 常熟理工學院學報， 2019， 33（2）： 6366.

[9] 劉旺， 侯路飛， 顧君杰， 等. 尾翼對汽車氣動性能的影響[J]. 汽車工程師， 2020（6）： 5657.

[10] 王毅剛， 楊超， 楊志剛， 等. 汽車外表面氣動噪聲特性分析[J]. 聲學技術， 2014， 33（1）： 5055.

[11] 方剛， 吳軍， 羅煜峰， 等. 汽車三維流場氣動網格的應用研究[J]. 湖南大學學報： 自然科學版， 2001， 28（3）： 3236.

[12] 梅元貴， 許建林， 趙鶴群， 等. 側風環境下高速列車外部流場數值模擬方法研究[J]. 工程力學， 2012， 29（6）： 253258， 278.

[13] 姬永魁. H型垂直軸風力機風輪氣動設計及氣動性能研究[D]. 蘭州： 蘭州理工大學， 2012.

[14] 文恒， 于夢閣， 盛旭高， 等. 沙塵環境下高速列車氣動特性分析[J]. 青島大學學報： 工程技術版， 2020， 35（4）： 5054.

[15] 楊金廣， 吳虎. 雙方程kω SST湍流模型的顯式耦合求解及其在葉輪機械中的應用[J]. 航空學報， 2014， 35（1）： 116124.

[16] Li H Q， Yu M G， Zhang Q， et al. A numerical study of the aerodynamic characteristics of a highspeed train under the effect of crosswind and rain[J]. Fluid Dynamics and Materials Processing， 2020， 15（4）： 7790.

[17] Morsi S A， Alexander A J. An investigation of particle trajectories in twophase flow systems[J]. Journal of Fluid Mechanics， 1972， 55（2）： 193208.

[18] Bset A C. Empirical formulae for the terminal velocity of water drops falling through the atmosphere[J]. Quarterly Journal of the Royal Meteorological Society， 1950， 76（329）： 302311.

[19] 周華， 楊志剛， 朱暉. 基于整車風洞試驗的MIRA車型數值計算[J]. 吉林大學學報： 工學版， 2019， 49（4）： 10431053.

[20] 龐加斌， 林志興， 余卓平， 等. TJ2風洞汽車模型試驗的修正方法[J]. 汽車工程， 2002， 24（5）： 371375.

Analysis of Vehicle Aerodynamic Characteristics under Heavy Rainfall Environment

YIN Shuo， YU Mengge， CHEN Huanming， WANG Yalun

（College of Mechanical and Electrical Engineering， Qingdao University， Qingdao 266071， China）

Abstract：? In order to study the aerodynamic characteristics of vehicles under heavy rainfall environment， the standard MIRA stepped back vehicle model is taken as the research object， and the aerodynamic numerical calculation model of vehicle outflow field is established. ANSYS fluent is used to analyze the influence of different rainfall intensity and yaw angle on the aerodynamic， velocity and pressure characteristics of car body， track the movement track of discrete phase raindrop and its contact with car body， study the distribution and change of raindrop concentration near the car body surface， and verify the accuracy of the model. The results show that： in the with rainfall environment， with the increase of vehicle speed， the coverage of positive pressure area and negative pressure area on the vehicle body surface increases obviously， and some raindrops rebound on the vehicle body surface obviously; Vehicle speed has great influence on aerodynamic drag， but has no substantial influence on aerodynamic drag coefficient; Under the same vehicle speed， the larger the rainfall intensity is， the larger the aerodynamic drag （coefficient） increases， and the relationship between aerodynamic drag （coefficient） and rainfall intensity is approximately linear; Under a certain rainfall intensity， with the increase of yaw angle， the aerodynamic resistance （coefficient） and aerodynamic side force （coefficient） of the whole vehicle increase compared with that without rain， the surface pressure of the left side gradually increases， the pressure of the rear side window edge area decreases， the coverage of the high velocity area on the roof increases， and extends to the right side of the vehicle body， and the coverage of the low velocity area on the right side of the vehicle body increases; Rainfall increases the surface pressure in the left side of the car， and changes the pressure distribution in the tailbox cover area. This study provides a reference for the research and application of vehicle aerodynamics in heavy rainfall environment.

Key words： numerical calculation; vehicle aerodynamic characteristics; speed; rainfall intensity; yaw angle