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

    不同輻射模型在太陽輻射數(shù)值模擬中的比較

    2015-12-30 02:07:11嚴(yán)寒,張鴻雁
    節(jié)能技術(shù) 2015年5期
    關(guān)鍵詞:太陽輻射數(shù)值模擬

    不同輻射模型在太陽輻射數(shù)值模擬中的比較

    嚴(yán)寒1,張鴻雁2

    (1.中國新時(shí)代國際工程公司,陜西西安710054;

    2.西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院,陜西西安710055)

    摘要:為了對太陽輻射傳熱過程進(jìn)行研究,在數(shù)值模擬計(jì)算中,分別比較了離散傳播(DT)輻射模型、基于球形諧波法的P-1輻射模型、羅斯蘭德(Rosseland)輻射模型、表面(S2S)輻射模型和離散坐標(biāo)(DO)輻射模型的特點(diǎn)和適用性。研究發(fā)現(xiàn)使用DO模型模擬太陽輻射傳熱的結(jié)果與實(shí)際吻合較好。

    關(guān)鍵詞:輻射模型;太陽輻射;數(shù)值模擬;CFD;DO模型

    中圖分類號:TK124

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1002-6339 (2015) 05-0428-04

    Abstract:In order to study the solar radiation heat transfer process,in the calculation by numerical simulation,the following radiation models were tested: Discrete Transfer (DT) Radiation Model, P-1 Radiation Model based on spherical-harmonics method, Rosseland Radiation Model, Surface-to-Surface (S2S) Radiation Model and Discrete Ordinates (DO) Radiation Model. It compared the characteristic and applicability of different radiation model. It was found that DO model was in good agreement with the practice for solar radiation.

    收稿日期2014-12-18修訂稿日期2015-03-20

    作者簡介:嚴(yán)寒(1980~),男,碩士,工程師,主要從事暖通設(shè)計(jì)和數(shù)值模擬研究。

    Comparison of Different Radiation Models for Numerical Simulation of Solar Radiation

    YAN Han1,ZHANG Hong-yan2

    (1.China New Era Intemational Engineering Corporation, Xi’an 710054,China;

    2.Xi’an University of Architecture and Technology, Xi’an 710055,China)

    Key words:radiation models; solar radiation; numerical simulation; CFD;DO model

    輻射傳熱是一種非接觸式換熱過程,它作為熱量傳遞的一種重要方式,在過程的機(jī)理上與導(dǎo)熱、對流換熱是根本不同的,它靠電磁波傳遞能量[1]。太陽能一般指太陽光的輻射能量,可以轉(zhuǎn)化為熱能、機(jī)械能、電能或者化學(xué)能,被人們所利用。2014年,澳大利亞研究人員已經(jīng)可以將太陽能轉(zhuǎn)化為電能的效率提高到40%以上。在中國,太陽能產(chǎn)業(yè)規(guī)模已位居世界第一,尤其是太陽能光伏發(fā)電系統(tǒng)和太陽能熱水系統(tǒng)[2]。同時(shí),太陽輻射是影響建筑室內(nèi)熱環(huán)境的重要外擾因素,在建筑室內(nèi)傳熱分析中必須加以準(zhǔn)確的考慮[3]。

    計(jì)算流體動力學(xué)(CFD)是建立在經(jīng)典流體動力學(xué)與數(shù)值計(jì)算方法基礎(chǔ)上的一門新型獨(dú)立學(xué)科[3]。與實(shí)驗(yàn)測量方法相比,CFD方法不受試驗(yàn)條件的限制,有較多的靈活性,可以拓寬試驗(yàn)研究的范圍,減少昂貴的試驗(yàn)成本。在眾多涉及傳熱過程的模擬計(jì)算中,都會用到輻射模型,比如根據(jù)文獻(xiàn)[4]和文獻(xiàn)[5]研究可知,在室內(nèi)熱環(huán)境模擬和擴(kuò)散爐內(nèi)部溫度場計(jì)算上,合理選擇輻射模型,均能得到較好的計(jì)算結(jié)果。本文對人工太陽輻射環(huán)境下的傳熱過程進(jìn)行模擬,通過模擬與試驗(yàn)的比較尋找最佳的輻射模型。

    1理論模型

    本文涉及到的問題是一個(gè)三維、定常、湍流流動和傳熱的過程。

    1.1 物理模型的數(shù)學(xué)描述

    流體流動要受質(zhì)量守恒定律、動量守恒定律和能量守恒定律的支配[6],它們的雷諾時(shí)均方程為

    (1)

    (2)

    (3)

    式中xi——代表直角坐標(biāo)下的三個(gè)坐標(biāo)分量;

    ui——代表直角坐標(biāo)系下瞬時(shí)速度的三個(gè)分量;

    ρ——密度;

    t——時(shí)間;

    μ——動力粘度;

    p——壓力;

    T——溫度;

    k——傳熱系數(shù);

    cp——比熱容;

    S——廣義源項(xiàng)。

    值得說明的是:

    (1)在Z方向動量方程的源項(xiàng)Sz中包括浮升力項(xiàng),采用Boussinesq假設(shè),將流體密度視為常數(shù)。

    (3)在能量方程的源項(xiàng)ST中包括了輻射熱量,由各種輻射模型求得。

    1.2 輻射模型

    輻射換熱研究一直以來都是一個(gè)科研難題,因?yàn)檩椛鋼Q熱問題有很高的非線性,描述輻射傳遞的方程為微分—積分方程,很難獲得解析解[7]。因此,工程實(shí)際中的輻射問題一般都作一定的簡化處理,形成不同的簡化模型。目前,主要的輻射模型包括:離散傳播(DT)輻射模型、基于球形諧波法的P-1輻射模型、羅斯蘭德(Rosseland)輻射模型、表面(S2S)輻射模型、離散坐標(biāo)(DO)輻射模型[8-10]。這些模型在模擬的精度、合理性和計(jì)算量上各有特點(diǎn)。

    DT模型的主要假設(shè)是用單一的輻射射線代替從輻射表面沿某個(gè)立體角的所有輻射效應(yīng)[11]。把體積微元向周圍的熱輻射均勻的離散成有限能束,求得他們的輻射強(qiáng)度變化之后,對所有的輻射變化求和,就得到對應(yīng)于每個(gè)流體單元內(nèi)的由于輻射所引起的能量源項(xiàng)。此模型的計(jì)算精度主要由所跟蹤射線的數(shù)目以及計(jì)算網(wǎng)格密度決定。

    P-1模型是P-N模型中最簡單的類型。P-N模型的出發(fā)點(diǎn)是把輻射強(qiáng)度展開成為正交的球諧函數(shù)。對于P-1模型,我們需要求解一個(gè)輻射輸運(yùn)方程并將所得的輻射熱量直接帶入能量方程的源項(xiàng)。

    Rosseland模型引入了與溫度成三次方的傳熱系數(shù)來計(jì)算輻射熱量。由于Rosseland模型不需要計(jì)算輻射強(qiáng)度的輸運(yùn)方程,所以它的計(jì)算量比P-1模型還小。

    S2S模型可計(jì)算出在封閉區(qū)域內(nèi)的漫灰表面之間的輻射換熱。某個(gè)表面接受到其余表面的入射輻射量用角系數(shù)來度量。角系數(shù)的含義就是離開表面1的輻射量被表面2所接收到的比例。該模型的主要假定是忽略了所有的輻射吸收、發(fā)射和散射,因此,模型中僅考慮表面之間的輻射傳熱。

    DO模型的主要思想是對輻射強(qiáng)度的方向變化進(jìn)行離散,通過求解覆蓋整個(gè)4π空間角的一套離散方向上的輻射輸運(yùn)方程而得到問題的解??臻g中某一位置的4π空間角的每個(gè)象限被分割成Nθ×Nφ個(gè)輻射立體角方向,θ、φ分別為經(jīng)/緯度角。有多少個(gè)立體角方向,DO模型就求解多少個(gè)輸運(yùn)方程。立體角的離散精度決定求解的精度,但是,增加立體角的精度會使計(jì)算量急劇增加。

    2物理模型與研究方法

    2.1 物理模型

    本文的研究對象是某沙漠環(huán)境模擬實(shí)驗(yàn)室,如圖1所示。其地面中心區(qū)域?yàn)?.28 m×5.04 m(X×Y)的沙坑,內(nèi)填沙漠原型沙。房間上方裝有全光譜日光模擬器,作為模擬太陽輻射的光源,其由188只反射型鏑燈組成,光源總功率為75.2 kW,平均功率為4.19 kW/m2。為了模擬自然風(fēng)環(huán)境,南側(cè)放置3臺軸流風(fēng)機(jī)送風(fēng),送風(fēng)溫度為5℃,送風(fēng)速度為1 m/s,并利用接近頂部的兩個(gè)排風(fēng)口自然排風(fēng)。

    模型采用分區(qū)混合網(wǎng)格技術(shù)建立,在模擬器周圍使用四面體非結(jié)構(gòu)化網(wǎng)格,其余區(qū)域使用六面體結(jié)構(gòu)化網(wǎng)格劃分。結(jié)構(gòu)化網(wǎng)格占網(wǎng)格總數(shù)的75.42%,圖2顯示了Y=3.475 m斷面的網(wǎng)格分布。這種網(wǎng)格劃分方法可以減少節(jié)點(diǎn),節(jié)省計(jì)算時(shí)間,同時(shí)保證主流計(jì)算區(qū)域的計(jì)算精度。

    圖1 物理模型示意圖

    圖2 Y=3.475 m截面網(wǎng)格分布圖

    2.2 研究方法

    試驗(yàn)采用SR5太陽輻射傳感器,其光譜范圍為300~3 000 nm,測量范圍為0~2000 W/m2,精度為±1 W/m2。如圖3所示,在日光模擬器投影到地面上的區(qū)域布置30(6×5)個(gè)測點(diǎn)進(jìn)行試驗(yàn),測量地面入射輻射值。

    圖3 地面輻射測點(diǎn)布置圖

    數(shù)值模擬采用FLUENT軟件,湍流模型采用Realizable k-ε湍流模型。為考慮溫差引起的浮升力的影響,采用Boussinesq假設(shè)。日光模擬器采用定熱流邊界條件,房間墻體設(shè)為絕熱表面,地面取為第三類邊界條件,傳熱系數(shù)取為k=0.58 W/(m2·℃),進(jìn)風(fēng)口設(shè)為速度入口,排風(fēng)口設(shè)為壓力出口。

    首先,將實(shí)驗(yàn)所得的地面入射輻射值與數(shù)值模擬結(jié)果進(jìn)行比較。其次,選取Y方向上的中點(diǎn)(Y=3.475 m)做XZ截面,輸出數(shù)值模擬的溫度云圖,分析該二維空間的溫度分布。

    3試驗(yàn)與模擬結(jié)果

    3.1 地面入射輻射

    將地面輻射的實(shí)測數(shù)據(jù)與數(shù)值模擬數(shù)據(jù)繪制成為以地面輻射測點(diǎn)為X軸,各測點(diǎn)的入射輻射值為Y軸的曲線圖,如圖4所示。由于Rosseland模型沒有入射輻射值的輸出,所以只比較其余四種模型。從試驗(yàn)曲線1至5點(diǎn)、6至10點(diǎn)等六段曲線可以看出:輻射分布呈拋物線狀,且13、18點(diǎn)輻射值最大,說明地面中心入射輻射最強(qiáng),向四周逐漸變?nèi)?。DT模型的模擬值均比實(shí)測值高,且曲線沒有規(guī)律性;P-1模型的模擬值也均比實(shí)測值高,曲線平行于X軸;S2S模型的模擬值比實(shí)測值偏低,曲線有微弱的起伏趨勢;DO模型的模擬值與實(shí)測值最接近,可見使用DO模型模擬地面入射輻射值較為可靠。

    圖4 地面入射輻射分布圖

    3.2 溫度分布

    圖5給出了房間Y=3.475 m截面的溫度云圖,由圖5(a)未考慮輻射時(shí)的情況可知,日光模擬器表面溫度最高;在日光模擬器上部溫度明顯高于下部,在上部能看到熱氣流上升現(xiàn)象;在下部溫度由上至下出現(xiàn)分層,但大部分區(qū)域溫度接近風(fēng)機(jī)送風(fēng)溫度,這與實(shí)際不符,可見輻射效應(yīng)不可忽略。

    由圖5(b)、圖5(c)可以看出,P-1模型中日光模擬器上部有微弱熱氣流上升現(xiàn)象,下部大部分區(qū)域溫度等于送風(fēng)溫度,這與未考慮輻射時(shí)的情況相似; Rosseland模型中日光模擬器周圍的溫度場呈橢圓形分布,這與實(shí)際情況完全不符。本模型空間的流體介質(zhì)為空氣,光學(xué)深度很小,而這兩個(gè)模型不適合光學(xué)深度小的介質(zhì)計(jì)算,尤其是Rosseland模型要求光學(xué)深度大于3。因此,這兩種模型不適合該問題的求解。

    由圖5(d)、圖5(e)、圖5(f)可以看出,使用了S2S模型、DT模型和DO模型的溫度場趨勢相似。在日光模擬器上部都能看到熱氣流上升現(xiàn)象,在下部溫度場均呈舌狀分布,一部分熱量以輻射的方式從日光模擬器下部散出。S2S模型、DT模型中上部與未考慮輻射時(shí)的情況相似,而DO模型上部溫度均低于未考慮輻射時(shí)的情況。這是因?yàn)榭諝饨橘|(zhì)中存在散射,而S2S模型和DT模型均忽略了散射項(xiàng),所以會造成溫度場失真。但是DO模型考慮了散射項(xiàng),同時(shí)它能夠求解所有光學(xué)深度區(qū)間的輻射問題,所以模擬結(jié)果更為合理。

    圖5 Y=3.475 m斷面溫度分布圖

    3.3 求解性能

    在模擬過程中,只有DO模型增加了收斂的難度,即需要更多的迭代步數(shù)才能收斂,其余模型均未有明顯增加。DO模型對電腦資源的使用量最大,因?yàn)閷τ?D問題,總共有8Nθ×Nφ個(gè)立體角方向,在默認(rèn)情況下,Nθ和Nφ的數(shù)目均為2,故要求解32個(gè)輸運(yùn)方程,所以計(jì)算量是未使用輻射模型時(shí)的四倍多;S2S模型雖然不用求解輸運(yùn)方程,但是要計(jì)算角系數(shù),在缺省情況下,有多少個(gè)網(wǎng)格表面,就要計(jì)算多少個(gè)角系數(shù),故計(jì)算量增加了一倍多;其次是DT模型,雖然也不用求解輸運(yùn)方程,但要計(jì)算射線每個(gè)行程的輻射強(qiáng)度;P-1模型增加了1個(gè)輸運(yùn)方程,計(jì)算量略有增加;Rosseland模型無方程增加,故無計(jì)算量增加。

    雖然DO模型對計(jì)算機(jī)的硬件要求較高,但是它的精度也最高?,F(xiàn)在雙核CPU、4GB內(nèi)存的計(jì)算機(jī)已經(jīng)可以滿足100萬左右網(wǎng)格、帶有DO模型的三維問題的計(jì)算,而且隨著計(jì)算機(jī)計(jì)算能力的迅速發(fā)展,更大型的計(jì)算也將不成問題。

    4結(jié)論

    經(jīng)過了上面的分析比較后,得到的結(jié)論是:

    (1)用DO模型模擬地面入射輻射與實(shí)驗(yàn)結(jié)果較為吻合;S2S模型的模擬值偏低;DT模型和P-1模型的模擬值偏高,且趨勢不同;Rosseland模型沒有入射輻射值的輸出。

    (2)DO模型最適合太陽輻射傳熱過程溫度場的模擬;在輻射強(qiáng)度不大或要求精度不高的前提下,為了節(jié)約計(jì)算機(jī)資源,可以使用DT模型或S2S模型,但必須保證散射項(xiàng)可以忽略;P-1模型和Rosseland模型不適合此類輻射傳熱模擬。

    參考文獻(xiàn)

    [1]楊世銘,陶文銓.傳熱學(xué)[M].3版.北京:高等教育出版社,1980.

    [2]陳一誠,王志新.適用于海島的獨(dú)立供電電源規(guī)劃設(shè)計(jì)[J].電網(wǎng)與清潔能源,2013,29(4):78-84.

    [3]渠滿,職更辰,屈國建.一種多功能光優(yōu)發(fā)電系統(tǒng)的可行性研究[J].電網(wǎng)與清潔能源,2013,29(1):93-96.

    [4]郭磊.對FLUENT輻射模型的數(shù)值計(jì)算與分析[J].制冷與空調(diào),2014,28(3):358-360.

    [5]侯立宏,劉祖明,張衛(wèi)東,等.太陽電池生產(chǎn)用管式擴(kuò)散爐的數(shù)值模擬:輻射模型的選擇[J].云南師范大學(xué)學(xué)報(bào),2012,32(2):25-28.

    [6]王福軍.計(jì)算流體動力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004.

    [7]姚磊,楊晶,王立.散射性介質(zhì)輻射傳熱的數(shù)值模擬[J].冶金能源,2006,25(12):26-30.

    [8]Fluent 6.2 User’s Guide. USA: Fluent Inc, 2005.

    [9]A. Habibi, B. Merci, G.J. Heynderickx. Impact of radiation models in CFD simulations of steam cracking furnaces[J].Computers and Chemical Engineering,2007(31):1389-1406.

    [10]李鐵,李偉力,袁竹林.用不同輻射模型研究下降管內(nèi)傳熱傳質(zhì)特性[J].中國電機(jī)工程學(xué)報(bào),2007,27(2):92-98.

    [11]N. W. Bressloff, J. B. Moss, P. A. Rubini. CFD Prediction of Coupled Radiation Heat Transfer and Soot Production in Turbulent Flames[J].Twenty-Sixth Symposium (International) on Combustion, The Combustion Institute,1996:2379-2386.

    猜你喜歡
    太陽輻射數(shù)值模擬
    邯鄲太陽輻射時(shí)空分布特征
    基于PCA 的太陽輻射觀測算法研究
    太陽輻射作用下鋼筒倉結(jié)構(gòu)溫度場分析研究
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設(shè)計(jì)與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    洛陽地區(qū)太陽輻射變化特征及影響因子分析
    河南科技(2015年11期)2015-03-11 16:25:00
    中国美女看黄片| 99久久中文字幕三级久久日本| 午夜福利欧美成人| 欧美日韩精品成人综合77777| 精华霜和精华液先用哪个| 亚洲专区中文字幕在线| 3wmmmm亚洲av在线观看| 99热只有精品国产| 国产精品自产拍在线观看55亚洲| 女人十人毛片免费观看3o分钟| 毛片一级片免费看久久久久 | 毛片女人毛片| 国产成人aa在线观看| 少妇人妻精品综合一区二区 | 日日夜夜操网爽| 波多野结衣巨乳人妻| 日韩在线高清观看一区二区三区 | 免费观看在线日韩| 亚洲av中文av极速乱 | 一进一出抽搐动态| 亚洲美女搞黄在线观看 | 久久久久久久久久久丰满 | av中文乱码字幕在线| 国产精品福利在线免费观看| 精华霜和精华液先用哪个| 久久精品国产亚洲av香蕉五月| 我的女老师完整版在线观看| 听说在线观看完整版免费高清| 午夜亚洲福利在线播放| 亚洲av免费高清在线观看| 一级黄片播放器| 91精品国产九色| 老司机福利观看| 欧美另类亚洲清纯唯美| 中文字幕久久专区| 中文字幕免费在线视频6| 国产aⅴ精品一区二区三区波| 成人综合一区亚洲| 国产精品国产高清国产av| av.在线天堂| 日日啪夜夜撸| 精品一区二区免费观看| 国产蜜桃级精品一区二区三区| 黄片wwwwww| 五月玫瑰六月丁香| 日日干狠狠操夜夜爽| 欧美性猛交黑人性爽| 偷拍熟女少妇极品色| 搡老妇女老女人老熟妇| 久久九九热精品免费| 两人在一起打扑克的视频| 成人欧美大片| 精品一区二区三区人妻视频| 全区人妻精品视频| 日日啪夜夜撸| АⅤ资源中文在线天堂| 狠狠狠狠99中文字幕| 亚洲最大成人手机在线| 成年女人毛片免费观看观看9| 国产高清三级在线| 一个人看视频在线观看www免费| 成人av一区二区三区在线看| 国产男靠女视频免费网站| 深夜a级毛片| 免费不卡的大黄色大毛片视频在线观看 | 久久天躁狠狠躁夜夜2o2o| 久久国内精品自在自线图片| 成人二区视频| 中文字幕人妻熟人妻熟丝袜美| 国产高清视频在线观看网站| 一进一出好大好爽视频| 国产一区二区三区av在线 | 不卡视频在线观看欧美| 国产成人a区在线观看| 国产精品永久免费网站| 男人和女人高潮做爰伦理| 琪琪午夜伦伦电影理论片6080| 欧美日韩乱码在线| 国产一区二区亚洲精品在线观看| 欧美日韩精品成人综合77777| 此物有八面人人有两片| 欧美区成人在线视频| 亚洲av二区三区四区| 亚洲无线观看免费| 真人做人爱边吃奶动态| 高清日韩中文字幕在线| 91久久精品国产一区二区成人| 波多野结衣高清作品| 黄色女人牲交| 婷婷精品国产亚洲av在线| 动漫黄色视频在线观看| netflix在线观看网站| 日韩一区二区视频免费看| 乱码一卡2卡4卡精品| 国产在视频线在精品| 亚洲精品456在线播放app | 国产69精品久久久久777片| 久久久色成人| 高清毛片免费观看视频网站| 九九热线精品视视频播放| 97超视频在线观看视频| 久久精品国产亚洲av香蕉五月| 乱人视频在线观看| 国产精华一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 色5月婷婷丁香| 国产高清视频在线播放一区| 哪里可以看免费的av片| 桃红色精品国产亚洲av| 韩国av一区二区三区四区| 亚洲第一区二区三区不卡| 色噜噜av男人的天堂激情| a级毛片免费高清观看在线播放| 国产一区二区三区在线臀色熟女| 非洲黑人性xxxx精品又粗又长| 淫妇啪啪啪对白视频| 色综合色国产| 精品人妻偷拍中文字幕| 亚洲乱码一区二区免费版| 亚洲 国产 在线| 在线国产一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 午夜免费激情av| 好男人在线观看高清免费视频| 三级毛片av免费| 国产淫片久久久久久久久| 88av欧美| 麻豆国产97在线/欧美| 久久久久久大精品| 日本免费a在线| 国产真实乱freesex| 春色校园在线视频观看| 久久6这里有精品| or卡值多少钱| 美女免费视频网站| 久久国产乱子免费精品| 久久久久久久久久成人| 国产成人a区在线观看| 久久人妻av系列| 99久久成人亚洲精品观看| 精品久久久久久久末码| 女同久久另类99精品国产91| 伦理电影大哥的女人| 联通29元200g的流量卡| 国内精品久久久久久久电影| 别揉我奶头~嗯~啊~动态视频| 最近视频中文字幕2019在线8| 亚洲最大成人av| 97超视频在线观看视频| 亚洲最大成人av| 九色成人免费人妻av| 无遮挡黄片免费观看| 搡老熟女国产l中国老女人| 窝窝影院91人妻| 亚洲人成伊人成综合网2020| 最近最新免费中文字幕在线| 色播亚洲综合网| 国产一区二区在线av高清观看| 日韩强制内射视频| 99热网站在线观看| 色在线成人网| 国内精品一区二区在线观看| 内射极品少妇av片p| 精品福利观看| 亚洲av二区三区四区| 免费在线观看影片大全网站| 国产毛片a区久久久久| 免费无遮挡裸体视频| 午夜福利在线观看免费完整高清在 | 黄色一级大片看看| a级一级毛片免费在线观看| 99热这里只有精品一区| 美女高潮的动态| 18禁黄网站禁片午夜丰满| 别揉我奶头~嗯~啊~动态视频| 午夜激情福利司机影院| 啦啦啦啦在线视频资源| 亚洲色图av天堂| 久久久久国产精品人妻aⅴ院| 天天一区二区日本电影三级| 日韩中字成人| 日韩,欧美,国产一区二区三区 | 在线观看美女被高潮喷水网站| 久久人人爽人人爽人人片va| 国产高清有码在线观看视频| 欧美一区二区亚洲| 在线观看美女被高潮喷水网站| 淫秽高清视频在线观看| 自拍偷自拍亚洲精品老妇| 欧美bdsm另类| 免费看a级黄色片| 成年女人看的毛片在线观看| 白带黄色成豆腐渣| 老司机深夜福利视频在线观看| av天堂中文字幕网| 日本-黄色视频高清免费观看| 免费在线观看影片大全网站| 搡老岳熟女国产| 国产精品免费一区二区三区在线| 最近视频中文字幕2019在线8| 亚洲aⅴ乱码一区二区在线播放| 亚洲va在线va天堂va国产| 给我免费播放毛片高清在线观看| 国内久久婷婷六月综合欲色啪| 成熟少妇高潮喷水视频| 国产精品久久电影中文字幕| 亚洲精品456在线播放app | 99久久精品一区二区三区| 18+在线观看网站| 亚洲va日本ⅴa欧美va伊人久久| 国产av麻豆久久久久久久| 成人永久免费在线观看视频| 一个人看的www免费观看视频| 看黄色毛片网站| 日日摸夜夜添夜夜添av毛片 | 18禁黄网站禁片免费观看直播| 大型黄色视频在线免费观看| 97碰自拍视频| 婷婷精品国产亚洲av在线| 美女免费视频网站| 丰满的人妻完整版| 一a级毛片在线观看| 国产男人的电影天堂91| 国产伦精品一区二区三区视频9| 22中文网久久字幕| 国产精品久久久久久亚洲av鲁大| 日本成人三级电影网站| 日日摸夜夜添夜夜添av毛片 | 女的被弄到高潮叫床怎么办 | 99热这里只有是精品在线观看| 午夜福利成人在线免费观看| 亚洲男人的天堂狠狠| 日本三级黄在线观看| 日本黄色片子视频| 亚洲欧美精品综合久久99| 久久精品国产鲁丝片午夜精品 | 动漫黄色视频在线观看| 嫩草影视91久久| 夜夜看夜夜爽夜夜摸| 亚洲乱码一区二区免费版| 99热网站在线观看| 国产私拍福利视频在线观看| 免费电影在线观看免费观看| 琪琪午夜伦伦电影理论片6080| 美女xxoo啪啪120秒动态图| 久久久久久久精品吃奶| 久久久久久久久久黄片| 国产一区二区三区av在线 | 99热精品在线国产| 午夜日韩欧美国产| 赤兔流量卡办理| 国内精品久久久久精免费| 国产精品98久久久久久宅男小说| 22中文网久久字幕| 国产一区二区在线观看日韩| 亚洲av.av天堂| 日本 欧美在线| 最近最新中文字幕大全电影3| 在线播放无遮挡| 亚洲精品一区av在线观看| 最后的刺客免费高清国语| 色尼玛亚洲综合影院| av.在线天堂| 国产一区二区在线av高清观看| 在线看三级毛片| 91狼人影院| 国产主播在线观看一区二区| a在线观看视频网站| 国产精品久久久久久久久免| 欧美日本亚洲视频在线播放| 精品久久久久久久末码| 精品一区二区三区av网在线观看| 91午夜精品亚洲一区二区三区 | 又爽又黄a免费视频| 人人妻人人澡欧美一区二区| 免费不卡的大黄色大毛片视频在线观看 | 人妻少妇偷人精品九色| 成人综合一区亚洲| 欧美一区二区精品小视频在线| 综合色av麻豆| 免费看av在线观看网站| 国产探花极品一区二区| 黄片wwwwww| 男插女下体视频免费在线播放| 波多野结衣巨乳人妻| 国内揄拍国产精品人妻在线| 露出奶头的视频| 精品人妻视频免费看| 亚洲精品日韩av片在线观看| 国产免费av片在线观看野外av| 舔av片在线| 嫩草影院入口| 久久久久国产精品人妻aⅴ院| 日韩一区二区视频免费看| 天堂影院成人在线观看| 伊人久久精品亚洲午夜| 亚洲真实伦在线观看| 国产高清有码在线观看视频| 国产麻豆成人av免费视频| 直男gayav资源| 亚洲成av人片在线播放无| 高清在线国产一区| 久久精品国产鲁丝片午夜精品 | 亚洲中文日韩欧美视频| 伦理电影大哥的女人| 亚洲av中文av极速乱 | 精品无人区乱码1区二区| 成熟少妇高潮喷水视频| 中文字幕高清在线视频| 变态另类成人亚洲欧美熟女| av在线观看视频网站免费| 久久99热这里只有精品18| 国产成年人精品一区二区| 最近最新免费中文字幕在线| 大又大粗又爽又黄少妇毛片口| 国产私拍福利视频在线观看| 免费人成视频x8x8入口观看| 天堂网av新在线| 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 一本精品99久久精品77| 精品一区二区免费观看| 午夜激情欧美在线| 黄色日韩在线| 舔av片在线| 99精品在免费线老司机午夜| 搡老岳熟女国产| 桃红色精品国产亚洲av| 国产免费男女视频| 国产成人福利小说| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区激情视频| 嫩草影院新地址| 99热这里只有精品一区| 一进一出抽搐gif免费好疼| 天天躁日日操中文字幕| a级一级毛片免费在线观看| 久久精品国产自在天天线| 精品欧美国产一区二区三| 可以在线观看毛片的网站| 国产爱豆传媒在线观看| 嫩草影院入口| 亚洲无线观看免费| 国产高清视频在线观看网站| 免费在线观看日本一区| 97超视频在线观看视频| 精品久久久噜噜| 亚州av有码| 美女高潮的动态| 全区人妻精品视频| 我的女老师完整版在线观看| 男女那种视频在线观看| 午夜福利18| 性色avwww在线观看| 成人欧美大片| 一a级毛片在线观看| 日韩 亚洲 欧美在线| 国产黄片美女视频| 久久人人爽人人爽人人片va| 亚洲人成网站高清观看| 亚洲精品日韩av片在线观看| 一本一本综合久久| 国产伦一二天堂av在线观看| 国产国拍精品亚洲av在线观看| 一级av片app| 亚洲国产高清在线一区二区三| 久久精品国产自在天天线| 久9热在线精品视频| 欧美激情国产日韩精品一区| 黄色女人牲交| 热99re8久久精品国产| 亚洲无线观看免费| 五月玫瑰六月丁香| 老司机福利观看| 欧美一区二区精品小视频在线| 成人av在线播放网站| 国产老妇女一区| 观看美女的网站| 91在线观看av| 美女高潮喷水抽搐中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产av在哪里看| 久久久久久久精品吃奶| 白带黄色成豆腐渣| 露出奶头的视频| 亚洲欧美激情综合另类| 久久热精品热| 一本久久中文字幕| 黄色一级大片看看| 色综合亚洲欧美另类图片| aaaaa片日本免费| 国产精品亚洲一级av第二区| 偷拍熟女少妇极品色| 直男gayav资源| 1000部很黄的大片| 免费观看的影片在线观看| 少妇熟女aⅴ在线视频| 亚洲精品色激情综合| 无人区码免费观看不卡| 国产一区二区三区视频了| 欧美在线一区亚洲| 亚洲在线观看片| 国产真实伦视频高清在线观看 | 中文在线观看免费www的网站| 成人高潮视频无遮挡免费网站| av在线天堂中文字幕| 欧美绝顶高潮抽搐喷水| 亚洲精品国产成人久久av| 黄色丝袜av网址大全| 村上凉子中文字幕在线| 能在线免费观看的黄片| 国产视频一区二区在线看| 舔av片在线| av在线老鸭窝| 国产精品无大码| 99在线视频只有这里精品首页| av在线蜜桃| 人妻少妇偷人精品九色| 日本与韩国留学比较| 91麻豆av在线| 亚洲av.av天堂| 淫妇啪啪啪对白视频| 麻豆国产97在线/欧美| 午夜免费男女啪啪视频观看 | 国产精品一区www在线观看 | 中文在线观看免费www的网站| 亚洲自偷自拍三级| 久久久久久久久中文| 我要看日韩黄色一级片| 不卡视频在线观看欧美| 天天躁日日操中文字幕| 不卡一级毛片| 韩国av在线不卡| 美女大奶头视频| 亚洲美女搞黄在线观看 | 毛片女人毛片| 波多野结衣高清无吗| xxxwww97欧美| 亚洲 国产 在线| 三级毛片av免费| 999久久久精品免费观看国产| 午夜亚洲福利在线播放| 国产黄片美女视频| av视频在线观看入口| 欧美日本视频| 97超视频在线观看视频| 久久久久九九精品影院| 婷婷色综合大香蕉| 白带黄色成豆腐渣| 九九久久精品国产亚洲av麻豆| 成人无遮挡网站| 亚洲欧美日韩高清在线视频| 亚洲国产高清在线一区二区三| 色哟哟·www| 级片在线观看| 亚洲一区二区三区色噜噜| 能在线免费观看的黄片| 日本黄大片高清| 国产真实乱freesex| 国产午夜福利久久久久久| 精品欧美国产一区二区三| 精品不卡国产一区二区三区| 观看美女的网站| av天堂中文字幕网| 欧美日韩乱码在线| 婷婷精品国产亚洲av| 91久久精品国产一区二区成人| 久久草成人影院| 无遮挡黄片免费观看| 亚洲成人久久爱视频| 亚洲真实伦在线观看| 校园春色视频在线观看| 在线免费观看的www视频| 男女那种视频在线观看| 国产精品,欧美在线| 啦啦啦观看免费观看视频高清| 免费在线观看影片大全网站| 欧美色视频一区免费| 欧美精品啪啪一区二区三区| 午夜亚洲福利在线播放| 波野结衣二区三区在线| 日韩av在线大香蕉| 亚洲人成网站高清观看| 久久精品国产清高在天天线| 精品日产1卡2卡| 深夜a级毛片| 99热这里只有是精品在线观看| 日韩中文字幕欧美一区二区| 成人毛片a级毛片在线播放| 久久午夜福利片| 国内精品一区二区在线观看| а√天堂www在线а√下载| 亚洲人成网站在线播| 婷婷精品国产亚洲av在线| 麻豆成人午夜福利视频| 国产伦精品一区二区三区四那| 国产一区二区在线观看日韩| 午夜激情福利司机影院| 国产爱豆传媒在线观看| 精品欧美国产一区二区三| av在线观看视频网站免费| 非洲黑人性xxxx精品又粗又长| 国内少妇人妻偷人精品xxx网站| 在现免费观看毛片| 色哟哟哟哟哟哟| 国模一区二区三区四区视频| 精品福利观看| 欧美日韩中文字幕国产精品一区二区三区| 日韩,欧美,国产一区二区三区 | 五月伊人婷婷丁香| 在线看三级毛片| 91久久精品电影网| 91麻豆av在线| 亚洲精品日韩av片在线观看| 麻豆精品久久久久久蜜桃| 精品一区二区三区视频在线观看免费| 国产亚洲精品久久久com| 在线a可以看的网站| 国产精品久久视频播放| 美女被艹到高潮喷水动态| 18禁黄网站禁片免费观看直播| 亚洲欧美精品综合久久99| 国产一区二区在线观看日韩| 亚洲一区高清亚洲精品| 欧美+亚洲+日韩+国产| 内地一区二区视频在线| 男女做爰动态图高潮gif福利片| 久久精品国产99精品国产亚洲性色| 国产精品98久久久久久宅男小说| 精品人妻熟女av久视频| 久久精品人妻少妇| 国产中年淑女户外野战色| 国产高清不卡午夜福利| 国产美女午夜福利| 国产av在哪里看| 午夜激情欧美在线| 日韩欧美一区二区三区在线观看| 麻豆一二三区av精品| 成人特级av手机在线观看| 中文字幕高清在线视频| 一a级毛片在线观看| 直男gayav资源| 国产精品久久久久久精品电影| bbb黄色大片| 一a级毛片在线观看| 欧美另类亚洲清纯唯美| 国产精品av视频在线免费观看| 国产精品av视频在线免费观看| 午夜视频国产福利| 午夜福利在线观看吧| 欧美又色又爽又黄视频| 老师上课跳d突然被开到最大视频| 九九久久精品国产亚洲av麻豆| 中文字幕av在线有码专区| 99在线人妻在线中文字幕| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 一个人看的www免费观看视频| 我要搜黄色片| 国产精品一区二区三区四区免费观看 | 黄色一级大片看看| 国国产精品蜜臀av免费| 国产在线男女| 在线观看美女被高潮喷水网站| 亚洲va日本ⅴa欧美va伊人久久| 久久久成人免费电影| 免费观看精品视频网站| 久久久久久久午夜电影| 中文字幕高清在线视频| 偷拍熟女少妇极品色| 久久久久免费精品人妻一区二区| 久久久久国内视频| 中文亚洲av片在线观看爽| 亚洲成a人片在线一区二区| 国产精品久久久久久久电影| 日本黄色片子视频| 欧美日韩精品成人综合77777| 亚洲三级黄色毛片| 久久精品综合一区二区三区| 欧美日本视频| 黄色视频,在线免费观看| 人妻制服诱惑在线中文字幕| 日本 欧美在线| 国语自产精品视频在线第100页| 欧美xxxx性猛交bbbb| 国产精品无大码| 成人亚洲精品av一区二区| 亚洲美女黄片视频| 国产综合懂色| 变态另类丝袜制服| 国产乱人伦免费视频| 久久九九热精品免费| 成人av一区二区三区在线看| 国产亚洲精品久久久com| 又粗又爽又猛毛片免费看| x7x7x7水蜜桃| 成人特级黄色片久久久久久久| 极品教师在线视频| 日本 av在线| 亚洲国产精品久久男人天堂| 国产精品精品国产色婷婷| 成人美女网站在线观看视频| 久久这里只有精品中国| 桃色一区二区三区在线观看| 国产午夜精品论理片| 自拍偷自拍亚洲精品老妇| 国产成人aa在线观看| 动漫黄色视频在线观看| 亚洲av成人av| 亚洲avbb在线观看| 嫩草影视91久久| 美女cb高潮喷水在线观看| 在线a可以看的网站| 91在线精品国自产拍蜜月| 久久精品国产亚洲av天美|