張忠強(qiáng) 范晉偉 張福建 程廣貴2) 丁建寧2)?

1) (江蘇大學(xué)， 智能柔性機(jī)械電子研究院， 鎮(zhèn)江 212013)2) (常州大學(xué)， 江蘇省光伏科學(xué)與工程協(xié)同創(chuàng)新中心， 常州 213164)(2020 年1 月16日收到; 2020 年3 月25日收到修改稿)

運(yùn)用分子動力學(xué)方法模擬研究了旋轉(zhuǎn)的黑磷納米管對管內(nèi)水流的軸向驅(qū)動特性， 研究結(jié)果表明: 手性黑磷納米管在旋轉(zhuǎn)時會驅(qū)動管內(nèi)水分子沿軸向運(yùn)動， 運(yùn)動方向由納米管轉(zhuǎn)向決定; 管內(nèi)水流的流速和驅(qū)動力會隨著黑磷管轉(zhuǎn)速的提高而增大. 采用黑磷雙壁Couette模型計算分析了水-黑磷界面的摩擦系數(shù)及滑移特性，闡明了黑磷表面天然的各向異性微結(jié)構(gòu)是旋轉(zhuǎn)黑磷管軸向驅(qū)動水流的本質(zhì)原因. 構(gòu)建了在雙層黑磷納米管間填充水分子的模型， 發(fā)現(xiàn)內(nèi)外黑磷管同時旋轉(zhuǎn)時， 管間水分子的軸向運(yùn)動會增強(qiáng). 納米管半徑也會對水分子的定向運(yùn)動產(chǎn)生影響， 具體表現(xiàn)為在相同轉(zhuǎn)速下， 隨著納米管半徑的增大， 管內(nèi)水分子在軸向上的運(yùn)動速度會減小， 而受力則會增大; 雙壁黑磷納米管在旋轉(zhuǎn)時管內(nèi)水分子的軸向運(yùn)動情況和單壁黑磷納米管模型差異很小， 證明黑磷管層數(shù)對水流驅(qū)動效果的影響不明顯; 溫度對水流驅(qū)動效果的影響規(guī)律取決于管內(nèi)壓強(qiáng)和溫度對流固界面摩擦系數(shù)的耦合作用， 當(dāng)溫度低于常溫時水分子在軸向上的速度和受力會隨著溫度的升高而增大， 當(dāng)溫度達(dá)到常溫時則趨于平穩(wěn). 研究結(jié)果可為基于黑磷納米管的流體傳動器件的設(shè)計和應(yīng)用提供理論基礎(chǔ).

1 引 言

近年來， 隨著微/納尺度技術(shù)的逐步成熟， 低維材料由于其在物理、化學(xué)方面的優(yōu)異性能[1，2]， 被廣泛應(yīng)用于微/納機(jī)電系統(tǒng)、柔性電子技術(shù)和傳感器等領(lǐng)域[3-5]. 其中， 微/納流體傳動器件作為新興的研究方向， 具有廣闊的應(yīng)用前景. 由于納米尺度下低維材料與流體在流-固界面上的獨(dú)特特性[6，7]，納米流體器件具有高通量、高流速和高精度等特點(diǎn). 目前學(xué)者們通過施加溫度梯度、壓力梯度、電場以及在通道內(nèi)構(gòu)建螺旋型結(jié)構(gòu)等方式對納米通道內(nèi)的流體進(jìn)行驅(qū)動[8-11]， 研究了流體在通道內(nèi)的運(yùn)動特性， 并以此設(shè)計出各類納米流體傳動器件.

黑磷(BP)是近幾年發(fā)現(xiàn)的二維材料之一， 是磷的同素異形體中熱力學(xué)性能較為穩(wěn)定的結(jié)構(gòu)， 相比于石墨烯和二硫化鉬等材料， 黑磷具有天然帶隙及高載流子遷移率， 是天然的半導(dǎo)體材料， 并且黑磷表面由于褶皺結(jié)構(gòu)的存在而呈現(xiàn)出各向異性[12-15]，針對上述這些特點(diǎn)， 學(xué)者們進(jìn)行了諸多研究. Yang等[16]利用分子動力學(xué)方法對單層黑磷在雙軸和單軸拉伸以及剪切變形下的溫度-應(yīng)力-應(yīng)變關(guān)系進(jìn)行了模擬， 并將結(jié)果與現(xiàn)有的實驗測量值和第一性原理計算值進(jìn)行對比. Zhao等[17]利用光學(xué)顯微鏡觀察了水蒸氣在黑磷表面的液化過程， 發(fā)現(xiàn)在黑磷表面結(jié)構(gòu)各向異性的影響下水滴呈現(xiàn)橢圓形結(jié)構(gòu).Hyun等[18]在可控的溫度和濕度下用原子力顯微鏡對黑磷表面原子的降解過程進(jìn)行了表征. 我們之前利用不同手性的層狀黑磷構(gòu)建Couette流和Poiseuille流模型， 研究了各向異性對黑磷納米通道內(nèi)水分子運(yùn)動特性的影響， 得出界面摩擦系數(shù)隨黑磷褶皺方向與流體運(yùn)動方向間夾角的增加而增加， 而水分子的黏度系數(shù)基本保持不變的結(jié)論[19，20].除此以外， 受到石墨烯與碳納米管結(jié)構(gòu)的啟發(fā)， 研究人員認(rèn)為可以通過卷繞單層黑磷的方式得到黑磷納米管(BPNT)的結(jié)構(gòu). Cai等[21]利用分子模擬方法研究了軸向壓縮條件下單壁黑磷納米管的強(qiáng)度和穩(wěn)定性， 明確了影響黑磷納米管臨界應(yīng)力和應(yīng)變的因素. Hao等[22]在分子動力學(xué)定性分析的基礎(chǔ)上， 研究了尺寸、應(yīng)變和空位對扶手椅型黑磷納米管傳熱性能的影響. 目前關(guān)于黑磷納米管內(nèi)流體運(yùn)動特性的研究鮮有報道， 因此本文構(gòu)建手性黑磷納米管模型， 利用黑磷表面的傾斜褶皺在管壁形成螺旋結(jié)構(gòu)， 從而實現(xiàn)黑磷納米管自轉(zhuǎn)時對水分子在軸向上的驅(qū)動作用， 并研究了各條件參數(shù)對水分子運(yùn)動特性的影響， 為黑磷納米管在納米尺度流體傳動領(lǐng)域的應(yīng)用提供理論基礎(chǔ).

2 模型和模擬方法

本文選用的材料為黑磷納米管， 黑磷是由磷原子通過sp3雜化形成褶皺狀的蜂窩結(jié)構(gòu)， 其特殊的褶皺結(jié)構(gòu)使黑磷表面呈現(xiàn)出各向異性. 通過對單層黑磷進(jìn)行幾何坐標(biāo)變換構(gòu)建黑磷納米管的模型， 并且參考碳納米管手性[23]的概念， 將納米管展開成的單層黑磷的褶皺方向與納米管軸向之間的夾角定義為手性角θ， 構(gòu)建具有不同手性的黑磷納米管，圖1(a)和圖1(b)分別是手性角度為23.4°的單層黑磷和黑磷納米管的模型圖. 為了研究旋轉(zhuǎn)的黑磷納米管中水分子的運(yùn)動特性， 在黑磷納米管中填充水分子， 模型如圖1(c)所示， 令黑磷納米管以一定的轉(zhuǎn)速繞軸線旋轉(zhuǎn)， 并將圖中納米管的旋轉(zhuǎn)方向定義為順時針方向. 由于單層黑磷的原子不在同一平面內(nèi)， 黑磷納米管內(nèi)外表面的原子會具有不同的半徑， 因此使用藍(lán)色和綠色兩種顏色分別代表黑磷納米管內(nèi)外表面的原子. 所構(gòu)建的黑磷納米管長度為4.3 nm， 半徑取黑磷納米管內(nèi)外管徑的平均值，用R表示， 并將黑磷納米管的軸向設(shè)置為周期性，而水分子則選用TIP4P-Ew模型[24].

本文采用分子動力學(xué)(MD)方法研究了黑磷納米管旋轉(zhuǎn)作用下管內(nèi)水分子軸向運(yùn)動的特性， 利用LAMMPS開源代碼完成模擬. 首先在正則(NVT)系綜下對系統(tǒng)進(jìn)行弛豫使系統(tǒng)達(dá)到穩(wěn)定狀態(tài)， 弛豫溫度為300 K. 模擬過程中采用Langevin恒溫器使水的溫度保持在300 K的恒溫. 整個模擬過程的時間為4 ns， 其中前1 ns用來進(jìn)行弛豫， 當(dāng)系統(tǒng)達(dá)到穩(wěn)定狀態(tài)后令黑磷納米管以一定的轉(zhuǎn)速繞軸線轉(zhuǎn)動， 持續(xù)時間為3 ns， 統(tǒng)計后2 ns的數(shù)據(jù)用于分析研究. 在勢函數(shù)選擇上， 采用Stillinger-Weber (SW)勢函數(shù)[21]描述黑磷原子之間的相互作用， 水分子模型選取改進(jìn)版的四位TIP4P-Ew模型， 水分子之間的長程庫侖勢作用通過particle-particle particle-mesh (PPPM)方法計算. 黑磷與水分子相互作用采用Lennard-Jones (LJ)勢函數(shù)進(jìn)行描述， 作用參數(shù)通過Lorentz-Berthelot混合算法計算得出， 具體參數(shù)見表1[25，26]. 其中，LJ作用和庫侖相互作用的截止距離為12 ?和10 ?，并且為了保持水分子的O—H鍵長和H—O—H鍵角的穩(wěn)定性， 本文使用SHAKE算法[27].

圖 1 (a)單層黑磷模型， 其中手性角度θ指黑磷褶皺方向與z軸方向(納米管軸向)的夾角; (b)手性角度為23.4°的黑磷納米管;(c)填充水分子的黑磷納米管旋轉(zhuǎn)模型圖Fig. 1. (a) Monolayer black phosphorus model， chiral angle θ is the intersection angle between the ripple direction of BP monolayer and z direction (the axial direction of the BPNT); (b) BPNT with a chiral angle of 23.4°; (c) model of the rotating BPNT filled with water molecules.

表 1 LJ勢能函數(shù)的參數(shù)值Table 1. Parameter values of LJ potential function.

3 模擬結(jié)果討論與分析

3.1 手性黑磷納米管旋轉(zhuǎn)時管內(nèi)水分子的軸向運(yùn)動特性

為了研究手性黑磷納米管旋轉(zhuǎn)作用對管內(nèi)水分子的軸向驅(qū)動特性， 構(gòu)建了具有不同手性角度的黑磷納米管， 黑磷納米管的半徑均設(shè)置為R =2.6 nm， 并且將黑磷管的轉(zhuǎn)動速度設(shè)置為繞順時針方向50 rad/ns， 研究黑磷納米管自轉(zhuǎn)作用下管內(nèi)水分子的運(yùn)動情況. 由于構(gòu)建納米管模型需要單層黑磷邊界處的原子相互匹配， 同時為了保證黑磷納米管管徑的一致性， 本文目前只構(gòu)造了0°， 23.4°，90° 三種手性黑磷納米管模型. 圖2所示為不同手性角度的黑磷納米管旋轉(zhuǎn)時管內(nèi)水分子的軸向速度以及水分子軸向上受到的驅(qū)動力隨時間的變化關(guān)系. 由圖易知， 手性角度為0°和90°時， 黑磷納米管內(nèi)水分子的軸向速度和受力均在0上下波動， 而手性角度為23.4°的黑磷納米管在旋轉(zhuǎn)時， 水分子的軸向速度和受力則始終為正， 且速度大小穩(wěn)定在1.5 m/s附近， 受力大小約為15 kcal·mol—1·?—1.由此可知， 相比于手性角度為0°和90°的黑磷納米管， 手性角度為23.4°的黑磷納米管在旋轉(zhuǎn)時管內(nèi)水分子出現(xiàn)了定向運(yùn)動的特性， 因此本文對此進(jìn)行重點(diǎn)討論.

令手性角度為23.4°的黑磷納米管繞軸線分別沿順時針(clockwise)和逆時針(anticlockwise)方向轉(zhuǎn)動， 轉(zhuǎn)動速度大小均為50 rad/ns， 統(tǒng)計水分子的軸向速度與受力隨時間的變化趨勢如圖3所示.為了減小計算誤差的影響， 每種情形通過改變隨機(jī)數(shù)的方式分別計算3次， 取3次數(shù)據(jù)的平均值并畫出了誤差帶.

圖 2 不同手性角度的黑磷納米管以50 rad/ns的轉(zhuǎn)速順時針旋轉(zhuǎn)時管內(nèi)水分子沿軸線方向的(a)速度和(b)受力隨時間的變化關(guān)系Fig. 2. For the angular velocity of the BPNT being 50 rad/ns， (a) the velocity in the axial direction of water molecules in BPNTs and (b) the resultant force in the axial direction of water molecules received from BPNTs with different chiral angles as a function of time.

由圖3易知， 當(dāng)黑磷納米管以相同轉(zhuǎn)動速度沿不同方向轉(zhuǎn)動時， 水分子軸向速度和驅(qū)動力均會保持穩(wěn)定， 且呈現(xiàn)出大小相等， 方向相反的現(xiàn)象. 由此可推斷出管內(nèi)水分子的運(yùn)動方向由黑磷納米管的轉(zhuǎn)動方向決定， 并且可進(jìn)一步證明手性角度為23.4°的黑磷納米管在繞自身軸線旋轉(zhuǎn)時會驅(qū)動管內(nèi)水分子沿軸線方向定向運(yùn)動.

為了研究黑磷納米管轉(zhuǎn)速對水分子運(yùn)動的影響， 令手性角度為23.4°的黑磷納米管繞軸線沿順時針方向轉(zhuǎn)動， 轉(zhuǎn)動速度大小分別為10， 20， 30，40， 50 rad/ns， 統(tǒng)計水分子運(yùn)動穩(wěn)定時的軸向速度與受力， 統(tǒng)計結(jié)果如圖4所示.

由圖4可以看出， 當(dāng)黑磷納米管以不同轉(zhuǎn)速進(jìn)行旋轉(zhuǎn)時， 水分子的軸向運(yùn)動速度和受力會隨著轉(zhuǎn)速的上升而增大， 并且增大的趨勢會逐漸變緩. 由此可知黑磷管對水分子的軸向驅(qū)動效果會隨著黑磷納米管轉(zhuǎn)速的上升而增強(qiáng)， 但增強(qiáng)的趨勢隨著黑磷納米管轉(zhuǎn)速的上升而逐漸減弱.

圖 3 手性角度為23.4°的黑磷納米管以50 rad/ns的轉(zhuǎn)速沿不同方向旋轉(zhuǎn)時管內(nèi)水分子沿軸向的(a)速度和(b)受力隨時間的變化關(guān)系Fig. 3. For the angular velocity of the BPNT being 50 rad/ns in different directions of rotation， (a) the velocity in the axial direction of water molecules in the BPNT and (b) the resultant force in the axial direction of water molecules received from the BPNT as a function of time when the chiral angle is 23.4°.

圖 4 手性角度為23.4°時， 黑磷納米管內(nèi)水分子的軸向速度與受力隨納米管轉(zhuǎn)動速度的變化關(guān)系Fig. 4. The velocity in the axial direction of water molecules in the BPNT and the resultant force in the axial direction of water molecules received from the BPNT as a function of the angular velocity of the BPNT when the chiral angle is 23.4°.

3.2 手性黑磷納米管旋轉(zhuǎn)時對管內(nèi)水流軸向驅(qū)動作用的機(jī)理分析

黑磷納米管旋轉(zhuǎn)時水分子的軸向運(yùn)動是由黑磷與水分子在界面處的相互作用引起的， 為了進(jìn)一步研究手性角度為23.4°的黑磷納米管旋轉(zhuǎn)作用下水分子定向運(yùn)動的原因， 利用相同手性角度的單層黑磷構(gòu)建納米通道并在其中填充水分子， 構(gòu)建出Couette流系統(tǒng)模型如圖5所示， 其中單層黑磷的尺寸為3.3 nm × 4.3 nm， 通道寬度為5.2 nm(上文中黑磷納米管的直徑). 令上下兩層黑磷分別以大小相等方向相反的速度對水分子進(jìn)行剪切， 剪切速度的大小分別為25， 50， 75， 100 m/s. 沿通道寬度的方向以0.3 nm的高度對水分子進(jìn)行分層， 統(tǒng)計水分子的速度分布如圖6所示.

圖 5 黑磷納米通道內(nèi)水分子的Couette流模型圖Fig. 5. Couette flow model diagram of water molecules flowing in BP nanochannel.

圖 6 黑磷納米通道寬度方向上水分子的速度分布Fig. 6. Velocity distribution of water molecules along the width of the BP nanochannel.

由圖6可知， 在黑磷的剪切作用下， 通道內(nèi)水分子的運(yùn)動速度在寬度方向上近似為線性分布， 并且隨著剪切速度的增加， 邊界處水分子的速度增大. 通過對曲線進(jìn)行擬合得到其斜率和截距， 斜率即為剪切應(yīng)變率， 并進(jìn)一步計算出不同剪切應(yīng)變率下水分子的邊界速度， 令剪切速度和邊界速度相減得到滑移速度， 同時通過模擬結(jié)果得到的單層黑磷對水分子的切向力， 計算出不同剪切應(yīng)變率下的剪切應(yīng)力， 統(tǒng)計滑移速度和剪切應(yīng)力隨剪切應(yīng)變率變化的趨勢如圖7所示.

由圖7可知， 水分子的邊界滑移速度和所受到的剪切應(yīng)力均會隨著剪切應(yīng)變率的增大而增大. 而界面摩擦系數(shù)f滿足公式， 根據(jù)圖7的結(jié)果計算出不同剪切應(yīng)變率下的摩擦系數(shù)并取平均， 得到剪切方向上的界面摩擦系數(shù)fx=0.691 MPa·s·m—1. 圖8(a)所示為水分子邊界處的微觀結(jié)構(gòu)圖， 可知黑磷的褶皺結(jié)構(gòu)會引起水分子近似沿褶皺方向分布， 進(jìn)而導(dǎo)致水分子的動量損失.已有研究證明納米管的傳動特性與其手性的關(guān)系取決于納米管勢能面的分布特征[28]， 因此通過計算黑磷與水分子之間的LJ勢能， 繪制出圖8(b)所示的黑磷-水交互界面的勢能分布情況， 圖8(b)顯示黑磷表面由于褶皺結(jié)構(gòu)形成粗糙勢能表面， 且表面勢能分布情況與褶皺的方向(即黑磷的手性)相符， 由此可知黑磷的手性會對黑磷與水分子在界面處的相互作用產(chǎn)生影響. 圖8(c)所示為黑磷對水分子作用力的示意圖， 當(dāng)黑磷的手性角度為23.4°時， 黑磷對水分子的作用力的方向為垂直于褶皺方向， 該作用力在與剪切速度平行和垂直的方向上會產(chǎn)生分力， 根據(jù)力的分解原理可得與剪切方向垂直方向上的摩擦系數(shù)fz= fx· tan23.4° =0.299 MPa·s·m—1. 由此可對黑磷納米管旋轉(zhuǎn)時水分子的軸向運(yùn)動進(jìn)行解釋: 當(dāng)黑磷納米管旋轉(zhuǎn)時，管壁可看作對管內(nèi)水分子進(jìn)行剪切， 由于黑磷和水分子的界面摩擦系數(shù)的存在， 水分子會因黑磷納米管的作用力而運(yùn)動， 當(dāng)手性角度為23.4°時， 黑磷納米管管壁會因其表面勢能的分布特性在軸向上具有一定的摩擦力分量， 從而對水分子產(chǎn)生軸向驅(qū)動力， 并且傾斜的褶皺使得黑磷納米管表面形成螺旋型結(jié)構(gòu)， 因此實現(xiàn)了管內(nèi)水分子的軸向運(yùn)動.

圖 7 (a)水分子的邊界滑移速度和(b)剪切應(yīng)力隨剪切應(yīng)變率的變化關(guān)系Fig. 7. (a) The boundary slip velocity of water molecules and (b) the shear stress as a function of the shear strain rate.

圖 8 (a)水分子邊界處的微觀構(gòu)型圖; (b)黑磷-水交互界面勢能分布圖; (c)黑磷對水分子的作用力示意圖Fig. 8. (a) Microstructure of the boundary of water molecules; (b) potential energy distribution cloud diagram of BP-water;(c) schematic diagram of the force of BP on water molecules.

3.3 雙層手性黑磷納米管旋轉(zhuǎn)時管間水分子的軸向運(yùn)動特性

由上文可知， 水分子能夠沿軸向運(yùn)動是由于手性黑磷納米管內(nèi)表面對水分子的作用引起的， 而黑磷納米管外表面具有和內(nèi)表面一樣的褶皺結(jié)構(gòu)， 推測手性黑磷納米管在旋轉(zhuǎn)時外表面也會促進(jìn)水分子的軸向運(yùn)動. 為了驗證這一想法， 構(gòu)建了手性角度為23.4°的雙層黑磷納米管， 管半徑分別為1.59 nm和3.19 nm， 并在兩管之間填充水分子，模型結(jié)構(gòu)如圖9所示. 分別令外管轉(zhuǎn)動、內(nèi)管轉(zhuǎn)動和內(nèi)外管同時同向轉(zhuǎn)動， 研究3種情形下水分子的軸向運(yùn)動特性. 3種情形下水分子的軸向運(yùn)動速度和受力如圖10所示.

圖 9 雙層黑磷納米管間填充水分子模型Fig. 9. Model of water molecules filling between two BPNTs.

由圖10可知， 3種情形下管間的水分子均產(chǎn)生定向運(yùn)動的現(xiàn)象， 并且當(dāng)內(nèi)管和外管同時旋轉(zhuǎn)時， 水分子的軸向速度和受力相比于另外兩種情形會有所增大， 由此可以看出手性黑磷納米管在旋轉(zhuǎn)時其外表面的褶皺結(jié)構(gòu)也會使水分子產(chǎn)生軸向運(yùn)動， 且內(nèi)外管同時轉(zhuǎn)動會使水分子的軸向運(yùn)動效果增強(qiáng).

3.4 條件參數(shù)對手性黑磷納米管旋轉(zhuǎn)時通道內(nèi)水分子軸向運(yùn)動特性的影響

圖 10 3種情形下黑磷納米管間水分子沿軸線方向的(a)速度和(b)受力隨轉(zhuǎn)速的變化關(guān)系Fig. 10. (a) The velocity in the axial direction of water molecules between BPNTs and (b) the resultant force in the axial direction of water molecules received from BPNTs as a function of the angular velocity of BPNTs in three cases.

管徑對于研究黑磷納米管內(nèi)水分子的運(yùn)動特性而言也是一個重要參數(shù). 為了研究手性角度為23.4°的黑磷納米管在旋轉(zhuǎn)時半徑對管內(nèi)水分子運(yùn)動規(guī)律的影響， 構(gòu)建了半徑分別為1.59， 2.66，3.72 nm的手性黑磷納米管模型， 在管內(nèi)填充水分子并使黑磷管以不同轉(zhuǎn)速轉(zhuǎn)動， 統(tǒng)計出水分子在軸向上的運(yùn)動速度和受力情況如圖11所示.

由圖11可以看出， 不同半徑的黑磷納米管在旋轉(zhuǎn)時， 管內(nèi)水分子在軸向上的速度和受力隨轉(zhuǎn)速的變化趨勢基本一致， 并且當(dāng)轉(zhuǎn)速相同時， 水分子軸向上的速度隨著黑磷納米管半徑的增大而減小，受力卻隨之增大. 這是由于半徑的增大一方面增加了水分子與黑磷納米管管壁的接觸面積， 造成水分子的受力隨之上升; 而另一方面也使得管內(nèi)的水分子數(shù)急劇上升， 從而增加了驅(qū)動的負(fù)擔(dān)， 因此水分子的軸向速度會有所下降.

參照多壁碳納米管的概念構(gòu)建手性角度為23.4°的多壁黑磷納米管， 討論黑磷納米管層數(shù)對水分子運(yùn)動的影響. 由于自然狀態(tài)下黑磷的層間距為0.53 nm， 而本文中選擇的LJ勢函數(shù)的截斷半徑為1.2 nm， 當(dāng)黑磷納米管的層數(shù)大于2時， 外層黑磷納米管的作用很小， 因此本節(jié)僅討論雙壁黑磷納米管(DWBPNT)的情形. 參考自然環(huán)境下黑磷層與層之間的間隔， 將兩層黑磷納米管之間的距離設(shè)置為0.53 nm構(gòu)建雙壁黑磷納米管模型. 在雙壁黑磷納米管內(nèi)填充水分子， 并使黑磷管沿順時針方向以不同轉(zhuǎn)速旋轉(zhuǎn)， 統(tǒng)計管內(nèi)水分子的軸向速度和受力并將其與單壁黑磷納米管(SWBPNT)對比，結(jié)果如圖12所示.

圖 11 手性角度為23.4°時， 不同半徑的黑磷納米管內(nèi)水分子沿軸線方向的(a)速度和(b)受力隨黑磷納米管轉(zhuǎn)速的變化關(guān)系Fig. 11. For different radius， (a) the velocity in the axial direction of water molecules in BPNTs and (b) the resultant force in the axial direction of water molecules received from BPNTs as a function of the angular velocity of BPNTs when the chiral angle is 23.4°.

由圖12可以看出， 單壁和雙壁黑磷納米管內(nèi)水分子在軸向上的速度和受力隨轉(zhuǎn)速的變化趨勢基本一致， 并且當(dāng)轉(zhuǎn)速相同時兩種情形下黑磷納米管內(nèi)水分子的軸向速度差異很小， 而受力則會因為層數(shù)的增加而有輕微的上升， 由此可以看出黑磷納米管在旋轉(zhuǎn)時， 其層數(shù)對管內(nèi)水分子的軸向運(yùn)動特性影響不明顯.

環(huán)境的溫度也可能會對黑磷納米管內(nèi)水分子的運(yùn)動特性產(chǎn)生影響， 因此將模擬溫度分別設(shè)置為200， 225， 250， 275， 300， 325， 350 K， 轉(zhuǎn)速設(shè)置為50 rad/ns， 探究溫度對水分子軸向運(yùn)動的影響.水分子運(yùn)動穩(wěn)定時的軸向速度與受力隨溫度的變化關(guān)系如圖13所示.

圖 12 手性角度為23.4°時， 不同層數(shù)黑磷納米管內(nèi)水分子沿軸線方向的(a)速度和(b)受力隨黑磷納米管轉(zhuǎn)速的變化關(guān)系Fig. 12. For different layers， (a) the velocity in the axial direction of water molecules in BPNTs and (b) the resultant force in the axial direction of water molecules received from BPNTs as a function of the angular velocity of BPNTs when the chiral angle is 23.4°.

圖 13 轉(zhuǎn)速為50 rad/ns時， 手性角度為23.4°的黑磷納米管內(nèi)水分子的軸向速度與受力隨溫度的變化關(guān)系Fig. 13. For the angular velocity of the BPNT being 50 rad/ns，the velocity in the axial direction of water molecules in the BPNT and the resultant force in the axial direction of water molecules received from the BPNT as a function of the temperature when the chiral angle is 23.4°.

由圖13可以看出， 溫度較低時， 水分子沿軸向的速度和受力會隨著溫度的升高而明顯增大， 而隨著溫度逐漸上升到常溫， 溫度對水分子的軸向速度和受力的影響逐漸減弱， 水分子的運(yùn)動趨于穩(wěn)定.

4 結(jié) 論

以黑磷納米管構(gòu)建納米通道， 研究了黑磷納米管自轉(zhuǎn)作用下通道內(nèi)水分子的定向運(yùn)動特性， 并從流-固界面角度對其原因進(jìn)行分析. 由研究結(jié)果可知， 手性黑磷納米管在旋轉(zhuǎn)時會驅(qū)動管內(nèi)水分子沿軸向運(yùn)動， 運(yùn)動方向由納米管轉(zhuǎn)向決定， 并且管內(nèi)水流的流速和驅(qū)動力會隨著納米管轉(zhuǎn)速的提高而增大. 基于單層黑磷構(gòu)建Couette流模型對水-黑磷界面的摩擦系數(shù)及滑移特性進(jìn)行分析， 證明黑磷表面各向異性的微結(jié)構(gòu)是旋轉(zhuǎn)黑磷管對水流進(jìn)行軸向驅(qū)動的本質(zhì)原因. 構(gòu)建了在雙層黑磷納米管間填充水分子的模型， 發(fā)現(xiàn)黑磷管外表面的微結(jié)構(gòu)也會驅(qū)動水流軸向運(yùn)動， 且內(nèi)外黑磷管同時旋轉(zhuǎn)時，管間水流的軸向運(yùn)動效果會增強(qiáng). 最后研究了黑磷管半徑、層數(shù)和環(huán)境溫度的影響， 結(jié)果表明: 相同轉(zhuǎn)速下， 隨著納米管半徑的增大， 管內(nèi)水流在軸向上的運(yùn)動速度會減小， 而受力則會增大; 雙壁黑磷納米管在旋轉(zhuǎn)時管內(nèi)水流的軸向運(yùn)動情況和單壁黑磷納米管模型差異很小， 證明納米管層數(shù)的影響不明顯; 溫度低于常溫時水流的速度和受力會隨著溫度的升高而增大， 當(dāng)溫度達(dá)到常溫后則趨于平穩(wěn).