• <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)尾焰輻射特性可視化研究
    亚洲国产欧美在线一区| 国产乱来视频区| 亚洲国产欧美在线一区| 日本-黄色视频高清免费观看| 99久久人妻综合| a级毛片免费高清观看在线播放| 久久人人爽人人片av| 久久综合国产亚洲精品| 中文精品一卡2卡3卡4更新| 成人一区二区视频在线观看| 久久 成人 亚洲| 精品久久久精品久久久| 国产精品无大码| 国产色爽女视频免费观看| 国产久久久一区二区三区| 99久久精品国产国产毛片| 老司机影院成人| 在线观看免费高清a一片| 人人妻人人添人人爽欧美一区卜 | 一区二区av电影网| 九九久久精品国产亚洲av麻豆| 一区二区av电影网| 天堂8中文在线网| 亚洲人与动物交配视频| 一级爰片在线观看| 人妻夜夜爽99麻豆av| 久久久久久久久大av| 亚洲欧美日韩卡通动漫| 狠狠精品人妻久久久久久综合| 熟女av电影| 久久久亚洲精品成人影院| 国产精品人妻久久久影院| 六月丁香七月| av在线app专区| 国产精品一区二区三区四区免费观看| av国产免费在线观看| 一级二级三级毛片免费看| 在线免费观看不下载黄p国产| 国产日韩欧美在线精品| 亚洲久久久国产精品| 大香蕉久久网| 午夜福利在线在线| 黄色怎么调成土黄色| 国产乱来视频区| 久久久欧美国产精品| 国产av码专区亚洲av| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人久久小说| 一级二级三级毛片免费看| 在线观看美女被高潮喷水网站| 有码 亚洲区| 777米奇影视久久| 最近最新中文字幕免费大全7| 最新中文字幕久久久久| 最新中文字幕久久久久| 国内揄拍国产精品人妻在线| 亚洲综合色惰| 日韩精品有码人妻一区| 爱豆传媒免费全集在线观看| 日韩强制内射视频| 视频中文字幕在线观看| 爱豆传媒免费全集在线观看| 草草在线视频免费看| 在线天堂最新版资源| 中国美白少妇内射xxxbb| 国产在线男女| a级一级毛片免费在线观看| 国产欧美亚洲国产| 黑丝袜美女国产一区| 91久久精品电影网| 亚洲av中文av极速乱| 久久精品人妻少妇| 22中文网久久字幕| 亚洲av欧美aⅴ国产| 久久精品人妻少妇| 午夜福利在线在线| 美女主播在线视频| 深夜a级毛片| 男女下面进入的视频免费午夜| 十分钟在线观看高清视频www | 成人亚洲欧美一区二区av| 国产女主播在线喷水免费视频网站| 一级毛片黄色毛片免费观看视频| 欧美高清成人免费视频www| 成人一区二区视频在线观看| 一级毛片aaaaaa免费看小| tube8黄色片| 久久国产精品大桥未久av | 色5月婷婷丁香| 身体一侧抽搐| 九九久久精品国产亚洲av麻豆| 亚洲色图综合在线观看| 精品久久国产蜜桃| 欧美日韩视频高清一区二区三区二| 亚洲精品日本国产第一区| 熟妇人妻不卡中文字幕| 欧美激情国产日韩精品一区| av.在线天堂| 人体艺术视频欧美日本| 亚洲真实伦在线观看| 一本色道久久久久久精品综合| 久久久久久伊人网av| 国产精品偷伦视频观看了| 少妇丰满av| 亚洲av.av天堂| 天堂8中文在线网| 性色avwww在线观看| 国产精品久久久久成人av| 狂野欧美白嫩少妇大欣赏| 国产欧美日韩一区二区三区在线 | 国产精品女同一区二区软件| 99久久综合免费| 精品人妻熟女av久视频| 熟妇人妻不卡中文字幕| 99久久精品一区二区三区| 国产精品一区二区三区四区免费观看| 男女免费视频国产| 免费高清在线观看视频在线观看| 五月玫瑰六月丁香| 欧美xxxx黑人xx丫x性爽| 国产在线一区二区三区精| 伦理电影大哥的女人| 日韩精品有码人妻一区| 好男人视频免费观看在线| 极品少妇高潮喷水抽搐| 精品人妻偷拍中文字幕| 亚洲欧洲日产国产| 99热这里只有是精品在线观看| 亚洲丝袜综合中文字幕| 免费高清在线观看视频在线观看| 亚洲,一卡二卡三卡| 精品99又大又爽又粗少妇毛片| 性高湖久久久久久久久免费观看| 日本wwww免费看| 一级毛片我不卡| 精品久久久噜噜| 中文字幕av成人在线电影| 久久久精品94久久精品| 免费黄网站久久成人精品| 99精国产麻豆久久婷婷| 亚洲欧洲日产国产| 欧美 日韩 精品 国产| 人人妻人人澡人人爽人人夜夜| 色吧在线观看| 性色av一级| 女性生殖器流出的白浆| 99热网站在线观看| 男人爽女人下面视频在线观看| a级一级毛片免费在线观看| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 国产日韩欧美在线精品| 午夜福利在线在线| 亚洲婷婷狠狠爱综合网| 日韩在线高清观看一区二区三区| 最新中文字幕久久久久| 毛片女人毛片| 人妻一区二区av| 少妇丰满av| 夜夜骑夜夜射夜夜干| 久久精品熟女亚洲av麻豆精品| 亚洲精品亚洲一区二区| 91精品国产国语对白视频| 精品99又大又爽又粗少妇毛片| 在线天堂最新版资源| 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 久久久久国产网址| 亚洲内射少妇av| 建设人人有责人人尽责人人享有的 | 久久精品国产亚洲av涩爱| 最新中文字幕久久久久| 日韩 亚洲 欧美在线| 欧美+日韩+精品| 少妇人妻久久综合中文| 日本黄色片子视频| 永久网站在线| 国精品久久久久久国模美| 国产色爽女视频免费观看| 亚洲av国产av综合av卡| 亚洲va在线va天堂va国产| 国产伦精品一区二区三区四那| 欧美日韩视频精品一区| 毛片女人毛片| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 国产精品蜜桃在线观看| 天堂中文最新版在线下载| 狂野欧美激情性xxxx在线观看| 欧美 日韩 精品 国产| 一级毛片黄色毛片免费观看视频| 一级二级三级毛片免费看| 国产在线一区二区三区精| 欧美精品人与动牲交sv欧美| 国产精品精品国产色婷婷| 国产精品蜜桃在线观看| 边亲边吃奶的免费视频| 国产精品精品国产色婷婷| av视频免费观看在线观看| 极品少妇高潮喷水抽搐| 国产欧美日韩一区二区三区在线 | 成人国产av品久久久| 久久99热6这里只有精品| 十八禁网站网址无遮挡 | 国产精品秋霞免费鲁丝片| 国产精品国产三级国产av玫瑰| 天堂8中文在线网| 亚洲欧美日韩另类电影网站 | 日韩欧美 国产精品| 日韩三级伦理在线观看| 亚洲人成网站高清观看| 成人免费观看视频高清| 2022亚洲国产成人精品| 多毛熟女@视频| 国产片特级美女逼逼视频| 国产日韩欧美亚洲二区| 亚洲在久久综合| 国语对白做爰xxxⅹ性视频网站| 欧美丝袜亚洲另类| 国产白丝娇喘喷水9色精品| 高清av免费在线| 婷婷色av中文字幕| 国产免费一级a男人的天堂| 久久综合国产亚洲精品| 在线观看人妻少妇| 国产精品无大码| 亚洲综合精品二区| 伊人久久精品亚洲午夜| 久久久久视频综合| 狠狠精品人妻久久久久久综合| 天美传媒精品一区二区| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕| 国产精品久久久久成人av| 久久久久久久久久成人| 免费久久久久久久精品成人欧美视频 | 免费黄频网站在线观看国产| 91aial.com中文字幕在线观看| av国产久精品久网站免费入址| 夜夜爽夜夜爽视频| 少妇人妻久久综合中文| 成人二区视频| 欧美精品一区二区大全| 春色校园在线视频观看| 国产精品.久久久| 多毛熟女@视频| 直男gayav资源| 日日啪夜夜撸| 男人爽女人下面视频在线观看| 免费大片黄手机在线观看| 亚州av有码| 晚上一个人看的免费电影| 午夜福利网站1000一区二区三区| 美女福利国产在线 | 高清日韩中文字幕在线| 99热这里只有是精品在线观看| 免费观看在线日韩| 18禁在线无遮挡免费观看视频| 国产午夜精品一二区理论片| 亚州av有码| 特大巨黑吊av在线直播| 在线观看免费日韩欧美大片 | 亚洲aⅴ乱码一区二区在线播放| 国产精品麻豆人妻色哟哟久久| av网站免费在线观看视频| 国产精品秋霞免费鲁丝片| 成人亚洲精品一区在线观看 | 少妇人妻 视频| 精品国产三级普通话版| 日韩av免费高清视频| 噜噜噜噜噜久久久久久91| 性色av一级| 色视频www国产| 99热这里只有精品一区| 亚洲精品国产av蜜桃| 高清毛片免费看| 精品亚洲成a人片在线观看 | 久久久久人妻精品一区果冻| 久久久久久久精品精品| 能在线免费看毛片的网站| 午夜精品国产一区二区电影| 国产伦精品一区二区三区四那| 色网站视频免费| 搡女人真爽免费视频火全软件| 国产av国产精品国产| 99热这里只有精品一区| 亚洲经典国产精华液单| 又大又黄又爽视频免费| 日韩欧美精品免费久久| 精品一品国产午夜福利视频| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区| 免费黄色在线免费观看| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 成人漫画全彩无遮挡| 极品教师在线视频| 国产白丝娇喘喷水9色精品| 在线观看一区二区三区| 男男h啪啪无遮挡| av在线app专区| 国产免费福利视频在线观看| 麻豆国产97在线/欧美| 高清不卡的av网站| 美女视频免费永久观看网站| 午夜精品国产一区二区电影| 亚洲精华国产精华液的使用体验| kizo精华| 老司机影院毛片| 日韩人妻高清精品专区| 亚洲精品自拍成人| 国产乱人视频| 精华霜和精华液先用哪个| 各种免费的搞黄视频| 男男h啪啪无遮挡| 亚洲图色成人| 亚洲成人av在线免费| av在线观看视频网站免费| 国产亚洲91精品色在线| 国产亚洲91精品色在线| 性色avwww在线观看| 国产大屁股一区二区在线视频| 男女免费视频国产| 欧美另类一区| 一级黄片播放器| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 啦啦啦啦在线视频资源| 欧美日本视频| 欧美日韩精品成人综合77777| 99热6这里只有精品| 伦理电影大哥的女人| 深爱激情五月婷婷| 国产av一区二区精品久久 | 如何舔出高潮| av播播在线观看一区| 少妇 在线观看| 夜夜看夜夜爽夜夜摸| 天美传媒精品一区二区| 国产精品不卡视频一区二区| 亚洲国产精品一区三区| 免费看不卡的av| 国产亚洲91精品色在线| 美女中出高潮动态图| av免费在线看不卡| 又爽又黄a免费视频| 亚洲av成人精品一二三区| 成年免费大片在线观看| freevideosex欧美| 99热这里只有是精品在线观看| 久久热精品热| 下体分泌物呈黄色| 纵有疾风起免费观看全集完整版| 色婷婷久久久亚洲欧美| 交换朋友夫妻互换小说| 99久国产av精品国产电影| 亚洲成色77777| 国精品久久久久久国模美| 91狼人影院| 夫妻午夜视频| 国产精品人妻久久久久久| 国产亚洲91精品色在线| 国产日韩欧美在线精品| 高清av免费在线| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 日日啪夜夜爽| 日韩精品有码人妻一区| 国产精品久久久久久久久免| av国产免费在线观看| 99热这里只有是精品在线观看| 91久久精品电影网| 少妇人妻精品综合一区二区| 男女免费视频国产| 日韩亚洲欧美综合| 国产免费一级a男人的天堂| 久久国内精品自在自线图片| 日本-黄色视频高清免费观看| 插阴视频在线观看视频| 免费av不卡在线播放| 精品少妇久久久久久888优播| 午夜免费男女啪啪视频观看| a级毛片免费高清观看在线播放| 国产色婷婷99| 一级毛片电影观看| 另类亚洲欧美激情| 中国国产av一级| 尾随美女入室| 亚洲国产精品一区三区| 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图 | 2021少妇久久久久久久久久久| 国产av码专区亚洲av| 18禁在线无遮挡免费观看视频| 国产白丝娇喘喷水9色精品| 久久韩国三级中文字幕| 成人亚洲精品一区在线观看 | 色5月婷婷丁香| 青春草亚洲视频在线观看| 午夜免费男女啪啪视频观看| 91久久精品国产一区二区三区| 国产成人午夜福利电影在线观看| 天堂8中文在线网| 欧美+日韩+精品| 日本av手机在线免费观看| 亚洲久久久国产精品| 免费观看性生交大片5| 一个人看视频在线观看www免费| 久久青草综合色| 久久av网站| 欧美日韩视频高清一区二区三区二| 91狼人影院| 久久人人爽av亚洲精品天堂 | 国产免费视频播放在线视频| 免费av不卡在线播放| 国模一区二区三区四区视频| 制服丝袜香蕉在线| 熟女人妻精品中文字幕| 这个男人来自地球电影免费观看 | 国语对白做爰xxxⅹ性视频网站| 国产高清不卡午夜福利| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站| 男人爽女人下面视频在线观看| 亚洲精品aⅴ在线观看| 狂野欧美白嫩少妇大欣赏| 国产黄色免费在线视频| 欧美丝袜亚洲另类| 黄片wwwwww| 特大巨黑吊av在线直播| 亚洲av.av天堂| 观看免费一级毛片| 日韩大片免费观看网站| h视频一区二区三区| 777米奇影视久久| 成年女人在线观看亚洲视频| 偷拍熟女少妇极品色| 在线播放无遮挡| 免费大片18禁| 波野结衣二区三区在线| 欧美成人一区二区免费高清观看| 国产女主播在线喷水免费视频网站| 人妻少妇偷人精品九色| 亚洲精品乱码久久久久久按摩| 99热这里只有是精品在线观看| 亚洲欧美日韩无卡精品| 晚上一个人看的免费电影| 亚洲国产av新网站| 啦啦啦啦在线视频资源| av网站免费在线观看视频| 午夜老司机福利剧场| 晚上一个人看的免费电影| 我要看黄色一级片免费的| 99re6热这里在线精品视频| 国产精品欧美亚洲77777| 尤物成人国产欧美一区二区三区| 亚洲av.av天堂| 亚洲精品国产av成人精品| www.av在线官网国产| 国产视频首页在线观看| videos熟女内射| 久久精品国产自在天天线| 少妇 在线观看| 女性生殖器流出的白浆| 两个人的视频大全免费| 久久97久久精品| 黄色配什么色好看| 国国产精品蜜臀av免费| 色婷婷久久久亚洲欧美| 夜夜看夜夜爽夜夜摸| 午夜老司机福利剧场| 一个人看的www免费观看视频| 人人妻人人添人人爽欧美一区卜 | 久久久精品免费免费高清| 一级a做视频免费观看| 免费黄频网站在线观看国产| 国产一区二区三区av在线| 亚洲精品第二区| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看| 久久国产精品大桥未久av | 日本欧美国产在线视频| 黑丝袜美女国产一区| 精品少妇久久久久久888优播| 18+在线观看网站| 有码 亚洲区| 夜夜爽夜夜爽视频| 最近的中文字幕免费完整| 日本-黄色视频高清免费观看| 王馨瑶露胸无遮挡在线观看| 女人久久www免费人成看片| 综合色丁香网| 18+在线观看网站| 久久ye,这里只有精品| 亚洲三级黄色毛片| 国产一区亚洲一区在线观看| 最近最新中文字幕大全电影3| 国产精品av视频在线免费观看| 国产精品蜜桃在线观看| 日本午夜av视频| 久久99热6这里只有精品| 日韩av在线免费看完整版不卡| 国产老妇伦熟女老妇高清| 九九在线视频观看精品| av卡一久久| 国产男女超爽视频在线观看| 在线观看免费高清a一片| 国产毛片在线视频| 国产精品女同一区二区软件| 少妇的逼水好多| 各种免费的搞黄视频| 伊人久久国产一区二区| 久久国产乱子免费精品| 国产69精品久久久久777片| 亚洲一区二区三区欧美精品| 国产 精品1| 成人毛片a级毛片在线播放| 国产精品久久久久成人av| 九九在线视频观看精品| av国产久精品久网站免费入址| 国产白丝娇喘喷水9色精品| 亚洲aⅴ乱码一区二区在线播放| 伊人久久精品亚洲午夜| 18禁在线播放成人免费| kizo精华| 九草在线视频观看| 亚洲,一卡二卡三卡| 黄色日韩在线| 水蜜桃什么品种好| 黄片wwwwww| 精品视频人人做人人爽| 嫩草影院新地址| 身体一侧抽搐| 国产精品伦人一区二区| xxx大片免费视频| 91狼人影院| 春色校园在线视频观看| 黄片无遮挡物在线观看| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| h日本视频在线播放| 亚洲人与动物交配视频| 日韩不卡一区二区三区视频在线| 人妻 亚洲 视频| 久久国产亚洲av麻豆专区| 午夜福利高清视频| 亚洲人成网站在线播| 国产亚洲午夜精品一区二区久久| 性高湖久久久久久久久免费观看| 中文乱码字字幕精品一区二区三区| 深夜a级毛片| 中文在线观看免费www的网站| 精品一区在线观看国产| 亚洲经典国产精华液单| 99热6这里只有精品| 国产精品秋霞免费鲁丝片| 亚洲人成网站在线观看播放| 亚洲aⅴ乱码一区二区在线播放| 我的老师免费观看完整版| 久热久热在线精品观看| 国产欧美另类精品又又久久亚洲欧美| 人人妻人人爽人人添夜夜欢视频 | 国产精品一区www在线观看| 日本午夜av视频| 精品一区二区三卡| 在线观看免费高清a一片| 欧美精品人与动牲交sv欧美| 精品99又大又爽又粗少妇毛片| 亚洲美女视频黄频| 丝袜喷水一区| 日韩av免费高清视频| 中文字幕久久专区| 亚洲内射少妇av| av专区在线播放| 中国国产av一级| 久久精品国产亚洲av涩爱| 久久6这里有精品| 99精国产麻豆久久婷婷| av国产久精品久网站免费入址| 男女免费视频国产| 97精品久久久久久久久久精品| 日本色播在线视频| 久久久久久人妻| 久久精品久久久久久噜噜老黄| 国产精品久久久久久久久免| 日本色播在线视频| 亚洲人成网站在线播| 99久久精品热视频| 国产色爽女视频免费观看| 又粗又硬又长又爽又黄的视频| 久久久精品免费免费高清| 青春草亚洲视频在线观看| 99re6热这里在线精品视频| 亚洲成人手机| 天堂中文最新版在线下载| 亚洲欧美中文字幕日韩二区| 在线观看免费日韩欧美大片 | 国产欧美另类精品又又久久亚洲欧美| 亚洲丝袜综合中文字幕| 搡老乐熟女国产| 黄色一级大片看看| 亚洲高清免费不卡视频| 蜜臀久久99精品久久宅男| 日本欧美视频一区| 欧美zozozo另类| 91午夜精品亚洲一区二区三区| 亚洲图色成人| 少妇裸体淫交视频免费看高清| 天堂中文最新版在线下载| 99热这里只有精品一区| 丰满迷人的少妇在线观看| 亚洲欧美精品自产自拍| 精品99又大又爽又粗少妇毛片| 亚洲精品国产色婷婷电影| 成年美女黄网站色视频大全免费 |