亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        宏觀速度對(duì)適應(yīng)系數(shù)的影響規(guī)律研究*

        2019-11-28 08:37:06張燁張冉2賴劍奇李樺
        物理學(xué)報(bào) 2019年22期
        關(guān)鍵詞:切向速度法向動(dòng)量

        張燁 張冉2)? 賴劍奇 李樺

        1) (國(guó)防科技大學(xué)空天科學(xué)學(xué)院,長(zhǎng)沙 410073)

        2) (西北核技術(shù)研究院,激光與物質(zhì)相互作用國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710024)

        適應(yīng)系數(shù)是稀薄氣體動(dòng)力學(xué)中的重要參數(shù),反映了氣體分子與固體表面的動(dòng)量和能量適應(yīng)程度,常常用于稀薄氣體數(shù)值模擬的邊界條件中.本文構(gòu)建了單個(gè)氣體分子Ar與金屬Pt表面相互作用的物理模型,通過Phontom模型來構(gòu)造固體壁面以反映真實(shí)物理特性.采用分子動(dòng)力學(xué)方法模擬了氣體分子在固體表面的碰撞過程,通過速度抽樣使得入射氣體具有宏觀速度特征,根據(jù)入射前和反射后氣體分子的平均動(dòng)量和能量計(jì)算出切向動(dòng)量適應(yīng)系數(shù)、法向動(dòng)量適應(yīng)系數(shù)和能量適應(yīng)系數(shù),分別在光滑和粗糙表面下分析了適應(yīng)系數(shù)隨宏觀切向速度和宏觀法向速度的變化規(guī)律.結(jié)果顯示,適應(yīng)系數(shù)對(duì)于氣體的宏觀速度和表面粗糙度均表現(xiàn)出了很強(qiáng)的相關(guān)性.切向速度的增大使得氣體分子與表面的作用時(shí)間縮短,適應(yīng)系數(shù)減小,對(duì)于氣固適應(yīng)有消極的作用.當(dāng)切向速度較大時(shí)會(huì)發(fā)生動(dòng)量由切向向法向轉(zhuǎn)移的情況.法向速度的增大對(duì)于光滑表面而言有助于氣體分子與表面的動(dòng)量和能量適應(yīng),對(duì)于粗糙表面而言卻呈現(xiàn)相反的規(guī)律.

        1 引 言

        由于稀薄氣體流動(dòng)在微納機(jī)電系統(tǒng)、真空技術(shù)、航天工程以及過程工藝比如多相催化、氣體分離薄膜等方面的廣泛應(yīng)用,有關(guān)稀薄氣體流動(dòng)機(jī)理的研究備受重視[1-5].氣體與固體表面的相互作用作為稀薄氣體流動(dòng)的關(guān)鍵問題,可以為物理吸附狀態(tài)、能量轉(zhuǎn)換機(jī)制以及分子碰撞的細(xì)節(jié)等稀薄氣體流動(dòng)邊界問題提供極有價(jià)值的信息[6].然而由于物理上的復(fù)雜性,目前對(duì)氣-固相互作用規(guī)律的研究還不充分,因此有必要對(duì)其展開深入研究.

        氣體的稀薄程度由克努森數(shù) (K n=λ/L,其中 λ 為氣體分子平均自由程,L為流動(dòng)特征長(zhǎng)度)來量化.隨著克努森數(shù)的增大,氣體稀薄程度增加,此時(shí)連續(xù)性假設(shè)不再適用,必須使用動(dòng)能理論來描述分子的輸運(yùn).許多基于玻爾茲曼方程的粒子類數(shù)值方法,比如直接模擬蒙特卡羅法(DSMC)、網(wǎng)格粒子法(PIC)以及格子玻爾茲曼法(LBM),可以很好地模擬稀薄氣體流動(dòng)[7-9].粒子類方法具備合適的處理分子間碰撞的模型,能夠模擬復(fù)雜的稀薄流動(dòng)系統(tǒng).然而,由于稀薄條件下氣體分子間碰撞頻率的減小,氣體與固體表面的相互作用則起主導(dǎo)作用[10].為了模擬氣體在固體邊界處的散射現(xiàn)象,粒子類方法中通常需引入散射核模型來描述氣體在表面的散射規(guī)律.Maxwell提出了最早的散射核模型[11],他認(rèn)為氣體分子與表面碰撞后,α 部分氣體分子在表面發(fā)生漫反射,即該部分氣體與表面完全適應(yīng),剩余的 1 -α 部分氣體分子在表面發(fā)生鏡面反射,反射氣體分子仍然保留了入射氣體分子的信息.另一個(gè)廣泛使用的散射核模型為CLL模型[12,13],由Cercignani和Lampis最先提出并由Lord進(jìn)一步發(fā)展.CLL模型是一種滿足歸一化條件和互易性原理并且恒為正的散射核模型,包含兩個(gè)獨(dú)立參數(shù),用于反映不同方向的散射特性,因此CLL模型在描述氣體散射規(guī)律方面比Maxwell模型更為靈活.上述散射核模型中使用的參數(shù)稱為適應(yīng)系數(shù),適應(yīng)系數(shù)包括三種不同類型,分別為切向動(dòng)量適應(yīng)系數(shù)(TMAC),法向動(dòng)量適應(yīng)系數(shù)(NMAC)以及能量適應(yīng)系數(shù)(EAC),分別表征氣體分子的切向動(dòng)量、法向動(dòng)量以及能量對(duì)表面的適應(yīng)程度.在稀薄氣體動(dòng)力學(xué)相關(guān)的各種模擬方法中,適應(yīng)系數(shù)是邊界條件的核心參數(shù)[14,15].不僅在散射核模型中適應(yīng)系數(shù)必須給定,在通過連續(xù)介質(zhì)方程求解滑移流動(dòng)時(shí)附加的滑移邊界條件和溫度階躍條件中,TMAC和EAC也是必須用到的參數(shù)[16,17].由此可見,適應(yīng)系數(shù)在稀薄氣體流動(dòng)的邊界條件中扮演著重要的角色,是描述氣-固相互作用的核心參數(shù),能否準(zhǔn)確地模擬稀薄氣體流動(dòng),很大程度上取決于適應(yīng)系數(shù)的取值能否符合氣體在表面的散射規(guī)律.然而,由于適應(yīng)系數(shù)的影響規(guī)律眾多,針對(duì)不同的問題,適應(yīng)系數(shù)的值難以確定,很多工程問題中只是經(jīng)驗(yàn)性地將適應(yīng)系數(shù)取值為1,這必然會(huì)導(dǎo)致流場(chǎng)模擬結(jié)果的不精確.為了構(gòu)造更加符合物理實(shí)際的邊界條件,使得稀薄氣體流動(dòng)的數(shù)值模擬更加精確,必須從氣體分子與表面相互作用的微觀過程展開研究,充分了解氣體分子在與表面碰撞過程中的動(dòng)量和能量交換機(jī)理,理清適應(yīng)系數(shù)與來流條件和固體表面條件之間的復(fù)雜關(guān)系.

        分子動(dòng)力學(xué)(MD)方法[18]是研究氣-固相互作用的有力手段[19,20],不引入任何假設(shè),許多學(xué)者都證實(shí)了分子動(dòng)力學(xué)的結(jié)果能夠很好地符合已有的實(shí)驗(yàn)結(jié)果[21,22].此外,分子動(dòng)力學(xué)方法還以可以對(duì)一些實(shí)驗(yàn)所不能實(shí)現(xiàn)的極端流場(chǎng)環(huán)境下的流動(dòng)機(jī)理開展研究[23].通過分子動(dòng)力學(xué)方法,前人開展了許多氣-固相互作用以及適應(yīng)系數(shù)方面的研究.Cao等[24]采用MD方法研究了微通道內(nèi)的氣體流動(dòng),發(fā)現(xiàn)TMAC隨著溫度升高的遞減規(guī)律符合指數(shù)衰減率,且在較低的溫度下,TMAC對(duì)于溫度的敏感性更大.Finger等[25]研究了He氣體分子在金屬銅表面的散射特性,研究了氣體入射角度、入射能量、表面晶體構(gòu)型以及表面吸附層對(duì)氣體分子散射特性和切向動(dòng)量輸運(yùn)的影響.Sun和Li[26-28]采用二維MD方法,在流動(dòng)過程中對(duì)氣體分子與表面的碰撞進(jìn)行跟蹤,統(tǒng)計(jì)得到氣體分子的入射及反射速度,從而直接計(jì)算得到氣體分子對(duì)表面的適應(yīng)系數(shù),并進(jìn)一步研究了表面構(gòu)型、表面溫度和物理吸附等因素對(duì)適應(yīng)系數(shù)的影響規(guī)律.Kovalev和Yakunchikov[22]通過MD方法研究了氫氣在石墨表面的散射過程,發(fā)現(xiàn)當(dāng)氣體分子溫度升高時(shí),表面溫度的影響會(huì)減弱,當(dāng)溫度較低時(shí),由于表面的物理吸附,適應(yīng)系數(shù)會(huì)增大.Prabha和Sathian[29]基于MD方法模擬光滑納米通道內(nèi)的氣體流動(dòng)來研究氣-固相互作用特性,發(fā)現(xiàn)EAC隨表面溫差的變化呈現(xiàn)出自然對(duì)數(shù)的關(guān)系,動(dòng)量的切向分量與分子入射速度表現(xiàn)出很強(qiáng)的相關(guān)性.Pham等[30]采用Kulginov勢(shì)能函數(shù)描述Ar原子與金屬Pt之間的勢(shì)能相互作用,進(jìn)一步研究了表面構(gòu)型、入射角度以及表面溫度對(duì)TMAC的影響規(guī)律.Liang和Keblinski[31]采用大規(guī)模MD方法模擬了單原子和雙原子氣體在Au表面的流動(dòng),研究了氣固相互作用強(qiáng)度、氣固原子質(zhì)量比、固體彈性剛度對(duì)熱適應(yīng)系數(shù)和動(dòng)量適應(yīng)系數(shù)的影響規(guī)律.Reinhold等[32]構(gòu)建了真實(shí)氣體分子Ar,N2和CO2在真實(shí)表面的散射系統(tǒng),研究了散射后的角度分布和TMAC的變化規(guī)律.Lim等[33]采用MD方法研究了He,Ne,Ar,Kr和Xe等惰性氣體的流動(dòng),發(fā)現(xiàn)隨著氣體分子質(zhì)量的增大,TMAC增大或減小的情況均會(huì)出現(xiàn).Yakunchikov等[34]通過惰性氣體在石墨表面的散射研究了氣體分子質(zhì)量、大小以及氣-固作用勢(shì)能對(duì)TMAC和EAC的影響規(guī)律.張冉等[35]模擬了過渡區(qū)納米通道內(nèi)剪切力驅(qū)動(dòng)的氣體流動(dòng),探究了通道高度和氣-固相互作用勢(shì)能系數(shù)對(duì)TMAC的影響規(guī)律.

        總結(jié)前人的工作,我們發(fā)現(xiàn),氣體分子在固體表面的適應(yīng)系數(shù)隨溫度、勢(shì)能系數(shù)、表面構(gòu)型、分子質(zhì)量、分子類型、入射角度等眾多因素影響,展現(xiàn)出了豐富的變化規(guī)律.然而對(duì)于具有一定宏觀速度的來流條件,其對(duì)適應(yīng)系數(shù)的影響規(guī)律尚未見系統(tǒng)深入的研究.事實(shí)上,在臨近空間稀薄大氣環(huán)境下,自由來流作用在航天飛行器的表面時(shí)不但有很大的宏觀速度,而且會(huì)與表面形成一定的夾角,導(dǎo)致法向速度分量的影響也不容忽視.而且真實(shí)情況下的固體表面總是具有一定的粗糙度[36,37],且粗糙度的尺寸相比于表面力場(chǎng)勢(shì)阱深度并不是小量,粗糙度的存在會(huì)大大改變表面的勢(shì)能分布[38],因此本文還探究了宏觀速度與粗糙度的耦合效應(yīng)的影響.本文構(gòu)建的物理模型為單個(gè)分子與表面的相互作用,采用該模型只關(guān)注氣體分子與表面的相互作用,避免計(jì)算主流區(qū)的流動(dòng)時(shí)占據(jù)大量計(jì)算時(shí)間,大大提高了計(jì)算效率,且氣體分子的入射條件可以靈活控制,氣-固相互作用信息易于獲得,方便對(duì)相關(guān)問題的研究.

        本文構(gòu)建了單分子氬氣與金屬鉑表面的相互作用模型,基于MD方法模擬了氣體分子在固體表面的散射過程,通過速度取樣法統(tǒng)計(jì)出氣體分子入射前和反射后的平均動(dòng)量和能量,計(jì)算出氣體分子的TMAC,NMAC以及EAC,進(jìn)一步分析了在光滑和粗糙表面下,來流宏觀切向和法向速度對(duì)氣體在固體表面的動(dòng)量和能量交換的影響規(guī)律.

        2 氣體-表面相互作用模型

        本文采用分子動(dòng)力學(xué)方法研究氣體分子Ar在金屬Pt表面的散射過程,其物理模型如圖1所示.固體壁面采用Phantom模型[39]來構(gòu)造,該模型在模擬壁面原子結(jié)構(gòu)、處理壁面原子與氣體分子間的能量交換等方面都比較真實(shí).該模型包含三層壁面原子: 第一層為真實(shí)原子層,與氣體分子直接發(fā)生相互作用,該層依靠虛擬彈簧力和勢(shì)能函數(shù)模擬壁面原子的運(yùn)動(dòng);第二層為緩沖層,將真實(shí)原子層與氣體分子相互作用的力、動(dòng)量、能量向基底原子傳遞,其原子的運(yùn)動(dòng)由連接第一層和第三層的虛擬彈簧力維持;第三層為基底,原子位置固定以穩(wěn)定壁面.

        圖1 氣體分子在固體表面散射的物理模型Fig.1.Physical model for scattering gas molecules on solid surfaces.

        壁面原子按照面心立方晶格點(diǎn)陣(FCC)排列,其(100)晶面與氣體分子相接觸,晶格常數(shù)A=1.15σ,X和Y方向上采用周期性邊界條件.為了考察粗糙度對(duì)氣體與表面相互作用的影響,本文通過一種典型的金字塔模型來構(gòu)造表面粗糙度[40],其幾何形狀示意圖見圖2.根據(jù)能量最低原理,塊狀粗糙度由于具有較大的表面自由能而不太穩(wěn)定,其粗糙單元的上層原子會(huì)逐漸遷移到較低位置,最終會(huì)演變成具有較小表面能的金字塔形結(jié)構(gòu).因此本文采取的粗糙度模型合乎物理實(shí)際,在結(jié)構(gòu)上具有很好的穩(wěn)定性.本文采用兩種不同大小的粗糙度單元,高度分別為 0.5 A 和 1.0 A.

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

        氣體分子與壁面原子之間的相互作用采用截?cái)郘ennard-Jones(L-J)6-12勢(shì)能函數(shù)[41]來模擬,其函數(shù)表達(dá)式為

        表1 分子動(dòng)力學(xué)模擬中的參數(shù)Table 1.Parametersin the MD simulations.

        分子間的作用力通過鄰位列表法來處理,運(yùn)動(dòng)方程通過velocity-Verlet算法來積分.在計(jì)算過程中,為了使壁面始終保持在給定的溫度下,我們對(duì)壁面原子進(jìn)行速度標(biāo)定,其公式如下:

        其中TC為期望溫度,本文取值為 300 K.計(jì)算時(shí)間步長(zhǎng)取為0.0005,前5萬(wàn)時(shí)間步為壁面弛豫時(shí)間.由于氣體的散射在各個(gè)切向方向是各向同性的[22],因此不妨將氣體分子始終沿著x軸正方向入射.每次模擬的初始時(shí)刻,氬原子處于坐標(biāo)(x,y,z=3.0σArPt)的位置(坐標(biāo)原點(diǎn)在z方向選取為表面最上層原子的位置),x和y是均勻隨機(jī)分布的變量.入射氣體分子的初始高度與表面最上層原子的距離剛好等于氣-固相互作用截?cái)喟霃絩c.MD模擬的分子運(yùn)動(dòng)通過短程力來控制,由(1)式可以看出,若氣體分子與表面的距離大于等于截?cái)喟霃絩c,氣-固相互作用力趨于零,可忽略不計(jì).只有當(dāng)氣體分子與壁面的距離小于截?cái)喟霃綍r(shí),氣-固相互作用才變得重要.因此這一入射高度值的選取保證了初始時(shí)刻氣體分子不受表面力場(chǎng)的作用,同時(shí)也避免了多余的計(jì)算時(shí)間步.

        單分子入射分子動(dòng)力學(xué)包含兩種不同的實(shí)施途徑[32].第一種為分子束法,受分子束實(shí)驗(yàn)啟發(fā),入射氣體分子的能量和角度固定以減小系統(tǒng)的復(fù)雜性,用于細(xì)致研究入射角和反射角的關(guān)系.第二種為采樣法,入射氣體分子的速度符合特定溫度下的Maxwell-Boltzmann分布,因?yàn)榭梢阅M自由來流條件下的氣-固相互作用問題.本文采用采樣法,入射氣體分子的初始速度可以表達(dá)為宏觀速度與符合Maxwell-Boltzmann分布的熱運(yùn)動(dòng)速度的疊加:

        式中,u′,w′分別為切向和法向熱運(yùn)動(dòng)速度;u0,w0分別為切向和法向宏觀速度;T為氣體溫度;kB=1.3806×10-23JK-1為玻爾茲曼常數(shù).氣體與壁面可能發(fā)生多次碰撞,判斷散射過程結(jié)束的標(biāo)志為氣體分子的高度大于初始時(shí)刻的高度值,因?yàn)榇藭r(shí)氣-固相互作用力消失,氣體分子不再受到表面力場(chǎng)的作用,從而離開表面區(qū)域,標(biāo)志著散射過程的結(jié)束.為了保證計(jì)算結(jié)果的準(zhǔn)確性,對(duì)于每個(gè)來流條件,均計(jì)算了30萬(wàn)次氣體分子在表面上的散射.

        氣體與表面的相互作用的切向動(dòng)量適應(yīng)系數(shù)σt(TMAC)、法向動(dòng)量適應(yīng)系數(shù) σn(NMAC)以及能量適應(yīng)系數(shù) α (EAC)的定義式分別為:

        其中,τ 為氣體分子的平均切向動(dòng)量,p為氣體分子的平均法向動(dòng)量,q為氣體分子的平均能量,下標(biāo)i,r和w分別代表氣體分子入射前、反射后以及與表面完全適應(yīng)后的狀態(tài).通過對(duì)與表面完全適應(yīng)后的反射氣體分子的速度分布函數(shù)進(jìn)行積分,可得

        如果氣體分子與表面的作用時(shí)間超過30 ps,則代表氣體分子在表面發(fā)生了吸附,其解吸附后的動(dòng)量遵從漫反射規(guī)律,即發(fā)生了完全的動(dòng)量適應(yīng),這種情況下的適應(yīng)系數(shù)均為1,所以應(yīng)當(dāng)對(duì)上述計(jì)算得到的適應(yīng)系數(shù)進(jìn)行吸附修正,即

        式中AC為三種適應(yīng)系數(shù)的統(tǒng)稱,ξ 表示氣體分子在表面的吸附概率.

        3 結(jié)果與討論

        3.1 光滑表面下宏觀速度的影響分析

        3.1.1 宏觀切向速度對(duì)適應(yīng)系數(shù)的影響規(guī)律

        首先,將法向宏觀速度分量設(shè)為0,僅改變切向宏觀速度分量的值,研究其對(duì)氣體分子表面適應(yīng)的影響規(guī)律.由于宏觀法向速度為0,則有 pi=pw,導(dǎo)致在計(jì)算NMAC時(shí)出現(xiàn)奇異值,且系統(tǒng)為溫度保持在300 K的等溫流動(dòng),因此,在這組算例中,只計(jì)算了TMAC.

        圖3為不同宏觀切向速度下的TMAC,由圖可知,TMAC隨著切向速度的增大而減小,TMAC的數(shù)值從 Vt=0.2σ/τ 時(shí) 的0.37降低到Vt=4.0σ/τ時(shí)的0.11,可見來流宏觀切向速度的增大對(duì)氣體分子與表面的動(dòng)量適應(yīng)具有明顯消極的影響,來流宏觀切向速度越大,則氣體對(duì)表面的切向動(dòng)量適應(yīng)情況越差.為了考察其中的原因,本文統(tǒng)計(jì)了不同來流切向速度大小下、氣體分子與表面的相互作用時(shí)間,如圖4所示.從圖中可知,由于來流切向速度的增大,氣體分子與表面作用時(shí)間隨之減小,進(jìn)而導(dǎo)致氣體分子與表面的適應(yīng)程度減弱,氣體分子仍能在一定程度上保持其原有的宏觀切向速度特征.

        圖3 不同宏觀切向速度下的切向動(dòng)量適應(yīng)系數(shù)Fig.3.TMAC values of gas molecular under different macroscopic tangential velocities.

        圖4 不同宏觀切向速度下氣體分子與表面作用時(shí)間Fig.4.Time of gas-solid interaction under different macroscopic tangential velocities.

        3.1.2 宏觀法向速度對(duì)適應(yīng)系數(shù)的影響規(guī)律

        為了考察法向速度對(duì)適應(yīng)系數(shù)的影響,在本小節(jié)中設(shè)置 4組算例,氣體的宏觀切向速度和3.1.1小節(jié)一樣,在 0 . 2σ/τ 到 4 . 0σ/τ 之間變化,4組算例的宏觀法向速度分別為為 0 . 5σ/τ,1 . 0σ/τ,1.5σ/τ,2 . 0σ/τ.通過計(jì)算,我們得到了不同法向速度和切向速度下的TMAC,NMAC和EAC,分別如圖5、圖6和圖7所示.

        圖5 不同宏觀切向速度及法向速度下的切向動(dòng)量適應(yīng)系數(shù)Fig.5.TMAC values of gas molecular under different macroscopic tangential and normal velocities.

        圖6 不同宏觀切向速度及法向速度下的法向動(dòng)量適應(yīng)系數(shù)Fig.6.NMAC values of gas molecular under different macroscopic tangential and normal velocities.

        圖7 不同宏觀切向及法向速度下的法向能量適應(yīng)系數(shù)Fig.7.EAC values of gas molecular under different macroscopic tangential and normal velocities.

        由圖5可知,在不同的法向速度下,TMAC隨切向速度的變化趨勢(shì)一致,切向速度的增大均會(huì)導(dǎo)致 TMAC 的減小.在切向速度小于 2 . 0σ/τ 時(shí),法向速度的變化對(duì)氣體分子TMAC的影響很小.當(dāng)切向速度大于 2 . 0σ/τ 后,法向速度的增大對(duì)氣體分子切向動(dòng)量與表面適應(yīng)的促進(jìn)作用開始顯現(xiàn)出來.總的來說,法向速度的變化對(duì)TMAC影響較小.當(dāng)來流切向速度較小時(shí),法向速度的增大對(duì)氣體分子與表面的切向動(dòng)量沒有太大影響,當(dāng)切向速度比較大時(shí),法向速度的增大對(duì)氣體分子與表面的切向動(dòng)量適應(yīng)有積極作用.

        由圖6可知,在不同的宏觀法向速度下,NMAC均隨著宏觀切向速度的增大而逐漸減小,這種趨勢(shì)隨著法向速度的增大而越來越不明顯.當(dāng)氣體的切向速度小于 1 . 0σ/τ 時(shí),法向速度對(duì)NMAC的影響較小.當(dāng)切向速度大于 1 . 0σ/τ 后,法向速度的增大導(dǎo)致NMAC明顯增大,法向速度對(duì)氣體在表面的法向動(dòng)量適應(yīng)的促進(jìn)作用開始顯現(xiàn)出來.同時(shí)我們注意到,當(dāng)來流切向速度較大且法向速度較小時(shí),NMAC出現(xiàn)了值為負(fù)的情況,這意味著氣體分子反射后的法向動(dòng)量大于入射時(shí)的法向動(dòng)量,說明氣固散射過程中發(fā)生了切向動(dòng)量向法向動(dòng)量的轉(zhuǎn)移,且這種動(dòng)量轉(zhuǎn)移程度隨著切向速度的增大以及法向速度的減小而變強(qiáng).

        圖7為不同來流速度下的EAC,從圖中可以明顯地看出,在相同的切向速度條件下,法向速度的增大有助于氣體與表面的能量適應(yīng),與此相反,在相同的法向速度條件下,切向速度的增大則不利于氣體與表面的能量適應(yīng).氣體分子與表面碰撞時(shí),來流切向速度的增大使得氣體分子與表面作用時(shí)間減少,導(dǎo)致氣-固能量交換不充分;而入射氣體分子的法向能量大小決定了氣體分子入射到表面勢(shì)能區(qū)域的深度,氣體分子進(jìn)入表面勢(shì)能區(qū)域越深,氣體分子與表面的能量交換就越充分,則EAC越大.

        3.2 粗糙表面下宏觀速度的影響分析

        上一節(jié)中對(duì)宏觀速度在光滑表面下對(duì)氣體分子動(dòng)量和能量與表面的適應(yīng)規(guī)律分別展開了分析,但大多數(shù)情況下,固體表面往往具有一定的粗糙度,這就需要我們對(duì)宏觀速度與粗糙度的綜合性影響展開研究.本節(jié)計(jì)算了以不同宏觀速度入射的氣體分子在不同表面上的TMAC,NMAC和EAC.

        具有不同宏觀切向速度及法向速度的氣體分子對(duì)粗糙度大小為0.5 A的粗糙表面的TMAC、吸附概率、NMAC及EAC分別如圖8—圖11所示.

        圖8 粗糙度為 0.5 A 的粗糙表面下,不同宏觀切向速度及法向速度下的切向動(dòng)量適應(yīng)系數(shù)Fig.8.TMAC values of gas molecular on rough surfaces of 0.5 A under different macroscopic tangential and normal velocities.

        圖9 粗糙度為 0.5 A 的粗糙表面下,不同宏觀切向速度及法向速度下的氣體分子吸附概率Fig.9.Sticking probability of gas molecular on rough surfaces of 0.5 A under different macroscopic tangential and normal velocities.

        圖10 粗糙度為 0.5 A 的粗糙表面下,不同宏觀切向速度及法向速度下的法向動(dòng)量適應(yīng)系數(shù)Fig.10.NMAC values of gas molecular on rough surfaces of 0.5 A under different macroscopic tangential and normal velocities.

        圖11 粗糙度為 0.5 A 的粗糙表面下,不同宏觀切向速度及法向速度下的能量適應(yīng)系數(shù)Fig.11.EAC values of gas molecular on rough surfaces of 0.5 A under different macroscopic tangential and normal velocities.

        從圖8可知,在不同的速度下,對(duì)粗糙度高度為0.5 A的表面的TMAC的變化范圍在0.7—0.85之間,TMAC的數(shù)值相比于在光滑表面下大幅度提高,同時(shí)TMAC的變化范圍相對(duì)于光滑表面來說卻是收窄的,這說明粗糙度的出現(xiàn)使氣體分子與表面的切向動(dòng)量適應(yīng)程度提高了,同時(shí)也降低了TMAC對(duì)宏觀速度變化的敏感度.在TMAC和宏觀速度的相關(guān)性上,TMAC隨著切向速度的增大而減小,隨著法向速度的增大而降低,這和氣體分子對(duì)光滑表面的動(dòng)量適應(yīng)規(guī)律有所不同.

        圖9為氣體分子對(duì)粗糙度為0.5 A的表面的吸附概率隨宏觀速度的變化.從圖中可知,吸附概率對(duì)宏觀速度比較敏感,來流的速度越大,則氣體分子的吸附概率就越小.當(dāng)來流切向與法向速度都較低時(shí),氣體分子在表面的吸附概率可以達(dá)到10%以上,相比于在光滑表面最大的吸附概率有很大的提高.當(dāng)來流速度逐漸增大后,氣體分子對(duì)表面的吸附概率降低,當(dāng)法向速度大于1.0 σ/τ,氣體分子對(duì)表面的吸附基本可以忽略.

        粗糙度為0.5 A時(shí),氣體分子在不同來流條件下對(duì)表面的NMAC如圖10所示,在圖中我們?nèi)匀话l(fā)現(xiàn)了NMAC為負(fù)數(shù)的情況,這說明對(duì)于粗糙表面,速度分量之間發(fā)生的動(dòng)量轉(zhuǎn)移現(xiàn)象依舊存在.此外,在同樣的速度條件下,NMAC 在粗糙表面條件下,相比于光滑表面要有所減小,表明粗糙度有助于氣體分子切向與法向動(dòng)量的轉(zhuǎn)移.

        粗糙度為0.5 A時(shí),氣體分子在不同來流條件下對(duì)表面的EAC如圖11所示,由圖可知,對(duì)于相同的來流法向速度,氣體分子對(duì)表面的EAC對(duì)切向速度的變化不敏感.在各法向速度下,EAC的值隨切向速度的增大而產(chǎn)生微小波動(dòng),其變化幅值均小于5%.可見,影響氣體分子與表面能量適應(yīng)的主要因素為法向速度的大小,法向速度越大,氣體分子與粗糙表面的能量適應(yīng)就越差(相比于光滑表面,粗糙表面下的EAC更大,說明粗糙度有助于氣固能量的適應(yīng)).

        圖12 粗糙度為 1.0 A 的粗糙表面下,不同宏觀切向速度及法向速度下的切向動(dòng)量適應(yīng)系數(shù)Fig.12.TMAC values of gas molecular on rough surfaces of 1.0 A under different macroscopic tangential and normal velocities.

        圖13 粗糙度為 1.0 A 的粗糙表面下,不同宏觀切向速度及法向速度下的氣體分子吸附概率Fig.13.Sticking probability of gas molecular on rough surfaces of 1.0 A under different macroscopic tangential and normal velocities.

        不同宏觀切向速度及法向速度下,氣體分子對(duì)粗糙度大小為1.0 A的粗糙表面的TMAC、吸附幾率、NMAC及EAC分別如圖12—圖15所示.由圖可知,各適應(yīng)系數(shù)在粗糙度為1.0 A的情況下,對(duì)不同速度的變化規(guī)律和表面粗糙度為0.5 A時(shí)的變化規(guī)律基本一致.但還是有一些小的變化,主要表現(xiàn)在以下幾個(gè)方面: 1)對(duì)于不同的來流速度,TMAC均在0.95以上,這說明隨著粗糙度的增大,TMAC對(duì)來流速度的敏感性進(jìn)一步降低;2)當(dāng)法向速度較大時(shí),氣體分子在表面的吸附概率提升,且吸附概率表現(xiàn)出隨切向速度的增大而明顯減小的特性;3)隨著粗糙度的增大,氣體分子在切向速度小于 2 . 0σ/τ 時(shí),NMAC對(duì)來流速度變化不敏感,當(dāng)切向速度大于 2 . 0σ/τ 后,NMAC 對(duì)來流速度的變化開始變得明顯起來;4)在不同的速度條件下,氣體分子的EAC隨著粗糙度的增大而顯著提升.

        圖14 粗糙度為 1.0 A 的粗糙表面下,不同宏觀切向速度及法向速度下的法向動(dòng)量適應(yīng)系數(shù)Fig.14.NMAC values of gas molecular on rough surfaces of 1.0 A under different macroscopic tangential and normal velocities.

        圖15 粗糙度為 1.0 A 的粗糙表面下,不同宏觀切向速度及法向速度下的能量適應(yīng)系數(shù)Fig.15.EAC values of gas molecular on rough surfaces of 1.0 A under different macroscopic tangential and normal velocities.

        4 結(jié) 論

        與宏觀連續(xù)流動(dòng)不同,稀薄條件下氣體分子之間的碰撞頻率顯著減小,氣體分子與固體表面的相互作用對(duì)于流動(dòng)的影響起主導(dǎo)作用.本文從微觀角度出發(fā),通過構(gòu)建合理的物理模型,選取恰當(dāng)?shù)臍怏w-表面間相互作用勢(shì)能函數(shù),實(shí)現(xiàn)了對(duì)單個(gè)氣體分子與表面碰撞的準(zhǔn)確模擬,研究方法真實(shí)可靠.采用速度取樣法獲得不同來流條件下氣體分子的入射速度,通過模擬得到反射速度,然后由入射和反射的平均動(dòng)量和能量計(jì)算得到了氣體分子對(duì)光滑和粗糙表面的 TMAC,NMAC,EAC以及氣體分子在表面的吸附概率,探討了切向速度、法向速度以及表面粗糙度對(duì)適應(yīng)系數(shù)的影響規(guī)律,并對(duì)其中的影響機(jī)理進(jìn)行分析,結(jié)論如下:

        1)在宏觀法向速度確定的情況下,對(duì)于相同狀況的表面,氣體分子的 TMAC,NMAC和EAC都隨著切向速度的增大而降低,表明切向速度的增大對(duì)氣體分子與表面的切向動(dòng)量、法向動(dòng)量及總的能量適應(yīng)都有著一定的消極作用;當(dāng)切向速度較大時(shí)會(huì)發(fā)生動(dòng)量由切向向法向轉(zhuǎn)移的情況,而這種轉(zhuǎn)移趨勢(shì)會(huì)隨著切向速度的增大而變得更加明顯;

        2)法向速度對(duì)動(dòng)量和能量適應(yīng)系數(shù)的影響較為復(fù)雜,對(duì)于光滑和粗糙的表面,適應(yīng)系數(shù)對(duì)法向系數(shù)的變化規(guī)律呈現(xiàn)出了相反的規(guī)律;對(duì)于光滑表面,法向速度的增大有助于氣體分子與表面的切向動(dòng)量和能量適應(yīng);而對(duì)于粗糙表面,氣體分子與表面的切向動(dòng)量和能量適應(yīng)程度卻隨著法向速度的增大而減弱;

        3)氣體分子對(duì)表面的適應(yīng)系數(shù)對(duì)來流宏觀速度以及表面粗糙度都表現(xiàn)出了一定的相關(guān)性,傳統(tǒng)散射核模型以及滑移邊界條件采用不變的常數(shù)來表征氣體分子與表面的適應(yīng)規(guī)律是不合理的;同時(shí),氣體分子在固體表面發(fā)生的動(dòng)量和能量的轉(zhuǎn)移現(xiàn)象也是改進(jìn)邊界條件的重點(diǎn).這些從微觀角度發(fā)現(xiàn)的現(xiàn)象對(duì)于理解稀薄氣體流動(dòng)機(jī)理,提高粒子類數(shù)值模擬方法的準(zhǔn)確性具有重要意義.

        猜你喜歡
        切向速度法向動(dòng)量
        動(dòng)量守恒定律在三個(gè)物體系中的應(yīng)用
        旋風(fēng)分離器內(nèi)氣相旋轉(zhuǎn)流不穩(wěn)定性的實(shí)驗(yàn)研究
        落石法向恢復(fù)系數(shù)的多因素聯(lián)合影響研究
        應(yīng)用動(dòng)量守恒定律解題之秘訣
        雙旋流耦合式旋流反應(yīng)器內(nèi)切向速度分布研究
        動(dòng)量相關(guān)知識(shí)的理解和應(yīng)用
        低溫狀態(tài)下的材料法向發(fā)射率測(cè)量
        落石碰撞法向恢復(fù)系數(shù)的模型試驗(yàn)研究
        不透明材料波段法向發(fā)射率在線測(cè)量方法
        重介質(zhì)微型旋流器內(nèi)切向速度的數(shù)值模擬
        日韩有码中文字幕在线观看 | 久久精品视频91| 中文字幕你懂的一区二区| 国产免费在线观看不卡| 日本久久精品视频免费| 日韩精品专区av无码| 猫咪免费人成网站在线观看| 亚洲一区二区久久青草| 精品亚洲一区二区在线观看| 狠狠色丁香婷婷综合潮喷| 精品国产av最大网站| 亚洲av成人一区二区三区网址| 成人爽a毛片免费网站中国| av天堂午夜精品一区| 黑人玩弄人妻中文在线| 中文无码日韩欧免费视频| av一区二区三区观看| 亚洲情综合五月天| 无码a∨高潮抽搐流白浆| 国内精品福利在线视频| 亚洲一区极品美女写真在线看| 国产在线91精品观看| 天天爽夜夜爱| 中文字幕国产欧美| 大量老肥熟女老女人自拍| 国产精品久久免费中文字幕| 久久综合国产乱子伦精品免费| 无码天堂亚洲国产av麻豆| 精品久久精品久久精品| 国产成人综合日韩精品无码| 日本午夜免费福利视频| 中文无码制服丝袜人妻AV| 国产亚洲av成人噜噜噜他| 亚洲人成影院在线观看| 国内精品一区二区2021在线| 亚洲精品国产成人久久av盗摄| 成人内射国产免费观看| 麻豆久久久9性大片| 亚洲欧美日韩国产综合久| 国内精品毛片av在线播放| 97碰碰碰人妻无码视频|