亚/超临界环境下气液界面性质的分子动力学模拟

发布时间：2021-03-23 03:55

　　采用OPLS-AA力场对正庚烷在亚/超临界环境下的气液界面性质进行分子动力学模拟。真空环境下,得到气液相密度、界面厚度及界面张力等性质随模拟分子数、截断半径及模拟温度的变化规律,并与实验值进行对比,验证了力场模型的合理性;在真空以及由氮气形成的亚/超临界环境下,给出流体由亚临界过渡到超临界状态时的温度分布情况,总结了气液界面性质随模拟环境的变化规律:随模拟环境由亚临界过渡到超临界,气、液相密度差及界面张力减小,界面厚度增大。在低超临界状态下,界面性质与亚临界状态尚无本质区别,仅当达到较高超临界状态时,流体才表现出明显的超临界特征。 

【文章来源】：工程热物理学报. 2016,37(08)北大核心EICSCD

【文章页数】：6 页

【部分图文】：

图I平衡模拟体系示意图(a)真空环境(b)亚临界环境Fig.1Equilibriumsimulationsystem:(a)vacuumsurroundings(h)subcriticalsurroundings

８期?邓磊等：亚／超临界环境下气液界面性质的分子动力学模拟?１８０５??了不同温度下界面张力７的模拟数值与实验值的对?表３所示？系统沿方向的温度分布如图５所示：由??比，二者基本吻合，说明ＯＰＬＳ－ＡＡ力场模型能较准?高温气相过渡到低温液相时温度未发生突变，气液??确合理地模拟正庚烷气液界面性质。?界面处较大的温度梯度导致温度逐渐连续的变化，??且温度梯度随模拟时间增加而减小，即正庚烷和氮??８００?ｒ￣￣．——■￣￣￣．——■￣．——■￣?Ｉ—??７〇〇?－?？你ｉｄ．?▼?ｒ＝２３７?Ｋ?－?３．０．，????３．０ｎ￣￣￣￣ｎ??＾１５：ｒ＝２７＾?ｉ?：?１５Ｔ＝２９３?ｆ?＾??，＃。?：ｕ?ｉ。。??１４〇〇－?？?ｉ?畏。。＿ｉ？珞厂：ｔ?ａ°—厂一??二＇ｔ?Ｗ?—１?５ｌ?■?■?＿￣￣Ｉ?—１?５＇一？一？?真■￣ｌ＿？＿ｉ—＞＿ｔ＿??３００?■?ｔ?＇?０?２?４?６?８?１０?０?２?４?６?８?１０??２〇〇?＾?Ｚ?：?ｘ／ｍｍ?ｘ／ｍｍ??１００「?」ｊ?＇?Ｔ?３?〇ｌ?ｆ＝３１３?Ｋ?Ｉ?－?３＇〇ｌ?Ｖ＝３３３?Ｋ??０－?Ｐ?１．５．?）?．：？?■?￡?１．５－?；，??－１００?〇?２?４?６?８?１０?＜?〇．〇?—－?｜＿?°－°?—??ｘ／ｍｍ?＿ｉ?５Ｌ???????一?１?５！?■?■?■?？??????＾?０?２?４?６?８?１０?＇?０?２?４?６?８?１０??图２不冋温度Ｆ密度分布曲线?；ｃ／ｍｍ?ａＶｉｔｕｔｉ??Ｆｉｇ．?２?Ｄｅｎｓｉｔｙ?ｄｉｓｔｒｉｂｕｔｉｏｎ?ｏｆ?ｖａｃｕｕｍ?ｓｕｒｒｏｕ

８期?邓磊等：亚／超临界环境下气液界面性质的分子动力学模拟?１８０５??了不同温度下界面张力７的模拟数值与实验值的对?表３所示？系统沿方向的温度分布如图５所示：由??比，二者基本吻合，说明ＯＰＬＳ－ＡＡ力场模型能较准?高温气相过渡到低温液相时温度未发生突变，气液??确合理地模拟正庚烷气液界面性质。?界面处较大的温度梯度导致温度逐渐连续的变化，??且温度梯度随模拟时间增加而减小，即正庚烷和氮??８００?ｒ￣￣．——■￣￣￣．——■￣．——■￣?Ｉ—??７〇〇?－?？你ｉｄ．?▼?ｒ＝２３７?Ｋ?－?３．０．，????３．０ｎ￣￣￣￣ｎ??＾１５：ｒ＝２７＾?ｉ?：?１５Ｔ＝２９３?ｆ?＾??，＃。?：ｕ?ｉ。。??１４〇〇－?？?ｉ?畏。。＿ｉ？珞厂：ｔ?ａ°—厂一??二＇ｔ?Ｗ?—１?５ｌ?■?■?＿￣￣Ｉ?—１?５＇一？一？?真■￣ｌ＿？＿ｉ—＞＿ｔ＿??３００?■?ｔ?＇?０?２?４?６?８?１０?０?２?４?６?８?１０??２〇〇?＾?Ｚ?：?ｘ／ｍｍ?ｘ／ｍｍ??１００「?」ｊ?＇?Ｔ?３?〇ｌ?ｆ＝３１３?Ｋ?Ｉ?－?３＇〇ｌ?Ｖ＝３３３?Ｋ??０－?Ｐ?１．５．?）?．：？?■?￡?１．５－?；，??－１００?〇?２?４?６?８?１０?＜?〇．〇?—－?｜＿?°－°?—??ｘ／ｍｍ?＿ｉ?５Ｌ???????一?１?５！?■?■?■?？??????＾?０?２?４?６?８?１０?＇?０?２?４?６?８?１０??图２不冋温度Ｆ密度分布曲线?；ｃ／ｍｍ?ａＶｉｔｕｔｉ??Ｆｉｇ．?２?Ｄｅｎｓｉｔｙ?ｄｉｓｔｒｉｂｕｔｉｏｎ?ｏｆ?ｖａｃｕｕｍ?ｓｕｒｒｏｕ

【参考文献】：
期刊论文
[1]内燃机跨临界/超临界燃料喷雾混合过程的机理与模型[J]. 解茂昭.  燃烧科学与技术. 2014(01)



本文编号：3095083