借助Matlab軟件，用蒙特卡洛方法對圓截面直角彎管在分子流態(tài)下的傳輸幾率進行了模擬計算。通過數學方法實現(xiàn)對氣體分子的虛擬約束和跟蹤。本文先通過模擬最簡單的圓直管道傳輸幾率，并與Clausing、Dushman方法進行比較，驗證了該方法的正確性和建模的合理性，進而對于圓截面直角彎管的傳輸幾率進行了模擬計算，將模擬結果與兩種等效算法比較得出：在分子流下計算圓截面直角彎管的傳輸幾率時，Davis的等效方法較傳統(tǒng)的等效方法更為合理。

　　在真空系統(tǒng)設計與計算中，為了表征稀薄氣體通過真空系統(tǒng)管路元件的流動，通常給出流導幾率(即傳輸幾率)，流導幾率是確定氣體流量的一個重要參數。圓截面直角彎管是真空系統(tǒng)常用的管道結構，其在分子流狀態(tài)下的傳輸幾率是衡量真空設備和真空性能的重要參數。

1、圓截面直角彎管傳輸幾率的模擬計算

　　1.1、基本假設及概率模型

　　基于管道中的氣流狀態(tài)為分子狀態(tài)，進行如下假設：

　　(1)氣流為穩(wěn)定氣流，氣體分子數守恒，即管壁無吸氣和放氣現(xiàn)象。這意味著射入管口的分子最終只有兩種可能：從出口逸出或者從入口逸出。兩者的幾率之和等于1。

　　(2)入射分子和反射分子都遵循余弦定律。

　　(3)分子在管道內的運動是相互獨立的，即分子之間互不碰撞，氣體分子只與管壁發(fā)生碰撞。

　　(4)評定參數C約等于1，忽略氣體分子在分子流態(tài)下通過直圓管道的位置束流效應。

　　因為氣體以分子流態(tài)流動，就每個分子而言，從分子飛入管道與管壁碰撞后產生漫反射直至分子逸出管道，分子的整個運動過程都是隨機的。故管道的傳輸幾率本身就是一種概率統(tǒng)計問題。每個分子的隨機運動都可以用一個隨機變量來表示，通常在計算機上采用(0，1)區(qū)間均勻分布的偽隨機數進行抽樣，用數學方法模擬每個分子的運動過程，根據計算機跟蹤每個分子，統(tǒng)計進入管道的分子總數N和逸出管道出口的分子數n，可以得到管道的傳輸幾率Pr

Pr=n/N(1)

　　N越大，Pr越準確。當N足夠大時，Pr就足夠準確。

　　1.2、圓直管道傳輸幾率

　　在計算圓截面直角彎管的傳輸幾率之前，先進行最簡單的圓直管道傳輸幾率計算并與Clausing、Dushman方法計算的傳輸幾率進行對比，驗證MonteCarlo法計算傳輸幾率的準確性。圓直管道傳輸幾率的計算方法與圓截面直角彎管的傳輸幾率計算中的橫管的部分相似，只不過不用判斷分子是否進入縱管，而是直接判斷其是否從出口飛出，具體方法在下文中詳述。

　　利用Matlab軟件進行編程模擬，計算出不同長徑比的圓直管道的傳輸幾率，將這些數據與Clausing、Dushman方法計算的相應傳輸幾率繪制在同一圖中，如圖1所示。

圖1 MonteCarlo法計算的傳輸幾率Pr和Clausing系數Kc、Dushman的對比

　　從圖1中可以看出，用MonteCarlo法模擬計算的傳輸幾率Pr與Clausing積分方程的近似解有很強的一致性。在L/R=1.5處出現(xiàn)最大殘差0.0105，造成約1.8%的相對誤差。相對誤差最大值出現(xiàn)在L/R=17處，達到了5.23%。而Dushman方法計算出來的傳輸幾率與其他兩者差別很大，在L/R=5.1處相對誤差達到了13.52%。隨著L/R逐漸增大，三條曲線逐漸趨于一致。用MonteCarlo法模擬計算的圓直管的傳輸幾率與用Clausing方程計算的近似解十分接近，這證明了該方法模擬計算管道的傳輸幾率的準確性。因此，可以用類似上述的模型對圓截面直角彎管的傳輸幾率進行模擬計算。

2、結論

　　本文主要采用MonteCarlo法模擬計算分子流態(tài)下管道的傳輸幾率，先驗證MonteCarlo法計算圓直管的傳輸幾率的準確性，再通過與Clausing積分方程的近似解和Dushman近似計算的結果進行對比，得到模擬結果與Clausing積分方程的近似解具有較好的吻合性，Dushman近似計算的結果與MonteCarlo模擬結果相差較大。然后采用該方法模擬計算圓截面直角彎管的傳輸幾率，并與目前常采用的傳統(tǒng)的等效長度法、Davis等效法進行對比，傳統(tǒng)的等效長度法與模擬結果相差較大，Davis等效法與模擬結果具有較好的吻合性。因此在理論上計算分子流態(tài)下的彎管的傳輸幾率時，可以根據Davis等效法來進行計算。