• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    水流在旋轉(zhuǎn)黑磷納米管內(nèi)軸向驅(qū)動特性*

    2020-06-30 12:12:08張忠強范晉偉張福建程廣貴2丁建寧2
    物理學報 2020年11期
    關鍵詞:黑磷納米管水分子

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

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

    運用分子動力學方法模擬研究了旋轉(zhuǎn)的黑磷納米管對管內(nèi)水流的軸向驅(qū)動特性, 研究結(jié)果表明: 手性黑磷納米管在旋轉(zhuǎn)時會驅(qū)動管內(nèi)水分子沿軸向運動, 運動方向由納米管轉(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)時, 管間水分子的軸向運動會增強. 納米管半徑也會對水分子的定向運動產(chǎn)生影響, 具體表現(xiàn)為在相同轉(zhuǎn)速下, 隨著納米管半徑的增大, 管內(nèi)水分子在軸向上的運動速度會減小, 而受力則會增大; 雙壁黑磷納米管在旋轉(zhuǎn)時管內(nèi)水分子的軸向運動情況和單壁黑磷納米管模型差異很小, 證明黑磷管層數(shù)對水流驅(qū)動效果的影響不明顯; 溫度對水流驅(qū)動效果的影響規(guī)律取決于管內(nèi)壓強和溫度對流固界面摩擦系數(shù)的耦合作用, 當溫度低于常溫時水分子在軸向上的速度和受力會隨著溫度的升高而增大, 當溫度達到常溫時則趨于平穩(wěn). 研究結(jié)果可為基于黑磷納米管的流體傳動器件的設計和應用提供理論基礎.

    1 引 言

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

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

    2 模型和模擬方法

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

    本文采用分子動力學(MD)方法研究了黑磷納米管旋轉(zhuǎn)作用下管內(nèi)水分子軸向運動的特性, 利用LAMMPS開源代碼完成模擬. 首先在正則(NVT)系綜下對系統(tǒng)進行弛豫使系統(tǒng)達到穩(wěn)定狀態(tài), 弛豫溫度為300 K. 模擬過程中采用Langevin恒溫器使水的溫度保持在300 K的恒溫. 整個模擬過程的時間為4 ns, 其中前1 ns用來進行弛豫, 當系統(tǒng)達到穩(wěn)定狀態(tài)后令黑磷納米管以一定的轉(zhuǎn)速繞軸線轉(zhuǎn)動, 持續(xù)時間為3 ns, 統(tǒng)計后2 ns的數(shù)據(jù)用于分析研究. 在勢函數(shù)選擇上, 采用Stillinger-Weber (SW)勢函數(shù)[21]描述黑磷原子之間的相互作用, 水分子模型選取改進版的四位TIP4P-Ew模型, 水分子之間的長程庫侖勢作用通過particle-particle particle-mesh (PPPM)方法計算. 黑磷與水分子相互作用采用Lennard-Jones (LJ)勢函數(shù)進行描述, 作用參數(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)水分子的軸向運動特性

    為了研究手性黑磷納米管旋轉(zhuǎn)作用對管內(nèi)水分子的軸向驅(qū)動特性, 構(gòu)建了具有不同手性角度的黑磷納米管, 黑磷納米管的半徑均設置為R =2.6 nm, 并且將黑磷管的轉(zhuǎn)動速度設置為繞順時針方向50 rad/ns, 研究黑磷納米管自轉(zhuǎn)作用下管內(nèi)水分子的運動情況. 由于構(gòu)建納米管模型需要單層黑磷邊界處的原子相互匹配, 同時為了保證黑磷納米管管徑的一致性, 本文目前只構(gòu)造了0°, 23.4°,90° 三種手性黑磷納米管模型. 圖2所示為不同手性角度的黑磷納米管旋轉(zhuǎn)時管內(nèi)水分子的軸向速度以及水分子軸向上受到的驅(qū)動力隨時間的變化關系. 由圖易知, 手性角度為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)了定向運動的特性, 因此本文對此進行重點討論.

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

    圖 2 不同手性角度的黑磷納米管以50 rad/ns的轉(zhuǎn)速順時針旋轉(zhuǎn)時管內(nèi)水分子沿軸線方向的(a)速度和(b)受力隨時間的變化關系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易知, 當黑磷納米管以相同轉(zhuǎn)動速度沿不同方向轉(zhuǎn)動時, 水分子軸向速度和驅(qū)動力均會保持穩(wěn)定, 且呈現(xiàn)出大小相等, 方向相反的現(xiàn)象. 由此可推斷出管內(nèi)水分子的運動方向由黑磷納米管的轉(zhuǎn)動方向決定, 并且可進一步證明手性角度為23.4°的黑磷納米管在繞自身軸線旋轉(zhuǎn)時會驅(qū)動管內(nèi)水分子沿軸線方向定向運動.

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

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

    圖 3 手性角度為23.4°的黑磷納米管以50 rad/ns的轉(zhuǎn)速沿不同方向旋轉(zhuǎn)時管內(nèi)水分子沿軸向的(a)速度和(b)受力隨時間的變化關系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)動速度的變化關系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ū)動作用的機理分析

    黑磷納米管旋轉(zhuǎn)時水分子的軸向運動是由黑磷與水分子在界面處的相互作用引起的, 為了進一步研究手性角度為23.4°的黑磷納米管旋轉(zhuǎn)作用下水分子定向運動的原因, 利用相同手性角度的單層黑磷構(gòu)建納米通道并在其中填充水分子, 構(gòu)建出Couette流系統(tǒng)模型如圖5所示, 其中單層黑磷的尺寸為3.3 nm × 4.3 nm, 通道寬度為5.2 nm(上文中黑磷納米管的直徑). 令上下兩層黑磷分別以大小相等方向相反的速度對水分子進行剪切, 剪切速度的大小分別為25, 50, 75, 100 m/s. 沿通道寬度的方向以0.3 nm的高度對水分子進行分層, 統(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)水分子的運動速度在寬度方向上近似為線性分布, 并且隨著剪切速度的增加, 邊界處水分子的速度增大. 通過對曲線進行擬合得到其斜率和截距, 斜率即為剪切應變率, 并進一步計算出不同剪切應變率下水分子的邊界速度, 令剪切速度和邊界速度相減得到滑移速度, 同時通過模擬結(jié)果得到的單層黑磷對水分子的切向力, 計算出不同剪切應變率下的剪切應力, 統(tǒng)計滑移速度和剪切應力隨剪切應變率變化的趨勢如圖7所示.

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

    圖 7 (a)水分子的邊界滑移速度和(b)剪切應力隨剪切應變率的變化關系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)時管間水分子的軸向運動特性

    由上文可知, 水分子能夠沿軸向運動是由于手性黑磷納米管內(nèi)表面對水分子的作用引起的, 而黑磷納米管外表面具有和內(nèi)表面一樣的褶皺結(jié)構(gòu), 推測手性黑磷納米管在旋轉(zhuǎ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種情形下水分子的軸向運動特性. 3種情形下水分子的軸向運動速度和受力如圖10所示.

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

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

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

    圖 10 3種情形下黑磷納米管間水分子沿軸線方向的(a)速度和(b)受力隨轉(zhuǎ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)水分子的運動特性而言也是一個重要參數(shù). 為了研究手性角度為23.4°的黑磷納米管在旋轉(zhuǎn)時半徑對管內(nèi)水分子運動規(guī)律的影響, 構(gòu)建了半徑分別為1.59, 2.66,3.72 nm的手性黑磷納米管模型, 在管內(nèi)填充水分子并使黑磷管以不同轉(zhuǎn)速轉(zhuǎn)動, 統(tǒng)計出水分子在軸向上的運動速度和受力情況如圖11所示.

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

    參照多壁碳納米管的概念構(gòu)建手性角度為23.4°的多壁黑磷納米管, 討論黑磷納米管層數(shù)對水分子運動的影響. 由于自然狀態(tài)下黑磷的層間距為0.53 nm, 而本文中選擇的LJ勢函數(shù)的截斷半徑為1.2 nm, 當黑磷納米管的層數(shù)大于2時, 外層黑磷納米管的作用很小, 因此本節(jié)僅討論雙壁黑磷納米管(DWBPNT)的情形. 參考自然環(huán)境下黑磷層與層之間的間隔, 將兩層黑磷納米管之間的距離設置為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)速的變化關系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)速的變化趨勢基本一致, 并且當轉(zhuǎn)速相同時兩種情形下黑磷納米管內(nèi)水分子的軸向速度差異很小, 而受力則會因為層數(shù)的增加而有輕微的上升, 由此可以看出黑磷納米管在旋轉(zhuǎn)時, 其層數(shù)對管內(nèi)水分子的軸向運動特性影響不明顯.

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

    圖 12 手性角度為23.4°時, 不同層數(shù)黑磷納米管內(nèi)水分子沿軸線方向的(a)速度和(b)受力隨黑磷納米管轉(zhuǎ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)水分子的軸向速度與受力隨溫度的變化關系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可以看出, 溫度較低時, 水分子沿軸向的速度和受力會隨著溫度的升高而明顯增大, 而隨著溫度逐漸上升到常溫, 溫度對水分子的軸向速度和受力的影響逐漸減弱, 水分子的運動趨于穩(wěn)定.

    4 結(jié) 論

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

    猜你喜歡
    黑磷納米管水分子
    電化學剝離黑磷制備納米黑磷研究進展
    多少水分子才能稱“一滴水”
    科教新報(2021年11期)2021-05-12 19:50:11
    最近鄰弱交換相互作用對spin-1納米管磁化強度的影響
    塊狀和少層黑磷的合成
    黑磷的礦化法制備及其形成機制研究進展
    磷肥與復肥(2020年2期)2020-05-06 06:55:14
    缺陷和硫摻雜黑磷的第一性原理計算
    為什么濕的紙會粘在一起?
    科學之謎(2016年9期)2016-10-11 08:59:04
    二氧化鈦納米管的制備及其應用進展
    應用化工(2014年11期)2014-08-16 15:59:13
    TiO2納米管負載Pd-Ag催化1,2-二氯乙烷的選擇性加氫脫氯
    你看到小船在移動了嗎?
    亚洲成人久久爱视频| 精品久久久久久成人av| 中国国产av一级| 国产精品一区二区三区四区免费观看| 成人综合一区亚洲| 国产黄色视频一区二区在线观看 | 综合色丁香网| 又粗又爽又猛毛片免费看| 欧美日本视频| 午夜精品一区二区三区免费看| 欧美丝袜亚洲另类| 国产片特级美女逼逼视频| 日韩欧美精品免费久久| 久久综合国产亚洲精品| 亚洲精品乱久久久久久| 免费看a级黄色片| 亚洲精品自拍成人| 91久久精品电影网| 欧美区成人在线视频| 天堂影院成人在线观看| 欧美3d第一页| 中文乱码字字幕精品一区二区三区 | 99在线视频只有这里精品首页| 超碰av人人做人人爽久久| 少妇的逼好多水| 高清日韩中文字幕在线| 国产三级在线视频| 免费观看在线日韩| 久久久久免费精品人妻一区二区| 啦啦啦观看免费观看视频高清| 亚洲久久久久久中文字幕| 日韩亚洲欧美综合| 国语对白做爰xxxⅹ性视频网站| 丰满乱子伦码专区| 免费电影在线观看免费观看| 免费播放大片免费观看视频在线观看 | 国产精品.久久久| 国产午夜精品一二区理论片| 波多野结衣巨乳人妻| 亚洲三级黄色毛片| 久久精品熟女亚洲av麻豆精品 | 六月丁香七月| 欧美zozozo另类| 精品免费久久久久久久清纯| 精品免费久久久久久久清纯| 久久综合国产亚洲精品| 亚洲av一区综合| 舔av片在线| 亚洲av电影在线观看一区二区三区 | 免费观看人在逋| 一级毛片电影观看 | 国产淫片久久久久久久久| 2021少妇久久久久久久久久久| 最后的刺客免费高清国语| 69人妻影院| 中国国产av一级| 狠狠狠狠99中文字幕| 深爱激情五月婷婷| 又粗又硬又长又爽又黄的视频| 久久99蜜桃精品久久| 亚洲精品亚洲一区二区| 纵有疾风起免费观看全集完整版 | 久久久精品大字幕| 国产成人精品一,二区| 又黄又爽又刺激的免费视频.| 别揉我奶头 嗯啊视频| 亚洲久久久久久中文字幕| 久久精品国产99精品国产亚洲性色| 久久99精品国语久久久| 黄色一级大片看看| 国产伦在线观看视频一区| 啦啦啦韩国在线观看视频| 日韩视频在线欧美| 国产在视频线在精品| 尤物成人国产欧美一区二区三区| 色尼玛亚洲综合影院| 三级国产精品片| 熟女人妻精品中文字幕| 成人无遮挡网站| 日韩av在线大香蕉| 99在线人妻在线中文字幕| 亚洲av成人精品一二三区| 久久久午夜欧美精品| 亚洲欧美精品自产自拍| 成人亚洲欧美一区二区av| 久久精品国产亚洲av涩爱| 国产高清三级在线| 亚洲精品456在线播放app| 国产成人aa在线观看| 春色校园在线视频观看| 午夜亚洲福利在线播放| 麻豆国产97在线/欧美| 亚洲18禁久久av| 亚洲av成人精品一二三区| 成年女人看的毛片在线观看| 国产成人freesex在线| 精品久久久久久久人妻蜜臀av| videos熟女内射| 久久精品国产亚洲av天美| 国产91av在线免费观看| 亚洲一区高清亚洲精品| 国产美女午夜福利| 国产麻豆成人av免费视频| 中文字幕精品亚洲无线码一区| 丰满少妇做爰视频| 国产精品久久久久久av不卡| 我要搜黄色片| 久久久久久国产a免费观看| 亚洲一级一片aⅴ在线观看| 亚洲精品乱久久久久久| 99在线人妻在线中文字幕| 亚洲av福利一区| 亚洲最大成人av| 最近中文字幕高清免费大全6| 狠狠狠狠99中文字幕| 一区二区三区乱码不卡18| 日韩亚洲欧美综合| 简卡轻食公司| 亚洲欧美一区二区三区国产| 三级国产精品片| 欧美另类亚洲清纯唯美| 久久久久久久久久久丰满| av在线老鸭窝| 国产色婷婷99| 成人午夜精彩视频在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲经典国产精华液单| 三级男女做爰猛烈吃奶摸视频| 美女脱内裤让男人舔精品视频| 国产极品精品免费视频能看的| 亚洲av日韩在线播放| 久久欧美精品欧美久久欧美| 成人三级黄色视频| 日日撸夜夜添| 欧美精品国产亚洲| 免费不卡的大黄色大毛片视频在线观看 | 性色avwww在线观看| 久久婷婷人人爽人人干人人爱| 精品国产三级普通话版| 日韩三级伦理在线观看| 久久精品影院6| 亚洲国产欧美在线一区| 亚洲最大成人av| 少妇裸体淫交视频免费看高清| 毛片一级片免费看久久久久| 亚洲欧美清纯卡通| 麻豆成人午夜福利视频| 色视频www国产| 色尼玛亚洲综合影院| 久久这里有精品视频免费| 亚洲在久久综合| 亚洲中文字幕日韩| 国产爱豆传媒在线观看| 日韩av在线免费看完整版不卡| 欧美又色又爽又黄视频| 草草在线视频免费看| 亚洲内射少妇av| 我要看日韩黄色一级片| 偷拍熟女少妇极品色| 久久韩国三级中文字幕| 精品免费久久久久久久清纯| 国产成人a∨麻豆精品| 两个人视频免费观看高清| 一区二区三区免费毛片| 国语对白做爰xxxⅹ性视频网站| 一二三四中文在线观看免费高清| 嫩草影院入口| 黄色一级大片看看| 黄色一级大片看看| 最近中文字幕2019免费版| 日本-黄色视频高清免费观看| 成人毛片a级毛片在线播放| 五月伊人婷婷丁香| 高清av免费在线| 日韩成人av中文字幕在线观看| 一区二区三区四区激情视频| 日本欧美国产在线视频| 国产精品蜜桃在线观看| 啦啦啦韩国在线观看视频| 国产精品乱码一区二三区的特点| 国产白丝娇喘喷水9色精品| 三级经典国产精品| 国产免费福利视频在线观看| www日本黄色视频网| 亚洲熟妇中文字幕五十中出| 日日干狠狠操夜夜爽| 日韩亚洲欧美综合| 99久国产av精品| 亚洲在线自拍视频| 91狼人影院| 国产av码专区亚洲av| 男人舔奶头视频| 我的女老师完整版在线观看| 亚洲av福利一区| 一区二区三区四区激情视频| 久久国内精品自在自线图片| 一级av片app| 精品久久久久久久人妻蜜臀av| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品乱码久久久久久按摩| 国产真实伦视频高清在线观看| 高清日韩中文字幕在线| 九九热线精品视视频播放| 午夜免费激情av| 国产69精品久久久久777片| 丝袜喷水一区| 国产成人一区二区在线| 一级黄色大片毛片| 国产精品综合久久久久久久免费| 亚洲精品日韩av片在线观看| 久久久久久久久久黄片| 乱人视频在线观看| 午夜老司机福利剧场| 久久久久九九精品影院| 麻豆乱淫一区二区| 真实男女啪啪啪动态图| 少妇熟女aⅴ在线视频| 亚洲精品aⅴ在线观看| 亚洲真实伦在线观看| 国产精品人妻久久久影院| 青青草视频在线视频观看| 一边摸一边抽搐一进一小说| 国产视频内射| 精品久久久久久久久久久久久| 特级一级黄色大片| 全区人妻精品视频| 亚洲激情五月婷婷啪啪| 久久这里只有精品中国| or卡值多少钱| 国产精华一区二区三区| 一级av片app| 夜夜爽夜夜爽视频| 看免费成人av毛片| 国产三级在线视频| 国产精品一区二区三区四区久久| av在线播放精品| 国产美女午夜福利| 国产精品永久免费网站| 久久精品国产亚洲av涩爱| 91午夜精品亚洲一区二区三区| 成人三级黄色视频| 嫩草影院精品99| 禁无遮挡网站| 亚洲av不卡在线观看| 国产午夜精品一二区理论片| 插逼视频在线观看| 亚洲电影在线观看av| 天天躁夜夜躁狠狠久久av| 亚洲精品日韩在线中文字幕| 亚洲激情五月婷婷啪啪| 亚洲自偷自拍三级| 日韩成人av中文字幕在线观看| 一本一本综合久久| 直男gayav资源| 国产探花极品一区二区| 亚洲真实伦在线观看| 国产精品一区二区三区四区久久| 久久精品久久精品一区二区三区| 欧美成人a在线观看| 日韩 亚洲 欧美在线| 国产综合懂色| 看十八女毛片水多多多| av在线观看视频网站免费| 久久人人爽人人爽人人片va| 欧美一区二区亚洲| 男插女下体视频免费在线播放| 国产精品国产三级国产专区5o | 久久久亚洲精品成人影院| 黄片无遮挡物在线观看| 亚洲av.av天堂| 深爱激情五月婷婷| 日本欧美国产在线视频| 国产乱来视频区| 亚洲图色成人| 国产淫语在线视频| 精品国产一区二区三区久久久樱花 | 麻豆国产97在线/欧美| 亚洲国产欧洲综合997久久,| 午夜福利在线观看免费完整高清在| 九草在线视频观看| 中文天堂在线官网| 在线免费观看的www视频| 免费看美女性在线毛片视频| 国产亚洲精品久久久com| 91精品国产九色| 一级毛片aaaaaa免费看小| 欧美不卡视频在线免费观看| 91午夜精品亚洲一区二区三区| 午夜福利在线在线| 精品不卡国产一区二区三区| 一本一本综合久久| 看十八女毛片水多多多| 亚洲aⅴ乱码一区二区在线播放| 国产高清三级在线| 国产综合懂色| 亚洲人成网站在线观看播放| 中文亚洲av片在线观看爽| 最新中文字幕久久久久| 久久人人爽人人爽人人片va| 午夜福利在线观看吧| 三级国产精品片| 可以在线观看毛片的网站| 精品熟女少妇av免费看| 久久99热这里只有精品18| 久久99热这里只频精品6学生 | 免费观看精品视频网站| 91午夜精品亚洲一区二区三区| 欧美日韩精品成人综合77777| 国产精品日韩av在线免费观看| 日韩欧美精品免费久久| 成人二区视频| 91狼人影院| 纵有疾风起免费观看全集完整版 | 三级男女做爰猛烈吃奶摸视频| 国产精品综合久久久久久久免费| 人人妻人人澡欧美一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 女人十人毛片免费观看3o分钟| 国产老妇伦熟女老妇高清| 在线播放国产精品三级| 伊人久久精品亚洲午夜| 中文字幕av成人在线电影| 久久久久久久久大av| 91在线精品国自产拍蜜月| 一级黄片播放器| 亚洲人成网站高清观看| 女人久久www免费人成看片 | 日日摸夜夜添夜夜爱| 国产黄色视频一区二区在线观看 | 久久热精品热| www.色视频.com| 亚洲国产最新在线播放| 日韩 亚洲 欧美在线| 十八禁国产超污无遮挡网站| 看黄色毛片网站| 免费观看在线日韩| 日韩精品有码人妻一区| 搡老妇女老女人老熟妇| 韩国av在线不卡| 国产激情偷乱视频一区二区| 插逼视频在线观看| 黄片wwwwww| 一级二级三级毛片免费看| 99久久精品国产国产毛片| 精品久久久久久成人av| 高清视频免费观看一区二区 | 日韩欧美在线乱码| 国产精品久久久久久久电影| 欧美日韩精品成人综合77777| 青青草视频在线视频观看| 麻豆精品久久久久久蜜桃| 国产精品野战在线观看| 日韩av不卡免费在线播放| 亚洲国产欧美在线一区| 天天躁日日操中文字幕| 欧美高清性xxxxhd video| 熟妇人妻久久中文字幕3abv| 在线观看av片永久免费下载| 亚洲自偷自拍三级| 最后的刺客免费高清国语| 在线播放国产精品三级| 少妇高潮的动态图| 热99在线观看视频| 69av精品久久久久久| a级毛色黄片| 97超碰精品成人国产| 日本-黄色视频高清免费观看| 美女黄网站色视频| 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 国产成人精品久久久久久| 婷婷六月久久综合丁香| 激情 狠狠 欧美| 亚洲成人久久爱视频| 天堂影院成人在线观看| 国产av码专区亚洲av| 五月伊人婷婷丁香| 日韩av在线大香蕉| 国产久久久一区二区三区| .国产精品久久| 亚洲综合色惰| 男插女下体视频免费在线播放| 久久精品国产鲁丝片午夜精品| 麻豆成人av视频| 2021天堂中文幕一二区在线观| 亚洲av电影在线观看一区二区三区 | 少妇熟女欧美另类| 亚洲丝袜综合中文字幕| 国产亚洲5aaaaa淫片| 国产精品人妻久久久影院| 99在线视频只有这里精品首页| 国产又黄又爽又无遮挡在线| 人体艺术视频欧美日本| 亚洲欧美日韩无卡精品| 久久精品久久精品一区二区三区| 精品人妻熟女av久视频| 亚洲在线自拍视频| 一个人看的www免费观看视频| 菩萨蛮人人尽说江南好唐韦庄 | 三级毛片av免费| 亚洲国产日韩欧美精品在线观看| 日韩欧美精品v在线| 国产高清有码在线观看视频| 18禁动态无遮挡网站| 内射极品少妇av片p| 亚洲av.av天堂| 国产亚洲精品久久久com| 联通29元200g的流量卡| 日本免费在线观看一区| 亚洲av福利一区| 在线免费十八禁| 亚洲精华国产精华液的使用体验| 亚洲精品久久久久久婷婷小说 | 性插视频无遮挡在线免费观看| 国产69精品久久久久777片| 国产毛片a区久久久久| 精品久久久久久电影网 | 日韩成人伦理影院| 日韩av不卡免费在线播放| 亚洲国产欧美人成| 超碰97精品在线观看| 搞女人的毛片| 欧美bdsm另类| 亚洲av一区综合| 直男gayav资源| 美女脱内裤让男人舔精品视频| 亚洲欧美成人综合另类久久久 | 国产在线男女| ponron亚洲| 中国国产av一级| 欧美丝袜亚洲另类| 欧美日韩国产亚洲二区| av免费观看日本| 亚洲五月天丁香| 美女被艹到高潮喷水动态| 婷婷色麻豆天堂久久 | 久久久久久久亚洲中文字幕| 看十八女毛片水多多多| 国产精品av视频在线免费观看| 亚洲va在线va天堂va国产| 国产 一区 欧美 日韩| 午夜福利网站1000一区二区三区| 成年女人永久免费观看视频| 国产精品久久久久久久久免| 白带黄色成豆腐渣| 亚洲综合色惰| 久久精品国产亚洲网站| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 99久久精品国产国产毛片| 欧美97在线视频| 2021天堂中文幕一二区在线观| 一级爰片在线观看| 69av精品久久久久久| 欧美97在线视频| 日本色播在线视频| 亚洲最大成人手机在线| 中国美白少妇内射xxxbb| 日韩制服骚丝袜av| 非洲黑人性xxxx精品又粗又长| 岛国毛片在线播放| 床上黄色一级片| 精品99又大又爽又粗少妇毛片| 国产免费一级a男人的天堂| 视频中文字幕在线观看| 精品人妻一区二区三区麻豆| 欧美一区二区精品小视频在线| 69人妻影院| 国产高清三级在线| 只有这里有精品99| 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 99久久中文字幕三级久久日本| 久久99热6这里只有精品| 日本免费在线观看一区| 少妇裸体淫交视频免费看高清| 国产成人精品一,二区| 2021天堂中文幕一二区在线观| 99久久精品国产国产毛片| 国产午夜精品论理片| 草草在线视频免费看| 人妻系列 视频| 全区人妻精品视频| 一个人免费在线观看电影| 久久这里有精品视频免费| 搡女人真爽免费视频火全软件| 51国产日韩欧美| 欧美日韩综合久久久久久| 一级二级三级毛片免费看| 亚洲丝袜综合中文字幕| 搡女人真爽免费视频火全软件| 免费黄色在线免费观看| 熟女电影av网| 一夜夜www| 亚洲无线观看免费| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美精品综合久久99| 久久热精品热| 禁无遮挡网站| 国产黄片美女视频| 青青草视频在线视频观看| 成人综合一区亚洲| 久久精品久久久久久噜噜老黄 | 亚洲欧美精品自产自拍| 卡戴珊不雅视频在线播放| 免费av观看视频| 亚洲熟妇中文字幕五十中出| 欧美极品一区二区三区四区| 午夜福利在线在线| 亚洲伊人久久精品综合 | 日韩欧美在线乱码| 亚洲精品色激情综合| 欧美极品一区二区三区四区| 欧美又色又爽又黄视频| 大香蕉97超碰在线| h日本视频在线播放| 91久久精品国产一区二区三区| 又黄又爽又刺激的免费视频.| 色噜噜av男人的天堂激情| 非洲黑人性xxxx精品又粗又长| 天堂中文最新版在线下载 | 99久久人妻综合| av播播在线观看一区| 国产精品美女特级片免费视频播放器| av福利片在线观看| 最近中文字幕高清免费大全6| 一级毛片aaaaaa免费看小| 中文精品一卡2卡3卡4更新| 国产精品蜜桃在线观看| 自拍偷自拍亚洲精品老妇| 久久久久网色| 亚洲成人av在线免费| 亚洲欧美日韩卡通动漫| 亚洲国产日韩欧美精品在线观看| 国产成人精品婷婷| 真实男女啪啪啪动态图| 国产私拍福利视频在线观看| 干丝袜人妻中文字幕| 国产乱人视频| av女优亚洲男人天堂| 寂寞人妻少妇视频99o| 国产在视频线精品| 九九在线视频观看精品| 亚洲国产精品成人综合色| 日韩,欧美,国产一区二区三区 | 久久久久久久久中文| 99在线人妻在线中文字幕| 在线播放国产精品三级| 亚洲欧美日韩无卡精品| 国产精品美女特级片免费视频播放器| 成人毛片a级毛片在线播放| 国产一区二区在线av高清观看| 亚洲av不卡在线观看| 国产极品天堂在线| 亚洲欧美一区二区三区国产| 天堂中文最新版在线下载 | 久久人人爽人人片av| 亚洲成色77777| 少妇猛男粗大的猛烈进出视频 | 亚洲综合色惰| 成人二区视频| 国产爱豆传媒在线观看| 能在线免费看毛片的网站| 色综合色国产| 18禁在线播放成人免费| 亚洲在线观看片| 少妇被粗大猛烈的视频| 日韩一区二区三区影片| 成年女人看的毛片在线观看| 日日干狠狠操夜夜爽| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 国产在视频线在精品| eeuss影院久久| 韩国av在线不卡| 91久久精品电影网| 午夜福利视频1000在线观看| 久久鲁丝午夜福利片| 国产精品三级大全| 国产日韩欧美在线精品| 日韩欧美三级三区| 国产91av在线免费观看| 日韩欧美精品v在线| 秋霞在线观看毛片| 91午夜精品亚洲一区二区三区| 国内精品宾馆在线| 国产精品1区2区在线观看.| 日日摸夜夜添夜夜爱| 亚洲欧美成人精品一区二区| 国国产精品蜜臀av免费| 精品久久久噜噜| 亚洲av.av天堂| 久久久久性生活片| 久久热精品热| 午夜爱爱视频在线播放| 麻豆久久精品国产亚洲av| 99久久精品国产国产毛片| 三级男女做爰猛烈吃奶摸视频| 一二三四中文在线观看免费高清| 尾随美女入室| 99九九线精品视频在线观看视频| 国产成人精品一,二区| 国产成年人精品一区二区| 国产成人a区在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产又黄又爽又无遮挡在线| 搡女人真爽免费视频火全软件| 亚洲美女搞黄在线观看| 欧美日韩精品成人综合77777| a级毛色黄片| 蜜臀久久99精品久久宅男| 亚洲成人久久爱视频| 国产三级在线视频| 久久精品久久久久久久性| 国产精品一区二区性色av|