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

    密閉金屬空腔中多孔介質(zhì)相變耦合傳熱模型研究

    2020-04-22 11:18:38徐漢中鄔潤輝邵春收單喜軍龔俊利
    關(guān)鍵詞:空腔腔體工質(zhì)

    徐漢中,鄔潤輝,2,邵春收,單喜軍,龔俊利,2

    密閉金屬空腔中多孔介質(zhì)相變耦合傳熱模型研究

    徐漢中1,鄔潤輝1,2,邵春收1,單喜軍1,龔俊利1,2

    (1. 北京航天長征飛行器研究所,北京,100076;2. 試驗(yàn)物理與計(jì)算數(shù)學(xué)國家級重點(diǎn)實(shí)驗(yàn)室,北京,100076)

    針對外加瞬間高能熱源,由密閉金屬空腔、腔體中多孔介質(zhì)、多孔介質(zhì)中可相變工質(zhì)組合成的耦合態(tài)目標(biāo)傳熱特性在航空航天領(lǐng)域的特殊應(yīng)用需求,基于金屬腔體、多孔介質(zhì)中不可壓縮流動相變工質(zhì)的能量守恒、動量守恒、質(zhì)量守恒等原理,以低氣壓環(huán)境下應(yīng)用的復(fù)雜物理結(jié)構(gòu)的密閉金屬腔體為研究對象,建立了由高能高效熱源-金屬空腔-多孔介質(zhì)-多孔介質(zhì)中可相變工質(zhì)的相變過程的傳熱特性的數(shù)值計(jì)算模型,通過典型物理參數(shù)狀態(tài)的仿真,獲取了密閉腔體內(nèi)溫度、壓力、相變等特征參數(shù)隨時間和空間變化的數(shù)值模擬結(jié)果,關(guān)鍵參數(shù)的仿真結(jié)果與已有試驗(yàn)結(jié)果對比一致性好,同時溫度、壓力、相變等特征參數(shù)的變化規(guī)律與傳熱傳質(zhì)理論相吻合。結(jié)果表明該建模方案思路的正確性以及仿真結(jié)果的有效性。該傳熱模型可推廣應(yīng)用于激勵熱源的優(yōu)化設(shè)計(jì)、密閉耦合態(tài)目標(biāo)的綜合性能優(yōu)化。

    密閉金屬腔體多孔介質(zhì);相變工質(zhì);傳熱模型;數(shù)值模擬

    0 引 言

    多孔介質(zhì)廣泛應(yīng)用于熱能存貯[1]、能源再生應(yīng)用[2]、消聲、減震、隔熱和電磁屏蔽、高溫導(dǎo)熱和傳熱等領(lǐng)域[3~6],其中多孔介質(zhì)及其工質(zhì)相變的傳熱傳質(zhì)模型研究是在各個領(lǐng)域應(yīng)用必不可少的研究環(huán)節(jié),因此通過采用理論建模和地面試驗(yàn)相結(jié)合的方法形成可準(zhǔn)確模擬多孔介質(zhì)傳熱傳質(zhì)過程的仿真模型是多孔介質(zhì)應(yīng)用領(lǐng)域的重要技術(shù)支撐,國內(nèi)外分別開展了理論模型和相關(guān)試驗(yàn)研究工作[7~9],開展了多孔介質(zhì)導(dǎo)熱系數(shù)計(jì)算模型,多孔介質(zhì)內(nèi)部的導(dǎo)熱、相變與傳質(zhì)計(jì)算模型,多孔介質(zhì)材料作為一種強(qiáng)化傳熱裝置的傳熱特性分析模型,多孔介質(zhì)的孔徑分布特性的優(yōu)化模型,多孔介質(zhì)傳熱物性參數(shù)的地面模擬試驗(yàn)等研究工作,形成了多孔介質(zhì)傳熱傳質(zhì)特性的數(shù)學(xué)與數(shù)值計(jì)算仿真模擬,為不同結(jié)構(gòu)的多孔介質(zhì)傳熱傳質(zhì)特性分析提供了重要的分析手段。

    而目前多孔介質(zhì)在航空航天領(lǐng)域的新應(yīng)用,需要針對外加瞬間高能熱源,建立由密閉金屬空腔、腔體中多孔介質(zhì)、多孔介質(zhì)中可相變工質(zhì)組合成的耦合態(tài)目標(biāo)傳熱模型,目前針對多組件耦合相變傳熱模型的建模和計(jì)算方法研究成果較少。本文在繼承已有多孔介質(zhì)傳熱模型算法的基礎(chǔ)上,針對新應(yīng)用需求提出了一種適用于模擬密閉耦合目標(biāo)中,瞬間高能熱源激勵條件下,由熱源-金屬空腔-多孔介質(zhì)-多孔介質(zhì)中可相變工質(zhì)的相變過程的傳熱特性的數(shù)值計(jì)算模型建模方法,通過建立多部件耦合傳熱特性分析模型,給出了密閉耦合態(tài)目標(biāo)的全目標(biāo)區(qū)域,溫度場、壓力場、速度場以及相變工質(zhì)的相變特征隨時間和空間變化的數(shù)值模擬結(jié)果,仿真模型的建立為激勵熱源的優(yōu)化設(shè)計(jì)、密閉耦合態(tài)目標(biāo)的綜合性能優(yōu)化提供重要的理論仿真平臺。

    1 傳熱及控制模型

    本文研究的對象包含金屬腔體、多孔介質(zhì)、多孔介質(zhì)中可相變可流動的工質(zhì)等,其傳熱過程是多個物理過程的耦合,因此建立傳熱模型時涉及到非定常動量方程,對流換熱與導(dǎo)熱方程,固-液-氣的相變模型,同時需要給出輻射換熱、變溫、變熱流邊界的邊界條件,整個傳熱過程滿足質(zhì)量守恒、動量守恒以及能量守恒。

    1.1 能量守恒方程

    傳熱模型遵循能量守恒方程,具體計(jì)算表達(dá)式為[10]

    1.2 動量和質(zhì)量守恒方程

    針對本文研究對象為金屬腔體、多孔介質(zhì)、多孔介質(zhì)中相變工質(zhì)、金屬空腔中工質(zhì)蒸汽流動等特性耦合的復(fù)雜特性,分別建立多孔介質(zhì)中工質(zhì)流動、金屬空腔中工質(zhì)流動的動量守恒和質(zhì)量守恒方程,具體如下:

    a)多孔介質(zhì)中,相變工質(zhì)流動的動量守恒方程為[10]

    由于多孔介質(zhì)中相變工質(zhì)流動的是不可壓縮流體,因此質(zhì)量守恒方程為[11]

    b)金屬空腔中工質(zhì)蒸汽流動的動量守恒為

    金屬空腔中工質(zhì)蒸汽流動的質(zhì)量守恒方程同 式(7)。

    1.3 工質(zhì)相變量判據(jù)方程

    1.4 邊界條件

    多孔介質(zhì)、相變工質(zhì)產(chǎn)生相變后由多孔介質(zhì)流動至金屬空腔的傳熱特性模擬的計(jì)算邊界條件表達(dá)式如下[10]:

    在確定計(jì)算邊界條件后,通過求解上述能量守恒方程、動量守恒方程、質(zhì)量守恒方程以及工質(zhì)相變判據(jù)方程,可以完成金屬空腔中多孔介質(zhì)相變耦合傳熱特性的仿真計(jì)算和分析。

    2 典型狀態(tài)數(shù)值模擬

    2.1 仿真模型

    為了驗(yàn)證提出的金屬空腔中多孔介質(zhì)相變耦合傳熱模型的數(shù)值模型計(jì)算功能,以某裝置為例(該裝置由金屬空腔、多孔介質(zhì)、可相變工質(zhì)等組成),開展密閉金屬腔體熱激勵作用下工質(zhì)的相變過程仿真。圖1為該環(huán)節(jié)傳熱仿真模型,傳熱過程為一定溫度的初始熱源(裝置底部)將能量傳導(dǎo)至金屬殼體,進(jìn)而至多孔介質(zhì),溫度達(dá)到工質(zhì)的汽化溫度后,工質(zhì)相變?yōu)闅鈶B(tài),氣體工質(zhì)向金屬空腔內(nèi)流動,為電弧放電提供氣源。

    在仿真計(jì)算時,以初始熱源加載時刻為零時刻,傳熱仿真時間段以熱源加載為起點(diǎn),以密閉腔體氣源開始供給于電弧放電為終點(diǎn),且典型的仿真時間段設(shè)置為5 s。典型的仿真結(jié)果包括不同時刻的溫度、壓力、相變等的二維空間分布。

    圖1 典型的二維仿真模型

    2.2 數(shù)值模擬主要輸入?yún)?shù)

    計(jì)算的主要輸入?yún)?shù)為:

    a)初始熱源:均勻分布的溫度場,且整個傳熱計(jì)算時間內(nèi)均勻分布的溫度場輸出值不變,均值為2000 ℃。

    b)工質(zhì)的相變溫度:700 ℃。

    c)多孔介質(zhì)的孔隙率:70%。

    d)多孔介質(zhì)的平均孔徑:10 μm。

    2.3 典型仿真結(jié)果

    以熱源加載為零時刻,圖2~4給出了0.1 s、2 s、4 s和5 s等4個時刻溫度、壓力、相變等二維空間分布結(jié)果。

    圖2 溫度場隨時間和空間的變化結(jié)果

    續(xù)圖2

    由圖2可以看出,熱源加載后,熱傳導(dǎo)作用下附近的溫度首先升高,隨著時間的推移,多孔介質(zhì)的溫升速率相對高于中心腔體的溫升速率;在5 s時刻密閉腔體內(nèi)一半以上區(qū)域的溫度都達(dá)到900 K,此時溫度滿足了裝置試驗(yàn)需求。

    圖3 蒸汽壓力隨時間和空間的變化結(jié)果

    由圖3可以看出,密閉腔體中的初始壓力為 1×105Pa,在熱源激勵下,隨著腔體內(nèi)溫度的不斷升高,工質(zhì)相變、汽化,腔內(nèi)壓力也不斷升高,特別是熱源附近的空腔中,在5 s時局部最大壓力達(dá)到4.5×105Pa,該仿真結(jié)果與某裝置工作時工質(zhì)的汽化條件相吻合。

    圖4為不同時刻的相變特征。

    圖4 工質(zhì)相變隨時間和空間的變化結(jié)果

    由圖4可以看出,無熱源激勵條件下,腔體內(nèi)多孔介質(zhì)吸附的工質(zhì)以固態(tài)或液態(tài)存在。在熱激勵條件下,隨著溫度的不斷升高,達(dá)到工質(zhì)由固-氣、液-氣相變的溫度時,多孔介質(zhì)中工質(zhì)不斷汽化為工質(zhì)蒸汽,蒸汽由多孔介質(zhì)脫附流向金屬空腔,隨著時間的推移,相變量不斷增大,空腔中蒸汽量也不斷增大,從而腔內(nèi)壓力也增大,腔內(nèi)的工質(zhì)質(zhì)量將為電弧放電提供源源不斷的放電用氣源。在5 s時刻多孔介質(zhì)中工質(zhì)相變可達(dá)到50%作用。

    為進(jìn)一步驗(yàn)證仿真模型和仿真結(jié)果,選取典型的仿真結(jié)果和已有的試驗(yàn)結(jié)果進(jìn)行定性對比(具體見表1),結(jié)果表明已建立的傳熱計(jì)算模型可有效應(yīng)用于計(jì)算開塞前密閉陰極部件熱傳遞問題,結(jié)果吻合較好。說明該仿真模型可以指導(dǎo)工程應(yīng)用。

    表1 主要性能參數(shù)仿真結(jié)果與試驗(yàn)結(jié)果對比

    Tab.1 Comparison of Parameter Simulation & Experiment Results

    序號參數(shù)名稱仿真結(jié)果已有試驗(yàn)結(jié)果對比對比結(jié)果 1加熱裝置反應(yīng)溫度/℃20001800~2000基本吻合

    續(xù)表1

    序號參數(shù)名稱仿真結(jié)果已有試驗(yàn)結(jié)果對比對比結(jié)果 2陰極腔內(nèi)壓力/Pa5s內(nèi)陰極管附近達(dá)到4×105陰極腔內(nèi)達(dá)到4×105~4.5×105,塞子頂開基本吻合 3陰極腔內(nèi)溫度/℃5s時刻陰極腔內(nèi)溫度為1100開塞前陰極腔內(nèi)溫度為1100基本吻合

    3 結(jié)束語

    本文針對復(fù)雜物理結(jié)構(gòu)的密閉金屬腔體,提出了一種適用于熱源-金屬空腔-多孔介質(zhì)-多孔介質(zhì)中可相變工質(zhì)的相變過程的傳熱特性的數(shù)值計(jì)算模型建模方法,建立了相變耦合傳熱計(jì)算模型,理論仿真與試驗(yàn)結(jié)果對比分析驗(yàn)證了建模方法的可行性、合理性,為該類型工程裝置的優(yōu)化設(shè)計(jì)提供重要的基礎(chǔ)理論支撐。

    [1] 趙長穎, 等. 多孔介質(zhì)的相變和熱化學(xué)儲熱性能[J]. 科學(xué)通報(bào), 2016(61): 1897-1911.

    Zhao Changyin, et al. Heat transfer of phase change materials (PCMs) and thermo chemical heat storage in porous materials[J]. Chinese Science Bulletin, 2016(61): 1897-1911.

    [2] 張一江, 等. 多孔過渡金屬氧化物材料在能源環(huán)境中的應(yīng)用進(jìn)展[J]. 中國陶瓷, 2017, 53(2): 1-12.

    Zhang Yijiang, et al. Research progress on applications of porous transition metal oxide materials in energy and environment[J]. China Ceramics, 2017, 53(2): 1-12.

    [3] 謝濤, 等. 氣凝膠納米多孔隔熱材料傳熱計(jì)算模型的研究[J]. 工程熱物理學(xué)報(bào), 2014, 35(2): 299-304.

    Xie Tao, et al. Study on theoretical model for the effective thermal conductivity of silica aerogel composite insulating materials[J]. Journal of Engineering Thermophysics, 2014, 35(2): 299-304.

    [4] 張虎, 等. 氣氛壓力對納米多孔材料等效熱導(dǎo)率的影響[J]. 工程熱物理學(xué)報(bào), 2013, 34(4): 756-759.

    Zhang Hu, et al. The influence of gas pressure on the effective thermal conductivity of nano-porous material[J]. Journal of Engineering Thermophysics, 2013, 34(4): 756-759.

    [5] 李慧燕, 等. 多孔材料傳熱特性的試驗(yàn)[J]. 低溫與超導(dǎo), 2008, 36(4): 15-17.

    Li Huiyan, et al. An experimental study on heat transfer characteristics of porous material[J]. Cryogenics and Superconductivity, 2008, 36(4): 15-17.

    [6] 吳志根, 等. 多孔介質(zhì)在高溫相變蓄熱中的強(qiáng)化換熱[J]. 化工學(xué)報(bào), 2012, 63(S1): 119-122.

    Wu Zhigen, et al. Heat transfer enhancement of high temperature thermal energy storage using porous materials[J]. CIESC Journal, 2012, 63(S1): 119-122.

    [7] 袁越錦, 等. 顆粒堆積多孔介質(zhì)干燥多尺度多層結(jié)構(gòu)傳熱傳質(zhì)模型及模擬[J]. 工程熱物理學(xué)報(bào), 2015, 36(12): 2726-2729.

    Yuan Yuejin, et al. Multi-Scale and multi-layer structural modeling and simulation of heat and mass transfer processes for drying of grain packing porous media[J]. Journal of Engineering Thermophysics, 2015, 36(12): 2726-2729.

    [8] 牛利嬌, 等. 具有預(yù)制孔隙多孔介質(zhì)冷凍干燥的多相傳遞模型[J]. 化工學(xué)報(bào), 2017, 68(5): 1833-1844.

    Niu Lijiao, et al. Multiphase transport model for freeze-drying of porous media with prefabricated porosity[J]. CIESC Journal, 2017, 68(5): 1833-1844.

    [9] 左遠(yuǎn)志, 等. 多孔介質(zhì)中熔融鹽流體高溫斜溫層蓄熱的熱特性[J]. 工程熱物理學(xué)報(bào), 2014, 35(2): 305-310.

    Zuo Yuanzhi, et al. Thermal analysis on the single-phase flow through porous media in the high-temperature molten salt thermocline thermal energy storage system[J]. Journal of Engineering Thermophysics, 2014, 35(2): 305-310.

    [10] Wang Haidong, Guo Zengyuan. Thermon gas as the thermal energy carrier in gas and metals[J]. Chinese Science Bulletin, 2010, 55(29): 3350-3355.

    [11] 朱杰. 多孔介質(zhì)內(nèi)的相變傳熱傳質(zhì)過程研究[D]. 大連: 大連理工大學(xué), 2007.

    Zhu Jie. Research on heat and mass transfer in the process of phase change of porous media[D]. Dalian: Dalian University of Technology, 2007.

    Study on the Coupled Heat Transfer Model of Porous Media inEnclosed Metallic Cavity

    Xu Han-zhong1, Wu Run-hui1,2, Shao Chun-shou1, Shan Xi-jun1, Gong Jun-li1,2

    (1. Beijing Institute of Space Long March Vehicle, Beijing, 100076;2. National Key Laboratory of Science and Technology on Test Physics & Numerical Mathematics, Beijing, 100076)

    Aiming at the heat source with external instantaneous high energy.According to the special requirements in the field of aerospace, the heat transfer characteristics of the variegated medium coupling state target’s are studied by considering the enclosed metal cavity, porous medium in the cavity, phase change in the porous media etc. The simulation model is obtained, based on the heat energy conservation, momentum conservation and mass conservation principle, and the heat-transfer mode is supposed from heat source to metal cavity, then to metal cavity porous medium, porous medium’s phase, its heat-transfer and so on. The parameters of temperature, pressure, velocity and phase change are computed, and the results are in conformity with tests also with theory of heat transfer medium as the incentive heat source is constant. This reveal the heat-transfer mode is right and the simulation model is valid. The heat transfer model can be popularized and applied to the optimization design of incentive heat source.

    enclosed metal cavity; porous foam materials; phase transition working medium; heat transfer model; numerical simulation

    O531

    A

    1004-7182(2020)02-0038-06

    10.7654/j.issn.1004-7182.20200208

    徐漢中(1975-),男,高級工程師,主要研究方向?yàn)榛鸸ぱb置技術(shù)。

    鄔潤輝(1976-),女,研究員,主要研究方向?yàn)榈入x子體技術(shù)。

    邵春收(1985-),男,高級工程師,主要研究方向?yàn)榭傮w技術(shù)。

    單喜軍(1970-),男,高級工程師,主要研究方向?yàn)闄C(jī)構(gòu)結(jié)構(gòu)。

    龔俊利(1962-),男,技師,主要研究方向?yàn)樵囼?yàn)技術(shù)。

    2018-12-24;

    2019-07-22

    國家自然科學(xué)基金項(xiàng)目資助(61302029, 61571031, 61571031,61871018)

    猜你喜歡
    空腔腔體工質(zhì)
    海洋溫差能發(fā)電熱力循環(huán)系統(tǒng)的工質(zhì)優(yōu)選
    基于邊光滑有限元法的二維復(fù)合彈性空腔聲振特性分析
    高鐵復(fù)雜腔體鑄造數(shù)值仿真及控制技術(shù)研究
    高鐵制動系統(tǒng)復(fù)雜腔體鑄造成形數(shù)值模擬
    采用R1234ze(E)/R245fa的非共沸混合工質(zhì)有機(jī)朗肯循環(huán)系統(tǒng)實(shí)驗(yàn)研究
    采用二元非共沸工質(zhì)的有機(jī)朗肯循環(huán)熱力學(xué)分析
    橡膠擠出裝置
    空腔參數(shù)對重力壩穩(wěn)定的影響分析
    前置污水去油池
    前置污水去油池
    国产精品麻豆人妻色哟哟久久| 国产亚洲一区二区精品| 伦理电影免费视频| e午夜精品久久久久久久| 免费久久久久久久精品成人欧美视频| 久久久久久久大尺度免费视频| 久久青草综合色| √禁漫天堂资源中文www| 色婷婷av一区二区三区视频| 精品国产乱码久久久久久小说| av超薄肉色丝袜交足视频| 免费观看av网站的网址| 精品人妻熟女毛片av久久网站| 美女午夜性视频免费| 999久久久国产精品视频| 在线观看免费视频日本深夜| 亚洲成人手机| 成人18禁高潮啪啪吃奶动态图| 最黄视频免费看| 日韩人妻精品一区2区三区| 久久久国产精品麻豆| 黄片小视频在线播放| 国产精品久久久av美女十八| 在线亚洲精品国产二区图片欧美| av免费在线观看网站| 黄色怎么调成土黄色| 国产亚洲午夜精品一区二区久久| 成人国语在线视频| e午夜精品久久久久久久| 日本av免费视频播放| 一级片免费观看大全| 老汉色av国产亚洲站长工具| 日本黄色日本黄色录像| 亚洲精品美女久久久久99蜜臀| 国产欧美日韩一区二区三| 国产日韩欧美在线精品| 黄色 视频免费看| 午夜福利影视在线免费观看| 黄色毛片三级朝国网站| 黄色丝袜av网址大全| 国产单亲对白刺激| 国产欧美日韩综合在线一区二区| 久久香蕉激情| a在线观看视频网站| 天堂中文最新版在线下载| 黄色视频,在线免费观看| 国产精品九九99| 宅男免费午夜| 91麻豆精品激情在线观看国产 | 国产精品av久久久久免费| 岛国毛片在线播放| 成人亚洲精品一区在线观看| 一区福利在线观看| 女人被躁到高潮嗷嗷叫费观| 久久国产精品影院| 色婷婷久久久亚洲欧美| 色老头精品视频在线观看| 丝袜美足系列| 亚洲自偷自拍图片 自拍| 精品国内亚洲2022精品成人 | 久久精品人人爽人人爽视色| 成年人免费黄色播放视频| 99re在线观看精品视频| 亚洲成人免费电影在线观看| 国产成人av激情在线播放| 色综合欧美亚洲国产小说| 99久久精品国产亚洲精品| 中文字幕人妻熟女乱码| 69av精品久久久久久 | 国产激情久久老熟女| 三级毛片av免费| 精品熟女少妇八av免费久了| 后天国语完整版免费观看| 久久亚洲真实| 人成视频在线观看免费观看| 一夜夜www| 叶爱在线成人免费视频播放| 国产精品自产拍在线观看55亚洲 | 国产淫语在线视频| 在线观看一区二区三区激情| 在线亚洲精品国产二区图片欧美| 亚洲专区中文字幕在线| 国产成人欧美| 最新的欧美精品一区二区| 国产主播在线观看一区二区| 久久精品aⅴ一区二区三区四区| 国产精品香港三级国产av潘金莲| 天天影视国产精品| 国产主播在线观看一区二区| kizo精华| 午夜福利一区二区在线看| 天天影视国产精品| 极品少妇高潮喷水抽搐| 日韩成人在线观看一区二区三区| 在线观看人妻少妇| 亚洲国产毛片av蜜桃av| 色在线成人网| 国产成人系列免费观看| 男人舔女人的私密视频| 日韩欧美国产一区二区入口| 精品亚洲成a人片在线观看| 波多野结衣av一区二区av| 极品少妇高潮喷水抽搐| 久久久久久久大尺度免费视频| 日韩一区二区三区影片| 国产日韩欧美视频二区| 国产一区二区三区视频了| 美女视频免费永久观看网站| 久久久精品区二区三区| 午夜两性在线视频| 动漫黄色视频在线观看| 欧美 日韩 精品 国产| 国产精品久久久久久精品古装| 丰满饥渴人妻一区二区三| 丝袜在线中文字幕| 欧美亚洲日本最大视频资源| 亚洲国产成人一精品久久久| 国产成人欧美| 精品一区二区三区视频在线观看免费 | 久久久久久亚洲精品国产蜜桃av| 国产麻豆69| 久久久久视频综合| 亚洲综合色网址| 午夜激情久久久久久久| 国产aⅴ精品一区二区三区波| 美女扒开内裤让男人捅视频| 精品少妇久久久久久888优播| 国产一区二区 视频在线| 国产精品一区二区免费欧美| 窝窝影院91人妻| 精品久久久精品久久久| 中文字幕精品免费在线观看视频| 在线十欧美十亚洲十日本专区| 美国免费a级毛片| 欧美乱码精品一区二区三区| 久久精品亚洲精品国产色婷小说| 精品福利观看| 999久久久精品免费观看国产| 丝袜在线中文字幕| 涩涩av久久男人的天堂| 国产真人三级小视频在线观看| 久久亚洲精品不卡| 国产免费av片在线观看野外av| 亚洲精品美女久久久久99蜜臀| 日本黄色日本黄色录像| 国产免费福利视频在线观看| 丝袜美腿诱惑在线| 母亲3免费完整高清在线观看| 如日韩欧美国产精品一区二区三区| 99热国产这里只有精品6| 久久这里只有精品19| 夜夜夜夜夜久久久久| 亚洲成a人片在线一区二区| 国产一区二区在线观看av| 欧美成人免费av一区二区三区 | 露出奶头的视频| 成人亚洲精品一区在线观看| 90打野战视频偷拍视频| 18禁美女被吸乳视频| 高清毛片免费观看视频网站 | 正在播放国产对白刺激| 亚洲午夜理论影院| 黑人操中国人逼视频| 99久久99久久久精品蜜桃| 免费日韩欧美在线观看| 精品国产超薄肉色丝袜足j| 韩国精品一区二区三区| 午夜福利视频在线观看免费| 成年动漫av网址| 精品少妇黑人巨大在线播放| 久久人妻福利社区极品人妻图片| 夜夜爽天天搞| 久久久国产成人免费| 两人在一起打扑克的视频| 午夜成年电影在线免费观看| 国产aⅴ精品一区二区三区波| 男男h啪啪无遮挡| 欧美日韩av久久| 黄网站色视频无遮挡免费观看| 亚洲伊人色综图| 亚洲欧美色中文字幕在线| 麻豆成人av在线观看| 不卡一级毛片| 丁香欧美五月| 亚洲第一青青草原| 俄罗斯特黄特色一大片| av网站免费在线观看视频| 免费在线观看黄色视频的| 国产一区二区在线观看av| 黑人猛操日本美女一级片| 成人精品一区二区免费| 午夜激情av网站| 伊人久久大香线蕉亚洲五| 久久精品熟女亚洲av麻豆精品| 亚洲av美国av| 色94色欧美一区二区| 汤姆久久久久久久影院中文字幕| tocl精华| 久久久久国产一级毛片高清牌| 青青草视频在线视频观看| 搡老熟女国产l中国老女人| 久久国产精品影院| 欧美亚洲日本最大视频资源| 搡老岳熟女国产| 美女高潮喷水抽搐中文字幕| 狠狠狠狠99中文字幕| 精品国产乱码久久久久久小说| 亚洲国产欧美在线一区| 中文字幕高清在线视频| 在线观看免费视频日本深夜| 国产精品久久久av美女十八| 岛国毛片在线播放| 久久免费观看电影| 精品国产一区二区久久| 亚洲欧洲精品一区二区精品久久久| 一级毛片女人18水好多| 久久国产精品人妻蜜桃| 午夜91福利影院| 久久久精品区二区三区| 久久精品aⅴ一区二区三区四区| 女性生殖器流出的白浆| 国产成人精品无人区| 18禁黄网站禁片午夜丰满| 国产单亲对白刺激| 亚洲精品av麻豆狂野| 亚洲精品一卡2卡三卡4卡5卡| 亚洲 国产 在线| 国产一卡二卡三卡精品| 久久人妻福利社区极品人妻图片| 麻豆国产av国片精品| 热99re8久久精品国产| 欧美大码av| 美国免费a级毛片| 欧美日韩亚洲高清精品| 国产日韩欧美在线精品| 欧美激情极品国产一区二区三区| 母亲3免费完整高清在线观看| 亚洲人成电影免费在线| 国产精品一区二区在线不卡| 女人爽到高潮嗷嗷叫在线视频| 超碰97精品在线观看| 日韩欧美一区二区三区在线观看 | 美女视频免费永久观看网站| 精品人妻1区二区| 亚洲熟妇熟女久久| 十八禁人妻一区二区| 久久精品国产99精品国产亚洲性色 | 美女午夜性视频免费| 亚洲精品美女久久久久99蜜臀| 麻豆av在线久日| 99久久精品国产亚洲精品| 啦啦啦视频在线资源免费观看| 一边摸一边做爽爽视频免费| tube8黄色片| 1024香蕉在线观看| 欧美在线一区亚洲| 日本一区二区免费在线视频| 我的亚洲天堂| 日韩一区二区三区影片| 天天影视国产精品| 国产免费现黄频在线看| tocl精华| 婷婷丁香在线五月| 9191精品国产免费久久| 三上悠亚av全集在线观看| 亚洲av美国av| 汤姆久久久久久久影院中文字幕| √禁漫天堂资源中文www| 天天躁日日躁夜夜躁夜夜| 亚洲国产欧美日韩在线播放| 欧美变态另类bdsm刘玥| 黄频高清免费视频| 极品人妻少妇av视频| 啦啦啦视频在线资源免费观看| 亚洲欧美激情在线| 美女午夜性视频免费| 免费观看人在逋| 别揉我奶头~嗯~啊~动态视频| 成人影院久久| 美女国产高潮福利片在线看| 国产野战对白在线观看| 91成年电影在线观看| 一区在线观看完整版| 一级片'在线观看视频| 国产精品国产高清国产av | 欧美日韩黄片免| 蜜桃在线观看..| 1024视频免费在线观看| 精品午夜福利视频在线观看一区 | 欧美亚洲日本最大视频资源| 国产主播在线观看一区二区| av网站免费在线观看视频| 色尼玛亚洲综合影院| 久久精品国产99精品国产亚洲性色 | 国产又爽黄色视频| 性高湖久久久久久久久免费观看| 久久天躁狠狠躁夜夜2o2o| 伦理电影免费视频| 免费一级毛片在线播放高清视频 | 精品国产亚洲在线| 男女无遮挡免费网站观看| 免费女性裸体啪啪无遮挡网站| 国产成人免费观看mmmm| 窝窝影院91人妻| 日韩免费av在线播放| 午夜福利在线免费观看网站| 成人黄色视频免费在线看| 天堂8中文在线网| 日韩中文字幕欧美一区二区| 国产精品久久久av美女十八| 久久久国产精品麻豆| 一边摸一边抽搐一进一小说 | 日本vs欧美在线观看视频| 免费一级毛片在线播放高清视频 | 国产老妇伦熟女老妇高清| 日本五十路高清| 国产精品一区二区在线观看99| 国产欧美日韩综合在线一区二区| 国精品久久久久久国模美| 三上悠亚av全集在线观看| 啦啦啦 在线观看视频| 久久九九热精品免费| 国产成人av激情在线播放| 女同久久另类99精品国产91| 18禁国产床啪视频网站| 亚洲五月婷婷丁香| 久久热在线av| 午夜视频精品福利| 国产一区二区激情短视频| 成年人黄色毛片网站| 高潮久久久久久久久久久不卡| 亚洲午夜精品一区,二区,三区| 久久精品国产亚洲av香蕉五月 | 日韩视频一区二区在线观看| 精品国产一区二区久久| 午夜福利视频精品| 少妇粗大呻吟视频| 国产欧美亚洲国产| 一边摸一边抽搐一进一小说 | 夫妻午夜视频| 在线 av 中文字幕| 少妇被粗大的猛进出69影院| 欧美亚洲 丝袜 人妻 在线| 大片电影免费在线观看免费| 国产精品98久久久久久宅男小说| 精品一品国产午夜福利视频| 麻豆av在线久日| 青青草视频在线视频观看| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 99热网站在线观看| 午夜福利影视在线免费观看| 亚洲国产看品久久| 免费观看av网站的网址| 国产区一区二久久| av不卡在线播放| 日韩精品免费视频一区二区三区| 电影成人av| 国产成人av激情在线播放| 亚洲精品中文字幕在线视频| 十八禁网站免费在线| 国产97色在线日韩免费| 大片免费播放器 马上看| 老熟女久久久| 狠狠婷婷综合久久久久久88av| 精品欧美一区二区三区在线| 老熟妇乱子伦视频在线观看| 欧美激情 高清一区二区三区| 国产黄频视频在线观看| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| videosex国产| 国产精品亚洲av一区麻豆| 国产精品国产高清国产av | 搡老岳熟女国产| 老鸭窝网址在线观看| 欧美中文综合在线视频| 九色亚洲精品在线播放| 国产免费av片在线观看野外av| 大香蕉久久成人网| 啦啦啦在线免费观看视频4| 久久久久久久国产电影| 18禁美女被吸乳视频| 国产精品一区二区在线不卡| 亚洲精品国产区一区二| 久久久久精品人妻al黑| e午夜精品久久久久久久| av不卡在线播放| kizo精华| 午夜老司机福利片| a级毛片黄视频| 久久精品国产亚洲av香蕉五月 | 国产精品av久久久久免费| 男女高潮啪啪啪动态图| 激情视频va一区二区三区| 女性生殖器流出的白浆| av又黄又爽大尺度在线免费看| 麻豆乱淫一区二区| 在线天堂中文资源库| 亚洲成国产人片在线观看| 精品少妇黑人巨大在线播放| av福利片在线| 亚洲一区中文字幕在线| 日韩视频在线欧美| 香蕉国产在线看| 国产亚洲精品一区二区www | 人人妻,人人澡人人爽秒播| 久久久久久久精品吃奶| 9色porny在线观看| 午夜视频精品福利| 久久人妻熟女aⅴ| 伦理电影免费视频| 国产欧美日韩一区二区三区在线| 香蕉丝袜av| 国产精品一区二区免费欧美| 韩国精品一区二区三区| 男女边摸边吃奶| 国产1区2区3区精品| 亚洲成人免费av在线播放| 国产黄频视频在线观看| 国产在视频线精品| 考比视频在线观看| 国产麻豆69| 欧美日韩av久久| 十八禁网站免费在线| e午夜精品久久久久久久| 国产亚洲av高清不卡| 性高湖久久久久久久久免费观看| 久久久久视频综合| www.自偷自拍.com| 亚洲av欧美aⅴ国产| 久久国产亚洲av麻豆专区| 一区二区日韩欧美中文字幕| 久久香蕉激情| 色老头精品视频在线观看| 捣出白浆h1v1| 久久久久视频综合| 在线观看免费视频日本深夜| 女人久久www免费人成看片| 亚洲欧洲日产国产| 日韩一区二区三区影片| videosex国产| 丰满人妻熟妇乱又伦精品不卡| 99久久人妻综合| 一边摸一边抽搐一进一小说 | 黑人猛操日本美女一级片| 精品国产亚洲在线| 麻豆成人av在线观看| 国产伦理片在线播放av一区| 两人在一起打扑克的视频| 精品国产亚洲在线| 男女之事视频高清在线观看| 天堂8中文在线网| 丁香六月天网| 人人妻人人澡人人看| 97人妻天天添夜夜摸| 色94色欧美一区二区| 国产精品一区二区在线不卡| 97在线人人人人妻| 欧美日韩成人在线一区二区| 一本综合久久免费| 日本黄色日本黄色录像| 亚洲精品国产一区二区精华液| 另类精品久久| 久久免费观看电影| 老司机深夜福利视频在线观看| 亚洲精品美女久久久久99蜜臀| 夫妻午夜视频| 十八禁高潮呻吟视频| 香蕉国产在线看| 亚洲一区中文字幕在线| 亚洲人成77777在线视频| 久久人妻av系列| 国产精品久久久av美女十八| 亚洲天堂av无毛| 热re99久久精品国产66热6| 精品免费久久久久久久清纯 | 国产一区二区 视频在线| 欧美在线黄色| 久久久精品国产亚洲av高清涩受| 国产视频一区二区在线看| 性色av乱码一区二区三区2| 国精品久久久久久国模美| 国产又色又爽无遮挡免费看| 久久精品国产99精品国产亚洲性色 | 又大又爽又粗| 老司机亚洲免费影院| 五月开心婷婷网| 久久人妻av系列| 在线观看一区二区三区激情| 国产97色在线日韩免费| 久久久久久久大尺度免费视频| 最近最新免费中文字幕在线| 夜夜骑夜夜射夜夜干| 国产午夜精品久久久久久| 正在播放国产对白刺激| 日本欧美视频一区| 99热网站在线观看| 成人av一区二区三区在线看| 久久天躁狠狠躁夜夜2o2o| 久久精品国产综合久久久| 美女高潮到喷水免费观看| 制服人妻中文乱码| avwww免费| 最新在线观看一区二区三区| 在线亚洲精品国产二区图片欧美| 最新的欧美精品一区二区| 欧美性长视频在线观看| 久久中文字幕人妻熟女| 啦啦啦在线免费观看视频4| 国产欧美日韩一区二区三| 久久久久久亚洲精品国产蜜桃av| 涩涩av久久男人的天堂| 麻豆国产av国片精品| 国产精品亚洲一级av第二区| 纵有疾风起免费观看全集完整版| 香蕉久久夜色| 久久国产精品大桥未久av| 日日夜夜操网爽| 国产精品二区激情视频| 伦理电影免费视频| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩综合在线一区二区| 午夜精品国产一区二区电影| 亚洲熟女精品中文字幕| 国产日韩欧美视频二区| 中亚洲国语对白在线视频| 老司机亚洲免费影院| 日韩欧美免费精品| 黑人操中国人逼视频| 国产三级黄色录像| 亚洲 欧美一区二区三区| kizo精华| 大片电影免费在线观看免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲免费av在线视频| 日韩欧美国产一区二区入口| 午夜精品久久久久久毛片777| 久久久精品区二区三区| 午夜福利在线观看吧| 九色亚洲精品在线播放| 高清黄色对白视频在线免费看| 亚洲自偷自拍图片 自拍| 黑丝袜美女国产一区| 午夜福利视频精品| 狠狠婷婷综合久久久久久88av| 亚洲欧美日韩高清在线视频 | 国产单亲对白刺激| 国产亚洲精品第一综合不卡| av片东京热男人的天堂| www.自偷自拍.com| 菩萨蛮人人尽说江南好唐韦庄| 久久人人97超碰香蕉20202| 99香蕉大伊视频| 妹子高潮喷水视频| 视频在线观看一区二区三区| 久久久精品区二区三区| 成人手机av| 日韩一区二区三区影片| 99riav亚洲国产免费| 国产视频一区二区在线看| 午夜福利在线观看吧| 亚洲精品成人av观看孕妇| 午夜福利在线观看吧| 侵犯人妻中文字幕一二三四区| 1024视频免费在线观看| 欧美精品亚洲一区二区| 极品教师在线免费播放| 中文字幕精品免费在线观看视频| 精品免费久久久久久久清纯 | 国产1区2区3区精品| 黄色视频不卡| 一区二区三区国产精品乱码| 日韩熟女老妇一区二区性免费视频| 精品高清国产在线一区| 日韩视频一区二区在线观看| 亚洲av成人一区二区三| 国产欧美亚洲国产| 美女福利国产在线| 久久人妻熟女aⅴ| 色婷婷久久久亚洲欧美| 在线看a的网站| 成年版毛片免费区| 一级片免费观看大全| 少妇 在线观看| 91国产中文字幕| 超碰97精品在线观看| 中文字幕人妻熟女乱码| 久久人人97超碰香蕉20202| 国产精品秋霞免费鲁丝片| 12—13女人毛片做爰片一| 色婷婷av一区二区三区视频| 国产精品麻豆人妻色哟哟久久| 成人国语在线视频| 亚洲第一av免费看| 大片电影免费在线观看免费| 久久精品aⅴ一区二区三区四区| 久久午夜综合久久蜜桃| 一二三四社区在线视频社区8| 精品少妇久久久久久888优播| 国产三级黄色录像| 色综合婷婷激情| 国产精品av久久久久免费| 欧美乱码精品一区二区三区| 国产免费福利视频在线观看| 国产亚洲av高清不卡| 亚洲av国产av综合av卡| 精品国产超薄肉色丝袜足j| 欧美日韩一级在线毛片| 国产深夜福利视频在线观看| 亚洲免费av在线视频| 精品视频人人做人人爽| 欧美大码av| 国产一区二区三区综合在线观看| 国产成人免费无遮挡视频|