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

    通風(fēng)及內(nèi)熱源參數(shù)的方腔內(nèi)混合對流模擬

    2021-03-16 09:28:44孫夢楠宋桂秋周世華董祉序
    關(guān)鍵詞:流口賽爾等溫線

    孫夢楠, 宋桂秋, 周世華, 董祉序

    (1.東北大學(xué) 機(jī)械工程與自動化學(xué)院,遼寧 沈陽 110819;2.沈陽工業(yè)大學(xué) 機(jī)械工程學(xué)院,遼寧 沈陽 110870)

    在針對換熱器[1]、空調(diào)系統(tǒng)[2]、電子設(shè)備冷卻[3]等眾多工程設(shè)計和能源問題的研究中,被分析流體往往受到外力和溫度的共同影響,多以混合對流形式存在。合理設(shè)計和布置設(shè)備結(jié)構(gòu)不僅可以有效加大能源利用率,還能夠減少重要零件的損耗。利用方腔模型可以分析多種復(fù)雜場內(nèi)流體的運(yùn)動特性和傳熱規(guī)律[4]。在針對混合對流的研究中,也多采用具有內(nèi)部加熱元件的開口方腔模型進(jìn)行模擬。Radhakrishnan等[5]通過實驗證明了模型的可靠性,并指出采用數(shù)值模擬方法有利于優(yōu)化熱流表現(xiàn)。孫猛等[6]在不同的格拉曉夫數(shù)、瑞利數(shù)以及縱橫比的條件下,研究了外界來流對方腔內(nèi)混合對流換熱的影響。Karimi等[7-8]建立了具有多個加熱單元的復(fù)雜混合對流模型,并對熱源的大小和位置進(jìn)行了討論。Chamkha等[9-10]借助于方腔模型對多種結(jié)構(gòu)因素進(jìn)行分析,研究結(jié)果表明,進(jìn)、出流口的位置對混合對流的影響很大,在改善換熱效率方面起著關(guān)鍵作用。

    隨著流體流動和換熱機(jī)理研究的深入,新的數(shù)值計算方法也在不斷涌現(xiàn)。格子Boltzmann方法(LBM)作為介觀模擬算法的代表,在近幾十年里發(fā)展迅速。與傳統(tǒng)計算方法相比,LBM使用簡化的動力學(xué)模型來處理微觀過程,這使其在具有出色并行計算能力的同時能夠更直觀地表現(xiàn)粒子間的相互作用[11]。多參數(shù)弛豫時間(MRT)的引入也進(jìn)一步提高了LBM的計算精度和穩(wěn)定性[12]。目前,LBM主要通過多速度法[13]、多分布法[14]和混合法[15]對多場問題進(jìn)行模擬,并在混合對流、輻射傳熱、雙擴(kuò)散對流等問題上得到了驗證。

    綜上可以發(fā)現(xiàn),在分析腔內(nèi)熱源對混合對流的影響時,研究工作主要集中在熱元件的布置方式和尺寸上,而對熱元件的形狀特征很少提及。此外,內(nèi)熱源結(jié)構(gòu)和外部方腔的開口位置共同變化下的流體和換熱表現(xiàn)也很少被討論。因此,本文采用結(jié)合MRT-LBM與有限差分法(FDM)的耦合算法[15],通過模擬流函數(shù)線圖、等溫線圖和努賽爾數(shù),探討了不同進(jìn)、出流口位置下熱源形狀變化對方腔內(nèi)混合對流的流動、溫度及換熱特性的影響。

    1 方腔內(nèi)混合對流模型的建立

    1.1 問題描述

    本文研究的含內(nèi)熱源的開口方腔結(jié)構(gòu)如圖1所示。計算基于笛卡爾坐標(biāo)系(x,y),坐標(biāo)原點(diǎn)位于計算域的左下角。正方形腔體與等長寬中心加熱元件的尺寸分別定義為L和d=0.3L。方腔內(nèi)壁為具有溫度Tc的恒溫壁面,熱源外表面的溫度為Th(Th>Tc)。

    具有速度uin和溫度Tin=Tc的流體通過左腔壁上尺寸為w=0.2L的進(jìn)流口流入,最終從對側(cè)壁面上相同尺寸的出流口以速度uout和溫度Tout流出,流動過程中所受重力加速度為g。以進(jìn)、出流口下邊緣高度l和l′ 來定義開口位置,并分別在底部、中部和頂部開口下進(jìn)行討論。同時,引入圓角半徑R來討論加熱元件的形狀變化,R的值取為0、0.25d和0.5d。由圖1中可見,R=0和R=0.5d分別表示熱源具有正方形截面和圓形截面。

    1.2 計算模型

    本文采用普朗特數(shù)Pr為0.71的空氣作為工作流體,認(rèn)為它是滿足Boussinesq近似的牛頓流體,并以穩(wěn)態(tài)、層流運(yùn)動。因此,可以得到混合對流的無量綱控制方程為:

    (1)

    (2)

    (3)

    (4)

    引入努賽爾數(shù)來表征熱源表面的對流換熱強(qiáng)度。由局部溫度梯度可求得的局部努賽爾數(shù)Nu和平均努賽爾數(shù)Nuav分別為:

    (5)

    (6)

    式中:n為沿?zé)岜砻嫱夥ň€方向的矢量;A為熱源廓線周長。

    用于求解速度場的MRT-LBM采用如圖2所示的9速度模型(D2Q9),其演化方程為:

    f(r+ciδt,t+δt)-f(r,t)=

    -M-1S(m(r,t)-m(eq)(r,t))+F

    (7)

    式中:f(r,t)=(f0,f1, …,f8)T為位置r處的速度分布函數(shù)矢量;m(r,t)、m(eq)(r,t) 為矩和矩對應(yīng)的平衡態(tài);M為正交轉(zhuǎn)換矩陣,并滿足m=Mf;S為馳豫系數(shù)矩陣;δt、δx為單位長度的離散時間步長和格子間距;c=δx/δt為格子速率。

    各個速度方向上的粒子速度集c,權(quán)重系數(shù)集ω和體積力F分別定義為:

    (8)

    (9)

    (10)

    在LBM中,宏觀物理量被視為分布函數(shù)的平均特性,因此流體密度ρ和速度u可以表示為:

    (11)

    (12)

    圖2 D2Q9模型Fig.2 Model of D2Q9

    在MRT-LBM的算式中,方腔內(nèi)壁面和進(jìn)、出流口位置采用非平衡態(tài)外推邊界條件[16]。對于曲線邊界與網(wǎng)格不完全一致的熱源表面,有些粒子在被壁面反彈后會脫離節(jié)點(diǎn),落到間隙位置(圖3中的xg位置),則采用線性插值邊界條件[17]。如圖3所示,該方法首先假設(shè)分布函數(shù)f1(xf,t)可以穿過邊界節(jié)點(diǎn)xw到達(dá)實體節(jié)點(diǎn)xs,由線性插值得到:

    f1(xw,t+δt)=(1-q)f1(xf,t+δt)+qf1(xs,t+δt)

    (13)

    式中q為流體節(jié)點(diǎn)到邊界間的距離與網(wǎng)格長度之比,q=|xf-xw|/|xf-xs|。再根據(jù)動量守恒和固體邊界的無滑移條件,最終得到在實體節(jié)點(diǎn)xs處的分布函數(shù)為:

    (1-q)f1(xff,t)]

    (14)

    式(14)具有二階精度,并可以避免按q≤1/2和q>1/2劃分時所引起的數(shù)值振蕩。

    在計算溫度場時,采用有限差分方法在與MRT-LBM相同的網(wǎng)格上對方程(4)進(jìn)行求解[18]。這里采用一階向前和二階中心差分方法來分別求解時間和空間的偏微分方程,即有:

    (15)

    (16)

    (17)

    圖3 曲面插值邊界條件Fig.3 Interpolation boundary condition for curved wall

    不同位置的無量綱熱邊界條件定義為:在方腔內(nèi)壁面和進(jìn)流口位置,θ=0;在熱源表面,θ=1;在出流口位置,?θ/?X=0。

    1.3 模型驗證

    為確認(rèn)所采用算法的可靠性和準(zhǔn)確性,以文獻(xiàn)[9]中所描述的具有中部進(jìn)流口、頂部出流口以及方形截面內(nèi)熱源的方腔混合對流模型作為算例進(jìn)行驗證。通過對比圖4中給出的流場和溫度場分布(Re=200,Ri=1)可以看出,采用本文方法所得到的模擬結(jié)果與文獻(xiàn)中的結(jié)果具有很好的一致性。同時,表1中的對比表明,當(dāng)Ri變化時,2種計算的Nuav具有相同的趨勢,最大百分比差異為3.52%,說明本文的計算方法是具有足夠精度的。

    圖4 本文和文獻(xiàn)[9]的流場和溫度場的算例驗證Fig.4 Comparison of streamlines and isotherms between the present work and that of ref.[9]

    表1 模擬結(jié)果對比Table 1 Comparison of simulation results

    考慮到算法對網(wǎng)格的敏感性,分別采用100×100、150×150、200×200和250×250的網(wǎng)格分布對本文的混合對流模型進(jìn)行計算。當(dāng)網(wǎng)格超過200×200時,繼續(xù)加密網(wǎng)格對模擬結(jié)果的改進(jìn)很小,因此在下述分析中均采用200×200的網(wǎng)格。

    2 對流換熱的計算結(jié)果與分析

    為保證混合對流換熱中始終存在強(qiáng)制對流和自然對流,且方腔內(nèi)的流體為層流狀態(tài),本文在Re=200,Ri=1的工況下對參數(shù)進(jìn)行討論。不同形狀的加熱元件在不同進(jìn)、出流口位置下的流場、溫度場分布如圖5~7所示。

    圖5給出了當(dāng)進(jìn)流口位于底部(l=0)時腔內(nèi)混合對流的流函數(shù)線圖和等溫線圖。在圖5(a)中,當(dāng)出流口也位于底部(l′=0)時,流體從方腔底部直接流過,慣性力使加熱元件被逆時針環(huán)流包圍。由于強(qiáng)迫對流和自然對流的共同作用,在加熱元件左側(cè)和方腔右上角還會出現(xiàn)2個渦流。R的增大會使環(huán)流增強(qiáng),渦流被擠壓。當(dāng)l′=0.4L時(圖5(b)),冷來流與加熱元件的右下角發(fā)生直接接觸,腔內(nèi)右部的流動發(fā)生變化。由于進(jìn)流口附近的流體受到出流口位置的影響較小,所以熱源左側(cè)與腔壁間仍存在渦流。當(dāng)R=0時,有小部分來流會從左側(cè)繞經(jīng)熱源后再流至腔外,但隨著R的增大,腔內(nèi)逐漸出現(xiàn)與圖5(a)中相似的環(huán)流結(jié)構(gòu)。當(dāng)l′=0.8L時(圖5(c)),熱源被整個包含在來流中,R的增大會減少上部繞流,但也為渦流的流動提供了更大的空間。在出流口位置由l′=0升高到l′=0.8L的過程中,冷來流的上移使腔內(nèi)右部的等溫線向上收縮,熱源底部的等溫線也變得越來越密集,這說明該位置的具有較大的溫差。R的增大會使熱源附近的等溫線發(fā)生變化,但對整體溫度場的影響不大。

    圖5 進(jìn)流口位于底部時,R和出流口位置對流函數(shù)線和等溫線的影響Fig.5 Effects of R and outlet location on the streamlines and isotherms for cavity with bottom inlet

    圖6給出的是當(dāng)進(jìn)流口位于中部(l=0.4L)時腔內(nèi)混合對流的流函數(shù)線圖和等溫線圖。此時,加熱元件正對進(jìn)流口,在進(jìn)行熱交換的同時也阻礙了流體運(yùn)動。流入的流體會在進(jìn)流口兩側(cè)產(chǎn)生上大下小的兩個低溫角渦,并從熱源上下兩側(cè)流過,最終在出流口處匯合。等溫線圖表明,腔體左側(cè)溫度主要由來流決定,而當(dāng)流體流至右側(cè)時熱源的作用開始變得明顯起來。隨著出流口位置l′的升高,右側(cè)腔體內(nèi)的流函數(shù)線和等溫線受到較大影響,但腔內(nèi)流體的整體流動模式變化很小。R的增大會使渦流的運(yùn)動區(qū)域變大,并使對應(yīng)位置的等溫線向右傾斜。

    進(jìn)流口位于頂部(l=0.8L)時腔內(nèi)混合對流的流函數(shù)線和等溫線如圖7所示。在圖7(a)和(b)中(l′=0和l′=0.4L),腔內(nèi)的繞流雖然占主導(dǎo)地位,但它們難以覆蓋進(jìn)流口下方的空間,因而產(chǎn)生了順時針渦流。等溫線聚集在熱源上表面,可以看出上部繞流對溫度場的影響要大于下部繞流。隨著R的增大,流體更多地從腔體上部流過,下部繞流減少,并逐漸脫離熱源表面。當(dāng)進(jìn)流口和出流口都位于圖7(c)頂部(l=l′=0.8L)時,流體不再受加熱元件阻礙,幾乎水平地流過腔體頂部。在冷來流未經(jīng)過的區(qū)域內(nèi),溫度也分布得較為均勻。R的增大會改變腔內(nèi)其他位置的流體運(yùn)動,使部分渦流開始向環(huán)流轉(zhuǎn)變。由于重力的存在,腔體內(nèi)的流函數(shù)線和等溫線分布與圖5(a)并不對稱。

    從上面的討論可以看出,進(jìn)、出流口位置選擇對腔內(nèi)流場的影響要大于加熱元件的形狀變化,R的變化對溫度場的影響并不明顯。實際上,熱源形狀的特征主要反映在加熱表面的對流換熱上。

    圖8給出了9種不同方腔結(jié)構(gòu)下的熱源表面局部努賽爾數(shù)Nu,圖中以熱源上表面與其中心線的交點(diǎn)為橫坐標(biāo)零點(diǎn),沿逆時針方向繪制了Nu曲線,R為0、0.25d和0.5d時所對應(yīng)的熱表面的無量綱長度分別為1.2、1.07和0.94。當(dāng)R=0時,在熱源的4個直角位置會出現(xiàn)尖銳的峰值,這表明該位置的溫度梯度很大,在相應(yīng)的溫度場中也會出現(xiàn)了密集的等溫線。隨著R的增大,4個峰值開始下降并變得平緩,至熱源截面變?yōu)閳A(R=0.5d)時,Nu曲線中只剩下一個明顯的主峰值。更改進(jìn)、出流口的位置可以改變峰值的大小和位置分布,相比而言,進(jìn)流口的影響要大于出流口。

    圖9中給出了不同方腔結(jié)構(gòu)下熱源表面的平均努賽爾數(shù)Nuav,圖中的大寫字母B、M和T分別表示底部、中部和頂部進(jìn)流口;B′、M′和T′分別表示底部、中部和頂部出流口。與局部努賽爾數(shù)的最大值出現(xiàn)在R=0時相反,平均努賽爾數(shù)是隨著R的增大而增大的,這說明雖然一些特殊位置的局部對流換熱能力有所下降,但整體表面的換熱卻是增強(qiáng)的,且Nuav的值并不與換熱表面的面積大小成正比。另外可以看出,當(dāng)方腔的進(jìn)、出流口在同一位置時,冷流體直接從頂部或底部通道流出腔體,Nuav較??;而當(dāng)方腔具有中部進(jìn)流口和頂部出流口時,熱源與來流正面接觸,換熱充分,Nuav取得最大值。

    3 結(jié)論

    1)流場和溫度場分布與熱源形狀關(guān)系不大,而主要受進(jìn)、出流口位置影響。當(dāng)進(jìn)流口位于腔體底部或頂部時,隨著出流口位置與進(jìn)流口位置之間距離的增大,熱源周圍的環(huán)流被擠壓,并逐漸被外部來流所取代;而當(dāng)進(jìn)流口位于中部時,外部來流始終流經(jīng)熱源的上、下表面,出流口位置不改變腔內(nèi)流體的流動模式。在溫度場中,由于受到外部強(qiáng)制對流的影響,等溫線會向進(jìn)流口附近的熱源表面以及出流口位置聚集。

    2)當(dāng)熱源截面為圓形時(R=0.5d),熱源表面的局部努賽爾數(shù)Nu曲線具有平緩變化的單峰值;隨著R的減小,Nu曲線上會逐漸出現(xiàn)多個尖銳峰值。Nu的最大值出現(xiàn)在熱源表面靠近進(jìn)流口位置的一側(cè),而出流口位置對Nu的影響不大。當(dāng)進(jìn)、出流口分別位于頂部和底部時,可在正方形熱源截面(R=0)的左上角獲得最佳局部對流換熱位置。

    3)當(dāng)進(jìn)、出流口位置固定時,熱源表面的平均努賽爾數(shù)Nuav隨著R的增大而增大。當(dāng)R不改變時,Nuav的最大值在方腔具有中部進(jìn)流口和頂部出流口時出現(xiàn)。 因此,在進(jìn)、出流口分別位于中部和頂部時,采用圓形截面的熱源(R=0.5d)可以使熱源整體表面的對流換熱最強(qiáng)。

    猜你喜歡
    流口賽爾等溫線
    萊賽爾產(chǎn)業(yè)力促優(yōu)化供給
    低溫吸附劑的低溫低壓吸附等溫線分布研究
    萊賽爾織物長車工藝探討
    轉(zhuǎn)套式配流系統(tǒng)配流口與泵腔壓力特性對比分析
    頌流口燕巖嶺
    可調(diào)底流口旋流器在某尾礦筑壩中的應(yīng)用
    轉(zhuǎn)套式配流系統(tǒng)配流口結(jié)構(gòu)及對工作脈動的影響
    《賽爾號大電影5:雷神崛起》
    電影故事(2015年29期)2015-02-27 09:02:53
    如何在新課改背景下突破等溫線判讀中的難點(diǎn)
    雷伊大戰(zhàn)孫悟空
    亚洲精品美女久久av网站| 18禁黄网站禁片午夜丰满| 亚洲精品国产区一区二| 久久 成人 亚洲| 国产成人欧美在线观看 | 可以免费在线观看a视频的电影网站| 天堂中文最新版在线下载| 亚洲熟妇熟女久久| www.999成人在线观看| 久久精品熟女亚洲av麻豆精品| 欧美日韩精品网址| av不卡在线播放| 一区二区三区国产精品乱码| 欧美久久黑人一区二区| 国内久久婷婷六月综合欲色啪| 精品无人区乱码1区二区| 国产不卡一卡二| 一区二区三区国产精品乱码| 可以免费在线观看a视频的电影网站| 精品国产一区二区三区久久久樱花| 一级黄色大片毛片| 91在线观看av| 黄片大片在线免费观看| 国产欧美亚洲国产| 精品人妻在线不人妻| 精品久久久久久,| videos熟女内射| 久久久久久人人人人人| 激情视频va一区二区三区| 国产精品乱码一区二三区的特点 | 免费黄频网站在线观看国产| 欧美不卡视频在线免费观看 | 侵犯人妻中文字幕一二三四区| 亚洲欧洲精品一区二区精品久久久| 最新的欧美精品一区二区| 久久久水蜜桃国产精品网| 国产又爽黄色视频| 自线自在国产av| 9色porny在线观看| 欧美精品一区二区免费开放| 一个人免费在线观看的高清视频| avwww免费| 欧美 日韩 精品 国产| 成人精品一区二区免费| 国产精品欧美亚洲77777| 日韩一卡2卡3卡4卡2021年| 精品福利永久在线观看| 最近最新中文字幕大全免费视频| 天堂√8在线中文| 久久久国产欧美日韩av| 精品国产亚洲在线| 欧美人与性动交α欧美软件| av天堂久久9| 女人被狂操c到高潮| 中文欧美无线码| 午夜福利欧美成人| 午夜视频精品福利| 国产精品99久久99久久久不卡| 久久精品国产综合久久久| 午夜免费鲁丝| 美女扒开内裤让男人捅视频| 日韩有码中文字幕| 一级片免费观看大全| 国产精品av久久久久免费| 一二三四社区在线视频社区8| 亚洲男人天堂网一区| 老汉色av国产亚洲站长工具| 欧美一区二区国产精品久久精品| 有码 亚洲区| 久久久久久久午夜电影| 日韩欧美三级三区| 久久人人精品亚洲av| 午夜福利在线观看免费完整高清在 | 久久国产精品人妻蜜桃| 久久草成人影院| 又爽又黄无遮挡网站| 国产精品亚洲av一区麻豆| 国产精品亚洲av一区麻豆| 蜜桃亚洲精品一区二区三区| 成熟少妇高潮喷水视频| 国产单亲对白刺激| 精品午夜福利视频在线观看一区| 欧美另类亚洲清纯唯美| 最近最新中文字幕大全电影3| 日日干狠狠操夜夜爽| eeuss影院久久| 午夜福利高清视频| 日本一本二区三区精品| 好男人在线观看高清免费视频| 国产成人福利小说| 亚洲av美国av| 宅男免费午夜| 国产精品一及| 日韩欧美一区二区三区在线观看| 男女下面进入的视频免费午夜| 日韩欧美精品v在线| 亚洲乱码一区二区免费版| 欧美+亚洲+日韩+国产| 亚洲人成网站在线播放欧美日韩| 日本熟妇午夜| 国产精品电影一区二区三区| or卡值多少钱| 在线观看免费视频日本深夜| 成年免费大片在线观看| 国产成人影院久久av| 69av精品久久久久久| 久久久精品欧美日韩精品| 久久久久久国产a免费观看| 日本精品一区二区三区蜜桃| 国产精品,欧美在线| www日本黄色视频网| 亚洲精品一区av在线观看| 久久欧美精品欧美久久欧美| 黄色日韩在线| 激情在线观看视频在线高清| 欧美日本视频| 日韩大尺度精品在线看网址| 黄色丝袜av网址大全| 蜜桃亚洲精品一区二区三区| 精品无人区乱码1区二区| 亚洲人与动物交配视频| 一本一本综合久久| 久久久久九九精品影院| 欧美成人a在线观看| 亚洲欧美日韩卡通动漫| 香蕉丝袜av| 淫秽高清视频在线观看| 亚洲国产欧美人成| 床上黄色一级片| 波多野结衣高清作品| 久久99热这里只有精品18| av天堂在线播放| 大型黄色视频在线免费观看| 欧美区成人在线视频| 欧美+亚洲+日韩+国产| 国产av在哪里看| 国产男靠女视频免费网站| 国产一区二区在线观看日韩 | 亚洲欧美日韩无卡精品| 精品国产美女av久久久久小说| 欧美黑人欧美精品刺激| a级毛片a级免费在线| 国模一区二区三区四区视频| 母亲3免费完整高清在线观看| 免费看美女性在线毛片视频| 亚洲av日韩精品久久久久久密| 综合色av麻豆| 免费看美女性在线毛片视频| 男女视频在线观看网站免费| а√天堂www在线а√下载| 一个人免费在线观看电影| 成人无遮挡网站| av女优亚洲男人天堂| 又爽又黄无遮挡网站| 一本精品99久久精品77| 麻豆一二三区av精品| 国产精品99久久99久久久不卡| 亚洲一区二区三区色噜噜| 人人妻,人人澡人人爽秒播| www日本黄色视频网| 亚洲av日韩精品久久久久久密| 91av网一区二区| 成年版毛片免费区| 日本熟妇午夜| 淫妇啪啪啪对白视频| 亚洲人成网站在线播| 成年免费大片在线观看| 国产精品永久免费网站| 国产又黄又爽又无遮挡在线| 亚洲不卡免费看| 看免费av毛片| 中文资源天堂在线| 亚洲精品一卡2卡三卡4卡5卡| 91字幕亚洲| 3wmmmm亚洲av在线观看| 国产亚洲欧美在线一区二区| 国产精品 国内视频| 毛片女人毛片| av片东京热男人的天堂| 国产精品久久久久久久久免 | 一个人看视频在线观看www免费 | 日韩欧美一区二区三区在线观看| 久久久久久久久大av| 99精品欧美一区二区三区四区| 变态另类丝袜制服| 一个人免费在线观看的高清视频| 一个人免费在线观看电影| 日韩高清综合在线| 成人av在线播放网站| 国产美女午夜福利| 国产精品久久久久久久久免 | a在线观看视频网站| 久久久国产成人精品二区| 91麻豆av在线| 一级黄片播放器| 性欧美人与动物交配| 一卡2卡三卡四卡精品乱码亚洲| a级一级毛片免费在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国内精品美女久久久久久| 最新美女视频免费是黄的| 欧美中文综合在线视频| 中文在线观看免费www的网站| 午夜视频国产福利| 在线观看66精品国产| 男人舔奶头视频| 国产真实伦视频高清在线观看 | 97超级碰碰碰精品色视频在线观看| 色尼玛亚洲综合影院| 亚洲av中文字字幕乱码综合| 午夜视频国产福利| 国产亚洲精品一区二区www| 天天一区二区日本电影三级| 久久精品国产亚洲av涩爱 | 欧美性猛交黑人性爽| 露出奶头的视频| 99国产精品一区二区蜜桃av| 一本久久中文字幕| 一进一出抽搐动态| 精品午夜福利视频在线观看一区| 久久人妻av系列| 午夜日韩欧美国产| 国产精品1区2区在线观看.| 日韩欧美 国产精品| 亚洲欧美日韩无卡精品| 国产私拍福利视频在线观看| 麻豆国产97在线/欧美| 制服丝袜大香蕉在线| 欧美在线一区亚洲| 欧美成狂野欧美在线观看| www国产在线视频色| 久久国产精品人妻蜜桃| 网址你懂的国产日韩在线| 老熟妇仑乱视频hdxx| 国产一区在线观看成人免费| 99热只有精品国产| 亚洲av熟女| 少妇人妻精品综合一区二区 | 久久国产精品影院| 亚洲精品亚洲一区二区| 在线观看一区二区三区| 老汉色av国产亚洲站长工具| 精品免费久久久久久久清纯| 久99久视频精品免费| 国产 一区 欧美 日韩| 成人永久免费在线观看视频| 国产极品精品免费视频能看的| 久久精品夜夜夜夜夜久久蜜豆| 亚洲18禁久久av| 亚洲aⅴ乱码一区二区在线播放| 一级黄片播放器| 人妻丰满熟妇av一区二区三区| 婷婷精品国产亚洲av| 搡老岳熟女国产| a级毛片a级免费在线| 一级作爱视频免费观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩 欧美 亚洲 中文字幕| 国产成人福利小说| 国产高清视频在线观看网站| 91麻豆精品激情在线观看国产| 18禁国产床啪视频网站| 久久久成人免费电影| 亚洲精品粉嫩美女一区| 午夜老司机福利剧场| 狂野欧美激情性xxxx| 精品一区二区三区av网在线观看| 少妇熟女aⅴ在线视频| 神马国产精品三级电影在线观看| 国产成年人精品一区二区| 午夜福利视频1000在线观看| 精品久久久久久成人av| 国产免费男女视频| 女人高潮潮喷娇喘18禁视频| 一个人免费在线观看的高清视频| 国产成人av教育| 亚洲av成人精品一区久久| 老司机在亚洲福利影院| 女人高潮潮喷娇喘18禁视频| 国产免费av片在线观看野外av| 亚洲狠狠婷婷综合久久图片| 嫩草影院精品99| 午夜福利在线观看吧| 亚洲人成电影免费在线| 在线播放国产精品三级| 99久久99久久久精品蜜桃| 精品一区二区三区视频在线观看免费| 成人三级黄色视频| 亚洲精品粉嫩美女一区| 男人和女人高潮做爰伦理| 亚洲久久久久久中文字幕| netflix在线观看网站| 一进一出抽搐gif免费好疼| 午夜a级毛片| 丝袜美腿在线中文| 99热这里只有精品一区| 国产探花极品一区二区| av在线蜜桃| eeuss影院久久| 偷拍熟女少妇极品色| 免费在线观看成人毛片| 在线国产一区二区在线| 好男人电影高清在线观看| 亚洲av成人av| 国产高清三级在线| 在线国产一区二区在线| 国产午夜精品久久久久久一区二区三区 | 亚洲精品在线美女| 国产成人系列免费观看| 久久久久久国产a免费观看| 久久久国产成人免费| 亚洲欧美日韩高清在线视频| 国产一级毛片七仙女欲春2| 欧美午夜高清在线| 91久久精品国产一区二区成人 | 亚洲久久久久久中文字幕| 又黄又粗又硬又大视频| 国产精品影院久久| 国产精品日韩av在线免费观看| 久久久色成人| 日韩欧美在线二视频| 91在线观看av| 久久久久久人人人人人| 国产精品,欧美在线| 一个人观看的视频www高清免费观看| aaaaa片日本免费| 日日夜夜操网爽| 身体一侧抽搐| 国产精品,欧美在线| 亚洲aⅴ乱码一区二区在线播放| xxxwww97欧美| 国产探花在线观看一区二区| 全区人妻精品视频| 午夜a级毛片| 岛国在线免费视频观看| 亚洲精品456在线播放app | www.色视频.com| 一个人观看的视频www高清免费观看| 久久精品国产99精品国产亚洲性色| 两性午夜刺激爽爽歪歪视频在线观看| 在线天堂最新版资源| 热99在线观看视频| 成人午夜高清在线视频| 久久精品国产99精品国产亚洲性色| 欧美+日韩+精品| 日韩成人在线观看一区二区三区| 亚洲无线在线观看| 国产麻豆成人av免费视频| 看片在线看免费视频| 久久欧美精品欧美久久欧美| 久久精品影院6| 国产麻豆成人av免费视频| 欧美乱码精品一区二区三区| 婷婷精品国产亚洲av在线| 久久精品综合一区二区三区| 国内精品久久久久久久电影| 香蕉久久夜色| 91av网一区二区| 午夜福利视频1000在线观看| 国产av一区在线观看免费| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| 亚洲国产精品成人综合色| 99久久综合精品五月天人人| 欧美不卡视频在线免费观看| 亚洲美女视频黄频| 国产亚洲精品久久久com| 亚洲人成网站在线播| 日韩国内少妇激情av| 久久精品夜夜夜夜夜久久蜜豆| 男女视频在线观看网站免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品嫩草影院av在线观看 | 一进一出抽搐gif免费好疼| 国产一区二区亚洲精品在线观看| 在线看三级毛片| 精品久久久久久久久久久久久| 熟妇人妻久久中文字幕3abv| 国产高潮美女av| 黄片小视频在线播放| 中亚洲国语对白在线视频| 欧美黄色淫秽网站| 欧美日本视频| 国产高清视频在线播放一区| 亚洲av成人av| e午夜精品久久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 久久九九热精品免费| 搡老熟女国产l中国老女人| avwww免费| 最近在线观看免费完整版| 欧美日本视频| 亚洲片人在线观看| 男人和女人高潮做爰伦理| 9191精品国产免费久久| 亚洲片人在线观看| 一本一本综合久久| 天堂av国产一区二区熟女人妻| ponron亚洲| 欧美性感艳星| 色尼玛亚洲综合影院| 动漫黄色视频在线观看| 女生性感内裤真人,穿戴方法视频| 久久久久九九精品影院| 精品久久久久久久久久久久久| 熟妇人妻久久中文字幕3abv| 国产成人系列免费观看| 色播亚洲综合网| 观看免费一级毛片| 欧美一区二区精品小视频在线| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品亚洲一级av第二区| 亚洲内射少妇av| aaaaa片日本免费| 成人国产综合亚洲| 国产精品免费一区二区三区在线| 成人av在线播放网站| 亚洲一区二区三区不卡视频| 九九久久精品国产亚洲av麻豆| 亚洲av免费高清在线观看| 中文字幕高清在线视频| 亚洲国产精品sss在线观看| 午夜老司机福利剧场| 一级毛片女人18水好多| 一区福利在线观看| 一夜夜www| 久久人妻av系列| 嫩草影视91久久| 免费人成视频x8x8入口观看| 亚洲欧美日韩高清专用| 亚洲中文日韩欧美视频| 国产不卡一卡二| 国产亚洲av嫩草精品影院| 看黄色毛片网站| 亚洲精品成人久久久久久| 日韩av在线大香蕉| 午夜精品一区二区三区免费看| 精品一区二区三区视频在线 | 国产午夜精品久久久久久一区二区三区 | 黄色日韩在线| 偷拍熟女少妇极品色| 一级毛片女人18水好多| 久久精品国产自在天天线| 最近在线观看免费完整版| 国产精品98久久久久久宅男小说| 成人特级av手机在线观看| 成人18禁在线播放| 成人国产一区最新在线观看| 一夜夜www| 黄色丝袜av网址大全| 国产男靠女视频免费网站| 成人无遮挡网站| 久久精品国产综合久久久| 色吧在线观看| 免费看日本二区| 日本五十路高清| 日韩有码中文字幕| 色哟哟哟哟哟哟| a在线观看视频网站| 熟妇人妻久久中文字幕3abv| 亚洲片人在线观看| 女警被强在线播放| 夜夜看夜夜爽夜夜摸| 我要搜黄色片| 国产探花在线观看一区二区| 男人舔女人下体高潮全视频| 非洲黑人性xxxx精品又粗又长| 久久香蕉精品热| 亚洲av第一区精品v没综合| 精品久久久久久久久久免费视频| 国产午夜精品久久久久久一区二区三区 | 欧美精品啪啪一区二区三区| 久久久久久国产a免费观看| 观看免费一级毛片| 亚洲真实伦在线观看| 最近最新中文字幕大全免费视频| 国产精品电影一区二区三区| 亚洲真实伦在线观看| 欧美zozozo另类| 国产精品电影一区二区三区| 亚洲真实伦在线观看| 桃红色精品国产亚洲av| 成人精品一区二区免费| 精华霜和精华液先用哪个| 桃红色精品国产亚洲av| 精品福利观看| 国产精品精品国产色婷婷| 国内精品久久久久精免费| 日本 av在线| 欧美黄色淫秽网站| 12—13女人毛片做爰片一| a在线观看视频网站| 日韩国内少妇激情av| 国产一区二区亚洲精品在线观看| 级片在线观看| 日韩欧美精品v在线| 亚洲国产色片| av天堂中文字幕网| 日本免费a在线| 看黄色毛片网站| 免费高清视频大片| 国产成人av激情在线播放| 丰满人妻一区二区三区视频av | 丁香六月欧美| 免费电影在线观看免费观看| 老司机午夜福利在线观看视频| 国产爱豆传媒在线观看| 毛片女人毛片| 国产精品国产高清国产av| 久久久国产精品麻豆| 日本一本二区三区精品| 国产一区二区在线观看日韩 | 久久国产精品人妻蜜桃| 99国产极品粉嫩在线观看| 午夜精品在线福利| av福利片在线观看| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 免费高清视频大片| 国产aⅴ精品一区二区三区波| 亚洲国产精品999在线| 国产淫片久久久久久久久 | 好男人在线观看高清免费视频| 国产免费男女视频| 国产一区二区在线av高清观看| 国产精品永久免费网站| 91av网一区二区| 精品无人区乱码1区二区| 日本黄大片高清| 一区二区三区免费毛片| 日韩国内少妇激情av| 久久久国产精品麻豆| 99精品在免费线老司机午夜| 国产私拍福利视频在线观看| 久久草成人影院| 久久欧美精品欧美久久欧美| 亚洲最大成人手机在线| 白带黄色成豆腐渣| а√天堂www在线а√下载| 人人妻人人看人人澡| 亚洲熟妇中文字幕五十中出| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 美女高潮的动态| 久久精品人妻少妇| 久久久久九九精品影院| 老司机在亚洲福利影院| 国产黄色小视频在线观看| 午夜免费男女啪啪视频观看 | 亚洲,欧美精品.| 高清在线国产一区| 在线观看一区二区三区| 成年女人毛片免费观看观看9| 国产精品影院久久| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 国产精华一区二区三区| 无遮挡黄片免费观看| 国产淫片久久久久久久久 | 成人18禁在线播放| 此物有八面人人有两片| 黄色片一级片一级黄色片| 午夜影院日韩av| 久久久国产成人精品二区| 精品久久久久久,| 成年女人毛片免费观看观看9| 亚洲美女黄片视频| 欧美丝袜亚洲另类 | 热99re8久久精品国产| 日韩人妻高清精品专区| 三级毛片av免费| 99国产极品粉嫩在线观看| 一级黄片播放器| 亚洲av成人不卡在线观看播放网| 亚洲国产色片| 国产一区二区三区视频了| 免费观看的影片在线观看| 亚洲 欧美 日韩 在线 免费| 嫁个100分男人电影在线观看| 国产精品亚洲一级av第二区| 欧美黑人巨大hd| 国产黄片美女视频| 一个人观看的视频www高清免费观看| 婷婷亚洲欧美| 久久久久久久精品吃奶| 激情在线观看视频在线高清| 91久久精品国产一区二区成人 | 此物有八面人人有两片| 中文字幕精品亚洲无线码一区| 午夜福利在线在线| 最新在线观看一区二区三区| 亚洲成人精品中文字幕电影| 久久草成人影院| 一a级毛片在线观看| 法律面前人人平等表现在哪些方面| 亚洲精品日韩av片在线观看 | 91九色精品人成在线观看| 日本免费a在线| 亚洲一区二区三区不卡视频| 亚洲av一区综合| 国产乱人伦免费视频| 69av精品久久久久久| 国产精品乱码一区二三区的特点| 叶爱在线成人免费视频播放| 国产伦在线观看视频一区| 欧美性猛交╳xxx乱大交人| 亚洲av中文字字幕乱码综合| 噜噜噜噜噜久久久久久91| av专区在线播放| 国产精品三级大全| 香蕉久久夜色| 免费大片18禁| 天天一区二区日本电影三级|