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

    鈉工藝間鈉霧火模型適用性分析

    2022-03-22 08:26:28俞之桐曹學(xué)武
    核技術(shù) 2022年3期
    關(guān)鍵詞:消耗量液滴預(yù)測值

    俞之桐 曹學(xué)武

    (上海交通大學(xué)機(jī)械與動力工程學(xué)院 上海 200240)

    針對鈉冷快堆安全研究,除了關(guān)注無保護(hù)失流等設(shè)計(jì)基準(zhǔn)事故外[1],鈉冷快堆一回路和二回路系統(tǒng)中鈉的泄漏[2]與燃燒也是核電廠的潛在威脅。鈉霧火是泄漏的鈉液滴在拋射狀態(tài)下發(fā)生的燃燒,燃燒速率大,難以撲滅[3],會破壞鈉工藝間安全設(shè)備和建 筑 結(jié) 構(gòu)[4]?,F(xiàn) 象 識 別 與 排 名 表(Phenomena Identification and Ranking Table,PIRT)認(rèn)為鈉霧火現(xiàn)象“非常重要”[5]。因此,評估鈉霧火模型適用性對研究鈉冷快堆鈉霧火事故后果具有重要意義。

    鈉霧火實(shí)驗(yàn)可以用于鈉霧火模型的建立和程序的開發(fā)驗(yàn)證。 其中,原子國際(Atomics International,AI)組織的J系列實(shí)驗(yàn)[6]、法國的小空間鈉霧火實(shí)驗(yàn)[7]、德國的FS 系列鈉霧火實(shí)驗(yàn)[8]以及美國桑迪亞實(shí)驗(yàn)室(Sandia National Laboratories,SNL)的T系列實(shí)驗(yàn)[9]等廣泛應(yīng)用于鈉霧火程序的開發(fā)和驗(yàn)證。鈉霧火程序可以分為兩類:計(jì)算流體力學(xué)方法,如日本的鈉泄漏火災(zāi)與氣溶膠分析程序[10](SOdium Leak Fire and Aerosol analysis code System,SOLFAS);集總參數(shù)法的鈉霧火模型,如基于Tsai 模型的開源鈉燃燒程序(Sodium Combustion,NACOM)與在Tsai 模型基礎(chǔ)上改進(jìn)的CONTAIN 程序鈉冷快堆版本(Liquid Metal Reactor version for CONTAIN code,CONTAIN-LMR)[11]以及多節(jié)點(diǎn)系統(tǒng)鈉火現(xiàn)象程序(Sodium fire Phenomenology IN multi-Cell System,SPHINCS)[12]等。Tsai 模型計(jì)算流程復(fù)雜,無法與鈉霧火期間其他重要現(xiàn)象,如鈉氣溶膠行為、鈉池火等耦合,對AI的J1與J2實(shí)驗(yàn)鈉消耗量模擬值比實(shí)驗(yàn)值高估30%~55%[13]。國內(nèi)在Tsai模型的基礎(chǔ)上對鈉霧火現(xiàn)象開展了較多模擬。張斌等[14]在Tsai鈉霧火模型基礎(chǔ)上加入流場模型,給出空間氣體溫度、氧氣濃度等物理量的三維分布。喻宏[15]在Tsai 模型的基礎(chǔ)上開發(fā)SSPRAY 鈉霧火程序,利用AI 的J3、J4 實(shí)驗(yàn)驗(yàn)證程序適用性。

    Tsai模型假設(shè)噴霧燃燒速率是鈉液滴燃燒速率之和,是固定粒徑鈉液滴隨時(shí)間變化的模型,不能反映下落過程中鈉液滴間相互碰撞融合作用以及噴霧區(qū)內(nèi)鈉液滴未與空氣接觸的現(xiàn)象[13],導(dǎo)致鈉液滴燃燒以及運(yùn)動預(yù)測過于保守[16],從而大大高估了鈉噴霧的燃燒速率,如Nashine 等[16]利用KAERI 的SP1與SP4鈉霧火實(shí)驗(yàn)對Tsai模型的驗(yàn)證發(fā)現(xiàn)氣溫與氣壓升高模擬值比實(shí)驗(yàn)值高了115%~1 550%。

    本文對Tsai 模型中鈉液滴運(yùn)動模型、噴霧燃燒模型進(jìn)行了修改,并討論了改進(jìn)過程對鈉霧火重要物理量(氣溫、氣壓峰值與鈉消耗量等)預(yù)測值的影響。針對SNL 的T3 鈉霧火實(shí)驗(yàn)以及AI 的J1-J4 鈉霧火實(shí)驗(yàn)進(jìn)行建模分析,通過對比模型預(yù)測結(jié)果與實(shí)驗(yàn)測量結(jié)果,評估鈉霧火階段和初始氧濃度等因素對改進(jìn)鈉霧火模型適用性的影響。

    1 鈉霧火模型

    1.1 液滴模型

    液滴模型包括鈉液滴尺寸和鈉液滴運(yùn)動兩部分。

    鈉液滴尺寸由平均直徑確定,按照Nukiyama-Tanasawa關(guān)系式將鈉液滴分為11組[13],如式(1):

    式中:D為鈉液滴直徑,m;Dˉ為鈉液滴平均直徑,m;R為所有直徑小于D的鈉液滴在噴霧中的體積份額。

    Tsai模型雖然能夠表征鈉液滴運(yùn)動的速度變化過程,但由于未考慮噴霧區(qū)內(nèi)鈉液滴與空氣非接觸以及液滴相互作用,會高估實(shí)際工況下鈉的燃燒速率和下落時(shí)間。因此,采用單一下落速度模型,將下落速度設(shè)定為鈉液滴落至地面的最終速度。改進(jìn)模型減少鈉液滴燃燒時(shí)間計(jì)算值,改進(jìn)了空間氣溫、氣壓峰值與鈉消耗量預(yù)測值的計(jì)算效果,并解決了Tsai模型開始下落時(shí)鈉液滴高加速度引起數(shù)值不穩(wěn)定性[13]的問題。

    1.2 單液滴燃燒模型

    自由下落鈉液滴燃燒速率基于D2規(guī)律,即鈉液滴直徑的平方隨下落時(shí)間線性減少,如式(2):

    式中:Di為鈉液滴初始直徑,m;t為鈉液滴下落時(shí)間,s;K為蒸發(fā)系數(shù),如式(3):

    式中:kg為氣體熱導(dǎo)率,W·(m·K)?1;Cp,g為氣體定壓比熱容,J·(kg·K)?1;ρNa為鈉液滴密度,kg·m?3;B為傳遞系數(shù),如式(4):

    式中:hfg為鈉液滴蒸發(fā)熱,J·kg?1;Tg為氣空間溫度,℃;TNa為鈉液滴溫度,℃;Hc為鈉燃燒熱,J·kg?1;Y為氧氣摩爾份額;i為化學(xué)計(jì)量率。

    靜態(tài)鈉液滴的質(zhì)量燃燒速率m?f,0如式(5):

    對于自由下落的鈉液滴的質(zhì)量燃燒速率m?f,考慮對流換熱的影響,引入增殖因子,如式(6):

    式中:Cf為經(jīng)驗(yàn)常數(shù),取0.3;Re為雷諾數(shù);Pr為Prandtl數(shù)。

    1.3 噴霧燃燒模型

    假定整個鈉工藝間的氣空間溫度壓力均勻分布??紤]下落過程鈉液滴尺寸不斷減少,采用Lagrange 方法跟蹤各組液滴從頂部到地板的燃燒過程。

    Tsai 模型中,對下落過程中不同直徑鈉液滴在不同高度處的總消耗量和放熱量進(jìn)行積分,積分步驟由時(shí)間步長確定,可能需要劃分出上百個積分步驟,計(jì)算復(fù)雜。改進(jìn)模型不考慮鈉液滴下落中的燃燒過程,將鈉液滴的積分放熱量直接加入氣空間,鈉燃燒對氣空間瞬發(fā)影響會使空間氣溫與氣壓達(dá)到峰值所需時(shí)間的預(yù)測值與噴霧持續(xù)時(shí)間相同,但是由于燃燒總放熱量確定,不會對鈉霧火導(dǎo)致的空間氣溫與壓力峰值以及鈉消耗量等預(yù)測值產(chǎn)生影響。噴霧燃燒模型的改進(jìn)在幾乎不影響鈉霧火重要物理量預(yù)測值的情況下,減少了模型復(fù)雜程度,方便鈉霧火模型與其他模型耦合。

    模型假設(shè)燃燒液滴無相互作用,各積分步驟中鈉噴霧的質(zhì)量燃燒速率如式(7):

    2 模型適用性驗(yàn)證

    2.1 鈉霧火階段對模型適用性影響

    SNL為了研究先進(jìn)反應(yīng)堆的金屬火災(zāi)以及其后果進(jìn)行了一系列的鈉火實(shí)驗(yàn)。其中,T3實(shí)驗(yàn)是在體積99 m3的Surtsey大尺寸鋼制裝置中進(jìn)行的鈉霧火實(shí)驗(yàn)。

    利用T3 實(shí)驗(yàn)對改進(jìn)模型進(jìn)行驗(yàn)證。實(shí)驗(yàn)中鈉噴霧的泄露高度為5.3 m,噴霧中的鈉液滴溫度為200 ℃;Surtsey 實(shí)驗(yàn)裝置的直徑為3.6 m,高為10.3 m,體積為99 m3,壁厚為0.01 m;氣空間的初始?xì)鉁貫?88.15 K,初始?xì)鈮簽?01.3 kPa,初始氧濃度為23.3%。

    實(shí)驗(yàn)利用光測高溫計(jì)測量發(fā)生鈉噴霧燃燒時(shí)噴嘴下方2 m 處鈉噴霧溫度[9]。由光測高溫計(jì)得到的鈉噴霧溫度峰值對應(yīng)的時(shí)間與實(shí)驗(yàn)中攝像頭拍攝到的兩次黃色閃光出現(xiàn)的時(shí)間相吻合,可以證明實(shí)驗(yàn)中鈉霧火分兩次進(jìn)行,持續(xù)時(shí)間為14 s[18]。第一次噴霧燃燒時(shí)間段為0~3 s,第二次噴霧燃燒時(shí)間段為9~14 s[19]。鈉霧火模擬參數(shù)如表1所示。

    表1 T3實(shí)驗(yàn)?zāi)P蛥?shù)Table 1 Parameter of T3 experiment

    實(shí)驗(yàn)在距離地板0.5~6.5 m共放置7個熱電偶用于測量氣空間氣體溫度[9],取7個熱電偶測量氣溫平均值作為氣空間氣體的平均溫度[19]。T3 鈉霧火實(shí)驗(yàn)早期階段(40 s)和長期階段(300 s)鈉霧火模型氣空間壓力、平均溫度預(yù)計(jì)值與實(shí)驗(yàn)測量值的對比分別如圖1、2 所示。其中,兩次鈉霧火與一次鈉霧火的模擬方式考慮了鈉噴霧總質(zhì)量相同,鈉霧火次數(shù)和鈉液滴質(zhì)量流量不同。

    從圖1可以看出,在早期工況下,改進(jìn)模型的氣壓預(yù)測結(jié)果與實(shí)驗(yàn)結(jié)果擬合良好;氣溫峰值預(yù)測值比測量值高了15 ℃,這可能是由于在發(fā)生鈉霧火時(shí),氣體與壁面的傳熱系數(shù)不斷減少[17],而模型中將傳熱系數(shù)視為定值,可能對實(shí)驗(yàn)初期傳熱系數(shù)存在低估。另外,由于實(shí)驗(yàn)裝置內(nèi)的氣體有流動換熱的過程,不同高度達(dá)到氣溫峰值的時(shí)間不同[9],導(dǎo)致平均氣溫峰值時(shí)間滯后于氣壓峰值時(shí)間約3 s。

    圖1 鈉霧火模型評估T3鈉霧火實(shí)驗(yàn)早期階段(前40 s)氣空間狀況 (a)氣壓,(b)平均氣溫Fig.1 Comparison of atmosphere condition between T3 experimental data and prediction data of sodium spray fire model in short term(the first 40 s) (a)Gas pressure,(b)Average temperature of atmosphere

    從圖2 可以看出,前100 s 模型與實(shí)驗(yàn)擬合較好;而200 s 后,改進(jìn)模型氣空間平均氣溫升高預(yù)測值、氣壓升高預(yù)測值分別在實(shí)驗(yàn)值的60%、72%以下,這是由于未燃燒的鈉液滴下落至地板形成鈉池,鈉池通過鈉池火等方式向氣空間放熱,說明事故長期工況的氣空間溫度壓力評估需要考慮未燃燒的鈉的影響。

    圖2 鈉霧火模型評估T3實(shí)驗(yàn)長期階段(前300 s)氣空間狀況 (a)氣壓,(b)平均氣溫Fig.2 Comparison of atmosphere condition between T3 experimental data and prediction data of sodium spray fire model in long term(the first 300 s) (a)Gas pressure,(b)Average temperature of atmosphere

    改進(jìn)模型中,鈉噴霧總質(zhì)量相同時(shí),不同的鈉液滴質(zhì)量流量對實(shí)驗(yàn)早期階段氣溫、氣壓峰值預(yù)測值偏差很小,約為2%,對實(shí)驗(yàn)長期階段氣溫、氣壓預(yù)測值幾乎無影響。

    2.2 初始氧濃度對模型適用性影響

    AI 在55.6 m3測試容器(Large Test Vessel,LTV)實(shí)驗(yàn)裝置中進(jìn)行了鈉霧火系列實(shí)驗(yàn),實(shí)驗(yàn)中鈉向上噴出至擋板形成雨?duì)顕婌F下落[13]。利用改進(jìn)后的鈉霧火模型對初始氧濃度為21%的J1、J2以及初始氧濃度小于2%的J3、J4 實(shí)驗(yàn)進(jìn)行模擬。實(shí)驗(yàn)參數(shù)如表2所示。

    表2 AI鈉霧火實(shí)驗(yàn)參數(shù)Table 2 Parameter of AI spray fire experiment

    改進(jìn)模型預(yù)測值與實(shí)驗(yàn)值的比較如表3 所示。由表3 可以看出,在模擬初始氧濃度為21%的J1 和J2 實(shí)驗(yàn)時(shí),改進(jìn)模型中鈉消耗預(yù)測值分別為實(shí)驗(yàn)值的118%和105%,氣空間最大升壓分別為實(shí)驗(yàn)值的105%和99%。改進(jìn)模型能夠很好地預(yù)測鈉泄漏事故下鈉工藝間中鈉噴霧在空氣中燃燒產(chǎn)生的氣壓峰值以及鈉消耗量。改進(jìn)模型將鈉液滴落地速度視為下落速度,液滴下落時(shí)間縮短,降低了鈉消耗量的預(yù)測值,減小了由于Tsai 模型高估鈉燃燒速率而造成高估鈉消耗量的預(yù)測偏差。

    表3 AI鈉霧火實(shí)驗(yàn)測量值與預(yù)測值比較Table 3 Comparison between experimental data and computational results for AI jet tests

    在模擬初始氧濃度分別為1.5%、0%的J3 和J4實(shí)驗(yàn)時(shí),鈉霧火模型的氣空間最大升壓預(yù)測值分別為實(shí)驗(yàn)值的為80%和71.5%。模擬低氧氣濃度的J3、J4 工況時(shí),鈉噴霧燃燒引起的氣溫氣壓升高較小,AI實(shí)驗(yàn)氣空間連接的鈉液滴產(chǎn)生裝置對氣空間氣溫氣壓的影響不可忽略,否則會導(dǎo)致模型氣溫氣壓預(yù)測值偏低。

    AI 實(shí)驗(yàn)中采用的鈉液滴產(chǎn)生裝置壓力恒為180 kPa,連接鈉液滴產(chǎn)生裝置和氣空間的噴管截面積 為5×10?4m2[13]???慮500 ℃下 鈉 的 密 度 為829 kg·m?3,由噴射速度和鈉質(zhì)量流量可以分別計(jì)算得到總體積流量和鈉體積流量。在J3、J4 實(shí)驗(yàn)中鈉體積流量分別約占總體積流量的60%和80%,所以分別設(shè)置鈉液滴與氣空間的流動路徑截面積為2×10?4m2、1×10?4m2。

    考慮鈉液滴產(chǎn)生裝置后,鈉霧火模型對J3、J4實(shí)驗(yàn)氣空間氣壓模擬結(jié)果如圖3所示。由圖3可得,考慮鈉液滴產(chǎn)生裝置的影響可以將鈉霧火模型對J3、J4實(shí)驗(yàn)低初始氧濃度工況的氣壓預(yù)測值偏差分別由?20.0%和?28.5%降低至?4.1%和?7.5%。

    圖3 鈉霧火模型氣空間氣壓升高模擬結(jié)果 (a)J3實(shí)驗(yàn),(b)J4實(shí)驗(yàn)Fig.3 Simulation results of atmosphere pressure rise by using sodium spray fire model for J3 experiment(a)and J4 experiment(b)

    3 結(jié)語

    通過對SNL 實(shí)驗(yàn)室開展的T3 實(shí)驗(yàn)以及AI 的J1-J4鈉霧火實(shí)驗(yàn)的模擬,評估了改進(jìn)后Tsai模型對鈉霧火不同階段以及不同初始氧濃度的鈉霧火現(xiàn)象的適用性。結(jié)果表明:改進(jìn)鈉霧火模型可以用于評估鈉工藝間鈉噴霧在空氣中燃燒引起的升溫與升壓現(xiàn)象;在模擬低初始氧濃度鈉霧火工況時(shí),最大升壓預(yù)測值與實(shí)驗(yàn)值偏差為?28%~?20%,考慮鈉液滴產(chǎn)生裝置對氣空間的影響時(shí),最大升壓預(yù)測值與實(shí)驗(yàn)值偏差降低至?7.5%~?4.1%。單一下落速度模型在一定程度上可以降低Tsai 模型高估鈉消耗量的程度,使鈉消耗量模擬值與實(shí)驗(yàn)值的偏差為5%~18%。但是,目前模型不能預(yù)測未消耗的鈉在長期工況中的池火對空間氣溫與壓力的影響,后續(xù)可進(jìn)一步耦合鈉池火模型。

    作者貢獻(xiàn)聲明俞之桐:模型研究,實(shí)施模擬研究,整理與分析數(shù)據(jù),起草文章;曹學(xué)武:提出研究思路、研究方法,對文章的知識性內(nèi)容作批評性審閱。

    猜你喜歡
    消耗量液滴預(yù)測值
    IMF上調(diào)今年全球經(jīng)濟(jì)增長預(yù)期
    企業(yè)界(2024年8期)2024-07-05 10:59:04
    路基石方爆破降低炸藥消耗量研究
    加拿大農(nóng)業(yè)部下調(diào)2021/22年度油菜籽和小麥產(chǎn)量預(yù)測值
    ±800kV直流輸電工程合成電場夏季實(shí)測值與預(yù)測值比對分析
    液滴間相互碰撞融合與破碎的實(shí)驗(yàn)研究
    噴淋液滴在空氣環(huán)境下的運(yùn)動特性
    法電再次修訂2020年核發(fā)電量預(yù)測值
    國外核新聞(2020年8期)2020-03-14 02:09:19
    有機(jī)化學(xué)反應(yīng)中試劑最大消耗量問題探析
    《輕型商用車輛燃料消耗量限值》強(qiáng)制性國家標(biāo)準(zhǔn)發(fā)布
    邁騰1.8TSI車燃油消耗量大
    欧美激情在线99| av.在线天堂| 青春草视频在线免费观看| 国产久久久一区二区三区| 久久久国产成人免费| 亚洲aⅴ乱码一区二区在线播放| 色综合站精品国产| 亚洲美女视频黄频| 欧美日本亚洲视频在线播放| 成人毛片a级毛片在线播放| 日本熟妇午夜| 精品少妇黑人巨大在线播放 | 色吧在线观看| 欧美bdsm另类| 2021天堂中文幕一二区在线观| 国产爱豆传媒在线观看| 日本熟妇午夜| 天天一区二区日本电影三级| 亚洲国产精品合色在线| 精品不卡国产一区二区三区| 欧美人与善性xxx| 国产av麻豆久久久久久久| 亚洲第一电影网av| 国产精品蜜桃在线观看 | 内地一区二区视频在线| 最后的刺客免费高清国语| 久久99蜜桃精品久久| 亚洲欧美成人精品一区二区| 夫妻性生交免费视频一级片| 日本撒尿小便嘘嘘汇集6| 久久精品影院6| 级片在线观看| 一区二区三区免费毛片| 淫秽高清视频在线观看| 色吧在线观看| 婷婷色综合大香蕉| 国内少妇人妻偷人精品xxx网站| 久久鲁丝午夜福利片| av专区在线播放| 级片在线观看| 色尼玛亚洲综合影院| 性欧美人与动物交配| 亚洲av免费高清在线观看| 亚洲国产精品成人综合色| 欧美性猛交黑人性爽| 国产精品久久久久久精品电影| 搡女人真爽免费视频火全软件| 国产精品久久久久久亚洲av鲁大| 国产熟女欧美一区二区| 天天一区二区日本电影三级| 人人妻人人看人人澡| 成人欧美大片| 精品久久久噜噜| 在线播放国产精品三级| 日本黄色片子视频| 人妻少妇偷人精品九色| 一本久久精品| 国产成人a∨麻豆精品| 日本av手机在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品欧美国产一区二区三| 最近视频中文字幕2019在线8| 26uuu在线亚洲综合色| 成人三级黄色视频| 日韩一本色道免费dvd| 国产午夜福利久久久久久| 亚洲第一区二区三区不卡| 变态另类成人亚洲欧美熟女| 在现免费观看毛片| 99久久无色码亚洲精品果冻| 97热精品久久久久久| 精品久久久久久久久av| 只有这里有精品99| 日本免费一区二区三区高清不卡| 美女脱内裤让男人舔精品视频 | 亚洲精品乱码久久久久久按摩| 自拍偷自拍亚洲精品老妇| 五月伊人婷婷丁香| 看非洲黑人一级黄片| 97人妻精品一区二区三区麻豆| 国产一区二区三区av在线 | 久久人人爽人人片av| 国产美女午夜福利| 51国产日韩欧美| 最近2019中文字幕mv第一页| 插阴视频在线观看视频| 99久国产av精品国产电影| 亚洲国产色片| 精品无人区乱码1区二区| 亚洲精品456在线播放app| 人妻久久中文字幕网| 国产又黄又爽又无遮挡在线| 国产av一区在线观看免费| 亚洲精品乱码久久久v下载方式| 高清毛片免费看| 国产成人a区在线观看| 99热精品在线国产| 97超碰精品成人国产| 不卡视频在线观看欧美| 在线观看66精品国产| 波野结衣二区三区在线| 久久久久久久久久久免费av| 久久这里有精品视频免费| 乱码一卡2卡4卡精品| 精品日产1卡2卡| 夜夜夜夜夜久久久久| 最近2019中文字幕mv第一页| 插阴视频在线观看视频| 能在线免费看毛片的网站| 少妇猛男粗大的猛烈进出视频 | АⅤ资源中文在线天堂| 国产精品福利在线免费观看| 久久久久久九九精品二区国产| 亚洲精品自拍成人| 国产av在哪里看| 成人av在线播放网站| 插阴视频在线观看视频| 美女内射精品一级片tv| 久久午夜福利片| 欧美zozozo另类| 一本久久中文字幕| 国产美女午夜福利| 晚上一个人看的免费电影| 特大巨黑吊av在线直播| 精品99又大又爽又粗少妇毛片| 最近最新中文字幕大全电影3| 亚洲丝袜综合中文字幕| 欧美日韩国产亚洲二区| 免费无遮挡裸体视频| 天堂中文最新版在线下载 | 欧美潮喷喷水| 亚洲国产精品合色在线| 成人特级黄色片久久久久久久| 两个人的视频大全免费| 亚洲欧美精品自产自拍| 边亲边吃奶的免费视频| 久久99精品国语久久久| 久久久成人免费电影| 又粗又硬又长又爽又黄的视频 | 国产精品野战在线观看| 99久国产av精品国产电影| a级毛色黄片| 国产乱人偷精品视频| 麻豆一二三区av精品| 日本熟妇午夜| 国产成人精品久久久久久| 自拍偷自拍亚洲精品老妇| 男女做爰动态图高潮gif福利片| 长腿黑丝高跟| 69人妻影院| 欧美成人一区二区免费高清观看| 国产精品嫩草影院av在线观看| 国产亚洲欧美98| 国产av一区在线观看免费| 精品久久久久久久久久免费视频| 国产伦精品一区二区三区四那| 欧美激情国产日韩精品一区| 三级男女做爰猛烈吃奶摸视频| 三级男女做爰猛烈吃奶摸视频| 天天躁夜夜躁狠狠久久av| avwww免费| 国产黄a三级三级三级人| 亚洲精品国产成人久久av| av在线播放精品| av女优亚洲男人天堂| 亚洲av电影不卡..在线观看| 国产综合懂色| 国产亚洲av片在线观看秒播厂 | 12—13女人毛片做爰片一| 日日啪夜夜撸| 午夜a级毛片| 亚洲国产欧美人成| 深爱激情五月婷婷| 日韩高清综合在线| 秋霞在线观看毛片| 国产视频内射| 亚洲不卡免费看| 免费一级毛片在线播放高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 少妇高潮的动态图| 我的老师免费观看完整版| 少妇人妻一区二区三区视频| 在线免费观看的www视频| 最近视频中文字幕2019在线8| 男插女下体视频免费在线播放| 男人狂女人下面高潮的视频| 尤物成人国产欧美一区二区三区| 别揉我奶头 嗯啊视频| 亚洲欧美成人精品一区二区| 日本-黄色视频高清免费观看| 国产乱人视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品一区二区三区四区免费观看| 精品一区二区三区人妻视频| 亚洲欧美精品专区久久| 免费电影在线观看免费观看| 99九九线精品视频在线观看视频| 亚洲欧美清纯卡通| 蜜桃亚洲精品一区二区三区| 又黄又爽又刺激的免费视频.| 97超碰精品成人国产| 夜夜看夜夜爽夜夜摸| 村上凉子中文字幕在线| www.色视频.com| 超碰av人人做人人爽久久| 男人和女人高潮做爰伦理| 欧美丝袜亚洲另类| 亚洲天堂国产精品一区在线| 国产精品不卡视频一区二区| 亚洲中文字幕日韩| 3wmmmm亚洲av在线观看| av又黄又爽大尺度在线免费看 | 精品人妻偷拍中文字幕| 18禁黄网站禁片免费观看直播| 亚洲精品色激情综合| 好男人视频免费观看在线| 久久久久性生活片| 国产探花在线观看一区二区| h日本视频在线播放| 蜜臀久久99精品久久宅男| 五月玫瑰六月丁香| 国产精品女同一区二区软件| 99久久精品国产国产毛片| 别揉我奶头 嗯啊视频| 在线国产一区二区在线| 久久精品91蜜桃| 五月伊人婷婷丁香| 欧美一区二区国产精品久久精品| 色播亚洲综合网| 亚洲欧美中文字幕日韩二区| 一本精品99久久精品77| 欧美人与善性xxx| 国产女主播在线喷水免费视频网站 | 国产亚洲5aaaaa淫片| 亚州av有码| 日日摸夜夜添夜夜爱| 国产黄色视频一区二区在线观看 | 婷婷色av中文字幕| 国产在线精品亚洲第一网站| 欧美人与善性xxx| 亚洲精品亚洲一区二区| 亚洲成人中文字幕在线播放| 久久精品国产自在天天线| 国产一区二区激情短视频| 免费观看的影片在线观看| 亚洲图色成人| 久久久久国产网址| 蜜臀久久99精品久久宅男| 欧美精品国产亚洲| 亚洲高清免费不卡视频| 又爽又黄无遮挡网站| 久久这里有精品视频免费| 国产高清不卡午夜福利| 久久久成人免费电影| 狠狠狠狠99中文字幕| 非洲黑人性xxxx精品又粗又长| 成人欧美大片| 亚洲第一区二区三区不卡| 国产不卡一卡二| 美女xxoo啪啪120秒动态图| 国产日本99.免费观看| 亚洲在线自拍视频| 婷婷精品国产亚洲av| 少妇被粗大猛烈的视频| or卡值多少钱| 在线免费观看不下载黄p国产| 成人综合一区亚洲| 欧美日韩综合久久久久久| av国产免费在线观看| 国产精品爽爽va在线观看网站| 男人和女人高潮做爰伦理| 中文字幕人妻熟人妻熟丝袜美| 欧美不卡视频在线免费观看| 美女大奶头视频| 超碰av人人做人人爽久久| 日韩 亚洲 欧美在线| 亚洲欧美日韩高清在线视频| 在线天堂最新版资源| 日韩制服骚丝袜av| 男女视频在线观看网站免费| 精品日产1卡2卡| 91狼人影院| 国产亚洲5aaaaa淫片| 国产黄a三级三级三级人| 能在线免费看毛片的网站| 午夜福利成人在线免费观看| 欧美+日韩+精品| 日本免费一区二区三区高清不卡| 久久久成人免费电影| 日日摸夜夜添夜夜爱| 亚洲图色成人| 国产高清不卡午夜福利| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品久久久久久婷婷小说 | 国产极品精品免费视频能看的| 婷婷六月久久综合丁香| 校园人妻丝袜中文字幕| 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 99热网站在线观看| 在线观看美女被高潮喷水网站| 国产真实乱freesex| 成人美女网站在线观看视频| 白带黄色成豆腐渣| 最好的美女福利视频网| 青春草亚洲视频在线观看| 麻豆乱淫一区二区| 国产激情偷乱视频一区二区| 亚洲精品乱码久久久久久按摩| 国产一区二区在线av高清观看| 一本一本综合久久| 亚洲在线观看片| 两个人的视频大全免费| 国产成人精品婷婷| 一本久久中文字幕| 日韩欧美国产在线观看| av福利片在线观看| 国产精品免费一区二区三区在线| 禁无遮挡网站| 久久久a久久爽久久v久久| 不卡一级毛片| 97在线视频观看| 午夜激情福利司机影院| 小蜜桃在线观看免费完整版高清| 亚洲成人久久性| 久久久久久久午夜电影| av女优亚洲男人天堂| 一级毛片aaaaaa免费看小| 国产亚洲精品久久久久久毛片| 久久久久久久久大av| 亚洲精品影视一区二区三区av| 观看免费一级毛片| 1000部很黄的大片| 天堂√8在线中文| 国产黄色视频一区二区在线观看 | 在线播放国产精品三级| 免费av不卡在线播放| 久久久久久伊人网av| 能在线免费看毛片的网站| 午夜精品国产一区二区电影 | 日本在线视频免费播放| 一个人看视频在线观看www免费| 亚洲中文字幕一区二区三区有码在线看| 成年女人看的毛片在线观看| 久久精品国产鲁丝片午夜精品| 我要搜黄色片| 午夜福利高清视频| 久久久久久久久久久丰满| 如何舔出高潮| 在线播放无遮挡| ponron亚洲| 天堂√8在线中文| 国产麻豆成人av免费视频| 欧美高清成人免费视频www| 国产精品无大码| av视频在线观看入口| 最近最新中文字幕大全电影3| 身体一侧抽搐| 国产成人福利小说| 美女被艹到高潮喷水动态| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利视频1000在线观看| 国产成人a区在线观看| 91久久精品电影网| 美女xxoo啪啪120秒动态图| 中出人妻视频一区二区| 大香蕉久久网| 国产蜜桃级精品一区二区三区| 禁无遮挡网站| 亚洲在线观看片| 12—13女人毛片做爰片一| 两个人视频免费观看高清| 成熟少妇高潮喷水视频| 亚洲成人av在线免费| АⅤ资源中文在线天堂| 少妇人妻精品综合一区二区 | 国产精品综合久久久久久久免费| 男女边吃奶边做爰视频| 成年免费大片在线观看| 免费无遮挡裸体视频| 99热只有精品国产| 国产黄色视频一区二区在线观看 | 97超碰精品成人国产| 在线观看一区二区三区| 春色校园在线视频观看| 久久这里有精品视频免费| 日本色播在线视频| 日日干狠狠操夜夜爽| 成年av动漫网址| 我的老师免费观看完整版| 国产精品爽爽va在线观看网站| 特大巨黑吊av在线直播| 高清午夜精品一区二区三区 | 又黄又爽又刺激的免费视频.| 日韩av不卡免费在线播放| 日韩av在线大香蕉| 午夜精品一区二区三区免费看| 国产极品天堂在线| 国产成人a∨麻豆精品| 精品久久久久久久久久免费视频| 免费看a级黄色片| 能在线免费看毛片的网站| 青春草视频在线免费观看| av免费观看日本| 久久国内精品自在自线图片| 可以在线观看毛片的网站| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 可以在线观看毛片的网站| 亚洲国产精品sss在线观看| 精品久久久久久久末码| 免费一级毛片在线播放高清视频| 日韩欧美精品v在线| 男人的好看免费观看在线视频| 97超碰精品成人国产| 免费av毛片视频| 秋霞在线观看毛片| 亚洲色图av天堂| 精品午夜福利在线看| 成熟少妇高潮喷水视频| 成人毛片a级毛片在线播放| 日韩成人av中文字幕在线观看| 99热这里只有是精品50| 亚洲自拍偷在线| 美女内射精品一级片tv| 国产成人a区在线观看| 91aial.com中文字幕在线观看| 男人的好看免费观看在线视频| 97超碰精品成人国产| 久久精品国产清高在天天线| 熟女人妻精品中文字幕| 99久久九九国产精品国产免费| 国产精品三级大全| 久久精品国产亚洲av天美| 又粗又爽又猛毛片免费看| 少妇的逼水好多| 熟女电影av网| 校园人妻丝袜中文字幕| 国产色爽女视频免费观看| 97人妻精品一区二区三区麻豆| 精品一区二区三区视频在线| 婷婷亚洲欧美| a级毛片a级免费在线| 中文在线观看免费www的网站| 欧美日本亚洲视频在线播放| 国产老妇伦熟女老妇高清| 老熟妇乱子伦视频在线观看| 国国产精品蜜臀av免费| 99久久人妻综合| 国产三级在线视频| 在线免费十八禁| 三级男女做爰猛烈吃奶摸视频| 青青草视频在线视频观看| 99热这里只有是精品50| 国产不卡一卡二| 中文字幕精品亚洲无线码一区| 美女脱内裤让男人舔精品视频 | 人妻少妇偷人精品九色| 国产黄片美女视频| 免费看美女性在线毛片视频| 天堂网av新在线| 热99在线观看视频| 久久国内精品自在自线图片| 少妇高潮的动态图| 国产精品精品国产色婷婷| 国产色爽女视频免费观看| 18禁黄网站禁片免费观看直播| 成人性生交大片免费视频hd| 三级经典国产精品| 免费大片18禁| 免费观看精品视频网站| 26uuu在线亚洲综合色| 国产精品一二三区在线看| 色综合色国产| 色吧在线观看| 免费看日本二区| 成年女人看的毛片在线观看| 99在线人妻在线中文字幕| 欧美另类亚洲清纯唯美| 日本五十路高清| 亚洲欧美日韩高清专用| 波野结衣二区三区在线| 嫩草影院新地址| 日韩欧美国产在线观看| 日韩中字成人| 国产亚洲欧美98| 又粗又爽又猛毛片免费看| 国产精华一区二区三区| 国产精品.久久久| 乱人视频在线观看| 高清午夜精品一区二区三区 | 久久这里有精品视频免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 夜夜爽天天搞| 亚洲在久久综合| 变态另类丝袜制服| 精品人妻一区二区三区麻豆| 18禁黄网站禁片免费观看直播| 51国产日韩欧美| 人体艺术视频欧美日本| 性插视频无遮挡在线免费观看| 日本熟妇午夜| 欧美性感艳星| 亚洲精华国产精华液的使用体验 | 午夜精品一区二区三区免费看| 九九久久精品国产亚洲av麻豆| 欧美成人免费av一区二区三区| 深夜a级毛片| 久久精品国产亚洲网站| 久久久久久九九精品二区国产| 一本一本综合久久| 国产熟女欧美一区二区| 91久久精品国产一区二区成人| 搞女人的毛片| 18禁在线播放成人免费| 国产在线男女| 亚洲内射少妇av| 国产成人a区在线观看| 日韩制服骚丝袜av| 亚洲久久久久久中文字幕| www日本黄色视频网| 成年免费大片在线观看| 亚洲av中文av极速乱| 亚洲四区av| 人体艺术视频欧美日本| 中国美女看黄片| 老熟妇乱子伦视频在线观看| 亚洲av第一区精品v没综合| 在线播放无遮挡| 日韩强制内射视频| 成人午夜精彩视频在线观看| 国产成年人精品一区二区| 精品久久久久久久久久久久久| www日本黄色视频网| 女的被弄到高潮叫床怎么办| 亚洲乱码一区二区免费版| 99热6这里只有精品| 亚洲精品亚洲一区二区| 国产黄色视频一区二区在线观看 | 一区二区三区高清视频在线| 在线免费十八禁| 12—13女人毛片做爰片一| 禁无遮挡网站| 日本黄大片高清| 禁无遮挡网站| 日本黄大片高清| 亚洲欧美中文字幕日韩二区| 身体一侧抽搐| 成人av在线播放网站| 久久热精品热| 熟女人妻精品中文字幕| ponron亚洲| 国产精品久久久久久精品电影小说 | av在线天堂中文字幕| 中国美白少妇内射xxxbb| 91久久精品国产一区二区三区| 国产在线男女| 欧美性猛交╳xxx乱大交人| 午夜免费男女啪啪视频观看| 久久久久久大精品| 亚洲精华国产精华液的使用体验 | 国产精品麻豆人妻色哟哟久久 | 精品久久久久久久人妻蜜臀av| 国产黄a三级三级三级人| 在线播放国产精品三级| 午夜a级毛片| 午夜久久久久精精品| 久久中文看片网| 午夜福利成人在线免费观看| avwww免费| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区免费观看| 国产乱人偷精品视频| 日韩制服骚丝袜av| 观看免费一级毛片| 国产精品人妻久久久影院| 国产在线精品亚洲第一网站| 亚洲av成人av| 国产日韩欧美在线精品| 国产成人福利小说| 日韩一区二区视频免费看| 久久久久久久午夜电影| 色5月婷婷丁香| 91在线精品国自产拍蜜月| 国产精品国产高清国产av| 成人无遮挡网站| 成年女人看的毛片在线观看| 精品欧美国产一区二区三| 精品一区二区三区人妻视频| 国产老妇伦熟女老妇高清| 亚洲国产精品久久男人天堂| 国产老妇女一区| 国产淫片久久久久久久久| 亚洲精品日韩在线中文字幕 | 午夜福利高清视频| 国产人妻一区二区三区在| 国产精品久久久久久亚洲av鲁大| 最近2019中文字幕mv第一页| 久久久午夜欧美精品| 夜夜夜夜夜久久久久| 美女大奶头视频| 午夜久久久久精精品| 日韩欧美三级三区| 99热精品在线国产| 亚洲无线在线观看| 日韩在线高清观看一区二区三区| 亚洲av第一区精品v没综合| 婷婷亚洲欧美| 亚洲婷婷狠狠爱综合网| 亚洲国产高清在线一区二区三| 欧美色视频一区免费| av在线观看视频网站免费| 特大巨黑吊av在线直播| 日韩欧美一区二区三区在线观看| 国产大屁股一区二区在线视频|