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

    扇形凹穴型微通道流動與傳熱的場協(xié)同與熵產(chǎn)分析

    2019-12-20 02:09:24馮振飛
    真空與低溫 2019年6期
    關(guān)鍵詞:摩阻雷諾數(shù)扇形

    張 穎,馮振飛 ,2

    (1.廣西大學(xué)機(jī)械工程學(xué)院,南寧 530004;2.廣西大學(xué)廣西石化資源加工及過程強(qiáng)化技術(shù)重點實驗室,南寧 530004)

    0 引言

    換熱器廣泛應(yīng)用于化工、航空航天和制冷等各個工業(yè)領(lǐng)域。隨著科技的不斷進(jìn)步,各行業(yè)對換熱器的換熱效果提出了更高的要求,包括要求換熱器向微型化方向發(fā)展。常規(guī)尺度通道的換熱器向微型化發(fā)展就形成了微細(xì)通道換熱器。通道直徑越小,傳熱系數(shù)就越高[1]。微細(xì)通道換熱器較高的傳熱系數(shù)促使其成為研究者們關(guān)注的熱點。為了使微通道換熱器滿足更大的熱負(fù)荷換熱需求,研究者們不斷探索其強(qiáng)化傳熱方法。被動式強(qiáng)化傳熱不需要額外提供能源,成為研究者們主要的研究方向。

    在微細(xì)通道內(nèi)設(shè)置凹穴結(jié)構(gòu)也屬于被動式強(qiáng)化傳熱方法。Pan等[2]實驗研究了扇形凹穴型微通道換熱器的傳熱特性。Alfellag等[3]在直細(xì)微通道內(nèi)加入梯形凹穴,并在凹穴內(nèi)設(shè)置橢圓形針肋,研究尺寸改變的最優(yōu)傳熱效果,結(jié)果表明,在研究范圍內(nèi)尺寸參數(shù)最優(yōu)時的最大綜合傳熱因子約1.37。Ma等[4]研究了在側(cè)壁設(shè)置周期性噴射節(jié)流凹穴結(jié)構(gòu)的微通道熱沉的換熱效果,研究發(fā)現(xiàn),節(jié)流型微通道熱沉的綜合性能因子約為噴射型微通道熱沉的2.5倍。Ghani等[5]研究了正弦型凹穴內(nèi)置矩形肋的流動傳熱特性。Liu等[6]對設(shè)置扇形凹穴的環(huán)形微通道進(jìn)行數(shù)值模擬,研究尺寸參數(shù)對微通道傳熱和流動特性的影響。馮振飛等[7]用對比方式研究了扇形凹穴結(jié)構(gòu)對直細(xì)微通道和螺旋形微通道的流動傳熱特性影響,結(jié)果表明,凹穴不能提高螺旋形微通道的綜合性能,但提高了直細(xì)微通道的綜合性能,最大強(qiáng)化傳熱因子達(dá)1.27。賈玉婷等[8]研究了等直徑段有不同夾角的水滴形凹穴微通道的流動傳熱特性。

    從上述研究發(fā)現(xiàn),研究者們往往基于熱力學(xué)第一定律對凹穴型微通道的流動與傳熱特性進(jìn)行研究,而基于熱力學(xué)第二定律的研究較少。此外,凹穴結(jié)構(gòu)對微通道流動和傳熱影響的內(nèi)在本質(zhì)尚未明晰。本工作是在前期扇形凹穴型微通道研究[7]的基礎(chǔ)上,對扇形凹穴型微通道流動與傳熱進(jìn)行場協(xié)同和熵產(chǎn)分析。具體是基于場協(xié)同理論,分析扇形凹穴結(jié)構(gòu)參數(shù)對微通道內(nèi)速度矢量與壓力梯度及速度矢量與溫度梯度的協(xié)同程度影響。并用熱力學(xué)第二定律,分析流體流動和傳熱過程中的不可逆損失,最后用熵產(chǎn)增大數(shù)來描述扇形凹穴結(jié)構(gòu)對各參數(shù)微通道能量損失的影響。研究成果可進(jìn)一步揭示凹穴結(jié)構(gòu)對微通道流動和傳熱影響的內(nèi)在本質(zhì)。

    1 模型描述

    1.1 物理模型

    扇形凹穴型微通道(簡稱凹穴微通道,Microchannel with Cavities,MCC)是在傳統(tǒng)的光滑微通道(Micro-channel,MC)基礎(chǔ)上,加入扇形凹穴結(jié)構(gòu)而形成的,其結(jié)構(gòu)尺寸如圖1所示。微通道的寬Wch=0.5 mm,高Hch=1 mm,微通道計算域?qū)扺=1.5 mm,高H=1.5 mm,半肋寬a=0.5 mm,微通道長L=40 mm?;谏鲜龀叽绲玫轿⑼ǖ喇?dāng)量直徑Dh=2WchHch/(Wch+Hch)=0.667 mm。根據(jù)不同的扇形凹穴尺寸R(整個扇形的圓半徑)和相鄰扇形凹穴間距S,可得到5種凹穴微通道結(jié)構(gòu),分別命名為MCC-0.5-1(R=0.5 mm,S=1 mm)、MCC-0.5-1.5、MCC-0.5-2、MCC-0.4-1.5、MCC-0.6-1.5,對應(yīng)有31、21、16、21和21對凹穴,進(jìn)出口與凹穴的距離b=5 mm。為了便于對比,建立了MC模型,其尺寸與上述一致。

    圖1 凹穴微通道結(jié)構(gòu)示意圖Fig.1 Schematic diagramof micro-channelwithcavities

    1.2 數(shù)值模型

    計算域模型如圖2所示。固體域的材料為銅,其導(dǎo)熱率為401.0 W/(m·K);流體域的材料為去離子水,其密度997.0 kg/m3、黏度為8.899×10-4Pa·s、比熱為4 181.7 J/(kg·K)和導(dǎo)熱率為0.606 9 W/(m·K)。為了防止發(fā)生出口回流現(xiàn)象,在微通道出口處設(shè)置了長為40 mm的過渡段。

    為了簡化計算,假設(shè)本研究的流體的流動狀態(tài)為連續(xù)且穩(wěn)定的層流狀態(tài);在熱傳遞中流體和固體的物性穩(wěn)定不變;忽略體積力、熱輻射影響和熱量損失。由此,流體域的連續(xù)性方程、動量方程和能量方程可以簡化為:

    圖2 計算域模型Fig.2 Camputational domain moldel

    由于固體域是沒有速度的,其能量方程為:

    式(1)~(4)中:下標(biāo)s和f分別代表固體和流體;U為速度矢量,m/s;ρ為該流體在某溫度下的密度,kg/m3;μ為該流體在某溫度下的粘度,Pa·s;cp為該流體在某溫度下的比熱,J/(kg·K);p為壓力,Pa;T為溫度,K;λ為導(dǎo)熱率,W/(m·K)。

    在數(shù)值模擬計算之前需進(jìn)行邊界條件設(shè)置。微通道進(jìn)口采用進(jìn)口速度邊界條件,進(jìn)口溫度Tin=293 K,進(jìn)口速度uin為0.4~2 m/s,根據(jù)式(5),可得到雷諾數(shù)Re為298.76~1 493.80;微通道出口設(shè)為壓力出口邊界條件,出口壓力pout=0 Pa(相對環(huán)境壓力);微通道的底面采用恒熱流邊界條件,熱流密度q=400 kW/m2;固體和流體域接觸部分設(shè)置為固液交界面邊界條件,沒有速度滑移和滲透;由于熱沉內(nèi)的多微通道具有周期性,因此研究的單個微通道計算域的兩側(cè)面采用周期性邊界條件;模型剩余的壁面均為絕熱邊界條件。

    上述控制方程組用CFD軟件進(jìn)行數(shù)值模擬求解,收斂殘差為10-5。在數(shù)值求解之前需對6組模型的各個計算域進(jìn)行網(wǎng)格劃分。鑒于模型的復(fù)雜程度,采用四面體和六面體的混合網(wǎng)格對其進(jìn)行劃分,所有模型的網(wǎng)格數(shù)維持在130萬左右。根據(jù)進(jìn)口速度的變化(0.4、0.6、0.8、1.0、1.2、1.4、1.6、1.8、2.0 m/s),每個模型數(shù)值模擬9個工況,6組模型,共54個算例。

    2 數(shù)據(jù)處理

    平均摩擦阻力系數(shù)f的計算公式為:

    式中:Δp為微通道段進(jìn)出口的壓降,Pa。

    Shah等[9]發(fā)現(xiàn),矩形直通道在正在發(fā)展的層流狀態(tài)下,其平均摩擦阻力系數(shù)的理論計算式為:

    其中:

    式中:Po為泊肅葉數(shù);α為微通道的寬高比。平均傳熱系數(shù)h的計算式為:

    式中:Aw和Aif分別為加熱面和流固傳熱面的面積,m2;Tw和Tout分別為加熱面和通道出口的平均溫度,K。

    平均努賽爾數(shù)Nu的計算式為:

    根據(jù)能量守恒,微通道進(jìn)出口溫差的理論計算公式為:

    式中:Ain為微通道進(jìn)口截面面積,m2。

    根據(jù)Gao等[10]提出的場協(xié)同理論,可得到流體速度矢量與溫度梯度之間的協(xié)同關(guān)系表達(dá)式為:

    式中:β為速度矢量和溫度梯度之間的協(xié)同角。協(xié)同角β越小,場協(xié)同效果越好,越有利于傳熱。

    同理可得到流體速度矢量與壓力梯度之間的協(xié)同關(guān)系表達(dá)式為[11]:

    式中:α為速度矢量與壓力梯度之間的協(xié)同角。協(xié)同角α越小,流體流動所消耗的泵功就越少。

    式(12)和式(13)中的協(xié)同角為0~180°,為了便于分析,取絕對值來說明協(xié)同角,其在0~90°之間。即α和β的計算式分別為:

    式中:u、v和w分別為x、y和z的速度,m/s。

    微通道熱沉內(nèi)部流體之間的對流傳熱包含流動過程和傳熱過程,按照熱力學(xué)第二定律,流體流動過程中壓降造成摩擦損失耗能,對流傳熱過程中因溫差存在造成傳熱能量損失。采用熵產(chǎn)原理分析微通道流動和傳熱過程的不可逆損失。流動過程中產(chǎn)生的流動熵產(chǎn)SF定義為[12]:

    式中:m為質(zhì)量流量,kg/s。

    傳熱過程中的傳熱熵產(chǎn)SH定義為[12]:

    總熵產(chǎn)Sg為[12]:

    為了更清晰的體現(xiàn)MCC相對MC的節(jié)能程度,引入熵產(chǎn)增大數(shù)評價各個微通道熱沉的綜合性能。熵產(chǎn)增大數(shù)定義為[12]:

    式中:Ns,a代表熵產(chǎn)增大數(shù),無量綱;Sg,0為光滑微通道的總熵產(chǎn)。Ns,a值小于1表明總熵產(chǎn)相對光滑微通道而言有所減少,Ns,a值越小,不可逆損失越小。

    3 數(shù)值模擬有效性驗證

    為了驗證所采用的數(shù)值模擬方法的有效性,選擇光滑微通道的數(shù)值模擬結(jié)果和理論計算結(jié)果進(jìn)行對比。圖3為光滑微通道摩阻系數(shù)和進(jìn)出口溫差的理論計算值與模擬值的對比結(jié)果。因為微通道內(nèi)的流體狀態(tài)是正在發(fā)展的層流,所以式(6)和式(7)中的L包含過渡段長度,即L=80 mm。

    政府應(yīng)該放寬出租車行業(yè)的準(zhǔn)入機(jī)制,不僅要鼓勵更多的企業(yè)參與到出租車行業(yè)的運營當(dāng)中,將出租車數(shù)量的控制權(quán)交給市場進(jìn)行管制,還可以通過更多的法律制度規(guī)范來約束整個出租車行業(yè),更多地實施監(jiān)管的責(zé)任,保證乘客的安全,提高人們對出租車的選擇率。

    由圖3可知,在所研究的雷諾數(shù)范圍內(nèi),光滑微通道摩阻系數(shù)的理論值和模擬值比較吻合,最大誤差為15.79%;光滑微通道的進(jìn)出口溫差的理論值和模擬值非常接近,最大相對誤差為1.55%。由此可見,所采用的數(shù)值模擬方法是可靠有效的,該方法也可以用于計算分析凹穴微通道的流動和傳熱特性。

    圖3 光滑細(xì)通道摩阻系數(shù)和進(jìn)出口溫差的模擬值與理論值的對比曲線Fig.3 Comparisons between simulation and theoretical data of friction factor and temperature difference between inlet and outlet in smooth micro-channels

    4 結(jié)果與討論

    4.1 流動與傳熱特性分析

    圖4為用式(6)計算得到的雷諾數(shù)范圍內(nèi)各結(jié)構(gòu)摩阻系數(shù)的變化趨勢。由圖4可見,隨著雷諾數(shù)增加,摩阻系數(shù)逐漸變小;在同一雷諾數(shù)下,各個微通道的摩阻系數(shù)各不相同,扇形凹穴型微通道的摩阻系數(shù)均高于光滑微通道。經(jīng)計算,MCC要比MC摩阻系數(shù)平均高約53.6%;相比MC,MCC-0.5-1.5阻力變化最大,摩阻系數(shù)平均高約72.6%;MCC-0.6-1.5摩阻系數(shù)變化最小,平均高出約36.9%。

    圖4 摩阻系數(shù)f隨雷諾數(shù)Re的變化曲線Fig.4 Variation of friction factor f with Reynolds number Re

    圖5 是由式(9)算得的傳熱系數(shù)隨雷諾數(shù)的變化趨勢。由圖5知,傳熱系數(shù)隨著雷諾數(shù)變大而變大;在多數(shù)工況下,相對于MC,MCC的傳熱系數(shù)均有所增加,計算得平均高約28.95%,說明扇形凹穴結(jié)構(gòu)提高了微通道的傳熱效果。其中MCC-0.5-1.5傳熱系數(shù)最高,相比MC平均高約40.02%。

    圖5 傳熱系數(shù)h隨著雷諾數(shù)Re的變化曲線Fig.5 Variation of heat transfer factor h with Reynolds number Re

    4.2 場協(xié)同分析

    圖6 給出了6種微通道的速度矢量和壓力梯度間的協(xié)同角云圖。由圖6(a)發(fā)現(xiàn)在MC整體區(qū)域內(nèi)協(xié)同角α較小,說明速度場和壓力場協(xié)同程度很高,這由于流體在直通道內(nèi)沿程壓降是沿著流動方向下降的,壓力梯度和速度矢量的方向就比較一致。所以一定程度上協(xié)同角α的大小,能夠體現(xiàn)能量的利用程度,即流體容易從壓力高的位置流向壓力低的位置。觀察圖6中5種MCC的云圖,發(fā)現(xiàn)協(xié)同角α較MC均有所增加。顯然在局部區(qū)域內(nèi)尤其是扇形凹穴區(qū)域內(nèi),由于流動擾動作用增強(qiáng),壓降突然增大,這與光滑微通道內(nèi)壓降是沿著流動方向逐漸減小的現(xiàn)象相反。因此,MCC通道局部位置的壓降梯度方向發(fā)生改變,導(dǎo)致協(xié)同角α變大。故而,相對MC,扇形凹穴結(jié)構(gòu)增加了微通道的阻力損失,降低了能量利用程度。

    圖6 6種微通道的速度矢量和壓力梯度間的協(xié)同角α云圖(進(jìn)口速度u in=1 m/s)Fig.6 Contours of synergistic angle α between velocity vectors and pressure gradients for six micro-channels(inlet velocity u in=1 m/s)

    圖7 是6種微通道的速度矢量和溫度梯度的協(xié)同角β云圖。對比圖7各圖可見,在整個MC內(nèi),速度矢量方向和溫度梯度的方向是相反的,這點在另外5種MCC的直通道部分也很明顯。而在扇形凹穴區(qū)域,發(fā)生層流邊界層中斷、流體擾動增強(qiáng)等作用,甚至產(chǎn)生不斷循環(huán)的二次流,致使凹穴區(qū)域的速度矢量方向隨著流動擾動和二次流作用發(fā)生改變,進(jìn)而減少了局部區(qū)域內(nèi)速度矢量和溫度梯度協(xié)同角β,甚至?xí)霈F(xiàn)β極小的局部區(qū)域。所以,扇形凹穴結(jié)構(gòu)的加入改善了速度場和溫度場的協(xié)同程度,減小了協(xié)同角β,強(qiáng)化了傳熱效果。

    但是,云圖只能體現(xiàn)局部協(xié)同角β的變化情況,不能定量的描述其大小,所以引入場協(xié)同數(shù)來計算不同進(jìn)口速度下的微通道速度場和溫度場的協(xié)同程度,場協(xié)同數(shù)定義為:

    式中:Pr為普朗特數(shù),無量綱。Fc為場協(xié)同數(shù),是一個小于1的無量綱數(shù),其值越大表示速度和溫度場協(xié)同程度越高,越利于傳熱。

    圖8為場協(xié)同數(shù)隨雷諾數(shù)的變化關(guān)系。由圖8可知,隨著雷諾數(shù)Re增加,場協(xié)同數(shù)Fc逐漸變小,這是因為微通道內(nèi)流體流速過高時會削弱扇形凹穴結(jié)構(gòu)的擾動作用,二次流現(xiàn)象減少,使場協(xié)同程度下降。在同一雷諾數(shù)Re下的大多數(shù)工況時,MCC的場協(xié)同數(shù)要高于MC。顯然,扇形凹穴結(jié)構(gòu)的加入改善了速度場和溫度場的協(xié)同程度。經(jīng)過計算,在所研究范圍內(nèi),MCC-0.5-1.5的場協(xié)同數(shù)最大,平均約為0.002 94。

    圖7 6種微通道的速度矢量和溫度梯度的協(xié)同角云圖(進(jìn)口速度u in=1 m/s)Fig.7 Contours of synergistic angles β between velocity vectors and temperature gradients for six micro-channels(inlet velocity u in=1 m/s)

    圖8 場協(xié)同數(shù)F c隨雷諾數(shù)Re的變化曲線Fig.8 Variation of field synergy number F c with Reynolds number Re

    4.3 熵產(chǎn)分析

    圖9 6組微通道的流動熵產(chǎn)S F和傳熱熵產(chǎn)S H隨雷諾數(shù)Re的變化曲線Fig.9 Variation of flow entropy generation S F and Heat entropy generation S H with Reynolds number Re in six microchannels

    圖10 為總熵產(chǎn)和熵產(chǎn)增大數(shù)隨雷諾數(shù)的變化趨勢。由圖可知,隨著雷諾數(shù)增大,總熵產(chǎn)逐漸減小,因傳熱熵產(chǎn)幾乎可以代表總熵產(chǎn),所以趨勢與傳熱熵產(chǎn)基本一致;相對MC,MCC的總熵產(chǎn)均有所減小,表現(xiàn)為熵產(chǎn)增大數(shù)均小于1,說明,扇形凹穴結(jié)構(gòu)雖帶來阻力損失,但也帶來更好的傳熱效果,總體上能夠減少微通道的不可逆能量損失。根據(jù)熵產(chǎn)最小原則,計算可知MCC-0.5-1.5的熵產(chǎn)增大數(shù)最小,平均約為0.74,在Re=747附近工況下的最小熵產(chǎn)增大數(shù)約為0.69。

    圖10 總熵產(chǎn)S g和熵產(chǎn)增大數(shù)N s,a隨雷諾數(shù)Re的變化曲線Fig.10 Variation of total entropy generation S g and augmentation entropy generation number N s,a with Reynolds number Re

    5 結(jié)論

    (1)對于光滑微通道而言,扇形凹穴結(jié)構(gòu)增大了流阻,MCC的摩阻系數(shù)平均增加53.6%;扇形凹穴結(jié)構(gòu)也提升了傳熱效果,MCC的傳熱系數(shù)平均提高了28.9%。

    (2)扇形凹穴結(jié)構(gòu)改變了MCC內(nèi)局部位置壓降梯度的方向,致使速度場和壓力場協(xié)同配合程度下降,協(xié)同角α變大,扇形凹穴結(jié)構(gòu)增加了微通道的阻力損失。

    (3)流體在凹穴位置產(chǎn)生擾動和二次流作用,改變了局部區(qū)域速度場方向,速度場和溫度場的協(xié)同角β有所減小,凹穴結(jié)構(gòu)給微通道帶來更好的傳熱效果。經(jīng)過場協(xié)同數(shù)Fc的計算,MCC-0.5-1.5場協(xié)同數(shù)最大,平均為0.002 94。

    (4)熵產(chǎn)分析表明,扇形凹穴結(jié)構(gòu)給微通道帶來的阻力損失遠(yuǎn)小于給微通道帶來的傳熱效果,流動熵產(chǎn)的增大量遠(yuǎn)小于傳熱熵產(chǎn)的減小量,MCC的熵產(chǎn)增大數(shù)均小于1。其中,MCC-0.5-1.5的熵產(chǎn)增大數(shù)最小,平均為0.74。

    猜你喜歡
    摩阻雷諾數(shù)扇形
    各種各樣的扇形
    探源拓思融會貫通
    ———《扇形的認(rèn)識》教學(xué)廖
    市政橋梁預(yù)應(yīng)力管道摩阻系數(shù)測試研究
    江西建材(2018年4期)2018-04-10 12:37:20
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    復(fù)扇形指標(biāo)集上的分布混沌
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    民機(jī)高速風(fēng)洞試驗的阻力雷諾數(shù)效應(yīng)修正
    計算隱式摩阻系數(shù)方程數(shù)值解的簡便方法
    考慮扶正器影響的套管摩阻計算方法研究
    97热精品久久久久久| 国内精品久久久久精免费| 99久久人妻综合| 一边摸一边抽搐一进一小说| 日韩高清综合在线| 天天躁日日操中文字幕| 99九九线精品视频在线观看视频| 一边亲一边摸免费视频| 麻豆国产97在线/欧美| 国产亚洲91精品色在线| 久久久久久久久久黄片| 国产精品一区二区三区四区久久| 国产精品人妻久久久影院| 一本精品99久久精品77| 一本精品99久久精品77| 日本一本二区三区精品| 日韩视频在线欧美| 成人午夜精彩视频在线观看| 精品久久久噜噜| 日韩大尺度精品在线看网址| 可以在线观看毛片的网站| 午夜视频国产福利| 久久这里有精品视频免费| 天堂网av新在线| 亚洲欧美日韩高清在线视频| 日韩大尺度精品在线看网址| a级一级毛片免费在线观看| 午夜视频国产福利| 日韩欧美在线乱码| 欧美变态另类bdsm刘玥| 国产精品不卡视频一区二区| 欧美+日韩+精品| 波野结衣二区三区在线| 不卡一级毛片| 国产精品综合久久久久久久免费| 寂寞人妻少妇视频99o| 久久人人精品亚洲av| 天堂影院成人在线观看| 十八禁国产超污无遮挡网站| 男女那种视频在线观看| 美女被艹到高潮喷水动态| 国产伦理片在线播放av一区 | 国产激情偷乱视频一区二区| a级一级毛片免费在线观看| 伦精品一区二区三区| h日本视频在线播放| 成年女人看的毛片在线观看| 此物有八面人人有两片| 亚洲精品成人久久久久久| 国产高清激情床上av| 亚洲av.av天堂| 国产蜜桃级精品一区二区三区| 91麻豆精品激情在线观看国产| 哪个播放器可以免费观看大片| 久久99热这里只有精品18| 亚洲国产欧洲综合997久久,| 欧美性猛交黑人性爽| 亚洲av一区综合| 色播亚洲综合网| 亚洲高清免费不卡视频| 日韩亚洲欧美综合| 韩国av在线不卡| 99久国产av精品| 中文字幕人妻熟人妻熟丝袜美| 国产精品电影一区二区三区| 亚洲成av人片在线播放无| 国产精品国产高清国产av| av.在线天堂| 亚洲激情五月婷婷啪啪| 国产午夜精品久久久久久一区二区三区| 亚洲成人精品中文字幕电影| 久久精品久久久久久噜噜老黄 | 99久久精品热视频| 天堂√8在线中文| 久久综合国产亚洲精品| 日韩欧美国产在线观看| 久久久精品大字幕| 女人被狂操c到高潮| 99热这里只有精品一区| 久久久久久久亚洲中文字幕| 亚洲成av人片在线播放无| 人人妻人人澡人人爽人人夜夜 | 狠狠狠狠99中文字幕| 春色校园在线视频观看| 秋霞在线观看毛片| 99热这里只有精品一区| 亚洲人成网站高清观看| 午夜免费男女啪啪视频观看| 人妻制服诱惑在线中文字幕| 你懂的网址亚洲精品在线观看 | 国产精品福利在线免费观看| www.色视频.com| 久久久久久久午夜电影| 男女视频在线观看网站免费| 毛片女人毛片| 久久精品国产亚洲网站| 日本三级黄在线观看| 看黄色毛片网站| h日本视频在线播放| 日韩欧美 国产精品| 久久久午夜欧美精品| 亚洲美女视频黄频| 丰满乱子伦码专区| 欧美xxxx黑人xx丫x性爽| 中文字幕精品亚洲无线码一区| 国产久久久一区二区三区| 成人美女网站在线观看视频| 国产午夜精品论理片| 在线观看免费视频日本深夜| 99国产极品粉嫩在线观看| 亚洲美女视频黄频| a级毛片a级免费在线| 男女视频在线观看网站免费| 三级毛片av免费| 日本三级黄在线观看| 中文字幕av在线有码专区| 变态另类成人亚洲欧美熟女| 日产精品乱码卡一卡2卡三| 亚洲综合色惰| 成人漫画全彩无遮挡| 亚洲最大成人手机在线| 国产亚洲欧美98| 天堂网av新在线| 国内精品一区二区在线观看| 12—13女人毛片做爰片一| 少妇丰满av| 亚洲国产欧美人成| 六月丁香七月| 国产视频内射| 国产色爽女视频免费观看| 91狼人影院| 亚洲第一区二区三区不卡| 美女大奶头视频| 久久人人爽人人爽人人片va| 日日摸夜夜添夜夜爱| 少妇熟女欧美另类| 尤物成人国产欧美一区二区三区| 中文资源天堂在线| 小说图片视频综合网站| 91av网一区二区| 有码 亚洲区| 好男人视频免费观看在线| 日本黄色片子视频| 久久久久久久亚洲中文字幕| 色综合色国产| 亚洲乱码一区二区免费版| 最新中文字幕久久久久| 久久久精品94久久精品| 国产淫片久久久久久久久| 校园春色视频在线观看| 国产精品野战在线观看| 婷婷精品国产亚洲av| 欧美日本亚洲视频在线播放| 亚洲自偷自拍三级| 欧美激情在线99| 成年免费大片在线观看| 黄色欧美视频在线观看| 一个人看的www免费观看视频| 国产女主播在线喷水免费视频网站 | 久久久久久久亚洲中文字幕| 不卡视频在线观看欧美| 日韩人妻高清精品专区| 99热只有精品国产| а√天堂www在线а√下载| 国产人妻一区二区三区在| 99久久成人亚洲精品观看| 久久精品国产亚洲av涩爱 | 欧美xxxx黑人xx丫x性爽| 女人被狂操c到高潮| 国产色爽女视频免费观看| 三级经典国产精品| 国产精品无大码| a级毛片免费高清观看在线播放| 日韩高清综合在线| 免费黄网站久久成人精品| kizo精华| 我的老师免费观看完整版| 亚洲综合色惰| 免费人成视频x8x8入口观看| 少妇熟女欧美另类| 欧美潮喷喷水| 网址你懂的国产日韩在线| 国产一区二区三区av在线 | 午夜视频国产福利| 国产极品精品免费视频能看的| 免费黄网站久久成人精品| 免费av观看视频| 日韩国内少妇激情av| 老熟妇乱子伦视频在线观看| 久久久午夜欧美精品| 欧美精品国产亚洲| 99久久精品热视频| 午夜亚洲福利在线播放| 91久久精品国产一区二区三区| videossex国产| 男人狂女人下面高潮的视频| 国产精品野战在线观看| av.在线天堂| www.色视频.com| 女人十人毛片免费观看3o分钟| 99热这里只有是精品在线观看| 成人一区二区视频在线观看| 高清午夜精品一区二区三区 | 国产精品伦人一区二区| 99久久精品一区二区三区| 成人特级黄色片久久久久久久| 少妇丰满av| 欧美日韩乱码在线| 国产精品一二三区在线看| 国产探花在线观看一区二区| 久久精品91蜜桃| 久久婷婷人人爽人人干人人爱| 一个人免费在线观看电影| 国产黄色视频一区二区在线观看 | 久久这里有精品视频免费| 色哟哟哟哟哟哟| 99热只有精品国产| 亚洲精品久久国产高清桃花| 晚上一个人看的免费电影| 国产伦一二天堂av在线观看| 国产精品综合久久久久久久免费| 亚洲精品色激情综合| 亚洲高清免费不卡视频| 国产精品久久电影中文字幕| 欧美不卡视频在线免费观看| 看非洲黑人一级黄片| 成人漫画全彩无遮挡| 99热全是精品| 赤兔流量卡办理| 日韩欧美 国产精品| 长腿黑丝高跟| 男插女下体视频免费在线播放| 国产乱人偷精品视频| 色综合色国产| 欧美bdsm另类| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕 | 精品国内亚洲2022精品成人| 日本黄色视频三级网站网址| 国产精品久久久久久精品电影| 久久韩国三级中文字幕| 蜜桃亚洲精品一区二区三区| 成人永久免费在线观看视频| 亚洲乱码一区二区免费版| 听说在线观看完整版免费高清| 99热6这里只有精品| 国产在线男女| 大型黄色视频在线免费观看| 少妇的逼好多水| 国产精品久久久久久精品电影小说 | 亚洲熟妇中文字幕五十中出| 国产精品美女特级片免费视频播放器| 啦啦啦观看免费观看视频高清| 久久精品91蜜桃| 国产高清视频在线观看网站| 日韩亚洲欧美综合| 日韩视频在线欧美| 亚洲成人久久性| 99在线人妻在线中文字幕| 亚洲国产色片| 欧美最黄视频在线播放免费| 国产女主播在线喷水免费视频网站 | 一级黄片播放器| 97超碰精品成人国产| 久久99精品国语久久久| 国产精品精品国产色婷婷| 欧美色欧美亚洲另类二区| 搞女人的毛片| 亚洲成人中文字幕在线播放| 丝袜美腿在线中文| 99久国产av精品| h日本视频在线播放| 亚洲一级一片aⅴ在线观看| 嫩草影院入口| 亚洲久久久久久中文字幕| 成年版毛片免费区| 热99在线观看视频| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3| 变态另类成人亚洲欧美熟女| 美女国产视频在线观看| 能在线免费看毛片的网站| 十八禁国产超污无遮挡网站| 亚洲第一区二区三区不卡| 长腿黑丝高跟| 狂野欧美白嫩少妇大欣赏| 久久久欧美国产精品| 成人三级黄色视频| 国产成人a∨麻豆精品| 国产在视频线在精品| 寂寞人妻少妇视频99o| 久久久久久国产a免费观看| 午夜激情福利司机影院| 国产伦一二天堂av在线观看| av免费观看日本| 夜夜爽天天搞| 亚洲国产精品久久男人天堂| 一级黄色大片毛片| 美女国产视频在线观看| 久久久色成人| 人人妻人人看人人澡| 日韩欧美精品v在线| 91在线精品国自产拍蜜月| 日本三级黄在线观看| 在线天堂最新版资源| 国产精品一区www在线观看| 99热精品在线国产| 日产精品乱码卡一卡2卡三| 日本免费a在线| 精品久久久久久成人av| 人体艺术视频欧美日本| 亚洲自偷自拍三级| 欧美xxxx性猛交bbbb| 国产伦理片在线播放av一区 | 精品免费久久久久久久清纯| 九九爱精品视频在线观看| 在线天堂最新版资源| 变态另类丝袜制服| 久久人人爽人人片av| 国产免费男女视频| 久久这里只有精品中国| 美女cb高潮喷水在线观看| 久久这里有精品视频免费| 人人妻人人澡欧美一区二区| 欧美激情在线99| 观看免费一级毛片| 一进一出抽搐gif免费好疼| 青春草亚洲视频在线观看| 青春草国产在线视频 | www日本黄色视频网| 一本久久精品| 成人午夜精彩视频在线观看| 91av网一区二区| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩高清在线视频| 日韩精品青青久久久久久| 国产黄片视频在线免费观看| 亚洲图色成人| 欧美三级亚洲精品| 国产精品爽爽va在线观看网站| 老司机福利观看| 午夜福利在线观看免费完整高清在 | 18禁裸乳无遮挡免费网站照片| 啦啦啦韩国在线观看视频| 国产av在哪里看| av卡一久久| 一级毛片aaaaaa免费看小| 国产av在哪里看| 精品久久久久久久末码| 亚洲欧美日韩高清在线视频| av天堂在线播放| 两个人视频免费观看高清| 国产亚洲av片在线观看秒播厂 | 亚洲欧美日韩卡通动漫| 日本一本二区三区精品| 波多野结衣高清无吗| 亚洲国产精品成人综合色| 亚洲欧美日韩高清在线视频| 日产精品乱码卡一卡2卡三| 91久久精品国产一区二区成人| 午夜精品一区二区三区免费看| 久久99精品国语久久久| 18禁在线无遮挡免费观看视频| 国产成人精品一,二区 | 免费在线观看成人毛片| 99久久九九国产精品国产免费| 亚洲av成人av| 婷婷色av中文字幕| 夜夜看夜夜爽夜夜摸| 最近视频中文字幕2019在线8| 亚洲精品乱码久久久v下载方式| 一区二区三区高清视频在线| 日本免费a在线| 丝袜喷水一区| av在线观看视频网站免费| 18禁在线无遮挡免费观看视频| 久久这里有精品视频免费| eeuss影院久久| 国产片特级美女逼逼视频| 欧美日韩国产亚洲二区| 国产蜜桃级精品一区二区三区| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 国产成人精品婷婷| 女的被弄到高潮叫床怎么办| 国产精品一区二区性色av| 超碰av人人做人人爽久久| 天堂影院成人在线观看| 2021天堂中文幕一二区在线观| videossex国产| 91精品一卡2卡3卡4卡| 久久这里只有精品中国| 51国产日韩欧美| 免费av不卡在线播放| 国产精品爽爽va在线观看网站| 国内久久婷婷六月综合欲色啪| 大型黄色视频在线免费观看| 99久久无色码亚洲精品果冻| 麻豆一二三区av精品| 婷婷六月久久综合丁香| 最近的中文字幕免费完整| 国产精品一区二区在线观看99 | 又爽又黄a免费视频| 日本五十路高清| 日韩欧美精品免费久久| 日韩欧美在线乱码| 色吧在线观看| 蜜桃亚洲精品一区二区三区| 国产免费一级a男人的天堂| 少妇猛男粗大的猛烈进出视频 | 国产成人freesex在线| 成人国产麻豆网| av女优亚洲男人天堂| 久久久久久久久久久免费av| 乱人视频在线观看| 高清在线视频一区二区三区 | 婷婷精品国产亚洲av| 中文字幕熟女人妻在线| 观看免费一级毛片| av免费观看日本| 国产美女午夜福利| 女同久久另类99精品国产91| 免费观看人在逋| 欧美最黄视频在线播放免费| 国产高清激情床上av| 99热这里只有是精品在线观看| 国产一级毛片七仙女欲春2| 亚洲欧美日韩卡通动漫| 小蜜桃在线观看免费完整版高清| 99久久中文字幕三级久久日本| 成人二区视频| 日韩欧美 国产精品| 美女内射精品一级片tv| 久久久久免费精品人妻一区二区| 亚洲精品色激情综合| 欧美成人免费av一区二区三区| 午夜福利高清视频| 高清毛片免费观看视频网站| 久久久久久久久久成人| av又黄又爽大尺度在线免费看 | 春色校园在线视频观看| 亚洲最大成人中文| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| avwww免费| 国产精品精品国产色婷婷| 干丝袜人妻中文字幕| 精品一区二区三区人妻视频| 国产精品久久久久久av不卡| 亚洲无线观看免费| 国产亚洲91精品色在线| 99久久人妻综合| 美女大奶头视频| 午夜激情欧美在线| 亚洲第一区二区三区不卡| 卡戴珊不雅视频在线播放| 在线播放无遮挡| 九色成人免费人妻av| 国产色爽女视频免费观看| 午夜a级毛片| 亚洲成a人片在线一区二区| 亚洲成人久久性| 18禁裸乳无遮挡免费网站照片| 日韩,欧美,国产一区二区三区 | 深爱激情五月婷婷| 成人一区二区视频在线观看| 美女国产视频在线观看| 国产精品蜜桃在线观看 | 综合色av麻豆| 国内精品久久久久精免费| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 国产精品福利在线免费观看| 乱系列少妇在线播放| 亚洲第一区二区三区不卡| 18禁在线无遮挡免费观看视频| 国产精品乱码一区二三区的特点| or卡值多少钱| 99精品在免费线老司机午夜| 国产精品国产三级国产av玫瑰| 一卡2卡三卡四卡精品乱码亚洲| 村上凉子中文字幕在线| 非洲黑人性xxxx精品又粗又长| 悠悠久久av| 久久国内精品自在自线图片| 变态另类成人亚洲欧美熟女| 欧美zozozo另类| 精品少妇黑人巨大在线播放 | 午夜亚洲福利在线播放| 18+在线观看网站| 色哟哟·www| 黄色欧美视频在线观看| 国产精品,欧美在线| 日本三级黄在线观看| 少妇丰满av| 99久国产av精品| 中文欧美无线码| 麻豆精品久久久久久蜜桃| 中文字幕av成人在线电影| 亚洲国产色片| 精品熟女少妇av免费看| 国产成人影院久久av| 免费av不卡在线播放| 亚洲国产精品成人综合色| 不卡视频在线观看欧美| 我要搜黄色片| 中国美白少妇内射xxxbb| 亚洲无线在线观看| 国产在线男女| 久久久久久国产a免费观看| 免费观看精品视频网站| 亚洲国产欧美在线一区| 最近最新中文字幕大全电影3| 亚洲综合色惰| 国产av一区在线观看免费| 99久久精品热视频| 亚洲图色成人| 热99在线观看视频| 九九爱精品视频在线观看| 久久精品国产自在天天线| 特大巨黑吊av在线直播| 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| 精品人妻视频免费看| 狂野欧美白嫩少妇大欣赏| 在线观看av片永久免费下载| 在线免费十八禁| 一区二区三区高清视频在线| 亚洲欧美中文字幕日韩二区| 国产一级毛片七仙女欲春2| 午夜精品在线福利| 国产午夜福利久久久久久| 久久九九热精品免费| 国产一级毛片七仙女欲春2| 日韩亚洲欧美综合| 悠悠久久av| 可以在线观看的亚洲视频| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 人妻夜夜爽99麻豆av| 人妻少妇偷人精品九色| 国内少妇人妻偷人精品xxx网站| 日韩 亚洲 欧美在线| 国产极品天堂在线| 美女 人体艺术 gogo| av在线观看视频网站免费| 老司机福利观看| 免费电影在线观看免费观看| 免费人成在线观看视频色| 国产一区二区三区在线臀色熟女| 日韩成人伦理影院| 国产高清有码在线观看视频| 亚洲国产精品合色在线| 日日摸夜夜添夜夜爱| 一级黄色大片毛片| 男人的好看免费观看在线视频| 亚洲经典国产精华液单| 久久久久九九精品影院| 高清午夜精品一区二区三区 | h日本视频在线播放| 日韩欧美精品免费久久| 久久久久久久久久久免费av| 一边摸一边抽搐一进一小说| 国产一级毛片七仙女欲春2| 99在线视频只有这里精品首页| 精品熟女少妇av免费看| 久久久国产成人免费| 亚洲精品乱码久久久久久按摩| 麻豆乱淫一区二区| 少妇人妻一区二区三区视频| 人妻久久中文字幕网| 99久久精品热视频| av在线老鸭窝| 99在线人妻在线中文字幕| 欧美变态另类bdsm刘玥| 美女xxoo啪啪120秒动态图| 国产探花极品一区二区| 日韩av在线大香蕉| 欧美日韩综合久久久久久| 亚洲四区av| 久久精品人妻少妇| 欧美色欧美亚洲另类二区| www日本黄色视频网| 天堂影院成人在线观看| 亚洲熟妇中文字幕五十中出| 麻豆一二三区av精品| 偷拍熟女少妇极品色| 国产精品一区二区三区四区免费观看| 久久午夜福利片| 免费观看在线日韩| 国产美女午夜福利| 黄片无遮挡物在线观看| 精品久久久久久成人av| 美女黄网站色视频| 成人欧美大片| 最近视频中文字幕2019在线8| 亚洲欧美日韩东京热| 精品久久久久久久末码| 国内精品美女久久久久久| 国产黄色视频一区二区在线观看 | 亚洲精品乱码久久久久久按摩| 秋霞在线观看毛片| 3wmmmm亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 哪里可以看免费的av片| 欧美日本视频| 两个人视频免费观看高清| 亚洲成人av在线免费| 国产淫片久久久久久久久| 在线观看一区二区三区| 级片在线观看| 色综合站精品国产| 国产精品99久久久久久久久|