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

    基于CFD的豎直圓管內(nèi)空氣-水逆向流動界面行為研究

    2023-12-16 05:24:02朱隆祥馬在勇張盧騰潘良明
    核科學與工程 2023年5期
    關鍵詞:波峰波幅氣相

    萬 潔,孫 皖,李 想,朱隆祥,馬在勇,張盧騰,潘良明

    基于CFD的豎直圓管內(nèi)空氣-水逆向流動界面行為研究

    萬潔1,2,孫皖1,*,李想1,2,朱隆祥1,2,馬在勇1,2,張盧騰1,2,潘良明1,2

    (1. 重慶大學低品位能源利用技術及系統(tǒng)教育部重點實驗室,重慶 400044; 2. 重慶大學核工程與核技術系,重慶 400044)

    為研究空氣-水逆向流動局部特性,本文基于VOF(volume of fluid)模型對=25 mm的圓管中氣液逆向流動過程中界面波的演變開展了數(shù)值計算,分析了流動過程中界面波波形的變化,計算結果與實驗中液膜厚度變化符合較好。對不同液相流速工況下的流場變化分析可得,流動過程中界面波前后產(chǎn)生的壓差引入了削弱界面波穩(wěn)定性的附加作用力,且界面波前后壓力變化隨液相流速的增大而增大;與此同時,氣相流速越大,界面波波幅沿軸向發(fā)展越迅速,表明界面波的不穩(wěn)定性隨氣相流速的增大而增大。由此可得,氣相流速及液相流速的增大都會使界面波的穩(wěn)定性降低,即更容易發(fā)生液泛。

    CFD;液泛;VOF;界面波;液膜厚度

    反應堆發(fā)生大破口事故時,應急堆芯冷卻系統(tǒng)(ECCS)投入運行,應急冷卻水與堆芯內(nèi)產(chǎn)生的蒸汽在堆芯孔板、下降段或熱管段等位置處發(fā)生液泛導致堆芯冷卻效率降低甚至無法冷卻從而影響反應堆安全。對于液泛的表面波動特性、氣-液相界面等的局部特性研究,即使是空氣-水逆向流動工況下都有很大不確定性,而數(shù)值模擬因其能較好反映局部特性而被廣泛使用。Jayanti等[1]研究發(fā)現(xiàn)施加在波上的剪切力由波附近的壓力變化引起且與流道尺寸有關,并基于管徑提出了一種區(qū)分液泛機理的方法。Shearer和Davidson[2]則提出了一種可以預測豎直表面上波形及波幅的理論。Da Riva和Del Col[3]采用VOF[4]界面重構法對一個軸對稱域進行了模擬,發(fā)現(xiàn)環(huán)狀波的形成頻率與Barbosa Jr.等[5]的實驗結果符合較好。Murase等[6],Minami等[7,8]和Utanohara等[9]在1/15 的PWR模型中對熱段進行了空氣-水三維逆向流動模擬,對界面摩擦、流型及逆向流動極限(CCFL)進行了研究并總結得出,采用具有合適的界面摩擦關系式的兩流體模型比VOF模型結果更好。兩流體模型雖在對逆向流動極限(CCFL)的模擬中表現(xiàn)較好,但本文關注界面波的界面行為,VOF模型對于相界面的識別更為準確,因此本文采用VOF模型展開模擬研究。

    本文基于VOF模型采用ANANSYS FLUENT研究了逆向流動過程中的界面波動特性,并將模擬結果與Wan等[10]的實驗結果進行了對比。本研究通過跟蹤界面波行為以得到液膜的演變,同時討論了氣相流速對波形變化的影響并對液泛過程中以及液泛后的流場進行了分析,得出氣液逆向流動過程中相間界面的相互作用,以更好地分析CCFL發(fā)生的機理。

    1 數(shù)值模擬

    1.1 幾何模型及邊界條件

    本文中數(shù)值計算幾何區(qū)域為一直徑25 mm,管長156 mm的豎直圓管??紤]到流動過程呈對稱分布所以本研究取管中心為對稱軸,以管左側(cè)區(qū)域的流動為研究對象如圖1所示,注水口長度為7.9 mm,采用常規(guī)光管注水,網(wǎng)格數(shù)量為109 950。管的下部為氣相速度入口,左側(cè)注水口為液相質(zhì)量入口,上部出口為兩相壓力出口,其余管道壁面為無滑移壁面條件。

    圖1 液泛數(shù)值計算幾何模型

    本文采用VOF模型對界面進行追蹤,壓力梯度耦合采用PISO算法,梯度項離散采用最小二乘法單元格方法。除空泡份額以外的變量都采用二階迎風格式進行離散,空泡份額采用幾何重構法離散,瞬態(tài)計算的時間離散采用一階隱式法。此外,當所有變量殘差小于10-4時認為計算達到收斂。本文共選取了包含液泛及液泛后過程的5個工況展開計算,其中液相流速根據(jù)Wan等[10]實驗中兩種液相流速下對應兩種不同的液泛機理分別選取了0.028 m/s及0.121 m/s,氣相流速分別設置為實驗中液泛起始點及液泛后對應的氣相流速工況,工況表如表1所示。所有物性參數(shù)均定義為常溫(25 ℃)常壓(0.1 MPa)下的數(shù)值。

    表1 液泛數(shù)值計算工況

    1.2 網(wǎng)格無關性驗證

    為驗證計算結果的網(wǎng)格無關性,本文采用不同的網(wǎng)格數(shù)量在l=0.028 m/s,g=6.379 m/s工況下進行數(shù)值計算并對比其界面波生長特性。從圖2中可看出,界面波波幅隨時間的變化隨著網(wǎng)格數(shù)量的增加變化不大,本文最終選取了網(wǎng)格數(shù)量109 950進行計算。

    圖2 不同網(wǎng)格數(shù)量下界面波變化特性(jl=0.028 m/s,jg=6.379 m/s)

    圖2 不同網(wǎng)格數(shù)量下界面波變化特性(jl=0.028 m/s,jg=6.379 m/s)(續(xù))

    2 結果與分析

    2.1 波形分析

    2.1.1小液相流速

    因界面波生長具有周期性,本文對某一周期內(nèi)變化進行了分析。如圖3所示,由于氣流的拖曳作用使界面產(chǎn)生的界面波波幅增大,在該工況下波沿徑向的增大也更為顯著。=0.031~0.056 s時,波峰沿著氣流方向不斷被拉長到一臨界點后其克服自身表面張力后被撕碎,分離出液滴隨氣流向上流動,此后界面波在氣流作用下不斷被撕碎使波幅逐漸變小,直至下一個小波產(chǎn)生從而進入到下一個周期的界面波生長過程。圖4中液膜隨時間的變化曲線與模擬結果變化一致,液膜厚度隨時間的推移先增大后逐漸減小,其達到最大值時對應模擬過程中波幅達到峰值。由=0.056 s及0.059 s可看出,波峰處液體被撕裂形成小液滴被氣相夾帶向上導致液膜變小,這是由于界面波在管內(nèi)徑向上不斷增大使得氣相在截面上所占流通面積減小而流速增大,氣液間的相間剪切增大,對應圖4中液膜減小階段。

    圖3 jl=0.028 m/s,jg=6.379 m/s時界面波形變化

    圖4 jl=0.028 m/s,jg=6.379 m/s時液膜在0.60 s內(nèi)的變化[10]

    2.1.2大液相流速

    如圖5所示,大液相流速下界面波的發(fā)展十分迅速。=0.019~0.026 s時,界面波在徑向上快速拉長與對稱側(cè)的界面波幾乎在管內(nèi)形成搭接,波峰向氣流流動方向彎曲直至=0.033 s時被氣相吹散,整個過程用時0.024 s,=0.057 s時波開始繼續(xù)生長并重復以上過程。由圖6可看出,大液相流速下管內(nèi)液膜厚度顯著增大,對比圖5也可看出管內(nèi)截面液相份額顯著增大,進而使氣相在管截面所占面積相應變小使得氣相流速增大,在氣流作用下易形成攪混流,對應圖6中液膜厚度的隨機增大或減小。攪混作用使液相入口處的波難以呈孤立波形式沿軸向及徑向發(fā)展。

    圖5 jl=0.121 m/s,jg=4.076 m/s時界面波形變化

    2.2 氣相流速的影響

    如圖7(a)所示,/0代表最大波幅所在軸向位置,其中為波峰位置在圓管豎直方向上的位置,0為最大波幅初始所在位置,這里假設界面波初始時是在液相入口處呈一個對稱分布的獨立波,因此將波峰最初位置定于液相入口中點處,即初始時0=68.95 mm。g=6.379 m/s時,最大波幅值隨時間變化有顯著增大,而最大波幅所在位置沿軸向變化很小,這說明在該氣相流速下界面波幾乎是在其初始產(chǎn)生的位置不斷增大且受氣流的拖曳力較小,界面波主要是沿徑向發(fā)展直至達到臨界波幅值時,波峰受到氣流作用沿氣流方向拉伸后被撕碎形成液滴進而被氣流帶走或重新在管壁上沉積。如圖7(b)及圖7(c)所示,g=7.000 m/s及7.721 m/s時,最大波幅較g=6.379 m/s工況減小且最大波幅所在軸向位置隨時間不斷變化,這說明隨著氣相流速的增大,界面波沿軸向的發(fā)展更加顯著且氣液界面更加平坦,波峰處被氣流撕碎形成液滴時對應的臨界波幅更小。而氣相流速越小界面波在徑向上的發(fā)展更顯著,即氣液界面波峰更顯著,波峰被撕碎時的臨界波幅更大。

    圖6 jl=0.0.121 m/s,jg=4.076 m/s時液膜在0.60 s內(nèi)的變化[10]

    2.3 流場分析

    2.3.1小液相流速

    如圖8所示,界面波的存在使得氣相流通面積變小從而導致氣相流速變大,因此界面波迎風側(cè)氣相流線較密,在波峰處速度分布變得平坦,而波峰后方為尾流低壓區(qū),這與Davide Del Col等[3]模擬結果一致,如圖9所示,壓力在界面波背風側(cè)突然下降,而在氣相流動的下游位置壓力則逐漸升高直至到達一常數(shù),這是因為氣流在遠離界面波的一段距離后擾動減小而不再受尾流區(qū)影響。界面波的背風側(cè)與迎風側(cè)之間的壓力差使界面波受到了向上的附加力作用,加速了液泛的發(fā)生。除此之外,隨著時間的推移,界面波在徑向上不斷生長使得壓差不斷增大,進而導致尾流區(qū)不斷增大直至波峰處被氣流撕碎,大尾流造成流動分離。

    圖7 jl=0.028 m/s時不同氣相流速下最大波幅值及其軸向位置隨時間變化圖

    圖8 jl=0.028 m/s時界面波附近流場圖:(a)jg=6.379 m/s(b)jg=7.000 m/s

    圖9 界面波達到最大波幅時的壓力變化[3]

    2.3.2大液相流速

    如圖10所示,界面波在初始形成時背風側(cè)也會形成尾渦,而在大液相流速下界面波在氣流作用下會在徑向上快速發(fā)展,使得背風側(cè)尾渦也快速增大進而導致波前后兩側(cè)壓差迅速變化,因此在大液相流速下界面波極易在管內(nèi)被氣流吹散形成攪混流,尾渦耗散后界面波重新生長并重復以上過程。大液相流速下由于管內(nèi)形成攪混流態(tài),管內(nèi)各截面上都有一定液相份額,隨著液相的不斷補充以及管內(nèi)擾動的存在,注水口處的波難以穩(wěn)定地呈孤立波生長。

    3 結論

    本文基于VOF模型,對=25 mm的圓管中氣液逆向流動過程中界面波的演變開展了CFD計算。得到如下結論:

    (1)小液相流速下界面波達到臨界波幅后,波峰在氣流作用下克服自身表面張力被撕碎;大液相流速下界面波發(fā)展迅速,難以孤立波穩(wěn)定生長且波幅較小,易形成攪混流。

    (2)小液相流速下隨著氣相流速的增大,界面波發(fā)展過程中界面的最大波幅不斷減小且最大波幅所在軸向位置沿氣流方向向下游發(fā)展。

    (3)小液相流速下界面波在受到氣流作用后波幅逐漸增大,波前后產(chǎn)生的壓差對界面波有附加力的作用,界面波在氣流及附加力雙重作用下變得不穩(wěn)定進而促使液泛發(fā)生。大液相流速下界面波沿徑向發(fā)展十分迅速,其背風側(cè)的尾渦快速增大從而導致界面波兩側(cè)的壓差也迅速增大,加大了波的不穩(wěn)定性。

    致謝

    本工作得到了國家重點研發(fā)計劃(2018YFB1900400)的支持。

    [1] Jayanti S,Tokarz A,Hewitt GF. Theoretical investigation of the diameter effect on flooding in countercurrent flow[J]. International Journal of Multiphase Flow,1996,22(2):307-324.

    [2] Shearer C J,Davidson. The investigation of a standing wave due to gas blowing upwards over a liquid film;its relation to flooding in wetted-wall columns[J]. Journal of fluid mechanics. 1965,22(2):321-335.

    [3] Enrico D R,Davide D C. Numerical simulation of churn flow in a vertical pipe[J]. Chemical engineering science,2009,64(17):3753-3765.

    [4] Hirt C W,Nichols B D. Volume of fluid(VOF)method for the dynamics of free boundaries[J]. Journal of computational physics,1981,39(1):201-225.

    [5] Barbosa J F,Lima E C,Zan?o R A,et al. The use of a W-Rh permanent modifier for direct determination of bismuth in urine and whole blood by electrothermal atomic absorption spectrometry[J]. Journal of Analytical Atomic Spectrometry,2001,16(8):842-846.

    [6] Murase M,Kinoshita I,Utanohara Y,et al. Numerical calculations for air-water tests on CCFL in different-scale models of a PWR hot leg;International Conference on Nuclear Engineering,2010.

    [7] Minami N,Murase M,Tomiyama A. Countercurrent Gas-Liquid Flow in a PWR Hot Leg under Reflux Cooling(Ⅱ)Numerical Simulation of 1/15 Scale Air-Water Tests[J]. Journal of Nuclear Science Technology,2010,47(2):149-155.

    [8] Minami N,Nishiwaki D,Nariai T,et al. Countercurrent gas-liquid flow in a PWR hot leg under reflux cooling(Ⅰ)air-water tests for 1/15-scale model of a PWR hot leg[J]. Journal of Nuclear Science Technology,2010,47(2):142-148.

    [9] Utanohara Y.,Kinoshita I.,Murase M.,et al. Numerical simulation using CFD software of countercurrent gas–liquid flow in a PWR hot leg under reflux condensation[J]. Nuclear Engineering Design,2011,241(5):1643-1655.

    [10] Wan J,Sun W,Deng J,et al. Experimental study on air-water countercurrent flow limitation in a vertical tube based on measurement of film thickness behavior[J]. Nuclear Engineering Technology,2021,53(6):1821-1833.

    Study on the Interfacial Behavior of the Air/Water Counter-current Flow in the Vertical Tube Based on the CFD

    WAN Jie1,2,SUN Wan1,*,LI Xiang1,2,ZHU Longxiang1,2,MA Zaiyong1,2,ZHANG Luteng1,2, PAN Liangming1,2

    (1. Key Laboratory of Low-Grade Energy Utilization Technologies and Systems,Ministry of Education,Chongqing University,Chongqing 400044,China;2. Department of Nuclear Engineering and Technology,Chongqing University,Chongqing 400044,China)

    In order to study the local characteristics of the air-water counter-current flow, a numerical simulation of interfacial wave in a vertical tube with 25 mm diameter has been carried out based on the VOF (volume of fluid) model. The change of interfacial wave was analyzed and showed a good agreement with the experimental film thickness. Through the analysis of the variation of the flow field under different liquid flow rates, it can be concluded that the pressure difference generated between the windward and leeward of the interfacial wave introduces an additional force, which weakens the stability of the interfacial wave. And the pressure difference increases with the increase of the liquid flow rate. Moreover, the greater the gas flow rate, the faster the amplitude of the interfacial wave develops along the axial direction, indicating that the instability of the interfacial wave increases with the increase of the gas flow rate. Consequently, the increase of the gas and liquid flow rate reduces the stability of the interfacial wave, namely, the flooding is more likely to occur.

    CFD; Flooding; VOF; Interfacial wave; Liquid film thickness

    TL334

    A

    0258-0918(2023)05-1105-08

    2022-08-04

    國家重點研發(fā)計劃資助項目(2018YFB1900400)

    萬 潔(1996—),女,貴州三穗人,博士研究生,現(xiàn)主要從事兩相逆向流動極限方面研究

    孫 皖,E-mail:sunwan@cqu.edu.cn

    猜你喜歡
    波峰波幅氣相
    開封市健康人群面神經(jīng)分支復合肌肉動作電位波幅分布范圍研究
    氣相過渡金屬鈦-碳鏈團簇的研究
    作用于直立堤墻與樁柱的波峰高度分析計算
    開不同位置方形洞口波紋鋼板剪力墻抗側(cè)性能
    價值工程(2018年25期)2018-09-26 10:14:36
    兒童標準12導聯(lián)T波峰末間期的分析
    考慮傳輸函數(shù)特性的行波幅值比較式縱聯(lián)保護原理
    新型釩基催化劑催化降解氣相二噁英
    頻率偏移時基波幅值計算誤差對保護的影響及其改進算法
    預縮聚反應器氣相管“鼓泡”的成因探討
    Dynamic Loads and Wake Prediction for Large Wind Turbines Based on Free Wake Method
    在线免费观看的www视频| 哪里可以看免费的av片| av在线观看视频网站免费| .国产精品久久| 亚洲一区二区三区不卡视频| 色尼玛亚洲综合影院| 热99re8久久精品国产| 色播亚洲综合网| 国产高清视频在线观看网站| 丰满的人妻完整版| 无遮挡黄片免费观看| 免费在线观看影片大全网站| 黄色日韩在线| 18+在线观看网站| 日本在线视频免费播放| 免费看日本二区| 热99re8久久精品国产| 国产亚洲精品久久久com| 亚洲av熟女| 综合色av麻豆| 有码 亚洲区| 久久精品国产亚洲av香蕉五月| 亚洲真实伦在线观看| 国内少妇人妻偷人精品xxx网站| 变态另类成人亚洲欧美熟女| 老司机深夜福利视频在线观看| 99久久精品一区二区三区| 国产 一区 欧美 日韩| 免费看日本二区| 亚洲人成伊人成综合网2020| 久久性视频一级片| 男人和女人高潮做爰伦理| 一级毛片久久久久久久久女| 麻豆国产97在线/欧美| 亚洲av成人av| 亚洲美女视频黄频| 3wmmmm亚洲av在线观看| 人人妻人人澡欧美一区二区| 久久精品人妻少妇| 12—13女人毛片做爰片一| 欧美区成人在线视频| 成人高潮视频无遮挡免费网站| 亚洲午夜理论影院| 欧美最黄视频在线播放免费| 欧美另类亚洲清纯唯美| 最新在线观看一区二区三区| 99在线视频只有这里精品首页| 欧美成人一区二区免费高清观看| 国产精品女同一区二区软件 | 男女下面进入的视频免费午夜| 一夜夜www| 床上黄色一级片| 久久中文看片网| 热99在线观看视频| 少妇人妻一区二区三区视频| 日韩成人在线观看一区二区三区| 欧美+亚洲+日韩+国产| 久久香蕉精品热| 免费av不卡在线播放| 国产乱人伦免费视频| 91九色精品人成在线观看| 午夜日韩欧美国产| 99精品在免费线老司机午夜| 久久国产精品人妻蜜桃| 亚洲avbb在线观看| 欧美中文日本在线观看视频| 精品一区二区三区视频在线| 岛国在线免费视频观看| 色吧在线观看| 男人舔女人下体高潮全视频| 亚洲国产精品999在线| 精品人妻1区二区| 不卡一级毛片| 三级国产精品欧美在线观看| 国产亚洲精品久久久久久毛片| 亚洲中文字幕日韩| 女同久久另类99精品国产91| 国产精品98久久久久久宅男小说| 国产白丝娇喘喷水9色精品| 757午夜福利合集在线观看| 久久久久久久久久成人| 欧美日韩国产亚洲二区| 国产熟女xx| 91狼人影院| 国产精品不卡视频一区二区 | 精品午夜福利在线看| 免费高清视频大片| 免费高清视频大片| 亚洲无线在线观看| 我要搜黄色片| 天美传媒精品一区二区| 欧美激情在线99| 欧美xxxx黑人xx丫x性爽| 亚洲最大成人手机在线| 亚洲人成伊人成综合网2020| 老司机午夜福利在线观看视频| 蜜桃久久精品国产亚洲av| 日韩欧美国产在线观看| 淫妇啪啪啪对白视频| 丁香六月欧美| 91九色精品人成在线观看| 黄色配什么色好看| 国产精品1区2区在线观看.| 熟女电影av网| 国产精品综合久久久久久久免费| 久久欧美精品欧美久久欧美| 亚洲欧美日韩无卡精品| 黄色女人牲交| 亚洲av美国av| 亚洲精品粉嫩美女一区| 国产精品自产拍在线观看55亚洲| 欧美性感艳星| 亚洲欧美日韩卡通动漫| 一区二区三区免费毛片| 成年人黄色毛片网站| 国产麻豆成人av免费视频| 成年免费大片在线观看| 色av中文字幕| 亚洲av成人精品一区久久| 脱女人内裤的视频| 9191精品国产免费久久| 国内精品一区二区在线观看| 国产黄片美女视频| 欧美三级亚洲精品| ponron亚洲| 精品人妻1区二区| 国产男靠女视频免费网站| 三级国产精品欧美在线观看| 91在线观看av| 欧美成人性av电影在线观看| 熟女电影av网| 日韩精品中文字幕看吧| 又黄又爽又免费观看的视频| .国产精品久久| 少妇高潮的动态图| 九色国产91popny在线| 国产综合懂色| 精品久久久久久久人妻蜜臀av| 高清在线国产一区| 久久精品国产99精品国产亚洲性色| 能在线免费观看的黄片| 免费观看的影片在线观看| 性插视频无遮挡在线免费观看| 国产精品爽爽va在线观看网站| 国产av不卡久久| 五月玫瑰六月丁香| 国产精品精品国产色婷婷| 88av欧美| 久久精品夜夜夜夜夜久久蜜豆| 国产av在哪里看| 我的老师免费观看完整版| 永久网站在线| 久久精品91蜜桃| 大型黄色视频在线免费观看| 又紧又爽又黄一区二区| 禁无遮挡网站| 我要看日韩黄色一级片| 久久亚洲精品不卡| 国产精品影院久久| 日本免费一区二区三区高清不卡| 在线天堂最新版资源| 免费无遮挡裸体视频| 在线十欧美十亚洲十日本专区| 欧美最黄视频在线播放免费| 99热这里只有精品一区| 亚洲精品456在线播放app | 国产高清视频在线播放一区| 国产野战对白在线观看| 欧美+亚洲+日韩+国产| 日韩免费av在线播放| 亚洲成人免费电影在线观看| 看十八女毛片水多多多| 成人精品一区二区免费| 99热这里只有是精品50| 亚洲avbb在线观看| 日本黄色片子视频| 首页视频小说图片口味搜索| 99久国产av精品| 亚洲精品456在线播放app | 国产野战对白在线观看| 久久国产乱子免费精品| 神马国产精品三级电影在线观看| 欧美一区二区精品小视频在线| 国产午夜精品论理片| 91九色精品人成在线观看| 在线看三级毛片| 国产色爽女视频免费观看| 欧美另类亚洲清纯唯美| 婷婷丁香在线五月| 国产精品一区二区免费欧美| av在线老鸭窝| 啦啦啦韩国在线观看视频| 亚洲在线自拍视频| 别揉我奶头 嗯啊视频| 又爽又黄a免费视频| 又爽又黄a免费视频| 欧美国产日韩亚洲一区| 日韩欧美在线二视频| 日韩欧美在线二视频| 欧美另类亚洲清纯唯美| 成熟少妇高潮喷水视频| 极品教师在线视频| 久久天躁狠狠躁夜夜2o2o| 欧美黄色片欧美黄色片| 五月玫瑰六月丁香| 国产精品,欧美在线| 亚洲激情在线av| 国产蜜桃级精品一区二区三区| 成人特级黄色片久久久久久久| 精品一区二区三区视频在线观看免费| 一级黄色大片毛片| 欧美日韩国产亚洲二区| 天天一区二区日本电影三级| 久久久久久久久久成人| 色5月婷婷丁香| 好看av亚洲va欧美ⅴa在| 国产精品一区二区免费欧美| 99久久精品一区二区三区| 99精品久久久久人妻精品| 亚洲一区高清亚洲精品| 国产精品99久久久久久久久| 亚洲中文字幕日韩| 精品久久久久久久末码| 丝袜美腿在线中文| 国产三级黄色录像| 一a级毛片在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲国产欧美人成| 有码 亚洲区| or卡值多少钱| 制服丝袜大香蕉在线| 嫩草影院新地址| 欧美成人一区二区免费高清观看| 欧美乱妇无乱码| 简卡轻食公司| 无人区码免费观看不卡| 国产成人aa在线观看| 久久久国产成人免费| 精品熟女少妇八av免费久了| 非洲黑人性xxxx精品又粗又长| 国产精品伦人一区二区| 蜜桃久久精品国产亚洲av| 丰满人妻熟妇乱又伦精品不卡| 蜜桃久久精品国产亚洲av| 99热这里只有精品一区| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av天美| 一级av片app| 亚洲国产精品久久男人天堂| 在线播放国产精品三级| 成人亚洲精品av一区二区| 国产v大片淫在线免费观看| 成熟少妇高潮喷水视频| 国内久久婷婷六月综合欲色啪| 亚洲国产精品久久男人天堂| 天堂影院成人在线观看| 国产成人a区在线观看| av专区在线播放| 亚洲va日本ⅴa欧美va伊人久久| av天堂中文字幕网| 熟女电影av网| 99国产极品粉嫩在线观看| 亚洲熟妇熟女久久| 黄色视频,在线免费观看| 亚洲欧美日韩卡通动漫| 高清毛片免费观看视频网站| 亚洲国产精品成人综合色| 色精品久久人妻99蜜桃| 免费观看的影片在线观看| 久99久视频精品免费| 国产精品伦人一区二区| 亚洲不卡免费看| 亚洲在线观看片| 国产成+人综合+亚洲专区| 在线a可以看的网站| 国产高清视频在线播放一区| 少妇人妻一区二区三区视频| 久久精品国产亚洲av香蕉五月| 国产成年人精品一区二区| 深爱激情五月婷婷| 97热精品久久久久久| bbb黄色大片| 午夜精品久久久久久毛片777| 国产成人av教育| 成人鲁丝片一二三区免费| 成熟少妇高潮喷水视频| 99久久成人亚洲精品观看| 黄色配什么色好看| ponron亚洲| 五月伊人婷婷丁香| 久久久久久久久久黄片| 少妇裸体淫交视频免费看高清| 午夜老司机福利剧场| 18禁黄网站禁片午夜丰满| 亚洲av不卡在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产成人av教育| 校园春色视频在线观看| 亚洲人成网站高清观看| 亚洲av成人精品一区久久| 伊人久久精品亚洲午夜| 亚洲人成网站在线播放欧美日韩| 日韩欧美国产一区二区入口| 九九久久精品国产亚洲av麻豆| 午夜a级毛片| 精品久久国产蜜桃| 狠狠狠狠99中文字幕| 国产欧美日韩一区二区精品| 精品无人区乱码1区二区| 欧洲精品卡2卡3卡4卡5卡区| 啦啦啦观看免费观看视频高清| 国产精品99久久久久久久久| 成人特级黄色片久久久久久久| 精品人妻1区二区| 色噜噜av男人的天堂激情| 亚州av有码| 琪琪午夜伦伦电影理论片6080| 免费人成视频x8x8入口观看| 一本一本综合久久| 午夜福利成人在线免费观看| 亚洲不卡免费看| 日韩成人在线观看一区二区三区| 亚洲精品色激情综合| 欧美性猛交╳xxx乱大交人| 免费大片18禁| 桃红色精品国产亚洲av| 国产69精品久久久久777片| 亚洲不卡免费看| 99久久99久久久精品蜜桃| av福利片在线观看| 午夜影院日韩av| 我的女老师完整版在线观看| 18禁在线播放成人免费| 在现免费观看毛片| 性欧美人与动物交配| 中文字幕久久专区| 亚洲av五月六月丁香网| 少妇被粗大猛烈的视频| 少妇的逼好多水| 中亚洲国语对白在线视频| 国产真实乱freesex| 久久草成人影院| 成人欧美大片| 变态另类成人亚洲欧美熟女| 亚洲真实伦在线观看| 91麻豆精品激情在线观看国产| 免费看美女性在线毛片视频| 国模一区二区三区四区视频| 欧美色欧美亚洲另类二区| 精品久久久久久久久久久久久| 搡老熟女国产l中国老女人| av在线蜜桃| 久久精品国产亚洲av香蕉五月| 九九热线精品视视频播放| 色播亚洲综合网| 国产精品嫩草影院av在线观看 | 淫妇啪啪啪对白视频| 色5月婷婷丁香| 欧美一区二区精品小视频在线| 我的女老师完整版在线观看| 日本免费a在线| 国产成人av教育| 久久久精品大字幕| 免费人成在线观看视频色| 免费av不卡在线播放| 国产野战对白在线观看| 精品午夜福利在线看| 亚洲精品影视一区二区三区av| 夜夜夜夜夜久久久久| 天堂网av新在线| 白带黄色成豆腐渣| 亚洲精品在线观看二区| 久久这里只有精品中国| 国内毛片毛片毛片毛片毛片| 欧美3d第一页| 日韩欧美在线乱码| 我的老师免费观看完整版| 热99在线观看视频| 日日干狠狠操夜夜爽| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久热精品热| 波多野结衣高清作品| 一个人观看的视频www高清免费观看| 中文字幕精品亚洲无线码一区| 精品久久久久久久久久久久久| 亚洲欧美日韩东京热| 国产三级中文精品| 香蕉av资源在线| 99热只有精品国产| 丝袜美腿在线中文| 极品教师在线视频| 小蜜桃在线观看免费完整版高清| 亚洲一区二区三区色噜噜| 国产私拍福利视频在线观看| 中亚洲国语对白在线视频| 亚洲欧美精品综合久久99| 99国产精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 天堂动漫精品| 国产日本99.免费观看| 99热精品在线国产| 亚洲成人免费电影在线观看| 久久这里只有精品中国| 欧美激情在线99| 亚洲成a人片在线一区二区| a级一级毛片免费在线观看| 亚洲中文字幕日韩| www日本黄色视频网| 国产精品爽爽va在线观看网站| 午夜精品在线福利| 尤物成人国产欧美一区二区三区| 成人一区二区视频在线观看| 在线观看舔阴道视频| 亚洲 欧美 日韩 在线 免费| 欧美在线黄色| 欧美性感艳星| 色播亚洲综合网| 久久草成人影院| 内射极品少妇av片p| 亚洲国产精品999在线| 国产精品自产拍在线观看55亚洲| av女优亚洲男人天堂| 一本精品99久久精品77| 婷婷精品国产亚洲av| 97人妻精品一区二区三区麻豆| 男女之事视频高清在线观看| 99久久精品国产亚洲精品| av专区在线播放| 国产精品一区二区性色av| 一区福利在线观看| 欧美zozozo另类| 99久国产av精品| 99国产精品一区二区三区| 亚洲美女黄片视频| 噜噜噜噜噜久久久久久91| 乱人视频在线观看| 熟女人妻精品中文字幕| 亚洲欧美日韩高清在线视频| 1000部很黄的大片| 男人狂女人下面高潮的视频| 在线十欧美十亚洲十日本专区| 精品国产三级普通话版| 啪啪无遮挡十八禁网站| 精品午夜福利视频在线观看一区| 日本三级黄在线观看| 少妇熟女aⅴ在线视频| 最好的美女福利视频网| 欧美性感艳星| 国产亚洲欧美98| 在线国产一区二区在线| 成人特级黄色片久久久久久久| 国产欧美日韩精品亚洲av| 青草久久国产| 九色国产91popny在线| 草草在线视频免费看| av在线天堂中文字幕| 亚洲人成伊人成综合网2020| 久久久久久久久中文| 欧美一区二区精品小视频在线| 亚洲av成人av| 最近在线观看免费完整版| 午夜视频国产福利| 日本五十路高清| 精品一区二区三区视频在线| 国产精品电影一区二区三区| 日韩人妻高清精品专区| 国产精品一区二区三区四区久久| 欧美日韩福利视频一区二区| 国产欧美日韩精品一区二区| 一进一出抽搐gif免费好疼| 国产一区二区激情短视频| 午夜激情欧美在线| 天天一区二区日本电影三级| 少妇被粗大猛烈的视频| 一个人观看的视频www高清免费观看| 婷婷丁香在线五月| 久久久久国产精品人妻aⅴ院| 国产精品av视频在线免费观看| 欧美成人一区二区免费高清观看| 日韩欧美在线二视频| 国产欧美日韩精品一区二区| 一进一出抽搐gif免费好疼| 亚洲内射少妇av| 欧美黄色淫秽网站| 观看免费一级毛片| 成人特级黄色片久久久久久久| 欧美午夜高清在线| www.www免费av| 久久久精品大字幕| 亚洲一区二区三区色噜噜| 国产av不卡久久| 国产精品永久免费网站| 成人特级黄色片久久久久久久| 亚洲自拍偷在线| 性色av乱码一区二区三区2| 国产色婷婷99| 久久久久亚洲av毛片大全| 一区二区三区四区激情视频 | 亚洲美女搞黄在线观看 | 一区二区三区高清视频在线| 制服丝袜大香蕉在线| 欧美极品一区二区三区四区| 国产精品98久久久久久宅男小说| 丰满的人妻完整版| 午夜精品久久久久久毛片777| 亚洲第一电影网av| 悠悠久久av| 免费电影在线观看免费观看| 少妇人妻精品综合一区二区 | 亚洲av五月六月丁香网| 丰满的人妻完整版| 久久九九热精品免费| 怎么达到女性高潮| 大型黄色视频在线免费观看| 亚洲av电影在线进入| 成人鲁丝片一二三区免费| 欧美精品啪啪一区二区三区| 99久久精品热视频| 中亚洲国语对白在线视频| 一本精品99久久精品77| 国产午夜福利久久久久久| 亚洲熟妇熟女久久| aaaaa片日本免费| 亚洲va日本ⅴa欧美va伊人久久| 高清日韩中文字幕在线| 欧美日韩福利视频一区二区| 国产成+人综合+亚洲专区| 99久久精品热视频| 婷婷丁香在线五月| 麻豆成人av在线观看| 在线免费观看的www视频| 99精品久久久久人妻精品| 日韩国内少妇激情av| 日韩欧美三级三区| 51国产日韩欧美| 亚洲中文日韩欧美视频| 精品久久久久久,| 国产成人影院久久av| 免费人成在线观看视频色| 搡老岳熟女国产| 国产精华一区二区三区| 午夜视频国产福利| 国产午夜精品久久久久久一区二区三区 | 深夜a级毛片| 欧美成人一区二区免费高清观看| 观看免费一级毛片| 婷婷精品国产亚洲av在线| 丝袜美腿在线中文| av黄色大香蕉| 久久欧美精品欧美久久欧美| 亚洲av第一区精品v没综合| 最近最新中文字幕大全电影3| 十八禁国产超污无遮挡网站| 国产成年人精品一区二区| 免费一级毛片在线播放高清视频| 哪里可以看免费的av片| 怎么达到女性高潮| 亚洲国产欧洲综合997久久,| 久久久久久久亚洲中文字幕 | 精品熟女少妇八av免费久了| 一区二区三区激情视频| 嫁个100分男人电影在线观看| 亚洲人成网站在线播放欧美日韩| 久久天躁狠狠躁夜夜2o2o| 网址你懂的国产日韩在线| 亚洲欧美日韩卡通动漫| 国产免费男女视频| 观看美女的网站| 国产探花极品一区二区| 中文字幕熟女人妻在线| 色综合婷婷激情| 淫秽高清视频在线观看| 国产在线男女| 久久久久久久午夜电影| 欧美3d第一页| 在线观看66精品国产| 国产精品日韩av在线免费观看| 人妻久久中文字幕网| 欧美一区二区精品小视频在线| 熟女人妻精品中文字幕| 国产精品98久久久久久宅男小说| 一个人观看的视频www高清免费观看| 看免费av毛片| 日韩中字成人| 99精品久久久久人妻精品| 亚洲欧美日韩东京热| 国产精品三级大全| 亚洲av不卡在线观看| 每晚都被弄得嗷嗷叫到高潮| 我要搜黄色片| 国内精品久久久久精免费| 国产精品1区2区在线观看.| 女生性感内裤真人,穿戴方法视频| 久久99热6这里只有精品| 婷婷精品国产亚洲av| 露出奶头的视频| 欧美国产日韩亚洲一区| 精品人妻偷拍中文字幕| 日韩av在线大香蕉| 国产精品电影一区二区三区| 深夜a级毛片| 欧美成人一区二区免费高清观看| 91狼人影院| 日韩欧美一区二区三区在线观看| 97超级碰碰碰精品色视频在线观看| 内射极品少妇av片p| 99热只有精品国产| 久久精品国产99精品国产亚洲性色| 小说图片视频综合网站| 90打野战视频偷拍视频| 老女人水多毛片| av在线观看视频网站免费| 一本久久中文字幕| 亚洲无线观看免费| 日韩欧美三级三区|