超聲速混合層中PIV粒子的湍流變動作用研究

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

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

超聲速混合層中PIV粒子的湍流變動作用研究

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

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

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

PIV；超聲速混合層；湍流變動；Stokes數；質量載荷

0 引 言

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

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

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

1 數值計算方法

1.1 氣相控制方程

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

分子粘性引起的張量σij和亞格子應力τijg分別表達為：

本文采用Smagorinsky-Lilly所提出的亞格子應力模型。Smagorinsky常數并不是一個固定值，Yoshizawa 和Horiuti[19]認為較大的Smagorinsky常數會導致耗散的增強，在均勻湍流中取Cs=0.2而剪切湍流中取Cs=0.065是合適的。氣相在空間上的數值離散是二階迎風格式，時間上的數值離散是二階隱式推進。

1.2 固相控制方程

單個粒子在拉格朗日框架下的運動方程為:

式中：xip、uip、mip、Tip和cip分別表示第i個粒子的位置坐標、速度、質量、溫度和比熱；Figp、Fipp和Qip分別代表第i個粒子受到的氣動力、粒子間相互作用力和傳熱量。因為本文中粒子的布撒體積分數均小于1×10-3，粒子之間的相互碰撞不明顯，所以不考慮離散相粒子之間的相互碰撞效應，故Fipp=0。

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

2 計算工況與驗證

2.1 氣相計算參數

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

表1 計算參數Table 1 Computational parameters

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

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

2.2 離散相計算參數

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

2.3 計算參數選擇

在PIV診斷測速技術中，示蹤粒子的直徑和布撒濃度是影響測試準確度的2個關鍵因素。粒徑和布撒濃度在兩相流中可以用參數St數和Φm來表征。在超聲速氣固兩相混合層中，粒徑導致的St數變化，直接改變示蹤粒子在大渦擬序結構中的分布特征從而產生湍流變動作用。而示蹤粒子Φm的增大會加劇粒子對混合層的影響。當Φm增大到一定程度之后，原始的連續相流場結構已經發生變化，此時的測量結果已經不具有可信性。

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

3 計算結果及分析

3.1 氣相計算結果驗證

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

3.2 不同St數時的湍流變動作用

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

Fig.1Normalizedstreamwisemeanandrootmeansquarevelocitydistribution

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

圖3～5分別是無量綱時均流向速度和溫度分布、無量綱流向和法向速度脈動和無量綱雷諾應力分布曲線。如圖3所示，3種粒徑的粒子在經過長時間平均后都能較好地逼近無負載時混合層的時均速度。St數對時均溫度影響較大，St=100時的粒子因為粒徑較大，兩相之間的傳熱過程與小粒子相比相對漫長。大粒徑粒子使氣相的時均溫度不能很好地跟隨無負載時的混合層。

Fig.3Normalizedmeanstreamwisevelocityandtemperaturedistribution

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

Fig.4Normalizedrootmeansquarestreamwiseandtransversevelocitydistribution

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

3.3 不同Φm時的湍流變動作用

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

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

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

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

Fig.7Normalizedmeanstreamwisevelocityandtemperaturedistribution

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

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

4 結 論

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

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

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

[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]林建忠, 石興, 余釗圣. 二維氣固兩相混合層中固粒對流場影響的研究[J]. 應用數學和力學, 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]林建忠, 石興. 氣固兩相混合層流場雙向耦合的數值研究[J]. 工程熱物理學報, 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]樊建人, 羅坤, 金晗輝, 等. 直接數值模擬三維氣固兩相混合層中顆粒與流體的雙向耦合[J]. 中國電機工程學報, 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]王兵, 劉毅, 王希麟. 顆粒粒徑和氣載比改變時湍流調制的變化[J]. 航空動力學報, 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]劉偉, 萬國新, 陳景兵, 等. 可壓縮氣固混合層中離散相與連續相的相互作用研究[J]. 計算力學學報, 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]劉偉, 陳景兵, 萬國新, 等. 空間模式下可壓縮氣固兩相氣固混合層流動特性研究[J]. 應用力學學報, 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]任兆新, 王兵, 張會強, 等. 超音速擬序旋渦中顆粒的彌散[J]. 工程熱物理學報, 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]張亞, 陳方, 劉洪, 等. 高速流動中PIV示蹤粒子松弛特性研究[J]. 實驗流體力學, 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]劉洪, 陳方, 勵孝杰, 等. 高速復雜流動PIV技術實踐與挑戰[J]. 實驗流體力學, 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

國家自然科學基金(11672183)

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

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

1672-9897(2017)06-0008-08

10.11729/syltlx20160144

V211.1+7, O359+.2

A

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

(編輯：李金勇)