• <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
    国产不卡一卡二| 亚洲 欧美 日韩 在线 免费| 国产视频一区二区在线看| 成熟少妇高潮喷水视频| 法律面前人人平等表现在哪些方面| 久久久久久亚洲精品国产蜜桃av| 国产麻豆成人av免费视频| 国产亚洲精品第一综合不卡| netflix在线观看网站| 欧美中文日本在线观看视频| 女生性感内裤真人,穿戴方法视频| 国产av在哪里看| 日韩欧美在线乱码| 大型av网站在线播放| 久久久精品欧美日韩精品| 九色成人免费人妻av| 国产三级在线视频| 国产午夜精品论理片| 国产激情久久老熟女| 99热6这里只有精品| 18禁裸乳无遮挡免费网站照片| 婷婷六月久久综合丁香| 黑人巨大精品欧美一区二区mp4| 天堂影院成人在线观看| 国产视频内射| 亚洲精品国产精品久久久不卡| 亚洲美女黄片视频| 国产精品亚洲一级av第二区| 啦啦啦观看免费观看视频高清| 免费在线观看亚洲国产| 久久精品国产亚洲av高清一级| 欧美一区二区国产精品久久精品 | 欧美高清成人免费视频www| 蜜桃久久精品国产亚洲av| 给我免费播放毛片高清在线观看| 国模一区二区三区四区视频 | 亚洲午夜理论影院| 国内久久婷婷六月综合欲色啪| 亚洲精品国产一区二区精华液| 久久久久久亚洲精品国产蜜桃av| 禁无遮挡网站| 亚洲无线在线观看| 国产午夜精品久久久久久| 亚洲精品在线观看二区| 午夜两性在线视频| 亚洲成人免费电影在线观看| 欧美日韩瑟瑟在线播放| 国产高清videossex| 精品久久久久久久末码| 亚洲国产欧美人成| 亚洲国产看品久久| 一进一出抽搐gif免费好疼| 波多野结衣巨乳人妻| 色综合欧美亚洲国产小说| 久久草成人影院| 亚洲中文av在线| 女人高潮潮喷娇喘18禁视频| 日韩精品中文字幕看吧| 国产精品乱码一区二三区的特点| 国产黄片美女视频| 精品不卡国产一区二区三区| www国产在线视频色| av片东京热男人的天堂| 黄片大片在线免费观看| 两个人视频免费观看高清| 日韩欧美免费精品| 国产成年人精品一区二区| 免费av毛片视频| 黄色女人牲交| 在线免费观看的www视频| 身体一侧抽搐| 啦啦啦韩国在线观看视频| 午夜激情福利司机影院| 91成年电影在线观看| 国产野战对白在线观看| 一二三四在线观看免费中文在| 久久 成人 亚洲| 男人舔女人下体高潮全视频| av免费在线观看网站| 国产午夜精品论理片| 精品国产乱子伦一区二区三区| 精品午夜福利视频在线观看一区| www国产在线视频色| 国产亚洲av嫩草精品影院| 婷婷六月久久综合丁香| 男人舔奶头视频| 麻豆av在线久日| 国产精华一区二区三区| 免费观看人在逋| 国产高清激情床上av| 欧美一区二区精品小视频在线| 12—13女人毛片做爰片一| 中文亚洲av片在线观看爽| 午夜精品一区二区三区免费看| 精品久久蜜臀av无| 麻豆国产av国片精品| 国产精品久久视频播放| 国产精品精品国产色婷婷| 美女高潮喷水抽搐中文字幕| 日本 av在线| 无遮挡黄片免费观看| 亚洲一区二区三区色噜噜| 又黄又粗又硬又大视频| 熟妇人妻久久中文字幕3abv| 午夜激情av网站| 黄色a级毛片大全视频| av片东京热男人的天堂| 免费在线观看黄色视频的| 热99re8久久精品国产| 亚洲欧美精品综合久久99| 在线观看免费午夜福利视频| 老汉色av国产亚洲站长工具| 免费在线观看完整版高清| 波多野结衣高清作品| 久久久精品欧美日韩精品| a级毛片在线看网站| 欧美乱妇无乱码| 男女之事视频高清在线观看| 久久精品国产亚洲av高清一级| 国产精品 国内视频| 成人午夜高清在线视频| 成在线人永久免费视频| 亚洲欧美一区二区三区黑人| 日本一二三区视频观看| 热99re8久久精品国产| 一本一本综合久久| 欧美zozozo另类| 叶爱在线成人免费视频播放| 久久精品综合一区二区三区| 十八禁网站免费在线| 88av欧美| 国产精品,欧美在线| 搡老岳熟女国产| 久久国产乱子伦精品免费另类| 国产高清视频在线播放一区| 久久精品91蜜桃| 动漫黄色视频在线观看| 中文资源天堂在线| 亚洲av片天天在线观看| 亚洲精品一卡2卡三卡4卡5卡| 日本黄色视频三级网站网址| 又黄又粗又硬又大视频| 国产亚洲精品综合一区在线观看 | 亚洲成人精品中文字幕电影| 午夜免费观看网址| 久久久国产精品麻豆| 国产伦人伦偷精品视频| 女生性感内裤真人,穿戴方法视频| 极品教师在线免费播放| 亚洲性夜色夜夜综合| cao死你这个sao货| av福利片在线| 99精品久久久久人妻精品| 两人在一起打扑克的视频| 一卡2卡三卡四卡精品乱码亚洲| 黄色a级毛片大全视频| 9191精品国产免费久久| 亚洲天堂国产精品一区在线| 天堂影院成人在线观看| 欧美国产日韩亚洲一区| 亚洲国产精品久久男人天堂| 亚洲av电影不卡..在线观看| 人妻久久中文字幕网| 国产v大片淫在线免费观看| 亚洲第一欧美日韩一区二区三区| 欧美黄色淫秽网站| 女人被狂操c到高潮| 在线播放国产精品三级| 国内精品一区二区在线观看| 国产亚洲av高清不卡| 欧美精品亚洲一区二区| 欧美不卡视频在线免费观看 | 丁香欧美五月| 国产精品一及| 欧美性猛交╳xxx乱大交人| 日韩大码丰满熟妇| 午夜老司机福利片| 国产精品久久视频播放| 中国美女看黄片| 久久精品91无色码中文字幕| 热99re8久久精品国产| 亚洲人成网站高清观看| 亚洲欧美日韩高清在线视频| 黑人欧美特级aaaaaa片| 久久精品成人免费网站| 母亲3免费完整高清在线观看| 久9热在线精品视频| 免费在线观看黄色视频的| 麻豆国产av国片精品| 99精品久久久久人妻精品| 久久这里只有精品19| 国产区一区二久久| 久久精品国产99精品国产亚洲性色| 久久草成人影院| 日本在线视频免费播放| 日本一本二区三区精品| 成人一区二区视频在线观看| 在线观看舔阴道视频| 久久久久国产精品人妻aⅴ院| 亚洲欧美精品综合一区二区三区| 欧美又色又爽又黄视频| 黄片小视频在线播放| 欧美黄色片欧美黄色片| 亚洲成人中文字幕在线播放| 草草在线视频免费看| 亚洲av成人av| 国产精品久久视频播放| 久久久国产成人免费| 很黄的视频免费| 午夜激情福利司机影院| 精品日产1卡2卡| 男男h啪啪无遮挡| 国产精品亚洲一级av第二区| 日韩精品免费视频一区二区三区| www.999成人在线观看| 视频区欧美日本亚洲| av超薄肉色丝袜交足视频| 非洲黑人性xxxx精品又粗又长| 久久精品国产99精品国产亚洲性色| 欧美中文综合在线视频| 亚洲自偷自拍图片 自拍| 亚洲 欧美 日韩 在线 免费| 国内揄拍国产精品人妻在线| 91成年电影在线观看| 午夜福利免费观看在线| 变态另类丝袜制服| 欧美一区二区国产精品久久精品 | 在线观看午夜福利视频| 最近在线观看免费完整版| 操出白浆在线播放| 免费看日本二区| 久久午夜综合久久蜜桃| 久久久久久久午夜电影| 亚洲精品久久国产高清桃花| e午夜精品久久久久久久| 国产成人一区二区三区免费视频网站| 国产成人一区二区三区免费视频网站| 法律面前人人平等表现在哪些方面| 国产精品永久免费网站| 久久久久免费精品人妻一区二区| 我的老师免费观看完整版| 国产高清激情床上av| 久久精品综合一区二区三区| 亚洲av电影在线进入| 日本一二三区视频观看| 欧美成人免费av一区二区三区| 黄色丝袜av网址大全| 国产激情偷乱视频一区二区| 亚洲午夜精品一区,二区,三区| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品999在线| videosex国产| 久9热在线精品视频| 亚洲av成人精品一区久久| 黑人操中国人逼视频| 在线观看舔阴道视频| 黄色成人免费大全| 日韩精品青青久久久久久| avwww免费| 19禁男女啪啪无遮挡网站| 亚洲一码二码三码区别大吗| tocl精华| 色综合欧美亚洲国产小说| 怎么达到女性高潮| 岛国在线观看网站| 成人国产一区最新在线观看| 日韩精品青青久久久久久| 精华霜和精华液先用哪个| 精品国产乱码久久久久久男人| 久久这里只有精品19| 国产片内射在线| 成人精品一区二区免费| 老司机深夜福利视频在线观看| 国产精品久久久久久久电影 | 可以在线观看的亚洲视频| 法律面前人人平等表现在哪些方面| 成人精品一区二区免费| 韩国av一区二区三区四区| 在线十欧美十亚洲十日本专区| 老司机靠b影院| 免费在线观看完整版高清| 毛片女人毛片| 91麻豆av在线| 观看免费一级毛片| 亚洲va日本ⅴa欧美va伊人久久| 在线观看www视频免费| 一级a爱片免费观看的视频| 日韩三级视频一区二区三区| 亚洲黑人精品在线| 精品一区二区三区四区五区乱码| 国产高清videossex| 在线观看www视频免费| 在线视频色国产色| 国产午夜精品久久久久久| 亚洲五月天丁香| 成年免费大片在线观看| 国内精品久久久久精免费| 精品福利观看| 久久人人精品亚洲av| 天堂√8在线中文| av超薄肉色丝袜交足视频| 国产区一区二久久| 啦啦啦韩国在线观看视频| 小说图片视频综合网站| 一级a爱片免费观看的视频| 香蕉丝袜av| 别揉我奶头~嗯~啊~动态视频| 国产精品一区二区免费欧美| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 国产精华一区二区三区| 99热这里只有精品一区 | 人成视频在线观看免费观看| 国产真实乱freesex| 夜夜躁狠狠躁天天躁| www.自偷自拍.com| 男女床上黄色一级片免费看| 国产伦人伦偷精品视频| 高清毛片免费观看视频网站| 成年免费大片在线观看| 亚洲一区中文字幕在线| 色综合亚洲欧美另类图片| 日本一二三区视频观看| 国产激情欧美一区二区| 此物有八面人人有两片| 久久久精品欧美日韩精品| a在线观看视频网站| 老司机福利观看| 国产亚洲精品久久久久久毛片| 国产欧美日韩一区二区三| 日韩欧美免费精品| 桃色一区二区三区在线观看| 国产精品一区二区三区四区免费观看 | 国产视频内射| 国产精品野战在线观看| 久久久久久免费高清国产稀缺| 一区福利在线观看| 91av网站免费观看| x7x7x7水蜜桃| 国产1区2区3区精品| 老汉色∧v一级毛片| 丰满人妻一区二区三区视频av | 久久草成人影院| 国产熟女xx| 午夜福利在线在线| 色综合站精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美久久黑人一区二区| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品sss在线观看| 国产精品av久久久久免费| 女人被狂操c到高潮| 国产探花在线观看一区二区| av福利片在线观看| 欧美黑人精品巨大| 欧美在线一区亚洲| 亚洲第一电影网av| 在线观看66精品国产| 亚洲中文字幕一区二区三区有码在线看 | 国产精品永久免费网站| 亚洲熟女毛片儿| 国产精品九九99| 欧美黄色淫秽网站| 女同久久另类99精品国产91| 午夜福利免费观看在线| 狂野欧美激情性xxxx| 国内精品久久久久精免费| 老熟妇仑乱视频hdxx| 久久久国产精品麻豆| 久久久国产成人精品二区| 天堂av国产一区二区熟女人妻 | 亚洲国产欧洲综合997久久,| 日韩av在线大香蕉| 久久伊人香网站| 日本成人三级电影网站| 长腿黑丝高跟| 久久久国产成人精品二区| 亚洲国产高清在线一区二区三| 久久香蕉国产精品| 久久久精品国产亚洲av高清涩受| 精品国内亚洲2022精品成人| 老汉色∧v一级毛片| 国产高清视频在线观看网站| 亚洲欧美日韩高清专用| 亚洲 欧美一区二区三区| 国产精品综合久久久久久久免费| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 亚洲国产日韩欧美精品在线观看 | 精品久久久久久久末码| 老司机午夜福利在线观看视频| 成人永久免费在线观看视频| 亚洲乱码一区二区免费版| 一个人观看的视频www高清免费观看 | 亚洲国产看品久久| 亚洲av第一区精品v没综合| www国产在线视频色| 在线观看舔阴道视频| 午夜免费成人在线视频| 听说在线观看完整版免费高清| 一卡2卡三卡四卡精品乱码亚洲| 在线a可以看的网站| 少妇被粗大的猛进出69影院| 亚洲国产精品成人综合色| 日本三级黄在线观看| 亚洲一区二区三区色噜噜| 亚洲性夜色夜夜综合| 国产精品免费一区二区三区在线| 免费人成视频x8x8入口观看| 婷婷丁香在线五月| 又紧又爽又黄一区二区| 五月玫瑰六月丁香| 国产麻豆成人av免费视频| 91九色精品人成在线观看| 美女大奶头视频| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久久久免费视频| 午夜视频精品福利| 最近视频中文字幕2019在线8| 精品久久蜜臀av无| 一二三四在线观看免费中文在| 美女高潮喷水抽搐中文字幕| 国产69精品久久久久777片 | 欧美一级a爱片免费观看看 | 俺也久久电影网| 亚洲av成人精品一区久久| www日本黄色视频网| 成人av在线播放网站| 欧美激情久久久久久爽电影| 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合久久99| 动漫黄色视频在线观看| 亚洲熟妇熟女久久| 又爽又黄无遮挡网站| 麻豆成人av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲男人的天堂狠狠| 久久精品人妻少妇| 久久精品成人免费网站| 国产激情欧美一区二区| 欧美av亚洲av综合av国产av| 免费看十八禁软件| 日日爽夜夜爽网站| 午夜a级毛片| 三级毛片av免费| 国产视频一区二区在线看| 精品国产乱码久久久久久男人| 亚洲在线自拍视频| 美女免费视频网站| 精品熟女少妇八av免费久了| 不卡av一区二区三区| 亚洲18禁久久av| 99久久综合精品五月天人人| 日韩免费av在线播放| 男女视频在线观看网站免费 | 色噜噜av男人的天堂激情| 成人av在线播放网站| 色老头精品视频在线观看| 久久久国产精品麻豆| 国产又黄又爽又无遮挡在线| 国产成人系列免费观看| 国产欧美日韩一区二区精品| 免费无遮挡裸体视频| av有码第一页| 久久久久国产一级毛片高清牌| 中文亚洲av片在线观看爽| 久久久久久国产a免费观看| www国产在线视频色| 一边摸一边做爽爽视频免费| 国产一区二区激情短视频| 欧美丝袜亚洲另类 | 十八禁人妻一区二区| 久久精品亚洲精品国产色婷小说| 99精品欧美一区二区三区四区| 精品欧美一区二区三区在线| 悠悠久久av| 国产亚洲精品久久久久5区| 两个人看的免费小视频| 亚洲精品一区av在线观看| 午夜福利在线观看吧| 国产亚洲精品av在线| 免费高清视频大片| 色尼玛亚洲综合影院| 国产在线观看jvid| 精品一区二区三区四区五区乱码| 欧美一区二区精品小视频在线| 亚洲精品久久国产高清桃花| 国产伦一二天堂av在线观看| 黄片小视频在线播放| www国产在线视频色| 两个人的视频大全免费| 免费搜索国产男女视频| 99国产精品99久久久久| 久久精品国产亚洲av香蕉五月| 两人在一起打扑克的视频| 午夜福利高清视频| 高清毛片免费观看视频网站| 国产伦在线观看视频一区| 19禁男女啪啪无遮挡网站| 怎么达到女性高潮| 国产精品一及| 久久久久精品国产欧美久久久| 国产真人三级小视频在线观看| 很黄的视频免费| 午夜福利高清视频| 在线观看午夜福利视频| 欧美绝顶高潮抽搐喷水| 中文字幕人妻丝袜一区二区| 亚洲第一欧美日韩一区二区三区| 悠悠久久av| 精品久久久久久久末码| 欧美日韩国产亚洲二区| 国产av不卡久久| 国产精品 欧美亚洲| 女人高潮潮喷娇喘18禁视频| 成人手机av| 午夜激情av网站| 欧美在线一区亚洲| 欧美日韩国产亚洲二区| 香蕉av资源在线| 正在播放国产对白刺激| 他把我摸到了高潮在线观看| 精品午夜福利视频在线观看一区| 午夜视频精品福利| 国产精品 国内视频| 女人爽到高潮嗷嗷叫在线视频| 国产午夜精品久久久久久| 久久午夜综合久久蜜桃| 国产一区二区在线观看日韩 | 一区二区三区国产精品乱码| 美女高潮喷水抽搐中文字幕| svipshipincom国产片| 男女午夜视频在线观看| 一本一本综合久久| 亚洲精品一卡2卡三卡4卡5卡| 两个人视频免费观看高清| 免费在线观看视频国产中文字幕亚洲| 成在线人永久免费视频| 黑人操中国人逼视频| 天天添夜夜摸| 国产熟女xx| 午夜福利免费观看在线| 国产精品电影一区二区三区| 精品乱码久久久久久99久播| 亚洲精品国产一区二区精华液| 成人国产一区最新在线观看| 最近在线观看免费完整版| 色哟哟哟哟哟哟| 1024手机看黄色片| 亚洲中文字幕一区二区三区有码在线看 | 国产视频一区二区在线看| 69av精品久久久久久| 无限看片的www在线观看| or卡值多少钱| 成人一区二区视频在线观看| 男女下面进入的视频免费午夜| 久久久国产成人精品二区| 此物有八面人人有两片| 久9热在线精品视频| 亚洲专区国产一区二区| cao死你这个sao货| 亚洲va日本ⅴa欧美va伊人久久| 成在线人永久免费视频| 中文字幕高清在线视频| 九色成人免费人妻av| 欧美一级a爱片免费观看看 | 婷婷丁香在线五月| 欧美av亚洲av综合av国产av| 免费在线观看成人毛片| 国产三级中文精品| 非洲黑人性xxxx精品又粗又长| 日韩欧美国产一区二区入口| 国产成人一区二区三区免费视频网站| 此物有八面人人有两片| 国产69精品久久久久777片 | 精品国产超薄肉色丝袜足j| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 两个人视频免费观看高清| 久久久久国产精品人妻aⅴ院| 亚洲黑人精品在线| 精品不卡国产一区二区三区| 男插女下体视频免费在线播放| 亚洲中文字幕日韩| 十八禁人妻一区二区| 中文亚洲av片在线观看爽| www.www免费av| 国产精品亚洲一级av第二区| 午夜福利在线观看吧| 三级男女做爰猛烈吃奶摸视频| 一夜夜www| 欧美精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久久久毛片| 成年人黄色毛片网站| 精品电影一区二区在线| 亚洲人成伊人成综合网2020| 午夜福利在线在线| 亚洲免费av在线视频| 久久国产乱子伦精品免费另类| 老司机午夜十八禁免费视频| 精品久久久久久,| 成人三级做爰电影| 亚洲美女视频黄频| 亚洲国产精品久久男人天堂| 99久久精品热视频| 精品第一国产精品| 欧美一区二区国产精品久久精品 | 在线观看www视频免费| 国产成人aa在线观看| 国产亚洲精品第一综合不卡| 老司机福利观看| 黄色女人牲交| 啦啦啦观看免费观看视频高清| 欧美日韩中文字幕国产精品一区二区三区|