• <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-二氯乙烷的選擇性加氫脫氯
    你看到小船在移動了嗎?
    女人十人毛片免费观看3o分钟| 亚洲精品乱码久久久v下载方式| 小蜜桃在线观看免费完整版高清| 国产一级毛片七仙女欲春2| 国产不卡一卡二| 18禁在线播放成人免费| xxxwww97欧美| 成人一区二区视频在线观看| av天堂中文字幕网| 18禁黄网站禁片免费观看直播| 久久精品综合一区二区三区| 最近在线观看免费完整版| 国产久久久一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 91午夜精品亚洲一区二区三区 | 91麻豆精品激情在线观看国产| 精品免费久久久久久久清纯| 中文字幕免费在线视频6| 国产成+人综合+亚洲专区| 嫩草影视91久久| 亚洲av成人av| 国产亚洲av嫩草精品影院| 亚洲精品乱码久久久v下载方式| 757午夜福利合集在线观看| 自拍偷自拍亚洲精品老妇| 五月玫瑰六月丁香| 99精品久久久久人妻精品| 亚洲人成网站高清观看| 999久久久精品免费观看国产| 久久久久久久久大av| 女人十人毛片免费观看3o分钟| 国产精品伦人一区二区| 国产精品乱码一区二三区的特点| 最新中文字幕久久久久| 岛国在线免费视频观看| 别揉我奶头 嗯啊视频| 国产毛片a区久久久久| 美女大奶头视频| 俄罗斯特黄特色一大片| 最近最新免费中文字幕在线| 禁无遮挡网站| 亚洲欧美精品综合久久99| 国产精品亚洲美女久久久| 自拍偷自拍亚洲精品老妇| 噜噜噜噜噜久久久久久91| 亚洲最大成人手机在线| 99热6这里只有精品| 亚洲av电影不卡..在线观看| 内地一区二区视频在线| 国内精品久久久久精免费| 欧美xxxx黑人xx丫x性爽| 中国美女看黄片| 国产精华一区二区三区| 三级毛片av免费| 成年女人毛片免费观看观看9| 亚洲精品乱码久久久v下载方式| 免费av观看视频| 日韩欧美在线乱码| 欧美3d第一页| 黄色丝袜av网址大全| 国产伦精品一区二区三区四那| 国产免费男女视频| 欧美日韩黄片免| 久久久久久久久久黄片| 熟女人妻精品中文字幕| 国产亚洲精品av在线| 久久久久亚洲av毛片大全| 露出奶头的视频| 国内精品久久久久久久电影| 日本熟妇午夜| 高清日韩中文字幕在线| 精品一区二区三区视频在线| 国产在线精品亚洲第一网站| 一区二区三区激情视频| 色哟哟·www| 99精品久久久久人妻精品| 国产欧美日韩精品一区二区| 国产午夜精品久久久久久一区二区三区 | 亚洲精华国产精华精| 日本三级黄在线观看| 国产精品98久久久久久宅男小说| 亚洲18禁久久av| 成年人黄色毛片网站| 久久久精品大字幕| 国产高潮美女av| 亚洲内射少妇av| 精品人妻熟女av久视频| 欧美激情国产日韩精品一区| 18禁黄网站禁片免费观看直播| 亚洲av熟女| 欧美中文日本在线观看视频| 一级黄色大片毛片| 国产91精品成人一区二区三区| 亚洲成av人片免费观看| 免费av观看视频| 国产午夜精品久久久久久一区二区三区 | 老熟妇乱子伦视频在线观看| 男女床上黄色一级片免费看| 美女大奶头视频| 国产一区二区激情短视频| 在线免费观看不下载黄p国产 | 又粗又爽又猛毛片免费看| 久久人人爽人人爽人人片va | 九九在线视频观看精品| 免费电影在线观看免费观看| 久9热在线精品视频| 窝窝影院91人妻| 此物有八面人人有两片| 69人妻影院| 蜜桃久久精品国产亚洲av| av福利片在线观看| 国产av不卡久久| 51国产日韩欧美| 亚洲国产日韩欧美精品在线观看| 日韩欧美精品v在线| 岛国在线免费视频观看| 高潮久久久久久久久久久不卡| 亚洲中文字幕一区二区三区有码在线看| 长腿黑丝高跟| 亚洲成a人片在线一区二区| 欧美乱色亚洲激情| 宅男免费午夜| 99在线视频只有这里精品首页| av专区在线播放| 99国产极品粉嫩在线观看| 午夜亚洲福利在线播放| 一卡2卡三卡四卡精品乱码亚洲| 国产乱人伦免费视频| 国产私拍福利视频在线观看| 亚洲人成伊人成综合网2020| 波野结衣二区三区在线| 欧美日韩中文字幕国产精品一区二区三区| 久9热在线精品视频| 久久草成人影院| 国产淫片久久久久久久久 | 色5月婷婷丁香| 12—13女人毛片做爰片一| 色综合亚洲欧美另类图片| 国产精品自产拍在线观看55亚洲| 成人av在线播放网站| 直男gayav资源| 熟妇人妻久久中文字幕3abv| 日韩中文字幕欧美一区二区| 午夜免费成人在线视频| 久久久久久久久大av| 成年女人永久免费观看视频| 一区二区三区免费毛片| 99热6这里只有精品| 国产精品永久免费网站| 在线国产一区二区在线| 三级男女做爰猛烈吃奶摸视频| 村上凉子中文字幕在线| 99热只有精品国产| 狠狠狠狠99中文字幕| 免费在线观看亚洲国产| 天堂网av新在线| 日韩精品青青久久久久久| 亚洲av电影在线进入| av女优亚洲男人天堂| 亚洲一区高清亚洲精品| 琪琪午夜伦伦电影理论片6080| 69人妻影院| 免费在线观看亚洲国产| 国产大屁股一区二区在线视频| 亚洲av中文字字幕乱码综合| 老熟妇乱子伦视频在线观看| 国产男靠女视频免费网站| 麻豆国产97在线/欧美| 久久6这里有精品| av天堂在线播放| 国产黄a三级三级三级人| 久久天躁狠狠躁夜夜2o2o| 免费观看的影片在线观看| 中文字幕免费在线视频6| 内射极品少妇av片p| 亚洲成人中文字幕在线播放| 九色成人免费人妻av| 国产精华一区二区三区| 真实男女啪啪啪动态图| 我的老师免费观看完整版| 日本五十路高清| 亚洲 国产 在线| 久久久久久国产a免费观看| 亚洲人成电影免费在线| 国产av麻豆久久久久久久| 亚洲一区二区三区色噜噜| 国产蜜桃级精品一区二区三区| 日本免费a在线| 欧美色视频一区免费| 欧美日韩亚洲国产一区二区在线观看| 日日夜夜操网爽| av福利片在线观看| 一级a爱片免费观看的视频| 床上黄色一级片| 午夜福利在线观看吧| 中出人妻视频一区二区| 在线国产一区二区在线| 永久网站在线| 国产免费av片在线观看野外av| 嫩草影院精品99| 国产麻豆成人av免费视频| 小蜜桃在线观看免费完整版高清| 成人国产一区最新在线观看| 欧美xxxx黑人xx丫x性爽| 国产激情偷乱视频一区二区| 一级a爱片免费观看的视频| 给我免费播放毛片高清在线观看| 成人高潮视频无遮挡免费网站| 亚洲,欧美精品.| 久久久久性生活片| 亚洲av二区三区四区| 女生性感内裤真人,穿戴方法视频| 亚洲av电影不卡..在线观看| 国产男靠女视频免费网站| 99riav亚洲国产免费| 亚洲中文日韩欧美视频| 五月玫瑰六月丁香| 国内久久婷婷六月综合欲色啪| 日韩人妻高清精品专区| 久久亚洲真实| 欧美乱色亚洲激情| 88av欧美| 免费看日本二区| 日本黄色视频三级网站网址| 男女视频在线观看网站免费| 久久久久久久久久成人| 亚洲av日韩精品久久久久久密| 亚洲精品乱码久久久v下载方式| 亚洲精品色激情综合| 一本一本综合久久| 伊人久久精品亚洲午夜| 成年女人永久免费观看视频| 精品久久久久久久久久免费视频| 精品久久久久久,| 精品欧美国产一区二区三| 内射极品少妇av片p| 十八禁网站免费在线| 久99久视频精品免费| 国产伦一二天堂av在线观看| 久久精品国产自在天天线| 欧美成人a在线观看| 国产成人a区在线观看| 亚洲第一区二区三区不卡| 欧美性猛交黑人性爽| 极品教师在线视频| .国产精品久久| 国产精品久久久久久亚洲av鲁大| 欧美日本视频| 99久久久亚洲精品蜜臀av| 欧美性猛交黑人性爽| 老司机福利观看| 亚洲av五月六月丁香网| 黄片小视频在线播放| 成人性生交大片免费视频hd| 波多野结衣巨乳人妻| 亚洲成人中文字幕在线播放| 精品久久国产蜜桃| 久久人人爽人人爽人人片va | 男女那种视频在线观看| 欧美区成人在线视频| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| 小说图片视频综合网站| 久久人人爽人人爽人人片va | 国产精品久久久久久久久免 | 日本 av在线| 18美女黄网站色大片免费观看| 99热这里只有精品一区| 国产精品国产高清国产av| 中文字幕人成人乱码亚洲影| 久久亚洲真实| 成人av在线播放网站| 欧美激情在线99| 成年女人毛片免费观看观看9| 美女被艹到高潮喷水动态| 性欧美人与动物交配| 国产欧美日韩精品一区二区| 亚洲最大成人中文| 少妇的逼水好多| 国内精品美女久久久久久| 欧美xxxx黑人xx丫x性爽| 99久久成人亚洲精品观看| 国产三级中文精品| 一边摸一边抽搐一进一小说| 欧美最黄视频在线播放免费| 在线观看av片永久免费下载| 国产一区二区激情短视频| 欧美中文日本在线观看视频| 欧美日韩综合久久久久久 | 亚洲 欧美 日韩 在线 免费| 欧美性感艳星| 国产久久久一区二区三区| 在线观看av片永久免费下载| 一个人观看的视频www高清免费观看| 国产在视频线在精品| 国内精品美女久久久久久| 69人妻影院| 丁香六月欧美| 他把我摸到了高潮在线观看| 亚洲人成伊人成综合网2020| 在线天堂最新版资源| 国产一区二区在线av高清观看| 男女之事视频高清在线观看| 免费观看精品视频网站| 日本五十路高清| 亚洲五月婷婷丁香| 熟妇人妻久久中文字幕3abv| 亚洲综合色惰| 亚洲成av人片免费观看| 又黄又爽又免费观看的视频| 亚洲熟妇熟女久久| 可以在线观看的亚洲视频| 搡老熟女国产l中国老女人| avwww免费| 国产成人aa在线观看| 在线a可以看的网站| 女生性感内裤真人,穿戴方法视频| 国产视频内射| 精品国产亚洲在线| 超碰av人人做人人爽久久| 久久伊人香网站| 日韩免费av在线播放| 色播亚洲综合网| 免费观看精品视频网站| 搡女人真爽免费视频火全软件 | 91午夜精品亚洲一区二区三区 | 亚洲av电影在线进入| 亚洲在线观看片| 国产三级在线视频| 久久久久久久精品吃奶| 中文字幕免费在线视频6| 国产一级毛片七仙女欲春2| 亚洲av一区综合| 亚洲av成人精品一区久久| 亚洲乱码一区二区免费版| 午夜福利视频1000在线观看| 少妇的逼水好多| 男女做爰动态图高潮gif福利片| 99热6这里只有精品| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| 级片在线观看| 91av网一区二区| 午夜精品久久久久久毛片777| 一本综合久久免费| 99热只有精品国产| 蜜桃亚洲精品一区二区三区| 免费av观看视频| 国产高清激情床上av| 99热6这里只有精品| 麻豆久久精品国产亚洲av| а√天堂www在线а√下载| 精品久久久久久成人av| 老熟妇仑乱视频hdxx| 免费搜索国产男女视频| 真实男女啪啪啪动态图| 好男人电影高清在线观看| 久久久久久久精品吃奶| 日日摸夜夜添夜夜添小说| 国产精品一区二区性色av| 亚洲第一区二区三区不卡| 国产免费av片在线观看野外av| 简卡轻食公司| 人妻制服诱惑在线中文字幕| 美女高潮喷水抽搐中文字幕| 99热6这里只有精品| 在线天堂最新版资源| 18禁在线播放成人免费| 男人的好看免费观看在线视频| 真实男女啪啪啪动态图| 两个人的视频大全免费| 欧美一区二区国产精品久久精品| 一进一出抽搐动态| 中文在线观看免费www的网站| 国内精品美女久久久久久| 九九久久精品国产亚洲av麻豆| 丝袜美腿在线中文| 国产在线精品亚洲第一网站| av天堂在线播放| 色播亚洲综合网| 天堂网av新在线| 日韩成人在线观看一区二区三区| 国模一区二区三区四区视频| 成人欧美大片| 日韩中字成人| 成人无遮挡网站| 9191精品国产免费久久| 香蕉av资源在线| 女人被狂操c到高潮| 身体一侧抽搐| 成人午夜高清在线视频| 国产精品电影一区二区三区| 一本久久中文字幕| 日日摸夜夜添夜夜添小说| 久久久成人免费电影| 亚洲欧美精品综合久久99| 麻豆久久精品国产亚洲av| 女人十人毛片免费观看3o分钟| 国产又黄又爽又无遮挡在线| 亚洲真实伦在线观看| 真人一进一出gif抽搐免费| 亚洲经典国产精华液单 | 国产探花在线观看一区二区| 国产精品av视频在线免费观看| 99热这里只有是精品在线观看 | 美女 人体艺术 gogo| 老熟妇仑乱视频hdxx| 亚洲欧美日韩高清专用| 欧美最新免费一区二区三区 | 91麻豆av在线| 亚洲第一区二区三区不卡| 日韩精品中文字幕看吧| 久久婷婷人人爽人人干人人爱| 午夜福利在线观看吧| 日本熟妇午夜| 狠狠狠狠99中文字幕| 一本综合久久免费| 成年女人永久免费观看视频| 日韩大尺度精品在线看网址| 日韩成人在线观看一区二区三区| 一本久久中文字幕| 久久国产精品影院| 观看美女的网站| 午夜视频国产福利| 国产69精品久久久久777片| 欧美乱色亚洲激情| 欧美黄色片欧美黄色片| 精品午夜福利视频在线观看一区| 动漫黄色视频在线观看| 性色av乱码一区二区三区2| 精品国内亚洲2022精品成人| 3wmmmm亚洲av在线观看| 麻豆成人午夜福利视频| 亚洲经典国产精华液单 | 国内精品美女久久久久久| 亚洲av一区综合| 国产毛片a区久久久久| 最近最新免费中文字幕在线| 国产一级毛片七仙女欲春2| 成人av在线播放网站| 亚洲精品在线观看二区| 成人亚洲精品av一区二区| 久久久成人免费电影| 国产成人a区在线观看| 国产一区二区亚洲精品在线观看| 欧美激情国产日韩精品一区| 欧美xxxx黑人xx丫x性爽| 999久久久精品免费观看国产| 简卡轻食公司| 亚洲欧美日韩无卡精品| 黄色丝袜av网址大全| 波多野结衣高清作品| 真人做人爱边吃奶动态| 性色av乱码一区二区三区2| 亚洲经典国产精华液单 | 欧美又色又爽又黄视频| 老司机福利观看| 91麻豆精品激情在线观看国产| 免费在线观看日本一区| 美女xxoo啪啪120秒动态图 | av女优亚洲男人天堂| 婷婷精品国产亚洲av在线| 亚洲精品456在线播放app | 国模一区二区三区四区视频| 悠悠久久av| 夜夜夜夜夜久久久久| 最新中文字幕久久久久| 国产又黄又爽又无遮挡在线| 丰满乱子伦码专区| 在线十欧美十亚洲十日本专区| 性插视频无遮挡在线免费观看| 国产午夜福利久久久久久| 三级毛片av免费| 18禁黄网站禁片免费观看直播| 尤物成人国产欧美一区二区三区| 国产一区二区在线av高清观看| 国产伦在线观看视频一区| 亚洲精华国产精华精| 国产精品一区二区性色av| 99久久精品国产亚洲精品| 少妇丰满av| 日本一本二区三区精品| 90打野战视频偷拍视频| 看十八女毛片水多多多| 一级作爱视频免费观看| 亚洲成人中文字幕在线播放| 十八禁国产超污无遮挡网站| 我要看日韩黄色一级片| 免费av不卡在线播放| 好男人在线观看高清免费视频| av中文乱码字幕在线| 内地一区二区视频在线| 床上黄色一级片| 久久久久久久久久成人| 亚洲欧美清纯卡通| 欧美又色又爽又黄视频| 别揉我奶头 嗯啊视频| 亚洲内射少妇av| 内射极品少妇av片p| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 亚洲avbb在线观看| 久久精品国产清高在天天线| 精品久久久久久久人妻蜜臀av| 日日干狠狠操夜夜爽| 亚洲第一区二区三区不卡| 日本撒尿小便嘘嘘汇集6| 国内揄拍国产精品人妻在线| 免费在线观看日本一区| 亚洲av第一区精品v没综合| 日韩中文字幕欧美一区二区| 99热这里只有是精品在线观看 | 国产激情偷乱视频一区二区| 欧美激情在线99| 人妻夜夜爽99麻豆av| 一级av片app| 丰满人妻一区二区三区视频av| 午夜免费男女啪啪视频观看 | 亚洲av一区综合| 午夜两性在线视频| 亚洲成a人片在线一区二区| 精品国产三级普通话版| a级一级毛片免费在线观看| 偷拍熟女少妇极品色| 熟女电影av网| 老司机午夜十八禁免费视频| 无遮挡黄片免费观看| 亚洲欧美精品综合久久99| 黄色女人牲交| 人人妻人人看人人澡| 夜夜爽天天搞| 国产三级中文精品| 久久久久精品国产欧美久久久| 亚洲无线在线观看| 一级黄色大片毛片| 在线看三级毛片| 高清日韩中文字幕在线| 日日摸夜夜添夜夜添av毛片 | 十八禁国产超污无遮挡网站| 国产精品女同一区二区软件 | 亚洲av成人av| 国产毛片a区久久久久| 成人性生交大片免费视频hd| 久久亚洲精品不卡| 午夜视频国产福利| 久久午夜福利片| 亚洲精品在线观看二区| 中文在线观看免费www的网站| 黄色视频,在线免费观看| 亚洲成人免费电影在线观看| 美女高潮的动态| 欧美黄色片欧美黄色片| 一个人看视频在线观看www免费| 免费观看的影片在线观看| 亚洲欧美日韩无卡精品| 看免费av毛片| 亚洲黑人精品在线| 亚洲av成人不卡在线观看播放网| 国产av一区在线观看免费| 国产精品自产拍在线观看55亚洲| 亚洲中文字幕日韩| 性欧美人与动物交配| 久久婷婷人人爽人人干人人爱| 波多野结衣巨乳人妻| 亚洲精品久久国产高清桃花| 自拍偷自拍亚洲精品老妇| 亚洲人成网站高清观看| av视频在线观看入口| 国产免费一级a男人的天堂| 人妻久久中文字幕网| 少妇的逼好多水| 亚洲经典国产精华液单 | 一区福利在线观看| 小蜜桃在线观看免费完整版高清| 久久久久国内视频| 久久国产乱子伦精品免费另类| 夜夜躁狠狠躁天天躁| 久久久久久久精品吃奶| 日韩欧美国产一区二区入口| 真实男女啪啪啪动态图| 国产高清有码在线观看视频| 好男人电影高清在线观看| 色尼玛亚洲综合影院| 听说在线观看完整版免费高清| 欧美区成人在线视频| 久久精品国产99精品国产亚洲性色| 国产欧美日韩精品亚洲av| 精品福利观看| 97超级碰碰碰精品色视频在线观看| 亚洲,欧美精品.| 天天一区二区日本电影三级| 可以在线观看的亚洲视频| 亚洲av.av天堂| 两性午夜刺激爽爽歪歪视频在线观看| 简卡轻食公司| 久久精品国产清高在天天线| 免费看光身美女| 国产麻豆成人av免费视频| 宅男免费午夜| 又黄又爽又免费观看的视频| 99热6这里只有精品| 88av欧美| 99久久九九国产精品国产免费| 久久久久性生活片| 欧美精品国产亚洲| 亚洲av五月六月丁香网| 亚洲成a人片在线一区二区| 长腿黑丝高跟| 又紧又爽又黄一区二区| 蜜桃亚洲精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 全区人妻精品视频| 欧美极品一区二区三区四区| 中文字幕熟女人妻在线| 精品一区二区三区人妻视频|