壁面效應(yīng)對納米尺度氣體流動的影響規(guī)律研究*

張燁 張冉 常青 李樺

(國防科技大學(xué)空天科學(xué)學(xué)院,長沙 410073)

1 引 言

近年來,隨著微納機電系統(tǒng)(MEMS/NEMS)的快速發(fā)展,納米尺度氣體流動問題受到廣泛關(guān)注,逐漸成為一個重要的研究領(lǐng)域[1-4].納米尺度氣體流動具有廣闊的應(yīng)用前景,比如生物樣本提取、集成電路冷卻以及流動主動控制[5].然而,納米尺度下的流動特性與宏觀流動相比有很大差異.為了理解納米尺度流動特性,有必要對其機理展開深入研究.

氣體的稀薄程度一般用克努森數(shù)(Kn=λ/L,其中λ為氣體分子平均自由程,L為流動特征長度)來表示.微納尺度氣體流動中,分子平均自由程與流動特征尺度相比不再是小量,連續(xù)介質(zhì)假設(shè)不再成立,氣體表現(xiàn)出明顯的稀薄效應(yīng).同時,尺度的減小會導(dǎo)致面積與體積比增大,與壁面相關(guān)的物理量對流動的影響更加突出.當Kn≤0.01 時,流動處于滑移流區(qū)域,Navier-Stokes (NS)方程在主流區(qū)域仍然有效,僅在壁面處出現(xiàn)速度滑移,附加具體的滑移邊界條件仍然可以求解.Maxwell滑移邊界條件[6]以及一些能夠反映固體界面處流體流動特性的滑移模型[7-9]已經(jīng)在工程中得到廣泛應(yīng)用,但是NS方程的線性本構(gòu)關(guān)系決定了滑移邊界條件對于壁面附近物理量非線性分布的描述不足.隨著克努森數(shù)的增大,平均自由程和特征長度處于同一量級,流動處于過渡流區(qū)域,克努森層變厚,壁面效應(yīng)影響更加突出,流場表現(xiàn)出強烈的非平衡和非線性特征,基于連續(xù)性假設(shè)的NS方程徹底失效.而描述稀薄氣體的最完備的玻爾茲曼方程由于自變量維數(shù)高,數(shù)學(xué)形式復(fù)雜,求解起來十分困難.對玻爾茲曼方程高階展開得到的近似方程比如Burnett方程、Super-Burnett方程,雖然較NS方程有一定優(yōu)勢,但仍然不能完全描述過渡區(qū)稀薄氣體的強非線性關(guān)系和非平衡效應(yīng).且由于穩(wěn)定性等諸多問題沒有解決,表達形式復(fù)雜,物理意義不明確,高階方程難以推廣應(yīng)用.因此,人們考慮采用數(shù)值模擬方法研究稀薄氣體流動.直接模擬蒙特卡羅(direct simulation Monte Carlo,DSMC)方法[10,11]是以分子混沌假設(shè)為前提的一種基于分子模型的方法,能夠有效模擬過渡流區(qū)域稀薄氣體流動.但是DSMC方法在處理邊界處氣體與壁面的相互作用時,仍然需要提供動量和能量適應(yīng)系數(shù).對于不同的氣體和壁面,其適應(yīng)系數(shù)相差很大,因此該方法不適合研究壁面占主導(dǎo)因素的納米尺度流動.

分子動力學(xué)(molecular dynamics,MD)方法[12]作為一種確定論的模擬方法,從原子運動的角度出發(fā),根據(jù)牛頓第二定律和分子間勢能函數(shù)追蹤每個分子的運動,不僅能獲得流場的流動細節(jié),還能獲得氣體分子在壁面處的散射規(guī)律,因此在研究微納尺度流動方面具有獨到的優(yōu)勢.基于MD模擬,許多學(xué)者針對微納通道流動開展了廣泛的研究.曹炳陽等[13,14]采用MD模擬方法對滑移區(qū)微納通道內(nèi)的氣體流動進行了研究,探討了速度滑移與勢能作用強度、系統(tǒng)溫度以及壁面粗糙度之間的關(guān)系.Priezjev[15]采用MD方法研究了具有納米尺度紋理的壁面構(gòu)成的通道內(nèi)的液體流動規(guī)律,發(fā)現(xiàn)近壁面第一層液體的有效滑移長度和耗散系數(shù)依賴于濕壁區(qū)域的面積.Sun和Li[16]采用非平衡分子動力學(xué)方法計算了光滑微納通道中的切向動量適應(yīng)系數(shù)和法向動量適應(yīng)系數(shù),通過等溫伯肅葉流動分析了溫度、氣體-壁面作用長度以及氣體分子的吸附現(xiàn)象對動量適應(yīng)系數(shù)的影響.Spijker等[17]采用MD方法模擬了納米通道內(nèi)的氣體流動,統(tǒng)計了氣體分子與壁面碰撞前后的速度分布,并在分析入射與反射速度間的關(guān)系時引入了速度相關(guān)性云圖.Xie和Liu[18]模擬了納米通道中混合氣體的Poiseuille流動,發(fā)現(xiàn)壁面作用力使得通道內(nèi)氣體混合物呈現(xiàn)出化學(xué)成分與物理結(jié)構(gòu)不再均勻的特殊現(xiàn)象.Kamali和Kharazmi[19]用MD方法研究了不同的壁面親疏水性和粗糙單元長寬比條件下,液氬納米尺度伯肅葉流動的密度和速度分布.Barisik和Beskok[20,21]引入了一種智能壁面分子動力學(xué)方法,假設(shè)壁面原子靜止,發(fā)現(xiàn)通道內(nèi)近壁區(qū)域氣體的密度、速度與壓力的變化僅由壁面力場確定,不受通道高度影響的規(guī)律.Noorian等[22]應(yīng)用MD方法研究了液氬在微納通道中的流動,結(jié)果顯示壁面吸引能和粗糙高度的增加,會加強壁面附近處的密度分層現(xiàn)象.Bao等[23,24]采用MD方法對有限長度納米通道內(nèi)的壓力驅(qū)動流動進行研究,得到了壓力差與溫度和克努森數(shù)無關(guān)而與通道縱橫比及壁面性質(zhì)相關(guān)的結(jié)論.To等[25]用MD模擬方法研究了微納通道中非平衡氣體的平均自由程分布情況,發(fā)現(xiàn)平均自由程在壁面附近顯著減小,而對流動類型不敏感.Liakopoulos等[26]采用非平衡分子動力學(xué)方法研究了壁面的親疏水性和粗糙度對納米通道內(nèi)摩擦系數(shù)的影響.Lim等[27]用MD方法模擬了He,Ne,Ar,Kr和Xe氣體的流動,研究了分子質(zhì)量和吸附效應(yīng)對切向動量適應(yīng)系數(shù)的影響,發(fā)現(xiàn)切向動量適應(yīng)系數(shù)隨分子質(zhì)量的變化并不單調(diào).王勝等[28]采用MD方法研究了液體在非對稱性浸潤性粗糙納米通道內(nèi)的流動與傳熱過程,發(fā)現(xiàn)隨著兩側(cè)壁面浸潤性不對稱程度的增加,流體反轉(zhuǎn)溫度分布更加明顯.

從前人的研究來看,壁面的影響在微納通道流動中占據(jù)著主導(dǎo)作用.為了進一步研究壁面效應(yīng)對納米尺度氣體流動的影響,我們在之前的工作中研究了通道高度和壁面勢能系數(shù)對納米通道內(nèi)氣體剪切流動的影響[29],但是對于具有一定剛度和粗糙度構(gòu)型的壁面,其對流場物理量的分布以及近壁區(qū)域流動特性的影響研究還不是很充分.實際上通道壁面的原子并不是靜止的,而是存在著分子熱運動.納米通道的壁面通常用Einstein固體理論[30]來構(gòu)造,該模型通過彈性系數(shù)來模擬壁面原子的熱運動,盡管有一些學(xué)者研究了彈性系數(shù)對納米通道內(nèi)液體流動的影響規(guī)律[31-33],但是在大部分微通道氣體流動研究中,由于壁面原子量巨大,出于節(jié)省計算量的目的,都將壁面原子設(shè)置為靜止不動[20,21],這種處理無疑忽略了壁面原子振動對氣體流動的影響,因此有必要研究壁面剛度對氣體流動的影響.此外,除了少數(shù)高精度加工的壁面外,絕大多數(shù)的壁面都是粗糙的.對于納米通道內(nèi)的氣體流動,粗糙度的尺寸相比于流動的特征長度并不是小量,粗糙度對氣體分子與壁面的能量和動量交換有著更加突出的影響,且粗糙度的存在增大了壁面力場的穿透深度,因此研究粗糙度對近壁區(qū)域氣體流動的影響具有重要意義.系統(tǒng)溫度的變化對納米通道稀薄氣體流動的影響主要通過兩方面來體現(xiàn).一方面,氣體溫度的變化會影響通道內(nèi)氣體分子之間的碰撞頻率,而碰撞頻率的變化又會導(dǎo)致氣體的黏性系數(shù)發(fā)生改變,進而對流動產(chǎn)生影響;另一方面,溫度的變化會影響氣體分子與壁面的動量和能量交換,從而導(dǎo)致氣體分子對壁面的切向動量適應(yīng)系數(shù)的變化.因此溫度也是納米通道流動中的一個重要變量.

本文采用三維MD模擬方法,針對處于過渡區(qū)的納米通道內(nèi)氣體剪切流動展開模擬,改變壁面剛度、壁面粗糙度以及系統(tǒng)溫度,獲得通道內(nèi)流場的密度、速度、正應(yīng)力及剪切應(yīng)力分布,以分析壁面條件和系統(tǒng)溫度對流動的影響規(guī)律,并探討近壁區(qū)域內(nèi)物理量的變化,進一步計算和分析切向動量適應(yīng)系數(shù)(TMAC)的大小和變化規(guī)律.

2 三維分子動力學(xué)模擬

2.1 物理模型及模擬細節(jié)

本文采用三維非平衡分子動力學(xué)方法來模擬微納通道中氣態(tài)氬Ar的剪切流動,其模型示意圖如圖1所示.單原子氣體Ar分布在由上下兩塊無限長平板構(gòu)成的高度為H=10.9 nm的納米通道內(nèi),X,Y方向采用周期性邊界條件,模擬區(qū)域的長寬尺寸為LX×LY=65.38 nm×65.38 nm.上下壁面以的速度沿相反方向運動,驅(qū)動氣體做剪切流動,其中M為馬赫數(shù),本文模擬時取為0.2;γ為比熱比(對于單原子氣體,γ=5/3);玻爾茲曼常數(shù)kB=1.3806×10-23J·K—1;T為氣體溫度(默認情況下T取298 K);m為氣體分子質(zhì)量.對于單原子氣體 Ar,其分子質(zhì)量mAr=6.63×10-26kg,分子直徑σAr=0.3405 nm,勢能參數(shù)εAr=1.67×10-27J.為了方便計算與分析,本文中壁面原子參數(shù)與氣體分子保持一致,即mwall=mAr,σwall=σAr,εwall=εAr.

圖1 納米通道內(nèi)氣體剪切流動示意圖Fig.1.Schematic of the shear-driven gas flow in the nanochannel.

氣體-氣體和氣體-壁面之間的相互作用采用截斷Lennard-Jones (L-J)6-12勢能函數(shù)[34]模擬,定義為

初始時刻,氣體分子均勻分布在通道內(nèi),且按照Maxwell分布產(chǎn)生初始隨機速度.速度積分采用Velocity-Verlet算法,時間步長Δ=0.002τ,系統(tǒng)運行 2×106步后達到穩(wěn)定,穩(wěn)定后繼續(xù)運行較長時間步,以確保所有宏觀量都達到熱力學(xué)穩(wěn)定狀態(tài).為了捕捉近壁區(qū)域的物理量變化,將模擬區(qū)域沿Z方向細分為若干等間距(0.2σAr)的統(tǒng)計層,以統(tǒng)計流場中各宏觀量的細致分布.

2.2 壁面剛度模型

壁面原子按照面心立方晶格點陣(FCC)排列,晶格常數(shù)為 1.6σAr,采用Einstein壁面模型[35],使壁面原子在其初始晶格點位置以相同的頻率做簡諧振動,簡諧振動由連接分子和晶格點的虛擬彈性力維持.其振動的彈性系數(shù)K表征了壁面的剛度,由Lindeman準則[36]〈取值〉,通常要避免壁面〈整體〉移動較大,限制條件為 dA2/d2＜0.023,其中dA2代表分子相對平衡位置振動時的均方位移,d為壁面中緊鄰的兩個分子之間的距離.根據(jù)壁面原子類型的不同,可以根據(jù)其對應(yīng)的Einstein溫度來確定虛擬彈簧力的彈性系數(shù)：

其中m為壁面原子的質(zhì)量,TE為固體的Einstein溫度,h為普朗克常數(shù).通過這種方法確定的彈性系數(shù)既滿足Lindemann準則,還可以模擬壁面原子熱運動及原子結(jié)構(gòu),能夠反映真實的壁面狀況.本文在默認情況下取K=1000ε/σ2,壁面剛度影響問題研究的計算條件如表1所列.

2.3 壁面粗糙度模型

本文構(gòu)造的粗糙度為一種典型的金字塔形狀的粗糙度模型[37],其幾何形狀示意圖見圖2.由于構(gòu)成粗糙壁面的原子均位于晶格點的位置上,因此該形狀的粗糙度具有合乎物理實際的傾斜度,在結(jié)構(gòu)上具有很好的穩(wěn)定性.本文設(shè)計了三種不同大小的粗糙度單元,高度分別為0.27,0.54及0.81 nm.

圖2 壁面納米粗糙度構(gòu)型Fig.2.Schematic diagram of the roughness geometry.

相比于光滑壁面,粗糙壁面的粗糙度具有不可忽視的物理尺寸,通道的特征長度會因此而發(fā)生一定的變化.圖3為粗糙通道的有效高度和光滑壁面的高度計算示意圖,在模擬過程中,為了保證各算例的克努森數(shù)保持相同,通道的高度必須一致.對于不同的粗糙壁面,統(tǒng)一地將粗糙通道的上邊界和下邊界分別放置在下壁面和上壁面粗糙度的中心位置.由于氣體分子會進入粗糙度內(nèi)部的間隙中,而在統(tǒng)計流場中的物理信息分布時,這部分氣體分子也是需要考慮的,因此,就會導(dǎo)致在后續(xù)的圖表中出現(xiàn)粗糙通道的真實高度(Z)為負值及約化高度(Z/H)小于—0.5的情況.粗糙度影響問題研究的計算條件如表1所列.

表1 計算條件設(shè)置Table 1.Parameters of the study on the influences of surface stiffness,roughness and temperature.

2.4 溫度控制方法

所有的系統(tǒng)都會與外界發(fā)生相互作用,比如受到外界溫度和壓力的影響.進行MD模擬時通常需要構(gòu)造一個等溫或等壓的系綜,一般來說,模擬系統(tǒng)的分子數(shù)目和體積比較容易保持恒定,而溫度和壓力則不同.由于系統(tǒng)中的粒子不停地運動與碰撞,并伴隨著動量、能量的交換,因此,溫度和壓力是不斷變化的.為了使系統(tǒng)溫度或壓力保持恒定,有必要對溫度及壓力進行適當?shù)目刂坪驼{(diào)節(jié),也稱作熱浴.本文采用Nose-Hoover熱浴方法[38],則模擬系統(tǒng)的運動方程可寫作：

式中F為粒子i所受到的作用力;η為熱浴阻尼系數(shù),

合作式治理路徑主張社區(qū)既非國家壟斷公共權(quán)力的領(lǐng)域，也非完全由社會自治的場所，而是國家權(quán)力和社會權(quán)力共同作用、合作治理的領(lǐng)域。研究者認為，社區(qū)治理是最接近于“治理”原義的治理形態(tài)，它是由多元主體共同參與、管理、決策和協(xié)商的合作治理結(jié)構(gòu)，具有基層草根性、利益相關(guān)性與主體多元性的特征。⑤中國的社區(qū)治理要完成從威權(quán)式治理向參與式治理體制轉(zhuǎn)型，賦予居民參與社區(qū)公共事務(wù)的權(quán)力，讓居民自己決定關(guān)系切身利益的社區(qū)公共事務(wù)，形成以社區(qū)需求為中心的權(quán)力和資源分配體制，發(fā)展社區(qū)居民自我治理的理念和能力，讓社區(qū)發(fā)展主體——居民能夠有效地介入到社區(qū)建設(shè)的決策和管理過程中。⑥

其中T0為目標 控 制 溫度,T(t)為當前時刻的 溫 度,通過阻尼系數(shù)使流場溫度保持在目標溫度下.本文在默認情況下溫度取為298 K,溫度影響問題研究的計算條件如表1所列.

3 結(jié)果與討論

3.1 壁面彈性系數(shù)對流動的影響分析

為了研究壁面彈性系數(shù)對納米通道內(nèi)氣體流動的影響規(guī)律,本節(jié)設(shè)置了一組算例,彈簧系數(shù)設(shè)置在 3 00ε/σ2—1 000ε/σ2之間,作為對比,同時計算了壁面原子固定的情況,也即彈性系數(shù)無窮大.

圖4為納米通道內(nèi)的速度分布曲線,由圖可知,壁面彈性系數(shù)對主流區(qū)的速度分布影響不大,但對近壁面區(qū)域的速度分布有一定影響.隨著壁面彈性系數(shù)的增大,氣體分子在近壁面區(qū)域的速度峰值 分 別 為 ：—51.17,—50.08,—49.72,—49.23,及—47.26 m/s.從原子運動的角度出發(fā),當氣體分子與壁面相互作用時,壁面原子的振動頻率越大,則氣體分子就越不容易在壁面停留,從而導(dǎo)致氣體分子與壁面的動量適應(yīng)性變差,氣體分子在近壁面區(qū)域的速度峰值隨之減小.

圖4 不同壁面彈性系數(shù)條件下流場速度分布 (a)無量綱速度分布;(b)近壁面區(qū)域速度分布Fig.4.Velocity profiles with different surface stiffness：(a)Normalized velocity profiles;(b)near-wall velocity profiles.

圖5—圖7分別為不同壁面彈性系數(shù)條件下,納米通道中的密度、正應(yīng)力及剪切應(yīng)力分布.由圖可知,壁面彈性系數(shù)對流場內(nèi)上述物理量的分布影響不是很明顯.不同壁面彈性系數(shù)條件下,流場中心區(qū)域的密度和剪切應(yīng)力如表2所列,通過剪切應(yīng)力與克努森數(shù)的關(guān)系[39,40],得到不同壁面彈性系數(shù)下的TMAC分別為0.740,0.736,0.731,0.731,0.726,總體上呈減小趨勢,但影響較小.

3.2 壁面粗糙度對流動的影響分析

為了研究粗糙度對納米尺度氣體流動的影響規(guī)律,改變通道壁面的粗糙高度,模擬了不同粗糙度大小情況下的納米尺度氣體流動,統(tǒng)計獲得了流場內(nèi)速度、密度、正應(yīng)力及剪切應(yīng)力分布.

圖5 不同壁面彈性系數(shù)條件下納米通道中的密度分布Fig.5.Density distribution with different surface stiffness.

圖6 不同壁面彈性系數(shù)條件下納米通道中的正應(yīng)力分布Fig.6.Normal stress distribution with different surface stiffness.

圖7 不同壁面彈性系數(shù)條件下納米通道中的剪切應(yīng)力分布Fig.7.Shear stress distribution with different surface stiffness.

圖8為不同粗糙通道內(nèi)氣體的無量綱速度分布,可知無論是主流區(qū)域還是近壁面區(qū)域,粗糙度對流動的影響都非常明顯.而且由于粗糙度的存在,在壁面出現(xiàn)了一個粗糙度影響區(qū)域,在這個區(qū)域內(nèi),也存在著一定的速度分布規(guī)律.為了更加細致地了解這一區(qū)域的速度分布情況,截取了距下壁面邊界2.0 nm以下范圍內(nèi)的無量綱速度分布曲線,如圖9所示.從圖9可以看出,在主流區(qū)域,氣體的速度隨著壁面粗糙度的增大而增大,當粗糙度高度為0.54和0.81 nm時,主流區(qū)域的速度分布趨于一致.此外,由于壁面粗糙度的增大,主流區(qū)的線性分布特性也受到了一定的影響,線性區(qū)域的范圍減小.在近壁面區(qū)域,隨著粗糙度的增大,流體速度迅速增大并接近于壁面速度的大小,當粗糙度足夠大(0.81 nm)時,幾乎整個粗糙單元縫隙中的氣體流速都與壁面速度相等.

表2 不同壁面彈性系數(shù)條件下納米通道中的剪切應(yīng)力及TMAC大小Table 2.Shear stress and TMAC values with different surface stiffness.

圖8 不同壁面粗糙度下的流場無量綱速度分布Fig.8.Normalized velocity profiles with different surface roughness.

圖9 不同壁面粗糙度下近壁區(qū)域流場無量綱速度分布Fig.9.Near-wall normalized velocity profiles with different surface roughness.

圖10為不同粗糙度高度下近壁面區(qū)域的密度分布曲線,可知隨著壁面粗糙度的增大,氣體的在近壁面附近的峰值隨之增大,同時粗糙度縫隙內(nèi)的氣體密度也比主流區(qū)域的密度要大,這說明氣體分子在粗糙壁面停留的時間較長,氣體分子在粗糙壁面的集聚現(xiàn)象明顯.圖11為不同粗糙度高度下近壁面區(qū)域的正應(yīng)力分布曲線,由于正應(yīng)力的分布和氣體密度直接相關(guān),因此正應(yīng)力分布曲線呈現(xiàn)出和密度分布相似的規(guī)律.

圖10 不同壁面粗糙度下近壁區(qū)域流場密度分布Fig.10.Near-wall density distribution with different surface roughness.

圖12為不同大小的粗糙度情形下剪切應(yīng)力分布曲線.在主流區(qū)域剪切應(yīng)力的大小趨于恒定,且隨著粗糙度高度的增大而增大.當粗糙度高度增大到一定程度時,剪切應(yīng)力不再繼續(xù)增大,比如當粗糙度高度為0.54和0.81 nm時,剪切應(yīng)力分別為18.85和18.99 kPa,二者幾乎相等.通過剪切應(yīng)力與切向動量適應(yīng)系數(shù)的關(guān)系,計算得到氣體分子對不同粗糙壁面的切向動量適應(yīng)系數(shù)分別為0.87,0.98,0.99,可以看到當粗糙度大于0.54 nm時,TMAC的數(shù)值已經(jīng)接近于1.0,表明氣體分子的切向動量與粗糙壁面完全適應(yīng).上述結(jié)果說明粗糙度越大,氣體與壁面的動量交換程度越高,氣體越容易在壁面發(fā)生漫反射.

圖11 不同壁面粗糙度下近壁區(qū)域流場正應(yīng)力分布Fig.11.Near-wall normal stress distribution with different surface roughness.

圖12 不同壁面粗糙度下近壁區(qū)域流場剪切應(yīng)力分布Fig.12.Near-wall shear stress distribution with different surface roughness.

3.3 系統(tǒng)溫度對流動的影響分析

為了研究系統(tǒng)溫度對納米通道內(nèi)氣體流動特性及切向動量適應(yīng)系數(shù)的影響規(guī)律,針對通道高度為10.90 nm的氣體剪切流動展開模擬,所有算例中的氣體分子數(shù)目都保持一致,待系統(tǒng)平衡后,統(tǒng)計獲得流場內(nèi)速度、密度、正應(yīng)力及剪切應(yīng)力分布.

圖13為不同溫度條件下流場的無量綱速度分布,可知溫度的變化明顯地改變了流場的速度分布.為了更好地觀察溫度對流場速度的影響,截取了距離壁面2.0 nm的區(qū)域內(nèi)的速度分布,如圖14所示.從圖14可以清晰地看到,整個流場的速度都隨著溫度的增大而減小.近壁面區(qū)域速度峰值大小的變化直觀印證了這一點,由248 K時的51.54 m/s降低到了373 K時的41.03 m/s.此外,主流區(qū)域的速度剪切率也隨著溫度的增大而降低.

圖13 不同溫度條件下的流場無量綱速度分布Fig.13.Normalized velocity profiles with different temperature conditions.

圖14 不同溫度條件下近壁區(qū)域流場速度分布Fig.14.Near-wall velocity profiles with different temperature conditions.

圖15為不同溫度條件下的流場密度分布,由于主流區(qū)域的密度趨于定值,本文只展示了距壁面2.0 nm范圍內(nèi)的密度分布,后面的正應(yīng)力分布及剪切應(yīng)力也采用了這種表示方法.由圖15可知,在壁面作用力影響的近壁面區(qū)域,氣體密度隨著溫度的降低而升高,由于在初始設(shè)置計算條件時,所有算例都布置了相同數(shù)目的氣體分子,因此,在溫度較低時,氣體分子在壁面的聚集會導(dǎo)致流場中心區(qū)域的密度有所降低.圖16為流場中的正應(yīng)力分布,正應(yīng)力的分布主要和氣體溫度相關(guān),流場的溫度越高,則氣體的正應(yīng)力越大,MD的計算結(jié)果也呈現(xiàn)了這種趨勢,即流場的正應(yīng)力隨著溫度的升高而升高.將MD計算的主流區(qū)域的正應(yīng)力結(jié)果和理想氣體理論預(yù)測的結(jié)果列于表3,通過對比發(fā)現(xiàn),對于所有溫度條件下MD計算的結(jié)果和理論預(yù)測的誤差不超過2%,以上結(jié)果也直接驗證了MD在計算宏觀量時的準確性.

圖15 不同溫度條件下近壁區(qū)域流場密度分布Fig.15.Near-wall density distribution with different temperature conditions.

圖16 不同溫度條件下近壁區(qū)域流場正應(yīng)力分布Fig.16.Near-wall normal stress distribution with different temperature conditions.

表3 不同溫度條件下主流區(qū)域正應(yīng)力與理論值對比Table 3.Comparisons of normal stress with its theoretic values in bulk flow under different temperature conditions.

圖17為不同溫度條件下的剪切應(yīng)力的分布,流場內(nèi)的剪切應(yīng)力在總體上是隨著溫度的升高而增大的.由氣體動理論可知,剪切應(yīng)力的變化不僅受到系統(tǒng)溫度的影響,它還受到氣體分子對壁面切向動量適應(yīng)系數(shù)的影響[39].計算得到的切向動量適應(yīng)系數(shù)如圖18所示,隨著溫度的升高,氣體分子對壁面的動量適應(yīng)隨之降低.溫度除了對TMAC有影響,對氣體分子在納米通道內(nèi)有效黏性系數(shù)也有明顯的影響,由氣體分子運動論可知,氣體溫度越高,則氣體分子間的碰撞越頻繁,則氣體的黏性越大.氣體的有效黏性系數(shù)可以通過下式計算：

圖17 不同溫度條件下近壁區(qū)域流場剪切應(yīng)力分布Fig.17.Near-wall normal stress distribution with different temperature conditions.

圖18 TMAC隨溫度的變化Fig.18.TMAC values with different temperature conditions.

其中τshear為流場中心區(qū)域剪切應(yīng)力,?u/?z為流場中心區(qū)域的速度剪切率.根據(jù)(5)式計算得到了不同溫度下氣體有效黏性系數(shù)如圖19所示,氣體分子的有效黏性系數(shù)隨著溫度的升高而增大.

圖19 有效黏性系數(shù)隨溫度的變化Fig.19.Effective viscosity values at different temperature.

4 結(jié) 論

采用三維非平衡分子動力學(xué)模擬方法研究了過渡區(qū)納米通道內(nèi)的氣體剪切流動,獲得了通道內(nèi)氣體的各物理量的分布,深入研究了近壁區(qū)域的流動特性,探討了壁面條件及系統(tǒng)溫度對納米尺度氣體流動的影響規(guī)律,并進一步分析了切向動量適應(yīng)系數(shù)的變化趨勢,發(fā)現(xiàn)了以下規(guī)律：

1)壁面原子的振幅相對于氣體分子平均自由程來說是小量,壁面彈性系數(shù)的變化對流動的影響不明顯,但是距壁面約一個分子直徑處的速度峰值以及切向動量適應(yīng)系數(shù)會隨著壁面彈性系數(shù)的增大而減小;

2)壁面粗糙度對流動的影響主要通過切向動量適應(yīng)系數(shù)的變化來體現(xiàn),隨著粗糙度的增大,氣體分子對壁面的切向動量適應(yīng)程度逐漸增強,直至完全適應(yīng);由于壁面粗糙度的存在會使壁面力場發(fā)生改變,進而造成近壁面區(qū)域速度、密度及應(yīng)力分布的波動;

3)系統(tǒng)溫度的變化對流動的影響是全局性的,無論是在近壁面區(qū)域還是在主流區(qū)域,氣體的速度、密度、正應(yīng)力和剪切應(yīng)力都隨著溫度變化而改變,溫度變化還會引起宏觀參數(shù)的改變,氣體分子的動量適應(yīng)系數(shù)隨著溫度的升高而降低,有效黏性系數(shù)隨著溫度的升高而增大;

4)對于納米尺度條件下的氣體流動,由于表面積/體積之比極大,壁面力場對流場的主導(dǎo)作用更加明顯,流動物理量對于壁面條件和流場溫度的改變也十分敏感;壁面剛度的改變雖然對主流區(qū)域沒有影響,但是對近壁區(qū)域的流動仍然有一定程度的干擾;而壁面粗糙度無論對主流區(qū)域還是近壁區(qū)域的流動都有較大影響,且氣體分子在壁面處的聚集現(xiàn)象明顯;系統(tǒng)溫度的影響則是全局性的,可以顯著改變氣體與氣體、氣體與壁面之間的動量交換程度,因此對主流區(qū)域和近壁面區(qū)域的影響也都較大.