陷窩誘導(dǎo)渦結(jié)構(gòu)數(shù)值模擬分析

劉 靜，李 杰，蔣勝矩

（1.西北工業(yè)大學(xué) 翼型葉柵空氣動(dòng)力學(xué)國(guó)防科技重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710072；2.中國(guó)兵器工業(yè)第203研究所，陜西 西安 710065）

0 引 言

旋渦是流動(dòng)的肌腱，對(duì)于流場(chǎng)中旋渦的準(zhǔn)確認(rèn)識(shí)是發(fā)展科學(xué)的旋渦控制和應(yīng)用技術(shù)的前提。陷窩作為一種渦流發(fā)生器，通過(guò)誘導(dǎo)產(chǎn)生渦結(jié)構(gòu)注入邊界層層動(dòng)量，在流動(dòng)分離控制、減阻及增強(qiáng)對(duì)流換熱等很多方面都有應(yīng)用研究［1－9］。陷窩具有在不帶來(lái)明顯的壓力損失的前提下，引起流場(chǎng)持續(xù)的擾動(dòng)的優(yōu)良特點(diǎn)，這得益于陷窩所誘導(dǎo)渦結(jié)構(gòu)奇特的流場(chǎng)擾動(dòng)方式。

國(guó)內(nèi)外用數(shù)值模擬和實(shí)驗(yàn)的手段進(jìn)行陷窩的有關(guān)研究。目前的數(shù)值研究方向分為兩類，第一類是以研究陷窩誘導(dǎo)渦結(jié)構(gòu)非定常脫落的基于URANS和LES的研究；第二類是基于RANS的穩(wěn)態(tài)假設(shè)下陷窩誘導(dǎo)渦結(jié)構(gòu)及其應(yīng)用的研究。URANS和LES方法著重研究陷窩誘導(dǎo)渦的非定常脫落性質(zhì)。RANS則側(cè)重于得到穩(wěn)態(tài)假設(shè)下的渦結(jié)構(gòu)、傳熱性能和流動(dòng)控制等的時(shí)均特征。Isaev［1－3］等人開(kāi)展基于 RANS的陷窩強(qiáng)化傳熱定常數(shù)值模擬研究，通過(guò)表面極限流線發(fā)現(xiàn)，深寬比0.2陷窩內(nèi)的渦結(jié)構(gòu)為雙單元類龍卷風(fēng)渦結(jié)構(gòu)，類龍卷風(fēng)渦從物面螺旋點(diǎn)升起，能引起流場(chǎng)的持續(xù)擾動(dòng)，強(qiáng)化換熱。Isaev［4］用 URANS和LES方法進(jìn)行了陷窩誘導(dǎo)渦結(jié)構(gòu)的非定常脫落頻率的研究。謝永慧［5］用LES研究陷窩控制邊界層分離研究中發(fā)現(xiàn)深寬比0.09陷窩內(nèi)誘導(dǎo)的渦結(jié)構(gòu)為馬蹄渦，馬蹄渦腿源自于物面焦點(diǎn)，馬蹄渦頭位于下游，馬蹄渦周期性脫落，并在陷窩后形成發(fā)夾渦排。劉高文［7－8］用實(shí)驗(yàn)測(cè)量了深寬比0.2陷窩內(nèi)的渦結(jié)構(gòu)，發(fā)現(xiàn)會(huì)在陷窩背風(fēng)面附近形成一個(gè)低速回流區(qū)，并發(fā)現(xiàn)氣流經(jīng)過(guò)陷窩后誘導(dǎo)形成一對(duì)尺寸較大的旋渦附著于尾跡區(qū)。何雅玲［9］用RANS研究陷窩強(qiáng)化換熱，發(fā)現(xiàn)陷窩的介入使得流動(dòng)中產(chǎn)生了旋渦死區(qū)，死區(qū)中的流體不容易從其中流出，降低了換熱效果。

現(xiàn)有研究主要通過(guò)物面極限流線和空間流線來(lái)研究陷窩內(nèi)渦結(jié)構(gòu)，也沒(méi)有對(duì)渦核處壓力和軸向速度的細(xì)節(jié)分析的文獻(xiàn)，為了澄清陷窩誘導(dǎo)渦結(jié)構(gòu)及其對(duì)尾流的擾動(dòng)方式，本文對(duì)布置單個(gè)陷窩兩平板間充分發(fā)展段流動(dòng)進(jìn)行RANS數(shù)值模擬。首次利用張涵信的旋渦沿軸線的非線性分叉理論［10－11］和沿渦軸的軸向速度和壓力分布等細(xì)致分析在穩(wěn)態(tài)假設(shè)下陷窩內(nèi)渦結(jié)構(gòu)發(fā)展、演變及破裂的規(guī)律。研究顯示本文模擬的陷窩內(nèi)渦結(jié)構(gòu)為失穩(wěn)破裂的半渦環(huán)，破裂后渦量擴(kuò)散形成弱縱向渦，形成對(duì)尾流的擾動(dòng)。陷窩內(nèi)的渦結(jié)構(gòu)決定了陷窩后流場(chǎng)的擾動(dòng)方式和強(qiáng)度。本文研究對(duì)陷窩的后續(xù)研究有指導(dǎo)意義。

1 物理模型的建立

物理模型為底板布置深寬比0.2陷窩的相距0.2m兩平板之間充分發(fā)展段流動(dòng)。計(jì)算區(qū)域取底板到兩平板等分面之間，等分面設(shè)為對(duì)稱面。計(jì)算模型如圖1所示，矩形計(jì)算區(qū)域長(zhǎng)0.406 4m，寬0.254m，高0.1m。流向?yàn)閤向，展向?yàn)閥向，高為z向。單個(gè)陷窩布置在底板上，陷窩中心和底板中心重合。陷窩表面直徑為0.025 4m，陷窩深度為0.005 08m。陷窩深寬比為0.2，陷窩深度和兩平板間距離之比為0.025。陷窩深度的設(shè)計(jì)依據(jù)過(guò)增元的場(chǎng)協(xié)同理論，即只有在粘性底層附近產(chǎn)生縱向渦，才能實(shí)現(xiàn)同功耗下湍流對(duì)流強(qiáng)化。本文計(jì)算入口速度型粘性底層和過(guò)渡層的總厚度約為0.004 2m，即陷窩深度約為粘性底層和過(guò)渡層的總厚度的1.2倍，見(jiàn)圖2。

圖1 計(jì)算模型示意圖Fig.1 Computational model

圖2 陷窩前速度型和陷窩截面圖Fig.2 Velocity profile before dimple

模型側(cè)面為對(duì)稱面，布置陷窩底板及陷窩為無(wú)滑移物面。為消除入口段影響，計(jì)算入口邊界設(shè)為充分發(fā)展邊界條件。入口邊界條件為速度入口。初始速度均勻分布為4m／s。通過(guò)流量守恒迭代得到流動(dòng)充分發(fā)展入口速度（圖3），湍流強(qiáng)度分布如圖4。出口邊界條件是壓力出口，出口相對(duì)靜壓為0Pa。基于入口平均速度4m／s和兩平板間距離0.2m的Re數(shù)為5.48×104。假設(shè)空氣不可壓縮。

圖3 x方向速度分布Fig.3 x－velocity distribution

圖4 湍動(dòng)能分布Fig.4 TKE distribution

2 數(shù)值計(jì)算方法

2.1 數(shù)值方法

本文用Fluent進(jìn)行兩平板間充分發(fā)展流動(dòng)定常RANS數(shù)值模擬。計(jì)算選擇基于壓力的求解器，速度壓力耦合采用Coupled方式。湍流模擬選擇k－ωSST湍流模型。k－ωSST湍流模型在壁面粘性模擬中采用增強(qiáng)壁面法，直接求解粘性底層的流動(dòng)。動(dòng)量方程，湍流強(qiáng)度k，比耗散率ω和能量方程的對(duì)流項(xiàng)離散采用三階MUSCL格式。

2.2 計(jì)算網(wǎng)格

空間離散采用ICEM生成多塊結(jié)構(gòu)網(wǎng)格。盡量保持網(wǎng)格方向和流動(dòng)方向的一致性，以及較好的正交性和較小的長(zhǎng)寬比。如圖5為整體和陷窩局部網(wǎng)格，為展示布置陷窩底板網(wǎng)格結(jié)構(gòu)，隱去頂面網(wǎng)格。因三個(gè)面均為對(duì)稱面，只需在底板法向加密網(wǎng)格，簡(jiǎn)化了計(jì)算模型。從圖中可見(jiàn)對(duì)無(wú)滑移底板邊界層網(wǎng)格進(jìn)行加密。無(wú)滑移物面的法向第一個(gè)網(wǎng)格滿足y＋＜1，并在粘性底層（y＋＜10）的區(qū)域內(nèi)有多于10個(gè)網(wǎng)格。采用O型網(wǎng)格拓?fù)浣Y(jié)構(gòu)對(duì)陷窩前后緣和陷窩內(nèi)部進(jìn)行局部加密。總網(wǎng)格單元數(shù)為169萬(wàn)。

圖5 模型整體和陷窩局部網(wǎng)格Fig.5 Model grid

在渦核處網(wǎng)格密度和空間離散格式精度會(huì)引起超出渦物理粘性耗散的強(qiáng)渦數(shù)值耗散，致使對(duì)旋渦結(jié)構(gòu)的模擬不準(zhǔn)確。為保證計(jì)算結(jié)果的準(zhǔn)確性，在原始169萬(wàn)網(wǎng)格的基礎(chǔ)上在研究關(guān)注的渦核附近區(qū)域加密得到269萬(wàn)網(wǎng)格。兩套網(wǎng)格渦核軸向速度對(duì)比如圖6，兩套不同網(wǎng)格模擬得到的渦結(jié)構(gòu)基本重合，即本文網(wǎng)格的節(jié)點(diǎn)布置和拓?fù)浣Y(jié)構(gòu)劃分能夠較好地模擬陷窩誘導(dǎo)的渦結(jié)構(gòu)。

圖6 不同網(wǎng)格計(jì)算渦核軸向速度對(duì)比Fig.6 Grid comparison of vortex core axial velocity

2.3 計(jì)算和實(shí)驗(yàn)對(duì)比驗(yàn)證

本文采用和文獻(xiàn)［8］的設(shè)置對(duì)陷窩強(qiáng)化對(duì)流傳熱進(jìn)行數(shù)值模擬，并與文獻(xiàn)［8］的實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，如圖7所示。光滑管道計(jì)算值為Nu0。計(jì)算與實(shí)驗(yàn)結(jié)果的Nu／Nu0相對(duì)誤差在6%內(nèi)，說(shuō)明了本文數(shù)值模擬方法的可靠性。

圖7 數(shù)值模擬和實(shí)驗(yàn)對(duì)比Fig.7 Comparison of numerical simulation and experiment

3 計(jì)算結(jié)果與分析

3.1 半渦環(huán)渦結(jié)構(gòu)分析

圖8所示為采用渦矢量方法所得的陷窩內(nèi)渦核俯視圖和側(cè)視圖。從圖中可以看出：渦核呈對(duì)稱分布；在兩側(cè)壁面起始段與物面垂直；上升到一定高度后，渦核沿流向發(fā)展至對(duì)稱面。結(jié)合圖9物面的極限流線可以看到物面極限流線流譜中奇點(diǎn)的分布為兩個(gè)鞍點(diǎn)＋兩個(gè)分離螺旋點(diǎn)。分離線上奇點(diǎn)分布為螺旋點(diǎn)／鞍點(diǎn)／螺旋點(diǎn)。分離屬于典型的閉式分離。物面兩個(gè)對(duì)稱穩(wěn)定的分離螺旋點(diǎn)對(duì)應(yīng)于空間對(duì)稱的類龍卷風(fēng)渦。類龍卷風(fēng)渦的形成條件為物面逆流的出現(xiàn)，分離螺旋點(diǎn)的形成，以及類龍卷風(fēng)的垂直渦升起。陷窩前緣的分離和迎風(fēng)面的再附形成了陷窩背風(fēng)面的逆流，從而形成了兩個(gè)對(duì)稱的分離螺旋點(diǎn)。物面分離螺旋點(diǎn)形成以后，匯聚了物面的部分渦量升起，在空間演化為垂直物面發(fā)展的類龍卷風(fēng)渦。閉式分離流動(dòng)具有明顯的禁區(qū)性特點(diǎn)，給分離渦輸送渦量的邊界層位于閉式分離線的下游。對(duì)稱類龍卷風(fēng)渦的渦量由分離渦面在分離線下游卷繞渦面提供。隨著渦沿流向的發(fā)展，渦核從垂直物面發(fā)展轉(zhuǎn)為沿流向向?qū)ΨQ面發(fā)展。從渦核的空間結(jié)構(gòu)圖可見(jiàn)，由對(duì)稱的類龍卷風(fēng)渦形成的雙渦結(jié)構(gòu)在對(duì)稱面閉合，形成半渦環(huán)渦結(jié)構(gòu)［12］。半渦環(huán)渦結(jié)構(gòu)指兩個(gè)渦軸端點(diǎn)終止于壁面分離螺旋點(diǎn)的半環(huán)型渦結(jié)構(gòu)，陷窩內(nèi)半渦環(huán)示意見(jiàn)圖10。

圖8 渦核線空間分布俯視圖和側(cè)視圖Fig.8 Vortex core vertical and side view

圖9 物面的極限流線Fig.9 Wall streamline

圖10 半渦環(huán)渦結(jié)構(gòu)示意圖Fig.10 Half ring vortex schematic drawing

從圖11可見(jiàn)沿渦核軸向流動(dòng)一直在逆壓梯度下。在渦的起始段，逆壓梯度較?。浑S著渦的發(fā)展，逆壓梯度逐漸增加。在x＝0.204至x＝0.205的渦軸位置，逆壓梯度增幅較大，對(duì)應(yīng)于圖12軸向速度減小至零，并出現(xiàn)軸向逆流的位置。旋渦對(duì)逆壓梯度很敏感，逆壓梯度的突然增加會(huì)引起旋渦的不穩(wěn)定性甚至破裂，同時(shí)會(huì)帶來(lái)軸向速度的減小和軸向逆流的出現(xiàn)。從圖12可見(jiàn)渦軸的軸向速度最大為0.1m／s，本文計(jì)算的入口平均速度為4m／s，相對(duì)于入口速度，渦軸的軸向速度值很小。閉式分離明顯的禁區(qū)性特點(diǎn)使得來(lái)流只能在分離線的下側(cè)回流進(jìn)入旋渦，動(dòng)量輸入小。

下面結(jié)合圖13垂直渦軸橫截面流態(tài)討論半渦環(huán)渦結(jié)構(gòu)沿流向的演變過(guò)程。圖13中陷窩外主流流場(chǎng)流向從左向右。在x＜0.202范圍內(nèi)，渦軸軸向速度加速度λ＞0，軸向速度增加，在渦軸附近，截面流態(tài)為穩(wěn)定的螺旋點(diǎn)形態(tài)，流線由外部指向旋渦中心，見(jiàn)圖13（a）。類龍卷風(fēng)渦結(jié)構(gòu)穩(wěn)定。渦核處因粘性耗散及粘性擴(kuò)散損失的渦量得到流向渦軸的渦量補(bǔ)給，渦量平衡，旋渦穩(wěn)定。在0.202＜x＜0.206范圍內(nèi)，λ＜0，軸向速度減小。沿渦軸線λ變號(hào)，由正變?yōu)樨?fù)，則從變號(hào)點(diǎn)起，其渦心附近的截面流線圖上出現(xiàn)穩(wěn)定的極限環(huán)，即出現(xiàn)Hopf分叉，見(jiàn)圖13（b）。穩(wěn)定極限環(huán)外部流線由外部指向極限環(huán)，流動(dòng)穩(wěn)定；穩(wěn)定極限環(huán)內(nèi)部為不穩(wěn)定的螺旋點(diǎn)形態(tài)，流線由旋渦中心指向旋渦外部。即旋渦的不穩(wěn)定最先出現(xiàn)在渦軸附近。隨著流動(dòng)沿著渦軸的發(fā)展，極限環(huán)逐漸變大，不穩(wěn)定區(qū)域也逐漸變大，見(jiàn)圖13（c）。在x＝0.204位置，極限環(huán)消失，不穩(wěn)定區(qū)域占據(jù)了全部的分離區(qū)域，見(jiàn)圖13（d）。在x＝0.205 4位置，截面流線拓?fù)浣Y(jié)構(gòu)保持不變，見(jiàn)圖13（e）。直至x＞0.206時(shí)，軸向速度加速度λ＞0，沿渦軸線λ變號(hào)，由負(fù)變?yōu)檎?，則從變號(hào)點(diǎn)起，其渦心附近的截面流線圖上開(kāi)始出現(xiàn)不穩(wěn)定的極限環(huán)，極限環(huán)內(nèi)部出現(xiàn)穩(wěn)定螺旋點(diǎn)。在0.206 6位置時(shí)，穩(wěn)定螺旋點(diǎn)已經(jīng)占據(jù)了較大的空間，見(jiàn)圖13（f）。

圖11 軸向壓力分布Fig.11 Vortex core pressure distribution

圖12 軸向速度分布Fig.12 Vortex core axial velocity distribution

圖13 垂直渦軸橫截面流線圖Fig.13 Slice streamline vertical to vortex core

數(shù)值模擬結(jié)果中的截面拓?fù)浣Y(jié)構(gòu)變化規(guī)律以及極限環(huán)的產(chǎn)生與張涵信的拓?fù)浞治龆ㄐ苑治隼碚撌窍喾?，說(shuō)明了本文數(shù)值模擬結(jié)果和分析結(jié)果的正確性。

3.2 半渦環(huán)渦破裂

如圖14空間流線圖可見(jiàn)在渦軸與x＝0.205相交處布置種子，計(jì)算所得空間流線圖?？梢?jiàn)在x＝0.205位置之前緊緊繞渦核旋轉(zhuǎn)流動(dòng)的流線突然擴(kuò)張，旋渦破裂。且由圖8渦核0.205處軸向速度降為零，同時(shí)壓力梯度增幅陡增。由圖9渦核空間俯視圖可見(jiàn)，半渦環(huán)渦核在對(duì)稱面出現(xiàn)拐折。渦環(huán)受到擾動(dòng)后，失去穩(wěn)定性，圓環(huán)變成多邊形形狀［13］。作為半渦環(huán)，渦核出現(xiàn)較明顯的夾角狀態(tài)意味著半渦環(huán)失穩(wěn)。從圖13垂直渦軸截面流線演變過(guò)程可見(jiàn)，在0.205位置不穩(wěn)定螺旋點(diǎn)已占據(jù)了絕大部分區(qū)域，旋渦不穩(wěn)定。由渦破裂理論知，此處渦破裂，其破裂的起點(diǎn)為駐點(diǎn)。渦破裂的主要特征是旋渦的軸向速度迅速下降，形成內(nèi)駐點(diǎn)，駐點(diǎn)后為回流區(qū)，渦核半徑突然擴(kuò)張，形成一個(gè)近似軸對(duì)稱的破裂泡［10－11］。

從λ＞0到λ＜0，由λ＝0的位置起，橫截面流線出現(xiàn)Hopf分叉產(chǎn)生穩(wěn)定的極限環(huán)。該極限環(huán)的邊界，將流動(dòng)分成內(nèi)外兩區(qū)，外區(qū)截面流線由外向內(nèi)轉(zhuǎn)而指向極限環(huán)；內(nèi)區(qū)是不穩(wěn)定的，截面流線由渦心向外轉(zhuǎn)而指向極限環(huán)。這就是說(shuō)，旋渦的失穩(wěn)，首先發(fā)生在渦軸附近。這樣，渦軸附近處于不穩(wěn)區(qū)的流體，在向前運(yùn)動(dòng)且不斷減速的過(guò)程中，就可能因失穩(wěn)而破裂。根據(jù)旋渦沿軸線的非線性分叉理論［11］，當(dāng)破裂點(diǎn)出現(xiàn)在軸向速度為零的位置時(shí)，其縱向截面流線為鞍點(diǎn)形態(tài)，如圖15所示，則說(shuō)明在x＝0.205破裂類型為泡型破裂。

圖14 過(guò)渦核駐點(diǎn)的空間流線Fig.14 Streamline through stagnation point

圖15 過(guò)渦核駐點(diǎn)縱向截面鞍點(diǎn)流態(tài)Fig.15 Slice streamline align to vortex core through stagnation point

一般旋渦的運(yùn)動(dòng)分區(qū)為粘性起主導(dǎo)作用的渦核區(qū)和渦核外近似無(wú)粘準(zhǔn)自由渦區(qū)［14］。準(zhǔn)自由渦區(qū)為渦核半徑以外的區(qū)域，此處流體的運(yùn)動(dòng)近似為無(wú)粘流，主要由壓力場(chǎng)決定，基本不受粘性力的作用，所以總壓變化幅度不大。即在準(zhǔn)自由渦區(qū)，存在壓力梯度，但總壓梯度基本為零。粘性主導(dǎo)渦核區(qū)為從最大渦核半徑處到軸心，切向速度從峰值迅速減小到零，就像邊界層的速度型，從邊界層頂部的速度很快降低到零。渦核區(qū)域是粘性高度集中的區(qū)域，由粘性力和壓力共同作用，由于粘性耗散，使渦核區(qū)域處總壓損失，有較大的總壓梯度。因此，從流場(chǎng)的低總壓分布區(qū)域和總壓梯度較大的區(qū)域可以準(zhǔn)確地劃定有強(qiáng)粘性耗散的渦核區(qū)域。從半渦環(huán)總壓切片圖16可以看到渦核處總壓從開(kāi)始有較大的梯度到梯度逐漸變小，直到0.204截面總壓最小值超出了所繪制云圖的上限，低總壓區(qū)逐漸消失，即強(qiáng)粘性耗散區(qū)逐漸消失。渦核尺寸在粘性擴(kuò)散的作用下逐漸增大，渦核強(qiáng)度在粘性耗散和粘性擴(kuò)散的共同作用下逐漸變?nèi)?，以致破裂?/p>

圖16 渦核處總壓等x切片圖Fig.16 Equal xslice total pressure distribution

以上種種證實(shí)了半渦環(huán)在x＝0.205位置渦破裂的事實(shí)。

從已發(fā)表文獻(xiàn)［1－9］中展示的流場(chǎng)結(jié)構(gòu)中可見(jiàn)，Isaev［1－3］和謝永惠［5］的數(shù)值模擬結(jié)果中都顯示出螺旋點(diǎn)／鞍點(diǎn)／螺旋點(diǎn)的閉式分離，兩個(gè)對(duì)稱的物面螺旋型焦點(diǎn)對(duì)應(yīng)于空間旋渦的物面起始點(diǎn)，和本文的數(shù)值模擬結(jié)果一致。

本文認(rèn)為同樣是物面兩個(gè)對(duì)稱的螺旋型焦點(diǎn)在空間里可能形成兩個(gè)完全不同的渦結(jié)構(gòu)。如果渦核在對(duì)稱面閉合，渦結(jié)構(gòu)應(yīng)為兩個(gè)渦軸端點(diǎn)終止于壁面的半渦環(huán)渦結(jié)構(gòu)；如果渦核在對(duì)稱面不閉合，沿流向延伸，則渦結(jié)構(gòu)應(yīng)為雙單元類龍卷風(fēng)渦。本文數(shù)值模擬的陷窩內(nèi)渦結(jié)構(gòu)定義為中間失穩(wěn)破裂的半渦環(huán)較合理。

3.3 旋渦對(duì)陷窩后流場(chǎng)影響

如圖17為陷窩后等x截面流線圖。從圖中可以看到陷窩后流體有縱向旋轉(zhuǎn)的趨勢(shì)。這是陷窩內(nèi)的半渦環(huán)破裂后渦量在下游流場(chǎng)中分散分布的結(jié)果。旋渦破裂后渦核內(nèi)的集中渦量分散開(kāi)來(lái)，隨著流體的向后流動(dòng)向下游輸運(yùn)，形成陷窩后的渦量分布，并誘導(dǎo)陷窩后流體橫向流動(dòng)。圖中貼近物面的截面流線旋渦是陷窩前緣分離流體在陷窩迎風(fēng)面再附后部分流體向后繞尖后緣流動(dòng)形成的對(duì)稱的邊渦對(duì)。主渦對(duì)破裂渦量在空間散開(kāi)，分布空間較大，渦量強(qiáng)度較小，誘導(dǎo)形成陷窩后弱縱向渦；邊渦對(duì)緊貼物面形成，分布集中，渦量強(qiáng)度較大。它們旋轉(zhuǎn)方向相同，共同形成對(duì)陷窩后流場(chǎng)的擾動(dòng)，加強(qiáng)了陷窩后流場(chǎng)的對(duì)流強(qiáng)度。劉高文［8］也試驗(yàn)研究發(fā)現(xiàn)氣流經(jīng)過(guò)陷窩后誘導(dǎo)形成一對(duì)尺寸較大的旋渦附著于尾跡區(qū)。

圖17 陷窩后縱向渦截面流線Fig.17 Wake longitudinal vortex slice streamline

4 結(jié) 論

本文對(duì)布置深寬比為0.2陷窩兩平板間充分發(fā)展段流動(dòng)進(jìn)行了數(shù)值模擬，利用張涵信的旋渦沿軸線的非線性分叉理論和沿渦軸的軸向速度和壓力分布等細(xì)致分析在穩(wěn)態(tài)假設(shè)下陷窩內(nèi)渦結(jié)構(gòu)發(fā)展、演變及破裂的規(guī)律。主要結(jié)論為：

（1）對(duì)于物面流線有兩個(gè)螺旋型焦點(diǎn)的陷窩內(nèi)流動(dòng)，當(dāng)渦核在對(duì)稱面閉合，渦結(jié)構(gòu)應(yīng)為半渦環(huán)；當(dāng)渦核在對(duì)稱面不閉合，渦結(jié)構(gòu)應(yīng)為雙單元類龍卷風(fēng)渦。

（2）本文數(shù)值模擬陷窩誘導(dǎo)的渦結(jié)構(gòu)為中間失穩(wěn)破裂的半渦環(huán)結(jié)構(gòu)。半渦環(huán)結(jié)構(gòu)由對(duì)稱的類龍卷風(fēng)渦在對(duì)稱面閉合形成。物面分離流譜為典型的螺旋點(diǎn)／鞍點(diǎn)／螺旋點(diǎn)分離。螺旋分離點(diǎn)為類龍卷風(fēng)渦在物面的起始點(diǎn)。物面分離螺旋點(diǎn)形成后，匯聚了物面的部分渦量升起，在空間演化為垂直物面發(fā)展的類龍卷風(fēng)渦。

（3）半渦環(huán)的發(fā)展經(jīng)歷了從垂直于物面穩(wěn)定升起，截面流線在渦心附近從穩(wěn)定的螺旋點(diǎn)形態(tài)；渦心附近出現(xiàn)不穩(wěn)定的螺旋點(diǎn)形態(tài)；再到渦泡型破裂的過(guò)程。半渦環(huán)經(jīng)過(guò)了穩(wěn)定升起，沿流向隨著渦量粘性擴(kuò)散和渦核的粘性耗散下變得不穩(wěn)定，最后在強(qiáng)逆壓梯度下泡型渦破裂的過(guò)程。

（4）半渦環(huán)渦破裂渦量散開(kāi)誘導(dǎo)形成陷窩后的弱縱向渦，和陷窩尖后緣繞流產(chǎn)生的邊渦旋轉(zhuǎn)同向，加強(qiáng)了陷窩后流場(chǎng)的對(duì)流強(qiáng)度。

［1］Isaev S A.The effect of rearrangement of the vortex structure on heat transfer under conditions of increasing depth of a spherical dimple on the wall of a narrow channel［J］.High Temperature，2003，41（2）：229－232.

［2］Isaev S A.Numerical simulation of vortex enhancement of heat transfer［J］.High Temperature，2003，41（5）：665－679.

［3］Isaev S A.Influence of the Reynolds number and the spherical dimple depth on turbulent heat transfer and hydraulic loss in a narrow channel［J］.International Journal of Heat and Mass Transfer，2010，53：178－197.

［4］Turnow J，Kornev N，Isaev S A.Vortex mechanism of heat transfer enhancement in a channel with spherical and oval dimples［J］.Heat and Mass Transfer，2011，47（3）：301－313.

［5］LAN Jibing，Xie Yonghui，Zhang Di.Large eddy simulation for control of boundary layer separation with dimples，Journal of Xi’an Jiaotong University，2010，44（9）：27－32.（in Chinese）藍(lán)吉兵，謝永惠，張荻.采用球窩控制邊界層分離流動(dòng)的大渦模擬［J］，西安交通大學(xué)學(xué)報(bào).2010，44（9）：27－32.

［6］Qiao Weiyang.Control of separation losses on turbine blade with low Reynolds number［J］，Journal of propulsion technology，2005，26（1）：42－45.（in Chinese）喬渭陽(yáng).低雷諾數(shù)渦輪流動(dòng)損失控制技術(shù) ［J］，推進(jìn)技術(shù).2005，26（1）：42－45.

［7］Liu Gaowen.Review of heat transfer enhancement for dimpled surface［J］.Journal of Aerospace Power，2007，22（11）：1785－1791.（in Chinese）劉高文.凹坑強(qiáng)化傳熱的研究進(jìn)展回顧［J］.航空動(dòng)力學(xué)報(bào)，2007，22（11）：1785－1791.

［8］Liu Gaowen.Transient infrared heat transfer measurements and surface flow visualization on a single dimpled wall［J］.Journal of Propulsion Technology，2009，30（1）：41－45.（in Chinese）劉高文.單凹坑壁面的瞬態(tài)紅外傳熱測(cè)量與流場(chǎng)顯示［J］.推進(jìn)技術(shù)，2009，30（1）：41－45.

［9］Li Rui，He Yaling.Numerical simulation of dimpled tubes for heat transfer enhancement［J］.Journal of Engineering Thermophysics，2008，29（11）：1947－1949.（in Chinese）李瑞，何雅玲.丁胞型強(qiáng)化換熱管的數(shù)值模擬［J］.工程熱物理學(xué)報(bào)，2008，29（11）：1947－1949.

［10］Zhang Hanxin，Deng Xiaogang.Analytic studies for three dimensional steady separated flows and vortex motion［J］.Acta Aerodynamica Sinica，1992，10（1）：8－20.（in Chinese）張涵信，鄧小剛.三維定常分離流和渦運(yùn)動(dòng)的定性分析研究［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，1992，10（1）：8－20.

［11］Zhang Hanxin.Bifurcation of vortex motion along its axis［J］.Acta Aerodynamica Sinica，1994，12（3）：242－251.（in Chinese）張涵信.旋渦沿軸線的非線性分叉［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，1994，12（3）：242－251.

［12］Wang Xiaoxin，Wu Ziniu.Study on wake vortex influence in insect flapping hovering flight［D］.［PhD Thesis］.Beijing：Tsinghua University，2012.（in Chinese）王曉欣，吳子牛.昆蟲(chóng)撲翼懸停飛行尾渦影響研究［D］.［博士學(xué)位論文］.北京：清華大學(xué)，2012.

［13］Tong Bingang.Unsteady flow and vortex motion［M］.Beijing：National Defence Industry Press，1993.（in Chinese）童秉綱.非定常流與渦運(yùn)動(dòng)［M］.北京：國(guó)防工業(yè)出版社，1993.

［14］Xia Xuejian，Deng Xueyin.Engineering separation dynamic flow［M］.Beijing：Beihang University Press，1991.夏學(xué)湔，鄧學(xué)鎣.工程分離流動(dòng)力學(xué)［M］.北京：北京航空航天大學(xué)出版社，1991.

［15］Wang Gang，Liang Xingang.Tornado－like vortices generation and its aerodynamic characteristics［J］.Acta Aerodynamica Sinica，2005，23（2）：167－172.（in Chinese）王剛，梁新剛.羊角渦結(jié)構(gòu)及其特性［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2005，23（2）：167－172.

［16］Wu Jiezhi.Introduction to vorticity and vortex dynamics［M］.Beijing：Higher Education Press，1993.（in Chinese）吳介之.渦動(dòng)力學(xué)引論［M］.北京：高等教育出版社，1993.

［17］Jeong J H，Hussain F.On the identification of a vortex［J］.J.Fluid Mech.，1995，（285）：69－94.