二氧化硅納米顆粒對冰漿中冰晶粒徑分布及存儲演化特性的影響

劉曦，林淑嫻，李歲，李學來

?

二氧化硅納米顆粒對冰漿中冰晶粒徑分布及存儲演化特性的影響

劉曦，林淑嫻，李歲，李學來

（福州大學石油化工學院，福建福州350116）

分別以乙二醇水溶液和氯化鈉水溶液為基液配制不同濃度的二氧化硅納米流體并以此制備冰漿，通過顯微裝置獲得冰晶圖像，將實驗得到的粒徑分布與正態分布、對數正態分布、Gamma分布和Weibull分布進行對比，探討納米二氧化硅對冰晶平均粒徑與分布特性的影響，同時觀測儲存過程中冰晶粒徑演化規律。結果表明：加入納米二氧化硅前后冰晶粒徑分布均可用Gamma分布描述；納米二氧化硅可起到細化晶粒的作用，而且添加濃度越高冰晶顆粒越小；當基液為乙二醇水溶液時，加入納米二氧化硅可較好地抑制儲存過程中的冰晶粒徑增長，但基液為氯化鈉水溶液時，納米二氧化硅濃度需達到0.75%，才可抑制冰晶增大。研究結果證明，一定濃度的納米二氧化硅流體可作為制冰溶液，起到減小冰晶粒徑并控制冰晶生長的作用，這對冰漿流動和傳熱性能的改善具有重要的應用價值。

冰漿；粒度分布；粒徑演化；納米粒子；二氧化硅

引 言

冰漿是指含有冰晶粒子的固液兩相混合物，其中冰晶顆粒的直徑不超過1 mm，具有傳熱效率高、釋冷速度快等特點。作為一種典型的相變蓄冷材料，冰漿已被應用于建筑物蓄冷、食品的加工與儲藏、礦井降溫、人造雪領域，潛在的應用領域包括管道清潔、醫療救助及消防滅火等[1-8]。冰漿中冰晶顆粒的形狀和粒徑是決定其流動性能、換熱性能及物性參數的關鍵[9-13]，同時也會影響系統設備（如泵、閥門）的使用[14]，過大的冰晶顆粒還可能導致冰漿流動過程管道堵塞并最終使系統發生故障。為了獲得小顆粒的冰晶，常常需要在制冰溶液中加入添加劑，由此控制冰晶的生長并改善冰漿的性能，有關添加劑對冰晶形態及粒徑分布的影響規律已有較多文獻可供參考。

Delahaye等[15]通過顯微觀察裝置研究了冰漿及水合物漿體中的晶體粒徑分布，并將實驗數據與正態分布、對數正態分布、Gamma分布及Weibull分布進行對比，研究表明冰漿和四氫呋喃水合物漿的晶體粒徑分布都基本符合Gamma分布和Weibull分布。劉志強等[16-18]建立了可用于反映實際儲存過程中冰晶演化規律的群體平衡模型，探討添加劑種類、濃度、耗損率、含冰率等因素對冰晶粒徑大小及分布的影響，研究認為添加劑種類及濃度是影響冰晶粒徑增大的主要因素，含冰率及耗損率的影響則較小，加入適量添加劑可抑制冰晶增大，而且當含冰率越低、耗損率越高時冰晶粒徑增長越慢。Peng等[19]通過液-液循環流化裝置制取冰晶，借助高分辨率的數碼攝像儀獲取冰晶圖像，分析了不同操作參數對冰晶粒徑分布及粒徑大小的影響。Pronk等[20-21]建立了冰漿懸浮液的Ostwald熟化動力學模型，用于模擬冰漿儲存過程中冰晶粒徑分布的變化，并研究了冰漿的儲存特性，研究表明，受Ostwald熟化影響，冰漿在恒溫儲存過程中冰晶粒徑顯著增大，而且生長速率隨溶液濃度升高而減緩。Inada等[22]研究了聚乙烯醇對恒溫冰漿中冰晶粒徑增大的抑制作用，并與聚乙二醇、聚乙烯基吡咯烷酮、聚丙烯酸、氯化鈉及Ⅰ類抗凍蛋白進行對比，實驗結果表明聚乙烯醇與Ⅰ類抗凍蛋白作用相似，低濃度下就可完全抑制冰晶粒徑的增大。

近年來，隨著納米技術的迅速發展，利用納米流體制備冰漿及其他相變蓄冷材料已引起國內外研究人員的廣泛興趣。納米流體是將納米顆粒懸浮于液體中形成的一種固液兩相混合物。通常添加納米顆粒可促進冰晶生長過程的非均質成核，有效降低溶液的過冷度，這對降低結晶所需能耗、提高系統COP具有重要意義[23-24]。同時，在水溶液中加入納米顆粒可提高溶液的熱導率并強化傳熱過程[25-28]。理論上，納米粒子的添加會促進異質成核，使水溶液中的晶核數增加，這必將影響冰晶的形態、粒徑及分布特性，故有必要對此進行深入研究。

本工作采用最常見的兩種制冰溶液乙二醇水溶液和氯化鈉水溶液為基液，分別加入納米二氧化硅并制備得到冰漿，通過顯微觀測系統及圖片處理技術獲得一系列冰晶粒徑數據，將實驗值與典型的4種概率分布即正態分布、對數正態分布、Gamma分布和Weibull分布進行對比，從而確定添加納米顆粒后冰晶的粒徑分布特點，同時研究了不同儲存時間下的粒徑變化規律，研究成果可為冰漿應用過程中冰晶粒徑的控制提供參考。

1 實驗部分

1.1 實驗原料

乙二醇，廣東光華科技股份有限公司，分析純；氯化鈉，國藥集團，分析純，含量≥99.5%；納米二氧化硅，深圳晶材化工有限公司，該產品采用溶劑熱法制得，納米顆粒的平均粒徑為30 nm，其水分散液具有良好的穩定性和分散性，而且納米顆粒不易團聚；去離子水，實驗室自制。

1.2 冰漿制備系統

冰漿制備系統如圖1所示，該系統主要由精密恒溫液浴槽、電動攪拌器、不銹鋼容器、溫度采集器、T型熱電偶等組成。液浴槽型號為RH-DC-V21-R20，其溫度可在-20～100℃之間無級調節，精度為0.05℃。溫度采集采用USB1608溫度記錄儀，誤差為±0.1℃，記錄儀的記錄間隔為1 s。攪拌裝置采用DJ型精密電動攪拌器，轉速在0～3000 r·min-1內可調。制冰漿容器直徑80 mm，高度60 mm，容積為300 ml，材質為不銹鋼。實驗過程中首先開啟液浴槽，待溫度降至-10℃，放入制冰漿容器，用支架固定，同時啟動攪拌器并開啟溫度采集系統，密切觀察溫度變化規律，當溶液結冰時可觀察到溫度會有明顯的躍變，待結冰時間達到120 s后停止實驗，取出冰漿。

1.3 冰晶粒徑觀測系統

常用的冰晶粒徑觀測方法包括聚焦光束反射測量法和顯微觀測法。聚焦光束反射測量法是一種較為先進的在線顆粒粒度測定法，其特點是無須取樣，可直接采用聚焦光束反射測量儀（FBRM）實時在線監測透明或不透明的懸浮液及高速流動懸浮液中的顆粒粒徑。目前已有研究人員利用該設備觀測水合物生成過程的顆粒變化[29]。顯微觀測則是一種較為傳統的冰晶粒徑觀測方法，雖無法實現實時觀測，但系統組成較為簡單且穩定，文獻[15，21-22]中均采用該方法對冰晶粒徑進行觀測，效果良好，故在本研究中采用顯微觀測法。

冰晶粒徑觀測系統由顯微鏡、低溫冷臺、工業相機、表面皿等組成，如圖2所示。顯微鏡采用SZM45體視顯微鏡，低溫冷臺采用KER5100-09S正負溫精密恒溫儀，工業相機型號為scA1390-17gc。通過顯微鏡對制得的冰漿進行觀察，觀察過程中采用低溫冷臺使表面皿底部保持恒溫，避免冰晶融化。利用工業相機拍攝若干照片，以保證至少獲取400個冰晶顆粒樣本。對每個樣品，圖像拍攝需在3 min內完成，因此可忽略拍攝過程中因Ostwald熟化、團聚、融化等動力學行為對冰晶粒徑的影響。

2 實驗數據處理

2.1 冰晶顆粒圖片處理

實驗所拍攝到的冰晶顆粒如圖3所示，從圖中可看出冰晶顆粒有著不同的輪廓，通常呈現橢圓形或近似圓形。由于所觀察的顯微照片中相鄰冰晶具有重疊性，冰晶的投影面積需在圖像分析軟件ImageJ中通過手工描繪冰晶的輪廓確定。

將冰晶的當量圓直徑定義為晶粒的特征尺寸Feret,i，可通過冰晶顆粒的投影面積計算得到[15,30][式（1）]；冰晶顆粒的平均粒徑則根據每個顆粒的直徑及顆粒總數獲得[式（2）]。

(2)

2.2 冰晶粒徑分布概率函數

常用的顆粒粒度概率分布函數包括正態分布、對數正態分布、Gamma分布和Weibull分布。本工作擬采用這4種經驗分布函數與實驗值進行對比，從而確定適用于描述冰晶顆粒粒徑分布的函數。

正態分布概率密度函數

式中，為冰晶顆粒粒徑的平均值，為冰晶顆粒粒徑的標準差。

對數正態分布概率密度函數

式中，為冰晶顆粒粒徑的對數平均值，為冰晶顆粒粒徑的對數標準差。

Gamma分布概率密度函數

式中，等于冰晶顆粒粒徑的平均值，2等于冰晶顆粒粒徑的方差，為Gamma函數。

圖4 納米二氧化硅濃度對乙二醇水溶液制得冰漿中冰晶粒徑分布的影響

Fig.4 Effects of nanosilica concentration on probability distributions for ice slurry made of aqueous ethylene glycol solutions

Weibull分布概率密度函數

，

式中，是形狀參數，是尺度參數，和分別為冰晶顆粒粒徑的平均值和標準方差。

為了比較實驗值和概率密度函數理論計算值之間的誤差，建立誤差評判指標如下

式中，為粒徑分布的區間分組數，PDFexp,為該分布區間的分布概率，PDFtheory,i為各經驗概率密度函數在對應分布區間內的分布概率理論值。

3 實驗結果與討論

3.1 納米二氧化硅對冰晶粒徑分布的影響

分別在150 ml 3%乙二醇水溶液和3%氯化鈉水溶液中添加納米二氧化硅，配備質量分數分別為0.25%、0.50%、0.75%的納米流體，并制得冰漿。觀察得到冰漿中的粒徑分布，并將冰晶粒徑概率分布實驗值與根據實驗值擬合的正態分布、對數正態分布、Gamma分布、Weibull分布曲線進行對比，如圖4和圖5所示。由式（7）計算得到冰晶粒徑分布實驗值與4種理論分布值間的誤差，見表1。

由圖4和表1可知：不添加納米二氧化硅時，乙二醇水溶液制得的冰漿中冰晶粒徑分布較符合正態分布和Gamma分布，冰晶粒徑分布在10～80 μm間；當納米二氧化硅的添加量為0.25%時，冰晶粒徑分布較符合對數正態分布和Gamma分布；當納米二氧化硅的添加量為0.50%時，冰晶粒徑分布較符合Gamma分布和對數正態分布；當納米二氧化硅的添加量為0.75%時，冰晶粒徑分布較符合對數正態分布和Gamma分布。總體來說，加入納米二氧化硅后，冰晶粒徑分布范圍基本不變，分布概率較符合對數正態分布和Gamma分布。

由圖5和表1可知：不添加納米二氧化硅時，氯化鈉水溶液制得的冰漿中冰晶粒徑分布較符合對數正態分布和Gamma分布，冰晶粒徑分布在20～120 μm間；加入0.25%、0.50%、0.75%的納米二氧化硅后，冰晶粒徑分布范圍分別為20～90 μm、10～90 μm、20～70 μm，均有不同程度的縮小。當納米二氧化硅的添加量為0.25%時，冰晶粒徑分布較符合Weibull分布和Gamma分布；當納米二氧化硅的添加量為0.50%時，冰晶粒徑分布較符合Gamma分布和對數正態分布；當納米二氧化硅的添加量為0.75%時，冰晶粒徑分布較符合Gamma分布和對數正態分布。

表1 冰晶粒徑分布實驗值與4種理論分布值間的誤差

基于以上分析可發現：當基液為乙二醇水溶液或氯化鈉水溶液時，加入不同劑量的納米二氧化硅后，兩種納米流體制得冰漿中的冰晶粒徑分布規律基本不變，均可用Gamma分布描述。

圖6所示為納米二氧化硅對乙二醇水溶液和氯化鈉水溶液中冰晶平均粒徑的影響。從圖中可看出：相比于不添加納米二氧化硅時制得的冰晶平均粒徑，在乙二醇溶液中加入0.25%、0.50%、0.75%納米二氧化硅時，平均粒徑從44.27 μm分別降至42.19、40.31、38.45 μm，降幅分別達4.70%、8.95%、13.15%；在氯化鈉溶液中加入0.25%、0.50%、0.75%納米二氧化硅時，平均粒徑從62.29 μm分別降至51.02、40.33、34.15 μm，降幅分別達18.09%、35.25%、45.17%。理論上，納米二氧化硅作為一種典型的成核劑，具有比表面積大、接觸角小、納米顆粒與水溶液間潤濕性好等特點，因此在乙二醇水溶液或氯化鈉水溶液中添加納米二氧化硅均可有效降低成核時的表面能壁壘，促進冰晶的異質成核。一般來說，成核劑可使晶粒尺寸細微化，溶液中含有的晶核數量越多形成的冰晶顆粒尺寸越小，故冰晶平均粒徑隨納米二氧化硅含量增大而減小。

觀察圖6還可發現，當基液為氯化鈉水溶液時納米二氧化硅的細化晶粒效果更明顯。分析其原因，主要是由于：乙二醇水溶液中，醇類物質自身所含的羥基可破壞水分子與水分子間形成的氫鍵，阻礙冰晶的生長與聚集，即乙二醇可起到一定的抑制冰晶生長作用，加入納米二氧化硅后冰晶粒徑進一步降低，但由于未添加成核劑時粒徑已較小，降幅不明顯；氯化鈉水溶液中雖含有一定量的氯離子，但氯離子半徑大、電負性低，形成的氫鍵較弱，故氯化鈉本身不能很好地抑制冰晶生長，加入納米二氧化硅后晶核數量顯著增加，使得冰晶粒徑明顯減小。

3.2 納米二氧化硅對冰晶粒徑演化規律的影響

冰漿在儲存過程中，冰晶易發生破碎、團聚、Ostwald熟化等動力學行為，這些動力學行為會使冰晶的粒徑發生變化，從而影響冰漿的流動和傳熱過程，因此有必要對冰晶粒徑的演化特性進行研究。

圖7給出了不同儲存時間下3%乙二醇水溶液制得冰漿中冰晶粒徑的演化曲線和添加0.25%納米二氧化硅后冰晶粒徑的演化曲線。圖8給出了不同儲存時間下3%氯化鈉水溶液制得冰漿中冰晶粒徑的演化曲線和添加0.25%納米二氧化硅后冰晶粒徑的演化曲線。從圖中可看出，添加納米前后冰漿中冰晶粒徑的演化規律大體上一致，而且粒徑分布基本與Gamma分布相吻合。隨著時間的推移，冰晶粒徑分布曲線向右移動，小直徑晶粒所占的比重減小，同時大直徑晶粒所占的比重增大，這表明冰晶在儲存過程中存在明顯的團聚和Ostwald熟化現象，團聚使得多個小直徑冰晶匯聚成大直徑冰晶，Ostwald熟化作用則使小直徑冰晶不斷消融而大直徑冰晶不斷增大，在這兩種動力學行為的聯合作用下粒度分布曲線不斷右移，小冰晶顆粒數減小而大冰晶顆粒數增大。此外，實驗所用的冰漿在儲存過程中未加攪拌，大直徑晶粒破碎成小直徑晶粒的概率極低，從而造成儲存過程中冰晶的團聚作用遠大于破碎作用，因此粒徑分布密度的峰值不斷減小，同時分布曲線的坡度逐漸變緩。

圖9所示為在乙二醇水溶液和氯化鈉水溶液中分別加入不同濃度的納米二氧化硅后所形成的冰漿在儲存過程中冰晶平均粒徑隨時間的演化曲線。觀察每條粒徑增長曲線可知：冰晶粒徑的增長速率基本上隨時間推移而降低。這是因為隨著時間推移冰晶顆粒總數不斷減少，團聚和Ostwald熟化現象的作用效果不斷削弱，因此粒徑的增長越來越不明顯。當納米二氧化硅質量分數為0、0.25%、0.50%、0.75%時，乙二醇水溶液中冰晶的初始平均粒徑分別為44.27、42.19、40.31、38.45 μm，儲存2 h后粒徑分別增至92.47、85.80、76.39、73.46 μm，增長率分別為108.88%、103.37%、89.51%、91.05%；氯化鈉水溶液中冰晶的初始平均粒徑分別為62.29、51.02、40.33、34.15 μm，儲存2 h后粒徑分別增至134.69、116.83、90.88、68.55 μm，增長率分別為116.23%、128.99%、125.34%、100.73%。以上數據表明，在乙二醇水溶液中加入納米二氧化硅可降低冰晶增長速率，抑制效果隨添加濃度而變，當添加量為0.50%時效果最明顯。在氯化鈉水溶液中加入少量的納米二氧化硅時，冰晶粒徑增長率有所上升，表明過低的添加濃度起不到抑制作用，只有當添加量達到0.75%時才可抑制冰晶生長。總體來說，當添加等量的納米二氧化硅時，乙二醇水溶液中冰晶平均粒徑的增長速率明顯低于氯化鈉中的增長速率，由此說明在乙二醇水溶液中添加納米二氧化硅抑制冰晶顆粒生長的效果優于氯化鈉水溶液。

4 結 論

（1）乙二醇水溶液和氯化鈉水溶液制得的冰漿中冰晶粒徑分布均可用Gamma分布描述。加入納米二氧化硅后，粒徑分布范圍略微變窄，但粒徑分布仍符合Gamma分布。因此，可統一采用Gamma分布描述添加納米二氧化硅前后乙二醇水溶液和氯化鈉水溶液中冰晶粒徑的分布規律。

（2）納米二氧化硅可起到細化冰晶顆粒的作用，在乙二醇水溶液和氯化鈉水溶液中加入納米二氧化硅均可獲得更小粒徑的冰晶，而且添加濃度越高冰晶顆粒越小，冰晶平均粒徑與納米二氧化硅濃度呈現良好的線性關系。

（3）受團聚和Ostwald熟化作用的影響，冰漿儲存過程中冰晶粒徑顯著增大，隨著時間的推移增大速率逐漸減緩。在乙二醇水溶液中加入納米二氧化硅可在一定程度上抑制顆粒粒徑的增長，而且當添加等量的納米二氧化硅時乙二醇水溶液中冰晶平均粒徑的增長速率均明顯低于氯化鈉中的增長速率，因此含納米二氧化硅的乙二醇水溶液較適用于冰漿制備，有望獲得具有良好的流動性能和傳熱性能的冰漿。

符 號 說 明

Ai——冰晶顆粒投影面積，m2 D——冰晶顆粒平均粒徑，m DFeret,i——第i個冰晶顆粒特征尺寸，m N——冰晶總顆粒數 x——冰晶顆粒粒徑變量，m a——形狀參數 b——尺度參數，m m——平均值，m s——標準差，m2

References

[1] SAITO A. Recent advances in research on cold thermal energy storage[J]. International Journal of Refrigeration, 2002, 25(2): 177-189.

[2] BELLAS I, TASSOU S A. Present and future applications of ice slurries[J]. International Journal of Refrigeration, 2005, 28(1): 115-121.

[3] FANG G Y, TANG F, CAO L. Dynamic characteristics of cool thermal energy storage systems—a review[J]. International Journal of Green Energy, 2016, 13(1): 1-13.

[4] EGOLF P W, KAUFFELD M. From physical properties of ice slurries to industrial ice slurry applications[J]. International Journal of Refrigeration, 2005, 28(1): 4-12.

[5] LI G, HWANG Y, RADERMACHER R. Review of cold storage materials for air conditioning application[J]. International Journal of Refrigeration, 2012, 35(8): 2053-2077.

[6] KAUFFELD M, WANG M J, GOLDSTEIN V,. Ice slurry applications[J]. International Journal of Refrigeration, 2010, 33(8): 1491-1505.

[7] WANG M J, KUSUMOTO N. Ice based storage in multifunctional buildings[J]. Heat Mass Transfer, 2001, 37(6): 594-604.

[8] QUARINI J. Ice-pigging to reduce and remove fouling and to achieve clean-in-place[J]. Applied Thermal Engineering, 2002, 22(7):747-753.

[9] YOUSSEF Z, DELAHAYE A, HUANG L,. State of the art on phase change material slurries[J]. Energy Conversion and Management, 2013, 65(1): 120-132.

[10] AYEL V, LOTTIN O, PEERHOSSAINI H. Rheology, flow behavior and heat transfer of ice slurries: a review of the state of the art[J]. International Journal of Refrigeration, 2003, 26(1): 51-59.

[11] DORON P, GRANICA D, BARNEA D. Slurry flow in horizontal pipes-experimental and modeling[J].Multiphase Flow, 1987, 13(4): 535-547.

[12] DORON P, BARNEA D. A three-layer model for solid-liquid flow in horizontal pipes[J]. Multiphase Flow, 1993, 19(6): 1029-1043.

[13] EGOLF P W, KITANOVSKI A, ATA-CAESAR D,. Thermodynamics and heat transfer of ice slurries[J]. International Journal of Refrigeration, 2005, 28(1): 51-59.

[14] MIKA L. Ice slurry flow in a poppet-type flow control valve[J]. Experimental Thermal and Fluid Science, 2013, 45(1): 128-135.

[15] DELAHAYE A, FOURNAISON L, GUILPART J. Characterisation of ice and THF hydrate slurry crystal size distribution by microscopic observation method[J]. International Journal of Refrigeration, 2010, 33(8): 1639-1647.

[16] 劉志強, 王肖肖, 王小倩, 等. 冰漿存儲過程中冰晶粒徑演化的影響因素研究[J]. 熱科學與技術, 2013, 12(4): 307-312. LIU Z Q, WANG X X, WANG X Q,. Analysis of influence factors of ice slurry in storage[J]. Journal of Thermal Science and Technology, 2013, 12(4): 307-312.

[17] 趙騰磊, 劉志強, 徐愛祥, 等. 冰漿存儲過程中冰晶粒徑演化數值模擬[J]. 中南大學學報, 2014, 45(10): 3651-3656. ZHAO T L, LIU Z Q, XU A X,. Numerical simulation of evolution of ice crystal size distribution during storage[J]. Journal of Central South University, 2014, 45(10): 3651-3656.

[18] 徐愛祥, 劉志強, 趙騰磊, 等. 冰漿存儲過程中冰晶粒徑動力學演化影響因素[J]. 中南大學學報, 2015, 46(8): 3138-3144. XU A X, LIU Z Q, ZHAO T L,. Factors influencing dynamics evolution of ice crystals during ice slurry storage[J]. Journal of Central South University, 2015, 46(8): 3138-3144.

[19] PENG Z B, YUAN Z L, LIANG K F,. Ice slurry formation in a cocurrent liquid-liquid flow[J]. Chinese Journal of Chemical Engineering, 2008, 16(4): 552-557.

[20] PRONK P, INFANTE FERREIRA C A, WITKAMP G J. A dynamic model of Ostwald ripening in ice suspensions[J]. Journal of Crystal Growth, 2005, 275(1/2): 1355-1361.

[21] PRONK P, HANSEN T M, INFANTE FERREIRA C A,. Time-dependent behavior of different ice slurries during storage[J]. International Journal of Refrigeration, 2005, 28(1): 27-36.

[22] INADA T, MODAK P R. Growth control of ice crystals by poly(vinyl alcohol) and antifreeze protein in ice slurries[J].Chemical Engineering Science, 2006,61(10): 3149-3158.

[23] YOSHI K T, AKIO S, SEIJI O. A study of supercooling-phenomenon and freezing probability of water inside horizontal cylinders[J]. International Journal of Refrigeration, 2004, 27(3): 242-247.

[24] FARID M M, KHUDHAIR A M, RAZACK S A K,. A review on phase change energy storage: materials and applications[J]. Energy Conversion and Management, 2004, 45(9): 1597-1615.

[25] HE Q B, WANG S F, TONG M W,. Experimental study on thermophysical properties of nanofluids as phase-change material (PCM) in low temperature cool storage[J]. Energy Conversion and Management, 2012, 64: 199-205.

[26] LEE J H, LEE S H, CHOI C J,. A review of thermal conductivity data, mechanisms and models for nanofluids[J]. International Journal of Micro-Nano Scale Transport, 2010, 1(4): 269-322.

[27] XIE H Q. Thermal conductivity enhancements of suspensions containing nanosized alumina particles [J]. Journal of Applied Physics, 2002, 91(7):4568-4572.

[28] WANG X, XU X F, CHOI S U S. Thermal conductivity of nanoparticle-fluid mixture[J]. Journal of Thermophysics and Heat Transfer, 1999, 13(4): 474-480.

[29] CLAIN P, NDOYE F T, DELAHAYE A,. Particle size distribution of TBPB hydrates by focused beam reflectance measurement (FBRM) for secondary refrigeration application[J]. International Journal of Refrigeration, 2015, 50: 19-31.

[30] 劉海紅, 李玉星, 王武昌, 等. 四氫呋喃水合物和一氟二氯乙烷水合物顆粒聚結特性[J]. 化工學報, 2014, 65(6):2049-2055. LIU H H, LI Y X, WANG W C,. Agglomeration characterization of THF and HCFC-141b hydrate particle[J]. CIESC Journal, 2014, 65(6): 2049-2055.

Effect of nanosilica on size distribution and evolution of ice crystal particles during storage of ice slurry

LIU Xi, LIN Shuxian, LI Sui, LI Xuelai

(School of Chemical Engineering, Fuzhou University, Fuzhou 350116, Fujian, China)

Ice slurries of nanosilica fluid at different concentrations were prepared by dispersing silica nanoparticles in aqueous solutions of eitherethylene glycol or sodium chloride.Images of ice crystal particles in ice slurries were obtained bya microscopic observation system. The resultingsize distributions of ice crystals were compared to normal, log-normal, Gamma, and Weibull distributions so as to investigate nanosilicaeffect on average diameter andsize distribution of ice particles andto study dimensionalchangeof ice particles during storageof ice slurries. The experimental results indicatedthatice particle sizeswerefitted well with Gamma distribution before and after addition of nanosilica. Nanosilica playedan important role in grain refinement such that sizes of ice particlesdecreasedlinearly with increasing concentrationsof nanosilica. When ethylene glycol aqueoussolutionwas base fluid, nanosilica addition effectively inhibitedgrowth of ice particlesduring slurry storage. When sodium chloride aqueous solutionwas base fluid, concentration of nanosilica hadto be increased from 0.25% to 0.75% for similar inhibition effect. Therefore, nanosilica couldbe usedto reduce particle size and to control growth of ice crystals, which has an important application forflow and heat transfer improvement of ice slurry.

ice slurry; particle size distribution; particle size evolution; nanoparticles; silica

10.11949/j.issn.0438-1157.20160937

TK 02

A

0438—1157（2017）03—0870—09

國家自然科學基金國家基礎科學人才培養基金項目（J1103303）；福建省中青年教師教育科研項目（JA14053）。

2016-07-06收到初稿，2016-11-24收到修改稿。

聯系人：李學來。第一作者：劉曦（1983—），女，博士研究生，講師。

2016-07-06.

LI Xuelai, lxl6632@sina.com

supported by the National Science Foundation for Fostering Talents in Basic Research of the National Natural Science Foundation of China (J1103303) and the Program of Young Teacher Education and Research of Fujian Province (JA14053).