熱水噴泉現(xiàn)象及其理論研究

李 強(qiáng) 覃志遠(yuǎn) 陳 聰 熊佳芬 馮公警

(南京航空航天大學(xué) 1理學(xué)院， 2航天學(xué)院,江蘇 南京 211100)

1 提出問(wèn)題

本文的研究?jī)?nèi)容是來(lái)源于2016年IYPT(International Youth Physics Tournament,國(guó)際青年物理學(xué)家競(jìng)賽)的一道賽題: Problem No.7 “Hot Water fountain” Partially fill a Mohr pipette with hot water. Cover the top of the pipette with your thumb. Turn the tip upwards and observe the fountain exiting the tip. Investigate the parameters describing the height of the fountain, and optimize them to get the maximum height. (用熱水部分裝填一個(gè)莫爾吸量管,用大拇指蓋住其上端并倒置,可觀察到從尖部噴出的熱水噴泉。請(qǐng)研究決定熱水的最大噴射高度的參量,并改變它們以獲得最大的高度。)[1]

熱水噴泉現(xiàn)象的研究在地?zé)釃娙拈_(kāi)發(fā)利用、噴水推進(jìn)技術(shù)研究等方面有潛在應(yīng)用價(jià)值。但到目前為止,國(guó)內(nèi)外對(duì)該現(xiàn)象定量研究的報(bào)道較少。劉建曉等人研究了莫爾吸量管中吸入部分熱水產(chǎn)生的噴泉現(xiàn)象,應(yīng)用伯努利方程得到最大噴射高度的表達(dá)式,討論了吸水高度、平衡溫度及噴嘴孔徑對(duì)最大噴射高度的影響。他們通過(guò)多項(xiàng)式擬合求得平衡溫度與吸水高度的關(guān)系式[2]。付偉娟等同樣應(yīng)用伯努利方程并具體考慮莫爾管管口形狀建立理論模型[6]。由于熱水噴泉現(xiàn)象屬于一種有多變量因素(如:熱水溫度、吸水高度以及莫爾管規(guī)格等)綜合作用的非定常流動(dòng),而伯努利方程只適用于定常流體的運(yùn)動(dòng)問(wèn)題,因此直接應(yīng)用伯努利方程處理該現(xiàn)象過(guò)于理想化。為此,我們從噴水過(guò)程中的熱水量守恒以及動(dòng)能定理出發(fā)，建立了一個(gè)既適用于定常流動(dòng)又適用于非定常流動(dòng)的更普遍的理論模型。在此基礎(chǔ)上,通過(guò)設(shè)計(jì)不同實(shí)驗(yàn)參量,從實(shí)驗(yàn)上對(duì)上述模型進(jìn)行了驗(yàn)證。結(jié)果表明:實(shí)驗(yàn)測(cè)量結(jié)果與理論模型的數(shù)值計(jì)算結(jié)果符合較好；熱水溫度、吸水高度以及莫爾管規(guī)格是影響最大噴射高度的主要因素。

2 實(shí)驗(yàn)研究與分析

2.1 實(shí)驗(yàn)裝置

實(shí)驗(yàn)中用到的器材裝置有:鐵架臺(tái)、量尺、溫度計(jì)以及莫爾管和高速攝像機(jī)。如圖1所示。

圖1 實(shí)驗(yàn)裝置圖

莫爾管是一種有上、下兩個(gè)端口的細(xì)長(zhǎng)玻璃儀器,常用于測(cè)量液體體積。管身上的刻度線用于標(biāo)記液體體積大小,兩個(gè)端口一大一小,便于液體進(jìn)出。

2.2 實(shí)驗(yàn)方法

為了和后面理論建模與分析保持一致,對(duì)實(shí)驗(yàn)中可能出現(xiàn)的物理量采用如下定義:H為莫爾管總長(zhǎng)度;m為口徑比,即莫爾管較細(xì)端口面積和較粗端口面積的比值；S2為莫爾管細(xì)口的面積;r為莫爾管粗口的半徑;h0為實(shí)驗(yàn)過(guò)程吸入的熱水高度;hmax為噴水過(guò)程中水柱達(dá)到的最大高度。

2.3 實(shí)驗(yàn)過(guò)程

實(shí)驗(yàn)準(zhǔn)備時(shí),先將莫爾管較小的管口朝下豎直浸在熱水池中,利用洗耳球從其上部的大管口吸氣,可使其下部的小管口能從熱水池中吸入一定量的熱水,隨后用橡膠塞堵住大管口。實(shí)驗(yàn)開(kāi)始時(shí),用旋轉(zhuǎn)夾將裝好熱水的莫爾管小管口豎直朝下采用可旋轉(zhuǎn)夾持器固定在鐵架臺(tái)上,隨后將莫爾管迅速轉(zhuǎn)動(dòng)180°(如圖2所示),隨即熱水柱從較小管口豎直向上噴出。為方便測(cè)量噴射高度,我們?cè)跓崴屑尤肓松倭亢谏?/p>

圖2 莫爾管夾持圖

2.4 實(shí)驗(yàn)結(jié)果

(1) 實(shí)驗(yàn)時(shí),首先使用高速攝像機(jī)觀察了從管口噴出去的熱水的形狀,發(fā)現(xiàn)噴出去的熱水呈完整的水柱狀。

(2) 在一定條件下實(shí)驗(yàn)，觀察到從莫爾管中噴出的水柱高度隨時(shí)間先從零開(kāi)始增加達(dá)到一個(gè)最大值后減小。實(shí)驗(yàn)結(jié)果如圖3所示。

圖3 噴射高度隨時(shí)間變化實(shí)驗(yàn)結(jié)果(吸水高度為0.17m,熱水溫度為73℃,環(huán)境溫度為19℃)

(3)其他變量(熱水溫度,莫爾管規(guī)格,環(huán)境溫度)相同，研究吸水高度對(duì)最大噴射高度的影響時(shí),實(shí)驗(yàn)得到的數(shù)據(jù)如表1、表2所示。

表1 最大噴射高度hmax隨吸水高度h0變化的實(shí)驗(yàn)數(shù)據(jù)

表2 最大噴射高度hmax隨吸水高度h0變化的實(shí)驗(yàn)數(shù)據(jù)

由表1和表2的實(shí)驗(yàn)數(shù)據(jù)可看出,在兩種不同的環(huán)境溫度下,莫爾管中熱水噴射的最大高度隨著初始時(shí)莫爾管中吸水高度逐漸增加都經(jīng)歷了先減小后增大的變化。

(4) 其他變量(吸水高度,莫爾管規(guī)格,環(huán)境溫度)相同,考慮熱水溫度對(duì)最大噴射高度的影響時(shí),實(shí)驗(yàn)得到的數(shù)據(jù)如表3所示。

表3 最大噴射高度hmax隨熱水溫度T變化的實(shí)驗(yàn)數(shù)據(jù)

3 理論建模與分析

為了解釋上述實(shí)驗(yàn)現(xiàn)象及結(jié)果,對(duì)噴射過(guò)程進(jìn)行了如下的理論分析和建模。

圖4 噴水裝置示意圖

此過(guò)程可看作是等容過(guò)程,因此有

(1)

(2)

其中,c與熱水-環(huán)境溫差有關(guān),還與初始時(shí)加入的熱水量以及管子規(guī)格有關(guān)。

將式(2)代入式(1),得

(3)

設(shè)噴水過(guò)程開(kāi)始時(shí)刻t0=0，t時(shí)刻,各部分熱水柱長(zhǎng)度分別為h(t),h1(t),噴水過(guò)程中任意時(shí)刻下部空氣壓強(qiáng)為p(t),體積為V(t)=πr2[H-h1(t)]。

由于熱水可認(rèn)為不可壓縮,因此水的總體積不變

在噴射過(guò)程中,可認(rèn)為管內(nèi)空氣經(jīng)歷絕熱過(guò)程:

(6)

在噴射過(guò)程中的t～t+dt時(shí)間段內(nèi),管中空氣柱的壓強(qiáng)和外界大氣壓以及水柱的重力做的總

功等于熱水柱總動(dòng)能的增加

(7)

將式(6)代入并在等式兩邊同時(shí)除以dt,得

(8)

式(8)即為描述熱水噴射高度隨時(shí)間變化的微分方程式。此式很難求出解析解,但可做數(shù)值求解,以便與實(shí)驗(yàn)結(jié)果進(jìn)行比較。

4 實(shí)驗(yàn)與理論比較

算得一系列c值,見(jiàn)表4。

表4 不同情況下的c值

4.1 噴射高度隨時(shí)間變化

考慮吸水高度為0.17m,熱水溫度為73℃,環(huán)境溫度為19℃的情況

根據(jù)式(8)代入?yún)?shù)做數(shù)值計(jì)算得出的噴射高度隨時(shí)間變化如圖4所示。

圖5 噴射高度隨時(shí)間變化理論結(jié)果(吸水高度為0.17m,熱水溫度為73℃,環(huán)境溫度為19℃)

由圖4可以看出,噴射高度先從零開(kāi)始增加達(dá)到一個(gè)最大值后逐漸減小。這與實(shí)驗(yàn)結(jié)果中(2)相符。

4.2 最大噴射高度與吸水高度的關(guān)系

考慮環(huán)境溫度為19℃，熱水溫度為83℃的情況

由圖6可得出,在兩種不同的環(huán)境溫度下,莫爾管中熱水噴射的最大高度隨著初始時(shí)莫爾管中吸水高度逐漸增加都經(jīng)歷了先減小后增大的變化。這是由于一方面當(dāng)加入熱水較少時(shí),其自身重力越小,從而對(duì)莫爾管內(nèi)空氣膨脹的阻力影響較小,最大噴射高度就較高；另一方面,熱水越多,空氣在預(yù)熱結(jié)束后溫度越高,壓強(qiáng)越大,從而對(duì)外做功能力越強(qiáng),這也可以從c值表4得到反映。

圖6 不同吸水高度對(duì)最大噴射高度的影響(環(huán)境溫度:19℃,熱水溫度:83℃)

從圖6可以看出,最大噴射高度隨初始時(shí)吸水高度變化很明顯,說(shuō)明吸水高度是一個(gè)重要的影響因素。

4.3 最大噴射高度與熱水溫度的關(guān)系

考慮初始時(shí)莫爾管吸水高度為0.12m,環(huán)境溫度為19℃的情況:

由圖7可得出,最大噴射高度隨著熱水溫度的升高快速地增加,這是由于熱水溫度越高,管內(nèi)空氣被充分預(yù)熱完成后的溫度就越高,從而壓強(qiáng)越大,對(duì)外做功的能力就越強(qiáng),水柱最大噴射高度越高。

圖7 不同熱水溫度對(duì)最大噴射高度的影響(吸水高度:0.12m,環(huán)境溫度:19℃)

由圖7可見(jiàn),最大噴射高度隨熱水溫度增加有明顯的變化,因此熱水溫度也是一個(gè)重要的影響因素。

5 結(jié)論

本文以莫爾管為基本實(shí)驗(yàn)裝置,采用控制變量法,通過(guò)對(duì)一系列實(shí)驗(yàn)數(shù)據(jù)的測(cè)量,探討了熱水噴泉現(xiàn)象中各種因素對(duì)水柱最大噴射高度的影響,得到了最大噴射高度隨熱水溫度、吸水高度等參量的變化規(guī)律。在此基礎(chǔ)上,基于噴水過(guò)程中的動(dòng)能定理,建立了一個(gè)最大噴射高度與相關(guān)量關(guān)系的理論模型,數(shù)值計(jì)算了上述實(shí)驗(yàn)條件下的最大噴射高度,并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較,結(jié)果表明:實(shí)驗(yàn)測(cè)量結(jié)果與理論模型的計(jì)算結(jié)果符合得較好。該研究進(jìn)一步揭示了熱水噴泉現(xiàn)象中的運(yùn)動(dòng)機(jī)理,研究結(jié)果對(duì)開(kāi)展熱水噴泉的相關(guān)應(yīng)用有一定的指導(dǎo)作用。

[1] IYPT. 2016年IYPT競(jìng)賽參考文獻(xiàn)[EB/OL]http://archive.iypt.org/problems. IYPT. 2016 IYPT competition reference[EB/OL]http://archive.iypt.org/problems.

[2] 劉建曉,鄭永春,史宮會(huì)，等.熱水噴泉現(xiàn)象的理論研究[J].物理實(shí)驗(yàn) 2016,36(4):23-26. LIU JianXiao, ZHENG Yongchun, SHI Gonghui, et al. Research on the hot water fountain[J]. Physical Experiment, 2016, 36(4): 23-26. (in Chinese)

[3] 植田辰洋,姚普明,何乃翔，等.氣、液二相流的流動(dòng)和傳熱[J].天津商學(xué)院學(xué)報(bào),1985(1):1-18. (in Chinese) ZHI Tianchenyang, YAO Puming, He Naixiang, et al. Flow and heat transfer in gas-liquid two-phase flow[J]. Journal of Tianjin University of Commerce, 1985(1): 1-18. (in Chinese)

[4] 趙近芳,王登龍.大學(xué)物理簡(jiǎn)明教程[M].2版.北京:北京郵電大學(xué)出版社,2013:67-76.

[5] 李復(fù).可壓縮流體的伯努利方程[J].大學(xué)物理,2008,27(8):15-18. LI Fu. Bernoulli’s equation for compressible flow[J]. College Physics. 2008, 27(8): 15-18. (in Chinese)

[6] 付偉娟,聶云漢,房振全，等.莫爾吸量管熱水噴泉現(xiàn)象的研究[J].物理實(shí)驗(yàn) 2017(6):37-6. FU Weijuan, NIE Yuanhan, FANG Zhenquan, et al. Study on the water fountain pipette moire phenomenon[J]. Physical Experiment, 2017(6): 37-6. (in Chinese)