微納尺度氣體流動和傳熱的Burnett方程研究.pdf

1、在微納尺度氣體流動中,通常用Knudsen數,Kn,來表征氣體的稀薄程度,Kn數定義為氣體分子的平均自由程與流場特征尺度之比。微納尺度氣體流動中,流動進入滑移過渡流區(qū),流場中缺乏足夠的分子碰撞,流體逐漸偏離熱平衡狀態(tài)。此時,基于線性本構關系的Navier-Stokes(N-S)方程得到的結果與實驗結果有較大的偏差,而基于顆粒模型的DSMC方法統(tǒng)計噪聲大,收斂慢,需要大量的計算時間。本論文采用Burnett方程來研究微納尺度通道里的氣體流

2、動和傳熱特性。Burnett方程由Chapman-Enskog展開從Boltzmann方程導出,能夠較好地描述稀薄氣體的非線性本質,近二十年來,Burnett方程逐漸受到關注。 首先采用線性小擾動理論,首次系統(tǒng)分析常規(guī)Burnett方程、增廣Burnett方程、BGK-Burner方程以及Woods方程的—維穩(wěn)定性問題。研究發(fā)現這四種方程在沒有引入對方程中隨體導數的近似時,對于小擾動都是不穩(wěn)定的,不穩(wěn)定的臨界Kn數分別是0.10

3、5,0.105,0.158和0.074。當采用Euler方程來近似這些方程里的隨體導數時,增廣Burner方程和BGK-Burnett方程對小擾動是無條件穩(wěn)定的,而常規(guī)Burnett方程和Woods方程是不穩(wěn)定的,臨界Kn數分別是0.499和0.184。當用N-S方程來近似這些方程里的隨體導數時,所有四個方程對于小擾動都是無條件穩(wěn)定的。 接著采用增廣Burnett方程結合二階滑移邊界條件研究微納尺度下的平板Couette流動,在

4、邊界上采用松弛方法,成功地獲得了任意網格尺度時所有Kn數下的收斂解。結果顯示Burnett方程的結果與DSMC以及IP方法的結果符合得較好,當Kn數較小時,N-S方程的計算結果與Burnett方程的結果基本符合,但是隨著Kn數的增大,N-S方程的結果開始偏離DSMC結果,而Burnett方程仍舊能夠較好地符合,Burnett方程更適合模擬高Kn數下的微納尺度流動。文中還分析比較了不同Kn數和馬赫數下的流場特性,給出了壁面上的滑移速度和溫

5、度躍變隨Kn數和馬赫數的變化。 最后采用Burnett方程結合高階通用滑移邊界條件研究微納尺度下的平板Poiseuille流動,首次在Burnett項中引入松弛方法,成功地拓寬了Burnett方程的求解范圍,獲得了Kn≤0.4時Poiseuille流的收斂結果。比較發(fā)現,Burnett方程的結果與理論分析結果、實驗結果以及DSMC結果符合得較好。當Kn數比較小時,Bumett方程和N-S方程的結果相符,但是當Kn數變大,達到過渡