流速梯度對懸浮顆粒脈動強度的影響

1、流速梯度對懸浮顆粒脈動強度的影響         06-02-21 17:15:00     作者：李丹勛，王興奎，王殿    編輯：studa9ngns摘要：從簡化的顆粒運動方程出發(fā), 分析了在剪切流場中顆粒脈動強度和流體脈動強度之間的關(guān)系。結(jié)果表明, 由于縱向時均流速的垂線分布梯度的作用, 顆粒在兩個方向上的脈動強度均可能超過相應(yīng)的流體脈動強度。 關(guān)鍵詞：流速梯度 Stokes數(shù) 脈動強度  1 引言對于細小懸浮顆粒在恒定

2、、均勻、各向同性紊動流場中的運動, 理論分析和試驗量測都表明: 顆粒的脈動強度總是小于相應(yīng)流體的脈動強度，而且隨顆粒Stokes數(shù)（Stokes數(shù)定義為=e/,e和分別是紊動場含能漩渦的特征頻率和顆粒的反應(yīng)頻率）的增大，顆粒脈動強度將減?。划斪銐虼髸r, 顆粒的脈動強度將趨于零1,2。隨著LDV(Laser Doppler Velocimetry)和PIV(Particle Image Velocimetry)技術(shù)的發(fā)展與應(yīng)用, 對顆粒在剪切流場中的運動有了較多的實測成果,不少文獻中發(fā)現(xiàn)了顆粒的脈動強度大于相應(yīng)流體的脈動強度的現(xiàn)象36。Lijegren 通過理論分析表明, 縱向時均流速的垂線分布

3、梯度的存在,使得顆?？v向脈動強度有可能超過流體脈動強度，而顆粒垂向脈動強度則不受流速梯度的影響7。本文在Lijegren工作的基礎(chǔ)上，考慮了流速梯度引起的Saffman力的影響, 進一步分析了顆粒脈動強度和流體脈動強度之間的關(guān)系。結(jié)果表明, 由于縱向時均流速垂向分布梯度的作用, 顆粒在兩個方向上的脈動強度均可能超過相應(yīng)的流體脈動強度。2 顆粒運動方程考慮圓球狀單顆粒在二維恒定均勻剪切流場中的運動。當粒徑較小時，阻力可用Stokes公式表示。 只考慮顆粒所受的Stokes阻力和垂向的Saffman力8, 忽略其它力的作用,則顆粒的運動方程可表示為(/6)D3p(dUp)/(dt)=3D(Uf-U

4、p)(1)(/6)D3p(dUp)/(dt)=3D(Vf-Vp)+(Uf-Up)(2)式中 D表示顆粒直徑，下標f,p分別代表流體和顆粒，U, V分別表示縱向和垂向流速，為流體動力粘性系數(shù)。, 其中G為縱向時均流的垂線分布梯度:G=d/dy. 為簡化分析, 這里取G為常數(shù).將式(1)、(2)整理得(dUp/dt)+Up=Uf(3)(dUp/dt)+Up=Uf +(Uf-Up)(4)其中=18/(D2p),=6/(pD3)。按照Lijegren7的方法對式(3)作進一步分析。引入U=Up-(5)其中Uf=+uf,uf為脈動速度。由式(5)可得(6)沿顆粒的運動軌跡觀察, 有(7)把式(5)、(6

5、)、(7)代入式(3)得(dU)/(dt)+U=-GVp+uf(8)不考慮體積力, 對細小顆粒有=0。對式(8)取平均可得(9)式(9)為簡單的一階線性微分方程，其漸近平穩(wěn)解=0, 即在擺脫初始條件的影響后,有p=f。結(jié)合式(5)可知U=up, 同時在二維恒定均勻流中有=0, 則式(8)、(4)可分別寫為(dup)/(dt)+up=-Gvp+uf(10)(dup)/(dt)+up=vf+(uf-up)(11)式(10)、(11)即為用脈動流速表示的簡化的顆粒運動方程。3 顆粒脈動強度分析對普通函數(shù)f(t), 只有當收斂時，其Fourier積分才存在。而對隨機函數(shù)而言, 任一個平穩(wěn)的隨機過程,

6、雖則當t為無窮大時并不趨于零, 但其具有明確物理意義的Fourier變換仍然是存在的9。定義,對式(10)、(11)進行Fourier變換(12)(13)整理可得(14)用 ()和()來表示()和()，則由式(14)、(15)可得(14)(15)記能譜密度:S()=|F()|2 , 則由式(16)、(17)分別可得(18)(19)其中(20)(21)(22)(23)(24)(25)對能譜密度進行積分,可得脈動速度的均方值。從式(18)(25)可得(26)(27)其中(28)(29)(30)(31)(32)(33)顆粒的縱向脈動強度由3項組成：M1項是顆粒脈動對流體脈動的響應(yīng), M2和M3兩項則

7、是由于流速梯度引起的附加項, 而且流速梯度越大, 二者的作用越明顯。易知 M2>0。在剪切流場中,的符號和流速梯度G的符號相反, 由此從式(22)、(30)知, M3的符號和(2-G)的符號相同, 當(2-G)>0時, M3>0。顆粒垂向的脈動強度也由3項組成：N1項是顆粒脈動對流體脈動的響應(yīng), N2和N3是流速梯度引起的附加項, 而且流速梯度越大, 二者的作用越明顯。易知N2>0, N3<0。下面定量考察流速梯度G的影響。記=e/,=G/B2(>0), 則有(2-2-G) 2+422=2+(1+)222+(1-)22(34)則Xi, Yi(i=1,2,3)

8、分別可以寫為(35)(36)其中記/e, 則/=, 從而有(37)(38)其中:1=/(1+),2=/(1-),1。下面的分析只考慮2>0, 對于2<0的情況, 可得到相同的結(jié)論。進一步的分析依賴于對能譜密度具體形式的認識。Lijegren7在其分析中采用如下能譜密度形式(39)u為流體的特征脈動速度, 在各向同性流場中, 有u2=,為Lagrangian長度積分比尺, 對于Xii, 有(40)d為截止頻率, C0為常數(shù):C0=e/u=3/2(2-e/d)本文采用同樣的能譜密度形式。記(41)把式(39)、(40)代入式(41), 在1很小而1d/e很大時，可得f(1)的近似值(4

9、2)同理在2很小而2d/e很大時, 可得f(2)的近似值為(43)由以上的分析可得(44)同理可得(45)(46)由此可得(47)從式(47)可知, 當下式成立時(48)顆粒的縱向脈動強度將大于相應(yīng)的流體的脈動強度。同時也可以看到，在不存在流速梯度時(G=0,=0)有：/=1-/4+O(2),即顆粒的脈動強度將小于流體的脈動強度。可見流速梯度對顆粒的縱向脈動強度有著直接而重要的影響。對于顆粒的垂向脈動強度, 按同樣的分析方法,可得(49)(50)(51)從而有(52)在<1時，有1/(1-)2>1，所以在較小時，總有可能使得(52)式的值大于1，即顆粒的垂向脈動強度大于流體的垂向脈

10、動強度。由上式也可知，在不存在流速梯度時(=0,=0)有：/=1-/4+O(2),即顆粒的垂向脈動強度總是小于流體的脈動強度?？梢娏魉偬荻葘︻w粒的垂向脈動強度大小也有著直接而重要的影響。(1)縱向(2)垂向 圖1 顆粒脈動強度和流體脈動強度的關(guān)系The relationship between turbulent intensity of two phases由以上分析可知，由于流速梯度的作用，顆粒脈動強度有可能超過流體的脈動強度, 下面給出一個具體的算例。假定流場中含能渦旋的特征頻率在100HZ以下, 表征紊動渦旋尺度的Komogorov比尺約為0.1mm的量級。對于0.010.1mm的細小輕質(zhì)塑料圓球顆粒(容重為1.056), 取e=5