王國俊　崔文詩　周華　朱暉　楊志剛

摘要： 為有效開展汽車氣動優(yōu)化設計和主動流動控制技術，采用正交分解（Proper Orthogonal Decomposition， POD）技術分析汽車非定常流場中POD模態(tài)的物理意義，并利用重構技術建立汽車外流場瞬態(tài)簡化模型.大渦模擬（Large Eddy Simulation， LES）仿真結果與實驗結果的對比表明該仿真方法有效.POD分析結果表明：低階模態(tài)主要對應分離泡的附著流動，包含流場中較高的能量；其他各個高階模態(tài)含能相對較低，表征小尺度渦的作用.當保留前15階POD模態(tài)時，通過重構技術可獲得70%的原流場能量，得到的降階模型可以較好地捕捉流動的特征.

關鍵詞： 汽車； 非定常流場； 正交分解； 大渦模擬； 簡化模型； Ahmed body； 流場重構； 渦能量耗散

中圖分類號： U463.51； TB115.1文獻標志碼： B

Abstract： To carry out the automobile aerodynamic optimization design and the active flow-control technology effectively， the Proper Orthogonal Decomposition （POD） technology is used to analyze the physical meaning of POD modes in unsteady flow field around automobile， and a simplified transient model is built for flow field of automobile outer flow field. The Large Eddy Simulation（LES） results are compared with the experimental results to verify that the simulation method is effective. The POD analysis results show that， the low-order modes correspond to the attached flow of separation bubbles， which contains high energy of the flow field. Other high-order modes contain relatively low energy， which represents the effect of small scale vortex. If the first 15 order POD modes are retained， 70% energy of the original flow field can be obtained through reconstuction technology， the low-dimensional model can well capture the characteristics of the flow.

Key words： automobile； unsteady flow field； proper orthogonal decomposition； large eddy simulation； simplified model； Ahmed body； flow field reconstruction； vortex energy dissipation

0引言

汽車外形復雜，其運動流場是一個高階非線性的復雜系統(tǒng)，存在著三維非定常分離流和近地面效應等復雜流動問題.汽車非定常流場對行駛穩(wěn)定性、氣動阻力、尾氣傳播和污染物夾帶等[1]有重要影響，因此有必要將復雜流體動力系統(tǒng)分解并簡化，得到主要流動結構，然后在此基礎上進行高精度的流場重構，建立汽車非定常流場簡化模型，進而實現高效的汽車氣動優(yōu)化設計和主動流動控制.

近年來，本征正交分解（Proper Orthogonal Decomposition，POD）技術被廣泛應用于解決實驗和數值模擬過程中的降維問題，如空氣動力學、最優(yōu)控制、結構動力學和熱傳導分析[2]等.POD技術是一種功能強大的降維和數據分解方法[3]，可在給出較高逼近度的情況下，有效降低物理問題的自由度，建立近似模型.因此，POD技術為高效分解和預測汽車非定常流場提供一種有效的工具.STH等[4]利用大渦模擬（Large Eddy Simulation， LES）對尾部傾角為90°的車體進行y向和x向速度計算，利用POD開展車體展向兩側渦脫落模式的研究.斜背后風窗對整車氣動阻力和噪聲影響較大[5]，但以往研究較少.國內外尚未開展基于渦量、壓力等參數直接表征車輛氣動性能的POD研究.

本文首先將數值仿真與實驗結果比較，驗證LES計算的有效性；然后，通過LES得到一系列快照的樣本后，基于POD技術得到不同模態(tài)的擬序結構；最后，基于主要的流動結構建立斜背后風窗非定常簡化流場（x和z向速度、渦量和壓力），并與LES結果進行對比驗證.

1數值仿真結果驗證

1.1模型與仿真方法

以尾部傾角為25°的Ahmed body模型[5]為研究對象，采用0.7倍縮比的Ahmed body模型計算.縮比模型長L=0.731 m，寬W=0.272 m，高H=0.202 m，斜背后風窗長Lr=0.155 m，來流風速為25 m/s，雷諾數為1.2×106.根據模型尺寸，確定計算域長為8L，寬為7W，高為4H.為保證仿真準確性，設置第一層網格尺寸y+小于1，網格尺寸為0.3～1.5 mm，采用六面體網格對計算域空間進行離散，在流場變化較為劇烈的斜背和尾部進行網格加密以提高計算精度，體網格總數約為2 444萬個，見圖1.使用商業(yè)軟件FLUENT 12.1進行LES計算，時間步長取0.000 1 s，迭代方法為SIMPLIC算法.

1.2仿真結果驗證

尾部平均壓力系數仿真結果見圖2.在z=0.17 m以下（車體尾部垂直背）和0.17 m

模型車身附近的幾個渦脫落頻率監(jiān)測點示意見圖3.對比實驗與仿真結果的斯特勞哈爾數Sr=Hf0U-10，其中H為模型高度，f0為自然渦脫落頻率.實驗獲得的Sr=0.2，LES得到的Sr=0.22，誤差較小，因此可以認為此次仿真是有效的.

通過風洞實驗，AHMED等[5]總結斜背傾角為25°的Ahmed body模型的氣動阻力來源，認為阻力主要源于斜背的負壓，而垂直背對阻力貢獻較小.外流沿著車頂壁面向下游流動（見圖4），氣流附著在斜背上，實驗中的氣流附著位置更靠近斜背下緣，LES同樣可捕捉到類似的流動現象.

2POD分析

2.1POD簡介及其理論

POD方法于1967年首先引入到湍流流動結構研究中.[8]該方法的基本思想是把原來在時間域和空間域上的連續(xù)物理量場分解成只與時間和空間相關的函數展開序列，且其在均方意義上最優(yōu)，在展開式中只需要少量的項數就可較準確地描述該物理過程，從而可以提供具有足夠精度而自由度又少的低維模型，以降低計算維數、節(jié)省計算時間和內存.

A，B，C和D的坐標為：A（1.64，0.236）B（0.171，1.781）

C（1.807，0.228）D（1.664，0.294）

2.3POD分析結果

分別基于速度、渦量和壓力的POD分析得到的模態(tài)數量皆為500，隨著模態(tài)階數的增多，含能降低.低階模態(tài)通常代表流場中具有較高動能的大尺度流動結構，是整個流動的載體.1階模態(tài)通常在流場發(fā)展中占支配地位，其含能最高，速度POD的1階模態(tài)占總動能的16%，渦量POD的1階模態(tài)占總動能的15%，壓力POD的1階模態(tài)占總動能的15%，皆遠高于其對應的2階模態(tài).各變量的2～15階模態(tài)含能占比皆超過1%，認為這些模態(tài)對流動仍有重要影響.速度、渦量和壓力的1階模態(tài)含能相差不到1%，2～7階模態(tài)能量占比差異不大，8～15階模態(tài)能量占比存在較小差異.速度、渦量和壓力的前100階模態(tài)各本征值占總動能的比例見圖5.由此可知：不同變量的15階模態(tài)含能總和皆在70%左右；隨著階數的增加，43階模態(tài)后的每個模態(tài)能量占比皆趨近于0.由于高階模態(tài)對總動能的貢獻低，對應的是小尺度結構的作用，對流場的影響較低.每階模態(tài)代表流動中不同尺度的流動結構.0階模態(tài)代表時均速度場.速度模態(tài)表征脈動速度的分布，見圖6.由此可知：1階模態(tài)為速度場中的主要含能結構，包含流場的主要流動信息，其中x向速度模態(tài)更為清楚地展現上述流動現象；x和z向速度的2階和3階模態(tài)中高速脈動和低速脈動集中區(qū)呈現交替分布狀態(tài)，z向速度模態(tài)中的交替結構更加明顯.定義兩高速脈動集中區(qū)（或低速脈動集中區(qū)）之間的距離為波長λ，通過圖6e和6f可以清楚的觀察到相鄰模態(tài)之間高低速區(qū)以λ/4進行交替并向下游運動，文獻[10]認為此現象與大尺度流動結構的耗散有關.

渦量的1階模態(tài)為流場中的主要含能結構，其1階模態(tài)見圖7.由此可知：1階模態(tài)對應的渦量分布形態(tài)對應斜背回流區(qū).應該指出的是，高強度的渦量本身意味著高能量的消耗.y=0面渦量云圖見圖8.由圖8可以看出：在斜背分離線處渦量大，隨著流動往下游進行，這些漩渦逐漸分離，由能量較高的渦分離出能量較小的渦，在到x=1.8左右的時候，漩渦逐漸破碎，能量進一步耗散.1階模態(tài)圖直觀地給出流場中能量消耗的主要來源.

基于靜壓的1階模態(tài)也是流場中的主要含能結構.在圖9a中，1階模態(tài)中的負壓區(qū)對應斜背回流區(qū)域，即為壓力脈動在斜背附近空間的分布規(guī)律.在圖9b和9c中，壓力的2階和3階模態(tài)的高、低壓力脈動集中區(qū)呈現交替分布，并向下游運動，這與基于z向速度的2階和3階模態(tài)結果類似，同樣與大尺度流動結構的耗散有關.

3流場重構與驗證

對含能貢獻較低的高階模態(tài)進行截斷.根據對流場能量的需求，認為對包含70%能量的前幾階模態(tài)進行重構，可獲得對流場較好的重構結果，即可簡化流場.利用基于速度的前15階模態(tài)進行重構，重構流場包含流場中70%的能量.速度場重構與驗證見圖10.

4結論

通過與實驗結果相比，驗證LES的有效性，進一步分析得到以下結論.

（1）相對于直接對非定常復雜流場進行分析而言，采用POD技術可以更加明確地辨別主導流動現象的主要流動結構.分別基于x和z向速度與渦量的POD分析指出：1階模態(tài)含能最高，主要對應回流區(qū)的附著流動，是脈動能量的主要耗散來源；高階模態(tài)含能相對較低，表征小尺度渦的耗散作用.

（2）在重構速度場中，當保留前15階POD模態(tài)時，通過重構技術可獲得70%的原場能量，得到的降階模型能很好地重現斜背分離流中復雜的大小附著渦結構以及二次分離渦結構等非定常流動現象.

（3）基于渦量和壓力POD的分析分別給出能量耗散的主要來源和壓力脈動分布，可作為評價參數指導工程應用.參考文獻：

[1]CHOI H， LEE J， PARK H. Aerodynamics of heavy vehicles[J]. Annual Review of Fluid Mechanics， 2014， 46： 441-68. DOI： 10.1146/annurev-fluid-011212-140616.

[2]SIROVICH L. Turbulence and the dynamics of coherent structures part I： coherent structures[J]. Quarterly of Applied Mathematics， 1987， 45（3）： 561-571.

[3]胡金秀， 鄭保敬， 高效偉. 基于特征正交分解降價模型的瞬態(tài)熱傳導分析[J]. 中國科學： 物理學 力學 天文學[J]. 2015， 45（1）： 73-84.

HU J X， ZHENG B J， GAO X W. Reduced order model analysis method via proper orthogonal decomposition for transient heat conduction[J]. Scientia Sinica（Physica， Mechanica & Astronomica）， 2015， 45（1）： 73-84.

[4]STH J， NOACK B R， KRAJNOVIC′ S， et al. On the need for a nonlinear subscale turbulence a term in POD models as exemplified for a high-Reynolds-number flow over an Ahmed body[J]. Journal of Fluid Mechanics， 2014， 747： 518-544. DOI： 10.1017/jfm.2014.168.

[5]AHMED S R， RAMM R， FALTIN G. Some salient features of the time-averaged ground vehicle wake[DB/OL]. （1984-02-01）[2016-04-11]. http：//papers.sae.org/840300/. DOI： 10.4271/840300.

[6]KOURTA A， LECLERC C. Characterization of synthetic jet actuation with application to Ahmed body wake[J]. Sensors and Actuators A： Physical， 2013， 192： 13-26. DOI： 10.1016/j.sna.2012.12.008.

[7]MINGUEZ M， PASQUETTI R， SERRE E. High-order large-eddy simulation of flow over the “Ahmed body” car model[J]. Physics of Fluids ， 2008， 20（9）： 095101. DOI： 10.1063/1.2952595.

[8]BEN CHIEKH M， MICHARD M， GUELLOUZ M S， et al. POD analysis of momentumless trailing edge wake using synthetic jet actuation[J]. Experimental Thermal and Fluid Science， 2012， 46： 89-102. DOI： 10.1016/j.expthermflusci.2012.11.024.

[9]王掩剛， 陳俊旭， 先松川. 基于POD方法的二維方柱低雷諾數繞流流場分析研究[J]. 西北工業(yè)大學學報， 2014， 32（4）： 612-617.

WANG Y G， CHEN J X， XIAN S C. Analysis of low Reynolds number flow around a square cylinder based on POD method[J]. Journal of Northwestern Polytechnical University， 2014， 32（4）： 612-617.

[10]THACKER A， AUBRUN S， LEROY A， et al. Experimental characterization of flow unsteadiness in the centerline plane of an Ahmed body rear slant[DB/OL]. （2013-02-20）[2016-04-11]. http：//link.springer.com/article/10.1007/s00348-013-1479-5？no-access=true. DOI： 10.1007/s00348-013-1479-5.

（編輯武曉英）