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

    結(jié)構(gòu)形式對擋板出口截面流體力學(xué)性能的影響

    2020-08-24 04:23:00張海春周圓圓吳劍華
    沈陽化工大學(xué)學(xué)報 2020年2期
    關(guān)鍵詞:渦量柱面示蹤劑

    張海春, 周圓圓, 龔 斌, 吳劍華, 張 靜

    (1.長春中車軌道車輛有限公司, 吉林 長春 130052;2.沈陽化工大學(xué) 遼寧省化工新技術(shù)轉(zhuǎn)移推廣中心, 遼寧 沈陽 110142)

    車輛清洗是提高軌道服務(wù)質(zhì)量的一項(xiàng)重要措施.洗車污水的濁度較高,含有大量的較細(xì)泥沙顆粒.一般要對沙水進(jìn)行重力分離,以達(dá)到污水排放標(biāo)準(zhǔn).分層器是一種應(yīng)用廣泛的沙水重力分離設(shè)備,通常會在分層器入口添加防沖擊結(jié)構(gòu),吸收流體介質(zhì)的動能,改善分層器內(nèi)已分層物料流場的分布狀況.在非均相重力分離的實(shí)際應(yīng)用中,進(jìn)料口處設(shè)置擋板可以明顯減小流場內(nèi)的回流面積,增大有效流動區(qū)域,提高分離效率,這為分層器的設(shè)計(jì)和運(yùn)行提供了參考[1].擋板的最佳位置是在回流區(qū)域,可用來破壞回流區(qū),提高分層器的運(yùn)行效率[2].在入口布置一塊徑向擋板和兩塊軸向擋板,以三塊擋板構(gòu)成的出口截面作為分層器主流場入口,其流動狀態(tài)對分離效率至關(guān)重要.張靜等[3]運(yùn)用模擬的方法對擋板局域內(nèi)沖擊射流受擋板曲率和沖擊間距的影響進(jìn)行分析,預(yù)測了滯點(diǎn)壓力和沖擊力,并利用速度場解釋了局域壓力場成因.然而,對于擋板局域和分層器內(nèi)關(guān)聯(lián)的重要環(huán)節(jié)——擋板出口截面沒有進(jìn)行詳盡分析.

    沖擊射流具有獨(dú)特的流動現(xiàn)象.Cooper[4]、Craft[5]和Ashforth-Frost等[6]對沖擊射流的流動結(jié)構(gòu)進(jìn)行了精確的實(shí)驗(yàn)測量和數(shù)值計(jì)算,發(fā)現(xiàn)在剪切層和噴嘴入口附近產(chǎn)生漩渦,其形成與射流出口雷諾數(shù)、是否受限和沖擊高度相關(guān).Ozmen等[7]利用熱線風(fēng)速儀研究發(fā)現(xiàn),近壁面徑向湍流強(qiáng)度存在兩個峰值,徑向擋板徑向壓力系數(shù)分布與沖擊寬度相關(guān).毛軍逵等[8]通過煙線流動可視化方法研究了柱面擋板,其中,凹柱面隨沖擊雷諾數(shù)Re的增加,在沖擊滯止區(qū)域兩側(cè)形成穩(wěn)定的旋流結(jié)構(gòu),高壓區(qū)逐漸增強(qiáng);而凹柱面上的流動不穩(wěn)定性降低并易形成回流[9-10]. Madarame 等[11]研究表明,對于向上平面射流撞擊自由面的自感振蕩,脈線法是CFD模擬的有效驗(yàn)證方法和可視化補(bǔ)充.Torré等[12-13]利用脈線法表示沖擊射流軌跡,實(shí)驗(yàn)和CFD結(jié)果吻合度較高,使用統(tǒng)計(jì)分析建立的相關(guān)性表明,對于橫流模型中的射流,脈線法能很好地描述射流軌跡.

    針對擋板局域范圍內(nèi)沖擊射流受限流動,本文采用脈線法與CFD模擬相結(jié)合,對擋板局域示蹤劑濃度及出口截面上流場分布特性進(jìn)行分析.通過示蹤劑對沖擊射流軌跡的可視化演示,探討擋板作用下受限沖擊射流流動特性.分析凹柱面(K>0)、平面(K=0)、凸柱面(K< 0)三種擋板形式對擋板局域受限射流出口截面上流體力學(xué)性能的影響,為立式圓筒體分層器內(nèi)非均相分離提供基礎(chǔ)研究數(shù)據(jù).

    1 數(shù)值模擬

    1.1 物理模型

    為分析擋板形式對擋板內(nèi)及出口局域示蹤劑濃度場的影響,本文的物理模型參照文獻(xiàn)[3]所示結(jié)構(gòu),附加示蹤劑加料管,如圖1所示.示蹤劑加料管為壁厚1 mm圓管,出口與主流體入口平齊.分層器抽象模型具體結(jié)構(gòu)參數(shù)見表1.其中,擋板相對曲率定義為

    表1 分層器抽象模型結(jié)構(gòu)參數(shù)

    圖1 分層器抽象模型簡圖

    K=d/2R.

    (1)

    1.2 網(wǎng)格分布

    本文是在文獻(xiàn)[3]的基礎(chǔ)上作進(jìn)一步研究,網(wǎng)格方案、邊界條件及模型選擇均沿用該文獻(xiàn)的結(jié)論.立式分層器內(nèi)流場整體劃分邊長為8 mm的四面體非結(jié)構(gòu)化網(wǎng)格.利用四面體網(wǎng)格邊上增加節(jié)點(diǎn)的方法,對入口局域長×寬×高=120 mm×120 mm×120 mm的正方體區(qū)域進(jìn)行網(wǎng)格加密.對平擋板結(jié)構(gòu)進(jìn)行網(wǎng)格無關(guān)性研究表明,未加密、一次加密、二次加密網(wǎng)格與三次加密網(wǎng)格的入口中心線上速度偏差分別為14.7 %、5.4 %、2.3 %.選擇二次加密網(wǎng)格作為網(wǎng)格劃分方案,三種擋板結(jié)構(gòu)水平對稱面上網(wǎng)格分布及對應(yīng)單元網(wǎng)格體積如圖2所示.計(jì)算域內(nèi)總網(wǎng)格為240~279萬,細(xì)化區(qū)域單元網(wǎng)格體積范圍在0.338~2.47 mm3.

    圖2 網(wǎng)格分析

    1.3 模型選擇及邊界條件

    研究介質(zhì)為不可壓縮流體水(ρ=998.2 kg/m3,μ=1.003×10-3Pa·s),湍流雷諾數(shù)范圍為Re=3 000~11 000,定性尺寸為入口管直徑d=30 mm.采用Realizablek-ε湍流模型對研究對象進(jìn)行模擬計(jì)算.入口為均勻速度入口,所有壁面為光滑無滑移壁面,圓筒體上下截面設(shè)置為壓力出口邊界,上下截面表壓分別為100 mm(979.2 Pa)和1 000 mm(9 792.3 Pa)水柱.

    為了與實(shí)驗(yàn)做對比研究,對Re=7 000(uin=0.234 m/s)條件下進(jìn)行示蹤劑擴(kuò)散數(shù)值模擬.其中,水循環(huán)入口中心位置設(shè)置直徑為1.0 mm的圓截面作為示蹤劑入口,在穩(wěn)態(tài)單相流計(jì)算過程中,該圓截面設(shè)置為壁面,數(shù)值計(jì)算的收斂極限設(shè)置為10-4.穩(wěn)態(tài)計(jì)算收斂結(jié)束,模型修改為非穩(wěn)態(tài)兩相流混合模型.基礎(chǔ)相和第二相(示蹤劑)均為不可壓縮水,直徑為1.0 mm的圓截面修改為速度入口,速度為0.234 m/s,第二相含率為1.第二相離散采用二階迎風(fēng)格式,時間步為0.001 s,最大迭代次數(shù)20次,每隔1 s保存一個計(jì)算結(jié)果.

    2 實(shí)驗(yàn)驗(yàn)證

    2.1 實(shí)驗(yàn)裝置

    為減小分層器出口對流體流動的影響,將圖1結(jié)構(gòu)放置于1 000 mm×1 000 mm×1 000 mm方箱體中,上下留出100 mm空間,如圖3(a)所示,使流體出口設(shè)置在圓筒體上下截面上,對擋板局域流體影響最小.為節(jié)約成本,三種形式徑向擋板利用UPBOX 3D打印機(jī)制作,并增加護(hù)板,擋板與護(hù)板如圖3(b)所示.在安裝三種徑向擋板時,通過護(hù)板保持出口截面寬度b=35 mm,而且保證徑向擋板受流體沖擊時位置穩(wěn)定.

    2.2 實(shí)驗(yàn)流程及方法

    實(shí)驗(yàn)流程由水循環(huán)系統(tǒng)、示蹤劑添加系統(tǒng)、拍攝系統(tǒng)3部分組成,如圖4所示.水循環(huán)系統(tǒng)中,先在箱體中注滿水,再啟動循環(huán)泵,閥門全開以排凈循環(huán)系統(tǒng)中空氣.調(diào)整閥門,控制入口流量.在方箱外側(cè)局域出口對應(yīng)位置,下軸向擋板外表面均貼附亞光白紙,降低背景對拍攝的影響,為后處理提供最佳圖像.同時,出口側(cè)上軸向擋板和護(hù)板側(cè)面貼附刻度尺[如圖3(b)],為定量判定示蹤劑分布提供輔助參考.

    圖4 實(shí)驗(yàn)流程裝置

    選擇酸性墨水藍(lán)溶液作為示蹤劑.示蹤劑添加系統(tǒng)由示蹤劑容器、蠕動泵、輸送管線組成.水循環(huán)系統(tǒng)流量穩(wěn)定后,啟動蠕動泵,流量控制在10 mL/min,與模擬示蹤劑速度保持一致,水循環(huán)系統(tǒng)流量計(jì)測量誤差為5‰.一定時間后可見示蹤劑流出,關(guān)閉蠕動泵,等待拍攝.

    拍攝系統(tǒng)由兩個攝像機(jī)及數(shù)據(jù)采集系統(tǒng)組成,相機(jī)1得到Z軸法向面圖像,相機(jī)2得到X軸法向面即擋板出口截面圖像,由軟件控制兩個攝像機(jī)同步.先啟動攝像機(jī),再啟動蠕動泵,連續(xù)拍攝.對拍攝圖像取數(shù)值模擬對應(yīng)時刻進(jìn)行后處理.

    2.3 模擬與實(shí)驗(yàn)比較

    圖5對入口Re=7 000,t=1 s,2 s,3 s,4 s四個時間點(diǎn)的實(shí)驗(yàn)和CFD模擬的示蹤劑濃度分布作出描述.從X軸法向面圖像可以看出:示蹤劑經(jīng)歷了噴出t=1 s→撞擊徑向擋板t=2 s→向上下軸向擋板流動,受軸向擋板限制沿平擋板表面向圓筒體壁面流動t=3 s→遇壁面再向入口流動,與入口流體匯合再流向徑向擋板t=4 s.從三種類型擋板t=4 s圖像可以看出:擋板內(nèi)流體以入口水平對稱面為基準(zhǔn)形成上下兩個對稱的旋流結(jié)構(gòu),在擋板沖擊區(qū)流動趨勢相同,壁面射流區(qū)差異不明顯.濃度分布規(guī)律與文獻(xiàn)[3]對速度場的研究結(jié)論一致.

    圖5 擋板出口截面示蹤劑濃度分布(Re=7 000)

    圖5中多個時間點(diǎn)示蹤劑濃度分布云圖與CFD模擬結(jié)果類似.運(yùn)用matlab對實(shí)驗(yàn)圖像進(jìn)行處理,不同濃度所占面積的數(shù)值模擬與實(shí)驗(yàn)數(shù)據(jù)吻合程度達(dá)到90 %,從側(cè)面說明了該模擬方案對擋板內(nèi)流場預(yù)測的可行性.

    3 擋板形式對出口截面流動特性的影響

    擋板出口截面是流體在擋板內(nèi)流動與分離器內(nèi)流動的分界面,該截面上的流動特性直接決定了分離器內(nèi)流體的流動特性.現(xiàn)針對不同結(jié)構(gòu)擋板進(jìn)行數(shù)值模擬分析,深入探究出口截面上流動參數(shù)的差異,為分離器入口擋板設(shè)計(jì)提供理論依據(jù).

    3.1 靜壓力

    對三種形式徑向擋板出口截面靜壓力進(jìn)行分析,定義靜壓力系數(shù)為

    (2)

    圖6描述了擋板形式對出口截面靜壓力系數(shù)分布的影響.在出口截面上,以水平對稱面為基準(zhǔn)形成上下兩個對稱的負(fù)壓區(qū)域.通過比較發(fā)現(xiàn):隨著徑向擋板曲率減小,負(fù)壓中心向上下兩側(cè)軸向擋板趨近.其中,K=0.1的Cp≤-0.1區(qū)域面積最大,占出口截面總面積的69.3 %;K=0面積最小,Cp≤-0.1區(qū)域面積占總面積的48.4 %;K=-0.1的Cp≤-0.1區(qū)域面積占總面積的57.5 %.出口截面壓力低于設(shè)備入口壓力的主要原因是流動受上下軸向擋板的限制,形成對稱旋流結(jié)構(gòu),出口截面上兩個對稱的低壓區(qū)成為回流形成的主要誘因.

    圖6 靜壓力系數(shù)分布(Re=7 000)

    3.2 回流比

    在立式圓筒體分層器中,物料在徑向擋板的作用下,進(jìn)入罐體后形成壁面射流.擋板出口截面上速度分布決定了罐體內(nèi)壁面射流流動結(jié)構(gòu),是分層效率的關(guān)鍵因素.不同形式擋板內(nèi)沖擊射流的差異決定了出口截面上的流場特性.

    圖7描述了三種形式擋板出口截面上的法向速度分布,以水平對稱面為基準(zhǔn)形成上下對稱結(jié)構(gòu).在徑向擋板附近有明顯的出流速度(ux>0)分布;在水平對稱面兩側(cè)形成回流(ux<0).

    圖7 法向速度分布(Re=7 000)

    在圖7中可看出明顯的回流范圍,盡管回流速度相對較低,卻增加了阻力和擾動,對后續(xù)流動有較大影響.同時,隨著K值減小,回流區(qū)向徑向擋板移動,對分層器內(nèi)中心區(qū)域的擾動也增大.

    圖8對Re=3 000~11 000范圍內(nèi)不同形式擋板出口截面上回流面積比S(回流面積與出口截面積比)進(jìn)行比較,K=-0.1的回流區(qū)域面積高出平板23 %,比K=0.1高出30 %.

    圖8 回流面積比例與雷諾數(shù)的關(guān)系

    3.3 渦量

    圖9為不同形式擋板出口截面上的二次流流線.

    圖9 二次流流線(Re=7 000)

    比照圖5、6、7可以看出:在水平對稱面兩側(cè)形成對稱漩渦結(jié)構(gòu),漩渦的中心處示蹤劑濃度(t=4 s)最低,局域靜壓力最低,并產(chǎn)生回流.

    為評價出口截面上二次流漩渦的強(qiáng)度,對該截面的法向渦量進(jìn)行計(jì)算,如圖10所示.渦量定義為

    (3)

    由圖10可以看出:忽略壁面影響,高渦量集中在靠近分層器壁面和上下軸向擋板區(qū)域,水平對稱面附近渦量較低.高渦量區(qū)域隨擋板曲率的減小而收縮,|ωx|≤2 s-1的區(qū)域相應(yīng)擴(kuò)大.K=0.1,0,-0.1三種擋板|ωx|≥10 s-1的面積分別占出口截面的30.0 %、25.3 %、18.3 %.

    圖10 渦量分布(Re=7 000)

    圖11為Re=3 000~11 000范圍內(nèi)三種形式出口截面渦量絕對值的均值比較.隨著雷諾數(shù)的增大,不同形式擋板渦量絕對值逐漸升高.平擋板出口截面的渦量最高,凹柱面擋板略低于平擋板,凸柱面擋板遠(yuǎn)低于其他兩種情況.K=0.1比K=0出口截面的渦量僅低2.3 %,而K=-0.1比K=0低19.8 %.

    圖11 絕對渦量平均值

    3.4 湍流強(qiáng)度

    擋板出口截面(分層器入口)流體的湍流強(qiáng)度是設(shè)備內(nèi)流體擾動的前提條件.圖12對Re=7 000條件下三種形式擋板出口截面湍流強(qiáng)度分布進(jìn)行描述.湍流強(qiáng)度定義為

    (4)

    由圖12可以看出:三種形式擋板出口截面上,徑向擋板近壁面湍流強(qiáng)度較高,出口截面中心靠近分層器壁面處偏低.凹柱面相對其他兩種擋板湍流強(qiáng)度分布更均勻,凸柱面最高和最低湍流強(qiáng)度差異最大.K=0.1,0,-0.1三種擋板I≥15 %的面積分別占出口截面的12.0 %、14.0 %、19.3 %.

    圖12 湍流強(qiáng)度分布

    圖13為出口截面上湍流強(qiáng)度均值計(jì)算結(jié)果.擋板形式對出口截面湍流強(qiáng)度均值影響較低,凹柱面板略高于其他形式擋板.雷諾數(shù)對湍流強(qiáng)度均值影響較大,Re=3 000~11 000范圍內(nèi)三種形式出口截面湍流強(qiáng)度隨雷諾數(shù)線性增長.

    圖13 湍流強(qiáng)度均值

    4 結(jié) 論

    在立式圓筒體分層器入口處設(shè)置凹柱面、平板和凸柱面三種擋板,在Re=7 000、分別采用三種結(jié)構(gòu)擋板的條件下,對示蹤劑分布隨時間變化進(jìn)行實(shí)驗(yàn)與數(shù)值模擬,并進(jìn)行比對分析.對分層器入口Re=3 000~11 000范圍內(nèi)的流場進(jìn)行數(shù)值模擬,研究擋板結(jié)構(gòu)對局域流場和出口截面流體力學(xué)性能的影響.具體結(jié)論如下:

    (1) 示蹤劑在擋板內(nèi)流體以水平對稱面為基準(zhǔn)形成上下兩個對稱的旋流結(jié)構(gòu),擋板出口截面水平對稱面兩側(cè)形成對稱漩渦和對稱的回流區(qū)域.

    (2) 擋板出口截面漩渦中心的示蹤劑濃度、靜壓力系數(shù)最低,并存在較低速度的回流;出口截面上水平對稱面附近渦量較低,湍流強(qiáng)度較低.

    (3) 凸柱面擋板出口截面上的回流面積比其他兩種擋板高出很多,絕對渦量平均值卻低很多,三種擋板出口截面湍流強(qiáng)度差異不大.

    (4) 凹柱面擋板使出口截面上的回流面積最小,對湍流強(qiáng)度也起到了限制作用,相對平擋板渦量也略有下降.而凸柱面盡管導(dǎo)致渦量大幅度下降,但回流面積卻大幅度上升.

    符號說明:

    CP—壓力系數(shù);

    I—湍流強(qiáng)度;

    K—相對曲率;

    P—靜壓力,Pa;

    Pin—入口靜壓力,Pa;

    t—時間,s;

    Vc—網(wǎng)格體積,m3;

    ρ—密度,kg/m3;

    uin—入口管內(nèi)平均流速,m·s-1;

    ux—出口截面法向速度,m·s-1;

    uy—出口截面y軸速度,m·s-1;

    uz—出口截面z軸速度,m·s-1;

    ωx—出口截面法向渦量,s-1.

    猜你喜歡
    渦量柱面示蹤劑
    含沙空化對軸流泵內(nèi)渦量分布的影響
    南海東部深水油田水平井產(chǎn)出剖面 示蹤劑監(jiān)測技術(shù)及應(yīng)用
    基于單攝像頭的柱面拼接
    Maple動畫功能在高等數(shù)學(xué)教學(xué)中的應(yīng)用示例(Ⅱ)
    自由表面渦流動現(xiàn)象的數(shù)值模擬
    矩形孔徑柱面鏡面形擬合基底多項(xiàng)式研究
    井間示蹤劑監(jiān)測在復(fù)雜斷塊油藏描述中的應(yīng)用
    錄井工程(2017年1期)2017-07-31 17:44:42
    多示蹤劑成像技術(shù)在腫瘤診斷方面的應(yīng)用研究
    航態(tài)對大型船舶甲板氣流場的影響
    基于節(jié)點(diǎn)剛度的柱面巨型網(wǎng)格結(jié)構(gòu)靜力性能研究
    中文字幕人妻丝袜一区二区| 久久婷婷人人爽人人干人人爱| 观看免费一级毛片| 麻豆国产av国片精品| 亚洲中文av在线| 美女大奶头视频| 男女床上黄色一级片免费看| 9191精品国产免费久久| 亚洲av成人一区二区三| 又黄又粗又硬又大视频| 一级a爱视频在线免费观看| 免费在线观看影片大全网站| av欧美777| 露出奶头的视频| 久久中文字幕人妻熟女| 亚洲五月婷婷丁香| 人人妻人人澡人人看| 99热这里只有精品一区 | 亚洲中文字幕日韩| 男人舔女人的私密视频| 国产一区在线观看成人免费| 久久精品国产清高在天天线| 免费看美女性在线毛片视频| 中文字幕另类日韩欧美亚洲嫩草| 性色av乱码一区二区三区2| av超薄肉色丝袜交足视频| 精品熟女少妇八av免费久了| 欧美色欧美亚洲另类二区| 99久久精品国产亚洲精品| 熟女少妇亚洲综合色aaa.| 最新美女视频免费是黄的| 我的亚洲天堂| 9191精品国产免费久久| 美女国产高潮福利片在线看| 国产午夜福利久久久久久| 国产男靠女视频免费网站| 婷婷六月久久综合丁香| 亚洲精品粉嫩美女一区| 精品卡一卡二卡四卡免费| 色播在线永久视频| 一本大道久久a久久精品| 婷婷精品国产亚洲av在线| 国产三级黄色录像| x7x7x7水蜜桃| 欧美性长视频在线观看| 一区二区日韩欧美中文字幕| 午夜免费激情av| 中文字幕最新亚洲高清| 999精品在线视频| 欧美激情 高清一区二区三区| 男女之事视频高清在线观看| www.www免费av| 欧美黑人巨大hd| 99国产综合亚洲精品| 99国产精品一区二区蜜桃av| 国产色视频综合| 中文字幕精品亚洲无线码一区 | 啦啦啦韩国在线观看视频| 中文字幕精品免费在线观看视频| 亚洲色图av天堂| 一区二区三区激情视频| 成人三级黄色视频| 亚洲熟妇熟女久久| 波多野结衣高清作品| 国产一区二区亚洲精品在线观看| 欧美一区二区国产精品久久精品| 国产一区二区亚洲精品在线观看| 日日干狠狠操夜夜爽| 日韩高清综合在线| 国产精品一二三区在线看| 久久久久久久久久久丰满| 午夜视频国产福利| 婷婷精品国产亚洲av| 22中文网久久字幕| 精华霜和精华液先用哪个| 久久人人爽人人爽人人片va| 卡戴珊不雅视频在线播放| 少妇人妻精品综合一区二区 | 亚洲欧美清纯卡通| 国产v大片淫在线免费观看| 亚洲av中文字字幕乱码综合| 白带黄色成豆腐渣| 欧美bdsm另类| a级毛片免费高清观看在线播放| 在线观看午夜福利视频| 又爽又黄无遮挡网站| 麻豆一二三区av精品| 亚洲电影在线观看av| 精品国产三级普通话版| 嫩草影院新地址| 亚洲精品国产av成人精品 | 色在线成人网| or卡值多少钱| 亚洲国产精品成人久久小说 | 亚洲在线观看片| 欧美日韩国产亚洲二区| 国产精品久久久久久久电影| 国产爱豆传媒在线观看| av天堂在线播放| 男女做爰动态图高潮gif福利片| 九九热线精品视视频播放| 97热精品久久久久久| 久久久久性生活片| 十八禁国产超污无遮挡网站| 国产高清三级在线| 国产人妻一区二区三区在| 免费黄网站久久成人精品| 久久久久免费精品人妻一区二区| 不卡视频在线观看欧美| 嫩草影视91久久| 午夜久久久久精精品| 日韩欧美精品免费久久| 国产av一区在线观看免费| 人人妻人人看人人澡| 久久久久久九九精品二区国产| 国产蜜桃级精品一区二区三区| 欧美在线一区亚洲| 极品教师在线视频| 亚洲性久久影院| 男人舔奶头视频| 波多野结衣高清作品| 午夜福利视频1000在线观看| 欧美国产日韩亚洲一区| 日韩精品中文字幕看吧| 亚洲乱码一区二区免费版| 色哟哟·www| 可以在线观看的亚洲视频| 亚洲三级黄色毛片| 大香蕉久久网| 久久精品国产自在天天线| 国产亚洲91精品色在线| 在线观看午夜福利视频| 亚洲图色成人| 国内久久婷婷六月综合欲色啪| 欧美最黄视频在线播放免费| 免费av不卡在线播放| 久久婷婷人人爽人人干人人爱| 91久久精品国产一区二区成人| 国产视频内射| 可以在线观看毛片的网站| 日韩在线高清观看一区二区三区| 免费在线观看影片大全网站| 国产伦精品一区二区三区视频9| 最近视频中文字幕2019在线8| 国产精品一区www在线观看| 精品99又大又爽又粗少妇毛片| 老女人水多毛片| 久久鲁丝午夜福利片| 亚洲自拍偷在线| 一级毛片电影观看 | 国产又黄又爽又无遮挡在线| 伦理电影大哥的女人| 看非洲黑人一级黄片| 亚洲欧美日韩卡通动漫| 极品教师在线视频| 久久精品国产亚洲av涩爱 | 国产熟女欧美一区二区| 日产精品乱码卡一卡2卡三| 长腿黑丝高跟| 国产亚洲av嫩草精品影院| 亚洲国产精品久久男人天堂| 伊人久久精品亚洲午夜| av中文乱码字幕在线| 麻豆国产97在线/欧美| 欧美在线一区亚洲| 欧美中文日本在线观看视频| 99久久精品一区二区三区| 亚洲一区二区三区色噜噜| 免费电影在线观看免费观看| 久久韩国三级中文字幕| 人妻久久中文字幕网| 日韩国内少妇激情av| 免费电影在线观看免费观看| 国产精品久久久久久久久免| 久久久久久九九精品二区国产| 搡老妇女老女人老熟妇| 丝袜喷水一区| 男人舔女人下体高潮全视频| 免费在线观看成人毛片| 亚洲中文字幕一区二区三区有码在线看| 舔av片在线| 亚洲第一电影网av| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久国产蜜桃| 亚洲天堂国产精品一区在线| 免费看a级黄色片| 国模一区二区三区四区视频| 日韩欧美在线乱码| 国产成年人精品一区二区| 成年免费大片在线观看| 国产高潮美女av| 日本a在线网址| 99热这里只有是精品50| 成人国产麻豆网| 久99久视频精品免费| 大香蕉久久网| 亚洲国产欧美人成| 亚洲不卡免费看| 日韩av不卡免费在线播放| 亚洲欧美日韩卡通动漫| av专区在线播放| 久久精品国产自在天天线| 久久欧美精品欧美久久欧美| 国产一区二区激情短视频| 国产亚洲欧美98| 成人综合一区亚洲| 久久草成人影院| 国产精品久久电影中文字幕| 欧美+日韩+精品| 日本五十路高清| 97在线视频观看| 亚洲人与动物交配视频| 联通29元200g的流量卡| 成人av一区二区三区在线看| 九色成人免费人妻av| 国产淫片久久久久久久久| 国产探花在线观看一区二区| 人妻夜夜爽99麻豆av| 变态另类丝袜制服| 亚洲欧美日韩高清专用| 免费人成在线观看视频色| 一区二区三区高清视频在线| 变态另类丝袜制服| 乱系列少妇在线播放| 日韩欧美在线乱码| 国产一区二区三区av在线 | 国产私拍福利视频在线观看| 51国产日韩欧美| 亚洲色图av天堂| 久久久久久大精品| 大又大粗又爽又黄少妇毛片口| 亚洲成人中文字幕在线播放| 欧美色视频一区免费| 天堂网av新在线| 国产探花极品一区二区| 国产黄色视频一区二区在线观看 | 人妻夜夜爽99麻豆av| 99热网站在线观看| 欧美成人a在线观看| 亚洲精品456在线播放app| 欧美激情久久久久久爽电影| 欧美日韩精品成人综合77777| 18禁在线播放成人免费| 女人被狂操c到高潮| 人妻制服诱惑在线中文字幕| 最新中文字幕久久久久| videossex国产| 日韩大尺度精品在线看网址| 一进一出抽搐gif免费好疼| 综合色av麻豆| 女人被狂操c到高潮| 黄色配什么色好看| 久久久久久伊人网av| 在现免费观看毛片| 免费一级毛片在线播放高清视频| 日日啪夜夜撸| 91在线精品国自产拍蜜月| 亚洲精品乱码久久久v下载方式| 成人亚洲欧美一区二区av| 美女黄网站色视频| 国产麻豆成人av免费视频| 国产伦在线观看视频一区| 日韩精品有码人妻一区| 性色avwww在线观看| 精品少妇黑人巨大在线播放 | 欧美日韩精品成人综合77777| 亚洲四区av| 免费人成在线观看视频色| 亚洲av中文av极速乱| 长腿黑丝高跟| 老司机影院成人| 老师上课跳d突然被开到最大视频| 国产极品精品免费视频能看的| 亚洲精品粉嫩美女一区| 欧美在线一区亚洲| 免费看av在线观看网站| 久久久久久久久中文| 熟女电影av网| 少妇的逼好多水| 最近2019中文字幕mv第一页| av在线亚洲专区| 国产不卡一卡二| 久久国内精品自在自线图片| 99热精品在线国产| 亚洲真实伦在线观看| 女人十人毛片免费观看3o分钟| 级片在线观看| 在线a可以看的网站| 久久久久久国产a免费观看| 国产69精品久久久久777片| 国产精品久久久久久久电影| 天堂动漫精品| 夜夜夜夜夜久久久久| 免费av不卡在线播放| 一级毛片电影观看 | 亚洲性久久影院| avwww免费| 日韩亚洲欧美综合| 搡老熟女国产l中国老女人| 日本 av在线| 国产精品一区www在线观看| 国产伦一二天堂av在线观看| www.色视频.com| 成人性生交大片免费视频hd| 国产精品亚洲美女久久久| 久久99热这里只有精品18| 99在线视频只有这里精品首页| 特级一级黄色大片| 亚洲综合色惰| 久99久视频精品免费| 久久久久国产精品人妻aⅴ院| 国产精品一区二区三区四区久久| 国产在线精品亚洲第一网站| 久久中文看片网| 直男gayav资源| 男人舔女人下体高潮全视频| 精品99又大又爽又粗少妇毛片| 国产一区二区三区av在线 | 亚洲,欧美,日韩| 狂野欧美激情性xxxx在线观看| 99热精品在线国产| 亚洲内射少妇av| 成人美女网站在线观看视频| 精品乱码久久久久久99久播| 国产成人a区在线观看| 99热网站在线观看| 最近最新中文字幕大全电影3| 久久婷婷人人爽人人干人人爱| 亚洲av电影不卡..在线观看| 午夜a级毛片| 亚洲欧美精品综合久久99| 亚洲乱码一区二区免费版| 久久人人爽人人片av| 大型黄色视频在线免费观看| 最近的中文字幕免费完整| 一进一出抽搐gif免费好疼| 久久午夜亚洲精品久久| 人妻夜夜爽99麻豆av| 亚洲av免费高清在线观看| 成人特级黄色片久久久久久久| 国模一区二区三区四区视频| 全区人妻精品视频| 日本a在线网址| 国产探花在线观看一区二区| 小蜜桃在线观看免费完整版高清| 亚洲色图av天堂| 国产av不卡久久| 99热这里只有是精品在线观看| 国语自产精品视频在线第100页| 变态另类丝袜制服| 色av中文字幕| 国产爱豆传媒在线观看| av视频在线观看入口| 国产伦精品一区二区三区四那| 最新中文字幕久久久久| 成人性生交大片免费视频hd| 久久久久国产精品人妻aⅴ院| 亚洲18禁久久av| 亚洲欧美中文字幕日韩二区| 国语自产精品视频在线第100页| 九色成人免费人妻av| 亚洲av不卡在线观看| 两个人的视频大全免费| 听说在线观看完整版免费高清| 别揉我奶头 嗯啊视频| 欧美+亚洲+日韩+国产| 亚洲五月天丁香| 午夜免费男女啪啪视频观看 | 少妇人妻精品综合一区二区 | 亚洲精华国产精华液的使用体验 | 男女边吃奶边做爰视频| 别揉我奶头 嗯啊视频| 婷婷精品国产亚洲av在线| 男人的好看免费观看在线视频| 亚洲,欧美,日韩| 国产黄a三级三级三级人| 99久久精品热视频| 久久久久性生活片| 身体一侧抽搐| 日韩av在线大香蕉| 亚洲国产日韩欧美精品在线观看| 日本 av在线| 麻豆国产av国片精品| 色5月婷婷丁香| 亚洲av.av天堂| 久久精品国产亚洲网站| 插逼视频在线观看| 午夜亚洲福利在线播放| 青春草视频在线免费观看| 久久精品久久久久久噜噜老黄 | 草草在线视频免费看| 亚洲成a人片在线一区二区| 国产一级毛片七仙女欲春2| 女人十人毛片免费观看3o分钟| 午夜视频国产福利| 久久中文看片网| 91精品国产九色| 91午夜精品亚洲一区二区三区| 给我免费播放毛片高清在线观看| 99久久九九国产精品国产免费| 欧美性感艳星| 国产av麻豆久久久久久久| 国产精品免费一区二区三区在线| 久久精品国产99精品国产亚洲性色| 日韩av不卡免费在线播放| 国产亚洲av嫩草精品影院| 国内精品美女久久久久久| 国产人妻一区二区三区在| 国产三级中文精品| 又爽又黄无遮挡网站| 色噜噜av男人的天堂激情| 高清毛片免费看| 国产私拍福利视频在线观看| 亚洲在线自拍视频| 亚洲精品国产成人久久av| 国产伦一二天堂av在线观看| 99久久无色码亚洲精品果冻| 九九在线视频观看精品| 久99久视频精品免费| av中文乱码字幕在线| 久久精品国产亚洲网站| 神马国产精品三级电影在线观看| 99久国产av精品| 99久久精品热视频| 黄色欧美视频在线观看| or卡值多少钱| 少妇高潮的动态图| 99久国产av精品国产电影| 天天躁夜夜躁狠狠久久av| 亚洲人与动物交配视频| 人妻夜夜爽99麻豆av| 午夜免费激情av| 国产午夜精品论理片| 日韩国内少妇激情av| 波多野结衣高清作品| 日本三级黄在线观看| 97人妻精品一区二区三区麻豆| 免费观看的影片在线观看| 国产白丝娇喘喷水9色精品| 卡戴珊不雅视频在线播放| 大香蕉久久网| 一级毛片久久久久久久久女| 亚洲成av人片在线播放无| 亚洲自偷自拍三级| 又黄又爽又刺激的免费视频.| 亚洲欧美成人综合另类久久久 | 亚洲精品粉嫩美女一区| 又黄又爽又刺激的免费视频.| 白带黄色成豆腐渣| 国产伦精品一区二区三区视频9| 午夜精品国产一区二区电影 | 亚洲av成人精品一区久久| 国产高清激情床上av| 老司机福利观看| 国产真实伦视频高清在线观看| 久久人人爽人人片av| 天天躁夜夜躁狠狠久久av| 99在线视频只有这里精品首页| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 少妇高潮的动态图| 日日干狠狠操夜夜爽| 精品少妇黑人巨大在线播放 | 亚洲激情五月婷婷啪啪| 国产人妻一区二区三区在| 99久久久亚洲精品蜜臀av| 亚洲国产色片| 老熟妇仑乱视频hdxx| 可以在线观看的亚洲视频| 久久久久久九九精品二区国产| 深爱激情五月婷婷| 1024手机看黄色片| h日本视频在线播放| 成人午夜高清在线视频| 亚洲性久久影院| 欧美性猛交╳xxx乱大交人| 亚洲精华国产精华液的使用体验 | 欧美精品国产亚洲| 黄色欧美视频在线观看| 国产精品人妻久久久影院| 最近最新中文字幕大全电影3| 国产成人91sexporn| 别揉我奶头 嗯啊视频| 久久精品夜夜夜夜夜久久蜜豆| 成人亚洲精品av一区二区| 欧美日韩精品成人综合77777| 国产v大片淫在线免费观看| 亚洲欧美日韩无卡精品| 久久精品国产99精品国产亚洲性色| 国产精品无大码| 欧美3d第一页| 久久国产乱子免费精品| 成人二区视频| 日本色播在线视频| 晚上一个人看的免费电影| 久久精品国产自在天天线| 亚洲久久久久久中文字幕| 日本 av在线| 直男gayav资源| 免费av毛片视频| 久久人妻av系列| 亚洲熟妇中文字幕五十中出| 成人综合一区亚洲| 精品久久久久久久久av| 99久久精品一区二区三区| 久久久久久久午夜电影| 亚洲人成网站在线观看播放| 最近视频中文字幕2019在线8| 国产亚洲精品av在线| 中国国产av一级| 成人鲁丝片一二三区免费| 看片在线看免费视频| 国产精品,欧美在线| 国产又黄又爽又无遮挡在线| 亚洲一级一片aⅴ在线观看| 精品久久久久久久久久久久久| 亚洲综合色惰| av在线老鸭窝| 长腿黑丝高跟| 国产伦精品一区二区三区四那| 欧美一区二区国产精品久久精品| 亚洲人成网站在线观看播放| 综合色av麻豆| 97超碰精品成人国产| 欧美性猛交╳xxx乱大交人| 亚洲,欧美,日韩| 日日干狠狠操夜夜爽| 久久久久国产精品人妻aⅴ院| 免费人成视频x8x8入口观看| 久久鲁丝午夜福利片| 卡戴珊不雅视频在线播放| 免费av不卡在线播放| 春色校园在线视频观看| 不卡一级毛片| 亚洲成人久久爱视频| 99国产极品粉嫩在线观看| 成年女人毛片免费观看观看9| 3wmmmm亚洲av在线观看| 直男gayav资源| 日韩三级伦理在线观看| 国产伦精品一区二区三区四那| 欧美日韩乱码在线| 日本成人三级电影网站| 老司机午夜福利在线观看视频| 日本色播在线视频| 亚洲久久久久久中文字幕| 久久99热这里只有精品18| 久久久a久久爽久久v久久| 波多野结衣高清作品| 真人做人爱边吃奶动态| 国产精品精品国产色婷婷| 国产精品99久久久久久久久| 亚洲aⅴ乱码一区二区在线播放| 欧美区成人在线视频| 久久久国产成人免费| 此物有八面人人有两片| 成人二区视频| 久久人人爽人人爽人人片va| 白带黄色成豆腐渣| 秋霞在线观看毛片| 亚洲av美国av| 黄色一级大片看看| 少妇人妻精品综合一区二区 | 国产精品av视频在线免费观看| 久久久久久久久中文| 午夜老司机福利剧场| 久久人人精品亚洲av| 精品久久久噜噜| 69av精品久久久久久| 亚洲激情五月婷婷啪啪| 国产大屁股一区二区在线视频| 免费观看人在逋| 日韩欧美 国产精品| 别揉我奶头 嗯啊视频| 成人漫画全彩无遮挡| 国产一级毛片七仙女欲春2| 午夜精品国产一区二区电影 | 亚洲婷婷狠狠爱综合网| 亚洲精品一区av在线观看| 精品一区二区三区av网在线观看| 99国产精品一区二区蜜桃av| 久久欧美精品欧美久久欧美| 十八禁国产超污无遮挡网站| av视频在线观看入口| 22中文网久久字幕| 亚洲天堂国产精品一区在线| 精品免费久久久久久久清纯| 在线观看av片永久免费下载| 乱码一卡2卡4卡精品| а√天堂www在线а√下载| 大型黄色视频在线免费观看| 国产大屁股一区二区在线视频| 嫩草影院新地址| 男女啪啪激烈高潮av片| 国产aⅴ精品一区二区三区波| 狠狠狠狠99中文字幕| 亚洲av中文av极速乱| 亚洲在线自拍视频| 99久久无色码亚洲精品果冻| 2021天堂中文幕一二区在线观| 青春草视频在线免费观看| videossex国产| 97人妻精品一区二区三区麻豆| 看片在线看免费视频| 日本一二三区视频观看| 在线观看一区二区三区| 日日摸夜夜添夜夜添av毛片| 国产淫片久久久久久久久| 国产成人91sexporn| 亚洲最大成人中文| 韩国av在线不卡| 欧美精品国产亚洲| 小说图片视频综合网站|