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

    液艙旋轉(zhuǎn)射流惰化模擬及優(yōu)化

    2021-01-05 03:03:16王恒遠(yuǎn)鄧佳佳盧金樹(shù)薛大文
    造船技術(shù) 2020年6期
    關(guān)鍵詞:惰化出氣口進(jìn)氣口

    王恒遠(yuǎn), 鄧佳佳, 盧金樹(shù), 薛大文

    (浙江海洋大學(xué) 船舶與海運(yùn)學(xué)院, 浙江 舟山 316022)

    0 引 言

    以石油、煤炭為主的燃料引發(fā)的環(huán)境污染問(wèn)題越來(lái)越受到關(guān)注[1]。天然氣憑借經(jīng)濟(jì)環(huán)保、安全性高等特點(diǎn)[2]已成為發(fā)展速度較快的燃料。我國(guó)天然氣資源并不豐富,天然氣消費(fèi)依賴于進(jìn)口,采用液化天然氣(Liquefied Natural Gas, LNG)運(yùn)輸是長(zhǎng)距離運(yùn)輸較經(jīng)濟(jì)的方式。當(dāng)采用新建或檢修后的液艙運(yùn)輸時(shí),需對(duì)液艙進(jìn)行惰化,將氧氣體積分?jǐn)?shù)降低至2%以下[3]。

    常用的惰性氣體置換方法主要包括混合惰化(推移式)、無(wú)混合惰化等4種惰化方式[4]。在4種惰化方式中混合惰化是使用率較高的方法。黃光容等[5]以某機(jī)用燃油箱為對(duì)象,通過(guò)數(shù)值模擬的方式模擬飛機(jī)燃油箱的惰化過(guò)程。KURLE等[6]建立LNG液艙模型,對(duì)液艙惰化過(guò)程及裝載過(guò)程進(jìn)行動(dòng)態(tài)模擬,發(fā)現(xiàn)增大流量可縮短惰化時(shí)間。王志偉[7]通過(guò)數(shù)值模擬的方法得出燃油箱惰化的進(jìn)氣方式、質(zhì)量分?jǐn)?shù)和載油率均影響惰化過(guò)程的進(jìn)行。宋洋[1]應(yīng)用Fluent軟件仿真液貨艙惰化的全過(guò)程,提出進(jìn)氣速度與氣體惰化效率的關(guān)系。但是,這些研究都較少得出惰化所需的具體的進(jìn)氣口大小、方式、位置和進(jìn)氣速度對(duì)惰化的影響,且均采用上進(jìn)下出的進(jìn)氣方式[8]進(jìn)行研究。

    采用數(shù)值仿真方法建立三維模型,分析采用臥式液艙端部、擴(kuò)張角較大的旋轉(zhuǎn)射流[9]等進(jìn)氣方式,探究氣體惰化機(jī)理,并提出優(yōu)化方案以實(shí)現(xiàn)推移式惰化。

    1 模 型

    1.1 物理模型

    實(shí)際工程中LNG液艙尺寸長(zhǎng)為25.5 m,寬為7.5 m,高為7.5 m,總體積約983 m3[10]。在液艙的右端設(shè)有射流噴嘴。按照15∶1的比例尺建立縮尺比模型,如圖1所示。具體尺寸如表1所示。

    圖1 簡(jiǎn)化后物理模型

    表1 簡(jiǎn)化模型參數(shù)

    在初始狀態(tài)下,液艙內(nèi)是溫度為300 K、壓力為0 MPa的空氣。采用純氮?dú)獬淙胍号撨M(jìn)行惰化,氮?dú)鉁囟葹?00 K??諝庵械?dú)饧把鯕獾奈锢韰?shù)如表2所示。

    表2 氮?dú)饧把鯕獾奈锢韰?shù)

    1.2 數(shù)學(xué)模型

    由于實(shí)際惰化過(guò)程復(fù)雜,為簡(jiǎn)化模型,對(duì)其進(jìn)行如下簡(jiǎn)化:

    (1) 將液艙內(nèi)空氣的成分簡(jiǎn)化為21%的氧氣和79%的氮?dú)狻?/p>

    (2) 控制方程主要包括質(zhì)量方程、動(dòng)量方程、組分控制方程和標(biāo)準(zhǔn)k-ε湍流模型。

    1.3 模型設(shè)置

    選擇軟件Gambit 2.4版本進(jìn)行網(wǎng)格劃分。空間內(nèi)流體區(qū)域設(shè)為流體。

    利用Fluent 15.0進(jìn)行模擬,選擇3D計(jì)算器[11-12]。模型方程包含質(zhì)量守恒、動(dòng)量守恒、能量守恒方程,湍流模型選擇標(biāo)準(zhǔn)k-ε湍流模型。設(shè)定y方向重力加速度為-9.81 m2/s。

    進(jìn)氣邊界采用速度進(jìn)口,組分為純氮?dú)猓怀隹诓捎脡毫Τ隹?;其余邊界均為壁面邊界?/p>

    Fluent求解器中壓力速度耦合選擇SIMPLE算法。取壓力松弛因子為0.3,取密度松弛因子為1.0,且動(dòng)量方程、湍流能量耗散、湍流黏度、能量等松弛因子分別取0.7、0.8、1.0、1.0。

    對(duì)所有區(qū)域進(jìn)行初始化,取速度為0 m/s,壓力設(shè)置為0 MPa,初始化溫度取300 K,氧氣平均體積分?jǐn)?shù)為0.231 7。不設(shè)置臨界收斂值,時(shí)間步長(zhǎng)設(shè)為0.01 s,計(jì)算步數(shù)設(shè)為6 000步。

    2 模型驗(yàn)證

    分別就計(jì)算模型、網(wǎng)格數(shù)量敏感度和時(shí)間步長(zhǎng)敏感度等3個(gè)方面開(kāi)展模型驗(yàn)證。

    2.1 計(jì)算模型驗(yàn)證

    建立與試驗(yàn)艙模型等效的物理模型,對(duì)其氣體惰化過(guò)程進(jìn)行數(shù)值模擬,將數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行比較以驗(yàn)證計(jì)算模型的準(zhǔn)確性和科學(xué)性。尺寸設(shè)置為2.20 m×0.92 m×1.22 m的矩形艙,液艙內(nèi)氧氣平均體積分?jǐn)?shù)與試驗(yàn)結(jié)果對(duì)比如圖2所示。由圖2可知:模擬結(jié)果與試驗(yàn)結(jié)果吻合較好。

    圖2 計(jì)算模型驗(yàn)證

    2.2 網(wǎng)格數(shù)量敏感度驗(yàn)證

    粗、中、細(xì)網(wǎng)格方案分別對(duì)應(yīng)網(wǎng)格數(shù)量為25萬(wàn)、53萬(wàn)、120萬(wàn),不同網(wǎng)格數(shù)據(jù)的計(jì)算結(jié)果對(duì)比如表3和圖3所示。由表3和圖3可知:三者的艙內(nèi)平均氧氣體積分?jǐn)?shù)非常接近;惰化完成時(shí)間對(duì)比中,粗、中網(wǎng)格之間的差值百分比為5.3%,而中、細(xì)網(wǎng)格之間的差值百分比僅為3.3%,小于5%。為節(jié)約計(jì)算資源并保證計(jì)算精度,選取中網(wǎng)格為后續(xù)計(jì)算網(wǎng)格。

    表3 網(wǎng)格參數(shù)及差值占比

    圖3 網(wǎng)格數(shù)量敏感性分析結(jié)果

    2.3 時(shí)間步長(zhǎng)敏感度分析

    3種不同時(shí)間步長(zhǎng)0.050 s、0.010 s和0.005 s的計(jì)算結(jié)果對(duì)比如表4和圖4所示。由圖4可以看出:每個(gè)工況的誤差較小。結(jié)合表4中不同網(wǎng)格完成惰化所需時(shí)間對(duì)比,選擇時(shí)間步長(zhǎng)為0.010 s。

    表4 時(shí)間步長(zhǎng)參數(shù)及差值占比

    圖4 時(shí)間步長(zhǎng)敏感性分析結(jié)果對(duì)比

    3 結(jié)果分析

    3.1 無(wú)量綱參數(shù)

    為對(duì)比不同工況惰化所需氮?dú)夂牧?,引入無(wú)量綱時(shí)間τ[8],τ代表輸入的氮?dú)怏w積與儲(chǔ)罐容積的體積比。

    3.2 基準(zhǔn)方案

    以直流射流為基準(zhǔn)方案,氮?dú)鈴膬?chǔ)罐右側(cè)端部進(jìn)入,從儲(chǔ)罐左側(cè)出氣,進(jìn)氣口直徑為0.1 m,進(jìn)氣速度為2 m/s。

    經(jīng)過(guò)數(shù)值模擬計(jì)算,直流射流儲(chǔ)罐內(nèi)氧氣平均體積分?jǐn)?shù)變化曲線如圖5所示,亦得到圖6中τ=0.3(t= 5 s)時(shí)液艙垂直縱剖面氧氣質(zhì)量分?jǐn)?shù)分布圖。由圖6可以看出:直流射流方案惰化的主要問(wèn)題是在噴嘴進(jìn)氣口區(qū)域,由于直流射流結(jié)構(gòu)的擴(kuò)張角較小,噴嘴附近死角區(qū)域惰化難度大,無(wú)法形成推移式惰化,且大量的噴入氮?dú)庵苯訌某鰵饪谂懦觯瑢?dǎo)致惰化時(shí)間長(zhǎng)、氮?dú)夂牧看蟆?/p>

    圖5 3種方案氧氣體積分?jǐn)?shù)變化曲線

    4 優(yōu)化及討論

    根據(jù)基準(zhǔn)方案結(jié)果,提出旋轉(zhuǎn)射流和混合射流兩種進(jìn)氣優(yōu)化方式。旋轉(zhuǎn)射流進(jìn)氣口結(jié)構(gòu)如圖7所示,進(jìn)氣口形狀為0.1 m×0.027 m的矩形,進(jìn)氣方向?yàn)樨Q直切入,從而產(chǎn)生旋轉(zhuǎn)氣流?;旌仙淞鬟M(jìn)氣口結(jié)構(gòu)如圖8所示,設(shè)有半徑為0.05 m的圓形進(jìn)氣口和形狀為0.1 m×0.027 m的矩形進(jìn)氣口,是直流射流與旋轉(zhuǎn)射流的結(jié)合。為了對(duì)比,兩種方案的進(jìn)氣體積流量與基準(zhǔn)方案相同。因此:將旋轉(zhuǎn)射流進(jìn)氣口速度設(shè)為5.82 m/s,使旋轉(zhuǎn)射流速度偏大;對(duì)混合射流的圓形進(jìn)氣口和矩形進(jìn)氣口按照12∶5的比例將速度分別設(shè)為3.62 m/s和1.51 m/s。

    圖7 旋轉(zhuǎn)射流進(jìn)氣口結(jié)構(gòu)

    圖8 混合射流進(jìn)氣口結(jié)構(gòu)

    圖5為3種液艙內(nèi)平均氧氣體積分?jǐn)?shù)和出氣口氮?dú)馄骄w積分?jǐn)?shù)變化曲線,由曲線可以看出:3種方案液艙內(nèi)惰化過(guò)程大體可以分成3個(gè)階段。

    0~4 s(0<τ<0.24)時(shí)間段為初始惰化階段,圖5顯示3種方案氧氣平均體積分?jǐn)?shù)的變化基本一致,儲(chǔ)罐內(nèi)氧氣體積分?jǐn)?shù)下降主要由出氣口附近氧氣被排出所致。結(jié)合t=4 s(τ=0.24)垂直縱剖面氮?dú)赓|(zhì)量分?jǐn)?shù)分布圖(見(jiàn)圖9)可知:直流射流在出氣口處氮?dú)赓|(zhì)量分?jǐn)?shù)最高;由圖5出氣口處氮?dú)怏w積分?jǐn)?shù)曲線看出:0~4 s(0<τ<0.24)時(shí)間段中,除直流射流逐漸增加外,其余方案無(wú)明顯變化。通過(guò)圖5中的出氣口處氮?dú)怏w積分?jǐn)?shù)曲線可知:直流射流方案出氣口氮?dú)饷黠@高于其他方案,使得惰化用氮?dú)饫速M(fèi),導(dǎo)致惰化效果弱化。

    圖9 t=4 s(τ=0.24)時(shí)垂直縱剖面氮?dú)赓|(zhì)量分?jǐn)?shù)分布圖

    4~12 s(0.24<τ<0.72)時(shí)間段為中間惰化階段,圖5顯示該階段3種方案的差異性體現(xiàn)出來(lái),旋轉(zhuǎn)射流和混合射流方案曲線下降速率較快,出氣口處氮?dú)怏w積分?jǐn)?shù)逐漸升高但是仍然小于直流射流。如圖9所示,旋轉(zhuǎn)射流和混合射流的惰化效果明顯優(yōu)于直流射流,主要是因?yàn)橹绷魃淞飨鄬?duì)于其他兩種方案進(jìn)氣擴(kuò)張角較小。但是由于混合射流中圓形進(jìn)氣口的推射作用,在該階段的后期,惰化氮?dú)獾竭_(dá)出氣口的時(shí)間略早于旋轉(zhuǎn)射流(見(jiàn)圖10),氮?dú)怏w積分?jǐn)?shù)增加的速率也快于旋轉(zhuǎn)射流(見(jiàn)圖5)。

    圖10 t=12 s(τ=0.72)時(shí)垂直縱剖面氮?dú)赓|(zhì)量分?jǐn)?shù)分布圖

    結(jié)合圖10可以得出:與中間惰化階段效果相同,由于直流射流形成較小的進(jìn)氣擴(kuò)張角,部分氮?dú)庵苯訌某鰵饪谂懦鲈斐衫速M(fèi),進(jìn)而減慢了惰化過(guò)程的推進(jìn);而旋轉(zhuǎn)射流和混合射流方案的氮?dú)饩黠@離出氣口較遠(yuǎn),基本形成推移式惰化。

    t=12 s(τ=0.72)后直至惰化過(guò)程結(jié)束時(shí)間段為后期惰化階段。選取惰化時(shí)間t=24 s(τ=1.44),此時(shí)旋轉(zhuǎn)射流惰化剛好結(jié)束。由圖11可以看出:在旋轉(zhuǎn)射流方案下儲(chǔ)罐內(nèi)氮?dú)赓|(zhì)量分?jǐn)?shù)最高;同時(shí)圖5中氧氣體積分?jǐn)?shù)下降更快。

    圖11 t=24 s(τ=1.44)時(shí)垂直縱剖面氮?dú)赓|(zhì)量分?jǐn)?shù)分布圖

    由圖5亦可得到:在12~24 s(0.72<τ<1.44)時(shí)間段,混合射流中直射進(jìn)氣口的推射作用被發(fā)揮,部分氮?dú)廨^旋轉(zhuǎn)射流更早地達(dá)到出氣口,但是會(huì)留下部分惰化死角。旋轉(zhuǎn)射流出氣口的氮?dú)怏w積分?jǐn)?shù)在24 s(τ=1.44)后變?yōu)樽罡?,這是由于惰化完成,旋轉(zhuǎn)射流推移式惰化效果較佳,出氣口氮?dú)怏w積分?jǐn)?shù)最高。

    圖12為τ=0.3(t= 5 s)即噴入的氮?dú)饬空純?chǔ)罐容積的30%時(shí)液艙垂直縱剖面流場(chǎng)速度流線圖。由圖12可以看出:進(jìn)入儲(chǔ)罐氮?dú)獾倪\(yùn)動(dòng)狀態(tài),旋轉(zhuǎn)射流比直流射流和混合射流有更大的擴(kuò)張角,使得罐內(nèi)氣體混合較為充分,形成推移式惰化。由τ=0.3(t= 5 s)時(shí)液艙垂直縱剖面氧氣質(zhì)量分?jǐn)?shù)分布圖(見(jiàn)圖13)看出:在進(jìn)入氮?dú)饬肯嗤那闆r下,旋轉(zhuǎn)射流比直流射流和混合射流的惰化效果明顯,噴嘴附近死角區(qū)域的惰化效果明顯改善,惰化進(jìn)程也更快,顯著提高了惰化效率,因此旋轉(zhuǎn)射流效果最好。

    3種方案完成惰化所需的時(shí)間對(duì)比如表5所示。由于進(jìn)氣口流量相同,因此完成時(shí)間長(zhǎng)短與耗氣量成正比。由表5可知:旋轉(zhuǎn)射流方案效果最好,相對(duì)于基準(zhǔn)方案可節(jié)省40.4%氮?dú)饬亢蜁r(shí)間,混合射流方案相對(duì)于基準(zhǔn)方案可節(jié)省19.4%氮?dú)饬亢蜁r(shí)間。

    表5 3種方案惰化完成時(shí)間 s

    綜上所述,旋轉(zhuǎn)射流優(yōu)于混合射流和直流射流,減少了惰化死角的存在,縮短了完成惰化的時(shí)間,節(jié)省了耗氮量。

    5 結(jié) 論

    采用計(jì)算流體力學(xué)技術(shù)對(duì)液艙的氮?dú)舛杌^(guò)程進(jìn)行仿真模擬,分析艙內(nèi)流場(chǎng)、氧氣體積分?jǐn)?shù)和氧氣質(zhì)量分?jǐn)?shù)分布特征,分析液艙不同射流氣體惰化機(jī)理,提出優(yōu)化方案,得到以下結(jié)論:

    (1) 3種射流方案中旋轉(zhuǎn)射流的推移式惰化效果最好。

    (2) 初始階段,由于氮?dú)膺M(jìn)入儲(chǔ)罐的時(shí)間較短,儲(chǔ)罐內(nèi)氧氣體積分?jǐn)?shù)下降較慢,3種方案的變化區(qū)別不大;中期階段,直流射流方案氮?dú)鈴街边M(jìn)入儲(chǔ)罐,進(jìn)氣擴(kuò)張角較小,留下大量惰化死角,減慢了惰化過(guò)程的推進(jìn),而旋轉(zhuǎn)射流和混合射流的氧氣體積分?jǐn)?shù)下降較快,艙內(nèi)氮?dú)怆x出氣口較遠(yuǎn),基本呈現(xiàn)推移式惰化;后期階段,旋轉(zhuǎn)射流方案中儲(chǔ)罐內(nèi)氧氣質(zhì)量分?jǐn)?shù)更低,氧氣體積分?jǐn)?shù)降低更快,氮?dú)赓|(zhì)量分?jǐn)?shù)更高,推移式惰化效果最好,大幅減少氮?dú)庀?,?jié)約時(shí)間。

    (3) 旋轉(zhuǎn)射流優(yōu)于直流射流與混合射流,具有更好的經(jīng)濟(jì)性。與直流射流方案相比,旋轉(zhuǎn)射流可減少40.4%氮?dú)饬亢投杌瘯r(shí)間;與混合射流方案相比,最多可減少26.2%氮?dú)饬亢秃富瘯r(shí)間。與直流射流方案相比,混合射流可節(jié)省19.4%氮?dú)饬亢投杌瘯r(shí)間。

    猜你喜歡
    惰化出氣口進(jìn)氣口
    抽吸氣流量對(duì)催化惰化系統(tǒng)性能影響
    閥門(mén)及氣相沉積設(shè)備
    低溫與特氣(2022年5期)2022-03-23 17:49:51
    延長(zhǎng)筒倉(cāng)存煤時(shí)間防止煤炭自燃的新型惰化技術(shù)及應(yīng)用
    基于仿真的空濾出氣口圓角尺寸對(duì)流阻的影響
    一種SBS橡膠尾氣收集處理裝置
    惰化知覺(jué)研究述評(píng)
    心理研究(2019年1期)2019-12-14 06:04:58
    飛機(jī)輔助進(jìn)氣口小流量工況的性能初探
    科技視界(2017年13期)2017-09-30 16:20:01
    北極熊的運(yùn)氣
    他山之石
    航空模型(2016年3期)2016-05-14 09:00:58
    埋入式進(jìn)氣口優(yōu)化設(shè)計(jì)
    亚洲精品国产色婷婷电影| 在线天堂最新版资源| 少妇被粗大的猛进出69影院 | 黄片无遮挡物在线观看| 国产精品一区二区性色av| 国产成人精品久久久久久| 精品国产一区二区三区久久久樱花| 亚洲经典国产精华液单| 夫妻性生交免费视频一级片| 日本-黄色视频高清免费观看| 插逼视频在线观看| 日韩人妻高清精品专区| 看非洲黑人一级黄片| 伦理电影免费视频| 少妇裸体淫交视频免费看高清| 伦理电影免费视频| 亚洲精品国产成人久久av| 亚洲欧洲日产国产| 国产成人91sexporn| 日韩不卡一区二区三区视频在线| 国产成人91sexporn| av又黄又爽大尺度在线免费看| 国产男人的电影天堂91| 免费久久久久久久精品成人欧美视频 | 色哟哟·www| 色哟哟·www| 久久久久人妻精品一区果冻| 97精品久久久久久久久久精品| 日韩强制内射视频| 少妇熟女欧美另类| 99久久精品一区二区三区| 永久免费av网站大全| 日韩三级伦理在线观看| 简卡轻食公司| 国产午夜精品一二区理论片| 国产精品成人在线| 亚洲色图综合在线观看| 极品教师在线视频| 少妇精品久久久久久久| 另类亚洲欧美激情| 日韩,欧美,国产一区二区三区| 亚洲精华国产精华液的使用体验| 国产成人一区二区在线| 九草在线视频观看| 国产熟女欧美一区二区| 在线看a的网站| 欧美亚洲 丝袜 人妻 在线| 亚洲精品,欧美精品| 国产真实伦视频高清在线观看| 如日韩欧美国产精品一区二区三区 | 欧美 日韩 精品 国产| 精品人妻一区二区三区麻豆| 中文字幕人妻丝袜制服| 日韩av免费高清视频| 女性被躁到高潮视频| 精品一区二区三卡| 午夜日本视频在线| 国产精品女同一区二区软件| 亚洲人成网站在线播| 五月开心婷婷网| 哪个播放器可以免费观看大片| 国产亚洲欧美精品永久| 精品久久久久久久久av| 一级二级三级毛片免费看| 久久国内精品自在自线图片| 嘟嘟电影网在线观看| 老司机亚洲免费影院| 成人18禁高潮啪啪吃奶动态图 | 亚洲av成人精品一二三区| 亚洲国产欧美在线一区| 97在线视频观看| 黄色怎么调成土黄色| 性色avwww在线观看| 老司机亚洲免费影院| 香蕉精品网在线| 久久精品国产a三级三级三级| 精品一品国产午夜福利视频| 亚洲真实伦在线观看| 亚洲国产精品一区三区| 蜜桃久久精品国产亚洲av| 国产白丝娇喘喷水9色精品| 免费播放大片免费观看视频在线观看| 日韩成人av中文字幕在线观看| 亚洲成人一二三区av| 两个人免费观看高清视频 | 国产免费一区二区三区四区乱码| 国产视频内射| 日韩电影二区| 女人久久www免费人成看片| 国产熟女欧美一区二区| 男人舔奶头视频| av福利片在线| 国产精品国产三级国产av玫瑰| 日韩强制内射视频| av在线播放精品| 国内精品宾馆在线| 亚洲av综合色区一区| 中文字幕久久专区| 人妻一区二区av| 99视频精品全部免费 在线| 亚洲激情五月婷婷啪啪| 日韩精品有码人妻一区| 精品少妇久久久久久888优播| 国产精品蜜桃在线观看| 国产精品一区二区在线观看99| 七月丁香在线播放| 麻豆乱淫一区二区| 国产成人精品福利久久| 亚洲国产欧美在线一区| 男人狂女人下面高潮的视频| 在线天堂最新版资源| 国产成人精品婷婷| 亚洲av综合色区一区| 三级经典国产精品| 五月玫瑰六月丁香| 热re99久久精品国产66热6| 国产成人freesex在线| 最新中文字幕久久久久| 国产淫片久久久久久久久| 成年美女黄网站色视频大全免费 | 一区二区三区免费毛片| 久久这里有精品视频免费| 汤姆久久久久久久影院中文字幕| 嫩草影院新地址| 亚洲av在线观看美女高潮| 你懂的网址亚洲精品在线观看| 又粗又硬又长又爽又黄的视频| 狂野欧美白嫩少妇大欣赏| 久久av网站| 3wmmmm亚洲av在线观看| 日韩亚洲欧美综合| 亚洲,欧美,日韩| 久久久久精品久久久久真实原创| 丝袜脚勾引网站| 伦精品一区二区三区| 亚洲av成人精品一二三区| 99久久精品国产国产毛片| 午夜免费观看性视频| 久久99蜜桃精品久久| 在线观看免费高清a一片| 乱码一卡2卡4卡精品| 永久免费av网站大全| 丰满迷人的少妇在线观看| 国产精品嫩草影院av在线观看| 久久久久久久大尺度免费视频| 国产成人精品福利久久| 精品久久久久久久久亚洲| 色视频在线一区二区三区| 国产深夜福利视频在线观看| 亚洲精品自拍成人| 99九九在线精品视频 | 黄色怎么调成土黄色| 久久久久久人妻| 国产在线视频一区二区| 亚洲欧美精品专区久久| 国产熟女午夜一区二区三区 | 久久久久视频综合| 人人澡人人妻人| 男人爽女人下面视频在线观看| freevideosex欧美| 国产一级毛片在线| av国产精品久久久久影院| 国产亚洲精品久久久com| 久久精品国产鲁丝片午夜精品| 国产一区二区三区综合在线观看 | 欧美精品国产亚洲| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 免费黄网站久久成人精品| 国产淫语在线视频| 夫妻性生交免费视频一级片| 午夜av观看不卡| 国国产精品蜜臀av免费| 黄色毛片三级朝国网站 | 汤姆久久久久久久影院中文字幕| 亚洲自偷自拍三级| 我要看黄色一级片免费的| 三上悠亚av全集在线观看 | 黄色毛片三级朝国网站 | 五月玫瑰六月丁香| 三上悠亚av全集在线观看 | 亚洲精品456在线播放app| 国产熟女欧美一区二区| 最近中文字幕高清免费大全6| videos熟女内射| 国产黄色免费在线视频| 91精品伊人久久大香线蕉| 成人综合一区亚洲| www.色视频.com| 51国产日韩欧美| 伦理电影大哥的女人| 一边亲一边摸免费视频| 如何舔出高潮| 午夜日本视频在线| 在现免费观看毛片| av天堂中文字幕网| 黄色怎么调成土黄色| 一级爰片在线观看| 精品国产一区二区久久| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人爽人人夜夜| 午夜久久久在线观看| 免费看日本二区| 亚洲av综合色区一区| 成人亚洲欧美一区二区av| av视频免费观看在线观看| 亚洲美女搞黄在线观看| 亚洲一区二区三区欧美精品| 亚洲国产最新在线播放| 精品人妻熟女毛片av久久网站| 97超视频在线观看视频| 熟妇人妻不卡中文字幕| 一级,二级,三级黄色视频| 美女内射精品一级片tv| 不卡视频在线观看欧美| 91久久精品国产一区二区三区| 亚洲精品乱码久久久v下载方式| 久久久久久久久久久久大奶| 一本久久精品| 国产黄片美女视频| 黑人猛操日本美女一级片| 丝瓜视频免费看黄片| 亚洲av中文av极速乱| 日韩欧美精品免费久久| 免费大片黄手机在线观看| 久久久国产一区二区| 嘟嘟电影网在线观看| 永久免费av网站大全| 久久国产亚洲av麻豆专区| 精品一品国产午夜福利视频| 18禁裸乳无遮挡动漫免费视频| 建设人人有责人人尽责人人享有的| 熟女电影av网| 精品少妇久久久久久888优播| 亚洲三级黄色毛片| av免费在线看不卡| 熟妇人妻不卡中文字幕| 一区在线观看完整版| 五月天丁香电影| 伦理电影大哥的女人| av一本久久久久| 精品人妻一区二区三区麻豆| 免费播放大片免费观看视频在线观看| 亚洲人与动物交配视频| 国产免费一级a男人的天堂| 高清不卡的av网站| 内射极品少妇av片p| 大片电影免费在线观看免费| 丝袜喷水一区| 国产一区二区在线观看av| 99热这里只有是精品在线观看| 国产成人精品久久久久久| 久久 成人 亚洲| 国产在线男女| 一边亲一边摸免费视频| 成人美女网站在线观看视频| 日本vs欧美在线观看视频 | 亚洲伊人久久精品综合| 在线亚洲精品国产二区图片欧美 | 成人影院久久| 夜夜爽夜夜爽视频| 日韩三级伦理在线观看| 91成人精品电影| 精品久久久久久久久亚洲| 在线免费观看不下载黄p国产| 国产有黄有色有爽视频| 好男人视频免费观看在线| 日日摸夜夜添夜夜添av毛片| 国产精品人妻久久久影院| 欧美日韩亚洲高清精品| av播播在线观看一区| 丝袜脚勾引网站| 精品国产露脸久久av麻豆| 免费观看性生交大片5| 国产精品国产av在线观看| 爱豆传媒免费全集在线观看| 97在线人人人人妻| 精品久久久久久电影网| 黑人巨大精品欧美一区二区蜜桃 | 男的添女的下面高潮视频| 精品国产国语对白av| 一个人看视频在线观看www免费| 欧美少妇被猛烈插入视频| 少妇精品久久久久久久| 色吧在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品亚洲一区二区| 2021少妇久久久久久久久久久| 岛国毛片在线播放| 久久久久久久精品精品| 在线观看免费高清a一片| 在线观看av片永久免费下载| 久久国产乱子免费精品| 亚洲欧美日韩卡通动漫| 内地一区二区视频在线| 亚洲成色77777| 久久久久久伊人网av| 中文字幕精品免费在线观看视频 | 久久久久视频综合| 国产 精品1| 亚洲,欧美,日韩| 人人妻人人澡人人看| 插阴视频在线观看视频| 午夜视频国产福利| 久久99蜜桃精品久久| 国产成人a∨麻豆精品| 欧美97在线视频| 国产一区二区三区av在线| 精品久久久精品久久久| 人人妻人人看人人澡| 亚洲欧美成人综合另类久久久| 老司机亚洲免费影院| 男人添女人高潮全过程视频| 夜夜爽夜夜爽视频| 汤姆久久久久久久影院中文字幕| 高清不卡的av网站| 黄色欧美视频在线观看| av免费观看日本| 亚洲国产最新在线播放| 国产精品国产av在线观看| 精品久久久久久久久亚洲| av不卡在线播放| 国产片特级美女逼逼视频| 黄色一级大片看看| 99热网站在线观看| 精品久久久久久久久亚洲| 亚洲欧美清纯卡通| 久久av网站| 国产一区亚洲一区在线观看| 自线自在国产av| 欧美精品一区二区免费开放| 熟女人妻精品中文字幕| 久久久久久久久久人人人人人人| 9色porny在线观看| av专区在线播放| 十分钟在线观看高清视频www | 日日摸夜夜添夜夜添av毛片| 国产国拍精品亚洲av在线观看| av线在线观看网站| 免费人妻精品一区二区三区视频| 色婷婷久久久亚洲欧美| av福利片在线观看| 日韩欧美精品免费久久| 肉色欧美久久久久久久蜜桃| 中国美白少妇内射xxxbb| 视频区图区小说| 老司机影院成人| 亚洲无线观看免费| 精品久久久久久电影网| 欧美xxxx性猛交bbbb| www.av在线官网国产| 少妇的逼好多水| 国产精品蜜桃在线观看| 国产免费一级a男人的天堂| 最近中文字幕高清免费大全6| 国产色婷婷99| 国产一区二区在线观看日韩| 成年人午夜在线观看视频| 成年人免费黄色播放视频 | av在线播放精品| 亚洲丝袜综合中文字幕| 亚洲国产av新网站| a 毛片基地| 久久99一区二区三区| 国内精品宾馆在线| 欧美精品人与动牲交sv欧美| 久久人人爽人人片av| 这个男人来自地球电影免费观看 | 99国产精品免费福利视频| 偷拍熟女少妇极品色| 一级毛片aaaaaa免费看小| 日韩精品有码人妻一区| 一区二区三区免费毛片| 啦啦啦视频在线资源免费观看| 成人美女网站在线观看视频| 欧美日韩精品成人综合77777| av不卡在线播放| 尾随美女入室| 亚洲av.av天堂| 香蕉精品网在线| a级毛片在线看网站| 午夜91福利影院| 久久av网站| 国产亚洲欧美精品永久| 特大巨黑吊av在线直播| 边亲边吃奶的免费视频| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 久久久久精品久久久久真实原创| 免费黄色在线免费观看| 黑丝袜美女国产一区| a级毛片免费高清观看在线播放| 韩国高清视频一区二区三区| 永久网站在线| 久久韩国三级中文字幕| 赤兔流量卡办理| 26uuu在线亚洲综合色| 亚洲va在线va天堂va国产| 国产伦在线观看视频一区| 精品午夜福利在线看| 欧美3d第一页| 在线看a的网站| 六月丁香七月| 国产淫语在线视频| 99久久精品热视频| 久久99热6这里只有精品| 美女脱内裤让男人舔精品视频| 午夜免费鲁丝| 亚洲av欧美aⅴ国产| 久久精品久久久久久噜噜老黄| 国产中年淑女户外野战色| 日本免费在线观看一区| 18禁在线播放成人免费| 日韩视频在线欧美| 日韩不卡一区二区三区视频在线| 另类精品久久| 你懂的网址亚洲精品在线观看| 欧美性感艳星| 两个人免费观看高清视频 | 午夜老司机福利剧场| 人妻制服诱惑在线中文字幕| 精品一区在线观看国产| 久久久久久久国产电影| 久久99蜜桃精品久久| 丝袜脚勾引网站| 欧美日韩一区二区视频在线观看视频在线| 国产高清国产精品国产三级| 又爽又黄a免费视频| 国产精品偷伦视频观看了| 男人和女人高潮做爰伦理| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看| 在线亚洲精品国产二区图片欧美 | 亚洲欧美精品专区久久| 婷婷色综合www| 伦精品一区二区三区| 色吧在线观看| 妹子高潮喷水视频| 在线观看av片永久免费下载| 男人和女人高潮做爰伦理| 九草在线视频观看| 日韩av免费高清视频| 91成人精品电影| 日韩不卡一区二区三区视频在线| 99热6这里只有精品| 日本-黄色视频高清免费观看| 国产成人a∨麻豆精品| 亚洲欧美日韩东京热| 成人18禁高潮啪啪吃奶动态图 | 亚洲国产成人一精品久久久| 国内少妇人妻偷人精品xxx网站| 国语对白做爰xxxⅹ性视频网站| 伦理电影大哥的女人| 午夜免费鲁丝| 国产精品人妻久久久久久| 制服丝袜香蕉在线| 国国产精品蜜臀av免费| 九九爱精品视频在线观看| 男人舔奶头视频| 一区二区三区乱码不卡18| 99re6热这里在线精品视频| 国产成人免费观看mmmm| 超碰97精品在线观看| 纵有疾风起免费观看全集完整版| 美女cb高潮喷水在线观看| 国产有黄有色有爽视频| 国产成人a∨麻豆精品| 内射极品少妇av片p| 久久久久人妻精品一区果冻| 人妻 亚洲 视频| 一本色道久久久久久精品综合| 中文字幕制服av| 国产永久视频网站| 久久午夜福利片| 精品久久久噜噜| 国产黄频视频在线观看| 精品一区二区三区视频在线| 午夜视频国产福利| 26uuu在线亚洲综合色| 97精品久久久久久久久久精品| 久久人妻熟女aⅴ| 王馨瑶露胸无遮挡在线观看| 黄色一级大片看看| 少妇裸体淫交视频免费看高清| 18+在线观看网站| 免费av中文字幕在线| 精品一品国产午夜福利视频| 欧美亚洲 丝袜 人妻 在线| 国产精品免费大片| 少妇精品久久久久久久| 女人精品久久久久毛片| 国产精品蜜桃在线观看| 男女边摸边吃奶| 国产亚洲精品久久久com| 欧美激情极品国产一区二区三区 | 麻豆成人av视频| av在线播放精品| 精品酒店卫生间| 日日摸夜夜添夜夜爱| 熟女电影av网| 韩国高清视频一区二区三区| 99热网站在线观看| 日韩在线高清观看一区二区三区| 精品国产国语对白av| 久久久亚洲精品成人影院| 在线亚洲精品国产二区图片欧美 | 简卡轻食公司| 26uuu在线亚洲综合色| 久久国产精品男人的天堂亚洲 | 国产精品欧美亚洲77777| 免费大片黄手机在线观看| 国产精品熟女久久久久浪| 亚洲美女黄色视频免费看| 精品亚洲成a人片在线观看| 日韩强制内射视频| 欧美日韩综合久久久久久| 在线免费观看不下载黄p国产| 亚洲av成人精品一二三区| 少妇精品久久久久久久| 777米奇影视久久| 亚洲图色成人| 欧美老熟妇乱子伦牲交| 欧美另类一区| 亚洲人成网站在线观看播放| 男女国产视频网站| 中文字幕久久专区| 九色成人免费人妻av| 秋霞伦理黄片| 亚洲欧洲精品一区二区精品久久久 | 日本免费在线观看一区| 久久国内精品自在自线图片| 一级毛片黄色毛片免费观看视频| 热re99久久国产66热| 久久女婷五月综合色啪小说| 国产免费一区二区三区四区乱码| 大陆偷拍与自拍| 欧美 亚洲 国产 日韩一| 国产成人精品久久久久久| 国产精品久久久久久精品电影小说| 久久免费观看电影| 麻豆乱淫一区二区| 看非洲黑人一级黄片| 欧美精品高潮呻吟av久久| 色视频在线一区二区三区| 日本黄大片高清| 国产精品久久久久久精品古装| 大码成人一级视频| 一个人看视频在线观看www免费| 色婷婷久久久亚洲欧美| 久久精品国产自在天天线| 男人添女人高潮全过程视频| 日日爽夜夜爽网站| 精品人妻偷拍中文字幕| 熟妇人妻不卡中文字幕| 伦理电影免费视频| 国产精品久久久久久久电影| 妹子高潮喷水视频| 又粗又硬又长又爽又黄的视频| 免费大片黄手机在线观看| 熟女人妻精品中文字幕| 亚洲欧美成人综合另类久久久| 国产精品麻豆人妻色哟哟久久| 夫妻午夜视频| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品乱码久久久v下载方式| 久久毛片免费看一区二区三区| 国产日韩一区二区三区精品不卡 | 极品少妇高潮喷水抽搐| 免费不卡的大黄色大毛片视频在线观看| 青春草视频在线免费观看| av一本久久久久| 十八禁网站网址无遮挡 | 久久女婷五月综合色啪小说| 国产精品熟女久久久久浪| 国产一区二区在线观看日韩| 国产av一区二区精品久久| 王馨瑶露胸无遮挡在线观看| 国产精品偷伦视频观看了| 亚洲第一av免费看| 少妇猛男粗大的猛烈进出视频| 国产黄色免费在线视频| 国产高清三级在线| 美女xxoo啪啪120秒动态图| av福利片在线观看| 超碰97精品在线观看| 一级a做视频免费观看| 成人亚洲欧美一区二区av| 婷婷色av中文字幕| 国产亚洲欧美精品永久| .国产精品久久| 天堂俺去俺来也www色官网| 菩萨蛮人人尽说江南好唐韦庄| 国产成人一区二区在线| 下体分泌物呈黄色| 七月丁香在线播放| 亚洲欧美日韩东京热| 国产精品嫩草影院av在线观看| 国内揄拍国产精品人妻在线| 日韩成人伦理影院| 中文乱码字字幕精品一区二区三区| 国产精品秋霞免费鲁丝片| 在线播放无遮挡| 日本欧美国产在线视频| 精品国产一区二区久久| 高清视频免费观看一区二区| 国产精品久久久久久精品电影小说| 极品人妻少妇av视频| 一个人看视频在线观看www免费| 亚洲欧美一区二区三区国产| 中国三级夫妇交换| 亚洲中文av在线| 国产乱来视频区| 高清av免费在线| 精品国产一区二区三区久久久樱花| 国产精品久久久久久久久免|