柴博森 王廣義 朱國仁? 閆東 陸振華 遲成芳

(1.吉林大學 機械與航空航天工程學院，吉林 長春 130022； 2.吉林大學 汽車仿真與控制國家重點實驗室，吉林 長春 130022； 3.長安大學 高速公路筑養(yǎng)裝備與技術教育部工程研究中心，陜西 西安 710061； 4.江蘇華隆興機械工程有限公司，江蘇 無錫 214023)

液力變矩器(TC)是自動變速汽車和工程機械領域廣泛應用的關鍵基礎部件[1]。目前高性能液力變矩器產品嚴重依賴國外進口，核心研發(fā)技術掌握在德國采埃孚(ZF)和日本愛信(Aisin)等國際知名企業(yè)中。液力變矩器先進設計理論與關鍵技術的缺失已成為制約我國自動變速汽車和工程機械領域發(fā)展的“卡脖子”問題，只有加強對液力變矩器基礎理論的研究和前瞻性技術的突破，才能滿足國家在自動變速汽車和工程機械領域發(fā)展的戰(zhàn)略需求，為國民經濟支柱產業(yè)汽車工業(yè)和工程機械的蓬勃發(fā)展強筋健骨。

液力變矩器是依靠流體傳動實現(xiàn)能量轉換的葉輪機械，其內部流動特性決定外部工作性能，深入研究其內部湍流時空演化規(guī)律與流動控制機理對于優(yōu)化設計其葉柵結構及性能改進具有重要的科學研究意義[2]。目前研究液力變矩器流場的主要手段是基于計算流體動力學(CFD)的數值模擬方法和采用非接觸測量方式的流場可視化技術[3]。在數值模擬方面，人們已經深刻地體會到發(fā)展高精度數值模擬方法對于精細刻畫湍流時空演化特征并準確地預測工作特性的重要性[4]?；诖鬁u模擬(LES)方法開展液力變矩器湍流高精度數值模擬已成為一項極具挑戰(zhàn)性的科研課題。國外，Odier等[5]提出了一種適用于復雜壁面模型大渦模擬的網格自適應方法，綜合考慮壁面邊界層y+值和平均動能耗散，預測了渦輪內部湍流特性。Yang等[6]提出一種新型的壁面建模方法，采用大渦模擬預測近壁面流場分布情況。日本北海道大學Tasaka等[7]采用大渦模擬仿真液力變矩器內部的三維流動，并進行了流場試驗對比分析。Fuente等[8]采用k-ε湍流模型仿真了液力變矩器在瞬態(tài)和穩(wěn)態(tài)時的流場特性，但該兩方程湍流模型仿真精度不足。國內，同濟大學王立軍等[9]基于大渦模擬技術，采用格子Boltzmann方法數值模擬了液力變矩器導輪流場，重現(xiàn)了導輪葉片尾跡區(qū)復雜的流動細節(jié)。北京理工大學劉城等[10]基于大渦模擬應力混合渦模型仿真了高功率密度液力變矩器內流場，揭示了流場能量損失機理。哈爾濱工業(yè)大學杜明杰采用大渦模擬方法仿真了渦輪內部流場，分析了葉片表面凹坑結構對壁面邊界層流場的影響規(guī)律[11]。在流動可視化試驗測量方面，激光多普勒測速(LDV)和粒子圖像測速(PIV)逐漸成為液力變矩器湍流實驗測量研究的主流技術。國外，美國亞拉馬巴州大學Lee等[12]基于LDV技術揭示了液力變矩器內部流動特性，研究了泵輪、渦輪之間轉速差對其內部復雜流動效應的影響規(guī)律。弗吉尼亞大學Flack等[13]采用激光測速技術測量了液力變矩器渦輪和導輪內部流場并分析了不同轉速比下流速的變化。韓國機械技術學會Yoo等[14]基于LDV技術對泵輪和渦輪葉片之間的流場區(qū)域進行測量，發(fā)現(xiàn)渦輪的瞬時位置對葉輪通道間隙區(qū)域內的流體流動有較大影響，且單點測量具有一定的局限性。日本福井大學Yamamoto等[15]基于PIV技術實現(xiàn)了液力變矩器流場可視化，提出了二進制圖像互相關算法進行流場計算。國內，北京理工大學祝自來等[16]基于LDV技術測量了液力變矩器流場，詳細分析了示蹤粒子直徑、濃度對其跟隨性的影響規(guī)律。北京理工大學李晉等[17]利用LDV技術測量了液力變矩器泵輪流場，發(fā)現(xiàn)導輪對泵輪流道上游的影響強于渦輪，而渦輪主要影響泵輪出口處的低流速區(qū)域。吉林大學柴博森等[18]基于PIV技術采集了液力變矩器渦輪流場圖像，發(fā)現(xiàn)提高轉速比后，二次流、反向流等復雜流動現(xiàn)象逐漸降低，流動能量損耗降低。吉林大學柴博森等[19]基于二維切面流場試驗測量系統(tǒng)采集了制動工況下液力偶合器徑向切面漩渦流場圖像，經圖像后處理識別剔除錯誤流速矢量，優(yōu)化了流動圖譜。綜上所述，近些年研究者越來越關注液力變矩器內部流場數值模擬與試驗測量，但是在流場高精度仿真與流動可視化精確測量等方面仍然存在很多亟待解決的科學問題。尤其是針對大渦模擬流場仿真，在湍流模型的合理選擇及其計算方法的正確應用等方面需要開展更深入的研究。

為了準確地揭示液力變矩器流場時空演化機理，基于大渦模擬技術采用5種亞格子湍流模型精細化仿真液力變矩器流場，對比分析流場仿真結果的差異性。通過PIV流場可視化試驗驗證數值模擬結果的可靠性，并對5種亞格子湍流模型的仿真適用性進行分析。

1 研究對象

1.1 計算模型

建立液力變矩器三維模型，其循環(huán)圓有效直徑為315 mm，泵輪葉片數28個，渦輪葉片數27個，導輪葉片數17個，如圖1(a)所示。為了簡化計算域，忽略各葉輪間的流體泄露，認為液力變矩器3個葉輪組成了封閉的耦合流動域，經過布爾運算提取全流道幾何模型作為流場計算域，如圖1(b)所示。

圖1 液力變矩器計算模型

1.2 網格模型

基于ICEM軟件采用六面體結構化網格對液力變矩器全流道計算模型進行精細化網格劃分。由于葉片近壁面邊界層流動情況對整體流道內部流場演化過程影響顯著，為了精確化求解葉片近壁面流場結構，對葉片近壁面區(qū)域進行局部網格加密處理，如圖2所示?？紤]網格無關性因素[20]，調整最終確定全流道計算域總體網格數為3427543。

圖2 全流道網格及葉片近壁面網格加密處理

2 數值計算方法

2.1 控制方程

采用大渦模擬(LES)方法數值計算液力變矩器流場。LES的控制方程是濾波函數對N-S方程進行濾波得到，基于網格尺寸篩選，比網格尺寸小的渦被忽略,大渦直接解析，小渦被?；?，濾波函數計算公式如下：

G(|x-x′|)=

(1)

式中，G為濾波函數，xi為張量形式的指標，Δxi為坐標系中沿i方向上的網格尺度。

經過濾波的大渦模擬控制方程如下：

(2)

式中，ρ為流體密度，t為時間，ui和uj為張量形式的時均速度，xj為張量形式的指標符號,p為壓強，μ為流體黏度，帶有上劃線的量為濾波后的場變量，ij表示濾波后得到的大小尺度脈動之間動量交換的亞格子應力項，可通過下式求解：

(3)

式中，Sij為應力張量比率，kk為各向同性的亞格子應力，μt為亞格子尺度的湍流黏度。

2.2 亞格子模型

亞格子湍流渦黏度μt定義不同則構成不同亞格子模型，文中采用渦黏性模型(SL)、局部渦黏度壁面自適應模型(WALE)、代數形式的壁面函數模型(WMLES)、改進的代數形式的壁面函數模型(WMLES S-Omega)和動態(tài)動能亞格子模型(KET)5種亞格子模型來數值模擬液力變矩器流場。

SL模型[21]滿足控制方程如下：

(4)

亞格子渦黏度由下式定義：

(5)

式中:δij為克羅內克函數,Smagorinsky系數Cs常取0.1；網格尺寸V決定濾波尺寸Δ，Δ=V1/3。

WALE模型[22]湍流渦黏度表達如下：

(6)

(7)

WMLES模型[23]是為了克服雷諾數的尺度限制而提出的，湍流渦黏度μt定義為

μt=(CSmagΔ)2S{1-exp[-(y+/25)3]}

(8)

式中，S為應變率,常數CSmag=0.2，y+為第一層網格質心到壁面的無量綱距離。

由于WMLES模型沒有考慮恒定剪切流中的零渦黏性，WMLES S-Omega模型被提出來解決恒定剪切流中的零渦黏性現(xiàn)象建模，在這個模型中只是將式(8)中的S替換成|S-Ω|，其中Ω為旋渦強度。

KET模型考慮了亞格子尺度湍動能的模擬，湍流渦黏度μt可通過以下方程組求得：

式中，Ck為模型系數，ksgs為湍動能，qsgs為亞格子湍動能。

2.3 計算方法

以純凈水作為液力變矩器工作介質，密度ρ=998.2 kg/m3，黏度μ=0.001 003 Pa·s，忽略流動介質溫度變化對密度和粘度的影響及流動泄漏的影響，設置各葉輪交界面為滑移網格，壁面處采用無滑移邊界條件?；贏NSYS-Fluent軟件采用大渦模擬亞格子湍流模型仿真制動工況(i=0)下液力變矩器流場，i為傳動比。泵輪輸入轉速為200 r/min，計算方法選取SIMPLE算法，壓力方程使用Standard格式離散，采用二階迎風格式求解動量方程與湍動能，設置瞬態(tài)求解時間步長0.01 s，迭代計算總步數為500，監(jiān)測殘差曲線與泵輪轉矩曲線作為判斷仿真結果是否收斂的依據。

3 流場仿真結果

3.1 三維渦識別與流場結構分析

流體力學科學家柯奇曼(Kucheman)認為旋渦是流體運動的肌腱[24]。液力變矩器工作時其內部充滿紊亂且不同尺度下的旋渦流動，多尺度旋渦流動的產生、發(fā)展以及彼此之間的相互作用支配著液力變矩器內部整體流動。準確地識別并提取液力變矩器三維渦旋結構特征，對于揭示其內部湍流時空演化規(guī)律及流動能量轉化與損耗機理具有重要的推動作用。在流場數值模擬的基礎上，利用CFD-Post軟件對流場仿真結果進行后處理，基于Q準則渦識別方法，選取合理閾值識別并提取液力變矩器內部整體三維渦旋結構，如圖3所示，其中渦輪葉片、泵輪葉片、導輪葉片分別采用紅色、白色和黑色區(qū)別表示，紫色箭頭表示從泵輪到渦輪的流動方向，黃色箭頭表示渦輪與導輪之間的流動趨勢。

圖3 三維渦時空結構

制動工況下，渦輪靜止不動，攜能高速流體由泵輪流道流出，沖擊渦輪葉片壓力面，由于渦輪葉片幾何扭曲程度大，導致沿葉片壁面處流體速度大小和方向急劇改變，形成多尺度渦旋結構，并由此驅動引起渦輪流道內部形成豐富多彩的不同結構的多尺度渦旋結構，如：通道渦(TD)、長條渦(CT)、發(fā)卡渦(FQ)、脫落渦(TL)等，如圖3所示。在離心力和葉片的雙重作用下，泵輪流道入口處小尺度分離渦經過旋轉、合并后逐漸形成較大尺度的混合渦結構，并加速沖擊進入渦輪流道內。由于渦輪葉片空間扭曲程度大，在扭曲葉片的作用下，逐漸在渦輪主流區(qū)域上形成明顯的通道渦，在這個過程中伴隨著多尺度渦結構的匯聚與合并，不斷強化通道渦結構的強度，并集聚大量能量，促成主流區(qū)域上的大尺度旋渦流動，考慮流體粘滯損耗，一部分流動動能轉化為流體內能，并最終以熱能形式耗散。隨著主流區(qū)域旋渦能量逐步向渦輪流道出口處傳遞，由于出口處流道空間變狹窄，流道出口處橫截面積變小，導致渦輪出口處流動加速，通道渦在拉伸、彎曲、扭轉、撕裂的作用下，逐漸由較大尺度渦旋結構轉變?yōu)樾〕叨认碌拈L條渦結構，如圖3渦輪出口處再次形成長條渦，且該處長條渦結構左側多呈現(xiàn)粗圓狀，右側多呈現(xiàn)尖細狀，整體結構形態(tài)類似于“子彈”形式，呈現(xiàn)加速狀態(tài)。在渦輪與導輪交界面無葉柵區(qū)域，出現(xiàn)部分小尺度脫落渦結構，伴隨一定的能量損耗。渦輪出口處的流體經過空間扭曲葉片的揚拋作用，以較大的入射角沖擊導輪葉片頭部，在導輪入口處靠近葉片吸力面位置形成發(fā)卡渦，并緊貼著葉片吸力面逐漸向后蔓延。為了揭示制動工況下渦輪流道內部非定常多尺度渦時空演化規(guī)律，基于Q準則渦識別方法精細提取5種亞格子湍流模型仿真結果中的多尺度渦旋結構，選取閾值55 842.1 s-2，提取結果如圖4所示。

圖4 五種湍流模型渦輪流道三維渦結構

圖4展示了5種亞格子模型渦輪流道瞬時三維渦空間分布形態(tài)。WMLES和WMLES S-Omega模型在葉片吸力面附近渦分布較密，流場情況更為復雜，而SL模型在該區(qū)域內渦結構識別不夠準確，沒有捕捉到較為詳細的流場狀態(tài)。主流區(qū)域上，WALE模型在該區(qū)域內渦結構破碎，這是由于攜能大尺度渦轉化為小尺度渦過程中，中間伴隨著渦的拉伸、撕裂、碰撞、合并等動力學作用，引起流場結構劇烈變化，引起能量傳遞損耗。SL模型沒有模擬出葉片壓力面處的小尺度渦旋結構，而其他模型均有體現(xiàn)。由于考慮亞格子尺度湍動能數值模擬，KET亞格子模型針對葉片尾部的脫落渦結構仿真結果相對于其他模型更加逼真。

3.2 瞬時流速場和渦量場分析

通過PIV試驗驗證仿真結果的準確性。截取數值模擬后的流速場和渦量場，如圖5所示，所截取的截面與PIV試驗測量中的二維截面保持一致。

(a)流速 (b)渦量

從流速場整體結構分布來看，5種亞格子湍流模型仿真結果周期性明顯，如圖5所示。從流速場局部分布定性來看，在渦輪流道入口處，外環(huán)與葉片吸力面的角隅區(qū)域存在明顯的高速區(qū)，這是由于來自泵輪流道內部的高速液流沖擊進入渦輪流道，經扭曲渦輪葉片的揚拋作用，進而加速流體運動而引起流速增大的結果。低速區(qū)主要集中在流道中間位置靠近吸力面處，這是由于液流沖擊渦輪壓力面導致流速大小和方向發(fā)生驟變，進而形成多尺度旋渦結構，且這些多尺度旋渦結構彼此之間發(fā)生沖擊、合并、分離相互作用，引起流動能量轉變和損耗。由于液流粘滯損耗的作用，部分損失的流體能量導致流體介質溫升，流體動能轉化為內能并最終以熱能形式耗散。在渦輪出口處由于流道變狹小，流速逐漸升高。

從渦量場整體分布來看，高渦量場主要集中在渦輪中間位置高、低流速交界處和近壁面區(qū)域附近。在高、低速局部流場交界處，流速梯度變化劇烈，引起大尺度旋渦流動，局部區(qū)域渦量值較大。小尺度渦旋主要分布在葉片近壁面，主流區(qū)域上大尺度旋渦流動引起的高速液流沖擊葉片表面，由于葉片幾何結構空間扭曲，在葉片近壁面區(qū)域形成較小尺度下的渦旋結構。主流區(qū)域大尺度渦旋攜帶的流動能量逐漸傳遞到小尺度渦旋，由于葉片壁面區(qū)域、角隅區(qū)域處所形成的小尺度渦運動方向與主流渦運動方向有所不同，方向相同的小尺度渦旋對于主流區(qū)域上整體渦旋運動起到促進作用，而方向相反的小尺度渦旋對于主流區(qū)域上整體渦旋運動起到阻礙作用，在流體摩擦的作用下，考慮流體的粘滯損耗作用導致流動能量損耗。

4 PIV試驗及數據處理

4.1 試驗測量系統(tǒng)

為了驗證仿真結果的準確性，基于粒子圖像測速技術設計并搭建流場可視化試驗臺。試驗臺主要由機械系統(tǒng)、光學系統(tǒng)和圖像采集系統(tǒng)組成，見圖6。機械系統(tǒng)主要由1.1 kW變頻調速電機、聯(lián)軸器、轉速測量儀及制動裝置組成，用于液力變矩器動力驅動和載荷施加，并記錄泵輪輸入軸轉速和渦輪輸出軸轉速。光學系統(tǒng)采用西安遠訊光電科技有限公司生產的型號為PT-532- 1.5 W-L的激光器，可提供最大功率1.5 W、波長532 nm、扇面張角為90度的激光片光，片光厚度可達2 mm，使用前先預熱，將液力變矩器待測流場區(qū)域照亮。圖像采集系統(tǒng)采用廣州市元奧儀器有限公司生產的型號為FR- 340- 10G圖像記錄系統(tǒng)及配套軟件，相機幀頻為338 fps，采集液力變矩器內部流動圖像。試驗樣機為循環(huán)圓有效直徑為315 mm的透明型液力變矩器，按實際液力變矩器大小以1:1比例制作，材質選擇為有機玻璃，試驗樣機的壁厚為2 mm，避免引入較大的光路折射誤差和圖像畸變等問題。

4.2 圖像采集與標定

以蒸餾水作為流動介質，選取直徑為1.5 μm的PSP示蹤粒子投入流場進行試驗測量，保證粒子濃度符合圖像采集要求[25]。通過變頻調速電機控制泵輪輸入轉速為200 r/min，合理布置高清CCD相機空間位置，保證相機光軸垂直激光平面，可采集清晰的渦輪流道內部二維切面流場視頻，經圖像分幀處理選取連續(xù)兩幀圖像，在不考慮丟幀的情況下，根據相機幀頻確定連續(xù)兩幀圖像之間的時間間隔為0.002 96 s。原始采集圖像如圖7所示，通過圖像灰度化增強、降噪處理等圖像預處理方法提高圖像質量以保證后續(xù)流場參數的計算精度。圖像標定精度決定流場參數提取精度。通過圖像外標定方法獲取圖像放大率，圖像預處理后采用Canny邊緣檢測算法識別圖像邊緣特征[26]，提取渦輪外環(huán)多段半徑數值，如圖8(a)所示，將檢測的半徑均值與實際外環(huán)半徑數值進行對比，獲得的圖像放大率為0.191 mm/px。標定結果見表1。

圖6 PIV試驗系統(tǒng)

(a)第1幀 (b)第2幀

(a)半徑檢測 (b)實際尺寸

基于歸一化遞歸互相關算法提取渦輪流場參數[27]。為了節(jié)省計算時間，只取渦輪單獨流道為計算區(qū)域，其他流場區(qū)域進行掩膜處理，提取流速場和渦量場，如圖9所示。為了定量驗證CFD仿真結果的準確性，以渦輪中心為圓心，對渦輪單獨流道進行圓弧等分劃分，在流道中間區(qū)域水平位置提取8個等分點，如圖9所示，提取等分點處流速值和渦量值，見表2和表3。

表1 標定結果

(a)流速等分點

(b)渦量等分點

表2 渦輪單流道流速場中間截面流速值

表3 渦輪單流道渦量場中間截面渦量值

5 結果與討論

5.1 流速場及渦量場

為了便于仿真與試驗結果對比分析，截取渦輪單流道流場仿真結果，如圖10所示。從流速場整體結構來看，不同湍流模型仿真結果中渦輪流場均呈現(xiàn)明顯的旋渦流動，但是旋渦結構特征及渦心位置均不同，差異性明顯。SL模型渦心分布在靠近葉片吸力面的中間位置處，呈現(xiàn)為“鞋墊狀”的旋渦結構；WALE模型渦心向右下方偏移，相較于SL模型，其呈現(xiàn)為“反鞋墊狀”的非包容型旋渦結構，且壓力面處的高流速區(qū)域相對于其他模型仿真效果更加明顯。WMLES模型渦心居中且靠近吸力面，呈現(xiàn)為縱向垂直狀的“心臟形”旋渦結構，其下方的高流速區(qū)域相較于其他模型更加明顯，且流速值較高。WMLES S-Omega模型渦心位置與WMLES模型接近，但是封閉狀的旋渦結構不明顯，低速區(qū)擴展到吸力面靠近出口位置。KET模型形成了兩個渦心，與其他模型仿真結果均不相同，兩個渦心在空間布局上形似“葫蘆狀”，較大的旋渦靠近主流區(qū)域中部，較小的旋渦位于其左上方靠近吸力面，它們的旋轉方向相同，都是逆時針方向。

由于來自泵輪的高速流體沖擊渦輪，渦輪入口處葉片吸力面與外環(huán)的交匯區(qū)域內出現(xiàn)大范圍高速區(qū)，WALE和KET模型仿真結果中所呈現(xiàn)的高流速區(qū)域范圍分布大，流速在3.8～4.6 m/s，如圖10(a)所示，這與PIV試驗結果中的B區(qū)域相對應，流速為1.7～3.0 m/s。SL和WMLES模型在低速區(qū)域仿真結果類似，流速在0.32～0.85 m/s，如圖10(a)所示。對比PIV試驗結果，如圖9(a)中D區(qū)域存在明顯低速區(qū)，流速為0.44～0.55 m/s。WMLES和WMLES S-Omega模型能夠捕捉到葉片近壁面區(qū)域小尺度渦旋結構，WALES模型仿真結果中在渦輪流道壓力面附近存在明顯的高流速區(qū)，流速在3.7～4.4 m/s，與PIV的試驗結果接近，如圖9(a)中C區(qū)域所示，流速為1.7～2.8 m/s，而WMLES S-Omega模型在該區(qū)域上流速過于平滑且數值偏低，仿真結果失真。

從渦量場結構分布來看，采用5種亞格子湍流模型對葉片近壁面流場仿真結果類似，但是主流區(qū)域上渦量場結構分布差異明顯，WALE模型在主流區(qū)域上的仿真結果存在兩個高渦量場結構，呈左右分布，但左側高渦量場區(qū)域分布范圍較大。WMLES模型高渦量場結構主要集中在主流區(qū)域偏下方位置處，呈橫向分布。KET模型仿真結果存在兩個高渦量場結構，分布范圍較大的高渦量場結構靠近吸力面一側，分布范圍較小的高渦量場結構位于主流核心區(qū)域，SL和WMLES S-Omega模型在主流區(qū)域上不存在明顯的高渦量場結構，如圖10(b)所示。對比PIV試驗結果，如圖9(b)中F、G、H區(qū)域上存在高渦量場結構，且渦量分別為582.17～793.42 s-1、425.41～596.79 s-1、471.48～612.37 s-1。針對F區(qū)域，數值模擬結果與試驗結果不對應，這是因為PIV試驗測量中液力變矩器外環(huán)與渦輪葉片交匯處出現(xiàn)明顯的反光現(xiàn)象，如圖7所示，因此導致PIV渦量場計算結果失真。針對G區(qū)域，WMLES模型仿真結果與PIV試驗結果相近，渦量為475.47～840.36 s-1。針對H區(qū)域，KET模型仿真結果與PIV試驗結果一致，渦量為426.37～782.46 s-1。

(a)流速 (b)渦量

為了量化對比數值模擬與試驗測量結果，參照試驗測量結果圖9，將渦輪單獨流道流場仿真結果進行等分圓弧劃分，在流道中間區(qū)域水平位置提取8個等分點，保證數值模擬結果與試驗測量結果比對點保持一致。由于文章篇幅有限，數據處理方法相同，以SL湍流模型仿真結果為例，提取流場8個等分點上的流速值和渦量值，如圖11所示，仿真數據提取結果見表2和表3。

(a)流速場等分點 (b)渦量場等分點

圖12為數值模擬結果與PIV試驗結果對比圖。從流速數值分布趨勢上來看，如圖12(a)所示，5種亞格子湍流模型仿真結果與PIV試驗結果分布走勢一致，但仿真數值均大于試驗數值，這與PIV流場圖像標定精度有關。綜合數據曲線整體趨勢和測量數據點數值兩個方面，WMLES S-Omega模型仿真結果與PIV試驗結果更加接近。從渦量數值分布趨勢來看，無明顯規(guī)律，尤其是數據點3～6上渦量數值發(fā)生巨大躍變，如圖12(b)所示，這是由于制動工況下渦輪流道主流區(qū)域上出現(xiàn)大尺度旋渦流動，流場梯度劇烈變化，在整體旋渦運動進程中同時伴隨著多尺度渦系結構的碰撞、交融和混合等動力學過程，導致渦量數值雜亂無章。從離散點的渦量數值比對來看，WMLES S-Omega模型仿真結果與試驗結果誤差最小。

(a)流速

(b)渦量

6 結論

(1)從流速場結構分布來看，采用SL和WMLES模型流場仿真結果類似，采用WMLES模型對葉片壓力面高流速區(qū)域仿真效果較好，與PIV試驗一致。從流速數值定量角度來看，WMLES S-Omega模型與PIV試驗結果最接近。

(2)從渦量場結構分布來看，5種亞格子湍流模型對于葉片近壁面渦量場仿真相近，但是主流區(qū)域上渦量場結構分布差異性較大。WMLES模型仿真結果所呈現(xiàn)的高渦量區(qū)域與PIV試驗結果最接近。從渦量數值定量角度來看，WMLES S-Omega模型仿真結果與試驗結果誤差最小。

(3)從三維渦時空結構分布角度來看，WMLES和WMLES S-Omega模型能夠時空再現(xiàn)葉片近壁面豐富的三維渦結構。KET模型能夠捕捉到渦輪葉片尾跡區(qū)域上明顯的脫落渦結構，對于泵輪與渦輪無葉柵區(qū)域上的三維渦結構仿真效果較好。