采用界面跟蹤法FTM（Front-Tracking Method）,并結合CSF（continuum surface force）模型,研究了在垂直方向上溫度分布不均勻的對稱流場中由Marangoni效應引起的氣泡上升運動問題。模擬了在不同的M a數(shù)和Pr數(shù)下單氣泡及同軸雙氣泡的運動。研究結果表明,在不同的M a數(shù)下氣泡的運動差異較大,M a數(shù)越大,氣泡運動至穩(wěn)態(tài)時的速度越大,且氣泡運動的最大速度值與M a數(shù)呈正相關關系;增大Pr數(shù)所造成的粘度增大或熱擴散率減小將削弱氣泡的遷移運動;Marangoni對流中雙氣泡的局部運動證實了溫度梯度和氣泡運動速度緊密相關。 

【文章來源】：計算力學學報. 2019,36(01)北大核心

【文章頁數(shù)】：7 頁

【部分圖文】：

圖１氣泡上升物理模型

下的熱毛細遷移距離較短。圖５相同工況下的溫度梯度分布能很好地解釋這一現(xiàn)象。Ｐｒ數(shù)增大意味著粘度增大或熱擴散系數(shù)減小，此時流體的傳熱能力減小，氣泡上下溫度梯度減小，Ｍａｒａｎｇｏｎｉ效應形成的驅動力減弱，故熱毛細遷移距離變短。圖６給出了Ｐｒ數(shù)分別為０．６４，１和４三種工況下速度隨時間的變化過程。與圖３的結果類似，氣泡的運動是一個連續(xù)平穩(wěn)的過程，沒有發(fā)生突變。氣泡加速至穩(wěn)定狀態(tài)后，其運動速度保持基本不變。比較發(fā)現(xiàn)，Ｐｒ數(shù)越大，氣泡的運動速度越校圖２不同時刻氣泡界面形態(tài)（Ｍａ＝１５３６）Ｆｉｇ．２ＳｈａｐｅｏｆｂｕｂｂｌｅｉｎｔｅｒｆａｃｅａｔｄｉｆｆｅｒｅｎｔｔｉｍｅｗｈｅｎＭａ＝１５３６圖３Ｍａ＝１２８０，１５３８和１７９２時的氣泡速度曲線Ｆｉｇ．３ＢｕｂｂｌｅｖｅｌｏｃｉｔｙｃｕｒｖｅｆｏｒＭａ＝１２８０，１５３８ａｎｄ１７９２圖４平衡速度與Ｍａ數(shù)的散點擬合Ｆｉｇ．４ＳｃａｔｔｅｒｐｌｏｔｓｏｆｅｑｕｉｌｉｂｒｉｕｍｖｅｌｏｃｉｔｙｗｉｔｈＭａｎｕｍｂｅｒａｎｄｉｔ’ｓｆｉｔｔｉｎｇｃｕｒｖｅ４５計算力學學報第３６卷

下的熱毛細遷移距離較短。圖５相同工況下的溫度梯度分布能很好地解釋這一現(xiàn)象。Ｐｒ數(shù)增大意味著粘度增大或熱擴散系數(shù)減小，此時流體的傳熱能力減小，氣泡上下溫度梯度減小，Ｍａｒａｎｇｏｎｉ效應形成的驅動力減弱，故熱毛細遷移距離變短。圖６給出了Ｐｒ數(shù)分別為０．６４，１和４三種工況下速度隨時間的變化過程。與圖３的結果類似，氣泡的運動是一個連續(xù)平穩(wěn)的過程，沒有發(fā)生突變。氣泡加速至穩(wěn)定狀態(tài)后，其運動速度保持基本不變。比較發(fā)現(xiàn)，Ｐｒ數(shù)越大，氣泡的運動速度越校圖２不同時刻氣泡界面形態(tài)（Ｍａ＝１５３６）Ｆｉｇ．２ＳｈａｐｅｏｆｂｕｂｂｌｅｉｎｔｅｒｆａｃｅａｔｄｉｆｆｅｒｅｎｔｔｉｍｅｗｈｅｎＭａ＝１５３６圖３Ｍａ＝１２８０，１５３８和１７９２時的氣泡速度曲線Ｆｉｇ．３ＢｕｂｂｌｅｖｅｌｏｃｉｔｙｃｕｒｖｅｆｏｒＭａ＝１２８０，１５３８ａｎｄ１７９２圖４平衡速度與Ｍａ數(shù)的散點擬合Ｆｉｇ．４ＳｃａｔｔｅｒｐｌｏｔｓｏｆｅｑｕｉｌｉｂｒｉｕｍｖｅｌｏｃｉｔｙｗｉｔｈＭａｎｕｍｂｅｒａｎｄｉｔ’ｓｆｉｔｔｉｎｇｃｕｒｖｅ４５計算力學學報第３６卷

【參考文獻】：
期刊論文
[1]自潤濕液體過冷核態(tài)沸騰射液流現(xiàn)象分析[J]. 周樂平,李媛園,魏龍亭,杜小澤,王補宣.  化工學報. 2014(S1)
[2]多相傳質過程中的Marangoni效應[J]. 汪洋,陳杰,王智慧,楊超,毛在砂.  化工學報. 2013(01)
[3]溫度和活性劑濃度梯度協(xié)同驅動的液滴鋪展特性[J]. 王松嶺,李春曦,葉學民.  化工學報. 2011(09)
[4]高濃度甘氨酸離子液體水溶液的降膜特性研究[J]. 彭璟,耿皎,王明明,張鋒,張志炳.  南京大學學報(自然科學版). 2010(05)
[5]豎直流動皂膜中解吸導致的Marangoni對流觀察[J]. 沙勇,李樟云,江桂仙,吐松,肖宗源,葉李藝.  化工學報. 2010(05)
[6]Marangoni效應對降膜加熱流動的影響[J]. 張鋒,耿皎,馬少玲,張志炳.  化工學報. 2005(09)
[7]伴有Marangoni效應的傳質動力學[J]. 沙勇,成弘,袁希鋼,余國琮.  化工學報. 2003(11)
[8]Marangoni效應與液膜振蕩[J]. 李冬梅,賀占博.  化學進展. 2003(01)



本文編號：3037151