• <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ù)處理
    国产精品无大码| 老熟妇乱子伦视频在线观看| 99热这里只有是精品在线观看| 五月玫瑰六月丁香| 午夜福利在线在线| 国产成人一区二区在线| 国产一区二区亚洲精品在线观看| 高清在线国产一区| 久久久久性生活片| 亚洲四区av| 日本欧美国产在线视频| 美女cb高潮喷水在线观看| 久久亚洲真实| 嫩草影院入口| 一区二区三区免费毛片| av视频在线观看入口| 91在线精品国自产拍蜜月| а√天堂www在线а√下载| 男女之事视频高清在线观看| 欧美成人性av电影在线观看| 亚洲aⅴ乱码一区二区在线播放| 男女之事视频高清在线观看| 91久久精品国产一区二区成人| 一夜夜www| 午夜福利在线观看吧| 久久午夜福利片| 一本一本综合久久| 久久午夜亚洲精品久久| 可以在线观看的亚洲视频| 男人和女人高潮做爰伦理| 国产欧美日韩一区二区精品| 国产人妻一区二区三区在| 欧美xxxx性猛交bbbb| 午夜激情福利司机影院| 非洲黑人性xxxx精品又粗又长| 久久精品国产亚洲网站| 久久久国产成人精品二区| 如何舔出高潮| 特大巨黑吊av在线直播| 亚洲美女搞黄在线观看 | 亚洲精品影视一区二区三区av| 免费搜索国产男女视频| 一边摸一边抽搐一进一小说| 俺也久久电影网| 久久九九热精品免费| 简卡轻食公司| 色视频www国产| 国产真实乱freesex| 国产人妻一区二区三区在| 日本撒尿小便嘘嘘汇集6| 在线免费观看不下载黄p国产 | 一个人看的www免费观看视频| 可以在线观看毛片的网站| 亚洲在线自拍视频| 精品久久久久久久久久免费视频| 美女被艹到高潮喷水动态| 可以在线观看毛片的网站| 国产激情偷乱视频一区二区| 91狼人影院| 白带黄色成豆腐渣| 国产老妇女一区| 国产午夜福利久久久久久| 一进一出抽搐gif免费好疼| 精品人妻1区二区| 婷婷色综合大香蕉| 国产伦精品一区二区三区视频9| 深夜精品福利| 成人影院久久| 久久久久久久精品精品| 欧美成人一区二区免费高清观看| 国产成人aa在线观看| 九九在线视频观看精品| 久久女婷五月综合色啪小说| 欧美区成人在线视频| 欧美精品一区二区大全| 男女啪啪激烈高潮av片| freevideosex欧美| 国产一区二区三区av在线| 少妇丰满av| 亚洲国产精品专区欧美| xxx大片免费视频| 男女国产视频网站| 欧美一区二区亚洲| 久久国产亚洲av麻豆专区| 毛片女人毛片| 人妻系列 视频| 99热这里只有精品一区| 亚洲美女搞黄在线观看| 男女无遮挡免费网站观看| 91精品一卡2卡3卡4卡| 国产伦在线观看视频一区| 亚洲中文av在线| 日本-黄色视频高清免费观看| 2022亚洲国产成人精品| 人妻夜夜爽99麻豆av| 亚洲色图av天堂| 亚洲欧美日韩东京热| 青春草视频在线免费观看| 赤兔流量卡办理| 久久久久久久久久成人| 国产男女超爽视频在线观看| 色网站视频免费| 日韩av不卡免费在线播放| 一区在线观看完整版| 国产精品国产三级国产av玫瑰| 一级毛片电影观看| 99久久综合免费| 少妇人妻精品综合一区二区| 久久久久久久久久久丰满| 18+在线观看网站| 波野结衣二区三区在线| 久久久亚洲精品成人影院| 国产精品秋霞免费鲁丝片| 久久ye,这里只有精品| av国产精品久久久久影院| 精品人妻熟女av久视频| 特大巨黑吊av在线直播| a级一级毛片免费在线观看| .国产精品久久| 精品人妻偷拍中文字幕| 久久99精品国语久久久| 亚州av有码| 久久6这里有精品| 国产黄色视频一区二区在线观看| 色视频www国产| 黄色欧美视频在线观看| 一级毛片黄色毛片免费观看视频| 国产精品无大码| 啦啦啦在线观看免费高清www| av国产精品久久久久影院| 国产人妻一区二区三区在| 国产日韩欧美亚洲二区| 91精品一卡2卡3卡4卡| 两个人的视频大全免费| 亚洲怡红院男人天堂| 美女福利国产在线 | 成人综合一区亚洲| 女性生殖器流出的白浆| 自拍偷自拍亚洲精品老妇| 中文精品一卡2卡3卡4更新| 久久 成人 亚洲| 国产淫片久久久久久久久| 国产91av在线免费观看| 欧美+日韩+精品| 午夜福利网站1000一区二区三区| 哪个播放器可以免费观看大片| 多毛熟女@视频| av不卡在线播放| 国产91av在线免费观看| 在线精品无人区一区二区三 | videossex国产| 久久精品国产亚洲网站| 日本欧美视频一区| 黄色配什么色好看| 婷婷色麻豆天堂久久| 我的女老师完整版在线观看| 天堂8中文在线网| 水蜜桃什么品种好| av黄色大香蕉| 日本黄色日本黄色录像| 亚洲精品日韩在线中文字幕| av又黄又爽大尺度在线免费看| 女人十人毛片免费观看3o分钟| 国产黄色视频一区二区在线观看| 蜜桃久久精品国产亚洲av| 日韩一区二区视频免费看| 国产白丝娇喘喷水9色精品| 亚洲电影在线观看av| 欧美日韩在线观看h| 大码成人一级视频| 亚洲婷婷狠狠爱综合网| 美女高潮的动态| 免费人成在线观看视频色| 成人亚洲精品一区在线观看 | 国产一区亚洲一区在线观看| 久久97久久精品| 三级国产精品片| 国产一区有黄有色的免费视频| 国产在线视频一区二区| av网站免费在线观看视频| 久久久久视频综合| 大片免费播放器 马上看| 亚洲精品成人av观看孕妇| 国产高清三级在线| 人人妻人人澡人人爽人人夜夜| 国内少妇人妻偷人精品xxx网站| 国产成人免费观看mmmm| 男男h啪啪无遮挡| 久久人人爽人人片av| 国产精品麻豆人妻色哟哟久久| 麻豆成人av视频| 亚洲激情五月婷婷啪啪| av免费在线看不卡| 免费观看性生交大片5| 人妻系列 视频| 中文字幕久久专区| 亚洲经典国产精华液单| 精品人妻视频免费看| 亚洲精品国产成人久久av| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 亚洲在久久综合| 亚洲精品自拍成人| 久久久久久久久大av| 一二三四中文在线观看免费高清| 少妇人妻精品综合一区二区| 国产 一区 欧美 日韩| 国产在线免费精品| 久久99精品国语久久久| 多毛熟女@视频| 亚洲va在线va天堂va国产| 中文字幕亚洲精品专区| 国产av一区二区精品久久 | 日韩中字成人| 国产亚洲午夜精品一区二区久久| 性色avwww在线观看| 晚上一个人看的免费电影| 欧美高清成人免费视频www| 亚洲欧美一区二区三区黑人 | 日韩视频在线欧美| 国产白丝娇喘喷水9色精品| 国产精品偷伦视频观看了| 久久久精品94久久精品| 我要看黄色一级片免费的| 国产爽快片一区二区三区| 国产探花极品一区二区| 久久久亚洲精品成人影院| 国产男女内射视频| 精品一区二区三区视频在线| av在线老鸭窝| 一个人看的www免费观看视频| 视频区图区小说| 超碰av人人做人人爽久久| 欧美精品亚洲一区二区| a级一级毛片免费在线观看| 国产女主播在线喷水免费视频网站| 日韩大片免费观看网站| 久久午夜福利片| 欧美成人午夜免费资源| 女性生殖器流出的白浆| 日韩中字成人| 国产精品精品国产色婷婷| 日韩一区二区三区影片| 少妇的逼水好多| 亚洲精品日韩在线中文字幕| 国产亚洲午夜精品一区二区久久| 免费观看的影片在线观看| 欧美成人午夜免费资源| 国产成人freesex在线| 美女主播在线视频| 超碰av人人做人人爽久久| 啦啦啦啦在线视频资源| 亚洲精华国产精华液的使用体验| 日韩av在线免费看完整版不卡| 久久ye,这里只有精品| 性色av一级| 国产成人freesex在线| av免费观看日本| av在线播放精品| 女的被弄到高潮叫床怎么办| 久久久久久久久久久免费av| 在线看a的网站| 亚洲国产av新网站| 国产免费视频播放在线视频| 中文字幕免费在线视频6| 免费大片黄手机在线观看| 又大又黄又爽视频免费| 久久久久久久精品精品| 夜夜看夜夜爽夜夜摸| 亚洲国产精品专区欧美| 天天躁日日操中文字幕| 国产精品蜜桃在线观看| 久久精品国产亚洲av涩爱| 3wmmmm亚洲av在线观看| 91精品伊人久久大香线蕉| 国产乱人偷精品视频| 国产熟女欧美一区二区| av又黄又爽大尺度在线免费看| 国产精品久久久久成人av| 亚洲av成人精品一区久久| 亚洲欧美一区二区三区黑人 | 亚洲久久久国产精品| 久久精品人妻少妇| 成人特级av手机在线观看| a 毛片基地| 午夜免费男女啪啪视频观看| 在线观看一区二区三区| 美女内射精品一级片tv| 久久精品国产a三级三级三级| 国产精品一区二区在线观看99| 久久久久久久久大av| 日本一二三区视频观看| 国产亚洲av片在线观看秒播厂| 色综合色国产| 色婷婷av一区二区三区视频| 婷婷色综合www| 久久久久久久亚洲中文字幕| 永久网站在线| 日韩av免费高清视频| 乱系列少妇在线播放| 99久久精品一区二区三区| 亚洲激情五月婷婷啪啪| 在线观看美女被高潮喷水网站| 亚洲成人一二三区av| 国产白丝娇喘喷水9色精品| 久久久成人免费电影| 18禁裸乳无遮挡免费网站照片| 丝袜喷水一区| 91久久精品电影网| 精品一区在线观看国产| 三级国产精品片| 国产精品一区二区三区四区免费观看| 亚洲伊人久久精品综合| 国产在视频线精品| 成人美女网站在线观看视频| 成人毛片a级毛片在线播放| 日本猛色少妇xxxxx猛交久久| 国产精品国产三级国产专区5o| 99国产精品免费福利视频| 国产一区二区三区av在线| 国产免费视频播放在线视频| 又黄又爽又刺激的免费视频.| 纵有疾风起免费观看全集完整版| 99九九线精品视频在线观看视频| 国产精品久久久久久久久免| 免费大片黄手机在线观看| 日韩中文字幕视频在线看片 | 高清午夜精品一区二区三区| 久久久久久伊人网av| 五月天丁香电影| 91狼人影院| 国产伦精品一区二区三区视频9| 日本一二三区视频观看| 亚洲欧美成人精品一区二区| 最黄视频免费看| videossex国产| 另类亚洲欧美激情| 亚洲精品,欧美精品| 午夜福利视频精品| 久久久久久久亚洲中文字幕| 欧美精品人与动牲交sv欧美| 成人特级av手机在线观看| 亚洲成人一二三区av| 九九在线视频观看精品| 国产亚洲91精品色在线| a级毛色黄片| 久久久久久久久久久免费av| 亚洲精品色激情综合| 国产美女午夜福利| 伦精品一区二区三区| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 小蜜桃在线观看免费完整版高清| 久久99蜜桃精品久久| 久久午夜福利片| 国产色婷婷99| 亚洲精品,欧美精品| 少妇猛男粗大的猛烈进出视频| 最近最新中文字幕大全电影3| 亚洲精品自拍成人| 久久青草综合色| 丝袜喷水一区| 国产精品国产三级专区第一集| 啦啦啦视频在线资源免费观看| 老女人水多毛片| 中文字幕av成人在线电影| 免费大片18禁| 尤物成人国产欧美一区二区三区| 国产成人精品久久久久久| 美女xxoo啪啪120秒动态图| 亚洲国产精品国产精品| 在线观看av片永久免费下载| 天美传媒精品一区二区| 久久久a久久爽久久v久久| 国产高潮美女av| 91久久精品国产一区二区成人| 国产成人a区在线观看| 在线观看免费日韩欧美大片 | 国产av精品麻豆| 少妇的逼好多水| 成年免费大片在线观看| 国产高清不卡午夜福利| 伦理电影免费视频| 免费黄网站久久成人精品| 免费不卡的大黄色大毛片视频在线观看| 下体分泌物呈黄色| 国产一区二区三区综合在线观看 | 青春草亚洲视频在线观看| 亚洲美女搞黄在线观看| 亚洲av欧美aⅴ国产| 国产免费视频播放在线视频| h视频一区二区三区| 少妇人妻 视频| 国产精品熟女久久久久浪| 久热这里只有精品99| 亚洲精品日韩av片在线观看| 国产精品久久久久成人av| 高清视频免费观看一区二区| 99久久中文字幕三级久久日本| 国产乱来视频区| 在线天堂最新版资源| 一级二级三级毛片免费看| 六月丁香七月| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 男人爽女人下面视频在线观看| 岛国毛片在线播放| 精品午夜福利在线看| 亚洲精品中文字幕在线视频 | 日韩视频在线欧美| 亚洲精品日韩av片在线观看| 日本一二三区视频观看| 亚洲,欧美,日韩| 亚洲精品色激情综合| 成人二区视频| 国产爱豆传媒在线观看| av卡一久久| 另类亚洲欧美激情| 九九久久精品国产亚洲av麻豆| 欧美xxxx性猛交bbbb| 亚洲中文av在线| 激情五月婷婷亚洲| 精品亚洲成国产av| 少妇人妻一区二区三区视频| 国产亚洲精品久久久com| 亚洲在久久综合| 国产精品免费大片| 中文字幕人妻熟人妻熟丝袜美| av免费观看日本| 51国产日韩欧美| 欧美三级亚洲精品| 久久综合国产亚洲精品| 精品一品国产午夜福利视频| 最新中文字幕久久久久| 亚洲精品第二区| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品国产色婷婷电影| 免费黄频网站在线观看国产| 五月伊人婷婷丁香| 乱系列少妇在线播放| 人人妻人人澡人人爽人人夜夜| 丰满乱子伦码专区| 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 91精品国产国语对白视频| 国产精品不卡视频一区二区| 色吧在线观看| 18禁动态无遮挡网站| 日韩人妻高清精品专区| 日本一二三区视频观看| 日本黄色日本黄色录像| 亚洲欧洲日产国产| 中文字幕人妻熟人妻熟丝袜美| 国产 一区精品| 丝袜脚勾引网站| 黄色日韩在线| 午夜免费鲁丝| 国产免费又黄又爽又色| 九九久久精品国产亚洲av麻豆| 久久韩国三级中文字幕| 亚洲精品456在线播放app| 日本黄大片高清| 亚洲av福利一区| a级毛色黄片| 人妻夜夜爽99麻豆av| 国产精品三级大全| 香蕉精品网在线| 亚洲精品日韩在线中文字幕| 国产亚洲精品久久久com| .国产精品久久| av国产精品久久久久影院| 18+在线观看网站| 国产精品国产三级国产av玫瑰| 久久久久久久国产电影| 亚洲精品,欧美精品| 国产精品一区www在线观看| 全区人妻精品视频| 嫩草影院入口| 亚洲国产日韩一区二区| 舔av片在线| 午夜视频国产福利| 亚洲国产欧美在线一区| 蜜桃久久精品国产亚洲av| 亚洲欧美精品专区久久| 久久精品夜色国产| 日本黄大片高清| 最近手机中文字幕大全| 一个人看的www免费观看视频| 日韩成人av中文字幕在线观看| 国产午夜精品一二区理论片| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品乱码久久久v下载方式| 国产黄频视频在线观看| 久久精品夜色国产| 交换朋友夫妻互换小说| 婷婷色综合www| 国产av精品麻豆| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av涩爱| 噜噜噜噜噜久久久久久91| 大又大粗又爽又黄少妇毛片口| 在线观看一区二区三区| 亚洲国产色片| 久久精品熟女亚洲av麻豆精品| 交换朋友夫妻互换小说| 一本色道久久久久久精品综合| 91在线精品国自产拍蜜月| 久久久久精品性色| 国产大屁股一区二区在线视频| 国产精品一区二区在线观看99| 中文天堂在线官网| 欧美日韩国产mv在线观看视频 | 国产精品欧美亚洲77777| 联通29元200g的流量卡| 一边亲一边摸免费视频| 少妇精品久久久久久久| 丰满人妻一区二区三区视频av| av线在线观看网站| 最新中文字幕久久久久| 精品99又大又爽又粗少妇毛片| 深夜a级毛片| 色视频www国产| 午夜日本视频在线| 少妇熟女欧美另类| 这个男人来自地球电影免费观看 | 国产午夜精品一二区理论片| 麻豆国产97在线/欧美| 久热这里只有精品99| 国产亚洲午夜精品一区二区久久| 久久精品熟女亚洲av麻豆精品| 少妇的逼好多水| 中文字幕久久专区| 大码成人一级视频| 国产成人精品一,二区| 日韩大片免费观看网站| 亚洲人成网站高清观看| 大码成人一级视频| 精品一品国产午夜福利视频| 欧美另类一区| 一级爰片在线观看| 一级毛片 在线播放| 精品人妻一区二区三区麻豆| 久久99热这里只有精品18| 高清在线视频一区二区三区| 亚洲伊人久久精品综合| 亚洲av在线观看美女高潮| 草草在线视频免费看| 久久精品国产鲁丝片午夜精品| 日本vs欧美在线观看视频 | 在线观看一区二区三区激情| 99久久中文字幕三级久久日本| 老师上课跳d突然被开到最大视频| 日日撸夜夜添| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 亚洲人成网站高清观看| 精品午夜福利在线看| 国产免费又黄又爽又色| 有码 亚洲区| 80岁老熟妇乱子伦牲交| 91久久精品国产一区二区三区| 久久久欧美国产精品| 久久综合国产亚洲精品| 日本黄色日本黄色录像| 国产av国产精品国产| h日本视频在线播放| 啦啦啦视频在线资源免费观看| 久久精品国产亚洲av涩爱| 国精品久久久久久国模美| 看十八女毛片水多多多| 毛片一级片免费看久久久久| 黄色日韩在线| 日韩欧美精品免费久久| 天堂中文最新版在线下载| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 久久精品夜色国产| av播播在线观看一区| 欧美 日韩 精品 国产| 少妇人妻 视频| av黄色大香蕉| 国产欧美另类精品又又久久亚洲欧美| 亚洲第一av免费看| 国国产精品蜜臀av免费| 久久毛片免费看一区二区三区| 久久影院123| 99国产精品免费福利视频| 国产精品久久久久久av不卡| 精品国产一区二区三区久久久樱花 | 亚洲av中文av极速乱| 韩国av在线不卡| 国产午夜精品久久久久久一区二区三区| 美女视频免费永久观看网站| 在线观看一区二区三区| 中文字幕人妻熟人妻熟丝袜美| xxx大片免费视频| 中文字幕久久专区| 男女国产视频网站| 国产伦在线观看视频一区| 在线观看免费高清a一片| 成人亚洲欧美一区二区av| 日本欧美国产在线视频| 日韩一本色道免费dvd| 干丝袜人妻中文字幕| 亚洲欧洲国产日韩| 在线观看人妻少妇| 国产精品三级大全| 久久精品国产亚洲av涩爱| 天天躁日日操中文字幕| 欧美日韩一区二区视频在线观看视频在线| 国产 精品1| 日产精品乱码卡一卡2卡三| 国产v大片淫在线免费观看| 久久99热这里只有精品18| 国产色爽女视频免费观看|