• <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車燃油消耗量大
    午夜精品国产一区二区电影 | 波野结衣二区三区在线| 一级黄片播放器| 日韩欧美三级三区| 草草在线视频免费看| 成人三级黄色视频| 婷婷色综合大香蕉| 18禁黄网站禁片免费观看直播| 男人舔奶头视频| 2021天堂中文幕一二区在线观| 少妇的逼水好多| 国产av麻豆久久久久久久| 欧美变态另类bdsm刘玥| 免费人成在线观看视频色| 精品日产1卡2卡| 久99久视频精品免费| 99在线视频只有这里精品首页| 一边摸一边抽搐一进一小说| 六月丁香七月| 一边摸一边抽搐一进一小说| 1024手机看黄色片| 精品人妻熟女av久视频| 免费观看的影片在线观看| 亚洲精品久久国产高清桃花| 中国美女看黄片| 亚洲精品影视一区二区三区av| 日韩欧美精品免费久久| 久久九九热精品免费| 在线免费十八禁| 桃色一区二区三区在线观看| 91精品国产九色| 给我免费播放毛片高清在线观看| 能在线免费看毛片的网站| 亚洲精品久久国产高清桃花| 好男人在线观看高清免费视频| 亚洲av一区综合| 亚洲精品国产成人久久av| 中文字幕熟女人妻在线| 成人漫画全彩无遮挡| av在线播放精品| 一级毛片我不卡| 男人和女人高潮做爰伦理| 禁无遮挡网站| av女优亚洲男人天堂| 久久99热6这里只有精品| a级毛色黄片| 桃色一区二区三区在线观看| 国产精品蜜桃在线观看 | 日韩,欧美,国产一区二区三区 | 熟妇人妻久久中文字幕3abv| 中文资源天堂在线| 国产探花极品一区二区| 国语自产精品视频在线第100页| 欧美区成人在线视频| 国产一级毛片七仙女欲春2| 精品熟女少妇av免费看| 日韩一本色道免费dvd| 亚洲精品自拍成人| 黄片无遮挡物在线观看| 久久精品国产自在天天线| 99热只有精品国产| 长腿黑丝高跟| 成年女人永久免费观看视频| 色综合站精品国产| 偷拍熟女少妇极品色| 在线观看美女被高潮喷水网站| 综合色丁香网| 亚洲国产精品成人综合色| 18禁裸乳无遮挡免费网站照片| 99久国产av精品国产电影| 国产人妻一区二区三区在| 成人性生交大片免费视频hd| 中文字幕精品亚洲无线码一区| 免费看av在线观看网站| 欧美极品一区二区三区四区| 尤物成人国产欧美一区二区三区| 男女那种视频在线观看| 成人午夜精彩视频在线观看| 亚洲精华国产精华液的使用体验 | 久久久午夜欧美精品| 国产亚洲精品久久久久久毛片| 人妻制服诱惑在线中文字幕| 欧美不卡视频在线免费观看| 久久精品91蜜桃| 国产极品天堂在线| 久久亚洲精品不卡| 国产一区二区在线av高清观看| 青青草视频在线视频观看| 国产蜜桃级精品一区二区三区| 黄色日韩在线| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久久免费av| 波野结衣二区三区在线| 国产精品不卡视频一区二区| 色噜噜av男人的天堂激情| 日本黄色视频三级网站网址| 国产成人freesex在线| 欧美变态另类bdsm刘玥| 乱系列少妇在线播放| 在线免费十八禁| 亚洲成a人片在线一区二区| 直男gayav资源| 小蜜桃在线观看免费完整版高清| 精品少妇黑人巨大在线播放 | 不卡视频在线观看欧美| 夜夜看夜夜爽夜夜摸| 免费看a级黄色片| 成人三级黄色视频| 99热网站在线观看| 丰满乱子伦码专区| 一夜夜www| 亚洲欧美日韩高清专用| 国产一级毛片在线| 精品久久久久久久久久久久久| 国内精品久久久久精免费| 久久久久免费精品人妻一区二区| 国产毛片a区久久久久| 国产一级毛片七仙女欲春2| 亚洲av中文av极速乱| 丰满的人妻完整版| 色综合站精品国产| 久久午夜福利片| 精品无人区乱码1区二区| 免费看av在线观看网站| 高清在线视频一区二区三区 | 午夜精品一区二区三区免费看| 欧美高清性xxxxhd video| 亚洲av成人精品一区久久| 人妻制服诱惑在线中文字幕| 亚洲欧美精品自产自拍| 精华霜和精华液先用哪个| 一边亲一边摸免费视频| 免费人成视频x8x8入口观看| 美女 人体艺术 gogo| 亚洲国产精品国产精品| 少妇的逼水好多| 日本黄色视频三级网站网址| 亚洲国产精品sss在线观看| 国产精品av视频在线免费观看| 国产成年人精品一区二区| 超碰av人人做人人爽久久| 在线观看午夜福利视频| 一级黄片播放器| 亚洲欧美日韩高清在线视频| 人妻制服诱惑在线中文字幕| 国产一级毛片七仙女欲春2| 国产不卡一卡二| or卡值多少钱| 一边摸一边抽搐一进一小说| 婷婷色av中文字幕| 成人午夜精彩视频在线观看| 欧美区成人在线视频| 久久99精品国语久久久| 丰满的人妻完整版| 少妇裸体淫交视频免费看高清| 国产亚洲av嫩草精品影院| 蜜桃亚洲精品一区二区三区| 一本精品99久久精品77| 日本一本二区三区精品| 国产精品99久久久久久久久| 午夜福利高清视频| 亚洲av免费在线观看| 91久久精品国产一区二区成人| 秋霞在线观看毛片| 黄色欧美视频在线观看| av天堂中文字幕网| 日产精品乱码卡一卡2卡三| 亚洲国产欧洲综合997久久,| 狠狠狠狠99中文字幕| 熟女人妻精品中文字幕| 国产av一区在线观看免费| 精品少妇黑人巨大在线播放 | 国产亚洲5aaaaa淫片| 我的女老师完整版在线观看| 一区二区三区四区激情视频 | 男人舔奶头视频| 久久久久久久久久黄片| 亚洲高清免费不卡视频| 色视频www国产| 午夜视频国产福利| 国产精品一区二区性色av| 丰满的人妻完整版| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产亚洲网站| 国产成人精品一,二区 | 能在线免费看毛片的网站| 午夜激情福利司机影院| 成人综合一区亚洲| 国产精品国产高清国产av| 男人舔奶头视频| 久久这里只有精品中国| 色综合站精品国产| 成人性生交大片免费视频hd| 老女人水多毛片| 波野结衣二区三区在线| 国产av在哪里看| 国产成人福利小说| 午夜福利视频1000在线观看| 女的被弄到高潮叫床怎么办| 欧美变态另类bdsm刘玥| 国产女主播在线喷水免费视频网站 | 欧美三级亚洲精品| www日本黄色视频网| 精品久久久久久久久av| 国产精品久久视频播放| 亚洲欧美日韩高清专用| 十八禁国产超污无遮挡网站| 国产精品99久久久久久久久| 国内少妇人妻偷人精品xxx网站| 日日撸夜夜添| 欧美zozozo另类| 久久人人爽人人片av| 亚洲av成人精品一区久久| 日本黄色视频三级网站网址| 国内揄拍国产精品人妻在线| 欧美性猛交╳xxx乱大交人| 久久久久网色| 亚洲欧美精品综合久久99| 亚洲三级黄色毛片| 岛国毛片在线播放| 亚洲精品日韩在线中文字幕 | 麻豆av噜噜一区二区三区| 国产av麻豆久久久久久久| 中文亚洲av片在线观看爽| 男女下面进入的视频免费午夜| av在线天堂中文字幕| 国产真实乱freesex| 亚洲国产精品国产精品| 亚洲欧美成人精品一区二区| 在线天堂最新版资源| 91aial.com中文字幕在线观看| av又黄又爽大尺度在线免费看 | 成人综合一区亚洲| 天堂影院成人在线观看| 在线播放无遮挡| 一进一出抽搐gif免费好疼| 日韩三级伦理在线观看| 国产免费男女视频| 国产伦一二天堂av在线观看| 在线观看66精品国产| av国产免费在线观看| 国产一区二区亚洲精品在线观看| 日本一二三区视频观看| 国内久久婷婷六月综合欲色啪| 国产成人freesex在线| 午夜福利在线在线| АⅤ资源中文在线天堂| 99久久精品国产国产毛片| 又爽又黄a免费视频| 国产毛片a区久久久久| 日韩中字成人| 最近中文字幕高清免费大全6| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久中文| 国产精品综合久久久久久久免费| 亚洲精品粉嫩美女一区| 亚洲va在线va天堂va国产| 久久久久久久亚洲中文字幕| 国产老妇女一区| 精品国产三级普通话版| 六月丁香七月| 亚洲第一电影网av| 岛国毛片在线播放| www.色视频.com| 国产精品久久视频播放| 亚洲国产精品国产精品| 亚洲国产高清在线一区二区三| 国产精品,欧美在线| 成人无遮挡网站| 中出人妻视频一区二区| 少妇熟女aⅴ在线视频| 少妇被粗大猛烈的视频| 一区二区三区四区激情视频 | 亚洲激情五月婷婷啪啪| 九九热线精品视视频播放| 久久精品国产亚洲网站| 国产探花极品一区二区| eeuss影院久久| 97热精品久久久久久| 日本一本二区三区精品| 身体一侧抽搐| 99在线人妻在线中文字幕| 最近2019中文字幕mv第一页| 一级毛片aaaaaa免费看小| 身体一侧抽搐| 精品久久国产蜜桃| 亚洲国产精品成人久久小说 | 国产精品人妻久久久影院| 春色校园在线视频观看| 丝袜喷水一区| 免费观看在线日韩| 亚洲熟妇中文字幕五十中出| 99在线视频只有这里精品首页| 亚洲丝袜综合中文字幕| 日韩一本色道免费dvd| 久久综合国产亚洲精品| 欧美成人精品欧美一级黄| 中国美白少妇内射xxxbb| 18禁黄网站禁片免费观看直播| 日本av手机在线免费观看| 啦啦啦啦在线视频资源| 色5月婷婷丁香| avwww免费| 美女cb高潮喷水在线观看| 天堂av国产一区二区熟女人妻| 国产av在哪里看| 欧美性猛交╳xxx乱大交人| 免费在线观看成人毛片| 欧美潮喷喷水| 日韩一本色道免费dvd| 精品免费久久久久久久清纯| 老司机福利观看| 日本成人三级电影网站| 亚洲国产欧美在线一区| 少妇人妻精品综合一区二区 | 不卡一级毛片| 97在线视频观看| 日本黄色片子视频| 国产毛片a区久久久久| 国产高清激情床上av| 国内精品宾馆在线| 欧美日本视频| 美女大奶头视频| 亚洲五月天丁香| 精品人妻一区二区三区麻豆| 你懂的网址亚洲精品在线观看 | 97超视频在线观看视频| 亚洲精品久久国产高清桃花| 亚洲在线观看片| 国产亚洲5aaaaa淫片| 一区二区三区四区激情视频 | 91麻豆精品激情在线观看国产| 99riav亚洲国产免费| 国产成人a区在线观看| 欧美一区二区精品小视频在线| 少妇裸体淫交视频免费看高清| 欧美高清性xxxxhd video| 国产极品精品免费视频能看的| 少妇的逼好多水| 国产中年淑女户外野战色| 黄色一级大片看看| 国产精品免费一区二区三区在线| 男人和女人高潮做爰伦理| 免费一级毛片在线播放高清视频| 日产精品乱码卡一卡2卡三| 女人十人毛片免费观看3o分钟| 亚洲国产精品久久男人天堂| 亚洲一区二区三区色噜噜| 黄色视频,在线免费观看| 国产91av在线免费观看| 不卡一级毛片| 黄色视频,在线免费观看| 日本三级黄在线观看| 日本免费一区二区三区高清不卡| 少妇裸体淫交视频免费看高清| 十八禁国产超污无遮挡网站| 久久热精品热| 亚洲av第一区精品v没综合| 国产亚洲av片在线观看秒播厂 | 岛国在线免费视频观看| 级片在线观看| ponron亚洲| 欧美在线一区亚洲| 亚洲国产精品国产精品| 中文亚洲av片在线观看爽| 最近中文字幕高清免费大全6| 丝袜喷水一区| 久久人人精品亚洲av| 午夜亚洲福利在线播放| 欧美一级a爱片免费观看看| 97超碰精品成人国产| 久久人人精品亚洲av| 亚洲va在线va天堂va国产| 中文资源天堂在线| 久久久久久久久大av| 非洲黑人性xxxx精品又粗又长| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看 | 国产成人a区在线观看| 村上凉子中文字幕在线| 日韩欧美精品免费久久| 麻豆一二三区av精品| 精品国产三级普通话版| 人妻少妇偷人精品九色| 色综合色国产| 日韩欧美一区二区三区在线观看| 51国产日韩欧美| 在线a可以看的网站| 深夜精品福利| 成人二区视频| 亚洲中文字幕一区二区三区有码在线看| 美女xxoo啪啪120秒动态图| 久久久久久久久久久丰满| 亚洲精品自拍成人| 久久午夜福利片| kizo精华| 亚洲美女视频黄频| 一级黄片播放器| 国产精品一区二区在线观看99 | 特大巨黑吊av在线直播| 亚洲中文字幕一区二区三区有码在线看| 少妇裸体淫交视频免费看高清| 中文在线观看免费www的网站| 2021天堂中文幕一二区在线观| 日本一二三区视频观看| 欧美+日韩+精品| 又粗又硬又长又爽又黄的视频 | www.av在线官网国产| 亚洲熟妇中文字幕五十中出| 亚洲精品456在线播放app| 国内精品美女久久久久久| 日日干狠狠操夜夜爽| 欧美最黄视频在线播放免费| 18禁黄网站禁片免费观看直播| 人妻少妇偷人精品九色| 亚洲经典国产精华液单| 我要看日韩黄色一级片| 一本精品99久久精品77| 久久欧美精品欧美久久欧美| 午夜福利在线观看吧| 国产欧美日韩精品一区二区| 久久精品国产自在天天线| 国产片特级美女逼逼视频| 欧美不卡视频在线免费观看| АⅤ资源中文在线天堂| 亚洲av电影不卡..在线观看| 在线观看av片永久免费下载| 99久久无色码亚洲精品果冻| 成人午夜精彩视频在线观看| eeuss影院久久| 男人和女人高潮做爰伦理| 久久久久国产网址| 亚洲精品粉嫩美女一区| 亚洲不卡免费看| 神马国产精品三级电影在线观看| 99久久无色码亚洲精品果冻| 午夜精品一区二区三区免费看| 亚洲精品成人久久久久久| 国产伦理片在线播放av一区 | 国产精品乱码一区二三区的特点| 欧美区成人在线视频| 国产精品精品国产色婷婷| 中国美白少妇内射xxxbb| 自拍偷自拍亚洲精品老妇| 能在线免费观看的黄片| 国产三级中文精品| 国产成人aa在线观看| 国产高潮美女av| 成人性生交大片免费视频hd| 午夜福利在线在线| 国产精品蜜桃在线观看 | 国产精品一及| 男人和女人高潮做爰伦理| 欧美潮喷喷水| 国产午夜精品久久久久久一区二区三区| 女人十人毛片免费观看3o分钟| 久久久a久久爽久久v久久| 男女啪啪激烈高潮av片| 免费大片18禁| 18+在线观看网站| 校园人妻丝袜中文字幕| 欧美最新免费一区二区三区| 亚洲国产精品成人综合色| 在线a可以看的网站| 国内精品宾馆在线| 2022亚洲国产成人精品| 国产精品蜜桃在线观看 | 国产一级毛片在线| eeuss影院久久| 亚洲高清免费不卡视频| 久久久国产成人精品二区| 国内精品宾馆在线| 国产精品永久免费网站| 白带黄色成豆腐渣| 男女做爰动态图高潮gif福利片| 美女xxoo啪啪120秒动态图| 一级黄片播放器| 亚洲三级黄色毛片| 三级国产精品欧美在线观看| 长腿黑丝高跟| 麻豆成人午夜福利视频| 精品日产1卡2卡| 免费无遮挡裸体视频| 少妇裸体淫交视频免费看高清| 色综合色国产| 精品久久久久久久久亚洲| 少妇的逼好多水| 亚洲欧美精品综合久久99| 中文字幕精品亚洲无线码一区| av黄色大香蕉| 亚洲av成人av| 日日摸夜夜添夜夜添av毛片| 欧美激情在线99| 人体艺术视频欧美日本| 国产精品美女特级片免费视频播放器| 国产精品无大码| 别揉我奶头 嗯啊视频| 深夜精品福利| 久久亚洲国产成人精品v| 成人午夜高清在线视频| 中文亚洲av片在线观看爽| 成人性生交大片免费视频hd| 亚洲精品日韩av片在线观看| 一本精品99久久精品77| 国产一区二区在线av高清观看| 麻豆成人av视频| 一边摸一边抽搐一进一小说| 91久久精品电影网| 国产精品伦人一区二区| 99久久中文字幕三级久久日本| 99久久精品一区二区三区| 超碰av人人做人人爽久久| 亚洲av二区三区四区| 精品一区二区三区视频在线| 亚洲成人精品中文字幕电影| 国产大屁股一区二区在线视频| 国产黄a三级三级三级人| 日韩成人伦理影院| 午夜老司机福利剧场| 亚洲精品乱码久久久久久按摩| 亚洲av不卡在线观看| 色吧在线观看| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 你懂的网址亚洲精品在线观看 | 欧美极品一区二区三区四区| 黑人高潮一二区| 国产亚洲精品av在线| 免费搜索国产男女视频| 国产女主播在线喷水免费视频网站 | 99久久九九国产精品国产免费| 春色校园在线视频观看| 99国产极品粉嫩在线观看| 精品久久久久久成人av| 哪里可以看免费的av片| 一个人看视频在线观看www免费| 在线观看午夜福利视频| 99国产极品粉嫩在线观看| 久久久久久久久久久丰满| 亚洲内射少妇av| 免费看光身美女| 禁无遮挡网站| 99国产极品粉嫩在线观看| 国产av一区在线观看免费| 久久久精品94久久精品| 国产精品国产高清国产av| 亚洲欧美日韩高清专用| 亚洲精品日韩av片在线观看| 午夜精品一区二区三区免费看| 免费大片18禁| 国产成人精品久久久久久| 色视频www国产| 国产私拍福利视频在线观看| 中国美白少妇内射xxxbb| 18禁裸乳无遮挡免费网站照片| 亚洲精品亚洲一区二区| 高清毛片免费看| a级毛片免费高清观看在线播放| 国产毛片a区久久久久| 国产片特级美女逼逼视频| 97人妻精品一区二区三区麻豆| 久久久久久久久久久丰满| 狂野欧美激情性xxxx在线观看| 国产黄片美女视频| 久99久视频精品免费| 婷婷色综合大香蕉| 亚洲三级黄色毛片| 色播亚洲综合网| 日韩欧美 国产精品| 欧美成人精品欧美一级黄| 精品午夜福利在线看| 特级一级黄色大片| 毛片女人毛片| www日本黄色视频网| 青春草视频在线免费观看| 91在线精品国自产拍蜜月| 寂寞人妻少妇视频99o| 久久精品国产99精品国产亚洲性色| 爱豆传媒免费全集在线观看| 特级一级黄色大片| 亚洲国产色片| 一区福利在线观看| 看片在线看免费视频| 97热精品久久久久久| 久久99蜜桃精品久久| 久久久a久久爽久久v久久| 搡老妇女老女人老熟妇| 国内精品美女久久久久久| 最近最新中文字幕大全电影3| av在线天堂中文字幕| 亚洲电影在线观看av| 国产一级毛片在线| 麻豆国产97在线/欧美| 亚洲经典国产精华液单| 国产亚洲av嫩草精品影院| 18禁黄网站禁片免费观看直播| 中文字幕久久专区| 欧美激情在线99| 此物有八面人人有两片| 中文字幕人妻熟人妻熟丝袜美| 嫩草影院精品99| 国产精品久久久久久久电影| 久久精品国产99精品国产亚洲性色| 天堂√8在线中文| 免费无遮挡裸体视频| 麻豆乱淫一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品国产三级普通话版| 国产在线男女| 日日摸夜夜添夜夜爱| 大型黄色视频在线免费观看| 国内久久婷婷六月综合欲色啪| 免费看美女性在线毛片视频| 身体一侧抽搐| 亚洲美女搞黄在线观看|