EN Menu
電(diàn)話(huà):191-1708-8493
郵箱:wynn.yang@vjingguan.com
更多分(fēn)類
技術分(fēn)享 | 基于分(fēn)子動理論的Fluent材料屬性定義
時間:2022年10月25日
浏覽:2131
AIYU.jpg

分(fēn)子動理論方式簡介

Fluent的流體(tǐ)材料屬性定義中(zhōng),有一(yī)種方式叫做kinetic theory(分(fēn)子動理論)。純氣體(tǐ)的比熱容、熱傳導率、粘度可設置爲分(fēn)子動理論;氣體(tǐ)混合物(wù)的組分(fēn)間質量擴散率可設置爲分(fēn)子動理論。

純氣體(tǐ)的設置

氣體(tǐ)混合物(wù)的設置

分(fēn)子動理論方式無需輸入各類宏觀的材料屬性數據,适用于氣體(tǐ)介質存在溫度變化很大(dà)、組分(fēn)數量極多等屬性數值變化較大(dà)、難以查找具體(tǐ)屬性數值的複雜(zá)情況。液體(tǐ)的分(fēn)子間作用機理遠比氣體(tǐ)複雜(zá),因此不推薦使用分(fēn)子動理論方法進行材料屬性定義。


一(yī)
FLUENT
輸入信息


采用分(fēn)子動理論定義材料屬性,前提條件爲密度采用理想氣體(tǐ)方式定義,ideal gas law和incompressible ideal gas law均可,這兩種密度定義方式的差異具體(tǐ)可參考Fluent幫助文檔。需要輸入的數據包括分(fēn)子量、L-J特征長度、L-J能量參數和分(fēn)子自由度。這幾個參數都是分(fēn)子微觀層面的參數,且爲常量,不受溫度、壓力等外(wài)部條件的影響,可查表獲取。

1

L-J參數


分(fēn)子間的作用力表現爲遠距離(lí)引力,近距離(lí)斥力。L-J勢是一(yī)種通過勢能對分(fēn)子間作用力描述的理論,由英國科學家John Lennard-Jones提出而得名。L-J勢描述的分(fēn)子間勢能和距離(lí)關系爲:
式中(zhōng),σ 爲零勢能距離(lí),ε 爲勢阱深度,均爲分(fēn)子的微觀常數。

L-J勢和距離(lí)的函數關系
(圖源:《Modified Lennard-Jones potentials for nanoscale atoms》,《Journal of Computational Chemistry》期刊)

在Fluent中(zhōng),L-J特征長度爲 σ(零勢能距離(lí))的數值,L-J能量參數爲 ε/k(勢阱深度/玻爾茲曼常數)的數值,其中(zhōng) k 爲玻爾茲曼常數。需要注意的是,L-J特征長度的單位爲angstroms(中(zhōng)文簡稱埃,簡寫爲瑞典字符 Å,或者特殊情況下(xià)簡寫爲英文字母A),其換算爲1Å=0.1nm。


2

分(fēn)子自由度


分(fēn)子自由度是指分(fēn)子内部原子的運動自由度,包括平動、轉動、振動三類自由度。不同分(fēn)子類型的自由度總結如表格所示,N表示原子數量。

在不太高的溫度下(xià),含有振動自由度的分(fēn)子比例極低,振動自由度可忽略不計,認爲僅包含平動和轉動兩類自由度。溫度升高時,含有振動自由度的分(fēn)子比例不斷上升,使自由度數量增加。随着溫度進一(yī)步上升,多數分(fēn)子均含有振動自由度,則達到分(fēn)子自由度的上限3N。

通常的氣體(tǐ)完全激發振動自由度的溫度(振動特征溫度,characteristic vibrational temperatures)極高。常見的雙原子氣體(tǐ),僅氯氣的振動特征溫度低于1000K,爲805K,其他常見氣體(tǐ)的振動特征溫度均遠高于1000K,例如氮氣爲3393K,氧氣2274爲K,氫氣爲6332K,一(yī)氧化碳爲3122K。含有多個振動自由度的氣體(tǐ),每個振動自由度均對應于一(yī)個振動特征溫度,不同溫度下(xià)可激發數量不同的振動自由度,情況更加複雜(zá)。

振動特征溫度的計算公式爲

式中(zhōng),h 爲普朗克常數,ν 爲分(fēn)子振動頻(pín)率,k 爲玻爾茲曼常數。


FLUENT
從微觀結構到宏觀數據


1

比熱容


Fluent中(zhōng)比熱容爲等壓比熱容。基于理想氣體(tǐ)假設,等壓比熱容和等體(tǐ)積比熱容關系爲:

式中(zhōng),R 爲通用氣體(tǐ)常數,M 爲氣體(tǐ)分(fēn)子量,f 爲分(fēn)子自由度

從公式中(zhōng)可以看出,比熱容僅受到分(fēn)子自由度影響,即受到溫度的間接影響。在室溫附近的較大(dà)範圍内,比熱容可以認爲是常數。

另外(wài),比熱容是唯一(yī)需要分(fēn)子自由度數據的材料屬性。


2

粘度


基于Chapman–Enskog理論,結合L-J勢理論,得到以下(xià)公式:

式中(zhōng), M 爲氣體(tǐ)分(fēn)子量,T 爲溫度,σ 爲零勢能距離(lí)(L-J特征長度,單位埃),Ω 爲關于T*的碰撞積分(fēn)函數(無初等函數表達式)。其中(zhōng)T* 表達式爲:

式中(zhōng),k/ε 爲L-J能量參數的倒數。 


3

熱傳導率


基于Chapman–Enskog理論,結合Arnold Eucken的工(gōng)作,熱傳導率和比熱容、粘度關系爲:

上式結合理想氣體(tǐ)的等壓比熱容和等體(tǐ)積比熱容關系,可得:

Fluent計算分(fēn)子動理論定義的熱傳導率,其引用的比熱容和粘度數據可以采用常數、多項式、UDF等其他方法定義,不需要都采用分(fēn)子動理論方式。 


4

質量擴散率


基于Chapman–Enskog理論,混合氣體(tǐ)中(zhōng)任意兩種組分(fēn)之間的擴散率爲:

式中(zhōng),p 爲絕對壓力, σij 爲組分(fēn)的L-J特征長度的函數,定義爲:

ΩD 爲擴散碰撞積分(fēn)函數,其表達式爲:

其中(zhōng),Eij 是組分(fēn)的L-J能量參數的函數,定義爲:

式中(zhōng),ε/k 爲組分(fēn)的L-J能量參數。



小(xiǎo)  結




對于基于分(fēn)子動理論方式進行定義的材料屬性,其均僅爲溫度的函數,不受壓力、密度的影響。具體(tǐ)的材料屬性是程序基于分(fēn)子微觀屬性進行的函數計算,因此可以減少搜集資(zī)料、評估數據準确性等工(gōng)作的耗時。






更多信息和資(zī)料請聯系 楊小(xiǎo)姐
19117088493
wynn.yang@vjingguan.com









上海華漁信息科技有限公司是一(yī)個以CAE軟件銷售、技術咨詢及服務,仿真咨詢及規劃布局爲一(yī)體(tǐ)的高科技公司,通用機械、新能源三電(diàn)、高科技、醫療、汽車(chē)有豐富的項目經驗和知(zhī)識積累。

作爲Ansys的合作夥伴,艾羽緻力于将Ansys推出的産品,通過業界性能頗佳、頗豐富的工(gōng)程仿真軟件産品組合幫助客戶解決複雜(zá)的仿真難題。力求與Ansys一(yī)起,共同爲中(zhōng)國制造業提供先進的仿真技術,通過仿真技術支撐中(zhōng)國2025。