尿素溶液制備冰漿的流動特性研究

張盛世 趙彥杰 李兆寧

(1 天津大學中低溫熱能高效利用教育部重點實驗室 天津 300350； 2 天津大學機械工程學院 天津 300350)

冰漿是指某類鹽溶液和冰晶粒子的共晶混合物，是一種呈泥漿狀的懸浮液。它具有良好的流動特性和較高的能量密度，既可作為相變冷庫材料，又可作為第二制冷劑來輸送冷能。在空調系統、食品加工、冷藏運輸、醫藥等領域應用廣泛[1]。

冰漿通常由鹽溶液而非純水制取產生。溶液中溶質的不同可以改變冰點溫度，降低結冰附著力并減小冰的粒徑。低濃度溶液的凝固點降低規律遵循溶液的“依數性”原則。溶液在凍結過程中溶質發生遷移，因此在固液兩相區，可以觀測到固相區中包含大量的小冰粒[2]，而且溶液中溶質的化學鍵具有阻止冰晶生長的作用[3]。所以冰層附著力低，生成的冰漿更加均勻，抑制了冰堵的發生和換熱的惡化。溶液中常見的溶質主要有氯化鈉、乙醇、乙二醇和丙二醇等。

在設計熱能存儲系統之前，首先應重點了解冰漿的流動特性。冰漿作為固液兩相流體，流動特性復雜，受多種因素影響，包括溶質變化、溶液濃度、流速、管道尺寸和含冰率(ice packing factor, IPF)等。國內外對冰漿的研究已有20多年的歷史，其中冰漿的流動特性已得到了廣泛的研究和應用。Tian Qiqi等[4]提出了3種管內冰漿的流動模式：懸浮流、移動床流和固定床流。懸浮流在高流動速率下發生，冰粒完全懸浮并均勻分布。移動床流在較低的流動速率下發生，流動過程中漿體被分為兩層，頂層為移動床層，而底層仍然是非均質流層。固定床流則在非常低的流動速率下發生，在該流動模式下，載流體因無法承載所有的移動床層而在管道頂部形成了固定床層。V. Ayel等[5]回顧了4種流變模型，即Bingham模型、 Casson模型、 Ostwald模型 和 Herschel Buckley模型，這些模型有助于預測不同條件下管中冰漿的流動行為。H. Kumano等[6]主要研究了乙醇溶液制取冰漿的流動特性，并建立了管道摩擦阻力系數與雷諾數之間的關系，文中指出，在層流條件下，可以將冰漿看做假塑性流體進行處理。H. Kumano等[7]分析了初始溶液濃度對冰漿流動特性的影響，結果表明，利用初始質量分數為5%和10%的乙醇溶液制取的冰漿，其流動特性同假塑性流體的流動特性相近，而由質量分數為2%的乙醇溶液制取冰漿則可看做牛頓流體。B. D. Knodel等[8]研究了乙二醇溶液制取冰漿的流動特性，通過實驗測量了冰漿流動過程中的再分層特性，指出冰漿流動過程中的摩擦壓降隨著含冰率的增加而降低。劉巨保等[9]研究了氯化鈉溶液制取的冰漿在水平管中的流動特性，基于兩相流和顆粒動力學理論，分別研究了管徑、管長、IPF、流速對冰漿流動過程中壓力損失的影響。

然而到目前為止，研究人員所選取的制取冰漿的溶液及實驗環境各不相同，因此對于冰漿流動特性的研究還未達成統一的共識。本課題組之前的研究[10]指出，尿素溶液的相變潛熱均高于同濃度的常規制備冰漿的溶液，如乙醇溶液和氯化鈉溶液。質量分數為10%的尿素溶液的相變潛熱高達311.2 J/g，略低于純水的潛熱，而且尿素溶液的相變潛熱隨著濃度的降低愈發接近于純水的潛熱。同時尿素溶液完全凍結后形成的冰層質地松軟，具有多孔結構，極易被液流的剪切力剝落，可以自發產生冰漿，并且目前還未發現有研究學者對尿素溶液制取冰漿的流動特性進行研究。因此，本文以尿素溶液所制取的冰漿作為研究對象，測量了不同初始質量濃度尿素溶液在常溫狀態和冰水混合物狀態的運動黏度，觀測了冰漿中冰粒分布情況及大小，分析了管徑、 IPF和流速對摩擦阻力系數的影響，并根據實驗結果闡明了冰漿在水平管中的流動特性。

1 實驗裝置及方法

1.1 實驗裝置

冰漿的高效制備和冰漿在水平管內的流動特性對于指導蓄冷系統的設計至關重要，為此我們沿用本課題組之前設計的冰漿制備和流動特性測試實驗臺[11]，實驗裝置如圖1所示。實驗系統主要包括3部分：冷卻液循環回路、冰漿制備循環回路和冰漿的流動特性實驗測試段。圖中虛線箭頭所示為冷卻液循環回路。恒溫箱為系統制備冰漿提供冷量，恒溫箱的最大制冷量為7.8 kW。恒溫箱中的冷卻液采用體積比為1∶1的乙二醇溶液，最大冷卻溫度可達 -33 ℃。在整個循環過程中，乙二醇溶液作為攜帶冷量的冷卻液，流經制漿罐的夾層，并將冷量轉移至罐中的尿素溶液中。

圖1 實驗裝置

在冰漿制備循環回路中，通過過冷法在制漿罐中制備冰漿。制漿罐設計為套筒式，內筒盛有質量分數為5%的尿素溶液，容器的容積為10 L；外筒循環流動乙二醇冷卻液，采用左右并進、下進上出的模式，以起到強化換熱的效果，并通過調節來自恒溫箱的乙二醇冷卻液溫度來控制制漿罐中溫度。同時在制漿罐內筒中安裝一個電動攪拌器，用于打破尿素溶液的過冷狀態并除去冰漿中的空氣。攪拌器的轉速范圍為0～300 r/min。攪拌罐頂部布置一個有機玻璃頂蓋，頂蓋上有兩個小孔，用于插入兩根T型熱電偶，以測量制備過程中溶液和冰漿的溫度變化，熱電偶的精度為0.001 ℃。系統中封閉的實線箭頭回路為溶液和冰漿的循環回路，整個循環管道均為DN15的不銹鋼直管和軟管，循環的動力裝置為帶有逆變器的小型離心泵，額定功率為370 W，揚程范圍為7.5～11.6 m，體積流量范圍為1.2～6 m3/h。通過對離心泵進行頻率轉換來控制流量。系統中采用Coriolis 質量流量計測量和記錄溶液和冰漿的質量流量和密度變化，流量測量量程為0～3 000 kg/h，測量精度為0.1%；密度測量范圍為0.3～3.000 g/cm3，測量精度為±0.002 g/cm3。流動速率用于計算冰漿流動過程中的雷諾數Re，冰漿的密度用于計算IPF。離心泵工作過程中內部產生的熱量對冰漿中IPF的降低無法避免也不能被忽略，所以將Coriolis質量流量計安裝在離心泵的出口處，以提高質量流量和密度數據的可信度。冰漿的溫度測量范圍為-200～350 ℃，測量精度為±0.01 ℃，工作壓力為0～4.0 MPa。在冰漿進入測試段之前首先會流經一段長度為1.5 m的充分發展段，以確保進入測試段的冰漿處于穩定狀態，減小測量誤差。測試段長1.2 m，并在測試段的兩端分別布置一個壓力傳感器，量程為0～0.6 MPa，精度為0.01%，輸出信號為4～20 mA，同時在壓力傳感器后并聯一個阻值為1 Ω的電阻，電阻精度為0.02%，用于將電流信號轉化成電壓信號，因為數據采集器僅接受電壓信號。又因4～20 mV電壓與0～0.6 MPa的壓力數據呈線性關系，所以可通過電壓信號線性計算出壓力數據，兩個壓力相減即可得到壓降Δp0。在壓降測量實驗結束后，通過溶液儲罐收集冰漿。測試段采用的不銹鋼管分為內徑4、6、8 mm三種，用于研究不同管徑內冰漿的流動特性。由于測試部分的長度有限，在測試過程中IPF的變化很小。此外，實驗設備中的制漿罐和所有管道均采用橡膠絕緣棉進行包裹，有效減少熱量損失。

1.2 制取冰漿

圖2 尿素溶液的過冷度的變化

本文采用循環和機械攪拌的方法制取冰漿，質量分數為5%的尿素溶液在套筒式攪拌罐內被不斷地降溫和攪拌，并為了防止冰堵和實現冰漿的制備，制取冰漿的過程中溶液在離心泵的作用下一直處于循環流動狀態。隨著溫度的不斷降低，溶液在攪拌和流動的作用下迅速接觸過冷態，之后冰漿大量生成。隨著時間的延續和冷卻液溫度的降低，冰漿的含冰率會不斷增加。本實驗主要制備含冰率為5%～30%的冰漿，并研究其流動特性。尿素溶液的過冷度主要受攪拌器轉速和冷卻速度的影響，如圖2所示。在冰漿的制備過程中，首先設定80、120、150 r/min三種不同的攪拌轉速，圖2(a)所示為制漿罐內壁面溫度Tw=-5.2 ℃時，質量分數為5%的尿素溶液在不同攪拌轉速下的冰漿的制取過程。可以看出，攪拌轉速的提高促進了冰漿的生成，并促進了過冷態的解除，三種轉速對應的過冷狀態的最低溫度分別為-1.74、-1.55、-1.34 ℃。圖2(b)所示為攪拌轉速為 150 r/min 時，質量分數為5%的尿素溶液在不同壁面冷卻溫度Tw下冰漿的制取過程，冷卻溫度的降低促進了冰漿的生成，但抑制了過冷態的解除，隨著Tc和Tw的降低，過冷態的最低溫度也會略有下降。而且，較低的壁面溫度Tw也會使壁面冰層附著加劇，從而會惡化傳熱，影響冰漿的持續生成。

圖3(a)所示為Tw=-8.6 ℃時出現冰漿30 min之后制漿罐內的狀態。制漿罐壁面附著了一層致密的冰層，占據了壁和攪拌器葉片之間的空間，傳熱惡化，不利于冰漿的制備。圖3(b)所示為Tw=-4.1 ℃時出現冰漿30 min之后制漿罐內的狀態。制漿罐內壁面幾乎無冰層附著現象，有利于冰漿的制備。通過實驗發現，在質量分數為5%的尿素溶液的冰漿制備過程中，由于攪拌和循環的作用，過冷度很低，僅為1～2 ℃，而且生成的冰漿更加均勻，可以認為利用鹽溶液的攪拌降溫制取冰漿是一種良好的制取蓄冷用冰漿的方法。

圖3 冰漿照片

冰漿的制備過程持續進行，直至制備出的冰漿滿足某一特定含冰率(IPF)的要求。IPF的計算式[6]為：

(1)

圖4 冰漿中的冰粒分布

式中：ρs為初始溶液的密度，g/cm3；ρsl為冰漿的密度，g/cm3，兩者均通過質量流量計測量獲得；ρi為純冰的密度，取0.9 g/cm3。

圖4(a)所示為不同含冰率的冰漿的形態實物圖，可以看出冰粒的形態類似于紅細胞，利用放大50倍的工業相機對含冰率為0.126、0.181和0.258的冰漿形態分別進行拍攝，并利用像素點測量的方法對冰粒的尺寸進行精確測量，取照片中每一個冰粒的最大尺寸和最小尺寸的平均值進行統計，繪制出冰漿中冰粒直徑的頻率分布圖，如圖4(b)所示。由圖4(b)可知，三種不同含冰率下冰粒的尺寸多分布在0.10～0.25 mm內，平均尺寸約為0.15 mm。對于不同含冰率的冰漿，冰粒的尺寸無較大區別。

1.3 尿素溶液的性質

隨著溫度的不斷降低，溶液將在其平衡凍結溫度下進行凍結，導致殘留溶液中溶質濃度提高，使冰和溶液之間的相平衡得以維持。在冰漿的流動特性的研究中，需要測量尿素溶液的各項性能，其中Re是一個非常重要的參數，而在確定Re的過程中，運動黏度是計算Re的關鍵參數。采用Ubbelohde黏度計對不同初始質量濃度尿素溶液在常溫狀態和冰水混合物狀態的運動黏度進行測量。

圖5所示為不同質量濃度尿素溶液在冰點狀態的運動黏度關系圖。整個實驗在過冷狀態下進行。隨著尿素溶液質量分數由0 %增至32.5 %，冰點溫度從0 ℃線性降至-11.5 ℃，冰點狀態的運動黏度從1.648 mm2/s增至6.736 mm2/s。通過對比可知，固液兩相的流體的運動黏度大于純液體的運動黏度。圖5中的公式為冰點狀態運動黏度與溫度之間的擬合公式。

1.4 實驗條件

實驗通過測量不同管徑、Re和IPF下的壓降來研究水平管中冰漿的流動特性。實驗條件如表1所示。

表1 實驗條件

2 結果與討論

2.1 流動壓降

冰漿本質上為一種固液兩相混合流體，不同含冰率的冰漿流動壓降較單相流體更為復雜，對冰漿流動壓降的分析是研究其流動特性的基礎。圖6所示分別為IPF為0.06、0.14和0.24的三種不同管徑的流動壓降與Re的關系。由圖6可知，冰漿的流動壓降與IPF呈正相關，壓降隨IPF的增加而增加，同時也隨Re的增加而增加，但隨著測試管徑的增加，壓降逐漸減小。冰漿流動壓降對參數的靈敏等級為：管徑>IPF>Re。圖7所示分別為在內徑為6 mm的測試段內，不同IPF的冰漿分別在層流和湍流狀態下流動壓降與Re的關系。在層流條件下，不同IPF的冰漿，其壓降隨Re的增加而增大。當從層流過度至湍流的過程中，不同IPF的冰漿壓降接近，在湍流條件下，不同IPF的冰漿，其壓降繼續隨Re的增加而增大。同時，對于高IPF的冰漿，湍流的壓降增長的速率較大，同時高IPF的流動壓降較大。流態決定了冰漿的基本流動特性。在層流和湍流條件下，IPF對高Re狀態下的壓降影響更為顯著。說明Re越大，冰漿的IPF對壓降的影響作用越大。

圖6 不同管的壓降與Re之間的關系

圖7 不同IPF冰漿在層流和湍流條件下的壓降對比

2.2 摩擦阻力系數

單純的壓降數據分析很難研究不同Re和不同管徑內的冰漿與單相流的流動特性之間的差別。而摩擦阻力系數λ是冰漿在管內流動過程中一個非常重要的參數。分析冰漿與純尿素溶液的管內流動摩擦阻力系數的比值偏離1的程度與IPF和Re的關系，能更加直觀地看出作為一種固液兩相流體的冰漿與牛頓流體的流動差別，對冰漿蓄冷系統漿體的高效輸送和能耗分析具有重要意義。本文研究了冰漿的實驗摩擦阻力系數結果與尿素溶液理論值之比，以闡明尿素溶液產生的冰漿的流動特性。

式(2)和式(3)分別為層流和湍流狀態下牛頓流體的摩擦阻力系數的理論公式，式中的Re通過式(4)計算得出。冰漿流動過程中管道摩擦阻力系數的實驗值λ由式(5)計算。

λlam=64/Re

(2)

λtur=0.316 4Re-1/4

(3)

Re=ud/μf

(4)

λ=2DiΔp/(Lρu2)

(5)

式中：u為冰漿的流速，m/s；Di為測試段管徑m；μf為冰漿的運動黏度，mm2/s；Δp為流動壓降，Pa；L為測試段管道的長度，m。

圖8所示為不同Re條件下，冰漿管內流動摩擦阻力系數的實驗值λ與把冰漿作為牛頓流體利用經驗公式得到的理論值λ0(λlam、λtur)的比值和IPF的關系。結果表明，在所有條件下，λ/λ0均隨IPF的增加而增大，λ/λ0偏離1的程度整體隨Re的增加而降低。層流狀態(Re=1 500)下，λ/λ0偏離1的程度遠大于湍流狀態(Re=3 500、6 000、8 500)。以管徑為8 mm為例，當Re=1 500且IPF=0.28時，在層流條件下λ/λ0為5.5。當Re=3 500且IPF=0.28時，在湍流條件下λ/λ0為2.2。當然這與λ0在計算過程中使用的經驗公式不同有直接關系，但根本原因可用Doron的顆粒流動理論[12]來解釋，當冰漿流速較大時，在強湍流狀態下，冰晶粒子保持充分的均勻懸浮狀態，冰漿接近單相牛頓流體；流速進一步降低，湍流度則降低，冰粒比水輕，故漂浮在管道頂部；隨著平均流速繼續減小，進入層流區域，冰晶粒子還會在管道頂層形成滑動的沉積物，稱為移動層，同時管道最底部形成固定層，冰漿偏離單相牛頓流體。即Re越大，冰晶粒子和溶液混合更均勻；而低流速下，由于分層現象的出現，冰晶粒子多分布在管道的底部和頂部，增加了冰漿中固相與管內壁的摩擦，黏性摩擦阻力相對更大，從而相對高流速、高Re的狀態更加偏離單相流。對比三種不同的管徑可知，最小管徑為4 mm中，不同Re的λ/λ0均小于其他管徑的結果。原因是在小管徑內，相同Re下，流速較大，冰粒在管內流動時相對大管徑內分布更加均勻，分層現象較弱，結合顆粒流動理論，所以小管徑內冰漿的流動相對大管徑更接近于單相流，故λ/λ0更接近1。

圖8 管道摩擦阻力系數與IPF間的關系

2.3 范寧摩擦阻力系數和修正雷諾數

分析摩擦阻力系數比，發現無法將冰漿作為牛頓流體去分析其流動特性，因此，本文將冰漿看做擬塑性流體來分析其流動特性。對于固液兩相的非牛頓流體，由于流動剪切應力τw和剪切速率γw之間非線性關系，所以黏度不恒定，而是隨著冰漿剪切速率的變化而變化。V. Ayel等[5]綜述了Bingham、Casson、Ostwald和Herschel- Buckley四種流變模型，用于預測不同流動狀態下非牛頓流體剪切速率和剪切應力之間的非線性關系。本文采用Ostwald模型，即冪律模型來描述冰漿的流動特性：

(6)

式中：τw為剪切應力，MPa；γw為剪切速率，s-1；n′為流動特征指數(無量綱數)，n′=1時為牛頓流體，n′>1時流體剪切增稠，n′<1時流體剪切變稀；K′為稠度系數，mm2/s，依賴于n′和運動黏度[13]。

在水平直管中，流體在管內作層流流動時，其壁面剪切應力τw和剪切速率γw[14]分別如式(7)和式(8)所示：

τw=DiΔp/(4L)

(7)

γw=8u/Di

(8)

流動特征指數n′和稠度系數K′滿足本構方程，如式(9)所示：

ln(DiΔp/(4L))=lnK′+n′ln(8u/Di)

(9)

因此通過測量流動時的壓降和流速，分別對ln(DiΔp/(4L))和ln(8u/Di)兩個數據進行微分擬合即可得到流體的特征指數n′，即式(10)：

(10)

圖9所示分別為管徑為4、6、8 mm的剪切應力與剪切速率的對數關系。從實驗數據的擬合曲線可以明顯看出，每條線的斜率是恒定的，說明冪律模型適用于冰漿。由式(10)可知，擬合曲線的斜率即為特征指數n′，截距即為稠度系數K′，分析了IPF分別為0.06、0.10、0.14、0.18、0.22和0.26的不同擬合結果，發現隨著IPF的增加，擬合直線的截距逐漸增大，斜率逐漸減小，反映在本構方程(式(9))中，即K′逐漸增大，n′逐漸減小。同時， IPF一定時，隨著管徑的增加，擬合直線的截距逐漸增大，斜率逐漸減小，即K′隨IPF的增大而增大，n′隨IPF的增大而減小。圖10所示為不同管徑的流動特征指數n′、稠度系數K′與IPF的關系。由圖10可知，n′隨著 IPF的增大而減小，當含冰率為0.06時，冰漿中主要為尿素溶液，此時流動特征指數接近于1，隨著IPF的增大，流動特征指數愈發小于1，且隨著IPF的增加更加偏離于1，說明低IPF下冰漿近似為牛頓流體，IPF越高，偏離牛頓流體的程度越大，為非牛頓流體，即擬塑性流體。而K′與IPF則呈正相關，IPF為0.06時，K′接近于0，隨著IPF的增加，稠度系數逐漸升高，這是因為冰漿中固相顆粒的增加使得流動的阻力變大，導致流體的稠度變大。

圖9 剪切應力與剪切速率間的關系

圖10 n′、K′和IPF之間的關系

分析不同管徑對特征指數n′和稠度系數K′的影響，發現在相同IPF下，隨著管徑的增大，n′更偏離于1，即管徑越大，流體偏離于牛頓流體的程度更大。隨著管徑的增大，K′也更加偏離于0，基本呈現規律是管徑越小，流體的流動粘稠度越低，小管徑內偏離牛頓流體的程度越小，這與小管徑內流速較大，溶液和冰粒混合更加均勻，流體更加接近單相流有關。分析IPF為0.26的狀態下的稠度系數，發現4 mm管徑下的稠度系數K′偏離0的程度甚至超過了6 mm管徑下的K′。這可能是因為IPF較高的情況下，雖然高流速下溶液和冰晶混合較均勻，但冰漿內冰晶的體積分數較大，且管徑較小，固相冰粒對管徑的碰撞作用相比于高流速對兩相均勻混合的影響更為顯著，所以在高含冰率(IPF=0.26)下表現出獨特的趨勢。

基于已有的冰漿的流動壓降數據，利用式(11)可以得到水平細管內冰漿流動的范寧摩擦阻力系數：

fF=λ/4=DiΔp/(2Lρu2)

(11)

通過分析摩擦阻力系數比可知，無法將冰漿作為牛頓流體分析其流動特性，因此，本文將冰漿看做一種固液兩相的非牛頓流體進行分析，其流動比單相流更加復雜，單相流體的摩擦阻力來自于黏性摩擦阻力，即流體在壁面附近的黏滯力；而對于包含固相顆粒的兩相流體而言，除了黏性摩擦阻力，還存在由固相顆粒同管壁間碰撞產生的運動摩擦阻力，對于額外的運動摩擦阻力，在層流流動狀態下的影響還較弱，但當流動狀態過渡到湍流狀態時，固相顆粒在徑向的流動也會增強，使固相顆粒同管壁間的相互作用增加，從而使運動摩擦阻力的影響增強[15]。為了更好地描述非牛頓冪律模型流體的復雜的流動特性，A. B. Metzner等[16]提出了修正雷諾數ReMR代替常規雷諾數Re，ReMR可以利用式(12)求得：

(12)

對于層流流動區域，因層流狀態下由冰粒和管壁之間的碰撞引起的運動摩擦阻力影響較小，仍采用牛頓流體的范寧摩擦阻力系數表達式，但在式中利用修正雷諾數ReMR代替常規雷諾數Re。如式(13)所示：

fF=16/ReMR

(13)

在湍流區域，由于冰粒與管壁之間的碰撞而引起的運動摩擦阻力較大而無法忽略，D. W. Dodge等[17]提出了采用充分發展階段光滑管內的半經驗關聯式計算摩擦阻力系數，如式(14)所示：

(14)

圖11所示為修正雷諾數與常規雷諾數的比值與IPF的關系。由圖11可知，比值均小于1，說明修正雷諾數均小于常規雷諾數，這主要是由于冰漿流動時固相黏滯力較大，而且隨著IPF的增大，雷諾數的比值更加偏離于1，修正雷諾數更小。尤其在弱湍流區域，如Re=3 500時，ReMR在三種管徑中均小于2 300，即越過了過渡區，實際狀態已經為層流。該現象可以用層流化[18]的理論來解釋，即在一定的流速下，固相質量分數的增加使固-液兩相的相互作用增強，更多能量由紊流狀態傳送到整個流動截面，導致漿體的整體紊流程度下降，弱紊流區域的流動甚至會從紊流流動狀態轉化為層流狀態。

圖11 雷諾數比值與IPF的關系

雖然雷諾數比值與IPF的變化呈負相關，但隨著常規雷諾數的增加，比值減小的程度逐漸減小，尤其在常規雷諾數Re=8 500時，IPF=0.18、0.22、0.26時的比值只略有減小，基本不變。而在低流速(Re=1 500、3 500)下，雷諾數比值減小的趨勢在IPF為0.3以內一直明顯，由前文分析可知，這兩種常規雷諾數的實際狀態其實均已經為層流狀態。說明冰漿在強湍流區域流動狀態更穩定，而層流和過渡區域的流動特性更為復雜，這主要是因高流速下的冰漿內液相和固相分布更均勻，而低流速下冰粒和溶液的分層現象的出現使流動規律更復雜。

圖12 范寧摩擦阻力系數與修正雷諾數間的關系

圖12所示為不同IPF的范寧摩擦阻力系數與修正雷諾數之間的關系。圖中實線為實驗數據的擬合曲線。當ReMR<2 300時，實驗值按照反比例函數擬合曲線，擬合公式遵循式(13)；當ReMR>2 300時，實驗值按照線性函數擬合曲線，擬合公式遵循式(14)。由圖12可知，隨著ReMR的增加，fF逐漸減小，這和摩擦阻力系數與常規雷諾數的關系是一致的。觀察擬合曲線可知，低雷諾數的曲線擬合和高雷諾數的曲線擬合在不同的IPF下有不同的交點。B. Niezgoda-zelasko等[19]將這種交點對應的修正雷諾數定義為冰漿從層流向紊流轉化的轉折修正雷諾數。M. Grozdek等[20]提出冰漿的轉折雷諾數范圍大多為1 700～2 500，而具體的轉折雷諾數的數值同冰漿的IPF和管徑有關。在本實驗中，隨著IPF的增加，轉折雷諾數的數值隨之增大，即會出現處在某一相同的轉折雷諾數時，低IPF的冰漿流動處于湍流狀態，而高IPF的冰漿流動仍處在過渡區或層流狀態。這就是固液兩相流體中經常出現的層流化現象，即IPF的增加使固液兩相間的相互作用增強，更多能量由湍流核心傳送到整個流動截面，導致漿體的整體湍流程度下降，從而弱湍流區域的流動狀態其實已經轉化為了層流狀態。分析不同管徑的結果，隨著IPF的增加，4 mm管道內的轉折雷諾數的范圍大約在2 500～3 200；6 mm管道內的轉折雷諾數的范圍約為1 600～2 500；8 mm管道內的轉折雷諾數的范圍約為1 500～1 900。所以冰漿的轉折雷諾數隨著管徑的增加而呈現減小的趨勢。

由圖12還可知，當實驗管徑為4 mm時，實驗值略大于理論分析值，主要是由于較小的管徑導致冰漿顆粒與管壁之間的碰撞和摩擦更為嚴重。當管徑為8 mm且ReMR在2 000～3 500之間時，實驗值大于理論分析值。因為此時冰漿流處于層流到湍流的過渡區域，并且充分發展的程度很低，而計算模擬值的半經驗公式是通過大量充分發展的流動狀態數據獲取的，因此此處的誤差較大。

3 結論

本文以質量分數為5%尿素溶液所制取的冰漿為研究對象，通過測量冰漿的運動黏度及冰粒尺寸，觀測冰漿中冰粒分布情況，研究冰漿制備及不同含冰率、不同雷諾數和不同管道內的流動特性，得到如下結論：

1)冰漿中冰粒的尺寸多分布在0.1～0.25 mm之間，平均尺寸約為0.15 mm，冰粒的平均尺寸與IPF之間無明顯關聯。

2)冰漿的壓降隨IPF或Re的增加而增大。隨管徑的增大而減小。冰漿流動壓降對參數的靈敏等級為：管徑>IPF>Re。對于高含冰率的冰漿，湍流狀態下壓降增長速率較大。流態決定了冰漿的基本流動特性。在層流和湍流條件下，IPF對高雷諾數狀態下的壓降影響更為顯著。

3)分析流動摩擦阻力系數的實驗值與利用牛頓流體的經驗公式求得的理論摩擦阻力系數的比值，發現偏離1的程度整體隨著雷諾數的增加而降低，在層流區域的比值較大，且遠大于湍流狀態。在小管徑內，由于流速較大，冰漿中固液兩相混合更加均勻，從而摩擦阻力系數的偏離理論值越小，即更接近于單相牛頓流體。

4)經過分析摩擦阻力系數比，發現無法將冰漿看作牛頓流體，將其看做非牛頓流體，利用冪律模型分析其流動特性。結果顯示，冰漿的流動隨著IPF的增加而愈發偏離牛頓流體。流動特征指數n′隨著IPF的增大而降低，IPF的增加使冰漿更加偏離牛頓流體。而稠度系數K′和IPF的增長呈正相關，當高含冰率(IPF>20%)時，K′的變化更加顯著。此外，當含冰率介于5%～30%范圍內，n′隨著管徑的增大略微減小，K′隨著管徑的增大而逐漸增大。

5)通過分析修正雷諾數和范寧摩擦阻力系數的關系可知，隨著IPF的增加，4 mm、6 mm和8 mm管道內的轉折修正雷諾數的范圍分別為2 500～3 200、1 600～2 300和1 500～1 900。所以冰漿的轉折修正雷諾數隨著管徑的增加而呈現減小的趨勢。