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

    漏斗狀儲熱水箱邊界斜度對熱量損失的影響

    2020-03-26 06:13:10苑中顯郭占全
    可再生能源 2020年3期
    關(guān)鍵詞:側(cè)壁對流壁面

    張 琦,苑中顯,郭占全

    (北京工業(yè)大學(xué) 環(huán)境與能源工程學(xué)院,北京 100124)

    0 引言

    太陽能是一種安全、清潔且分布廣泛的可再生能源。為了充分地利用太陽能,在沒有熱量需求時,對太陽能進行收集并加以儲存,從而提高太陽能的全年利用率。目前,太陽能長周期跨季節(jié)儲熱技術(shù)為研究熱點[1]~[3]。

    帶有儲熱子系統(tǒng)的太陽能供熱系統(tǒng)如圖1所示。

    圖1 帶有儲熱子系統(tǒng)的太陽能供熱系統(tǒng)Fig.1 Trans-seasonal solar heating system with a heat storage tank

    該系統(tǒng)中的儲熱水箱埋于地下。在非供暖季,該系統(tǒng)將太陽能集熱器所獲得的熱量以顯熱的形式存儲于儲熱水箱內(nèi);在供暖季,將儲熱水箱中的熱水送至建筑末端,從而實現(xiàn)跨季節(jié)供暖。儲熱水箱的側(cè)壁斜度會影響自身的熱量儲存效率,這是由于貼近壁面處的水受壁面冷卻,密度隨之增大,當(dāng)水在水箱側(cè)壁自上而下流動時,若水箱側(cè)壁為豎直面,則水流不會受到壁面的阻礙,流速較大,產(chǎn)生的熱量損失也較大;若水箱側(cè)壁為傾斜面,則水流會受到傾斜壁面的阻礙,流速相對減小,產(chǎn)生的熱量損失也較小,因此傾斜的壁面有利于提高儲熱水箱的熱量儲存效率。

    目前,鮮有研究儲熱水箱壁面傾斜度對自身自然對流熱損失的影響。本文主要研究埋于地下的漏斗狀儲熱水箱(以下簡稱為儲熱水箱)側(cè)壁斜度對其自然對流換熱所產(chǎn)生的熱損失的影響,并將儲熱水箱簡化為由土壤包裹的封閉腔體,以便進行模擬研究。

    1 物理模型

    儲熱水箱周圍的土壤溫度遠低于其內(nèi)部的熱水溫度。儲熱水箱中的熱水會與土壤發(fā)生熱交換,導(dǎo)致儲熱水箱內(nèi)的水溫隨之降低,自身溫度場也變得不均勻,使得水箱內(nèi)的水產(chǎn)生流動,形成自然對流。當(dāng)儲熱水箱側(cè)壁存在一定的傾斜度時,該壁面會阻礙流體的流動,壁面傾斜度不同,此阻礙作用也不同。改變側(cè)壁與地面之間的夾角(水箱側(cè)壁傾角)θ,可得到不同的自然對流換熱強度。

    儲熱水箱的模型圖見圖2。

    圖2 儲熱水箱模型圖Fig.2 Numerical domain of the water tank of funnel shape

    由圖2可知,本文的儲熱水箱呈四棱臺結(jié)構(gòu)。令該水箱上、下底面的邊長分別為L1,L2,高為H;土壤的長、寬、高分別為a,b,h。本文將儲熱水箱的下底面幾何尺寸和高設(shè)為定值,通過改變θ,構(gòu)造出不同側(cè)壁傾角的水箱。

    模擬過程中,儲熱水箱的相關(guān)參數(shù)見表1。

    表1 儲熱水箱的相關(guān)參數(shù)Table 1 Related parameters of hot water storage tank

    計算區(qū)域中,水和土壤的物性參數(shù)隨溫度的變化很小,故視為定值。其中,水的密度為979.8 kg/m3,定壓比熱為 4 214 J/(kg·K),導(dǎo)熱系數(shù)為0.682 W/(m·K),動力粘度η為0.000 3 Pa·s,容積膨脹系數(shù)為0.000 5;土壤的密度為2 000 kg/m3,定壓比熱為 1 010 J/(kg·K),導(dǎo)熱系數(shù)為 1.16 W/(m·K)[4]。

    2 數(shù)學(xué)模型

    式中:ρ0為土壤密度;為土壤被加熱后的溫度;TS為土壤初始溫度;β為土壤的體積膨脹系數(shù)。

    2.1 控制方程

    在上述假設(shè)的基礎(chǔ)上,整個計算區(qū)域內(nèi)的流體流動以及傳熱過程滿足以下微分方程。

    連續(xù)性方程為

    式中:u為流體流速;t為時間;ρ為流體密度。

    動量方程為

    式中:p為流體壓力;T為流體溫度;T0為參考溫度;μ為流體運動粘度;β為流體容積膨脹系數(shù);μt為流體湍動粘度。

    能量方程為

    式中:σT為湍流普朗特數(shù),取0.9。

    湍動能方程為

    式中:k為湍流動能;σk為湍流動能k相應(yīng)的普朗特數(shù),取1;ε為湍流耗散率;Gk為流體平均速度梯度所引起的湍流動能的產(chǎn)生項。

    式中:C1ε為小渦拉伸項,取 1.44;C2ε為黏性破壞項,取1.92;σε為ε相應(yīng)的普朗特數(shù),取1.3。

    計算區(qū)域內(nèi)流體湍動粘度μt的計算式為

    對于計算區(qū)域內(nèi)流體平均速度梯度所引起的湍動能k,其產(chǎn)生項Gk的表達式為

    2.2 初始條件

    儲熱水箱中水的初始溫度T0=368K,外圍土壤的初始溫度TS=295 K,水與土壤之間耦合換熱。

    2.3 邊界條件

    遠離水箱處的土壤溫度不會受到水箱影響,故設(shè)置土壤模型的邊界為絕熱邊界。

    3 數(shù)值求解

    3.1 網(wǎng)格劃分

    本文用CFD前處理軟件Gambit建立儲熱水箱和外圍土壤的三維模型并生成三維實體計算網(wǎng)格。整個計算區(qū)域劃分為水箱和土壤兩部分:采用非均分網(wǎng)格在水體和土壤交界處進行局部細(xì)化,遠離水箱的土壤受水箱影響很小,因此,該處的網(wǎng)格較為稀疏。由于控制了水箱的底部邊長和高,各傾角水箱的體積不同,生成的網(wǎng)格數(shù)也不同。

    當(dāng)儲熱水箱側(cè)壁傾角θ為60°時,儲熱水箱和外圍土壤三維模型的網(wǎng)格劃分如圖3所示。

    圖3 當(dāng)θ為60°時,儲熱水箱和外圍土壤三維模型的網(wǎng)格劃分情況Fig.3 When θ=60 °,meshing of 3D model of hot water storage tank and surrounding soil

    3.1.1 網(wǎng)格無關(guān)性驗證

    當(dāng) θ為 90°,網(wǎng)格數(shù)分別為 28 061,36 286,50 419,75 993時,通過計算得到儲熱水箱的平均溫度,而后根據(jù)計算結(jié)果進行網(wǎng)格無關(guān)性驗證。網(wǎng)格無關(guān)性驗證結(jié)果見圖4。

    圖4 網(wǎng)格無關(guān)性驗證結(jié)果圖Fig.4 The average tank temperature changing with the grid number

    由圖4可知,當(dāng)網(wǎng)格數(shù)由36 286逐漸增加至50 419時,儲熱水箱平均溫度的變化量小于0.01數(shù)量級,因此,可認(rèn)為當(dāng)網(wǎng)格數(shù)為50 419時,已達到網(wǎng)格無關(guān)。

    本文還分別對 θ為 15,30,45,60,75,90 °條件下的模型進行網(wǎng)格無關(guān)性驗證。驗證結(jié)果表明,當(dāng) θ分別為 15,30,45,60,75,90 °時,本文模型的網(wǎng)格數(shù)分別應(yīng)取 1 113 027,3 649 938,416 826,207 819,98 221,50 419。

    3.1.2 時間步長無關(guān)性驗證

    通過計算發(fā)現(xiàn),當(dāng)時間步長分別為1,0.5,0.1,0.05s時,儲熱水箱的平均溫度分別為367.407,367.526,367.655,367.657 K。由此可見,當(dāng)時間步長為由0.1 s逐漸減小至0.05 s時,儲熱水箱的平均溫度變化量在0.01數(shù)量級,因此,本文取時間步長為0.1 s。

    3.2 求解器設(shè)置

    本文采用的數(shù)值方法為基于原始變量的SIMPLEC方法,壓力項采用PRESTO!方案進行離散,動量方程和能量方程采用二階迎風(fēng)格式進行離散。

    土壤與水體采取耦合換熱,這樣能夠避免將壁面設(shè)置為恒溫或恒熱流邊界條件時,土壤對流換熱模擬結(jié)果的誤差。

    4 數(shù)值模擬結(jié)果與分析

    圖5為當(dāng)θ為60°時,儲熱水箱及其周圍土壤的溫度分布情況。由圖5可知,土壤的溫度低于儲熱水箱,導(dǎo)致儲熱水箱的壁面溫度較低,儲熱水箱壁面處的水被冷卻,密度變大并沿壁面向下流動,遠離壁面處水的溫度較高,密度也較小,并填補壁面處向下流動的水,所以水箱中存在自然對流換熱。被冷卻的水沉積于底部,使儲熱水箱出現(xiàn)溫度分層現(xiàn)象,水箱底部水的溫度最低。當(dāng)儲熱水箱底部的水逐漸遠離該處,溫度逐漸升高。綜上可知,本文的儲熱水箱可以實現(xiàn)分層取水,具體的取水方式:需要低溫水時,從儲熱水箱的下部區(qū)域取水;需要高溫水時,從儲熱水箱的上部區(qū)域取水。

    圖5 當(dāng)θ為60°時,儲熱水箱及其周圍土壤的溫度分布情況Fig.5 Temperature distribution of the water tank and the surrounding soil as θ=60 °

    4.1 傾角對流場的影響

    流體沿儲熱水箱壁面的自然對流換熱強度會受到自身運動高度的影響。儲熱水箱側(cè)壁傾角會影響流體沿壁面的運動高度。建立儲熱水箱剖面的坐標(biāo)系如圖6所示。

    圖6 儲熱水箱剖面的坐標(biāo)系Fig.6 The coordinate system of the water tank for the numerical simulation

    圖7為當(dāng)儲熱水箱側(cè)壁傾角θ為15~90°時,該水箱內(nèi)部流體的流動情況。由圖7可知,當(dāng)θ較?。é?15,30 °)時,接近壁面處出現(xiàn)多個渦旋(當(dāng)θ=15°時,水箱兩側(cè)各存在3個渦旋;當(dāng)θ=30°時,水箱兩側(cè)各存在2個渦旋)。這是由于當(dāng)θ較小時,水箱壁面斜度較平緩,當(dāng)流體沿壁面向下流動時,受到傾斜壁面的阻礙作用較大,因此生成渦旋較多。渦旋越多,渦旋附近水的流速越小,導(dǎo)致水箱壁面的對流換熱強度減弱;隨著θ逐漸增大(θ=45~90°),水箱兩側(cè)的渦旋各減少至1個,且渦旋逐漸向壁面靠近。由圖6,7可知,當(dāng)θ=45,60,75,90 °時,儲熱水箱右側(cè)渦旋中心的坐標(biāo)分別為 (9.48,0.34),(8,-0.26),(6.85,-0.83),(5.9,-0.71),該渦旋與儲熱水箱右側(cè)壁面之間的垂直距離分別為2,1.5,1.1,1.1 m。這是由于隨著θ逐漸增大,水箱壁面對流動的阻礙作用逐漸減弱,導(dǎo)致渦旋與儲熱水箱壁面之間的距離大體上逐漸減小。由圖7還可看出,渦旋距離壁面越近,壁面處水的速度梯度和流動速度越大,從而導(dǎo)致水箱壁面的對流換熱強度越大。綜上可知,隨著θ逐漸增大,水箱壁面的對流換熱強度逐漸增大。

    圖7 當(dāng)θ為15~90°時,該水箱內(nèi)部流體的流動情況Fig.7 The streamline pattern of the natural convection flow in the tank as θ varies from 15 °to 90 °

    4.2 溫度、速度分布情況

    本文通過數(shù)值模擬計算得到,不同θ條件下,儲熱水箱中心截面(x=0)處溫度場的分布情況如圖8所示。

    圖8 不同θ條件下,儲熱水箱中心截面處(x=0)的溫度場分布情況Fig.8 Distribution of the temperature in the tank at x=0 palne as θ varies

    由圖8可知,儲熱水箱頂部、底部區(qū)域溫度較低,其他區(qū)域水的溫度大體一致。由于儲熱水箱頂部和底部區(qū)域受土壤冷卻作用明顯,導(dǎo)致這兩個區(qū)域中水的溫度較低。此外,由于儲熱水箱底部存在溫度分層情況,導(dǎo)致該區(qū)域溫度梯度較大。

    由圖8還可看出,當(dāng)θ不同時,儲熱水箱底部區(qū)域的溫度也不同。當(dāng)θ較大(60~90°)時,儲熱水箱底部的溫度較低,為362.7~363.5 K,則儲熱水箱底部區(qū)域的溫度梯度較大,當(dāng)θ較?。?5~45°)時,儲熱水箱底部的溫度較高,為365.7~366.6 K,則儲熱水箱底部區(qū)域的溫度梯度較小。

    不同θ條件下,儲熱水箱中心截面(x=0)處速度場的分布情況如圖9所示。

    圖9 不同θ條件下,儲熱水箱中心截面處(x=0)的速度場分布情況Fig.9 Distribution of the velocity in the tank at x=0 palne as θ varies

    由圖 9 可知,當(dāng) θ較?。?5,30 °)時,水箱內(nèi)水的流動速度自上而下(不包含底部區(qū)域)逐漸增大,這是由于當(dāng)水從儲熱水箱頂部向下流動時,水箱頂部(y=5 m處)水的流速較小,為 0~0.015 m/s,隨后在重力加速度的影響下,水的流動速度逐漸增大。當(dāng)水流至水箱底部時,由于底面阻擋,速度有所降低。當(dāng)θ由45°逐漸增大至90°時,水箱內(nèi)水的流動速度整體上逐漸增大,此外,水箱內(nèi)水的流動速度受渦旋影響較大,與高度之間的關(guān)系不明顯。隨著θ逐漸增大,儲熱水箱壁面對流動的阻礙能力逐漸減弱,則流體的流動速度大體上隨傾角θ的增大而增大。

    由圖9還可看出,當(dāng)θ為75,90°時,儲熱水箱底部區(qū)域 (y=-5 m附近)水的流動速度為0~0.005 m/s;當(dāng)θ為15~60°時,儲熱水箱底部水的流動速度為0.006~0.033 m/s。總體來看,對于θ為75,90°的儲熱水箱,其底部區(qū)域水的流速明顯小于其他傾角條件下的流速,這是由于當(dāng)水流到儲熱水箱底面后,流動方向會發(fā)生改變,底面對流動產(chǎn)生堵塞作用,導(dǎo)致水的流速減小,并且這種阻礙能力會隨著θ的增大而增強,當(dāng)θ大于60°時,此影響尤為明顯。這一現(xiàn)象也會減弱儲熱水箱內(nèi)壁的自然對流換熱強度。綜上可知,改變θ,對儲熱水箱內(nèi)水的流速分布情況影響較大,進而影響儲熱水箱內(nèi)的對流換熱情況。

    4.3 儲熱水箱側(cè)壁傾角對換熱強度的影響

    儲熱水箱的頂部、底部和側(cè)壁的平均Nu數(shù)隨θ的變化情況如圖10所示。

    圖10 儲熱水箱的頂部、底部和側(cè)壁的平均Nu數(shù)隨θ的變化情況Fig.10 Variation of the average Nu number of the top,bottom and side walls of the hot water storage tank with θ

    由圖10可知,儲熱水箱頂部的平均Nu數(shù)隨θ變化不大,說明θ的改變對儲熱水箱頂部的自然對流換熱強度影響較?。划?dāng)θ由15°逐漸增大至60°時,儲熱水箱側(cè)壁的平均Nu數(shù)逐漸增大,當(dāng)θ由60°逐漸增大至90°時,儲熱水箱側(cè)壁的平均Nu數(shù)逐漸減小。由上文可知,θ越小,儲熱水箱側(cè)壁自然對流換熱強度越小,則儲熱水箱側(cè)壁的平均Nu數(shù)越小,當(dāng)其他條件相同時,儲熱水箱因自然對流產(chǎn)生的熱損失就越小。隨著θ逐漸增大,儲熱水箱側(cè)壁的自然對流強度逐漸增強,導(dǎo)致儲熱水箱側(cè)壁的平均Nu數(shù)逐漸增大,儲熱水箱的散熱量也隨之增加;當(dāng)θ大于60°時,熱水箱側(cè)壁的平均Nu數(shù)有所減小,這是因為在有限空間的自然對流中,當(dāng)流體流動到水箱底部時,流動方向改變,且θ越大,流向改變得越明顯,水箱底面對流動的堵塞作用越大,流體的流速隨之減小,導(dǎo)致水箱底面的自然對流換熱強度變?nèi)酢?/p>

    儲熱水箱邊界面的平均Nu數(shù)隨θ的變化情況如圖11所示。

    將水箱側(cè)壁、底面、頂面的Nu數(shù)求平均值,得到儲熱水箱整體的平均Nu數(shù)。由圖11可知,當(dāng)θ由0°逐漸增大至60°時,儲熱水箱邊界面的平均Nu數(shù)逐漸增大,這是由于隨著θ逐漸增大,水箱側(cè)壁處的自然對流換熱強度逐漸增強,導(dǎo)致儲熱水箱邊界面的平均Nu數(shù)隨之增大;當(dāng)θ由60°逐漸增大至90°時,由于水箱底面對流動的阻礙能力越發(fā)明顯,因此水箱底面的自然對流換熱強度逐漸減小,導(dǎo)致儲熱水箱邊界面的平均Nu數(shù)隨之減小。

    圖11 儲熱水箱邊界面的平均Nu數(shù)隨θ的變化情況Fig11 Variation of the average Nu number of the boundary surface of the hot water storage tank with θ

    4.4 側(cè)壁傾角對儲熱水箱散熱過程的影響

    不同的θ條件下,儲熱水箱平均溫度隨時間的變化情況如圖12所示。

    圖12 不同的θ條件下,儲熱水箱平均溫度隨時間的變化情況Fig.12 Changes of the average temperature of the water tank with time for different θ

    由圖12可知,當(dāng)θ為75°時,儲熱水箱的散熱速度最快;當(dāng)θ為45°時,儲熱水箱的散熱速度最慢。由此可以看出,不同側(cè)壁傾角條件下,儲熱水箱的散熱量與其Nu數(shù)的變化規(guī)律不完全相符,即散熱量不僅與換熱強弱有關(guān),可能還與其比表面積相關(guān)。

    儲熱水箱的比表面積為其表面積與體積之間的比值。儲熱水箱比表面積S的計算式為

    式中:S′為儲熱水箱的表面積;V為儲熱水箱的體積。

    儲熱水箱比表面積S隨θ的變化情況如圖13所示。

    圖13 儲熱水箱比表面積S隨θ的變化情況Fig.13 The specific surface area S of the tank changes with θ

    由圖13可知,隨著θ逐漸增大,儲熱水箱比表面積S呈現(xiàn)出先減小再增大的變化趨勢。當(dāng)θ為60°時,儲熱水箱的比表面積達到最小值;當(dāng)θ為15°時,儲熱水箱的比表面積達到最大值,這與儲熱水箱各面及整體的平均Nu數(shù)的變化趨勢恰好相反。

    綜上可知,儲熱水箱的比表面和對流換熱強度對自身散熱情況的影響結(jié)果:當(dāng)θ為15°時,雖然儲熱水箱內(nèi)的對流換熱強度較弱,但比表面積較大,導(dǎo)致儲熱水箱單位時間的散熱量較大,平均溫度隨時間的下降速率較大,可見,此時比表面積對儲熱水箱的散熱量影響較大,40 min時,儲熱水箱的平均溫度降至367.1 K;當(dāng)θ為45°時,儲熱水箱的比表面積顯著減小,并接近于最小值,并且其對流換熱強度相對較弱,因此,儲熱水箱溫度下降得比較緩慢,當(dāng)τ為40 min時,儲熱水箱平均溫度為367.4 K。;當(dāng)θ為60°時,雖然儲熱水箱的比表面積較小,但是其對流換熱強度較強,所以儲熱水箱溫度下降得比較迅速,此時對流換熱強度對儲熱水箱散熱量的影響較大,當(dāng)τ為40 min時,儲熱水箱的平均溫度降低至367 K??傮w來看,在對流換熱強度和比表面積的雙重作用下,當(dāng)θ為45°時,儲熱水箱的溫度下降得最慢,散熱量也最少。

    4.5 比表面積對散熱的影響

    基于不同θ下,儲熱水箱溫度的變化情況可以看出,比表面積對儲熱水箱的散熱強弱有一定的影響。下面進一步分析儲熱水箱的比表面積與其平均溫度之間的關(guān)系。

    本文分別以邊長為 10,12,15,20,25,30 m 的正方體水箱(θ為90°)為研究對象,對應(yīng)的比表面積分別為 0.6,0.5,0.4,0.3,0.24,0.2 m-1。水的初始溫度T0為368 K,散熱時間為20 min。

    儲熱水箱的比表面積與平均溫度的關(guān)系如圖14所示。

    圖14 儲熱水箱的比表面積與平均溫度的關(guān)系Fig.14 Changes of the average temperature of the water tank with time for the cubic shape tank

    由圖14可知,比表面積越大,儲熱水箱的平均溫度越低,這是由于比表面積越大,單位體積儲熱水箱與外界之間的接觸面積越大,即換熱面積越大,因此儲熱水箱的換熱速率越大,最終導(dǎo)致儲熱水箱的平均溫度越低。

    5 結(jié)論

    本文基于數(shù)值模擬結(jié)果分析了埋于土壤中的儲熱水箱的側(cè)壁傾角θ對其內(nèi)部自然對流換熱強度的影響,分析結(jié)論如下。

    ①儲熱水箱側(cè)壁傾角越小,側(cè)壁對水箱內(nèi)部流體流動的阻礙能力越大,流體流速越小,從而影響水箱內(nèi)部的自然對流換熱強度:當(dāng)儲熱水箱側(cè)壁傾角由15°逐漸增大至60°時,該水箱的自然對流換熱強度逐漸增強,當(dāng)儲熱水箱側(cè)壁傾角由60°逐漸增大至90°時,由于水箱底面對流動的阻礙作用愈發(fā)顯著,使得該水箱內(nèi)的對流換熱強度開始緩慢減弱。

    ②當(dāng)儲熱水箱側(cè)壁傾角大于60°時,水箱底面對流體流到的阻礙作用較大,從而導(dǎo)致水箱內(nèi)的換熱強度減弱。

    ③隨著比表面積逐漸大,儲熱水箱的散熱強度逐漸增大。

    ④由自然對流換熱強度和比表面積對儲熱水箱散熱量的影響可知,當(dāng)儲熱水箱側(cè)壁傾角為45°時,其內(nèi)部水溫的下降速率較慢,散熱量較小。

    ⑤儲熱水箱側(cè)壁傾角較大時,水箱底部的溫度分層比較明顯,需要分層取水時,應(yīng)使得儲熱水箱的側(cè)壁傾角較大。

    猜你喜歡
    側(cè)壁對流壁面
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    二維有限長度柔性壁面上T-S波演化的數(shù)值研究
    高速公路隧道洞口淺埋段施工中的雙側(cè)壁導(dǎo)坑法
    柳鋼6號高爐應(yīng)對爐缸側(cè)壁溫度超高的護爐實踐
    四川冶金(2019年4期)2019-11-18 09:30:50
    越野車側(cè)壁加強筋板工藝改進
    模具制造(2019年3期)2019-06-06 02:10:52
    壁面溫度對微型內(nèi)燃機燃燒特性的影響
    基于ANSYS的自然對流換熱系數(shù)計算方法研究
    層合板上層建筑側(cè)壁抗空爆性能研究
    二元驅(qū)油水界面Marangoni對流啟動殘余油機理
    顆?!诿媾鲎步Ec數(shù)據(jù)處理
    日韩一区二区视频免费看| 国产伦一二天堂av在线观看| 国产老妇女一区| 日韩欧美精品v在线| 亚洲国产欧美人成| 熟女人妻精品中文字幕| 男女那种视频在线观看| 日本黄色视频三级网站网址| 久久精品国产自在天天线| 亚洲最大成人av| 99久久精品国产国产毛片| 亚洲不卡免费看| 两性午夜刺激爽爽歪歪视频在线观看| 村上凉子中文字幕在线| 亚洲av一区综合| 99精品久久久久人妻精品| 黄色丝袜av网址大全| 国产精品国产高清国产av| 无遮挡黄片免费观看| 淫秽高清视频在线观看| 亚洲成人久久爱视频| 国产色爽女视频免费观看| 欧美日韩综合久久久久久 | 我要搜黄色片| 精品日产1卡2卡| 国产精品99久久久久久久久| 人人妻人人澡欧美一区二区| 在线观看美女被高潮喷水网站| 午夜激情欧美在线| 久久久久免费精品人妻一区二区| 亚洲欧美日韩无卡精品| 欧美zozozo另类| 一级av片app| 搡老熟女国产l中国老女人| 亚洲 国产 在线| 精品久久国产蜜桃| 能在线免费观看的黄片| 国产精品人妻久久久久久| 久久久久久九九精品二区国产| 亚洲精品色激情综合| 最后的刺客免费高清国语| 亚洲精品乱码久久久v下载方式| 亚洲最大成人中文| avwww免费| 成人午夜高清在线视频| 国产精品伦人一区二区| 亚洲av五月六月丁香网| 久久精品国产清高在天天线| 好男人在线观看高清免费视频| 观看免费一级毛片| 999久久久精品免费观看国产| 色吧在线观看| 麻豆成人午夜福利视频| 免费av观看视频| 国产三级中文精品| 嫩草影院精品99| 亚洲人成伊人成综合网2020| 少妇的逼水好多| 91午夜精品亚洲一区二区三区 | 一个人看视频在线观看www免费| 色在线成人网| а√天堂www在线а√下载| 国产视频内射| 天天一区二区日本电影三级| 国产精品久久久久久精品电影| 久久午夜福利片| 亚洲图色成人| 欧美日韩黄片免| 国产男人的电影天堂91| 亚洲欧美日韩高清专用| 国产亚洲精品av在线| 精品人妻1区二区| 亚洲狠狠婷婷综合久久图片| 国产精品伦人一区二区| 免费av不卡在线播放| 三级男女做爰猛烈吃奶摸视频| 久久久久久九九精品二区国产| 男女边吃奶边做爰视频| 免费看av在线观看网站| 一个人看的www免费观看视频| 色综合亚洲欧美另类图片| 99riav亚洲国产免费| 精品久久久久久成人av| 国产高清不卡午夜福利| 国产91精品成人一区二区三区| 亚洲va在线va天堂va国产| 观看免费一级毛片| 两个人的视频大全免费| 深爱激情五月婷婷| 黄色丝袜av网址大全| 欧美人与善性xxx| 国产精品乱码一区二三区的特点| 久久久久免费精品人妻一区二区| 国产精品不卡视频一区二区| 嫩草影视91久久| or卡值多少钱| 别揉我奶头 嗯啊视频| 一级黄色大片毛片| 国国产精品蜜臀av免费| 国产精品久久久久久久久免| 日本成人三级电影网站| 十八禁网站免费在线| 亚洲成人久久爱视频| 欧美高清成人免费视频www| 国产av在哪里看| 免费观看人在逋| 国产精品久久久久久久久免| 国产极品精品免费视频能看的| 免费高清视频大片| 久久精品国产99精品国产亚洲性色| 亚洲va日本ⅴa欧美va伊人久久| 无遮挡黄片免费观看| 日韩av在线大香蕉| 免费人成在线观看视频色| 午夜福利高清视频| 欧美zozozo另类| 国产 一区 欧美 日韩| 麻豆久久精品国产亚洲av| 欧美性猛交黑人性爽| 色视频www国产| 熟妇人妻久久中文字幕3abv| 99久久中文字幕三级久久日本| 男人和女人高潮做爰伦理| 99久久精品一区二区三区| 亚洲av中文av极速乱 | 伦理电影大哥的女人| 99久久中文字幕三级久久日本| 国产真实乱freesex| 在线观看舔阴道视频| 精华霜和精华液先用哪个| 夜夜爽天天搞| 精品久久久久久久久久免费视频| 成人国产一区最新在线观看| 少妇裸体淫交视频免费看高清| 一卡2卡三卡四卡精品乱码亚洲| 在线a可以看的网站| 亚洲国产欧洲综合997久久,| 亚洲成a人片在线一区二区| 免费看光身美女| 草草在线视频免费看| 99久久精品国产国产毛片| 赤兔流量卡办理| 精品人妻视频免费看| 欧美成人a在线观看| 亚洲在线自拍视频| eeuss影院久久| 午夜老司机福利剧场| 精品久久久久久久久久免费视频| 亚洲欧美日韩无卡精品| 最近在线观看免费完整版| 国产精品人妻久久久久久| 黄片wwwwww| a在线观看视频网站| 窝窝影院91人妻| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 欧美日韩国产亚洲二区| avwww免费| 久久午夜亚洲精品久久| 久久久久免费精品人妻一区二区| 国产精品av视频在线免费观看| 成年人黄色毛片网站| 在线播放无遮挡| 3wmmmm亚洲av在线观看| 日韩亚洲欧美综合| 婷婷亚洲欧美| 女生性感内裤真人,穿戴方法视频| 一级av片app| 在线天堂最新版资源| 久久久久久久午夜电影| 国产久久久一区二区三区| 精品一区二区三区人妻视频| 一区二区三区四区激情视频 | 久久久国产成人精品二区| 三级男女做爰猛烈吃奶摸视频| 亚洲精品粉嫩美女一区| 中文字幕免费在线视频6| 熟女电影av网| 免费观看人在逋| 中文在线观看免费www的网站| 22中文网久久字幕| 久久午夜福利片| 亚洲中文日韩欧美视频| 国产蜜桃级精品一区二区三区| 亚洲四区av| 亚洲国产日韩欧美精品在线观看| 深夜a级毛片| 91av网一区二区| 免费无遮挡裸体视频| 欧美日本亚洲视频在线播放| 亚洲最大成人av| 熟女人妻精品中文字幕| 国产av麻豆久久久久久久| 自拍偷自拍亚洲精品老妇| 日本a在线网址| 免费黄网站久久成人精品| 在线观看66精品国产| 亚洲综合色惰| 日韩高清综合在线| 高清在线国产一区| 少妇高潮的动态图| 日韩欧美国产一区二区入口| 亚洲最大成人手机在线| 亚洲av五月六月丁香网| 久久久精品欧美日韩精品| 精品久久久久久久久久久久久| 中文字幕av成人在线电影| 欧美在线一区亚洲| 91在线观看av| 床上黄色一级片| 亚洲av电影不卡..在线观看| 欧美又色又爽又黄视频| 给我免费播放毛片高清在线观看| 亚洲色图av天堂| 五月玫瑰六月丁香| 久久精品国产鲁丝片午夜精品 | 国产欧美日韩一区二区精品| 国产一区二区激情短视频| 午夜激情欧美在线| 99视频精品全部免费 在线| av在线天堂中文字幕| 国产高清三级在线| 久久久成人免费电影| 亚洲欧美激情综合另类| 国产女主播在线喷水免费视频网站 | 超碰av人人做人人爽久久| 国产综合懂色| 成人国产麻豆网| 1024手机看黄色片| 国产老妇女一区| 免费看美女性在线毛片视频| 99久久精品一区二区三区| 国产真实伦视频高清在线观看 | 蜜桃亚洲精品一区二区三区| 午夜亚洲福利在线播放| 97热精品久久久久久| 国产单亲对白刺激| 韩国av一区二区三区四区| 欧美又色又爽又黄视频| 国产精品女同一区二区软件 | 亚洲人成网站在线播放欧美日韩| 日韩欧美三级三区| 亚洲欧美日韩高清专用| 99精品在免费线老司机午夜| 老熟妇仑乱视频hdxx| 国产精品日韩av在线免费观看| 国产中年淑女户外野战色| 国产伦一二天堂av在线观看| 最新中文字幕久久久久| 欧美成人一区二区免费高清观看| 久久国产精品人妻蜜桃| 99国产极品粉嫩在线观看| 午夜福利视频1000在线观看| 看十八女毛片水多多多| 午夜亚洲福利在线播放| 国产亚洲精品综合一区在线观看| 国产av不卡久久| 日本在线视频免费播放| 男女下面进入的视频免费午夜| x7x7x7水蜜桃| 中文字幕av成人在线电影| 高清毛片免费观看视频网站| av天堂中文字幕网| 听说在线观看完整版免费高清| 精品久久久噜噜| 国产一区二区三区av在线 | 欧美成人一区二区免费高清观看| 少妇高潮的动态图| 最近视频中文字幕2019在线8| 成人美女网站在线观看视频| 在线免费观看的www视频| 99久国产av精品| 大型黄色视频在线免费观看| 在线免费十八禁| 波野结衣二区三区在线| av视频在线观看入口| 成人精品一区二区免费| 亚洲熟妇熟女久久| 男女那种视频在线观看| 美女黄网站色视频| 免费黄网站久久成人精品| 2021天堂中文幕一二区在线观| 一本一本综合久久| 国产午夜福利久久久久久| 性插视频无遮挡在线免费观看| 欧美+日韩+精品| 亚洲内射少妇av| 午夜久久久久精精品| 日本色播在线视频| 欧美中文日本在线观看视频| 村上凉子中文字幕在线| 日本黄色片子视频| 午夜影院日韩av| bbb黄色大片| 亚洲电影在线观看av| 亚洲五月天丁香| 欧美精品啪啪一区二区三区| 免费观看人在逋| 色综合婷婷激情| 99在线人妻在线中文字幕| 黄色丝袜av网址大全| 精品人妻视频免费看| x7x7x7水蜜桃| 一进一出抽搐gif免费好疼| 久久久久九九精品影院| 国语自产精品视频在线第100页| av在线天堂中文字幕| 波多野结衣巨乳人妻| 欧美黑人欧美精品刺激| 高清在线国产一区| 一进一出好大好爽视频| 一本精品99久久精品77| 高清日韩中文字幕在线| 久久精品国产亚洲av天美| 日韩中字成人| 一级av片app| av天堂在线播放| 日韩 亚洲 欧美在线| 1000部很黄的大片| 国产亚洲av嫩草精品影院| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩高清在线视频| 亚洲欧美精品综合久久99| av在线亚洲专区| 村上凉子中文字幕在线| 麻豆成人午夜福利视频| 99久久久亚洲精品蜜臀av| 午夜免费激情av| 欧美+日韩+精品| 久久久久久大精品| 中文字幕高清在线视频| 亚洲无线在线观看| 久久久久久久久中文| 哪里可以看免费的av片| 国产激情偷乱视频一区二区| 成年女人永久免费观看视频| 国产伦精品一区二区三区四那| 日韩欧美三级三区| 国产色婷婷99| АⅤ资源中文在线天堂| a级一级毛片免费在线观看| 波多野结衣高清作品| 高清在线国产一区| 丰满人妻一区二区三区视频av| 国产久久久一区二区三区| avwww免费| 成人亚洲精品av一区二区| 91av网一区二区| 免费看日本二区| 日本与韩国留学比较| 一级黄片播放器| 日韩精品青青久久久久久| 51国产日韩欧美| 一进一出抽搐gif免费好疼| 日韩亚洲欧美综合| 在线观看av片永久免费下载| 精品无人区乱码1区二区| 精品人妻偷拍中文字幕| 久久久久性生活片| 国产 一区精品| 成人二区视频| 婷婷色综合大香蕉| 国产久久久一区二区三区| 久久这里只有精品中国| 精品人妻一区二区三区麻豆 | 少妇被粗大猛烈的视频| 别揉我奶头 嗯啊视频| 搡女人真爽免费视频火全软件 | 日日摸夜夜添夜夜添小说| 免费高清视频大片| 少妇的逼水好多| or卡值多少钱| 一个人看视频在线观看www免费| 一个人免费在线观看电影| 欧美人与善性xxx| 校园人妻丝袜中文字幕| 国产白丝娇喘喷水9色精品| av视频在线观看入口| 国产毛片a区久久久久| 少妇猛男粗大的猛烈进出视频 | 非洲黑人性xxxx精品又粗又长| 免费在线观看影片大全网站| 在线免费十八禁| 亚洲专区中文字幕在线| 在线看三级毛片| 国产精品福利在线免费观看| 成年女人毛片免费观看观看9| 国产精品无大码| 国产v大片淫在线免费观看| 国产一区二区在线观看日韩| 又爽又黄a免费视频| 久久久久久伊人网av| 午夜久久久久精精品| 亚洲av免费在线观看| 国产在线精品亚洲第一网站| 国产一区二区亚洲精品在线观看| 久久久久免费精品人妻一区二区| 亚洲欧美日韩卡通动漫| 成人国产综合亚洲| 欧美成人性av电影在线观看| 国产精品电影一区二区三区| 午夜a级毛片| 国产三级在线视频| 欧美又色又爽又黄视频| 亚洲精品粉嫩美女一区| 国产精品久久电影中文字幕| 成人鲁丝片一二三区免费| .国产精品久久| 偷拍熟女少妇极品色| 啦啦啦韩国在线观看视频| 日本 av在线| 欧美3d第一页| 两个人的视频大全免费| 在线观看午夜福利视频| 日韩国内少妇激情av| 亚洲精品国产成人久久av| 国产色爽女视频免费观看| 久久精品国产清高在天天线| 欧美极品一区二区三区四区| 嫩草影视91久久| 美女 人体艺术 gogo| 成人国产一区最新在线观看| 国产黄a三级三级三级人| 国产午夜福利久久久久久| 春色校园在线视频观看| 永久网站在线| 国产成人av教育| 亚洲成人精品中文字幕电影| 免费黄网站久久成人精品| 国产私拍福利视频在线观看| 国产欧美日韩精品一区二区| 综合色av麻豆| 不卡视频在线观看欧美| 日本 欧美在线| 亚洲中文字幕日韩| 乱系列少妇在线播放| 亚洲av中文av极速乱 | 日韩精品青青久久久久久| 欧美3d第一页| 国产精品久久久久久精品电影| 淫妇啪啪啪对白视频| 色哟哟·www| 国产成人一区二区在线| 日本 欧美在线| 亚洲中文字幕一区二区三区有码在线看| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| 99久久成人亚洲精品观看| 亚洲熟妇中文字幕五十中出| 在线播放无遮挡| av视频在线观看入口| 九九在线视频观看精品| 女人被狂操c到高潮| av中文乱码字幕在线| 亚洲av不卡在线观看| 国产私拍福利视频在线观看| 国产在线男女| а√天堂www在线а√下载| 日本成人三级电影网站| 日韩 亚洲 欧美在线| 深夜a级毛片| 久久久久精品国产欧美久久久| АⅤ资源中文在线天堂| 欧美人与善性xxx| 国产精品久久视频播放| 久久久久久九九精品二区国产| 国产69精品久久久久777片| 男女边吃奶边做爰视频| 五月玫瑰六月丁香| 成熟少妇高潮喷水视频| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 久久久久久九九精品二区国产| 给我免费播放毛片高清在线观看| 精品不卡国产一区二区三区| 毛片一级片免费看久久久久 | 在线观看av片永久免费下载| 麻豆国产av国片精品| 全区人妻精品视频| www.www免费av| 在线观看66精品国产| 中亚洲国语对白在线视频| 婷婷亚洲欧美| 丝袜美腿在线中文| 又黄又爽又免费观看的视频| а√天堂www在线а√下载| 日韩国内少妇激情av| 嫩草影院新地址| 婷婷亚洲欧美| 国产淫片久久久久久久久| 国产午夜精品久久久久久一区二区三区 | 一区二区三区免费毛片| 91久久精品电影网| 老司机福利观看| 国产精品一区二区三区四区久久| 中文资源天堂在线| 一区二区三区免费毛片| 国产三级中文精品| 日本在线视频免费播放| 欧美性猛交黑人性爽| 白带黄色成豆腐渣| 久99久视频精品免费| 91av网一区二区| 国产在线精品亚洲第一网站| 偷拍熟女少妇极品色| 国产视频一区二区在线看| 日韩 亚洲 欧美在线| 久久精品久久久久久噜噜老黄 | 制服丝袜大香蕉在线| 一级黄色大片毛片| 韩国av在线不卡| 99久久精品热视频| 日韩欧美在线乱码| 97热精品久久久久久| 97碰自拍视频| 高清在线国产一区| 嫩草影院精品99| 欧美最新免费一区二区三区| 非洲黑人性xxxx精品又粗又长| 日本 av在线| 日韩一区二区视频免费看| 亚洲欧美日韩卡通动漫| 国产高潮美女av| 成人亚洲精品av一区二区| 长腿黑丝高跟| 国产精品一区二区三区四区免费观看 | 亚洲人成伊人成综合网2020| 日韩欧美一区二区三区在线观看| 在线国产一区二区在线| 日韩欧美精品免费久久| 女同久久另类99精品国产91| 18禁裸乳无遮挡免费网站照片| 真人一进一出gif抽搐免费| 亚洲最大成人中文| 夜夜爽天天搞| 搡老熟女国产l中国老女人| 波多野结衣高清作品| 欧美最黄视频在线播放免费| 亚洲精品色激情综合| 极品教师在线视频| 日韩精品青青久久久久久| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av在线| 很黄的视频免费| 99国产极品粉嫩在线观看| 一级av片app| 日本 av在线| 九九爱精品视频在线观看| 国产私拍福利视频在线观看| 亚洲欧美日韩高清在线视频| 亚洲精品国产成人久久av| 久久久久免费精品人妻一区二区| 亚洲av第一区精品v没综合| 日韩,欧美,国产一区二区三区 | 男人的好看免费观看在线视频| 在线a可以看的网站| 亚洲精品久久国产高清桃花| 国产亚洲精品久久久久久毛片| 日韩高清综合在线| 国产亚洲精品久久久久久毛片| 国产精华一区二区三区| av在线蜜桃| x7x7x7水蜜桃| 91av网一区二区| 亚洲美女搞黄在线观看 | 床上黄色一级片| 日韩在线高清观看一区二区三区 | 亚洲avbb在线观看| 久久久久久伊人网av| 黄色欧美视频在线观看| 欧美一区二区国产精品久久精品| 搡老熟女国产l中国老女人| 欧美成人性av电影在线观看| 亚洲经典国产精华液单| 男插女下体视频免费在线播放| 亚洲成人久久性| 99热6这里只有精品| 亚洲欧美日韩无卡精品| 香蕉av资源在线| 欧洲精品卡2卡3卡4卡5卡区| 日本 av在线| 黄色丝袜av网址大全| 一进一出好大好爽视频| 十八禁国产超污无遮挡网站| 亚洲精品粉嫩美女一区| 亚洲国产日韩欧美精品在线观看| 岛国在线免费视频观看| 国产精品98久久久久久宅男小说| 桃红色精品国产亚洲av| 22中文网久久字幕| 97超级碰碰碰精品色视频在线观看| 亚洲精品成人久久久久久| 亚洲第一电影网av| 国产精品一及| a级一级毛片免费在线观看| 一级毛片久久久久久久久女| 伦理电影大哥的女人| 毛片女人毛片| 久久欧美精品欧美久久欧美| 女的被弄到高潮叫床怎么办 | x7x7x7水蜜桃| 午夜福利欧美成人| 观看免费一级毛片| 成年女人永久免费观看视频| 91狼人影院| 好男人在线观看高清免费视频| 久久这里只有精品中国| 伦精品一区二区三区| 亚洲成人久久爱视频| 丰满人妻一区二区三区视频av| 国产精品国产高清国产av| 国产视频内射| 麻豆国产97在线/欧美| 欧美日韩综合久久久久久 | 精品人妻1区二区|