2模擬結果與討論

2.1溫度對密度分布的影響

在模擬分子數N=256和截斷半徑rc=0.9498 nm的條件下，當溫度T=400、450、500、550和610 K時，模擬得到的密度分布如圖3所示。統計得到的汽相主體密度ρV、液相主體密度ρL及汽液界麵厚度d如表2所示。由圖3及表2可見，汽相主體密度和汽液界麵厚度隨溫度的提高而增加，而液相主體密度隨溫度的提高而減小。

液相主體密度與汽相主體密度之差(ρL－ρV)與溫度T的關係如圖4所示。可見，液、汽相主體密度之差隨溫度的升高而降低;從理論上講，在臨界點處，其差值應該趨近於零，這與圖3所示的規律一致。液、汽相主體密度之差與溫度的關係可以擬合成式(14)的形式。

式中水臨界溫度Tc=647.3 K，利用表2數據對式(14)進行擬合，得到參數ρ0=1545.8 kg/m3，指數因子x=0.5516。

2.2溫度對界麵張力的影響

在模擬分子數N=256和截斷半徑rc=0.9498 nm的條件下，當溫度T=400、450、500、550和610 K時，水汽液界麵張力的模擬結果見表3。

圖5為局部界麵張力的分布曲線(500 K)。由圖5可見，汽相主體的局部界麵張力基本為零;從汽相主體向液相主體的過渡過程中，界麵張力值逐漸增加，在汽液界麵區達到峰值;在液相主體又在零值附近波動。水汽液界麵張力模擬值隨溫度變化規律如圖6所示。

由圖6可以看出，隨著溫度的提高，界麵張力降低，模擬值與實驗值之間的誤差逐漸減小。界麵張力與溫度的關係可以擬合得到方程(15)。

將表3的數據對式(15)進行擬合，得到的參數γ0=254.3 mN·m－1，指數因子k=1.305。

2.3溫度對勢能分布的影響

在模擬分子數N=256和截斷半徑rc=0.9498 nm的條件下，當溫度T=400、450、500、550和610 K時，汽相主體總勢能UV、液相主體總勢能UL及總勢能勢阱深度ΔU的模擬結果如表4所示。圖7為水分子的勢能分布曲線(500 K)，圖8為液相主體區域的勢能隨溫度的變化趨勢。

圖8液相主體區域的勢能隨溫度的變化趨勢

前已述及，水的勢能分為L－J勢能和靜電勢能。由圖7可以看出，L－J勢能均為正值，在液相區形成勢壘，勢壘高度ΔULJ為液相主體L－J勢能與汽相主體L－J勢能之差;靜電勢能均為負值，在液相區形成勢阱，勢阱深度ΔUe為汽相主體靜電勢能與液相主體靜電勢能之差;由於靜電勢能起主導作用，總勢能也為負值，同樣在液相區形成勢阱，分子之間主要為吸引作用。從圖8可以看出，汽相主體勢能作用不明顯，勢壘高度隨溫度升高而降低，液相主體勢能的勢阱深度隨體係溫度的升高而減小。

2.4模擬分子數對模擬結果的影響

在溫度500 K和截斷半徑rc=0.9498 nm的條件下，當模擬分子數N=108、256、500和864時，模擬得到的密度分布見圖9。統計得到的汽相主體密度ρV、液相主體密度ρL及汽液界麵厚度d見表5。

圖9水分子數對密度分布的影響

表5不同水分子數下界麵性質的模擬結果

由表5和圖9可見，隨著模擬分子數的增加，液相主體密度有所增加，液相主體區域寬度加大，汽液界麵厚度稍有增大，汽相主體密度有所波動。

2.5截斷半徑對模擬結果的影響

在溫度為500 K和模擬分子數為864的條件下，當截斷半徑rc=0.7915、0.9498、1.2660 nm時，模擬得到的密度分布如圖10所示。統計平均得到的汽相主體密度ρV、液相主體密度ρL及汽液界麵厚度d如表6所示。從表6和圖10可以看出，隨著截斷半徑的增加，液相主體密度增大，汽相主體密度減小，汽液界麵厚度變化不大。

3結論

采用SPC模型，對水汽液界麵特性的分子動力學模擬研究結果表明，隨著溫度的升高，汽相主體密度增加，汽液界麵厚度增大，液相主體密度降低，界麵張力逐漸減小，液相主體區域勢能的勢阱深度也逐漸降低。隨著模擬分子數的增加，液相主體密度增加，汽液界麵厚度稍有增大。隨著截斷半徑的增加，液相主體密度增加，汽液界麵厚度變化不大。