堵頭管中軸對稱紊動(dòng)射流的實(shí)驗(yàn)研究

堵頭管中軸對稱紊動(dòng)射流的實(shí)驗(yàn)研究Experimentalstudyofanaxisymmetricturbulentjetdischargedintoadead-endtunnel1前言紊動(dòng)射流是流體運(yùn)動(dòng)的一種重要類型，在航天航空、動(dòng)力機(jī)械、化工設(shè)備、自動(dòng)控制、水力采掘、民用消防、農(nóng)業(yè)噴灌、環(huán)境工程、通風(fēng)空調(diào)等許多工程技術(shù)部門中都有大量的射流問題。按照射流周圍環(huán)境條件，射流可分為受納空間無限大或遠(yuǎn)比射流特征尺寸大的自由射流與受納空間有限的非自由射流或稱受限射流。在通風(fēng)空調(diào)工程中，廣泛采用射流送風(fēng)。進(jìn)入室內(nèi)的射流大部分屬于受限射流，如圖1所示。射流由噴口噴入房間，只要噴口尺寸斷面尺寸相比不是很小，其擴(kuò)展就會(huì)受到房間圍護(hù)結(jié)構(gòu)（頂棚、墻壁或地板等）的限制和影響，形成回流?；亓鞯拇嬖谑鞘芟奚淞髋c自由射流的主要區(qū)別，對自由射流的實(shí)驗(yàn)和理論研究均已比較成熟，國內(nèi)外系統(tǒng)介紹這方面研究的文獻(xiàn)很多[1～6]，但對實(shí)際工程中更常見的受限射流流動(dòng)特性研究卻不多見，所以，從工程實(shí)際出發(fā)，研究受限射流具有重要意義。圖1房間通風(fēng)射流示意圖從理論上說，受限射流是紊動(dòng)射流的一種，而紊動(dòng)射流作為典型的剪切紊流，盡管對它的研究已有幾十年的歷史，但至今仍然是一個(gè)非?；钴S的領(lǐng)域，因?yàn)檫@些研究有助于加深人們對率流特性的認(rèn)識(shí)，如對自由射流中擬序結(jié)構(gòu)的研究就使人們對紊流結(jié)構(gòu)有了新的認(rèn)識(shí)和理解[7，8]，所以對受限射流的研究也可以增進(jìn)人們對紊流現(xiàn)象的認(rèn)識(shí)。本課題采用熱膜流速儀研究堵頭圓管中不同受限度等溫軸對稱射流流場的時(shí)均和脈動(dòng)特性，探求分析計(jì)算方法進(jìn)行校核、補(bǔ)充或修改，并對射流內(nèi)部紊動(dòng)特性和運(yùn)動(dòng)機(jī)理進(jìn)行探討。2實(shí)驗(yàn)裝置2．1實(shí)驗(yàn)設(shè)備實(shí)驗(yàn)設(shè)備如圖2所示，由風(fēng)機(jī)、送風(fēng)管道、穩(wěn)流箱和實(shí)驗(yàn)段--堵頭圓管等幾部分組成。選用圓形截面堵頭管來研究受限射流，是考慮到與矩形斷面模型相比，圓管的軸對稱性可大大減少測點(diǎn)數(shù)目，而且文獻(xiàn)[9]中比較了矩形斷面堵頭管和圓形斷面堵頭管中受限射流的實(shí)驗(yàn)成果，指出以斷面面積為特征尺度選擇無量綱數(shù)，模型的斷面形狀對流速、流量的大小影響不大，只是對射流和回流的位置有影響。圖2實(shí)驗(yàn)裝置示意圖1.吸風(fēng)口2.送風(fēng)機(jī)3.軟接頭4.鐵皮送風(fēng)管道5.風(fēng)量調(diào)節(jié)閘板6.穩(wěn)流方箱7.格柵與紗網(wǎng)8.連接段9.射流噴嘴10.堵頭圓管2．2量測斷面及測點(diǎn)布置模型實(shí)驗(yàn)段如圖3所示，堵頭圓管內(nèi)徑D=240mm，全長L=170mm，等于模型直徑的7倍多。根據(jù)羅津別爾格[9]的實(shí)驗(yàn)，射流作用的最大長度不超過堵頭管直徑的4倍數(shù)，故該模型長度充分有余。圓管制一端被封死，開口的一端裝有送風(fēng)噴嘴和連接段，送風(fēng)噴嘴用螺栓固定在噴嘴支板上。圖3模型實(shí)驗(yàn)段1.連接段2.噴嘴支板3.噴嘴4.堵頭圓管由于模型的軸對稱性，每個(gè)斷面只須在一條直徑上量測，在堵頭管的上部和下部各開有一排圓孔，上部圓孔直徑為8.5mm，其中標(biāo)有"c"的圓孔在位置為最測斷面，直探針（熱膜探針和邊界探針）從這些圓孔引入模型內(nèi)部矛盾測量；標(biāo)有"f"的圓孔是輔助量測孔，供彎頭長50mm的皮托管引入模型對量測斷面進(jìn)行測量表中4．1．2各組實(shí)驗(yàn)量測結(jié)果實(shí)驗(yàn)量測的時(shí)均流速及脈動(dòng)流速相對值的典型結(jié)果見圖6。圖中的軸向流速脈動(dòng)相對值是指，um為射流軸線時(shí)均流速。射流邊界位置由正反流邊界探針確定。邊界附近的時(shí)均流速參考皮托管量測進(jìn)行了修正（圖6a給出修正后的時(shí)均流速分布）。距噴嘴出口距離增大后，由于流速減小，正反流邊界探針分辨不出邊界，所以后面的斷面沒有給出時(shí)均流速分布，但給出了脈動(dòng)流速分布。圖6時(shí)均流速與軸向流速脈動(dòng)相對值沿程分布圖4．2堵頭管中受限射流的一般特性及分區(qū)根據(jù)對實(shí)驗(yàn)成果的整理分析，受限射流流場可分為四個(gè)區(qū)：I自由擴(kuò)展區(qū)，Ⅱ受限擴(kuò)展區(qū)，Ⅲ收縮區(qū)，Ⅳ零均流區(qū)，如圖7所示。圖7堵頭管中受限射流流場示意圖I自由擴(kuò)展區(qū)，Ⅱ受限擴(kuò)展區(qū)，Ⅲ收縮區(qū)，Ⅳ零均流區(qū)在第一臨界斷面C1之前為自由擴(kuò)展區(qū)，C1斷面之后，射流的擴(kuò)展受到邊界和回流的限制，其卷吸周圍流體作用減弱，射流半徑流量的增加速率逐漸減慢，但總的趨勢仍是半徑增大，流量增加，故稱為受限擴(kuò)展區(qū)，直到第二臨界斷面C2。此后，射流運(yùn)動(dòng)發(fā)生根本轉(zhuǎn)折，射流流線開始越出邊界產(chǎn)生回流，射流區(qū)流量沿程減少，在第二臨界斷面上射流區(qū)流量達(dá)到最大值，而射流半徑在C2斷面稍后達(dá)到最大值，然后逐漸縮小，稱為收縮區(qū)，收縮區(qū)末端位置（第三臨界斷面）由射流最遠(yuǎn)作用長度決定?？筛鶕?jù)實(shí)際工程中對軸線流速的限值給出（詳見4.3.1）。第三臨界斷面和圓管堵頭之間是第Ⅳ區(qū)--零均流區(qū)，或稱尾部旋渦區(qū)，該區(qū)氣流的特點(diǎn)是處理"阻滯"和"堆積"狀態(tài)，形成一個(gè)或多個(gè)不穩(wěn)定的大尺度旋渦，旋渦誘導(dǎo)的瞬時(shí)流甚，時(shí)間平均流速為零，壓力均勻。堵頭管中的射流，在第I和第Ⅱ區(qū)所卷吸的流體來自本身的回流，管內(nèi)空氣產(chǎn)生多次循環(huán)，排除的只是從噴口噴入的那部分空氣量。4．3受限射流的時(shí)均流動(dòng)特性4．3．1受限射流軸線流速變化規(guī)律射流進(jìn)入堵頭管后，由于動(dòng)量的衰減和橫向傳遞，流速不斷降低。圖8是射流軸線上時(shí)均流速隨無量綱距離x/d0的沿程變化。從圖中可見，軸線流速的變化取決于射流的受限度。當(dāng)受限度很大時(shí)，軸線流速急劇下降，射流相對射程x/d0隨著受限度減小而增大。而且隨著受限程度減小，軸線流速的變化曲線互相趨于接近，當(dāng)受限度很時(shí)，趨近于圖8中用虛線表示的自由射流軸線流速變化曲線，射流受了程度越大，其軸線流速變化曲線um/u0=f(x/d0)離自由射流的曲線越遠(yuǎn)。圖8軸線流速隨x/d0的變化圖9是用無量綱距離表示的軸線流速沿程變化。圖9a中曲線表示巴哈列夫經(jīng)驗(yàn)公式計(jì)算值。從圖中可見，實(shí)測值與經(jīng)驗(yàn)公式計(jì)算值在趨勢上是一致的，但在數(shù)值上，前者大于后者。巴哈列夫的受限射流軸線流速經(jīng)驗(yàn)公式是把受限度引入自由射流軸線流速公式之后擬合實(shí)驗(yàn)數(shù)據(jù)得到的。他采用的自由射流軸線流速公式為（1）式中<1時(shí)，≥1時(shí)，A=6.5a（2）圖9軸線流速隨的變化巴哈列夫把自由射流看作受限射流的極限情況，將受限度以e的指數(shù)形式引入式（1），得以受限射流軸線流速表示式：（3）式中的指數(shù)m即反映射流受限程度的影響，根據(jù)實(shí)驗(yàn)數(shù)據(jù)，巴氏得到（4）由此可見1自由射流的軸線流速公式在巴氏受限射流軸線流速經(jīng)驗(yàn)公式中起著重要作用。圖10比較了式（1）和Tollmien的理論解[1]，即，發(fā)現(xiàn)巴氏所采用的自由射流軸線流速公式（1）的計(jì)算值低于其它學(xué)者的研究成果。本文采用Tollmien理論解的數(shù)據(jù)擬合式（1）中的因數(shù)A，得到當(dāng)<1時(shí)（5a）A=6.85a當(dāng)≥1時(shí)（5b）于是經(jīng)過修正的受限射流軸線流速經(jīng)驗(yàn)公式變?yōu)椋?）式中A按式（5）選取。圖10自由射流線流速沿程變化圖8、圖9b中的曲線均是根據(jù)式（6）計(jì)算得到的。從圖中可見，實(shí)測值與之基本相符?？紤]到射流作用的末端，斷面上流速分布趨于均勻，軸線流速與斷面平均流速相差不大，在實(shí)際工程中可把軸線流速達(dá)到給定值ud時(shí)的長度定義為射流最遠(yuǎn)作用長度。ud根據(jù)需要選定，在通風(fēng)空調(diào)中ud一般為0.07～0.04m/s[9][13]，代入式（6）可以得到受限射流的最遠(yuǎn)作用長度，也就是第三臨界斷面的位置4．3．2受限射流流量的變化規(guī)律根據(jù)實(shí)測斷面流速分布，可得受限射流的流量Qj。圖10是受限射流流量隨時(shí)無量綱距離的沿程變化。由于管壁邊界限制了射流的發(fā)展，受限射流的流量不像自由射流那樣沿程增加，而是先逐漸增大，到達(dá)第二臨界斷面之后，由于氣流的回轉(zhuǎn)，射流流量逐漸減小。在圖10中，從射流噴口到第一臨界斷面的自由擴(kuò)展區(qū)，流量增加很快；隨后的受限擴(kuò)展區(qū)，流量增加速度逐漸減慢并在第二臨界斷面處達(dá)到最大值；第二臨界斷面以后，射流流量逐漸減小。射流流量的沿程變化在一定程度上反映了射流的卷吸性質(zhì)，從圖11a可見，射流受限越大，相對流量Qj/Q0最大值越小，卷吸量越小。圖11b是各個(gè)受限度射流無量流量的綜合結(jié)果。從圖中可見，不同受限度的實(shí)驗(yàn)點(diǎn)比較接近，可用一條曲線表示。圖中曲線1是根據(jù)巴哈列夫經(jīng)驗(yàn)公式計(jì)算得到的。從圖中可見實(shí)驗(yàn)點(diǎn)與之偏差較大，歸因于巴哈列夫量測流速的儀器--電風(fēng)速儀在紊動(dòng)強(qiáng)度較高時(shí)精度較差。仍選用巴氏經(jīng)驗(yàn)公式的形式，擬合實(shí)驗(yàn)數(shù)據(jù)，得到（7）如圖11b中曲線2所示。圖11a中的曲線也是根據(jù)式（7）給出的。兩圖中小受限度的一些斷面與其它實(shí)驗(yàn)點(diǎn)偏離的原因可能是小受限度時(shí)，噴口直徑很小，相應(yīng)的噴口流量Q0也小，這樣流速量測時(shí)較小的誤差可能導(dǎo)致Qj/Q0較大的誤差，故在擬合式（7）時(shí)沒有包括明顯有誤差之點(diǎn)。圖11射流流量隨的沿程變化根據(jù)實(shí)驗(yàn)結(jié)果，射流在各種受限度下，第二臨界斷面均在同一無量綱距離c2=0.21處,在此斷面上，射流量達(dá)到最大值：（8）4．3．3受限射流射流區(qū)平均流速變化規(guī)律根據(jù)受限射流流量及射流半徑的變化規(guī)律，可得到受限射流射流區(qū)平均流速的沿程變化規(guī)律，如圖12所示。從圖12b可見，各受限度下射流平均流速無量綱綜合值的分布很有規(guī)律。由上述實(shí)驗(yàn)成果，可得到射流區(qū)平均流速變化的經(jīng)驗(yàn)公式（9）根據(jù)式（9）所作的曲線與實(shí)驗(yàn)數(shù)據(jù)基本相符。圖12射流區(qū)平均流速沿程變化4．3．4受限射流回流區(qū)平均流速變化規(guī)律如前所述，堵頭管中受限射流，由于邊界的限制，射流的外部存在回流。在通風(fēng)空調(diào)工程中，工作區(qū)通常位于回流區(qū)，這是因?yàn)榛亓鲄^(qū)流動(dòng)的流速和溫差已充分衰減，具有比較均勻穩(wěn)定流速場和溫度場，所以回流區(qū)平均流速是通風(fēng)空調(diào)設(shè)計(jì)的重要指標(biāo)，在實(shí)際工程中，往往要限制回流區(qū)平均流速最大值不得超過工作區(qū)風(fēng)速的允許值。目前，在通風(fēng)空調(diào)工程設(shè)計(jì)中，巴哈列夫?qū)嶒?yàn)得到的回流區(qū)平均流速經(jīng)驗(yàn)公式（10）被廣泛采用[13，14]，但在文獻(xiàn)[9]中這個(gè)經(jīng)驗(yàn)公式里的回流區(qū)平均流速定義不明確。斷面平均流速是相應(yīng)流量和相應(yīng)面積的比值，但文獻(xiàn)[9]卻先給出回流區(qū)平均流速在堵頭管長度上的變化。擬合得到式（10），然后由連續(xù)原理得到回流流量Qr（等于同一斷面的射流流量Qj）和式（10），給出了回流面積，未免令人費(fèi)解。本文定義每個(gè)斷面上的回流區(qū)面積Fr等于圓管面積Fn與射流區(qū)面積Fj之差，即Fr=Fn-Fj（11）再由（12）得到各斷面的回流區(qū)平均流速。圖13表示回流區(qū)平均流速隨的沿程變化。從圖中可見，噴口直徑越大，越小，回流區(qū)平均流速無量綱vr/u0越大。并且受限度不同的射流，在第二臨界為面c2=0.21處回流區(qū)平均流速均達(dá)到最大值。圖13回流區(qū)平均流速隨的沿程變化圖13b中曲線是不同受限度實(shí)驗(yàn)數(shù)據(jù)的擬合曲線，所得的擬合公式為（13）回流區(qū)平均流速從噴口開始逐漸增加，至第二臨界斷面c2=0.21達(dá)到最大值，然后又逐漸減小。在第二臨界斷面上，回流區(qū)平均流速最大值等于（14）4．4受限射流的紊動(dòng)特性根據(jù)圖6中給出的堵頭管中受限射流軸向流速及相對值的沿程分布，可以對受限射流的紊動(dòng)特性進(jìn)行初步探討。圖14是不同受限度下射流軸線上隨x/d0的變化（在射流軸線上，ε1與相對紊動(dòng)強(qiáng)度相同）。圖中虛線取自參考文獻(xiàn)[15][16]中自由射流數(shù)據(jù)。圖14射流軸線上相對紊動(dòng)強(qiáng)度沿程變化從圖中可知，在自由擴(kuò)展區(qū)，不同受限度射線軸線上相對紊動(dòng)強(qiáng)度變化與自由射流基本一致；而自由擴(kuò)展區(qū)之后，受限射流細(xì)線上紊動(dòng)強(qiáng)度沿程持續(xù)增加，并不是趨于常數(shù)。分析其原因，無論是自由射流還是堵頭管中的受限射流，在射流軸線上ui/xj=0，因而軸線上不產(chǎn)生紊動(dòng)，軸線上的紊動(dòng)是由上游的輸運(yùn)和附近的擴(kuò)散引起的。在受限條件下，由于回流的存在，射流斷面上流速分布更不均勻，軸線附近紊動(dòng)加大，因而擴(kuò)散到軸線上的紊動(dòng)也相應(yīng)地增加；同時(shí)，主流的波動(dòng)也是導(dǎo)致相對紊動(dòng)強(qiáng)度增加的原因。從圖6可見，不同受限度射流軸向流速脈動(dòng)ε1在斷面上的徑向分布情況，氣流從噴口噴入堵頭管中，在勢流核心區(qū)沒有完全消失的斷面，核心區(qū)內(nèi)ε1很?。浑S著徑向坐標(biāo)r的增加，ε1先是逐漸增大，增大到一定程度后開始減小，直至射流邊界附近；此后r再增大，ε1基本保持為常數(shù)。結(jié)合斷面上時(shí)均流速的分布，可以得到，從射流軸線開始，隨著r的增加，相對紊動(dòng)強(qiáng)度先是和ε1一樣逐漸增大，但當(dāng)ε1達(dá)到最大值并開始減小時(shí)，由于時(shí)均流速u的減小，ε將趨近于∞。而在射流邊界以外的回流區(qū)，雖然ε1保持在一個(gè)較低水平，但由于時(shí)均流速同樣較低，相對紊動(dòng)強(qiáng)度ε值仍然較大。從圖6中還可以看出，隨著軸向距離的增加，整個(gè)斷面的軸向流速脈動(dòng)均增大，并逐漸趨于均勻。圖15給出的是d0=20mm，=10.6的射流各斷面上回流區(qū)軸向流速脈動(dòng)相對于回流區(qū)平均速度vr的沿程變化，從圖中可以看出，隨著回流越接近出口，ε2值越大，這主要是vr迅速減小的結(jié)果。圖15回流區(qū)ε2沿程變化5結(jié)語5．1以往對堵頭管中射流時(shí)均特性的研究成果定性上是正確的。但是，由于實(shí)驗(yàn)多在早期完成，受到量測手段的局限，定量存在一定的誤差。本文采用了熱膜流速儀以及專門設(shè)計(jì)的正反流邊界探針，使得研究成果具有較高的可靠性。5．2堵頭管中射流的發(fā)展可分為四個(gè)區(qū)：自由擴(kuò)散區(qū)、受限擴(kuò)展區(qū)、收縮區(qū)和零均流區(qū)。分區(qū)的第一臨界斷面和第二臨界斷面分別在c1=0.09和c2=0.21處，第三臨界斷面則由射流最大作用長度決定。5．3對巴哈列夫關(guān)于受限射流軸線流速沿程變化的經(jīng)驗(yàn)公式進(jìn)行了修正，見式（6）受限射流、流量射流區(qū)平均速度沿程變化可分別按經(jīng)驗(yàn)公式（7）和（9）計(jì)算。5．4通過定義回流區(qū)面積等于堵頭管面積與射流區(qū)面積之差，給出了通風(fēng)空調(diào)設(shè)計(jì)的重要指標(biāo)--回流區(qū)平均速度的經(jīng)驗(yàn)公式（13），澄清了巴哈列夫?qū)亓鲄^(qū)平均速度定義的不明確之處。5．5在實(shí)驗(yàn)所包括的協(xié)諾數(shù)范圍內(nèi)(Re0>104),受限射流的流動(dòng)與雷諾數(shù)無關(guān)，這主要是由于在高雷諾數(shù)情況下，旋渦結(jié)構(gòu)的變化主要在于最小旋渦尺度的減小，而大旋渦尺度基本不變，由于卷吸強(qiáng)度受尺度旋渦的影響，故受限射流的卷吸強(qiáng)度亦不隨Re0而變化。但是，受限度影響流動(dòng)的發(fā)生發(fā)展以及卷吸強(qiáng)度的大小，受限度越大，射流所卷吸的流量越小。5．6在受限條件下，由于回流的存在，射流的紊動(dòng)在整個(gè)堵頭管中都很劇烈，尤其是在收縮區(qū)和零均流區(qū)。預(yù)示依據(jù)研究堵頭管中的射流，無論是理論分析還是數(shù)值計(jì)算，必須對運(yùn)動(dòng)方程中流速脈動(dòng)項(xiàng)給予充分的重視。6參考文獻(xiàn)1AbramovichGN.TheTheoryofTurbulentJets.Massachusetts:MITPress,1963.2RajaratnamN.TurbulentJets.Elsevier,Amsterdam.1976.3RodiW.TurbulentBuoyantJetsandPlumes.Oxford:PergramonPress,1982.4HinzeJO.Turbulence.2ndedition,NewYork:McGraw-Hill,1975.5余常昭，環(huán)境流體力學(xué)導(dǎo)論，北京：清華大學(xué)出版社，1992。6余常昭，紊動(dòng)射流，北京：高等教育出版社，1993。7KovasznayLGG.StructureandMechanis