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

    復(fù)燃對氫氧火箭發(fā)動(dòng)機(jī)尾焰流場及輻射特性影響數(shù)值研究

    2016-06-05 09:34:02聶萬勝豐松江吳高楊
    關(guān)鍵詞:尾焰氫氧流場

    喬 野,聶萬勝,豐松江,吳高楊

    ?

    復(fù)燃對氫氧火箭發(fā)動(dòng)機(jī)尾焰流場及輻射特性影響數(shù)值研究

    喬 野,聶萬勝,豐松江,吳高楊

    (中國人民解放軍裝備學(xué)院航天裝備系,北京,101416)

    為深入研究復(fù)燃對氫氧火箭發(fā)動(dòng)機(jī)尾焰流場及輻射特性的影響,以氫氧發(fā)動(dòng)機(jī)喉部截面參數(shù)為入口條件,采用耦合Realizable k-ε湍流模型的三維N-S方程,考慮尾焰復(fù)燃反應(yīng)影響,利用PISO算法求解得到尾焰流場參數(shù)。在此基礎(chǔ)上,通過氣體輻射傳輸方程和大氣透過率計(jì)算模型SLG對尾焰輻射特性進(jìn)行計(jì)算,對比復(fù)燃反應(yīng)對尾焰流場及其輻射特性的影響。結(jié)果表明,復(fù)燃反應(yīng)對氫氧發(fā)動(dòng)機(jī)尾焰流場計(jì)算影響較大,使溫度場以及燃燒產(chǎn)物的質(zhì)量分?jǐn)?shù)大幅增加,從而導(dǎo)致尾焰的輻射特性增強(qiáng),因而在氫氧發(fā)動(dòng)機(jī)尾焰流場和輻射計(jì)算中,考慮復(fù)燃反應(yīng)是極為必要的。

    氫氧火箭發(fā)動(dòng)機(jī);尾焰;復(fù)燃反應(yīng);紅外輻射;數(shù)值仿真

    0 引 言

    火箭發(fā)動(dòng)機(jī)的燃?xì)鈱儆诟蝗既細(xì)?,不能完全燃燒的燃?xì)鈬娙氪髿夂?,?huì)和空氣中的氧氣摻混燃燒發(fā)生復(fù)燃反應(yīng),從而影響尾焰的溫度場和燃?xì)饨M分分布[1]?;鸺l(fā)動(dòng)機(jī)尾焰具有高溫、高速、大流量的特點(diǎn),其在飛行階段高溫尾焰會(huì)產(chǎn)生強(qiáng)烈的輻射特性,這對實(shí)現(xiàn)紅外追蹤與預(yù)警具有重要意義。但是,火箭發(fā)動(dòng)機(jī)尾焰流動(dòng)是耦合強(qiáng)烈化學(xué)反應(yīng)的復(fù)雜流動(dòng),在計(jì)算中定量分析復(fù)燃對尾焰的影響極為必要。

    Leone等[2]研究了燃燒尾焰的化學(xué)成分在高空發(fā)生復(fù)燃反應(yīng)對同溫層產(chǎn)生的影響;王雁鳴等[3]對多噴管發(fā)動(dòng)機(jī)在低空的紅外輻射特性進(jìn)行研究,得到多噴管發(fā)動(dòng)機(jī)尾焰流場在2~5 μm波段光譜的紅外特性;聶萬勝等[4、5]對液體火箭尾部噴焰紅外輻射特性進(jìn)行數(shù)值仿真計(jì)算,并分析燃燒室內(nèi)一步反應(yīng)和兩步反應(yīng)、無燃燒室3種情形下尾焰流場與輻射特點(diǎn);王偉臣等[6]研究了固體火箭發(fā)動(dòng)機(jī)尾焰復(fù)燃的輻射效應(yīng),在復(fù)燃尾焰中耦合紅外輻射傳輸模型;劉尊洋等[7~9]針對復(fù)燃對固體火箭及液體火箭尾焰紅外輻射特性影響分別進(jìn)行研究,分析了不同飛行參數(shù)對火箭尾焰輻射特性的影響。上述文獻(xiàn)從不同角度對火箭發(fā)動(dòng)機(jī)尾焰流場及輻射特性進(jìn)行了研究,為火箭發(fā)動(dòng)機(jī)尾焰研究奠定了重要基礎(chǔ)。

    以液氫液氧作為推進(jìn)劑的液體火箭發(fā)動(dòng)機(jī),反應(yīng)能量遠(yuǎn)遠(yuǎn)大于液氧煤油發(fā)動(dòng)機(jī)和常規(guī)的偏二甲肼/四氧化二氮發(fā)動(dòng)機(jī)。同時(shí),由于氫的分子量極低,氫氧發(fā)動(dòng)機(jī)比推力比其他任何推進(jìn)劑發(fā)動(dòng)機(jī)的比推力都高,并具有無毒、無污染等特點(diǎn),因而也是國外爭相研究的對象,如美國的J-2發(fā)動(dòng)機(jī)、日本的LE-7/A發(fā)動(dòng)機(jī)及蘇聯(lián)的RD-0120發(fā)動(dòng)機(jī)[10]等。目前,針對復(fù)燃反應(yīng)對氫氧發(fā)動(dòng)機(jī)尾焰流場及輻射特性的影響還不完善,需要進(jìn)一步研究。

    本文以氫氧發(fā)動(dòng)機(jī)為研究對象,對發(fā)動(dòng)機(jī)內(nèi)流場以及尾焰外流場進(jìn)行計(jì)算,并以內(nèi)流場喉部截面參數(shù)作為尾焰外流場計(jì)算的入口條件,從而達(dá)到簡化計(jì)算的目的。尾焰輻射特性以尾焰流場參數(shù)為基礎(chǔ),并利用氣體輻射傳輸方程和大氣透過率計(jì)算模型SLG對尾焰輻射特性進(jìn)行計(jì)算。通過以上方法,可以得到氫氧發(fā)動(dòng)機(jī)內(nèi)流場、尾焰外流場以及輻射特性參數(shù),從而定量分析復(fù)燃反應(yīng)對尾焰流場和輻射特性的影響。

    1 物理模型與計(jì)算方法

    1.1 幾何模型

    氫氧發(fā)動(dòng)機(jī)結(jié)構(gòu)具有對稱性,可以取圓周面的三分之一進(jìn)行對稱處理,從而簡化計(jì)算。發(fā)動(dòng)機(jī)幾何模型如圖1所示。發(fā)動(dòng)機(jī)主要由短噴管推力室和噴管延伸段組成,如圖1a。由于在發(fā)動(dòng)機(jī)內(nèi)流動(dòng)處于超聲速流動(dòng)狀態(tài),流動(dòng)擾動(dòng)無法逆向傳播,因而可通過只計(jì)算短噴管推力室的內(nèi)流場,獲取喉部截面參數(shù),提高計(jì)算效率。

    a)整體

    b)噴注面板

    圖1 氫氧發(fā)動(dòng)機(jī)幾何模型

    1.2 網(wǎng)格劃分與邊界條件

    經(jīng)對稱簡化后,噴注面板共有120個(gè)噴注單元,每個(gè)噴注單元都為氣液同軸直流式噴嘴。由于液氫是作為發(fā)動(dòng)機(jī)管壁再生冷卻的冷卻劑,經(jīng)過冷卻套后,液氫被加熱轉(zhuǎn)化為氣態(tài),因而將H2作為氣相,O2作為液相進(jìn)行計(jì)算。由于噴注面板噴嘴數(shù)量過多,增加了網(wǎng)格劃分的復(fù)雜程度,因此選用非結(jié)構(gòu)網(wǎng)格對幾何模型進(jìn)行貼合,提高計(jì)算精度,降低網(wǎng)格劃分難度。

    氫氧發(fā)動(dòng)機(jī)短噴管推力室計(jì)算網(wǎng)格如圖2所示。從圖2可知,由發(fā)動(dòng)機(jī)入口可給定內(nèi)流場計(jì)算邊界條件(見表1);發(fā)動(dòng)機(jī)壁面采用標(biāo)準(zhǔn)壁面函數(shù),無滑移邊界條件;對稱面上徑向速度為零,所有變量徑向梯度為零。

    表1 氫氧發(fā)動(dòng)機(jī)內(nèi)流場計(jì)算邊界條件

    氫氧發(fā)動(dòng)機(jī)尾焰外流場仿真采用三維180°對稱計(jì)算網(wǎng)格,如圖3所示。

    a)全流場

    b)局部放大圖

    圖3 氫氧發(fā)動(dòng)機(jī)尾焰流場計(jì)算網(wǎng)格

    仿真模型由噴管擴(kuò)張段內(nèi)流場及尾焰外流場組成。噴管擴(kuò)張段面積比為49∶1,模型以噴管喉部為入口,入口條件通過氫氧發(fā)動(dòng)機(jī)內(nèi)流場計(jì)算結(jié)果給定。噴管壁面采用標(biāo)準(zhǔn)壁面函數(shù)求解,無滑移邊界條件。環(huán)境遠(yuǎn)場邊界及出口邊界給定環(huán)境壓力0.1 MPa、溫度300 K以及空氣組分質(zhì)量分?jǐn)?shù),其中來流速度忽略不計(jì)。對稱面上徑向速度為零,所有變量徑向梯度為零。

    1.3 物理模型

    氫氧發(fā)動(dòng)機(jī)內(nèi)氣相流動(dòng)過程采用帶化學(xué)反應(yīng)的三維N-S方程[11]描述,采用Realizable k-ε湍流模型封閉流動(dòng)方程組,液相流動(dòng)過程采用離散相模型進(jìn)行描述,采用Rosin-Rammler分布[12]確定噴射液滴初始尺寸分布,采用Wave模型[13]和O’Rourke模型[14]對液滴破碎及碰撞過程進(jìn)行計(jì)算,以獲得更為精確的計(jì)算結(jié)果。燃燒過程化學(xué)反應(yīng)速率CH采用湍流脈動(dòng)機(jī)制EBU和Arrhenius機(jī)制Arr控制[15],氫氧發(fā)動(dòng)機(jī)化學(xué)反應(yīng)采用氫氧單步化學(xué)反應(yīng),即2H2+O2=2H2O。流動(dòng)方程對流項(xiàng)采用QUICK格式進(jìn)行離散,采用PISO算法進(jìn)行壓力-速度耦合求解。輻射模型采用氣體輻射傳輸方程[16]和大氣透過率計(jì)算模型SLG[17]進(jìn)行計(jì)算。

    2 計(jì)算結(jié)果分析

    2.1 氫氧發(fā)動(dòng)機(jī)內(nèi)流場計(jì)算

    發(fā)動(dòng)機(jī)內(nèi)流場計(jì)算結(jié)果的各項(xiàng)參數(shù)分布如圖4所示。

    a)壓力

    b)溫度

    c)H2O質(zhì)量分?jǐn)?shù)

    d)H2質(zhì)量分?jǐn)?shù)

    從圖4中可以看到,發(fā)動(dòng)機(jī)內(nèi)每個(gè)截面的壓力、溫度和各個(gè)燃?xì)饨M分質(zhì)量分?jǐn)?shù)的分布數(shù)據(jù)。在壓力云圖中,發(fā)動(dòng)機(jī)在燃燒室中達(dá)到最大壓力為8.78 MPa,并在噴管內(nèi)壓力逐步降低,其中燃燒室最大壓力與該工況下試車測得的燃燒室壓力8.87 MPa極為接近,證明了算法的有效性;發(fā)動(dòng)機(jī)的溫度分布呈現(xiàn)出相同的變化趨勢,只是溫度的變化較小。由圖4還可以看出,在發(fā)動(dòng)機(jī)燃燒室中,H2與O2占主要成分,并會(huì)發(fā)生劇烈的燃燒反應(yīng);在發(fā)動(dòng)機(jī)噴管中,H2O和H2占主要成分,O2因反應(yīng)而消耗殆盡,說明噴管噴出的燃?xì)鉃楦蝗既細(xì)狻?/p>

    2.2 復(fù)燃對尾焰流場特性影響

    為方便對比,將復(fù)燃反應(yīng)流與凍結(jié)流流場參數(shù)進(jìn)行對比分析。圖5給出復(fù)燃反應(yīng)對溫度場的影響。

    a)全流場

    b)軸線上溫度分布

    c)發(fā)動(dòng)機(jī)出口不同位置處縱向溫度分布

    從圖5a中可見,在=0~14 m區(qū)間,反應(yīng)流高溫區(qū)主要存在于尾焰邊界區(qū)域;在=14~26 m區(qū)間,反應(yīng)流高溫區(qū)則存在于尾焰中心區(qū)域。主要是由于在=0~14 m范圍內(nèi),尾焰流速較快,燃?xì)馀c周圍空氣的摻混主要發(fā)生在尾焰邊界,隨即發(fā)生復(fù)燃反應(yīng),導(dǎo)致該區(qū)域溫度升高;隨著流動(dòng)的進(jìn)行,尾焰流速逐漸減慢,燃?xì)馀c周圍空氣的摻混加劇,復(fù)燃反應(yīng)隨之加強(qiáng),造成中心區(qū)域溫度升高。圖5b、5c中,軸線上溫度場的最大相對增幅達(dá)40.48%,發(fā)動(dòng)機(jī)出口=5 m、10 m和20 m位置處縱向溫度最大增幅均達(dá)30%以上。由此可知,復(fù)燃反應(yīng)對氫氧發(fā)動(dòng)機(jī)尾焰流場計(jì)算影響很大。

    圖6給出了計(jì)算結(jié)果同文獻(xiàn)[18]提供的實(shí)驗(yàn)試車與仿真數(shù)據(jù)的對比。從圖6中可以看出,計(jì)算結(jié)果同文獻(xiàn)中的結(jié)果符合較好,證明了算法的有效性和正確性。圖7給出了復(fù)燃對尾焰流場軸線上各組分質(zhì)量分?jǐn)?shù)分布影響。

    a)文獻(xiàn)實(shí)驗(yàn)結(jié)果

    b)文獻(xiàn)仿真結(jié)果

    c)本文仿真結(jié)果

    圖6 氫氧發(fā)動(dòng)機(jī)尾焰流場計(jì)算結(jié)果同文獻(xiàn)[18]結(jié)果對比

    a)H2O

    b)H2

    從圖7中可以看出,在尾焰流場中,氣態(tài)H2O是尾焰的主要成分。氫氧發(fā)動(dòng)機(jī)尾焰紅外輻射計(jì)算應(yīng)圍繞H2O的輻射進(jìn)行。從復(fù)燃的影響來看,復(fù)燃反應(yīng)會(huì)使尾焰中H2O的質(zhì)量分?jǐn)?shù)明顯增加,而H2和O2相應(yīng)減少,軸線上H2O質(zhì)量分?jǐn)?shù)最大相對增幅達(dá)43.75%。

    2.3 復(fù)燃對尾焰輻射特性影響

    將氫氧發(fā)動(dòng)機(jī)尾焰流場計(jì)算得到的壓力、溫度以及燃?xì)饨M分的質(zhì)量分?jǐn)?shù)帶入尾焰輻射計(jì)算模型中進(jìn)行計(jì)算,計(jì)算結(jié)果如圖8所示。對于氫氧發(fā)動(dòng)機(jī)其燃燒產(chǎn)物主要是H2O,也是主要輻射體,其輻射主要體現(xiàn)在2.7 μm波段上。從圖8中可以看出,復(fù)燃反應(yīng)會(huì)使整個(gè)波段內(nèi)的光譜輻射照度大幅增加,在=2.67 μm處增幅最大。這主要是由于復(fù)燃反應(yīng)會(huì)對尾焰溫度場和燃?xì)赓|(zhì)量分?jǐn)?shù)分布產(chǎn)生影響,從而影響尾焰的輻射特性。由此可知,復(fù)燃反應(yīng)對氫氧發(fā)動(dòng)機(jī)尾焰輻射計(jì)算的影響同樣不容忽視(光譜輻射照度表示波長為λ時(shí),單位波長間隔內(nèi)的入射物體單位面積上接收到的輻射功率[19])。

    3 結(jié) 論

    本文通過求解氫氧發(fā)動(dòng)機(jī)尾焰流場與輻射計(jì)算模型,定量分析了復(fù)燃反應(yīng)對發(fā)動(dòng)機(jī)尾焰流場和輻射特性的影響,所得結(jié)論如下:

    a)復(fù)燃反應(yīng)會(huì)改變尾焰溫度場和燃?xì)饨M分分布。對于溫度場,復(fù)燃反應(yīng)會(huì)使尾焰溫度大幅升高,最大增幅可達(dá)40.48%。對于燃?xì)饨M分分布,復(fù)燃反應(yīng)會(huì)使H2O的含量增加,相應(yīng)反應(yīng)物的含量減少,H2O質(zhì)量分?jǐn)?shù)的最大增幅達(dá)43.75%。

    b)氫氧發(fā)動(dòng)機(jī)尾焰主要輻射波段在2.7 μm波段,復(fù)燃反應(yīng)會(huì)使整個(gè)波段內(nèi)的光譜輻射照度大幅增加,在=2.67μm處增幅最大。

    c)復(fù)燃反應(yīng)會(huì)對氫氧發(fā)動(dòng)機(jī)尾焰流場和輻射特性計(jì)算產(chǎn)生較大影響,因而在相應(yīng)計(jì)算中應(yīng)考慮復(fù)燃反應(yīng)以獲得較為精確的結(jié)果。

    [1] 苗瑞生. 發(fā)射氣體動(dòng)力學(xué)[M]. 北京:國防工業(yè)出版社, 2006.

    [2] Leone D M, Turns S R. Active Chlorine and Nitric Oxide Formation From Chemical Rocket Plume Afterburning[R]. AIAA 94-0788, 1994.

    [3] 王雁鳴, 談和平, 董士奎, 等. 低空多噴管發(fā)動(dòng)機(jī)噴焰紅外輻射特性研究[J]. 固體火箭技術(shù), 2009, 32(6): 634-637.

    [4] 聶萬勝, 楊軍輝, 何浩波, 等. 液體火箭發(fā)動(dòng)機(jī)尾噴焰紅外輻射特性[J]. 國防科技大學(xué)學(xué)報(bào), 2005, 27(5):91-94.

    [5] 豐松江, 聶萬勝, 解慶紛, 等. 燃燒室內(nèi)燃燒模型對尾焰流場及其輻射的影響[J]. 火箭推進(jìn), 2006, 32(2):6-10.

    [6] WANG W C, WEI Z J, ZHANG Q, et al. Study on infrared signature of solid rocket motor afterburning exhaust plume[R]. AIAA 2010-6847, 2010.

    [7] 劉尊洋, 邵立, 汪亞夫, 等. 復(fù)燃對固體火箭尾焰紅外輻射特性的影響[J]. 光學(xué)學(xué)報(bào), 2013, 33(6): 1-8.

    [8] 劉尊洋, 邵立, 汪亞夫, 等. 復(fù)燃對液體火箭尾焰紅外輻射特性的影響[J]. 光子學(xué)報(bào), 2013, 42(4): 480-485.

    [9] 劉尊洋, 邵立, 汪亞夫, 等. 飛行參數(shù)對液體火箭尾焰紅外輻射特性的影響[J]. 光學(xué)學(xué)報(bào), 2013, 33(4): 1-7.

    [10] 丁兆波, 孫繼國, 路曉紅. 國外典型大推力氫氧發(fā)動(dòng)機(jī)推力室技術(shù)方案綜述[J]. 導(dǎo)彈與航天運(yùn)載技術(shù), 2012(4): 28-30.

    [11] 陶文銓. 數(shù)值傳熱學(xué)[M]. 西安: 西安交通大學(xué)出版社, 1988.

    [12] 莊逢辰. 液體火箭發(fā)動(dòng)機(jī)噴霧燃燒的理論、模型及應(yīng)用[M]. 長沙: 國防科技大學(xué)出版社, 1995.

    [13] Reitz R D. Mechanisms of atomization processes in high-pressure vaporizing sprays[J]. Atomization and Spray Technology, 1987(3): 309-337.

    [14] O’Rourke P J. Collective drop effects on vaporizing liquid sprays[D]. New Jersey: Princeton University, 1981.

    [15] FENG S J, CHENG Y F, NIE W S. Combustion instability analysis in a Hydrogen-Oxygen rocket engine[R]. AIAA 2009-4864, 2009.

    [16] 楊軍輝. 液體火箭燃燒過程及其尾部噴焰輻射特性研究[D]. 北京: 裝備指揮技術(shù)學(xué)院,2004.

    [17] Ludwig C B, et al. Handbook of infrared radiation from combustion gases[R]. NASA SP-3080, 1973.

    [18] 李茂, 王占林. 氫氧火箭發(fā)動(dòng)機(jī)射流仿真與試驗(yàn)臺(tái)熱防護(hù)[J]. 載人航天, 2014, 20(5):421-426.

    [19] 徐南榮, 卞南華. 紅外輻射與制導(dǎo)[M]. 北京: 國防工業(yè)出版社, 1997.

    Numerical Research on Influence Exerted by Afterburning on Flow Field and Radiation Characteristics of LH2/LOX Rocket Engine

    Qiao Ye, Nie Wan-sheng, Feng Song-jiang, Wu Gao-yang

    (Department of Space Equipment, Equipment Academy of PLA, Beijing, 101416)

    To have a deep understanding of the influence exerted by afterburning on the flow field and radiation characteristics of LH2/LOX rocket engine plume, parameters at the throat section of the engine are set as the inlet condition to consider the afterburning reaction. The N-S equation coupled with Realizable k-ε model is used as well. Plume flow field parameters are obtained through PISO algorithm. On this basis, plume radiation characteristics are calculated with gas radiation transmission equation and SLG model. The influence of afterburning reaction on the plume flow field and radiation characteristics are compared as well. The results indicate that the afterburning reaction exerts a significant effect on plume flow field calculation, making the mass fraction of combustion products go up greatly, which enhances plume’s radiation characteristics. In conclusion, it is necessary to take account afterburning reaction in plume flow field and radiation calculations.

    LH2/LOX rocket engine; Plume; Afterburning; Infrared radiation; Numerical simulation

    1004-7182(2016)02-0022-05

    10.7654/j.issn.1004-7182.20160205

    V43

    A

    2015-08-21;

    2015-12-25

    國家自然科學(xué)基金(51206185,91441123)

    喬野(1991-),男,碩士研究生,研究領(lǐng)域?yàn)橐后w火箭發(fā)動(dòng)機(jī)

    猜你喜歡
    尾焰氫氧流場
    基于粒子系統(tǒng)的尾焰紅外圖像實(shí)時(shí)仿真技術(shù)
    氫氧燃料電池演示實(shí)驗(yàn)的改進(jìn)
    自制液壓儲(chǔ)氣式氫氧燃料電池
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計(jì)算
    固體火箭尾焰等離子體特性影響因素?cái)?shù)值仿真
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    連鑄坯氫氧切割應(yīng)用與碳排放研究
    多噴管液體火箭動(dòng)力系統(tǒng)尾焰輻射特性研究
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計(jì)分析
    多噴管火箭動(dòng)力系統(tǒng)尾焰輻射特性可視化研究
    看十八女毛片水多多多| 一边摸一边做爽爽视频免费| 久久精品国产自在天天线| 亚洲在久久综合| 久久久久国产精品人妻一区二区| 日韩一卡2卡3卡4卡2021年| 岛国毛片在线播放| 国产老妇伦熟女老妇高清| 亚洲av成人精品一二三区| 大香蕉久久网| 国产一级毛片在线| 久久99热这里只频精品6学生| 看非洲黑人一级黄片| 国产白丝娇喘喷水9色精品| 国产黄色视频一区二区在线观看| 久久久久久伊人网av| 欧美日韩综合久久久久久| 99国产精品免费福利视频| 久久影院123| 国产极品天堂在线| 少妇的逼水好多| 国产亚洲一区二区精品| 看免费av毛片| 久久这里只有精品19| 成人免费观看视频高清| 免费黄网站久久成人精品| 亚洲国产最新在线播放| 色吧在线观看| 在线天堂最新版资源| 日本欧美视频一区| 下体分泌物呈黄色| 午夜久久久在线观看| 亚洲一区二区三区欧美精品| 看免费成人av毛片| 乱人伦中国视频| 国语对白做爰xxxⅹ性视频网站| 成人国产av品久久久| 你懂的网址亚洲精品在线观看| 老熟女久久久| 精品国产国语对白av| 亚洲成国产人片在线观看| 国产精品无大码| 黄色配什么色好看| 成年av动漫网址| 日本欧美国产在线视频| 色哟哟·www| 日韩人妻精品一区2区三区| 一区福利在线观看| 日本欧美国产在线视频| 国产一区二区 视频在线| 欧美人与性动交α欧美精品济南到 | 久久久久国产一级毛片高清牌| 午夜福利网站1000一区二区三区| 午夜影院在线不卡| 97人妻天天添夜夜摸| 久久久久人妻精品一区果冻| 日本av免费视频播放| 女性生殖器流出的白浆| 国产日韩欧美亚洲二区| 人妻少妇偷人精品九色| xxxhd国产人妻xxx| 日韩制服骚丝袜av| 纵有疾风起免费观看全集完整版| 亚洲欧美成人精品一区二区| 亚洲国产成人一精品久久久| 夜夜骑夜夜射夜夜干| 国产精品人妻久久久影院| 中文天堂在线官网| 国产淫语在线视频| 晚上一个人看的免费电影| 美女福利国产在线| 欧美精品国产亚洲| 只有这里有精品99| 国产成人一区二区在线| 午夜精品国产一区二区电影| 国产av国产精品国产| 精品久久久久久电影网| 精品国产乱码久久久久久男人| 亚洲精品一区蜜桃| av国产精品久久久久影院| 免费播放大片免费观看视频在线观看| 免费高清在线观看日韩| 欧美精品国产亚洲| 91午夜精品亚洲一区二区三区| 蜜桃在线观看..| 国产精品久久久久久久久免| 亚洲国产欧美网| 久久精品久久久久久噜噜老黄| av.在线天堂| 日韩视频在线欧美| 日韩 亚洲 欧美在线| 中文乱码字字幕精品一区二区三区| 美女福利国产在线| 亚洲精品久久午夜乱码| 夫妻午夜视频| 在线观看国产h片| 在线观看美女被高潮喷水网站| 日韩人妻精品一区2区三区| 免费大片黄手机在线观看| 国产成人午夜福利电影在线观看| 国产日韩欧美亚洲二区| 久久久久久人人人人人| 少妇 在线观看| a 毛片基地| 黄色毛片三级朝国网站| 亚洲 欧美一区二区三区| 欧美中文综合在线视频| 亚洲一级一片aⅴ在线观看| 日本wwww免费看| 大香蕉久久成人网| 国产成人精品无人区| 久久人妻熟女aⅴ| 十八禁网站网址无遮挡| 另类亚洲欧美激情| 男人添女人高潮全过程视频| 丝袜喷水一区| www.自偷自拍.com| 精品少妇一区二区三区视频日本电影 | 亚洲欧美成人综合另类久久久| 丝袜人妻中文字幕| 久热这里只有精品99| 最近最新中文字幕大全免费视频 | 美女脱内裤让男人舔精品视频| 午夜福利乱码中文字幕| 日韩中文字幕欧美一区二区 | 免费黄频网站在线观看国产| 久久国内精品自在自线图片| 久久影院123| 啦啦啦在线观看免费高清www| 香蕉国产在线看| 久久免费观看电影| 国产精品一区二区在线观看99| 精品一品国产午夜福利视频| 久久久久久久国产电影| 三级国产精品片| 国产精品久久久久久精品电影小说| 99久久综合免费| 久久久久久久精品精品| 熟女少妇亚洲综合色aaa.| 亚洲av成人精品一二三区| 亚洲精品视频女| 亚洲人成电影观看| 精品酒店卫生间| 男人添女人高潮全过程视频| 久久久久久久国产电影| 99热全是精品| 午夜福利影视在线免费观看| 国产成人av激情在线播放| 日本猛色少妇xxxxx猛交久久| 中国三级夫妇交换| 五月伊人婷婷丁香| 亚洲av欧美aⅴ国产| 人体艺术视频欧美日本| 国产欧美亚洲国产| 深夜精品福利| 亚洲人成网站在线观看播放| 精品亚洲成国产av| 老熟女久久久| 女的被弄到高潮叫床怎么办| 在线免费观看不下载黄p国产| 最近手机中文字幕大全| 国产av国产精品国产| 精品第一国产精品| 男女国产视频网站| 亚洲欧洲日产国产| 青青草视频在线视频观看| 国产精品 欧美亚洲| 有码 亚洲区| 大片电影免费在线观看免费| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美成人综合另类久久久| 成人手机av| 亚洲人成电影观看| 最新的欧美精品一区二区| 大话2 男鬼变身卡| 久久久久久久久久人人人人人人| 国产不卡av网站在线观看| 80岁老熟妇乱子伦牲交| 美女中出高潮动态图| 欧美成人精品欧美一级黄| 免费黄色在线免费观看| 男人舔女人的私密视频| 中国三级夫妇交换| 日本-黄色视频高清免费观看| 精品人妻熟女毛片av久久网站| 在线观看人妻少妇| 国产成人欧美| 在线 av 中文字幕| 一本大道久久a久久精品| 91精品国产国语对白视频| 日本爱情动作片www.在线观看| 久久久国产一区二区| 大码成人一级视频| 成人国产麻豆网| 精品一区二区免费观看| 美女国产视频在线观看| 亚洲伊人久久精品综合| 国产欧美日韩综合在线一区二区| 一区二区av电影网| 久久久久视频综合| 大片电影免费在线观看免费| 亚洲内射少妇av| 一级毛片电影观看| 丰满乱子伦码专区| 丰满迷人的少妇在线观看| 成人影院久久| 多毛熟女@视频| 精品第一国产精品| 精品少妇久久久久久888优播| 黄片播放在线免费| 午夜免费鲁丝| 热re99久久国产66热| 午夜福利视频精品| 色播在线永久视频| 午夜福利影视在线免费观看| 亚洲国产色片| 老汉色∧v一级毛片| 97在线人人人人妻| 美女中出高潮动态图| 91精品三级在线观看| 国产免费现黄频在线看| 国产精品无大码| 成人国产av品久久久| 久久精品国产综合久久久| 一本大道久久a久久精品| a级毛片黄视频| 亚洲欧美日韩另类电影网站| 精品99又大又爽又粗少妇毛片| 久久精品国产a三级三级三级| 美女午夜性视频免费| 国产伦理片在线播放av一区| 精品酒店卫生间| 国产av码专区亚洲av| 成人毛片60女人毛片免费| 亚洲av中文av极速乱| 最黄视频免费看| 亚洲成国产人片在线观看| 亚洲欧洲日产国产| av国产精品久久久久影院| 咕卡用的链子| 曰老女人黄片| 三级国产精品片| 999久久久国产精品视频| 国产在线一区二区三区精| 亚洲国产毛片av蜜桃av| 18+在线观看网站| 国产高清不卡午夜福利| 免费女性裸体啪啪无遮挡网站| 国产97色在线日韩免费| 国产在视频线精品| 少妇人妻精品综合一区二区| 香蕉精品网在线| 999久久久国产精品视频| 国产亚洲精品第一综合不卡| 成年动漫av网址| 少妇的逼水好多| 日韩熟女老妇一区二区性免费视频| 免费高清在线观看日韩| 亚洲精品aⅴ在线观看| 69精品国产乱码久久久| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 免费黄网站久久成人精品| 久久久精品免费免费高清| 国产一区二区三区综合在线观看| 欧美日韩一区二区视频在线观看视频在线| 美女大奶头黄色视频| 欧美日韩精品成人综合77777| 国产成人免费观看mmmm| 夫妻午夜视频| 国产精品久久久久久av不卡| av女优亚洲男人天堂| 女的被弄到高潮叫床怎么办| 久久久a久久爽久久v久久| av视频免费观看在线观看| 日本黄色日本黄色录像| 亚洲av日韩在线播放| 观看av在线不卡| 午夜久久久在线观看| 免费日韩欧美在线观看| 亚洲精品在线美女| 男人添女人高潮全过程视频| 国产成人精品福利久久| 只有这里有精品99| 人人妻人人添人人爽欧美一区卜| 考比视频在线观看| 午夜影院在线不卡| 亚洲中文av在线| 一级片免费观看大全| 菩萨蛮人人尽说江南好唐韦庄| 色哟哟·www| 国精品久久久久久国模美| 国产精品熟女久久久久浪| 欧美亚洲日本最大视频资源| 人妻 亚洲 视频| 色94色欧美一区二区| 2022亚洲国产成人精品| 色吧在线观看| 校园人妻丝袜中文字幕| 水蜜桃什么品种好| 老熟女久久久| 在线精品无人区一区二区三| 满18在线观看网站| 极品人妻少妇av视频| 久久影院123| 中文乱码字字幕精品一区二区三区| 天天躁日日躁夜夜躁夜夜| 不卡av一区二区三区| 一本大道久久a久久精品| 一级a爱视频在线免费观看| 日韩制服丝袜自拍偷拍| 蜜桃在线观看..| 欧美激情高清一区二区三区 | 美女脱内裤让男人舔精品视频| 超色免费av| 国产亚洲午夜精品一区二区久久| 国产成人a∨麻豆精品| 精品99又大又爽又粗少妇毛片| 久久久国产精品麻豆| 啦啦啦在线观看免费高清www| 午夜福利视频精品| 亚洲欧美日韩另类电影网站| 精品亚洲成国产av| 日本黄色日本黄色录像| 欧美xxⅹ黑人| 人人妻人人澡人人看| 多毛熟女@视频| 精品久久久精品久久久| 中文字幕亚洲精品专区| 考比视频在线观看| 色网站视频免费| 日韩av在线免费看完整版不卡| 99热网站在线观看| 一区二区三区四区激情视频| 久久人人爽人人片av| 亚洲伊人色综图| 国产精品二区激情视频| 精品福利永久在线观看| 桃花免费在线播放| 美女福利国产在线| 欧美日韩亚洲高清精品| 亚洲视频免费观看视频| 制服丝袜香蕉在线| 女人精品久久久久毛片| 91在线精品国自产拍蜜月| 青春草国产在线视频| 午夜免费观看性视频| 成年美女黄网站色视频大全免费| 久久精品国产自在天天线| 亚洲国产精品一区三区| 久久人人爽人人片av| 日韩电影二区| 日本免费在线观看一区| 最黄视频免费看| 亚洲av男天堂| 亚洲av成人精品一二三区| 天堂中文最新版在线下载| 欧美 日韩 精品 国产| 午夜激情久久久久久久| 少妇的丰满在线观看| 国产一区二区在线观看av| 爱豆传媒免费全集在线观看| 欧美精品av麻豆av| 美女xxoo啪啪120秒动态图| 欧美精品亚洲一区二区| 国产又爽黄色视频| 久久99蜜桃精品久久| 人妻人人澡人人爽人人| 欧美人与善性xxx| xxx大片免费视频| 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 黑丝袜美女国产一区| 亚洲国产av新网站| 高清不卡的av网站| 熟妇人妻不卡中文字幕| 五月开心婷婷网| 国产熟女午夜一区二区三区| 久久ye,这里只有精品| 国产精品99久久99久久久不卡 | 晚上一个人看的免费电影| 高清视频免费观看一区二区| 亚洲精品久久久久久婷婷小说| 高清av免费在线| 美女国产高潮福利片在线看| 欧美97在线视频| 永久免费av网站大全| 各种免费的搞黄视频| 五月天丁香电影| 国产又爽黄色视频| 大话2 男鬼变身卡| 日本猛色少妇xxxxx猛交久久| 九色亚洲精品在线播放| 欧美97在线视频| 交换朋友夫妻互换小说| 亚洲综合色网址| 人妻人人澡人人爽人人| 极品少妇高潮喷水抽搐| 香蕉丝袜av| 一本色道久久久久久精品综合| 欧美精品亚洲一区二区| 美女中出高潮动态图| 久久狼人影院| 永久网站在线| 国产成人免费观看mmmm| 我要看黄色一级片免费的| 日韩电影二区| 亚洲一区二区三区欧美精品| 久久久久久久久免费视频了| 成年av动漫网址| 美女脱内裤让男人舔精品视频| xxx大片免费视频| 黄色视频在线播放观看不卡| 18禁国产床啪视频网站| 91aial.com中文字幕在线观看| 日韩成人av中文字幕在线观看| 欧美日韩亚洲高清精品| 女人久久www免费人成看片| 看非洲黑人一级黄片| 亚洲欧美一区二区三区久久| videos熟女内射| 欧美精品一区二区大全| 久久午夜综合久久蜜桃| 飞空精品影院首页| 伦理电影大哥的女人| 日韩大片免费观看网站| 成年人午夜在线观看视频| 秋霞伦理黄片| 男女边摸边吃奶| av在线播放精品| 久久精品国产鲁丝片午夜精品| 亚洲欧美精品自产自拍| 91久久精品国产一区二区三区| 在线观看www视频免费| 女性生殖器流出的白浆| 少妇的丰满在线观看| 欧美精品亚洲一区二区| 久久精品国产亚洲av涩爱| 国产一区二区在线观看av| 国产精品国产三级国产专区5o| www.av在线官网国产| av免费在线看不卡| 一级a爱视频在线免费观看| 国产欧美亚洲国产| 亚洲av.av天堂| 久久人人爽人人片av| 精品酒店卫生间| 中文精品一卡2卡3卡4更新| 国产亚洲一区二区精品| 国产免费现黄频在线看| www.精华液| 在线亚洲精品国产二区图片欧美| 香蕉丝袜av| 黑人欧美特级aaaaaa片| 亚洲精品一区蜜桃| 日韩人妻精品一区2区三区| 人人妻人人添人人爽欧美一区卜| 老女人水多毛片| 在线观看免费日韩欧美大片| 天天躁日日躁夜夜躁夜夜| 赤兔流量卡办理| 亚洲欧美成人综合另类久久久| 亚洲欧美中文字幕日韩二区| 亚洲国产日韩一区二区| 各种免费的搞黄视频| 亚洲精品在线美女| 久久人妻熟女aⅴ| 久久综合国产亚洲精品| 国产亚洲精品第一综合不卡| freevideosex欧美| 午夜91福利影院| 黄片小视频在线播放| 最近最新中文字幕免费大全7| 国产精品久久久久久精品电影小说| 亚洲在久久综合| 亚洲精品美女久久久久99蜜臀 | 午夜av观看不卡| 丰满少妇做爰视频| a级片在线免费高清观看视频| 亚洲欧美成人综合另类久久久| 免费黄色在线免费观看| 亚洲av电影在线进入| 男女高潮啪啪啪动态图| 精品99又大又爽又粗少妇毛片| 成人毛片a级毛片在线播放| 啦啦啦在线免费观看视频4| 波多野结衣av一区二区av| av免费观看日本| 国产老妇伦熟女老妇高清| 久久久久国产网址| www.熟女人妻精品国产| 国产成人欧美| 最近最新中文字幕大全免费视频 | 欧美老熟妇乱子伦牲交| 最近最新中文字幕免费大全7| 999精品在线视频| 亚洲欧美一区二区三区黑人 | 少妇的逼水好多| 久久久久精品久久久久真实原创| 欧美日韩av久久| 99精国产麻豆久久婷婷| 国产人伦9x9x在线观看 | 大香蕉久久成人网| 中文字幕精品免费在线观看视频| 捣出白浆h1v1| 国产精品一区二区在线不卡| videos熟女内射| 欧美精品一区二区免费开放| 成年人午夜在线观看视频| 一级毛片黄色毛片免费观看视频| 观看美女的网站| 美女中出高潮动态图| 免费播放大片免费观看视频在线观看| 久久人人爽人人片av| 亚洲第一区二区三区不卡| 黄色毛片三级朝国网站| 国产伦理片在线播放av一区| 久久久国产一区二区| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品久久久久久婷婷小说| 午夜老司机福利剧场| 亚洲精品美女久久av网站| 国产精品免费大片| 三上悠亚av全集在线观看| av网站免费在线观看视频| 亚洲国产最新在线播放| 日韩,欧美,国产一区二区三区| 精品第一国产精品| 久久久国产欧美日韩av| 三上悠亚av全集在线观看| av在线观看视频网站免费| 美国免费a级毛片| 免费大片黄手机在线观看| 中文字幕最新亚洲高清| 亚洲国产看品久久| 极品人妻少妇av视频| 精品国产超薄肉色丝袜足j| 母亲3免费完整高清在线观看 | 下体分泌物呈黄色| 你懂的网址亚洲精品在线观看| 亚洲,一卡二卡三卡| 日韩在线高清观看一区二区三区| 亚洲成国产人片在线观看| 韩国精品一区二区三区| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 边亲边吃奶的免费视频| 午夜日韩欧美国产| av有码第一页| 久久青草综合色| 亚洲,欧美,日韩| 伦理电影免费视频| 午夜免费观看性视频| 免费不卡的大黄色大毛片视频在线观看| 日韩制服丝袜自拍偷拍| 欧美bdsm另类| 黄色 视频免费看| 欧美精品一区二区大全| 男人爽女人下面视频在线观看| 精品国产露脸久久av麻豆| av卡一久久| 99re6热这里在线精品视频| 中国国产av一级| 午夜福利乱码中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩一区二区视频在线观看视频在线| 看免费av毛片| 欧美变态另类bdsm刘玥| 有码 亚洲区| 国产淫语在线视频| 久久久久国产一级毛片高清牌| 精品亚洲成国产av| 一级毛片 在线播放| 国产亚洲av片在线观看秒播厂| av片东京热男人的天堂| 精品亚洲乱码少妇综合久久| 久久久精品国产亚洲av高清涩受| 国产爽快片一区二区三区| 中文字幕人妻熟女乱码| 国产一级毛片在线| 国产 一区精品| 精品国产超薄肉色丝袜足j| 在线观看免费日韩欧美大片| 亚洲欧美一区二区三区国产| 精品人妻偷拍中文字幕| 丝袜脚勾引网站| 亚洲av欧美aⅴ国产| 国产精品99久久99久久久不卡 | 叶爱在线成人免费视频播放| av国产久精品久网站免费入址| 国产一区二区激情短视频 | 亚洲色图综合在线观看| 免费大片黄手机在线观看| 免费观看性生交大片5| 免费在线观看完整版高清| 美女xxoo啪啪120秒动态图| 伦精品一区二区三区| 午夜免费鲁丝| 不卡视频在线观看欧美| 国精品久久久久久国模美| 伦理电影大哥的女人| 亚洲国产欧美网| 在线天堂最新版资源| 国产精品熟女久久久久浪| 久久99一区二区三区| 少妇 在线观看| av在线观看视频网站免费| 久久久久久免费高清国产稀缺| 哪个播放器可以免费观看大片| 国产精品蜜桃在线观看| 日本wwww免费看| 老女人水多毛片| 在线观看三级黄色| 精品人妻在线不人妻| 中国国产av一级|