王 超,周艷麗,陳英才

(臺州學院物理系,浙江臺州318000)

重力場中理想氣體密度分布的Monte Carlo模擬

王 超,周艷麗,陳英才

(臺州學院物理系,浙江臺州318000)

利用Monte Carlo方法模擬了重力場中理想氣體的密度分布,直觀展現(xiàn)了重力場中氣體分子位置的改變和分布特點,討論了分子質量和系統(tǒng)溫度對氣體密度分布曲線以及重力勢能零點處密度n0的影響.模擬結果與玻耳茲曼分布律完全吻合.另外,模擬結果表明玻耳茲曼分布律不僅對純的理想氣體成立,而且對混合理想氣體中各成分氣體也成立.

理想氣體;密度分布;玻耳茲曼分布;Monte Carlo模擬

1 引 言

氣體分子運動論是大學物理熱學部分的重要知識,其中涉及玻耳茲曼分布律.考慮到玻耳茲曼分布律難于用實驗驗證,很多學校在講授時常采用純粹的數(shù)學推導,這明顯不利于學生對知識的掌握.為此,我們利用Monte Carlo方法[1-2]在計算機上對重力場中理想氣體系統(tǒng)進行模擬,直觀展現(xiàn)氣體分子空間分布圖像,驗證玻耳茲曼分布律.計算機模擬在現(xiàn)代教學過程中發(fā)揮著日益重要的作用[3-5],通過模擬可以幫助學生增加對所學知識的感性認識,鞏固所學理論.

在重力場中,理想氣體系統(tǒng)達到平衡時,分子分布會呈現(xiàn)上疏下密的特征.結合理想氣體狀態(tài)方程和氣柱模型[6-7],可以推出重力場中理想氣體分子數(shù)密度n隨高度h的變化規(guī)律,即玻耳茲曼分布律:

其中m為氣體分子的質量,g為重力加速度,kB為玻耳茲曼常量,T為系統(tǒng)溫度,n0為高度h=0 (零重力勢能)處分子數(shù)密度.n0與分子質量和系統(tǒng)溫度都有關.假定h=0處的壓強為p0,則由理想氣體狀態(tài)方程可得p0=n0kBT;若氣柱的底面積為A,氣柱內的分子數(shù)目為N,則p0又可表示

2 模型的建立及模擬方法

對于理想氣體系統(tǒng)而言,分子之間大量的隨機碰撞使系統(tǒng)平衡時分子具有特定的速度和數(shù)密度分布.考慮到分子速度分布和數(shù)密度分布是相互獨立的[6-7],因此可以撇開分子速度來討論分子的密度分布.分子之間的碰撞使分子空間位置的改變出現(xiàn)隨機性和無規(guī)性,因此在模擬中可用分子位置的隨機改變來描述分子之間的碰撞.重力場的加入使得分子選擇重力勢能小的位置的概率增大,選擇重力勢能大的位置的概率減小.在模擬中用Metropolis算法[1]來刻畫重力場對分子位置選擇概率的改變.另外,分子間的碰撞還反映了分子之間的不可侵入性,因此在模擬中采用排除體積作用來描述分子之間的相互作用,即2個或多個分子不能同時占據(jù)空間同一位置.具體的模型、模擬方法如下:

模型系統(tǒng)建立在二維正方格子點陣上(格子常量為a),重力加速度方向沿z軸負方向,如圖1所示.模擬區(qū)域在x和z方向的尺寸分別為L和H+1.在z方向選取2個邊界z=-1和z=H,而在x方向上采用周期性邊界條件.理想氣體分子只能分布在高度z=-1到z=H區(qū)域內.氣體分子用質點代替,每個分子只能占據(jù)1個格點,分子與分子或分子與邊界不能同時共用1個格點(排除體積作用).假定每個分子的質量為m,并選取z=0處重力勢能為零,那么處在高度z上的分子重力勢能可表示為E=mgz.

圖1 模型系統(tǒng)示意圖

在模擬中,分子的位置改變是通過分子嘗試運動來實現(xiàn)的:假定氣體系統(tǒng)共有N個分子,對于每次嘗試運動,先隨機選擇1個分子,然后再隨機選擇該分子周圍8個近鄰格點中的1個格點作為其新位置.假定新位置沒有被其他分子占據(jù),那么根據(jù)Metropolis算法[1]該嘗試運動成功的概率為min(1,e-ΔE/kBT),其中ΔE代表分子嘗試運動所伴隨的重力勢能的增加量.定義1個Monte Carlo步長為時間單位,在1個Monte Carlo步長內共有N次分子嘗試運動.當系統(tǒng)經過一定時間演化達到平衡后,便終止抽樣模擬并對模擬空間每一高度(z)上的分子數(shù)目(NI)進行統(tǒng)計,進而計算得到每一高度上的分子數(shù)密度n(z)=NI/L.本文所給出的結果都由200次獨立抽樣模擬平均得到.

本文分別選取正方格子常量a、溫度T0= 300 K和氫原子質量m0=1.67×10-27kg為長度、溫度和質量單位,并用z,T＊和m＊表示高度、溫度和質量.假定格子常量a=100 m,同時取重力加速度g=10 m/s2、玻耳茲曼常量kB=1.38× 10-23J/K,則式(1)可作如下變換:其中無量綱常數(shù)C=4.0×10-4.選擇模擬區(qū)域尺寸L=1 000,H=499.如果不做特別說明,模擬區(qū)域的理想氣體分子數(shù)目均為N=1 000.

3 模擬結果及討論

首先模擬了溫度T＊=1.000時氧氣分子的空間分布.在無重力情形下,模擬區(qū)域內的每個分子向各個位置運動具有相同的概率,因此系統(tǒng)平衡時分子在模擬區(qū)域內呈均勻分布,如圖2所示.這與實際氣體系統(tǒng)在無重力條件下分子空間分布特點相一致.

圖2 無重力時氧氣分子的空間平衡分布

加入重力場后,分子向重力勢能減小方向運動的概率增大,而向勢能增大方向運動的概率減小,原來的均勻分布平衡被破壞.經過一定時間的演化,體系將達到新的平衡,分子的空間位置也將出現(xiàn)新的分布.圖3給出了加入重力場后模擬區(qū)域內幾個高度上的氧氣分子數(shù)目隨時間的演化(t=0對應剛加入重力場),其中虛線對應無重力場時各高度上的平均分子數(shù).由圖3可以看出,加入重力場后,模擬區(qū)域下部(z值小)各高度上分子數(shù)不斷增加,而上部(z值大)各高度分子數(shù)不斷減小.隨著下部分子數(shù)增多,分子之間的排除體積作用越來越強.當重力與分子之間的排除體積作用相平衡時,各高度上的分子數(shù)達到飽和,系統(tǒng)也達到新的平衡,在整個模擬區(qū)域氧氣分子形成了上疏下密的空間分布,如圖4所示.系統(tǒng)達到平衡并不意味各高度上的分子數(shù)目恒定不變.由圖3還可以看出,各高度上的分子數(shù)達到飽和后,分子數(shù)仍存在一定的漲落,這說明氣體系統(tǒng)的平衡是一種動態(tài)平衡.為了進一步說明這一點,當系統(tǒng)平衡時,隨機選擇1個氧氣分子,記錄其高度隨時間的演化,如圖5所示.由圖5可以看出,即便整個體系處于平衡,但就每個分子而言其高度仍可在很大的范圍內變化,只不過分子處在低位置的概率要比高位置的大.當溫度T＊=1.000時,由式(3)可以得到氧氣(m＊=32)的密度隨高度的變化關系:n/n0= exp(-0.012 8z).圖6給出了氧氣分子數(shù)密度分布的模擬結果,其中曲線是由上述理論公式得到的結果.由圖6可以看出,模擬結果與理論曲線幾乎完全吻合.

圖3 不同高度上分子數(shù)目隨時間的演化

圖4 有重力時氧氣分子的空間平衡分布

圖5 平衡態(tài)下氧氣分子的高度隨時間的演化

由式(1)可以看出分子質量以及系統(tǒng)溫度都會對氣體密度分布產生影響.為了驗證這一問題,分別模擬了相同分子數(shù)(N=1 000)的氧氣(O2)、二氧化碳(CO2)和二氧化硫(SO2)3種氣體在一定溫度下的密度分布,同時模擬了一定分子數(shù)(N=1 000)的氧氣(O2)在不同溫度下的密度分布,圖7和圖8分別給出了對應的模擬結果.可以直觀看到:質量和溫度都會對密度分布圖線的陡峭程度產生影響.分子質量越大,密度分布圖線越陡峭;系統(tǒng)溫度越高,密度分布圖線越平緩.處在重力場中的氣體分子一方面做無規(guī)則熱運動,一方面受重力作用.質量越大,分子受到的重力越大,向下運動的概率就越大,從而使模擬區(qū)域下部的分子數(shù)越多,密度分布曲線越陡峭;溫度越高,分子的熱運動就越劇烈,克服重力作用而達到模擬區(qū)域上部的分子數(shù)目也就越多,從而使得密度分布曲線越平緩.

圖7 不同質量的氣體分子的密度分布圖線

圖6 氧氣分子的密度分布圖線

圖8 不同溫度下氧氣分子的密度分布圖線

由圖7和圖8容易看出,質量和溫度不僅會對密度分布圖線的陡峭程度產生影響,而且對n0也有影響.質量越大,處在z=0高度的分子越多,n0就越大;溫度越高,處在z=0高度的分子越少,n0就越小.式(2)表明n0與質量成正比,與溫度成反比.為了驗證這一關系,我們分別模擬研究了n0隨m＊以及T＊的變化,如圖9和圖10所示.由圖可以看出,隨著分子質量(m＊)增大或溫度(T＊)減小,n0不斷增大.另外,由模擬結果可以直觀得到:n0∝m＊,n0∝1/T＊,這就很好地驗證了式(2).

圖9 lgn0與lgm＊的關系

圖10 lgn0與lgT＊的關系

上面討論了純的理想氣體在重力場中的密度分布.很明顯,模擬不僅直觀展現(xiàn)了氣體分子位置改變和分布特點,而且定量驗證了玻耳茲曼分布律.接下來模擬研究混合理想氣體在重力場中的密度分布.在實際模擬中,同時將氧氣(O2)、二氧化碳(CO2)和二氧化硫(SO2)3種氣體各1 000個分子放入模擬區(qū)域,考查溫度T＊=1.000時各種成分氣體以及混合氣體平衡態(tài)分子密度分布,模擬結果如圖11所示,其中插圖代表混合氣體總的分子數(shù)密度分布.由圖可以看出:當系統(tǒng)處于平衡時,混合氣體和各種成分氣體分子數(shù)密度分布都保持上疏下密的特點,這說明在重力場中不管是純的氣體還是混合氣體,氣體分子都會向重力勢能小的地方聚集.另外,通過對比圖11和圖3發(fā)現(xiàn)混合氣體中各成分氣體分子數(shù)密度分布與對應的純氣體分子數(shù)密度分布幾乎完全相同,這說明玻耳茲曼分布關系不僅適用于重力場中純的理想氣體,而且適用于混合氣體中各成分氣體.

圖11 混合氣體中各成分氣體密度分布圖線(插圖為混合氣體總的密度分布圖線)

4 結束語

本文用Monte Carlo方法模擬了重力場中理想氣體的密度分布.模擬直觀展現(xiàn)了重力場中氣體分子位置的演化,揭示了氣體分子位置的改變和分布特點.模擬得到的氣體密度分布曲線與玻耳茲曼分布律完全一致.模擬結果表明,分子質量(m＊)和系統(tǒng)溫度(T＊)對分子密度分布產生重要影響,m＊和T＊不僅影響分布曲線的陡峭程度而且影響重力勢能零點處分子密度n0的大小. n0與分子質量成正比,與溫度成反比.模擬結果與理論結果描述的關系完全吻合.另外,我們還討論了混合理想氣體密度分布,模擬結果顯示混合氣體中各成分氣體的密度分布與純的理想氣體的密度分布相一致,這說明玻耳茲曼分布律不僅對純的理想氣體成立,而且對混合理想氣體中各成分氣體也成立.

[1] Metropolis N,Rosenbluth A W,Rosenbluth M N, et al.Equation of state calculations by fast computing machines[J].J.Chem.Phys.,1953,21: 1 087-1 092.

[2] 王道光,周紅生,李松嶺.Monte Carlo法模擬油滴的理論分布[J].物理實驗,2008,28(9):26-29.

[3] 吳俊,吳本科,謝莉莎,等.微波單縫衍射的實驗研究及數(shù)值模擬[J].物理實驗,2008,28(6):39-41.

[4] 蔡星漢,繆騰飛,周進.字母屏的夫瑯禾費衍射分析及計算機模擬[J].物理實驗,2008,28(12):38-42.

[5] 孫存英,喬衛(wèi)平,沈學浩.偏振光的干涉場分布特征數(shù)值模擬及實驗研究[J].物理實驗,2009,29(9): 42-45.

[6] 費曼R P,萊登R B,桑茲M.費曼物理學講義(第一卷)[M].上海:上?？茖W技術出版社,1983:390-392.

[7] 趙凱華,羅蔚茵.熱學[M].北京:高等教育出版社,2000:79-81.

Monte Carlo simulation of density distribution of ideal gas in gravity field

WANG Chao,ZHOU Yan-li,CHEN Ying-cai
(Department of Physics,Taizhou University,Taizhou 318000,China)

The density distribution of ideal gas in gravity field is simulated using the Monte Carlo method.The characteristics of the moving process and distribution of ideal gas molecules in gravity field are displayed.The influence of molecule mass and system temperature on the density distribution and the densityn0at the position where the gravity potential energy is defined zero are discussed.The simulation results are in good agreement with the Boltzmann distribution.In addition,the simulation results indicate that the Boltzmann distribution hold true not only for pure ideal gas,but also for each component in mixed ideal gas.

ideal gas;density distribution;Boltzmann distribution;Monte Carlo simulation

O4-39;O552.3

A

1005-4642(2010)10-0024-05

[責任編輯:郭 偉]

2009-12-28;修改日期:2010-04-05

王 超(1981-),男,河南南陽人,臺州學院物理系講師,碩士,主要從事物理教學及高分子模擬研究工作.