導(dǎo)語：紊動射流研究論文一文來源于網(wǎng)友上傳，不代表本站觀點，若需要原創(chuàng)文章可咨詢客服老師，歡迎參考。

提要介紹了采用熱膜流速儀量測堵頭管中不同受限度軸對稱射流的時均流速和流速脈動。根據(jù)實驗，把堵頭管中射流分為自由擴展區(qū)、受限擴展區(qū)、收縮區(qū)和零均流區(qū)四個區(qū)，并給出了臨界斷面位置，對現(xiàn)有計算方法進(jìn)行了改

進(jìn)。

關(guān)鍵詞受限射流堵頭圓管熱膜測速

AbstractMeasuresthetimemeanvelocityandvelocityfluctuationusingahotfilmanemometerwithdifferentconfinementratiofortheaxisymmetricturbulentjitIadead-endtunnel.Byanalysingthedatafromthemeasurement,proposesthedivisionoffourflowregions-thefreespread,confinedspread,contractionandzero-time-averaged-flowregions,andthecriterionandflowcharacteristicsforeachofthem,withafewnewempiricalformulasderived.

Keywordsconfinedaxisymmetricjetdead-endtunnelHWAHFAvelocitymeasurement

1前言

紊動射流是流體運動的一種重要類型，在航天航空、動力機械、化工設(shè)備、自動控制、水力采掘、民用消防、農(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所示。射流由噴口噴入房間，只要噴口尺寸斷面尺寸相比不是很小，其擴展就會受到房間圍護(hù)結(jié)構(gòu)（頂棚、墻壁或地板等）的限制和影響，形成回流。回流的存在是受限射流與自由射流的主要區(qū)別，對自由射流的實驗和理論研究均已比較成熟，國內(nèi)外系統(tǒng)介紹這方面研究的文獻(xiàn)很多[1～6]，但對實際工程中更常見的受限射流流動特性研究卻不多見，所以，從工程實際出發(fā)，研究受限射流具有重要意義。

圖1房間通風(fēng)射流示意圖

從理論上說，受限射流是紊動射流的一種，而紊動射流作為典型的剪切紊流，盡管對它的研究已有幾十年的歷史，但至今仍然是一個非常活躍的領(lǐng)域，因為這些研究有助于加深人們對率流特性的認(rèn)識，如對自由射流中擬序結(jié)構(gòu)的研究就使人們對紊流結(jié)構(gòu)有了新的認(rèn)識和理解[7，8]，所以對受限射流的研究也可以增進(jìn)人們對紊流現(xiàn)象的認(rèn)識。

本課題采用熱膜流速儀研究堵頭圓管中不同受限度等溫軸對稱射流流場的時均和脈動特性，探求分析計算方法進(jìn)行校核、補充或修改，并對射流內(nèi)部紊動特性和運動機理進(jìn)行探討。

2實驗裝置

2．1實驗設(shè)備

實驗設(shè)備如圖2所示，由風(fēng)機、送風(fēng)管道、穩(wěn)流箱和實驗段--堵頭圓管等幾部分組成。選用圓形截面堵頭管來研究受限射流，是考慮到與矩形斷面模型相比，圓管的軸對稱性可大大減少測點數(shù)目，而且文獻(xiàn)[9]中比較了矩形斷面堵頭管和圓形斷面堵頭管中受限射流的實驗成果，指出以斷面面積為特征尺度選擇無量綱數(shù)，模型的斷面形狀對流速、流量的大小影響不大，只是對射流和回流的位置有影響。

圖2實驗裝置示意圖

1.吸風(fēng)口2.送風(fēng)機3.軟接頭4.鐵皮送風(fēng)管道5.風(fēng)量調(diào)節(jié)閘板

6.穩(wěn)流方箱7.格柵與紗網(wǎng)8.連接段9.射流噴嘴10.堵頭圓管

2．2量測斷面及測點布置

模型實驗段如圖3所示，堵頭圓管內(nèi)徑D=240mm，全長L=170mm，等于模型直徑的7倍多。根據(jù)羅津別爾格[9]的實驗，射流作用的最大長度不超過堵頭管直徑的4倍數(shù)，故該模型長度充分有余。圓管制一端被封死，開口的一端裝有送風(fēng)噴嘴和連接

段，送風(fēng)噴嘴用螺栓固定在噴嘴支板上。

圖3模型實驗段

1.連接段2.噴嘴支板3.噴嘴4.堵頭圓管

由于模型的軸對稱性，每個斷面只須在一條直徑上量測，在堵頭管的上部和下部各開有一排圓孔，上部圓孔直徑為8.5mm，其中標(biāo)有"c"的圓孔在位置為最測斷面，直探針（熱膜探針和邊界探針）從這些圓孔引入模型內(nèi)部矛盾測量；標(biāo)有"f"的圓孔是輔助量測孔，供彎頭長50mm的皮托管引入模型對量測斷面進(jìn)行測量。

3量測儀器

3．1熱膜流速儀

實驗采用美國TSI-1050A恒溫型熱膜/熱線流速儀量測流速，探針選用長1mm，直徑51μmr1210-20型圓柱熱膜探針。熱膜流速儀以熱平衡原理為基礎(chǔ)，利用置于流場中由電流加熱的敏感元件來量測測點速度。TSI-1050A熱膜流速儀框圖如圖4所示，液體速度的變化使探針產(chǎn)生的電訊號經(jīng)1050A放大處理后，輸入IFA-200進(jìn)行A-D轉(zhuǎn)換，然后輸入計算機主機（HP-100）。智能型終端（2623A型）由鍵盤控制，并通過它來選擇各種采樣參數(shù)，進(jìn)行數(shù)據(jù)存儲與計算分析處理，結(jié)果可用表格和圖形在屏幕上顯示，也可由熱敏打印機輸出。

圖4TSI-1050A熱膜流速儀框圖

3．2正反流分界線探針

噴入堵頭管中的受限射流，其流場特點是內(nèi)部存在零流速點，且上下流動方向相反，而采用標(biāo)準(zhǔn)的熱膜探針不論流速方向如何，輸出的電壓值總是正值，所以必須采用其它方法判定流速方向。根據(jù)文獻(xiàn)[11][12]設(shè)計了正反流分界線探針（圖5），探針頭部是兩根相對的直徑為1mm的細(xì)管，兩細(xì)管彼此錯開約1mm，可感受兩個方向的壓力。將此探針由各測斷面的圓孔伸入堵頭管內(nèi)，只有當(dāng)軸向流速為零時，兩細(xì)管感受的壓力才相等，將探針在各斷面上下移動，即可測得正反流分界位置，即射流邊界。文獻(xiàn)[12]用這種探針測量了通道內(nèi)穩(wěn)定器后帶回流區(qū)流場正反流分界線，通過分析認(rèn)量測結(jié)果是可靠的。

5正反流分界線探針

4實驗成果與分析

4．1實驗條件與量測結(jié)果

4．1．1實驗基本參數(shù)

為了比較全面地研究受限射流，實驗采用了6種不同的噴嘴出口直徑，對應(yīng)著較大的受限度范圍（=4.25～53.17，F(xiàn)n為堵頭圓管橫截面面積，do為噴嘴出口直徑）。噴嘴出口流速的選擇，保證了出口雷諾數(shù)處于自模區(qū)（Re0>104）。各組實驗的基本參數(shù)見表1。

表1受限射流實驗基本參數(shù)

實驗組次123456

噴口直徑d0/mm4510204050

堵頭圓管直徑/mm240240240240240240

受限度53.1742.521.310.65.324.25

噴嘴出口溫度/℃22.3±0.222．5±0.223.4±0.327.0±0.424.5±0.323.6±0.2

噴嘴出口流速/u0/m/s57.4747.057.8629.3629.337.84

噴嘴出口雷諾數(shù)Re0/1041.441.473.623.677.335.03

表中

4．1．2各組實驗量測結(jié)果

實驗量測的時均流速及脈動流速相對值的典型結(jié)果見圖6。圖中的軸向流速脈動相對值是指，um為射流軸線時均流速。射流邊界位置由正反流邊界探針確定。邊界附近的時均流速參考皮托管量測進(jìn)行了修正（圖6a給出修正后的時均流速分布）。距噴嘴出口距離增大后，由于流速減小，正反流邊界探針分辨不出邊界，所以后面的斷面沒有給出時均流速分布，但給出了脈動流速分布。

圖6時均流速與軸向流速脈動相對值沿程分布圖

4．2堵頭管中受限射流的一般特性及分區(qū)

根據(jù)對實驗成果的整理分析，受限射流流場可分為四個區(qū)：I自由擴展區(qū)，Ⅱ受限擴展區(qū)，Ⅲ收縮區(qū)，Ⅳ零均流區(qū)，如圖7所示。

圖7堵頭管中受限射流流場示意圖

I自由擴展區(qū)，Ⅱ受限擴展區(qū)，Ⅲ收縮區(qū)，Ⅳ零均流區(qū)

在第一臨界斷面C1之前為自由擴展區(qū)，C1斷面之后，射流的擴展受到邊界和回流的限制，其卷吸周圍流體作用減弱，射流半徑流量的增加速率逐漸減慢，但總的趨勢仍是半徑增大，流量增加，故稱為受限擴展區(qū)，直到第二臨界斷面C2。此后，射流運動發(fā)生根本轉(zhuǎn)折，射流流線開始越出邊界產(chǎn)生回流，射流區(qū)流量沿程減少，在第二臨界斷面上射流區(qū)流量達(dá)到最大值，而射流半徑在C2斷面稍后達(dá)到最大值，然后逐漸縮小，稱為收縮區(qū)，收縮區(qū)末端位置（第三臨界斷面）由射流最遠(yuǎn)作用長度決定。可根據(jù)實際工程中對軸線流速的限值給出（詳見4.3.1）。第三臨界斷面和圓管堵頭之間是第Ⅳ區(qū)--零均流區(qū)，或稱尾部旋渦區(qū)，該區(qū)氣流的特點是處理"阻滯"和"堆積"狀態(tài)，形成一個或多個不穩(wěn)定的大尺度旋渦，旋渦誘導(dǎo)的瞬時流甚，時間平均流速為零，壓力均勻。

堵頭管中的射流，在第I和第Ⅱ區(qū)所卷吸的流體來自本身的回流，管內(nèi)空氣產(chǎn)生多次循環(huán)，排除的只是從噴口噴入的那部分空氣量。

4．3受限射流的時均流動特性

4．3．1受限射流軸線流速變化規(guī)律

射流進(jìn)入堵頭管后，由于動量的衰減和橫向傳遞，流速不斷降低。圖8是射流軸線上時均流速隨無量綱距離x/d0的沿程變化。

從圖中可見，軸線流速的變化取決于射流的受限度。當(dāng)受限度很大時，軸線流速急劇下降，射流相對射程x/d0隨著受限度減小而增大。而且隨著受限程度減小，軸線流速的變化曲線互相趨于接近，當(dāng)受限度很時，趨近于圖8中用虛線表示的自由射流軸線流速變化曲線，射流受了程度越大，其軸線流速變化曲線um/u0=f(x/d0)離自由射流的曲線越遠(yuǎn)。

圖8軸線流速隨x/d0的變化

圖9是用無量綱距離表示的軸線流速沿程變化。圖9a中曲線表示巴哈列夫經(jīng)驗公式計算值。從圖中可見，實測值與經(jīng)驗公式計算值在趨勢上是一致的，但在數(shù)值上，前者大于后者。巴哈列夫的受限射流軸線流速經(jīng)驗公式是把受限度引入自由射流軸線流速公式之后擬合實驗數(shù)據(jù)得到的。他采用的自由射流軸線流速公式為

（1）

式中<1時，≥1時，A=6.5a（2）

圖9軸線流速隨的變化

巴哈列夫把自由射流看作受限射流的極限情況，將受限度以e的指數(shù)形式引入式（1），得以受限射流軸線流速表示式：

（3）

式中的指數(shù)m即反映射流受限程度的影響，根據(jù)實驗數(shù)據(jù)，巴氏得到

（4）

由此可見1自由射流的軸線流速公式在巴氏受限射流軸線流速經(jīng)驗公式中起著重要作用。圖10比較了式（1）和Tollmien的理論解[1]，即，發(fā)現(xiàn)巴氏所采用的自由射流軸線流速公式（1）的計算值低于其它學(xué)者的研究成果。本文采用Tollmien理論解的數(shù)據(jù)擬合式（1）中的因數(shù)A，得到當(dāng)<1時（5a）A=6.85a當(dāng)≥1時（5b）于是經(jīng)過修正的受限射流軸線流速經(jīng)驗公式變?yōu)?/p>

（6）

式中A按式（5）選取。

圖10自由射流線流速沿程變化

圖8、圖9b中的曲線均是根據(jù)式（6）計算得到的。從圖中可見，實測值與之基本相符。

考慮到射流作用的末端，斷面上流速分布趨于均勻，軸線流速與斷面平均流速相差不大，在實際工程中可把軸線流速達(dá)到給定值ud時的長度定義為射流最遠(yuǎn)作用長度。ud根據(jù)需要選定，在通風(fēng)空調(diào)中ud一般為0.07～0.04m/s[9][13]，代入式（6）可以得到受限射流的最遠(yuǎn)作用長度，也就是第三臨界斷面的位置。

4．3．2受限射流流量的變化規(guī)律

根據(jù)實測斷面流速分布，可得受限射流的流量Qj。圖10是受限射流流量隨時無量綱距離的沿程變化。由于管壁邊界限制了射流的發(fā)展，受限射流的流量不像自由射流那樣沿程增加，而是先逐漸增大，到達(dá)第二臨界斷面之后，由于氣流的回轉(zhuǎn)，射流流量逐漸減小。在圖10中，從射流噴口到第一臨界斷面的自由擴展區(qū)，流量增加很快；隨后的受限擴展區(qū)，流量增加速度逐漸減慢并在第二臨界斷面處達(dá)到最大值；第二臨界斷面以后，射流流量逐漸減小。

射流流量的沿程變化在一定程度上反映了射流的卷吸性質(zhì)，從圖11a可見，射流受限越大，相對流量Qj/Q0最大值越小，卷吸量越小。圖11b是各個受限度射流無量流量的綜合結(jié)果。從圖中可見，不同受限度的實驗點比較接近，可用一條曲線表示。圖中曲線1是根據(jù)巴哈列夫經(jīng)驗公式計算得到的。從圖中可見實驗點與之偏差較大，歸因于巴哈列夫量測流速的儀器--電風(fēng)速儀在紊動強度較高時精度較差。

仍選用巴氏經(jīng)驗公式的形式，擬合實驗數(shù)據(jù)，得到

（7）

如圖11b中曲線2所示。圖11a中的曲線也是根據(jù)式（7）給出的。兩圖中小受限度的一些斷面與其它實驗點偏離的原因可能是小受限度時，噴口直徑很小，相應(yīng)的噴口流量Q0也小，這樣流速量測時較小的誤差可能導(dǎo)致Qj/Q0較大的誤差，故在擬合式（7）時沒有包括明顯有誤差之點。

圖11射流流量隨的沿程變化

根據(jù)實驗結(jié)果，射流在各種受限度下，第二臨界斷面均在同一無量綱距離c2=0.21處,在此斷面上，射流量達(dá)到最大

值：（8）

4．3．3受限射流射流區(qū)平均流速變化規(guī)律

根據(jù)受限射流流量及射流半徑的變化規(guī)律，可得到受限射流射流區(qū)平均流速的沿程變化規(guī)律，如圖12所示。從圖12b可

見，各受限度下射流平均流速無量綱綜合值的分布很有規(guī)律。由上述實驗成果，可得到射流區(qū)平均流速變化的經(jīng)驗公式

（9）

根據(jù)式（9）所作的曲線與實驗數(shù)據(jù)基本相符。

圖12射流區(qū)平均流速沿程變化

4．3．4受限射流回流區(qū)平均流速變化規(guī)律

如前所述，堵頭管中受限射流，由于邊界的限制，射流的外部存在回流。在通風(fēng)空調(diào)工程中，工作區(qū)通常位于回流區(qū)，這是因為回流區(qū)流動的流速和溫差已充分衰減，具有比較均勻穩(wěn)定流速場和溫度場，所以回流區(qū)平均流速是通風(fēng)空調(diào)設(shè)計的重要指標(biāo)，在實際工程中，往往要限制回流區(qū)平均流速最大值不得超過工作區(qū)風(fēng)速的允許值。

目前，在通風(fēng)空調(diào)工程設(shè)計中，巴哈列夫?qū)嶒灥玫降幕亓鲄^(qū)平均流速經(jīng)驗公式

（10）被廣泛采用[13，14]，但在文獻(xiàn)[9]中這個經(jīng)驗公式里的回流區(qū)平均流速定義不明確。斷面平均流速是相應(yīng)流量和相應(yīng)面積的比值，但文獻(xiàn)[9]卻先給出回流區(qū)平均流速在堵頭管長度上的變化。擬合得到式（10），然后由連續(xù)原理得到回流流量Qr（等于同一斷面的射流流量Qj）和式（10），給出了回流面積，未免令人費解。

本文定義每個斷面上的回流區(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ù)據(jù)的擬合曲線，所得的擬合公式為

（13）回流區(qū)平均流速從噴口開始逐漸增加，至第二臨界斷面c2=0.21達(dá)到最大值，然后又逐漸減小。在第二臨界斷面上，回流區(qū)平均流速最大值等于（14）

4．4受限射流的紊動特性

根據(jù)圖6中給出的堵頭管中受限射流軸向流速及相對值的沿程分布，可以對受限射流的紊動特性進(jìn)行初步探討。

圖14是不同受限度下射流軸線上隨x/d0的變化（在射流軸線上，ε1與相對紊動強度相同）。圖中虛線取自參考文獻(xiàn)[15][16]中自由射流數(shù)據(jù)。

圖14射流軸線上相對紊動強度沿程變化

從圖中可知，在自由擴展區(qū)，不同受限度射線軸線上相對紊動強度變化與自由射流基本一致；而自由擴展區(qū)之后，受限射流細(xì)線上紊動強度沿程持續(xù)增加，并不是趨于常數(shù)。分析其原因，無論是自由射流還是堵頭管中的受限射流，在射流軸線上ui/xj=0，因而軸線上不產(chǎn)生紊動，軸線上的紊動是由上游的輸運和附近的擴散引起的。在受限條件下，由于回流的存在，射流斷面上流速分布更不均勻，軸線附近紊動加大，因而擴散到軸線上的紊動也相應(yīng)地增加；同時，主流的波動也是導(dǎo)致相對紊動強度增加的原因。

從圖6可見，不同受限度射流軸向流速脈動ε1在斷面上的徑向分布情況，氣流從噴口噴入堵頭管中，在勢流核心區(qū)沒有完全消失的斷面，核心區(qū)內(nèi)ε1很小；隨著徑向坐標(biāo)r的增加，ε1先是逐漸增大，增大到一定程度后開始減小，直至射流邊界附近；此后r再增大，ε1基本保持為常數(shù)。結(jié)合斷面上時均流速的分布，可以得到，從射流軸線開始，隨著r的增加，相對紊動強度先是和ε1一樣逐漸增大，但當(dāng)ε1達(dá)到最大值并開始減小時，由于時均流速u的減小，ε將趨近于∞。而在射流邊界以外的回流區(qū)，雖然ε1保持在一個較低水平，但由于時均流速同樣較低，相對紊動強度ε值仍然較大。從圖6中還可以看出，隨著軸向距離的增加，整個斷面的軸向流速脈動均增大，并逐漸趨于均勻。

圖15給出的是d0=20mm，=10.6的射流各斷面上回流區(qū)軸向流速脈動相對于回流區(qū)平均速度vr的沿程變

化，從圖中可以看出，隨著回流越接近出口，ε2值越大，這主要是vr迅速減小的結(jié)果。

圖15回流區(qū)ε2沿程變化

5結(jié)語

5．1以往對堵頭管中射流時均特性的研究成果定性上是正確的。但是，由于實驗多在早期完成，受到量測手段的局限，定量存在一定的誤差。本文采用了熱膜流速儀以及專門設(shè)計的正反流邊界探針，使得研究成果具有較高的可靠性。

5．2堵頭管中射流的發(fā)展可分為四個區(qū)：自由擴散區(qū)、受限擴展區(qū)、收縮區(qū)和零均流區(qū)。分區(qū)的第一臨界斷面和第二臨界斷面分別在c1=0.09和c2=0.21處，第三臨界斷面則由射流最大作用長度決定。

5．3對巴哈列夫關(guān)于受限射流軸線流速沿程變化的經(jīng)驗公式進(jìn)行了修正，見式（6）受限射流、流量射流區(qū)平均速度沿程變化可分別按經(jīng)驗公式（7）和（9）計算。

5．4通過定義回流區(qū)面積等于堵頭管面積與射流區(qū)面積之差，給出了通風(fēng)空調(diào)設(shè)計的重要指標(biāo)--回流區(qū)平均速度的經(jīng)驗公式（13），澄清了巴哈列夫?qū)亓鲄^(qū)平均速度定義的不明確之處。

5．5在實驗所包括的協(xié)諾數(shù)范圍內(nèi)(Re0>104),受限射流的流動與雷諾數(shù)無關(guān)，這主要是由于在高雷諾數(shù)情況下，旋渦結(jié)構(gòu)的變化主要在于最小旋渦尺度的減小，而大旋渦尺度基本不變，由于卷吸強度受尺度旋渦的影響，故受限射流的卷吸強度亦不隨Re0而變化。但是，受限度影響流動的發(fā)生發(fā)展以及卷吸強度的大小，受限度越大，射流所卷吸的流量越小。

5．6在受限條件下，由于回流的存在，射流的紊動在整個堵頭管中都很劇烈，尤其是在收縮區(qū)和零均流區(qū)。預(yù)示依據(jù)研究堵頭管中的射流，無論是理論分析還是數(shù)值計算，必須對運動方程中流速脈動項給予充分的重視。

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余常昭，紊動射流，北京：高等教育出版社，1993。

7KovasznayLGG.StructureandMechanismsofTurbulence.EdbyFiedlerH.ProcSymponTurbulence,Berlin,1997,vol1:1～18.

8LauferJ.AnnRevFluidMech.1975,7:307～326.

9巴哈列夫BA，特羅楊諾夫斯基著，寧德平譯，集中送風(fēng)式采暖通風(fēng)設(shè)計計算原理。北京：中國工業(yè)出版社，1965。

10HusssainAKMFandRamjeeV.Effectsoftheaxisym-metricshapeonincompressibleturbulentflowJFluidEng,1976,98:55～69.

11王世蛤，陳睿，姚志化，二元火焰穩(wěn)定器冷態(tài)流場測量，工程熱物理學(xué)報，1983，4（1）：74～79。

12ГорбуноваГМ.СгaбициэaцияПлaмeнииРaзвитиеПроцeccaСгopaиявТурбуцнентномПотоке.ГоcyцapcтвeннoeНаучно-ТexничeckoeИэцательство,Оборонтиэ,1961.

13陸耀慶主編，供暖通風(fēng)設(shè)計手冊，北京：中國建筑工業(yè)出版社，1986。

14周謨?nèi)手骶帲黧w力學(xué)泵與風(fēng)機，第二版，北京：中國建筑工業(yè)出版社，1985。

15WygnanskiIandFiedlerH.Somemeasurementsinself-re-servingjet.JFluidMech,1969,38:577～612.

16李永祥，苑明順，李春華，熱膜流速儀在水流和摻氣水流中應(yīng)用研究，清華大學(xué)水電系研究報告，199。