基于FLUENT的LNG儲罐瞬態蒸發率仿真實驗

管 官， 林 焰，b， 楊 蕖， 周 帥

(大連理工大學 a.船舶CAD工程中心； b.工程裝備結構分析國家重點實驗室，遼寧 大連 116024)

·專題研討——虛擬仿真實驗(49)·

基于FLUENT的LNG儲罐瞬態蒸發率仿真實驗

管 官a， 林 焰a，b， 楊 蕖a， 周 帥a

(大連理工大學 a.船舶CAD工程中心； b.工程裝備結構分析國家重點實驗室，遼寧 大連 116024)

為使學生接觸專業科研前沿，更好地培養其科研與創新能力，將LNG儲罐瞬態蒸發率研究內容引入到實驗教學當中，開展了基于FLUENT的LNG儲罐瞬態蒸發率仿真實驗研究。利用FLUENT進行仿真實驗，可直觀呈現LNG儲罐在不同充滿率、不同環境溫度、不同內徑尺寸下LNG蒸發率的變化規律，便于學生理解LNG蒸發規律。通過該教學和科研相結合的仿真實驗，提高了實驗教學效果，有利于提高學生的仿真能力，激發學生的科學探索精神，培養學生運用實驗手段解決科研問題的能力。

液化天然氣儲罐； 瞬態蒸發率； 仿真實驗； FLUENT

0 引 言

科研驅動教學模式是研究型大學教學改革的重要內容之一[1]，具有教學方法和教學內容上的創新性。科學研究是高質量教學的有力保證[2]。在教學改革中，將專業前沿的科研課題融入到教學環節，使學生有機會接觸到實際的工程實踐，有助于提升學生的學習積極性，培養學生的科學素養和實際工程能力以至有所創新[3]。

隨著液化天然氣(Liquefied Natural Gas,LNG)船型的開發和廣泛應用，LNG的安全存儲已經成為LNG船能否安全使用一個關鍵問題。由于LNG易燃易爆的特性，必須保證其不能泄漏，LNG儲罐的絕熱性就顯得至關重要，而蒸發率是評價絕熱性的重要指標，因此，為保證船舶LNG儲罐的安全使用，對LNG蒸發率進行研究，無論是從工程應用還是科研教學來講，都具有著重要的意義[4-7]。

為了加強培養學生的科研能力與創新能力，本文將LNG儲罐瞬態蒸發率研究內容引入到教學實驗中，設計了“基于FLUENT的LNG儲罐瞬態蒸發率仿真實驗”，給學生創建一個接觸科研前沿、應用專業知識的平臺，創建科學研究的情境，激發學生的科學探索熱情。FLUENT作為當今計算仿真技術中最有效、最實用的工具之一，為仿真實驗教學提供了工具[8]。將FLUENT仿真融入LNG儲罐瞬態蒸發率實驗教學中，不僅可以促進學生對蒸發規律的理解，還能增強學生的仿真實踐能力。這種科研與實踐相結合的教學模式能有效提高教學質量，培養學生從中學到新的思維方式和運用實驗手段解決科研問題的能力[9-11]。

1 實驗計算模型

1.1計算模型分析

本實驗針對大連理工大學船舶制造國家工程中心研制的1.5 m3獨立C型LNG儲罐(已授權發明專利)[12]，對其瞬態蒸發過程采用均相表面蒸發率模型進行仿真，將環境漏熱量均用于LNG蒸發，初始假定罐內所有介質溫度都相等，環境漏熱量作為LNG蒸發唯一熱源，計算采用MIXture模型[13]，因此，整個儲罐內的溫度場分為氣相，液態相以及氣液混合相。此外，由于C型罐外形簡單，為圓柱形且邊界條件軸對稱，在求解C型罐內LNG自然蒸發率問題時，可以把求解區域看做C—H—C型區域，考慮減少計算量，本實驗將LNG自然蒸發作為一個二維的熱傳導問題，計算模型如圖1所示，由于C型罐實物尺寸不大，采用實際尺度計算仿真。

圖1 獨立C型LNG儲罐

1.2模型建模及網格

采用FLUENT前處理軟件Gambit建模及網格劃分，由于LNG儲罐模型形狀相對規則，故采用全結構化網格劃分策略，并定義相應邊界條件類型以及網格流域，此外，為了確保計算精確，在漏熱區網格進行加密。網格劃分圖如圖2所示，網格數約為3萬個。

(a) 0.2充滿率

(b) 0.5充滿率

本仿真實驗探討影響LNG蒸發率各種因素，包括LNG儲罐充滿率、LNG儲罐外部環境溫度及LNG儲罐內徑。為能清楚反應各影響因子對儲罐內LNG蒸發率的影響，本實驗采用單參數變化影響因子，考慮充滿率變化0.1～0.9，間隔0.1；內徑0.5～2.5 m，間隔0.5 m；環境漏熱溫度118～122 K，間隔1 K。

1.3模型邊界條件設置

本實驗計算LNG瞬態蒸發率，求解時均連續性方程、動量方程、時均湍流方程以及能量方程，采用基于壓力求解器以及SIMPLE格式，多相流模型使用Mixture模型，并結合標準兩方程湍流方程k-ε計算湍流場。此外在計算LNG蒸發率時考慮相應的模型簡化，即認為液態LNG為均質不可壓縮流體，并忽略流體粘性耗散，且不考慮LNG蒸發時物性參數變化。本實驗結合均相表面蒸發模型，模型邊界均設置為無滑移壁面邊界條件，初始計算C型罐內氣液兩相溫度均為109.15 K，漏熱區壁面溫度高于飽和溫度111.15 K。LNG蒸發過程是由液態相達到飽和溫度時，向氣態進行轉換，存在相變過程，本實驗基于C語言，編寫相變UDF函數，通過FLUENT中的UDF模塊，將相變函數加載到FLUENT求解源相，從而模擬LNG蒸發相變過程。

UDF函數主要分為四段分別為：C語言頭文件及定義相變條件；定義液相轉變條件及算法，即液相單元溫度若達到飽和溫度時，液相向氣相進行質量轉移；定義氣相轉變算法，即定義氣液混合區單元溫度若達到飽和溫度時，液相向氣相進行質量轉移；定義混合模型能量源相，即定義混合區單元溫度達到飽和溫度時，液態LNG相氣相轉換，并吸收熱量[14-15]。

2 實驗結果與分析

2.1充滿率影響

對1.5 m3LNG儲罐在不同充滿率下的蒸發進行仿真，得到從0.1～0.9充滿率下蒸發量隨時間變化的曲線，漏熱溫度為120 K，內徑為0.5 m，以0.5充滿率為例，圖3給出了0.5充滿率蒸發量變化曲線。

圖3 0.5充滿率蒸發量變化曲線

0.1～0.9充滿率下蒸發率計算結果為：39.7%,26.6%,19.9%,16.1%,13.1%,11.5%,9.5%,8.5%,4.4%。為了預測此LNG儲罐在不同充滿率下的蒸發率，得到多項式函數的擬合:

y=-1.606 9x3+2.966 7x2-1.940 1x+0.558

R2=0.997 4

繪制0.1～0.9不同充滿率下蒸發率的計算曲線，如圖4所示。

圖4 蒸發率隨充滿率變化曲線

為了更直觀地表達罐內LNG液態與氣態的分布，給出罐內氣液態分布云圖如圖5所示。可見，隨著充滿率的增大，LNG蒸發率會不斷變小，在0.9的充滿率下，LNG的日蒸發率為4.4%。環境進入的熱量一方面作用在蒸發一方面作用在為過冷的液相來提高內能，當充滿率不斷增大的時候，需要為過冷的液相增加的內能就會增多，用于蒸發的就會減少，所以其蒸發率會降低。

圖5 充滿率0.5時的蒸汽比例分布云圖

2.2環境溫度影響

對1.5 m3LNG儲罐在不同環境溫度下的蒸發進行仿真，設由外界環境溫度熱傳導形成的壁面溫度的不同來分析，其充滿率為0.9，環境溫度產生的壁面溫度分別取118，119，120，121，122 K 5種不同的情況進行分析。圖6所示為118 K蒸發量變化曲線。

圖6 118 K蒸發量變化曲線

不同的環境漏熱溫度下蒸發率的計算結果如表1所示。

表1 不同環境漏熱溫度下的蒸發率計算結果

可見，當環境溫度增加時，LNG蒸發率也在不斷增大。此結論是由于由外界進入其內部的熱量是其內部與外界溫度大小的差值與保溫層材料自身的導熱系數成正比，是由這兩個已知量所構成的函數。如果環境溫度不斷升高，保溫層內外兩邊的溫度差距的數值也會增大，這樣會導致進入罐內的熱量變多，使氣體的蒸發加快，也就提高了蒸發率。所以，冬天LNG蒸發的速率將低于夏天LNG蒸發的速率，所以在設計過程中，為了安全起見，應當按照夏天的溫度標準來進行LNG儲罐的設計。

2.3儲罐內徑尺寸影響

對1.5 m3LNG儲罐在不同儲罐內徑下的蒸發進行仿真，在充滿率為0.5時，對內徑為0.5,1.0,1.5，2.0，2.5 m 5種情況下不同的蒸發率進行計算。圖7給出了內徑為1.5 m蒸發量變化曲線。

圖7 內徑為1.5 m蒸發量變化曲線

5種不同的內徑下蒸發率的計算結果如表2所示。可見，隨著內徑的增大，LNG的蒸發率在不斷變大(用于小罐的計算)。這是由于隨著內徑的增大，其對應的傳熱面積不同，內徑增大，傳熱面積會變大，繼而漏入的熱量會變多，蒸發就會變快，蒸發率也會升高。

表3 不同內徑下的蒸發率計算結果

3 結 語

本文將LNG儲罐瞬態蒸發率研究內容引入到教學實驗中，設計了基于FLUENT的LNG儲罐瞬態蒸發率仿真實驗，通過仿真實驗，獲得不同充滿率、不同環境溫度、不同內徑尺寸下LNG蒸發率，通過對實驗結果分析，得出：隨著充滿率的增大，LNG的蒸發率在不斷變小；隨著環境溫度的增大，LNG的蒸發率在不斷變大；隨著內徑尺寸的增大，LNG的蒸發率在不斷變大。

通過該仿真實驗，可以給學生創建一個接觸科研前沿、提高實驗技能的平臺，學生通過操作實驗過程、分析實驗結果和撰寫實驗報告，可以從中學習FLUENT的建模、網格劃分、仿真分析等方法，培養學生利用先進仿真軟件解決實際科研問題的能力。將FLUENT仿真融入LNG儲罐瞬態蒸發率實驗教學中，不僅可以促進學生對蒸發規律的理解，還能激發學生的科學探索熱情。這種科研與實踐相結合的教學模式能有效提高實驗教學的質量，培養學生運用實驗手段解決科研問題的能力。

[1] 許邁進，楊行昌．教學與科研并重：研究型大學和諧發展戰略的重要選擇[J]．中國高教研究，2007，27(4)：49-51.

[2] 張德高．科研教學結合為人才培養提供強力支撐[J]．中國高等教育，2013，33(17)：44-45, 54．

[3] 吳再生，吳有訓．大學教育思想及其在清華的實踐:以高水平科學研究支撐的高質量大學教育[J]．清華大學教育研究,2013，33(3)：112-118.

[4] 邵小耿. 1.5 m3獨立C型LNG儲罐的蒸發率研究[D].大連:大連理工大學，2015.

[5] Dimopoulos C A F. A dynamic model for liquefied natural gas evaporation during marine transportation[J]. Interna-tional Journal of Thermodynamics. 2008, 11(3): 123-131.

[6] 張湘鳳，郭 雷，姚 云. LNG低溫儲罐冷量損失在設計中的計算方法[J]. 化工設計, 2012, 21(6): 40-42.

[7] 劉文華，陸 晟. 中小型LNG船C型獨立液貨艙蒸發率計算[J]. 船舶設計通訊, 2012, 25(1): 25-28.

[8] 林建忠. 流體力學[M]. 北京: 清華大學出版社, 2013.

[9] 錢 潔,石嘉豪,費 儉．模擬科研情境開展研究型實驗教學[J]．實驗室研究與探索,2010，29(1)：132-134.

[10] 周立亞,龔福忠,蘭宇衛,等．構建研究型實驗教學法培養學生創新能力[J]．實驗室研究與探索,2011，30(5)：127-129.

[11] 馮根生,馮 婷,杜春榮．基于學科特點的研究型實驗教學模式建設[J]．實驗技術與管理, 2012，29(3)：239-241.

[12] 陳 明,于雁云,李 楷,等. 一種獨立LNG實驗罐裝置:中國, ZL201410290979.0[P].2015.11.18.

[13] 趙斌娟，袁壽其，劉厚林，等. 基于Mixture多相流模型計算雙流道泵全流道內固液兩相湍流[J]. 農業工程學報, 2008，23(1): 7-12.

[14] 肖時雄. 小型LNG船液罐儲運方案中的傳熱與升壓過程研究[D]. 上海：上海交通大學, 2012.

[15] 陳 罕，周昆穎. 應用有限單元法計算瞬時溫度場的基本原理[J]. 北京化工學院學報, 1987，21(2): 63-72.

StudyonSimulationExperimentofTransientEvaporationRateofLNGStorageTankBasedonFLUENT

GUANGuana,LINYana,b,YANGQua,ZHOUShuaia

(a.Ship CAD Engineering Center; b. State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, Liaonin, China)

In order to make students touch professional scientific frontier and better train their scientific research and innovation abilities, the research of transient evaporation rate of LNG storage tank is introduced into experimental teaching, the study on simulation experiment of transient evaporation rate of LNG storage tank based on FLUENT is carried out. The simulation experiment with FLUENT can present the changing rules of the evaporation rate of LNG under fill depth, temperature, variation of the inner diameter of LNG storage tank intuitively. The experiment can promote students’ understanding. The teaching effect can be improved by the teaching experiment combined the teaching with scientific research. It is beneficial to improve students’ ability of simulation, motivate students’ scientific curiosity and cultivate students’ ability to solve scientific research problems by experiment.

LNG storage tank; transient evaporation rate; simulation experiment; FLUENT

U 677.2；G 642.423

A

1006-7167(2017)09-0079-04

2016-12-01

國家自然科學基金資助項目(51609036)；中國博士后科學基金資助項目(2014M561234、2015T80256)；遼寧省博士啟動基金(201501176)資助項目；中央高校基本科研業務費專項資金資助(DUT16RC(4)26)

管 官(1983-),男,滿族,遼寧丹東人，博士，講師，主要從事船舶與海洋工程專業相關教學與科研工作。Tel.:13610926011; E-mail： guanguan@dlut.edu.cn