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

    壁面催化條件對(duì)氣動(dòng)熱環(huán)境的影響研究

    2020-03-27 08:23:52李旭東史增民
    關(guān)鍵詞:駐點(diǎn)馬赫數(shù)激波

    李旭東,張 賦,史增民,王 欣

    壁面催化條件對(duì)氣動(dòng)熱環(huán)境的影響研究

    李旭東1,2,張 賦1,2,史增民1,2,王 欣1,2

    (1. 北京航天長(zhǎng)征飛行器研究所,北京,100076;2. 高超聲速飛行器防隔熱與熱控技術(shù)中心,北京,100076)

    為考察再入過程中壁面催化效應(yīng)對(duì)氣動(dòng)熱環(huán)境的影響,針對(duì)高超聲速熱化學(xué)非平衡流動(dòng)條件下不同壁面條件(全催化壁面/非催化壁面)對(duì)氣動(dòng)加熱的影響規(guī)律展開了研究。基于熱化學(xué)非平衡流動(dòng)特征構(gòu)造了數(shù)值算法,并通過激波-膨脹風(fēng)洞中開展的高焓高超聲速地面試驗(yàn)結(jié)果進(jìn)行了驗(yàn)證,然后考察了不同飛行高度下和典型幾何外形下不同壁面條件對(duì)氣動(dòng)熱環(huán)境的影響并分析了其影響機(jī)制。研究結(jié)果表明在50 km的飛行高度上催化效應(yīng)明顯,此時(shí)全催化/非催化壁面條件下熱流存在明顯差異,這種差異隨著高度降低而減小。

    高超聲速流動(dòng);真實(shí)氣體效應(yīng);催化壁面;數(shù)值方法

    0 引 言

    飛行器再入地球或者進(jìn)入火星大氣層時(shí)面臨非常復(fù)雜的流動(dòng)環(huán)境,由于強(qiáng)激波壓縮,激波層內(nèi)氣體溫度急劇上升,高溫引起氣體分子振動(dòng)能激發(fā),產(chǎn)生離解甚至電離以及電子激發(fā)和輻射等一系列復(fù)雜的物理、化學(xué)現(xiàn)象[1],而當(dāng)氣體高溫離解后流動(dòng)到具有催化能力的壁面時(shí),已經(jīng)離解的原子會(huì)重新復(fù)合成為分子而在壁面上釋放出大量的化學(xué)熱,導(dǎo)致壁面加熱量上升。當(dāng)前飛行器飛行時(shí)間越來越長(zhǎng),如果不能準(zhǔn)確評(píng)估壁面催化效應(yīng)的影響將造成氣動(dòng)熱環(huán)境預(yù)測(cè)精度偏差較大。長(zhǎng)時(shí)間飛行下這種氣動(dòng)熱預(yù)測(cè)偏差將導(dǎo)致防隔熱系統(tǒng)設(shè)計(jì)方案偏保守,嚴(yán)重影響和制約了飛行器所能達(dá)到的水平。

    本文以典型球錐外形——ELECTRE球錐為基礎(chǔ),通過多組分?jǐn)?shù)值方法對(duì)不同飛行高度下催化/非催化壁面條件氣動(dòng)加熱的影響和加熱機(jī)理展開研究。

    1 研究方法

    1.1 數(shù)值方法

    對(duì)于球錐這種簡(jiǎn)單外形的非平衡可壓縮流動(dòng),國(guó)內(nèi)外研究中一般采用多組分N-S方程[2,3]進(jìn)行計(jì)算,本文針對(duì)壁面催化問題選用了5組分17基元反應(yīng)模型。在零攻角情況下球錐流場(chǎng)是軸對(duì)稱的,因此可以用二維軸對(duì)稱方程來模擬。軸對(duì)稱的N-S控制方程為

    在計(jì)算中控制方程的離散采用基于密度的有限體積法。對(duì)流項(xiàng)的離散采用AUSMPW+格式[4]。該格式在高超聲速的模擬中有廣泛的應(yīng)用,其計(jì)算量較??;而與矢通量分裂格式(Flux Vector Splitting,F(xiàn)VS)相比,耗散較低,對(duì)激波間斷和邊界層的捕捉有一定優(yōu)勢(shì)。本文應(yīng)用該格式時(shí),將通量分為對(duì)流項(xiàng)和壓力通量項(xiàng),根據(jù)馬赫數(shù)來確定分裂的形式。在離散方程中離散界面通量通過左、右兩側(cè)的物理量來獲得,但聲速選擇相同的值。時(shí)間積分采用隱式LU-SGS方法來加快收斂速度。

    1.2 熱化學(xué)反應(yīng)模型

    本文中化學(xué)反應(yīng)模型選擇Dung & Kang的高溫空氣化學(xué)反應(yīng)模型[5],用到其中5組分17基元反應(yīng)模型,如表1所示。高溫情況下流場(chǎng)中會(huì)出現(xiàn)振動(dòng)能激發(fā)的現(xiàn)象,本文采用Park[6]的雙溫度模型,雙溫度模型的具體細(xì)節(jié)可以參閱文獻(xiàn)[6]。

    表1 Dung & Kang 高溫空氣化學(xué)反應(yīng)模型

    Tab.1 Dung & Kang's Reaction Model for High-temperature Air

    序號(hào)反應(yīng)方程式序號(hào)反應(yīng)方程式 1O2+N=2O+N10N+N2=2N+N 2O2+NO=2O+NO11O2+O=2O+O 3N2+O=2N+O12O2+O2=2O+O2 4N2+NO=2N+NO13O2+N2=2O+N2 5N2+O2=2N+O214N2+N2=2N+N2 6NO+O2=N+O+O215NO+O=N+2O 7NO+N2=N+O+N216NO+N=O+2N 8O+NO=N+O217NO+NO=N+O+NO 9O+N2=N+NO——

    在表1的反應(yīng)體系中,非平衡流動(dòng)下每個(gè)基元化學(xué)反應(yīng)都可以朝著正、逆方向進(jìn)行,但是正、逆反應(yīng)速率都是獨(dú)立的而不受平衡常數(shù)約束。正、逆反應(yīng)速率均采用三系數(shù)的Arrhenius定律,其中的指前因子、溫度系數(shù)和活化能,以及多組分下的輸運(yùn)系數(shù)(粘性系數(shù)、二元擴(kuò)散系數(shù)等)均在文獻(xiàn)[5]附表中列出。

    1.3 數(shù)值方法的驗(yàn)證

    本文采用=25 mm的球頭開展了多次試驗(yàn),獲取了試驗(yàn)氣流條件下模型駐點(diǎn)的溫度變化曲線以及熱流曲線,并與本程序的計(jì)算結(jié)果進(jìn)行了對(duì)比。地面試驗(yàn)在正向爆轟驅(qū)動(dòng)激波-膨脹風(fēng)洞JF-16E中展開,其結(jié)構(gòu)如圖1所示。該風(fēng)洞利用正向爆轟大功率驅(qū)動(dòng)技 術(shù)[7],以及擴(kuò)容環(huán)腔結(jié)構(gòu)進(jìn)一步改善驅(qū)動(dòng)性能[8],并應(yīng)用高速運(yùn)動(dòng)激波在激波管-膨脹加速管內(nèi)的非定常膨脹過程,實(shí)現(xiàn)高超聲速試驗(yàn)氣流。目前試驗(yàn)氣流速度可在6~10 km/s之間自由選擇。試驗(yàn)過程中首先在JF-16E開展了一系列典型模型的高超聲速測(cè)熱試驗(yàn)包括球頭、尖楔和尖錐(如圖1c所示)。通過上述試驗(yàn)表明試驗(yàn)設(shè)備具備產(chǎn)生離解氣體的能力。

    圖1 JF-16E的爆轟驅(qū)激波-膨脹風(fēng)洞以及典型模型試驗(yàn)(V∞=8.1 km/s)

    本研究開展了3個(gè)典型試驗(yàn)狀態(tài)的測(cè)熱試驗(yàn),試驗(yàn)狀態(tài)分別為0=35、45、47 MJ/kg,試驗(yàn)件前來流條件如表2所示。

    表2 來流條件

    Tab.2 Inflow Condition

    工況總焓/(MJ?kg-1)速度/(km?s-1)溫度/K壓力/Pa 1358.02390952.46 2458.92914607.95 3479.2950709.28

    測(cè)熱過程中使用了鎳鉻-康銅同軸熱電偶傳感器。同軸熱電偶是基于熱電效應(yīng)的溫度傳感器,利用塞貝克效應(yīng)將2種不同材料的導(dǎo)體或半導(dǎo)體A和B連接成一個(gè)閉合回路。結(jié)點(diǎn)的溫度可以根據(jù)輸出的電勢(shì)確定,然后通過表面溫度變化經(jīng)過積分計(jì)算得到表面熱流率,其直徑為1.4 mm。在球頭駐點(diǎn)中心打孔并安裝傳感器之后,熱電偶與模型表面一起打磨拋光并鍍銅膜導(dǎo)通內(nèi)外2種材料構(gòu)成回路;鍍膜前后球頭表面的變化如圖1d所示。

    計(jì)算中分別使用了全催化和非催化2種壁面化學(xué)模型,結(jié)果如表3所示。試驗(yàn)中測(cè)量得到了高達(dá) 37 MJ/m2的超高熱流。對(duì)比結(jié)果表明,非催化壁面條件的計(jì)算結(jié)果與試驗(yàn)結(jié)果差別很大,而全催化壁面模型條件下的計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合很好,在總焓 35 MJ/kg條件的兩者差別小于10%,而總焓47 MJ/kg情況下兩者差別接近17%。試驗(yàn)?zāi)P统跗诒砻驽冦~膜,可以認(rèn)為是全催化壁面。通過上述對(duì)比結(jié)果可知,雖然數(shù)值方法直接給出熱流數(shù)據(jù)比較困難,但仍然能處于合理區(qū)間內(nèi)。

    表3 JF-16E高超聲速測(cè)熱試驗(yàn)與計(jì)算數(shù)據(jù)對(duì)比

    Tab.3 Comparison of Computation and Hypersonic Measurement on JF-16E

    序號(hào)總焓MJ/kg測(cè)熱試驗(yàn)熱流MJ/m2Fay-Riddell公式數(shù)值計(jì)算熱流 熱流MJ/m2誤差全催化MJ/m2誤差非催化MJ/m2 13521.3922.203.8%23.178.3%10.05 221.921.3%5.7% 34532.5430.067.6%29.589.1%13.28 44737.3933.879.4%31.1016.8%14.23

    2 壁面催化條件影響研究

    2.1 典型外形與工況

    計(jì)算模型外形選取ELECTRE頭部縮比的球錐模型,頭部半徑為35 mm,總長(zhǎng)為70 mm,錐角為4.6°,如圖2a所示。計(jì)算中采用的網(wǎng)格為100×100,壁面處首層網(wǎng)格為0.001 mm,網(wǎng)格高度沿法向以1.05的比例遞增,如圖2b所示。

    圖2 球錐模型以及計(jì)算網(wǎng)格

    2.2 40 km高空馬赫數(shù)10~30氣動(dòng)熱規(guī)律

    來流條件選擇40 km高空,飛行馬赫數(shù)分別為10、20、25、30,計(jì)算得到模型壁面的熱流分布如圖3所示。從圖3可以看出,模型頭部駐點(diǎn)區(qū)域?yàn)闊崃鞣逯?,沿著壁面向尾部熱流逐漸減小,對(duì)于球錐模型,頭部球形區(qū)域?yàn)橹饕母邿崃鲄^(qū)域,隨著飛行馬赫數(shù)提高,熱流值變大。對(duì)比研究2種催化條件對(duì)熱流的影響,=10時(shí),2種催化條件得到的熱流分布差別不明顯,隨著馬赫數(shù)增加,2種催化條件差別逐漸增大,全催化高于非催化條件熱流值。

    圖3 全催化&非催化條件下壁面熱流分布

    續(xù)圖3

    提取5組馬赫數(shù)工況下駐點(diǎn)熱流值,如圖4所示。圖4驗(yàn)證了之前的分析,在低馬赫數(shù)時(shí),2種催化條件下駐點(diǎn)熱流率幾乎重合,但在高馬赫數(shù)時(shí)差異明顯。

    圖4 Ma=10~30駐點(diǎn)熱流

    全催化與非催化條件下駐點(diǎn)熱流的比值如圖5所示。圖5中,全催化條件在馬赫數(shù)較高時(shí),對(duì)壁面熱流的影響有顯著增加的效應(yīng),且隨著馬赫數(shù)不斷提高,這種效應(yīng)在加劇。分析原因認(rèn)為,隨著飛行馬赫數(shù)的提高,頭部激波的強(qiáng)度不斷提高,導(dǎo)致波后氣體離解程度加大,離解出的原子成分比重增加,造成壁面處有更多的原子進(jìn)行催化復(fù)合釋熱,因此全催化條件的效應(yīng)得到加強(qiáng),對(duì)熱流的提升也更為顯著。

    采用SPSS 22.0軟件對(duì)數(shù)據(jù)進(jìn)行分析處理,牙周袋探診深度、臨床附著水平、牙齦指數(shù)、牙齒松動(dòng)度、牙槽骨高度等指標(biāo)為計(jì)量資料,采用(±s)表示,滿足正態(tài)分布和方差齊性后,組間相較采用獨(dú)立t檢驗(yàn),以P<0.05表示差異具有統(tǒng)計(jì)學(xué)意義。

    圖5 駐點(diǎn)熱流率比值(全催化/非催化)

    提取=20時(shí)駐點(diǎn)線氣流粒子質(zhì)量分?jǐn)?shù)分布,如圖6所示。

    圖6 Ma=20駐點(diǎn)線組分質(zhì)量分?jǐn)?shù)

    波后O2基本完全離解成O,N2也有較大的離解程度,在壁面附近,全催化條件下原子完全復(fù)合成分子狀態(tài),非催化條件則保留了一部分原子成分。表4為非催化條件駐點(diǎn)處粒子質(zhì)量分?jǐn)?shù)表。

    表4 非催化條件駐點(diǎn)粒子質(zhì)量分?jǐn)?shù)表

    Tab.4 Mass-fraction of Particles on Stagnant Point in Non-catalytic Condition

    質(zhì)量分?jǐn)?shù)馬赫數(shù) 1015202530 N20.76490.75690.69180.67090.6534 O20.22580.03250.02290.02230.0145 N0.00030.00630.07560.10690.1449 O0.00840.20180.20760.19740.1855 NO0.00060.00260.00210.00250.0016

    由表4可以看出,馬赫數(shù)較低時(shí),離解效應(yīng)微弱,因此壁面催化條件對(duì)熱流的影響不大。隨著馬赫數(shù)提高,駐點(diǎn)處的原子質(zhì)量比重增加,而相對(duì)應(yīng)的全催化條件下原子質(zhì)量比重為0,2種催化條件的差異在不斷加大,催化效應(yīng)對(duì)熱流的影響隨馬赫數(shù)增加而加劇。

    2.3 馬赫數(shù)為20的40~70 km高空氣動(dòng)熱規(guī)律

    選擇固定飛行馬赫數(shù)20,來流條件分布選擇 40 km、50 km、60 km、70 km高空氣流參數(shù),模型外形不變,研究不同高空氣流參數(shù)對(duì)氣動(dòng)熱的影響規(guī)律。

    圖7分別為高度為50 km、60 km、70 km時(shí)2種催化條件下壁面熱流的結(jié)果,并加入2.2節(jié)中已經(jīng)得到的40 km結(jié)果一起進(jìn)行比較分析。

    圖7 全催化&非催化條件下壁面熱流分布

    由圖7可以看出,隨著飛行高度的增加,熱流值減小。對(duì)比催化條件對(duì)熱流的影響,70 km高空2種催化條件對(duì)熱流的影響不大,壁面熱流曲線幾乎重合,但是在其余3組高度均表現(xiàn)出顯著的差異。提取4組工況駐點(diǎn)熱流如圖8所示,并計(jì)算得到全催化對(duì)比非催化條件駐點(diǎn)熱流比值如圖9所示。

    圖8 40~80 km高空駐點(diǎn)熱流

    圖9 駐點(diǎn)熱流率比值(全催化/非催化)

    由圖8可知,隨著飛行高度增加,2種催化條件駐點(diǎn)熱流的絕對(duì)差值減小,在70 km高空處趨近于重合。圖9中駐點(diǎn)熱流的比值線出現(xiàn)一個(gè)拐點(diǎn),40 km的值小于50 km,50 km高空以上該比值隨高度增加遞減。分析具體原因,提取50 km、70 km高空駐點(diǎn)線氣流粒子組分質(zhì)量分?jǐn)?shù),如圖10所示。

    a)50 km

    b)70 km

    續(xù)圖10

    表5 非催化條件壁面粒子質(zhì)量百分比表

    Tab.5 Mass-fraction of Particles on Wall Surface in Non-catalytic Condition

    質(zhì)量分?jǐn)?shù)高度/km 4050 60 70 N20.69180.61480.69250.7608 O20.02290.00400.00840.2124 N0.07560.15470.08180.0068 O0.20760.22450.21440.0188 NO0.00210.00200.00300.0011

    從表5中可以看出,40 km原子組分質(zhì)量比重小于50 km,50 km高空以上原子組分質(zhì)量比重隨高度增加而減小,驗(yàn)證了之前的分析。

    3 結(jié) 論

    本文首先利用地面試驗(yàn)開展小尺度模型高焓流動(dòng)條件下的非平衡氣動(dòng)熱測(cè)量,并驗(yàn)證了本文所用的數(shù)值算法。然后開展高超聲速條件下全催化/非催化壁面的氣動(dòng)熱環(huán)境影響研究,在研究中發(fā)現(xiàn)高焓條件下壁面催化條件對(duì)熱流影響比較大。

    a)在高馬赫數(shù)下,激波后的分子離解更強(qiáng),原子質(zhì)量分?jǐn)?shù)占比較高,因此在壁面復(fù)合的時(shí)候產(chǎn)生了較多的熱流。

    b)來流密度隨高度減小導(dǎo)致在較高的高度上粒子碰撞不足,更多的原子運(yùn)動(dòng)到壁面時(shí)才復(fù)合,而低空中原子在流動(dòng)過程中就得到足夠的碰撞而復(fù)合。而當(dāng)氣體過于稀薄時(shí)(粒子碰撞數(shù)更少)導(dǎo)致分子離解很少,因此基本不受壁面條件影響。

    [1] Bertin J J, Cummings R M. Critical hypersonic aerothermodynamics phenomena[J]. Annual Review of Fluid Mechanics, 2006(38): 129-157.

    [2] 潘沙. 高超聲速氣動(dòng)熱數(shù)值模擬方法及大規(guī)模并行計(jì)算研究[D]. 長(zhǎng)沙:國(guó)防科技大學(xué), 2010.

    Pan Sha. Hypersonic aerothermal numerical simulation method and massive parallel computation research[D]. Changsha: National University of Defense Technology, 2010.

    [3] 李康, 等. 高焓熱化學(xué)非平衡流動(dòng)中再壓縮效應(yīng)對(duì)球錐氣動(dòng)力的影 響[J]. 力學(xué)學(xué)報(bào), 2015, 47(2): 346-350.

    Li Kang, et al. The effect of re-compression to aerodynamic force in hypersonic thermo-chemical nonequilibrium flow[J].Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(2): 346-350.

    [4] Kim Kyu Hong, Kim Chongam, Rho Oh-Hyun. Methods for the accurate computations of hypersonic flows I[J]. AUSMPW+ SchemeJournal of Computational Physics, 2001(174): 38-80.

    [5] Gnoffo P A, Gupta R N, Shinn J L. Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium[R]. NASA TP-2867, 1989.

    [6] Park C. Assessment of two-temperature kinetics model for dissociation and weakly-inozing nitrogen[J]. Journal of Thermophysics And Heat Transfer, 1988, 2(1): 8-16.

    [7] 俞鴻儒, 趙偉. 氫氧爆轟驅(qū)動(dòng)激波風(fēng)洞性能[J]. 氣動(dòng)實(shí)驗(yàn)與測(cè)量控制, 1993, 7(3): 3, 38-42.

    Yu Hongru, Zhao Wei. Performance of shock tunnel with H2-O2detonation driver[J].Aerodynamic Experiment and Measurement & Control, 1993, 7(3): 3, 38-42.

    [8] 姜宗林, 等.長(zhǎng)試驗(yàn)時(shí)間爆轟驅(qū)動(dòng)激波風(fēng)洞技術(shù)研究[J]. 力學(xué)學(xué)報(bào), 2012, 44(5): 824-831.

    Jiang Zonglin, et al. Approaches to extend the test time for Detonation-driven Shock Tunnel[J].Chinese Journal of Theoretical and Applied Mechanics, 2012, 44(5): 824-831.

    Research on the Influence and Mechanisms of Catalytic Wall Condition on Aerodynamic-Heating

    Li Xu-dong1,2, Zhang Fu1,2, Shi Zeng-min1,2, Wang Xin1,2

    (1. Beijing Institute of Space Long March Vehicle, Beijing, 100076; 2. Hypersonic Vehicle Research Center of Thermal Protection and Insulation, Beijing, 100076)

    To investigate the influence of catalytic wall condition to aerodynamic-heating,the effects of different wall conditions (Fully catalytic wall/ Non-catalytic wall) on aerodynamic heating for hypersonic thermal-chemical non-equilibrium flow are studied. Numerical method, which is constructed according to the thermal-chemical non-equilibrium flow characteristics and verified by high enthalpy & hypersonic ground experiment results, are used to investigate the influence of different catalytic wall conditions on aerodynamic thermal environment and its intrinsic mechanism for some typical geometrical configurations and different flight heights. Study shows that there’s a significant catalytic effect when flight height is above 50km, and heat-flux under FCW/NCW differs apparently. But the difference between heat-flux und FCW/NCW reduces with flight altitude decreasing.

    hypersonic flow; real-gas effect; catalytic wall; numerical method

    TP333

    A

    1004-7182(2020)01-0112-06

    10.7654/j.issn.1004-7182.20200120

    2018-06-28;

    2018-09-22

    李旭東(1983-),男,博士,高級(jí)工程師,主要研究方向?yàn)楦咚倩瘜W(xué)反應(yīng)流動(dòng)與氣動(dòng)熱。

    張 賦(1992-),男,博士,工程師,主要研究方向?yàn)榧げㄏ嗷プ饔靡约霸囼?yàn)流體力學(xué)。

    史增民(1970-),男,研究員,主要研究方向?yàn)闅鈩?dòng)熱環(huán)境。

    王 欣(1992-),男,工程師,主要研究方向?yàn)闅鈩?dòng)熱環(huán)境。

    猜你喜歡
    駐點(diǎn)馬赫數(shù)激波
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    載荷分布對(duì)可控?cái)U(kuò)散葉型性能的影響
    一種基于聚類分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    基于游人游賞行為的留園駐點(diǎn)分布規(guī)律研究
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    利用遠(yuǎn)教站點(diǎn),落實(shí)駐點(diǎn)干部帶學(xué)
    利用遠(yuǎn)教站點(diǎn),落實(shí)駐點(diǎn)干部帶學(xué)
    videos熟女内射| 精品视频人人做人人爽| 日韩制服骚丝袜av| 99视频精品全部免费 在线| 亚洲成人一二三区av| 午夜老司机福利剧场| 各种免费的搞黄视频| freevideosex欧美| 人人妻人人澡人人爽人人夜夜| 亚洲第一av免费看| 色婷婷久久久亚洲欧美| 91午夜精品亚洲一区二区三区| 免费人成在线观看视频色| 香蕉精品网在线| 国产精品久久久久久精品电影小说| 一级毛片久久久久久久久女| av免费观看日本| 纵有疾风起免费观看全集完整版| 黄色一级大片看看| 久久国产精品男人的天堂亚洲 | 国产av一区二区精品久久| 久热久热在线精品观看| 亚洲va在线va天堂va国产| 少妇裸体淫交视频免费看高清| 亚洲精品色激情综合| 国产精品秋霞免费鲁丝片| 男女免费视频国产| 国产白丝娇喘喷水9色精品| 蜜桃在线观看..| 国产爽快片一区二区三区| 国产 精品1| av卡一久久| 亚洲欧洲精品一区二区精品久久久 | 超碰97精品在线观看| 插阴视频在线观看视频| 国产视频首页在线观看| 国产亚洲91精品色在线| 蜜臀久久99精品久久宅男| 亚洲中文av在线| 美女中出高潮动态图| 蜜桃久久精品国产亚洲av| 狂野欧美激情性bbbbbb| 久久国产亚洲av麻豆专区| 国国产精品蜜臀av免费| 七月丁香在线播放| 亚洲熟女精品中文字幕| 亚洲欧美精品自产自拍| 国产探花极品一区二区| 免费看不卡的av| av线在线观看网站| 国产视频内射| 国产片特级美女逼逼视频| 亚洲,一卡二卡三卡| 各种免费的搞黄视频| 街头女战士在线观看网站| 高清黄色对白视频在线免费看 | 亚洲内射少妇av| 又大又黄又爽视频免费| 亚洲av欧美aⅴ国产| 免费看光身美女| 日韩欧美一区视频在线观看 | 22中文网久久字幕| 美女大奶头黄色视频| 看免费成人av毛片| 国产淫语在线视频| 国产欧美日韩一区二区三区在线 | 2018国产大陆天天弄谢| 午夜av观看不卡| 肉色欧美久久久久久久蜜桃| 97超碰精品成人国产| 少妇精品久久久久久久| 九九爱精品视频在线观看| 国产精品久久久久久精品古装| 99热6这里只有精品| 插阴视频在线观看视频| 一级片'在线观看视频| av免费观看日本| 亚洲中文av在线| 秋霞伦理黄片| 97在线人人人人妻| 黑人巨大精品欧美一区二区蜜桃 | 欧美97在线视频| 久久久久久人妻| 午夜激情福利司机影院| 在线 av 中文字幕| 日韩免费高清中文字幕av| 一级二级三级毛片免费看| 精品国产乱码久久久久久小说| 久久久欧美国产精品| 国产无遮挡羞羞视频在线观看| 最黄视频免费看| 三级经典国产精品| 国产伦精品一区二区三区四那| 超碰97精品在线观看| 一本久久精品| 日本av手机在线免费观看| 男女国产视频网站| 一级毛片 在线播放| 高清午夜精品一区二区三区| 久久狼人影院| 熟妇人妻不卡中文字幕| 狂野欧美激情性bbbbbb| 国产成人a∨麻豆精品| 青青草视频在线视频观看| 少妇的逼好多水| 日韩制服骚丝袜av| 日韩精品有码人妻一区| 一区在线观看完整版| 这个男人来自地球电影免费观看 | 日本猛色少妇xxxxx猛交久久| 菩萨蛮人人尽说江南好唐韦庄| 午夜日本视频在线| 寂寞人妻少妇视频99o| 精品国产露脸久久av麻豆| 99国产精品免费福利视频| 亚洲精品一二三| kizo精华| 国产精品人妻久久久久久| 国产免费又黄又爽又色| 亚洲图色成人| 国产无遮挡羞羞视频在线观看| 热99国产精品久久久久久7| 美女福利国产在线| 青春草亚洲视频在线观看| 欧美少妇被猛烈插入视频| 成人国产麻豆网| 亚洲国产av新网站| 欧美日韩精品成人综合77777| 日日啪夜夜撸| 亚洲熟女精品中文字幕| 久久99蜜桃精品久久| 亚洲激情五月婷婷啪啪| 国产精品国产三级专区第一集| 久久久亚洲精品成人影院| 日韩av不卡免费在线播放| 大陆偷拍与自拍| 少妇高潮的动态图| 国产欧美亚洲国产| 伦理电影大哥的女人| 亚洲成人av在线免费| 久久女婷五月综合色啪小说| 三级经典国产精品| 亚洲成色77777| 精品酒店卫生间| 久久久国产精品麻豆| 人妻一区二区av| 国产在线免费精品| h日本视频在线播放| 黄色一级大片看看| 午夜福利,免费看| 精品人妻熟女毛片av久久网站| 国产伦在线观看视频一区| 亚洲欧美日韩东京热| 久久久久久久久久久免费av| 三级经典国产精品| 插逼视频在线观看| 日韩av不卡免费在线播放| 亚洲真实伦在线观看| av播播在线观看一区| 成人黄色视频免费在线看| a级毛片在线看网站| 久久久亚洲精品成人影院| 日产精品乱码卡一卡2卡三| 久久热精品热| 国产精品偷伦视频观看了| 久久99精品国语久久久| 成人毛片a级毛片在线播放| 久久久久久久久久久久大奶| 国产成人精品久久久久久| 男人舔奶头视频| 最黄视频免费看| 精品一区二区免费观看| 99久久精品热视频| 亚洲av男天堂| 十分钟在线观看高清视频www | 精品一区在线观看国产| 国产成人精品婷婷| 日韩强制内射视频| 99久久综合免费| 一区二区三区免费毛片| 一本大道久久a久久精品| 最近手机中文字幕大全| 亚洲美女视频黄频| 亚洲av成人精品一二三区| 国产欧美另类精品又又久久亚洲欧美| 欧美少妇被猛烈插入视频| 免费av中文字幕在线| 国产日韩一区二区三区精品不卡 | 日本免费在线观看一区| 69精品国产乱码久久久| 国产一区二区在线观看av| 18禁裸乳无遮挡动漫免费视频| 免费av中文字幕在线| 啦啦啦中文免费视频观看日本| 国产精品久久久久久av不卡| 99久久精品国产国产毛片| 国产 精品1| 另类精品久久| av在线app专区| 欧美一级a爱片免费观看看| 97超碰精品成人国产| 久久久久久久久久成人| 大陆偷拍与自拍| av网站免费在线观看视频| 观看av在线不卡| 日本vs欧美在线观看视频 | 国产视频内射| 久久精品国产鲁丝片午夜精品| 亚洲第一av免费看| 人妻一区二区av| 久久女婷五月综合色啪小说| 一级,二级,三级黄色视频| 少妇人妻一区二区三区视频| 色吧在线观看| 日本爱情动作片www.在线观看| 热re99久久国产66热| 免费黄色在线免费观看| 亚洲精品国产av蜜桃| 亚洲av在线观看美女高潮| 人妻 亚洲 视频| 国产探花极品一区二区| 免费观看av网站的网址| 亚洲激情五月婷婷啪啪| 国内揄拍国产精品人妻在线| 亚洲精品乱久久久久久| 搡老乐熟女国产| 在线观看美女被高潮喷水网站| 亚洲欧美日韩卡通动漫| 久久久久网色| 在线免费观看不下载黄p国产| 精品99又大又爽又粗少妇毛片| 丝袜喷水一区| 人妻一区二区av| 99视频精品全部免费 在线| 国产成人精品无人区| 最黄视频免费看| 一级毛片电影观看| 国产免费视频播放在线视频| 日韩欧美精品免费久久| a级毛色黄片| 久久精品国产鲁丝片午夜精品| 国产av一区二区精品久久| 人人妻人人爽人人添夜夜欢视频 | 夜夜看夜夜爽夜夜摸| 最新的欧美精品一区二区| 少妇丰满av| 一级毛片 在线播放| 成人特级av手机在线观看| 国产精品女同一区二区软件| 大片免费播放器 马上看| 51国产日韩欧美| 欧美 日韩 精品 国产| 夜夜爽夜夜爽视频| 久久久久国产网址| 国产精品一区二区性色av| 日本与韩国留学比较| 最新中文字幕久久久久| 亚洲国产欧美在线一区| 多毛熟女@视频| 久久韩国三级中文字幕| 久久久久国产精品人妻一区二区| 亚洲av在线观看美女高潮| 欧美变态另类bdsm刘玥| 亚洲欧美日韩东京热| 国产成人精品久久久久久| 国内少妇人妻偷人精品xxx网站| 亚洲精品亚洲一区二区| 日韩三级伦理在线观看| 天美传媒精品一区二区| 亚洲欧美清纯卡通| 国产无遮挡羞羞视频在线观看| 亚洲国产毛片av蜜桃av| 国产又色又爽无遮挡免| 欧美最新免费一区二区三区| 久久99一区二区三区| 日韩中字成人| av一本久久久久| 国产亚洲精品久久久com| 丝袜在线中文字幕| 大又大粗又爽又黄少妇毛片口| 国产一区有黄有色的免费视频| 免费大片黄手机在线观看| 亚洲精品日韩在线中文字幕| 免费观看a级毛片全部| 国产一区亚洲一区在线观看| 一区二区三区四区激情视频| 日韩一区二区三区影片| av在线app专区| 免费看光身美女| 精品酒店卫生间| 日韩,欧美,国产一区二区三区| 我的老师免费观看完整版| 草草在线视频免费看| 中文字幕av电影在线播放| 午夜激情福利司机影院| 久久久a久久爽久久v久久| 亚洲经典国产精华液单| 国产 一区精品| 一级二级三级毛片免费看| 97精品久久久久久久久久精品| 另类亚洲欧美激情| 伦理电影大哥的女人| 在线观看免费日韩欧美大片 | 男人舔奶头视频| 中国三级夫妇交换| 国精品久久久久久国模美| 在线观看www视频免费| 精品视频人人做人人爽| 在线天堂最新版资源| 亚洲欧美成人综合另类久久久| 男女国产视频网站| 国产精品嫩草影院av在线观看| 国产爽快片一区二区三区| 色5月婷婷丁香| 亚洲av国产av综合av卡| 中文精品一卡2卡3卡4更新| 好男人视频免费观看在线| 国产日韩欧美视频二区| 一级毛片aaaaaa免费看小| 人妻夜夜爽99麻豆av| 少妇人妻久久综合中文| 国产av国产精品国产| 妹子高潮喷水视频| 国产伦理片在线播放av一区| 久久久久久久久久久久大奶| 国产国拍精品亚洲av在线观看| 春色校园在线视频观看| 99热这里只有精品一区| 国产成人一区二区在线| 亚洲av国产av综合av卡| 成人18禁高潮啪啪吃奶动态图 | 成人影院久久| 一本大道久久a久久精品| av黄色大香蕉| 日本免费在线观看一区| 中文精品一卡2卡3卡4更新| 简卡轻食公司| 日韩一区二区视频免费看| 日本-黄色视频高清免费观看| 我的女老师完整版在线观看| 插阴视频在线观看视频| 成人美女网站在线观看视频| 永久网站在线| 各种免费的搞黄视频| 男人和女人高潮做爰伦理| 国产成人精品一,二区| 在线看a的网站| 欧美bdsm另类| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 最近中文字幕2019免费版| 国产成人免费观看mmmm| 五月天丁香电影| 日韩一区二区视频免费看| 久久热精品热| 在线观看国产h片| 啦啦啦啦在线视频资源| 午夜视频国产福利| 久久婷婷青草| 99九九线精品视频在线观看视频| 三级国产精品欧美在线观看| 国产精品一区二区三区四区免费观看| 另类精品久久| 国产乱人偷精品视频| 亚洲成人av在线免费| 极品人妻少妇av视频| 久久综合国产亚洲精品| 亚洲成人一二三区av| 国产精品不卡视频一区二区| 国产精品福利在线免费观看| 欧美日韩视频高清一区二区三区二| 极品教师在线视频| 久久久国产精品麻豆| 一区二区av电影网| 中文字幕av电影在线播放| 狠狠精品人妻久久久久久综合| 女性生殖器流出的白浆| 国产成人91sexporn| 最后的刺客免费高清国语| 中文字幕制服av| 三级经典国产精品| 亚洲一级一片aⅴ在线观看| 日韩免费高清中文字幕av| 男人添女人高潮全过程视频| 亚洲精品久久久久久婷婷小说| 人人妻人人添人人爽欧美一区卜| 亚洲av不卡在线观看| 国产黄片美女视频| 九九久久精品国产亚洲av麻豆| av国产久精品久网站免费入址| 在线观看免费视频网站a站| 亚洲激情五月婷婷啪啪| 日韩强制内射视频| 色视频www国产| 伦理电影大哥的女人| 亚州av有码| 简卡轻食公司| 久久久久久久久久成人| 国产 一区精品| 男男h啪啪无遮挡| 精品亚洲乱码少妇综合久久| 国产av精品麻豆| 一级a做视频免费观看| 国产亚洲一区二区精品| kizo精华| 美女主播在线视频| 午夜91福利影院| 在线观看一区二区三区激情| 边亲边吃奶的免费视频| 多毛熟女@视频| 国产精品蜜桃在线观看| 插逼视频在线观看| 免费久久久久久久精品成人欧美视频 | 中文精品一卡2卡3卡4更新| 男人爽女人下面视频在线观看| 国产成人午夜福利电影在线观看| 国产极品天堂在线| 中文字幕av电影在线播放| 精品亚洲乱码少妇综合久久| 大陆偷拍与自拍| 国产精品国产三级国产av玫瑰| 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 国产精品秋霞免费鲁丝片| 一区二区三区免费毛片| 久久狼人影院| 亚洲一区二区三区欧美精品| 一二三四中文在线观看免费高清| 伊人久久国产一区二区| 制服丝袜香蕉在线| 久久影院123| 国产无遮挡羞羞视频在线观看| 十八禁网站网址无遮挡 | 免费黄色在线免费观看| 亚洲欧美日韩东京热| 美女国产视频在线观看| 下体分泌物呈黄色| 人人妻人人爽人人添夜夜欢视频 | 777米奇影视久久| 国产精品熟女久久久久浪| 丝袜喷水一区| 亚洲av中文av极速乱| 久久国产亚洲av麻豆专区| 午夜免费观看性视频| 亚洲精品日韩av片在线观看| 国产在视频线精品| a级一级毛片免费在线观看| 少妇人妻一区二区三区视频| 国产精品久久久久成人av| 亚洲精华国产精华液的使用体验| 国产精品偷伦视频观看了| 美女福利国产在线| av视频免费观看在线观看| 国产黄频视频在线观看| 最近手机中文字幕大全| 亚洲国产av新网站| 久久ye,这里只有精品| 国产爽快片一区二区三区| 我的老师免费观看完整版| 国产精品女同一区二区软件| 又爽又黄a免费视频| 日韩av在线免费看完整版不卡| av卡一久久| 秋霞伦理黄片| av免费观看日本| 中文乱码字字幕精品一区二区三区| 麻豆成人av视频| 777米奇影视久久| 99热这里只有是精品在线观看| 国产男女超爽视频在线观看| 成人特级av手机在线观看| 内地一区二区视频在线| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 在线播放无遮挡| 日韩中文字幕视频在线看片| 亚洲国产成人一精品久久久| 大话2 男鬼变身卡| 久久久久久久精品精品| 国产色爽女视频免费观看| 边亲边吃奶的免费视频| 91精品伊人久久大香线蕉| 18+在线观看网站| 最近最新中文字幕免费大全7| 狂野欧美白嫩少妇大欣赏| 国产av精品麻豆| 午夜免费鲁丝| 大又大粗又爽又黄少妇毛片口| 最新的欧美精品一区二区| 亚洲国产最新在线播放| 久久亚洲国产成人精品v| 777米奇影视久久| 亚洲精品乱码久久久v下载方式| 亚洲精品自拍成人| 插阴视频在线观看视频| 国产伦在线观看视频一区| 国产精品国产三级国产av玫瑰| 精品人妻熟女av久视频| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 免费看日本二区| 两个人的视频大全免费| 在线观看免费高清a一片| 成人二区视频| 亚洲国产精品一区二区三区在线| 日韩av免费高清视频| 少妇 在线观看| 天堂中文最新版在线下载| 毛片一级片免费看久久久久| 国产成人午夜福利电影在线观看| 69精品国产乱码久久久| 丝袜脚勾引网站| 欧美日韩综合久久久久久| 亚洲精品久久午夜乱码| 精品一品国产午夜福利视频| 亚洲av.av天堂| 亚洲精品456在线播放app| 国产永久视频网站| 久久午夜福利片| 亚洲,欧美,日韩| 国国产精品蜜臀av免费| 久久精品国产鲁丝片午夜精品| 亚洲欧美精品专区久久| 高清欧美精品videossex| 黄色一级大片看看| 精品一区在线观看国产| 中文字幕亚洲精品专区| 国产永久视频网站| 欧美日韩国产mv在线观看视频| 欧美+日韩+精品| 免费av中文字幕在线| 熟妇人妻不卡中文字幕| 久久久午夜欧美精品| 中文乱码字字幕精品一区二区三区| 熟女电影av网| 精品少妇内射三级| 亚洲精品一二三| 欧美国产精品一级二级三级 | 日韩伦理黄色片| 啦啦啦啦在线视频资源| 久久免费观看电影| 99久国产av精品国产电影| 亚洲成色77777| 啦啦啦中文免费视频观看日本| 91精品一卡2卡3卡4卡| 少妇猛男粗大的猛烈进出视频| 亚洲,欧美,日韩| 久久ye,这里只有精品| 赤兔流量卡办理| 欧美少妇被猛烈插入视频| 亚洲图色成人| 人妻夜夜爽99麻豆av| 国语对白做爰xxxⅹ性视频网站| 国产精品麻豆人妻色哟哟久久| 老司机影院毛片| 国产综合精华液| 欧美97在线视频| 国产男人的电影天堂91| 欧美日韩亚洲高清精品| 少妇人妻一区二区三区视频| 一级av片app| 国产亚洲5aaaaa淫片| 免费在线观看成人毛片| 日本免费在线观看一区| 卡戴珊不雅视频在线播放| 人妻一区二区av| 99九九线精品视频在线观看视频| 中文资源天堂在线| 免费大片18禁| 亚洲,欧美,日韩| 久久99热这里只频精品6学生| 成年人免费黄色播放视频 | 青春草国产在线视频| 在线观看av片永久免费下载| a 毛片基地| 一级,二级,三级黄色视频| 精品人妻熟女av久视频| 久久热精品热| 97精品久久久久久久久久精品| 日韩av不卡免费在线播放| 国产日韩欧美视频二区| 国产精品无大码| 亚洲国产欧美日韩在线播放 | 欧美3d第一页| 69精品国产乱码久久久| 晚上一个人看的免费电影| 欧美bdsm另类| 亚洲av在线观看美女高潮| 亚洲av.av天堂| 亚洲美女搞黄在线观看| 亚洲av.av天堂| 夜夜爽夜夜爽视频| 亚洲精品国产成人久久av| 人妻一区二区av| 久久综合国产亚洲精品| 少妇高潮的动态图| 日韩欧美精品免费久久| 成年美女黄网站色视频大全免费 | 老司机亚洲免费影院| 国产女主播在线喷水免费视频网站| 乱系列少妇在线播放| 国产亚洲5aaaaa淫片| 国产精品蜜桃在线观看| 亚洲av成人精品一二三区| 久久久国产欧美日韩av| 少妇的逼好多水| av视频免费观看在线观看| 一级二级三级毛片免费看| 亚洲在久久综合| 在线观看免费高清a一片| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人久久小说| 一级毛片黄色毛片免费观看视频| 啦啦啦中文免费视频观看日本| 精品亚洲成国产av|