謝 龍,靳思宇,2,王玉璋,于建國(guó)
(1.上海交通大學(xué)動(dòng)力機(jī)械與工程教育部重點(diǎn)試驗(yàn)室,上海 200240;2.中國(guó)船舶重工集團(tuán)公司第704研究所,上海 200031)
彎管流動(dòng)受到諸多因素的影響,尤其是在彎管前端存在閥門(mén)時(shí),彎管內(nèi)的流場(chǎng)呈現(xiàn)出十分復(fù)雜的流動(dòng)特性,如在管內(nèi)壁附近形成二次流動(dòng)分離區(qū)、管外側(cè)產(chǎn)生大尺度渦街,這些現(xiàn)象不利于流體輸運(yùn)的平穩(wěn)性能,造成流體總壓和能量的損失,并且特殊的流動(dòng)結(jié)構(gòu)還對(duì)管道噪聲特性有著重要影響。
彎管流動(dòng)的研究一直受到內(nèi)流研究者的關(guān)注[1-7]。文獻(xiàn)[2-3]中分別采用基于重整化群方法提出的RNG k-ε湍流模型和大渦模擬(Large Eddy Simulation,LES)對(duì)彎管內(nèi)部流場(chǎng)進(jìn)行數(shù)值模擬,揭示了管道橫截面上的二次流動(dòng)及沿流向的壓力分布。文獻(xiàn)[4-6]中分別采用旋轉(zhuǎn)探針技術(shù)、激光多普勒測(cè)速儀、CTA熱線風(fēng)速儀和五孔探針對(duì)不同曲率直徑比Rc/D(Rc、D分別代表曲率平均半徑和彎管的水動(dòng)力直徑)的圓形或方形截面彎管內(nèi)部流場(chǎng)進(jìn)行了試驗(yàn)研究。已有文獻(xiàn)的研究主要關(guān)注無(wú)閥門(mén)、大曲率直徑比90°彎管的流動(dòng)特性,對(duì)彎管前段存在閥門(mén),或Rc/D<1.0的圓截面彎管流動(dòng),無(wú)論試驗(yàn)還是理論研究都進(jìn)行得比較少,而這類(lèi)彎管在工業(yè)中應(yīng)用卻十分廣泛,因此,已有研究對(duì)工程應(yīng)用貢獻(xiàn)有限。
針對(duì)這種情況,結(jié)合某工程實(shí)際問(wèn)題,對(duì)前端配置蝶閥,Rc/D=0.75、D=200mm的90°圓截面彎管內(nèi)部流場(chǎng)進(jìn)行了PIV測(cè)量及POD分析。研究結(jié)果對(duì)進(jìn)一步揭示閥體后彎管內(nèi)流動(dòng)特性具有重要參考價(jià)值;為以后繼續(xù)研究管道噪聲與流場(chǎng)渦結(jié)構(gòu)特性等之間的關(guān)系打下了基礎(chǔ),對(duì)于彎管降噪的研究具有重要意義。
試驗(yàn)裝置主要包括管路系統(tǒng)、試驗(yàn)段和PIV測(cè)量系統(tǒng)[7]。管路系統(tǒng)為試驗(yàn)搭建的低速循環(huán)水槽。試驗(yàn)工質(zhì)水由磁力泵驅(qū)動(dòng),通過(guò)變頻儀調(diào)節(jié)水泵轉(zhuǎn)速,以達(dá)到對(duì)試驗(yàn)段流速的控制。循環(huán)水從水箱引流,由磁力泵驅(qū)動(dòng)并經(jīng)過(guò)整流進(jìn)入試驗(yàn)測(cè)量段,通過(guò)試驗(yàn)段后,用電磁流量計(jì)測(cè)量流速,控制流量為226.08~452.16m3/h,對(duì)應(yīng)于來(lái)流速度為2~4m/s。本次試驗(yàn)工質(zhì)水含有大量粒子,試驗(yàn)前經(jīng)測(cè)試,粒子具有良好的跟隨性、光散射性且分布均勻,滿足示蹤要求,故試驗(yàn)中不需要額外添加示蹤粒子。
試驗(yàn)段結(jié)構(gòu)如圖1所示。其中,Rc/D=0.75,D=200mm。在試驗(yàn)段前端,安裝控制開(kāi)度的一次偏心蝶閥(其閥桿軸心偏離了蝶板中心,從而使得蝶板上下端不再成為回轉(zhuǎn)軸心),蝶閥偏心距30mm,閥桿軸心與彎管入口平面距離L=112mm,試驗(yàn)中蝶閥保持全開(kāi)的開(kāi)度。均勻來(lái)流在蝶閥以及彎管的約束作用下產(chǎn)生復(fù)雜流動(dòng),試驗(yàn)測(cè)量區(qū)取在流速變化劇烈的彎管段,如圖1中陰影部分所示。以蝶閥所處位置為參考將測(cè)量區(qū)分為外側(cè)區(qū)域和內(nèi)側(cè)區(qū),取區(qū)域間的交界線(虛線)為特征線。試驗(yàn)段有機(jī)玻璃窗口處于彎管上部以及外側(cè)面,分別用于拍攝與激光照明。圖2為試驗(yàn)測(cè)量區(qū)域截面示意圖,試驗(yàn)測(cè)量段的拍攝截面以管中心為基準(zhǔn)面,選取在z向彼此相距10mm的9個(gè)水平測(cè)試截面。
高分辨率PIV測(cè)量系統(tǒng)基于Dantec Dynamics公司的二維粒子圖像測(cè)速系統(tǒng)構(gòu)建。主要組成為:雙諧振脈沖式Nd:YAG激光器、高分辨率工業(yè)相機(jī)FlowSense 4MMKII(分辨率為2048pixel×2048pixel)、同步時(shí)序控制器Timer Box 80N75及PIV應(yīng)用軟件。其中,激光器片光源鏡頭及相機(jī)通過(guò)二維坐標(biāo)架實(shí)現(xiàn)縱向位移的精確控制,相機(jī)工作頻率7.4Hz,跨幀延時(shí)100μs,所拍攝試驗(yàn)段圖像與實(shí)際尺寸之間的轉(zhuǎn)換通過(guò)拍攝標(biāo)尺圖像進(jìn)行圖像標(biāo)定來(lái)完成?;ハ嚓P(guān)圖像分析采用先進(jìn)的多重網(wǎng)格迭代計(jì)算、多重自適應(yīng)變形窗口算法、判讀區(qū)交疊查詢及互相關(guān)峰值確定的高斯擬合方法,計(jì)算結(jié)果可達(dá)±0.1pixel的亞像素精度[8],并對(duì)處理出的矢量場(chǎng)進(jìn)行進(jìn)一步的相干性過(guò)濾及局部有效性判定,所分析出的速度場(chǎng)真實(shí)反映了高分辨率流動(dòng)瞬態(tài)圖像所記錄的實(shí)際流動(dòng)信息。
圖1 試驗(yàn)段彎管結(jié)構(gòu)示意圖Fig.1 Schematic diagram of the bend of measured section
圖2 拍攝截面Fig.2 Shooting section
圖3為閥體后彎管流場(chǎng)的瞬態(tài)速度矢量圖及特征線上的速度分布。流動(dòng)的典型特征為流體在彎管圓弧內(nèi)壁的作用下順時(shí)針轉(zhuǎn)向與來(lái)流垂直的方向。因蝶閥對(duì)來(lái)流的有效阻滯作用面積較小,閥體后形成的分離區(qū)內(nèi)速度與來(lái)流流速相差不大,從而不能造成強(qiáng)烈的剪切流動(dòng)并形成彎管段內(nèi)可以觀測(cè)的較大尺度的渦旋結(jié)構(gòu)。內(nèi)側(cè)壁面渦流的生成因?yàn)榇蟪叨染鶆騺?lái)流的存在而受到抑制,因此圖3(a)中彎管內(nèi)側(cè)速度矢量分布均勻;外側(cè)矢量方向沿流向存在不規(guī)則擺動(dòng),且幅度逐漸變小直至最后出口段所有矢量互相平行,對(duì)比已有彎管研究結(jié)果[1],可以推論外側(cè)有小尺度的渦旋從蝶閥后的分離流中生成并脫落。利用瞬態(tài)速度場(chǎng)減去時(shí)均速度場(chǎng)得到脈動(dòng)速度場(chǎng),提取脈動(dòng)速度場(chǎng)中蝶閥所處位置流線上的速度分布,如圖3(b)所示,彎管外側(cè)流體沿流向存在逐漸減弱的速度起伏波動(dòng),表明該特征線上存在沿流動(dòng)方向逐漸減弱的渦列。
圖4為時(shí)均速度矢量圖及以速度云圖為背景的流線圖。彎管內(nèi)側(cè)為流動(dòng)高速區(qū)[1],原因是彎管段存在由外側(cè)指向內(nèi)側(cè)的壓力梯度,在此壓力梯度作用下在垂直于流向的橫截面上產(chǎn)生一組上下反對(duì)稱的渦流二次流,且渦核位置隨著流動(dòng)的發(fā)展逐漸偏向彎管內(nèi)側(cè),截面上下兩側(cè)近壁面流動(dòng)由外側(cè)指向內(nèi)側(cè),中間部分由內(nèi)側(cè)指向外側(cè)[2-3],該二次流與彎管主流的矢量和形成了內(nèi)側(cè)高速區(qū),且流向以小角度向彎管外側(cè)偏離該處彎管同心圓的切線方向。彎管外側(cè)矢量分布均勻,說(shuō)明瞬態(tài)場(chǎng)中觀察到的渦街的時(shí)均流動(dòng)效果為0??梢缘贸鼋Y(jié)論,閥體后擾流在彎管外側(cè)區(qū)域產(chǎn)生的渦流結(jié)構(gòu)具有宏觀不可測(cè)性,但閥體后擾流在向下游輸運(yùn)的過(guò)程中會(huì)引起當(dāng)?shù)厮俣让}動(dòng),而所有脈動(dòng)信息中渦結(jié)構(gòu)對(duì)流場(chǎng)貢獻(xiàn)權(quán)重大于流動(dòng)噪聲,因此可以從速度脈動(dòng)中提取旋渦結(jié)構(gòu)的信息。
圖3 瞬態(tài)流場(chǎng)Fig.3 Transient flow field
圖4 時(shí)均流場(chǎng)Fig.4 Time-averaged flow field
本征正交分解(POD)基于主元分析(PCA)思想,根據(jù)Karhunen-Loève展開(kāi)將流場(chǎng)進(jìn)行分解,按能量比計(jì)算出各階模態(tài),識(shí)別流動(dòng)中的大尺度能級(jí)結(jié)構(gòu)。POD之于PIV相當(dāng)于LES之于CFD,POD對(duì)流動(dòng)系統(tǒng)的識(shí)別旨在獲得多維度系統(tǒng)的低維度近似描述,準(zhǔn)確把握住主導(dǎo)的流動(dòng)結(jié)構(gòu),對(duì)復(fù)雜流動(dòng)進(jìn)行合理的降階研究。本征正交分解已經(jīng)被廣泛用于各種湍流速度場(chǎng)、渦量場(chǎng)、壓力場(chǎng)等的研究中[9-11]。本次試驗(yàn)采用由Sirovich提出的快照POD法(snapshot POD)[12],將速度場(chǎng)分解為時(shí)均速度場(chǎng)和脈動(dòng)速度場(chǎng),即u(x)=ū+u′(x),然后再對(duì)脈動(dòng)速度場(chǎng)進(jìn)行本征正交分解,全部流動(dòng)信息便分解為時(shí)均流動(dòng)(0階模態(tài))與脈動(dòng)流動(dòng)(分解為1,2,…階模態(tài))。給出0mm截面位置在3m/s來(lái)流速度工況下的POD分析結(jié)果。
圖5為脈動(dòng)速度場(chǎng)各模態(tài)能量分布(不包含0階模態(tài))。彎管速度脈動(dòng)場(chǎng)的1階模態(tài)占據(jù)了脈動(dòng)動(dòng)能的14.6%,在脈動(dòng)流場(chǎng)中占據(jù)支配地位,適于用來(lái)描述典型的脈動(dòng);2~3階模態(tài)的能量貢獻(xiàn)次之,所占脈動(dòng)動(dòng)能比例分別為9.2%和6.8%,可認(rèn)為這些模態(tài)對(duì)脈動(dòng)流場(chǎng)起了重要作用。而后曲線的下降及后續(xù)的平整意味著高階的模態(tài)對(duì)總動(dòng)能的貢獻(xiàn)很低,可以簡(jiǎn)單描述為流動(dòng)噪聲,如短暫的外部擾動(dòng)、大尺度渦引起的高頻小尺度結(jié)構(gòu)或者耗散渦等。由于彎管中的渦街是脈動(dòng)流的載體,因此脈動(dòng)場(chǎng)的低階模態(tài)中包含了彎管中的渦結(jié)構(gòu)信息。
圖5 POD模態(tài)能量分布Fig.5 POD modes energy distribution
1~3階模態(tài)包含了閥體后擾流在彎管流場(chǎng)中產(chǎn)生的旋渦結(jié)構(gòu),除了簡(jiǎn)單的流線識(shí)別外,渦的識(shí)別方法有Calilean分解方法、LES分解方法、Reynolds分解方法、小渦分解方法、渦量方法和旋渦強(qiáng)度(swirling strength)方法[13]等。圖6給出了1~4階模態(tài)流動(dòng)相干結(jié)構(gòu)以渦量為背景的速度矢量圖,其中表征時(shí)均流動(dòng)的0階模態(tài)前文已給出,圖6中可見(jiàn)4階模態(tài)已不具備明顯的流動(dòng)結(jié)構(gòu)特征,再次說(shuō)明4階及其以上模態(tài)表征流動(dòng)噪聲。1~3階模態(tài)中渦列中渦核間距離較為恒定,渦核出現(xiàn)位置近似處于蝶閥安裝位置所在流線,且渦量數(shù)值分布為一正一負(fù)交替排列,說(shuō)明沿流線旋渦旋轉(zhuǎn)方向相反,這表現(xiàn)為典型的卡門(mén)渦街特征。過(guò)程為:在蝶閥前緣分離的剪切層,由于流體的粘性剪切作用使得剪切層卷起形成旋渦,而閥體兩側(cè)剪切層內(nèi)渦對(duì)的相互吸引和相互作用,最終導(dǎo)致旋渦從閥體兩側(cè)剝離脫落,形成卡門(mén)渦街,以較高頻率脫落的渦隨流體進(jìn)入彎管段,彎管橫截面上的渦流二次流將卡門(mén)渦向外側(cè)擠壓,彎管外側(cè)輪廓對(duì)渦流的向外側(cè)發(fā)展產(chǎn)生束縛,同時(shí)旋渦本身從剪切層脫落以后相互吸引、相互作用,3種原因共同作用下使得旋渦的運(yùn)動(dòng)軌跡不斷向閥體所在位置流線靠近,因此彎管中渦街的各個(gè)渦出現(xiàn)的位置近似在同一流線上;彎管橫截面上二次流向外側(cè)的偏移擠壓及外側(cè)壁面效應(yīng)使渦外圍能量隨著流線逐漸削弱,同時(shí)隨著卡門(mén)渦向下游的不斷輸運(yùn),渦的旋轉(zhuǎn)與周?chē)黧w產(chǎn)生摩擦也耗散了部分渦能量,因此可清晰觀測(cè)到渦量沿流線逐漸減弱的過(guò)程。根據(jù)渦量的定義ω=▽×u,渦量方法除了識(shí)別出流場(chǎng)中的渦核,同時(shí)也反映了彎管流場(chǎng)中大量剪切運(yùn)動(dòng)信息。除脈動(dòng)場(chǎng)中的渦結(jié)構(gòu)以外,模態(tài)1中可以觀測(cè)到由彎管內(nèi)側(cè)指向外側(cè)的二次流動(dòng)[2-3];模態(tài)2中在彎管出口段還觀測(cè)到大尺度不規(guī)則回流區(qū),該回流由出口管外壁對(duì)出口斜射流的反作用逆壓梯度引起;在彎管內(nèi)壁出口位置也存在小尺度的渦流,這是由于內(nèi)壁出口位置的局部低速區(qū)壓力較高,與高速區(qū)形成的低壓之間存在壓力梯度,從而引起局部回流;流場(chǎng)中渦流區(qū)與二次流區(qū)以外為大范圍的低速區(qū)。
圖6 低階POD模態(tài)相干結(jié)構(gòu)Fig.6 Coherent structures of low POD modes
將流動(dòng)狀態(tài)投影到各POD模態(tài),可以用若干基模態(tài)的線性組合來(lái)表示流場(chǎng),即POD模態(tài)重構(gòu),模態(tài)重構(gòu)可以還原主要流動(dòng)結(jié)構(gòu)。模態(tài)在重構(gòu)具體流場(chǎng)時(shí)所承擔(dān)的主要意義所占的比重由其對(duì)應(yīng)的POD系數(shù)表征,圖7給出了所研究時(shí)間范圍內(nèi)模態(tài)1的POD系數(shù)的時(shí)變圖。POD系數(shù)隨時(shí)間以類(lèi)似正弦曲線的方式正負(fù)交替波動(dòng)代表了該模態(tài)對(duì)應(yīng)的流動(dòng)結(jié)構(gòu)隨時(shí)間以正反作用交替的方式影響流動(dòng)特性。由于該模態(tài)代表的是渦結(jié)構(gòu),因此其在流場(chǎng)中的作用體現(xiàn)為渦結(jié)構(gòu)的順時(shí)針與逆時(shí)針旋向交替。這是由于卡門(mén)渦街中渦旋總是以成對(duì)的旋轉(zhuǎn)方向相反的反對(duì)稱渦形式交錯(cuò)排列,流場(chǎng)中的固定空間點(diǎn)隨時(shí)間受到不同旋向渦的作用。本次試驗(yàn)采用的高分辨率相機(jī)沒(méi)有很好的時(shí)間分辨特性,只采樣到流場(chǎng)中渦流旋向真實(shí)交替過(guò)程的幾個(gè)片段,因此觀測(cè)到的旋向交替過(guò)程中表現(xiàn)為正負(fù)值之間的大幅度波動(dòng),故試驗(yàn)不能反映實(shí)際流動(dòng)中渦流旋向的交替頻率,即卡門(mén)渦街脫落頻率,但可以預(yù)計(jì)渦流旋向交替隨時(shí)間的變化形態(tài)表現(xiàn)為一個(gè)周期很短的正弦函數(shù)[10],且該頻率遠(yuǎn)大于試驗(yàn)觀測(cè)到的頻率。
圖7 模態(tài)1POD系數(shù)時(shí)變圖Fig.7 Time-varying graph of POD coefficient value of mode 1
1~3階模態(tài)中彎管外側(cè)渦結(jié)構(gòu)的能量逐漸減弱至很低的水平,4階及以上模態(tài)已不具備明顯流動(dòng)結(jié)構(gòu)特征,說(shuō)明這3個(gè)模態(tài)中包含了閥體后彎管流場(chǎng)中由蝶閥擾流引起的大部分渦結(jié)構(gòu)信息。排除時(shí)均流動(dòng)與流動(dòng)噪聲,利用1~3階模態(tài)重構(gòu)流場(chǎng)可以還原彎管脈動(dòng)場(chǎng)中由閥體后卡門(mén)渦引起的流動(dòng)結(jié)構(gòu),得到無(wú)噪聲情況下脈動(dòng)流場(chǎng)的合理壓縮。由于卡門(mén)渦以很高的頻率從閥體后脫落,可以預(yù)計(jì)重構(gòu)出的渦街中渦的旋向與強(qiáng)度隨時(shí)間變化存在較大差異。針對(duì)還原的渦結(jié)構(gòu)信息進(jìn)行進(jìn)一步分析可以揭示渦的流動(dòng)特性,利用旋渦強(qiáng)度分布識(shí)別重構(gòu)脈動(dòng)場(chǎng)中的渦結(jié)構(gòu),用笛卡爾坐標(biāo)下速度梯度張量的復(fù)特征值虛部分析能夠使渦結(jié)構(gòu)顯示出來(lái),并量化局部旋渦運(yùn)動(dòng)的強(qiáng)度,這一復(fù)特征值虛部定義為局部旋渦運(yùn)動(dòng)的旋渦強(qiáng)度。正特征值意味著流場(chǎng)中存在剪切流動(dòng),而不存在渦結(jié)構(gòu);負(fù)的旋渦強(qiáng)度表明該處有渦結(jié)構(gòu)存在,而局部最小值則可以用來(lái)識(shí)別渦核的位置。圖8給出彎管脈動(dòng)流場(chǎng)信息POD重構(gòu)出的典型渦街結(jié)構(gòu),其中計(jì)算結(jié)果只返回流場(chǎng)中的負(fù)特征值,而將正特征值全部設(shè)為0(紅色區(qū)域),以排除剪切作用的影響。圖中可以清晰觀察到處于脫落過(guò)程中的卡門(mén)渦,渦列中相鄰渦旋向交替、渦核位置近似處于同一流線、脫落過(guò)程中渦能量逐漸衰減。此外,“淹沒(méi)”在下游主流中的旋渦也能從云圖中觀察到。
圖8 旋渦強(qiáng)度分布Fig.8 Distribution of swirling strength
進(jìn)一步的研究可以基于利用脈動(dòng)主導(dǎo)模態(tài)重構(gòu)的脈動(dòng)場(chǎng),對(duì)流場(chǎng)的降階研究可以降低問(wèn)題的復(fù)雜度,增強(qiáng)研究的針對(duì)性,且低階模型具有較高精度。對(duì)固定來(lái)流流速、變觀察截面的研究表明:由于閥體圓形結(jié)構(gòu)在靠近上下端時(shí)尺寸逐漸變小,閥體對(duì)來(lái)流的有效阻滯作用減小,故相對(duì)于中間截面,上下兩側(cè)截面位置上閥體后渦流強(qiáng)度減小且其分離點(diǎn)更靠前,并以較小能量匯入主流,在彎管前的直管段內(nèi)其能量已被大部分消耗,因此在彎管中僅有微弱的表現(xiàn);在彎管橫截面上渦流二次流的擠壓下彎管中渦出現(xiàn)的位置在徑向產(chǎn)生較大的指向外輪廓的偏移;同時(shí)觀察到兩側(cè)截面上渦脫落頻率小于中間截面,說(shuō)明蝶閥后卡門(mén)渦街脫落頻率與該截面上蝶閥長(zhǎng)度正相關(guān)。對(duì)固定觀察截面、變來(lái)流流速的研究表明:試驗(yàn)中阻流體迎面寬度恒定不變,卡門(mén)渦街脫落頻率隨來(lái)流流速成正比線性變化,因此在所研究流速范圍以內(nèi),流場(chǎng)中的渦結(jié)構(gòu)不隨流速產(chǎn)生顯著改變,隨著流速的變化旋渦強(qiáng)度做近似線性變化。綜上所述,彎管中的脈動(dòng)大尺度結(jié)構(gòu)為位于彎管外側(cè)壁面附近的渦街,因此如果能改進(jìn)該區(qū)域的脈動(dòng)流,使該區(qū)域大尺度結(jié)構(gòu)消除或弱化,將有利于流體輸運(yùn)的平穩(wěn)性能在整體上的提高,進(jìn)而有效降低流動(dòng)噪聲。改進(jìn)方式包括調(diào)整蝶閥安裝位置;將閥體背離來(lái)流部分補(bǔ)齊,使其表面接近流線型,從而消除旋渦的脫落;改進(jìn)彎管外側(cè)面結(jié)構(gòu)或配合脈動(dòng)緩沖器使用等。
(1)流體流經(jīng)蝶閥在彎管中形成卡門(mén)渦街,渦流結(jié)構(gòu)宏觀不可測(cè),可以從脈動(dòng)速度場(chǎng)中提取旋渦結(jié)構(gòu)信息;
(2)彎管中渦結(jié)構(gòu)信息包含于低階POD模態(tài),低階模態(tài)流動(dòng)相干結(jié)構(gòu)中渦列處于蝶閥安裝位置所在流線,渦核間距恒定,相鄰渦旋向相反,渦脫落過(guò)程中渦能量逐漸衰減;
(3)脈動(dòng)場(chǎng)的模態(tài)重構(gòu)可以準(zhǔn)確還原彎管中的渦結(jié)構(gòu),彎管中間截面渦流強(qiáng)度大于上下兩側(cè)截面,蝶閥后卡門(mén)渦街脫落頻率與該截面上蝶閥長(zhǎng)度正相關(guān),并隨來(lái)流流速做近似線性變化。消除或減弱外側(cè)壁面附近的大尺度渦結(jié)構(gòu)將有利于流體輸運(yùn)平穩(wěn)性能在整體上的提高,并降低流動(dòng)噪聲。
[1] 江山,張京,偉吳崇,等.基于FLUENT的90°圓形彎管內(nèi)部流場(chǎng)分析[J].中國(guó)艦船研究,2008,3(1):37-41.
[2] 丁玨,翁培奮.90°彎管內(nèi)流動(dòng)的理論模型及流動(dòng)特性的數(shù)值研究[J].計(jì)算力學(xué)學(xué)報(bào),2004,21(3):314-321.
[3] 梁開(kāi)洪,曹樹(shù)良,陳炎,等.入流角對(duì)圓截面90°彎管內(nèi)高雷諾數(shù)流動(dòng)的影響[J].清華大學(xué)學(xué)報(bào),2009,49(12):1971-1975.
[4] SUDO K,SUMIDA M,HIBARA H.Experimental investigation on turbulent flow in a circular-sectioned 90-degree bend[J].Experiments in Fluids,1998,25(1):42-49.
[5] WILLIAMY,BENNY K.Experimental investigation of dilute turbulent particulate flow inside a curved 90°bend[J].Chemical Engineering Science,2006,61(11):3593-3601.
[6] 尚虹,王尚錦,席光,等.90°圓截面彎管內(nèi)三維紊流場(chǎng)試驗(yàn)研究[J].航空動(dòng)力學(xué)報(bào),1994,9(3):263-266.
[7] 謝龍,靳思宇,于建國(guó),等.閥體后90°圓形彎管內(nèi)部流場(chǎng)PIV分析[J].上海交通大學(xué)學(xué)報(bào),2011,45(9):1395-1399.
[8] 萬(wàn)津津,施鎏鎏,余俊,等.貼壁方柱湍流場(chǎng)TR-PIV實(shí)驗(yàn)研究[J].實(shí)驗(yàn)流體力學(xué),2009,23(2):31-35.
[9] 陽(yáng)祥,丁鵬,李增耀,等.基于POD方法的湍流大尺度信息提取與分析[J].工程熱物理學(xué)報(bào),2010,31(6):1019-1022.
[10]SHI Liu-liu,LIUI Ying-zheng,WAN Jin-jin.Influence of wall proximity on characteristics of wake behind a square cylinder:PIV measurements and POD analysis[J].Experimental Thermal and Fluid Science,2010,34(1):28-36.
[11]O'DONNELL B J,HELENBROOK B T.Proper orthogonal decomposition and incompressible flow:An application to particle modeling[J].Computers and Fluids,2007,36(7):1174-1186.
[12]SIROVICH L.Turbulence and the dynamics of coherent structures[J].Quarterly of Applied Mathematics,1987,45(3):561-590.
[13]ADRIAN M J,CHRISTENSEN K T,LIU Z C.Analysis and interpretation of instantaneous turbulent velocity fields[J].Experiments in Fluids,2000,29(3):275-290.