陳小虎, 陳 方,*, 劉 洪, 沙 莎, 逯雪鈴, 張慶兵

(1. 上海交通大學(xué) 航空航天學(xué)院, 上海 200240; 2. 北京電子工程總體研究所, 北京 100854)

超聲速混合層中PIV粒子的湍流變動(dòng)作用研究

陳小虎1, 陳 方1,*, 劉 洪1, 沙 莎2, 逯雪鈴2, 張慶兵2

(1. 上海交通大學(xué) 航空航天學(xué)院, 上海 200240; 2. 北京電子工程總體研究所, 北京 100854)

對(duì)二維超聲速氣固兩相混合層進(jìn)行雙向耦合，研究了粒子圖像測(cè)速技術(shù)(PIV)中示蹤粒子對(duì)超聲速混合層的湍流變動(dòng)作用。超聲速氣固兩相混合層的氣相采用大渦模擬，離散相采用拉格朗日顆粒軌道模型求解。結(jié)果表明：與無負(fù)載示蹤粒子時(shí)的超聲速混合層相比，小Stokes數(shù)示蹤粒子在超聲速混合層中的布撒減弱了流向湍流，而強(qiáng)化了法向湍流，使雷諾應(yīng)力峰值增大了9.68%；大Stokes數(shù)示蹤粒子對(duì)混合層的湍流脈動(dòng)起到了一定的削弱作用，最大雷諾應(yīng)力值只有無負(fù)載時(shí)的41.74%。大質(zhì)量載荷時(shí)，大量示蹤粒子的運(yùn)動(dòng)尾跡抹平了部分法向速度脈動(dòng)，使最大法向速度脈動(dòng)只有無負(fù)載粒子時(shí)的38.63%；中等質(zhì)量載荷時(shí)，超聲速混合層的法向速度脈動(dòng)和雷諾應(yīng)力峰值與無負(fù)載粒子時(shí)相近；而小質(zhì)量載荷時(shí)，超聲速混合層中心線及其附近的法向速度脈動(dòng)得到較小的增強(qiáng)，而最大流向速度脈動(dòng)卻被削弱了19.29%。小Stokes數(shù)和中等質(zhì)量載荷示蹤粒子對(duì)原始無負(fù)載粒子時(shí)的流場(chǎng)影響相對(duì)較小，研究結(jié)論對(duì)高速流動(dòng)PIV測(cè)試有著重要的參考價(jià)值。

PIV；超聲速混合層；湍流變動(dòng)；Stokes數(shù)；質(zhì)量載荷

0 引 言

PIV測(cè)試技術(shù)作為一種非接觸光學(xué)測(cè)量方式，在流場(chǎng)診斷測(cè)量的應(yīng)用中展現(xiàn)出強(qiáng)大的功能。示蹤粒子Stokes(St)數(shù)表征粒子的松弛時(shí)間尺度與流場(chǎng)特征時(shí)間尺度之比，PIV示蹤粒子直徑(dp)越小，St數(shù)越小。質(zhì)量載荷Φm衡量單位質(zhì)量氣相中離散相粒子的質(zhì)量百分?jǐn)?shù)。根據(jù)Mie氏散射理論，粒徑越小散射強(qiáng)度越弱。此時(shí)，需要通過增大PIV示蹤粒子的質(zhì)量載荷以獲得到更多的流場(chǎng)信息。但相較于無負(fù)載粒子時(shí)的被測(cè)流場(chǎng)，較大的Φm會(huì)給原始流場(chǎng)帶來劇烈的干擾，影響PIV實(shí)驗(yàn)的可靠性[1]。因此，在高速復(fù)雜流場(chǎng)PIV測(cè)量中，探究粒子St數(shù)和Φm對(duì)超聲速流場(chǎng)的湍流變動(dòng)顯得尤為重要。

隨著PIV等實(shí)驗(yàn)測(cè)試和數(shù)值模擬方法的快速發(fā)展，研究者對(duì)超聲速氣固兩相流動(dòng)的研究也越發(fā)深入。特別是在較大Φm的氣固兩相流動(dòng)中，粒子對(duì)氣相湍流變動(dòng)的研究成為多相流領(lǐng)域的研究熱點(diǎn)之一。Gore, Crowe和Kenning[2-4]總結(jié)了之前的研究，認(rèn)為粒子直徑與流體湍流尺度的比值(dp/l)可以作為增強(qiáng)或削弱湍流的判斷準(zhǔn)則。Hetsroni[5]研究了粒子雷諾數(shù)(Rep)對(duì)湍流變動(dòng)的影響。Rep大于400時(shí)，湍流增強(qiáng)；Rep小于400時(shí)，湍流減弱。Elghobashi[6-7]認(rèn)為質(zhì)量攜帶率和顆粒尺寸會(huì)影響連續(xù)相湍流。當(dāng)顆粒體積分?jǐn)?shù)處于[10-6,10-3]之間，St數(shù)大于1時(shí)，顆粒使湍流增強(qiáng)，反之則減弱。Savolainen和Karvinen[8]以及Lain和Sommerfeld[9-10]采用Lagrange方法研究了豎直向上和水平管流中的粒子對(duì)湍流的影響并進(jìn)行了對(duì)比試驗(yàn)的驗(yàn)證，認(rèn)為粒子使低速區(qū)域的湍流得到增強(qiáng)，高速流動(dòng)區(qū)域的湍流被削弱。Druzhinin[11]和Ling[12]采用直接數(shù)值模擬研究了氣固雙向耦合對(duì)湍流的影響。林建忠等[13-14]研究表明，顆粒的濃度與St數(shù)對(duì)流場(chǎng)的影響不可以忽略，但空間發(fā)展的超聲速流動(dòng)中是否有相似規(guī)律未做說明。樊建人等[15]用擬譜法對(duì)時(shí)間發(fā)展的氣固混合層雙向耦合，顆粒使混合層大渦卷起延遲，且質(zhì)量攜帶率越大氣相湍流能量衰減越厲害。王兵等[16]通過多普勒顆粒分析儀實(shí)驗(yàn)研究了粒徑和氣載比對(duì)圓湍射流的湍流調(diào)制。劉偉等[17-18]研究發(fā)現(xiàn)對(duì)流馬赫數(shù)0.8時(shí)，加入粒子會(huì)使流場(chǎng)中的小激波結(jié)構(gòu)消失，但未解釋?duì)祄對(duì)流場(chǎng)的反作用。

基于以上研究現(xiàn)狀，現(xiàn)有研究很少有以實(shí)際高速流動(dòng)PIV測(cè)試為研究對(duì)象,且對(duì)空間發(fā)展的超聲速氣固混合層中示蹤粒子影響湍流變動(dòng)的判斷準(zhǔn)則以及作用機(jī)理研究仍不充分。因此本文著重分析不同St數(shù)和Φm的PIV示蹤粒子對(duì)空間發(fā)展的二維超聲速混合層的湍流變動(dòng)作用，并與無負(fù)載粒子時(shí)的流場(chǎng)參數(shù)對(duì)比，預(yù)期為PIV測(cè)試中示蹤粒子粒徑和布撒濃度的選擇提供參考依據(jù)。

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

1.1 氣相控制方程

為研究超聲速混合層中PIV示蹤粒子的湍流變動(dòng)作用，對(duì)氣相求解Navier-Stokes(N-S)方程。其中，氣相物理量用下標(biāo)“g”表示，離散相物理量用下標(biāo)“p”表示，如下：

分子粘性引起的張量σij和亞格子應(yīng)力τijg分別表達(dá)為：

本文采用Smagorinsky-Lilly所提出的亞格子應(yīng)力模型。Smagorinsky常數(shù)并不是一個(gè)固定值，Yoshizawa 和Horiuti[19]認(rèn)為較大的Smagorinsky常數(shù)會(huì)導(dǎo)致耗散的增強(qiáng)，在均勻湍流中取Cs=0.2而剪切湍流中取Cs=0.065是合適的。氣相在空間上的數(shù)值離散是二階迎風(fēng)格式，時(shí)間上的數(shù)值離散是二階隱式推進(jìn)。

1.2 固相控制方程

單個(gè)粒子在拉格朗日框架下的運(yùn)動(dòng)方程為:

式中：xip、uip、mip、Tip和cip分別表示第i個(gè)粒子的位置坐標(biāo)、速度、質(zhì)量、溫度和比熱；Figp、Fipp和Qip分別代表第i個(gè)粒子受到的氣動(dòng)力、粒子間相互作用力和傳熱量。因?yàn)楸疚闹辛Ｗ拥牟既鲶w積分?jǐn)?shù)均小于1×10-3，粒子之間的相互碰撞不明顯，所以不考慮離散相粒子之間的相互碰撞效應(yīng)，故Fipp=0。

在PIV試驗(yàn)中，示蹤粒子在運(yùn)動(dòng)的過程中受到諸多的作用力，如虛假質(zhì)量力、馬格努斯力和熱泳力等。因本文所涉及的粒子均在微米量級(jí)及以下，運(yùn)動(dòng)阻力對(duì)粒子運(yùn)動(dòng)的影響較其他作用力大得多。粒子運(yùn)動(dòng)阻力的計(jì)算，如式(9)和(10)所示。離散相在空間各個(gè)時(shí)刻和位置上的彌散位置由離散隨機(jī)游走模型確定。

2 計(jì)算工況與驗(yàn)證

2.1 氣相計(jì)算參數(shù)

本文對(duì)空間發(fā)展的超聲速氣固兩相混合層進(jìn)行雙向耦合，探究了超聲速混合層中PIV示蹤粒子的湍流變動(dòng)作用。選取Mac=0.20的Goebel-Dutton[20]可壓縮超聲速混合層實(shí)驗(yàn)Case1，來流條件如表1所示。對(duì)流馬赫數(shù)Mac定義為：

表1 計(jì)算參數(shù)Table 1 Computational parameters

本文選取的計(jì)算區(qū)域[21]為Lx×Ly=0.5m×0.04m，對(duì)計(jì)算域進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分并加密。采用3套網(wǎng)格(1001×201、1601×251、2001×301)計(jì)算超聲速混合層，并進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。1001×201的網(wǎng)格與1601×251的網(wǎng)格相比，混合層厚度在發(fā)展前期相差較大。而1601×251的網(wǎng)格與2001×301的網(wǎng)格相比，混合層厚度相差不大。為減少計(jì)算量選取1601×251的網(wǎng)格進(jìn)行計(jì)算。為了使計(jì)算符合物理事實(shí)，以雙曲正切型速度剖面初始化入口速度u(y)。設(shè)置上下邊界為無反射邊界，以過濾邊界對(duì)混合層的影響。

式中：U1、U2為上下來流初始速度；θ0為初始時(shí)刻的邊界層厚度。根據(jù)Goebel-Dutton的實(shí)驗(yàn)數(shù)據(jù)在此取θ0=0.0004m。

2.2 離散相計(jì)算參數(shù)

選取上海交通大學(xué)多馬赫數(shù)高超聲速風(fēng)洞PIV測(cè)試[22-24]所用TiO2納米粒子作為離散相。TiO2納米粒子的材料密度為4230kg/m3，摩爾質(zhì)量為79.9g/mol。St數(shù)和Φm在兩相流中都是重要的無量綱評(píng)價(jià)參數(shù)，可以描述為：

2.3 計(jì)算參數(shù)選擇

在PIV診斷測(cè)速技術(shù)中，示蹤粒子的直徑和布撒濃度是影響測(cè)試準(zhǔn)確度的2個(gè)關(guān)鍵因素。粒徑和布撒濃度在兩相流中可以用參數(shù)St數(shù)和Φm來表征。在超聲速氣固兩相混合層中，粒徑導(dǎo)致的St數(shù)變化，直接改變示蹤粒子在大渦擬序結(jié)構(gòu)中的分布特征從而產(chǎn)生湍流變動(dòng)作用。而示蹤粒子Φm的增大會(huì)加劇粒子對(duì)混合層的影響。當(dāng)Φm增大到一定程度之后，原始的連續(xù)相流場(chǎng)結(jié)構(gòu)已經(jīng)發(fā)生變化，此時(shí)的測(cè)量結(jié)果已經(jīng)不具有可信性。

為探究超聲速混合層中示蹤粒子St數(shù)和Φm的湍流變動(dòng)作用，選擇St=0.01、1和100，對(duì)應(yīng)dp分別為50nm、525nm和5μm，以及Φm=0.0121、0.1092和0.5507的PIV粒子作為超聲速氣固兩相混合層雙向耦合的離散相進(jìn)行數(shù)值模擬。

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

3.1 氣相計(jì)算結(jié)果驗(yàn)證

如圖1所示是超聲速混合層無量綱流向時(shí)均速度和脈動(dòng)速度曲線。通過對(duì)比發(fā)現(xiàn)，按照實(shí)驗(yàn)參數(shù)計(jì)算的超聲速混合層的無量綱流向時(shí)均速度剖面呈雙曲正切型分布，與Goebel-Dutton實(shí)驗(yàn)數(shù)據(jù)基本吻合。因?yàn)閷?shí)驗(yàn)中的湍流流動(dòng)是三維的，而本文計(jì)算的是二維超聲速混合層，不能對(duì)展向湍流進(jìn)行很好地預(yù)測(cè)，所以脈動(dòng)速度分布與實(shí)驗(yàn)值很接近，最大速度脈動(dòng)值減小了約7.33%。本文利用UDF自定義程序嵌入Fluent求解器進(jìn)行大渦模擬，能模擬出合理的超聲速混合層流場(chǎng)結(jié)構(gòu)。能很好的支持下文討論超聲速混合層中不同St數(shù)和Φm的PIV示蹤粒子的湍流變動(dòng)作用。

3.2 不同St數(shù)時(shí)的湍流變動(dòng)作用

PIV示蹤粒子的直徑影響兩相流動(dòng)的St數(shù)大小，而研究表明，不同St數(shù)的粒子在混合層大渦擬序結(jié)構(gòu)中存在傾向性的富集現(xiàn)象。如圖2所示是Φm=0.1092時(shí)，St=0.01、1和100的粒子與氣相雙向耦合下的氣相渦量場(chǎng)和粒子分布圖。粒子在某一區(qū)域的大密度富集勢(shì)必對(duì)原始的流場(chǎng)產(chǎn)生劇烈干擾，因而研究不同St數(shù)的粒子對(duì)超聲速混合層的湍流變動(dòng)具有重要意義。

Fig.1Normalizedstreamwisemeanandrootmeansquarevelocitydistribution

質(zhì)量載荷Φm=0.1092對(duì)應(yīng)的初始粒子布撒體積分?jǐn)?shù)是Φv=1×10-5,而計(jì)算發(fā)現(xiàn)，選擇以質(zhì)量載荷作為研究變量比體積分?jǐn)?shù)更具有參考意義。一定的質(zhì)量載荷下，隨著St數(shù)從0.01增大到100，粒子傾向性分布加劇。St=100時(shí)粒子跟隨性不好，無法隨渦旋轉(zhuǎn),此時(shí)的大粒徑粒子運(yùn)動(dòng)尾跡產(chǎn)生的湍流效應(yīng)比小粒徑粒子明顯，氣相流場(chǎng)會(huì)受到一定的影響。St=0.01的粒子對(duì)氣相渦量場(chǎng)的影響大于St=1和100的粒子。從氣相渦量等值線可以看到，St=0.01時(shí)粒子帶來的擾動(dòng)較大，導(dǎo)致等值線出現(xiàn)小波動(dòng)。St=1時(shí)氣相渦量等值線趨近橢圓形，與無負(fù)載粒子時(shí)充分發(fā)展的渦量場(chǎng)相似。

圖3～5分別是無量綱時(shí)均流向速度和溫度分布、無量綱流向和法向速度脈動(dòng)和無量綱雷諾應(yīng)力分布曲線。如圖3所示，3種粒徑的粒子在經(jīng)過長時(shí)間平均后都能較好地逼近無負(fù)載時(shí)混合層的時(shí)均速度。St數(shù)對(duì)時(shí)均溫度影響較大，St=100時(shí)的粒子因?yàn)榱捷^大，兩相之間的傳熱過程與小粒子相比相對(duì)漫長。大粒徑粒子使氣相的時(shí)均溫度不能很好地跟隨無負(fù)載時(shí)的混合層。

Fig.3Normalizedmeanstreamwisevelocityandtemperaturedistribution

超聲速混合層的流向速度脈動(dòng)呈現(xiàn)出無負(fù)載粒子時(shí)大于氣固雙向耦合時(shí)的情形，而不同St數(shù)粒子對(duì)法向速度脈動(dòng)有加強(qiáng)也有削弱，如圖4所示。示蹤粒子St數(shù)的大小對(duì)流向速度脈動(dòng)的影響相較于法向速度脈動(dòng)更為明顯。對(duì)于流向的速度脈動(dòng)，與無負(fù)載的超聲速混合層相比，加入粒子后摩擦帶來的流向流動(dòng)阻力增大，流向速度脈動(dòng)減弱。St=100時(shí)的速度脈動(dòng)是無負(fù)載時(shí)的88%，大于St=1和0.01時(shí)的速度脈動(dòng)。對(duì)于法向速度脈動(dòng)，與無負(fù)載的超聲速混合層相比，0

Fig.4Normalizedrootmeansquarestreamwiseandtransversevelocitydistribution

超聲速混合層的無量綱雷諾應(yīng)力的湍流變動(dòng)如圖5所示。由于dp的增大，粒子的傾向性聚集現(xiàn)象明顯，對(duì)混合層的湍流影響顯著。相同質(zhì)量載荷下，St=1和100與St=0.01時(shí)相比粒徑大而粒子數(shù)較少，布撒PIV粒子后湍流正應(yīng)力比湍流切應(yīng)力變化更明顯。對(duì)于局部的流體微團(tuán)，粒子的加入減小了湍流脈動(dòng)造成的湍流正應(yīng)力，從而使混合層中心線附近的雷諾應(yīng)力被削弱。St=100的氣固混合層的最大雷諾應(yīng)力只有無負(fù)載時(shí)的41.74%，粒子與氣體的摻混能力減弱。St=0.01的示蹤粒子因?yàn)榫哂辛己玫碾S流性，粘性引起的湍流切應(yīng)力主導(dǎo)雷諾應(yīng)力的變化。湍流脈動(dòng)帶來的流層之間的剪切增強(qiáng)，湍流切應(yīng)力隨之增大，雷諾應(yīng)力值大于無負(fù)載粒子時(shí)的混合層。小St數(shù)的示蹤粒子雖然使超聲速氣固兩相混合層中心線處雷諾應(yīng)力增大了9.68%，但此時(shí)氣固摻混能力強(qiáng)并且與大St數(shù)相比對(duì)原始流場(chǎng)影響相對(duì)較小，適合作為PIV示蹤粒子。

3.3 不同Φm時(shí)的湍流變動(dòng)作用

在超聲速氣固兩相混合層的雙向耦合中，PIV示蹤粒子的Φm大小會(huì)不同程度地影響氣相流場(chǎng)。相同St數(shù)，不同Φm的粒子與超聲速混合層雙向耦合，氣相作為載體輸運(yùn)示蹤粒子做跟隨運(yùn)動(dòng)。但是，當(dāng)Φm太大時(shí)，兩相間的動(dòng)能和內(nèi)能交換劇烈，此時(shí)示蹤粒子對(duì)氣相的湍流作用顯著并產(chǎn)生明顯的湍流變動(dòng)作用。圖6所示是不同Φm時(shí)氣固雙向耦合下的氣相渦量場(chǎng)和粒子運(yùn)動(dòng)。

從圖6可知，對(duì)于粒子St=0.01的氣固混合層，氣相渦量等值線隨Φm的增大越來越不規(guī)則。Φm=0.0121時(shí)，高速來流粒子產(chǎn)生“順流下卷”運(yùn)動(dòng)，低速來流的粒子隨流動(dòng)做“逆流上拋”運(yùn)動(dòng)。示蹤粒子此時(shí)對(duì)氣相渦量的影響微小，幾乎可以忽略。而當(dāng)質(zhì)量載荷增大到Φm=0.1092時(shí)，粒子運(yùn)動(dòng)形態(tài)差別不大但氣相渦量等值線已經(jīng)不再呈橢圓形或圓弧形而開始出現(xiàn)抖動(dòng)。質(zhì)量載荷繼續(xù)增大到Φm=0.5507時(shí)，氣固兩相強(qiáng)烈耦合致使大量示蹤粒子的運(yùn)動(dòng)尾跡嚴(yán)重干擾氣相流場(chǎng)。此時(shí)的氣相渦量等值線已經(jīng)完全呈現(xiàn)鋸齒狀。

在做高速復(fù)雜流場(chǎng)PIV測(cè)量時(shí)，示蹤粒子的布撒尤為重要。質(zhì)量載荷Φm的增大會(huì)對(duì)氣相流場(chǎng)的影響顯著，所以探究如何選取合適的布撒濃度對(duì)PIV測(cè)量具有指導(dǎo)意義。當(dāng)粒子St=0.01時(shí)，分別選取質(zhì)量載荷Φm=0.0121、0.1092和0.5507作為混合層初始的布撒濃度進(jìn)行氣固雙向耦合。不同Φm示蹤粒子對(duì)超聲速混合層的湍流變動(dòng)如圖7～10所示。

當(dāng)St=0.01時(shí)，大質(zhì)量載荷Φm=0.5507的混合層的時(shí)均流向速度和溫度曲線嚴(yán)重偏離無負(fù)載粒子時(shí)，如圖7所示。小質(zhì)量載荷Φm=0.0121和中質(zhì)量載荷Φm=0.1092混合層的時(shí)均速度和溫度曲線與無負(fù)載時(shí)差別不大，但有稍許的遲滯現(xiàn)象。當(dāng)St數(shù)一定時(shí)，大質(zhì)量載荷的氣固混合層中粒子基數(shù)很大，氣固兩相之間的動(dòng)量和熱量交換強(qiáng)烈，這就造成高速來流和低速來流之間的分層明顯，混合層內(nèi)的法向速度梯度和溫度梯度大。

Fig.7Normalizedmeanstreamwisevelocityandtemperaturedistribution

從圖8和9無量綱流向和法向速度脈動(dòng)曲線可以看到，Φm=0.5507時(shí)大量粒子所產(chǎn)生的尾跡抹平了部分流向和法向速度脈動(dòng)。此時(shí)的最大流向速度脈動(dòng)比無負(fù)載時(shí)減小了79.85%，最大法向速度脈動(dòng)只有無負(fù)載粒子時(shí)的38.63%。當(dāng)Φm=0.1092時(shí)，在-0.125

圖10是不同Φm的示蹤粒子對(duì)超聲速混合層雷諾應(yīng)力的影響。Φm=0.0121和0.1092時(shí)粒子布撒到混合層中，增強(qiáng)了混合層中心線附近的雷諾應(yīng)力。而大質(zhì)量載荷Φm=0.5507時(shí)，添加示蹤粒子后的整個(gè)混合層區(qū)域的雷諾應(yīng)力值小于無負(fù)載時(shí)。小質(zhì)量載荷時(shí)在渦核及其周圍粒子分布較為均勻，湍流切應(yīng)力的增大導(dǎo)致雷諾應(yīng)力的增大。中等Φm=0.1092時(shí)粒子在低渦量區(qū)域的渦邊緣和渦辮區(qū)域濃度增大。與無負(fù)載時(shí)相比較，-0.125

4 結(jié) 論

(1)Φm=0.0121時(shí)，St=100粒子的粒徑大，熱量交換緩慢，粒子不能很好地跟隨無負(fù)載粒子時(shí)的時(shí)均溫度。St=100和1的粒子削弱了湍流脈動(dòng)，St=100時(shí)雷諾應(yīng)力峰值只有無負(fù)載時(shí)的41.74%。St=0.01的粒子減弱了流向湍流而強(qiáng)化了法向湍流，使混合層中心線處的雷諾應(yīng)力值增大了9.68%。但St=0.01的粒子能較好地跟隨超聲速氣固混合層時(shí)均速度和溫度，摻混能力強(qiáng)且對(duì)原始流場(chǎng)影響較小，適合作為PIV示蹤粒子；

(2)St=0.01時(shí)，Φm=0.5507時(shí)氣固混合層中的粒子基數(shù)很大，大量粒子的運(yùn)動(dòng)尾跡抹平了部分法向速度脈動(dòng)，使流向和法向速度脈動(dòng)的峰值分別減小了約28.83%和61.37%。Φm=0.0121時(shí)的混合層渦核的附近聚集了少量的示蹤粒子，增強(qiáng)了中心線及其附近的法向速度脈動(dòng)，最大流向速度脈動(dòng)被削弱了19.29%。而Φm=0.1092時(shí)的法向速度脈動(dòng)與無負(fù)載粒子時(shí)差別不大，只在低速流動(dòng)區(qū)域略微減小，適合作為PIV初始布撒濃度。粒子Φm對(duì)雷諾應(yīng)力的影響與速度脈動(dòng)相似；

(3) 對(duì)比超聲速混合層中不同St數(shù)和Φm的湍流變動(dòng)作用，PIV示蹤粒子St=0.01和Φm=0.1092時(shí)，粒子給原始無負(fù)載流場(chǎng)帶來的湍流變動(dòng)作用最小，在高速流動(dòng)PIV診斷測(cè)量中可以作為示蹤粒子粒徑和布撒濃度的選擇依據(jù)。

[1]Elsinga G E, Van Oudheusden B W, Scarano F. Evaluation of aero-optical distortion effects in PIV[J]. Experiments in Fluids, 2005, 39(2): 246-256.

[2]Gore R A, Crowe C T. Effect of particle size on modulating turbulent intensity[J]. Multiphase Flow, 1989, 15: 297-285.

[3]Kenning V M, Crowe C T. On the effect of particles on carrier phase turbulence in gas-particle flows[J]. International Journal of Multiphase Flow, 1997, 23(2): 403-408.

[4]Crowe C T. On models for turbulence modulation in fluid-particle flows[J]. International Journal of Multi-phase Flow, 2000, 26(5): 719-727.

[5]Hetsroni G. Particle-turbulence interaction[J]. Multiphase Flow, 1989, 15: 735-746.

[6]Truesdell G C, Elghobashi S. On the two-way interaction between homogeneous turbulence and dispersed solid particles(II): Particle dispersion[J]. Physics of Fluids, 1994, 6(3): 1405-1407.

[7]Elghobashi S, Truesdell G C. On the two-way interaction between homogeneous turbulence and dispersed solid particles I: Turbulence modification[J]. Physics of Fluids A: Fluid Dynamics (1989-1993), 1993, 5(7): 1790-1801.

[8]Lain S, Broder D, Sommerfeld M, et al. Modelling hydrodynamics and turbulence in a bubble column using the Euler-Lagrange procedure[J]. International Journal of Multiphase Flow, 2002, 28(8): 1381-1407.

[9]Lain S, Sommerfeld M. Turbulence modulation in dispersed two-phase flow laden with solids from a La-grangian perspective[J]. International Journal of Heat and Fluid Flow, 2003, 24(4): 616-625.

[10]Savolainen K, Karvinen R. The effect of particles on gas turbulence in a vertical upward pipe flow[C]. Third International Conference on Multiphase Flow, ICMF, 1998.

[11]Druzhinin O A. The influence of particle inertia on the two-way coupling and modification of isotropic turbu-lence by microparticles[J]. Physics of Fluids, 2001, 13(12): 3738-3755.

[12]Ling W, Chung J N, Troutt T R, et al. Direct numerical simulation of a three-dimensional temporal mixing layer with particle dispersion[J]. Journal of Fluid Mechanics, 1998, 358: 61-85.

[13]林建忠, 石興, 余釗圣. 二維氣固兩相混合層中固粒對(duì)流場(chǎng)影響的研究[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2000, 21(8): 771-776. Lin J Z, Shi X, Yu Z S. Research on the effect of partical of two-dimensional shear flow[J]. Applied Mathematics and Mechanics, 2000, 21(8): 771-776.

[14]林建忠, 石興. 氣固兩相混合層流場(chǎng)雙向耦合的數(shù)值研究[J]. 工程熱物理學(xué)報(bào), 2001, 22(4): 496-499. Lin J Z, Shi X. Numerical study on the gas-solid two-phase mixing layer by two-way coupling model[J]. Journal of Engineering Thermophysics, 2001, 22(4): 496-499.

[15]樊建人, 羅坤, 金晗輝, 等. 直接數(shù)值模擬三維氣固兩相混合層中顆粒與流體的雙向耦合[J]. 中國電機(jī)工程學(xué)報(bào), 2003, 23(4): 153-157. Fan J R, Luo K, Jin H H, et al. Direct numerical simulation of the two-way coupling effects between particles and fluid in the three-dimensional partical-laden mixing layer[J]. Proceedings of the CSEE, 2003, 23(4): 153-157.

[16]王兵, 劉毅, 王希麟. 顆粒粒徑和氣載比改變時(shí)湍流調(diào)制的變化[J]. 航空動(dòng)力學(xué)報(bào), 2009, 24(8): 1818-1823. Wang B, Liu Y, Wang X L. Turbulence modulation caused by change of particle size and particle mass loading ratio[J]. Journal of Aerospace Power, 2009, 24(8): 1818-1823.

[17]劉偉, 萬國新, 陳景兵, 等. 可壓縮氣固混合層中離散相與連續(xù)相的相互作用研究[J]. 計(jì)算力學(xué)學(xué)報(bào), 2009, 26(1): 8-14. Liu W, Wan G X, Chen J B, et al. Study on the interaction between the continuous and the dispersed phase in compressible gas-solid mixing layer[J]. Chinese Journal of Computation Mechanics, 2009, 26(1): 8-14.

[18]劉偉, 陳景兵, 萬國新, 等. 空間模式下可壓縮氣固兩相氣固混合層流動(dòng)特性研究[J]. 應(yīng)用力學(xué)學(xué)報(bào), 2010, 27(2): 8-14. Liu W, Chen J B, Wan G X, et al. Characteristics of the flow of spatial developing gas-solid mixing layer[J]. Chinese Journal of Applied Mechanics, 2010, 27(2): 8-14.

[19]Yoshizawa A, Horiuti K. A statistically-derived subgrid-scale kinetic energy model for the large-eddy simulation of turbulent flows[J]. Journal of the Physical Society of Japan, 1985, 54(8): 2834-2839.

[20]Goebel S G, Dutton J C, Krier H, et al. Mean and turbulent velocity measurements of supersonic mixing layers[J]. Experiments in Fluids, 1990, 8: 263-272.

[21]任兆新, 王兵, 張會(huì)強(qiáng), 等. 超音速擬序旋渦中顆粒的彌散[J]. 工程熱物理學(xué)報(bào), 2015, 36(9): 1957-1960. Ren Z X, Wang B, Zhang H Q, et al. Dispersion of particles in the supersonic coherent vortex[J]. Journal of Engineering Thermophysics, 2015, 36(9): 1957-1960.

[22]Rong Z, Liu H, Chen F. Development and application of PIV in supersonic flows[C]. AIAA Aerospace Sciences Meeting Including the New Horizons Forum & Aerospace Exposition. American Institute of Physics, 2011.

[23]張亞, 陳方, 劉洪, 等. 高速流動(dòng)中PIV示蹤粒子松弛特性研究[J]. 實(shí)驗(yàn)流體力學(xué), 2013, 27(6): 70-75. Zhang Y, Chen F, Liu H, et al. Research on the relaxation characteristics of PIV tracer particles in supersonic flow[J]. Journal of Experiments in Fluid Mechanics, 2013, 27(6): 70-75.

[24]劉洪, 陳方, 勵(lì)孝杰, 等. 高速復(fù)雜流動(dòng)PIV技術(shù)實(shí)踐與挑戰(zhàn)[J]. 實(shí)驗(yàn)流體力學(xué), 2016, 30(1): 28-42. Liu H, Chen F, Li X J, et al. Practices and challenges on PIV technology in high speed complex flows[J]. Journal of Experiments in Fluid Mechanics, 2016, 30(1): 28-42.

InvestigationofturbulencemodificationbyPIVtracerparticlesinasupersonicmixinglayer

Chen Xiaohu1, Chen Fang1,*, Liu Hong1, Sha Sha2, Lu Xueling2, Zhang Qingbing2

(1. School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China; 2. Beijing Institute of Electronic System Engineering, Beijing 100854, China)

The turbulence modification by Particle Image Velocimetry (PIV) tracer particles was investigated with the two-way coupling of 2D spatial development supersonic gas-solid two-phase mixing layer. The gas phase coupled with the dispersion phase were simulated by the large eddy simulation and the Lagrangian trajectory model, respectively. It is found that the streamwise turbulence of the mixing layer is weakened, the transverse turbulence of the mixing layer strengthened and the Reynolds stress peak value increased by 9.68% than that of the unladen mixing layer due to the small Stokes particles. However, the root mean square velocity of the mixing layer is weakened by the large Stokes particles, and the Reynolds stress peak value is only 41.74% of that of the unladen mixing layer. The root mean square transverse velocity with the large mass loading is 38.63% of that of unladen mixing layer, and it is partly counteracted by the motion of particles. The root mean square transverse velocity and the Reynolds stress with the middle mass loading are almost equal to those of the unladen mixing layer. The largest root mean square streamwise velocity is reduced by 19.29% whereas the root mean square transverse velocity near the centerline of the mixing layer is increased by the small mass loading. The modification to the turbulence of the supersonic mixing layer under the small Stokes number and middle mass loading condition is smaller than other cases. The study provides important reference for improving PIV experiment in high speed flows.

PIV; supersonic mixing layer; turbulence modification; Stokes number; mass loading

2016-09-20;

2016-10-30

國家自然科學(xué)基金(11672183)

*通信作者 E-mail: fangchen@sjtu.edu.cn

ChenXH,ChenF,LiuH,etal.InvestigationofturbulencemodificationbyPIVtracerparticlesinasupersonicmixinglayer.JournalofExperimentsinFluidMechanics, 2017, 31(6): 8-14,21. 陳小虎, 陳 方, 劉 洪, 等. 超聲速混合層中PIV粒子的湍流變動(dòng)作用研究. 實(shí)驗(yàn)流體力學(xué), 2017, 31(6): 8-14,21.

1672-9897(2017)06-0008-08

10.11729/syltlx20160144

V211.1+7, O359+.2

A

陳小虎(1990-)，男，河南周口人，碩士研究生。研究方向：空氣動(dòng)力學(xué)和兩相流。通信地址：上海市閔行東川路800號(hào)，上海交通大學(xué)航空航天學(xué)院(200240)。E-mail: chenxiaohu@sjtu.edu.cn

(編輯：李金勇)