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

    液體火箭發(fā)動機針栓噴注器霧化燃燒技術(shù)研究進展

    2020-08-18 06:40:10朱國強李進賢
    火箭推進 2020年4期
    關鍵詞:發(fā)動機研究

    石 璞,朱國強,李進賢,侯 曉

    (1.西北工業(yè)大學 燃燒、熱結(jié)構(gòu)與內(nèi)流場重點實驗室,陜西 西安710072;2.中國航天科技集團有限公司,北京 100048)

    0 引言

    變推力液體火箭發(fā)動機具有廣闊的應用范圍和良好的應用前景,目前主要應用于衛(wèi)星入軌與降落、空間交會、軌道機動、空間定位和穩(wěn)定、運載火箭回收與重復利用等方面[1-2]。此外,在發(fā)動機推力控制、彈道防御導彈軌道控制和組合發(fā)動機轉(zhuǎn)級等方面也有一定的需求[3]。目前實現(xiàn)發(fā)動機推力調(diào)節(jié)的技術(shù)途徑主要包括[4]:調(diào)節(jié)固定噴注器壓降(如控制閥門等);在噴注器上游管路注入惰性氣體調(diào)節(jié)推進劑流量;通過多路歧管的開閉調(diào)節(jié)有效噴注面積;通過可移動部件調(diào)節(jié)噴注面積;調(diào)節(jié)噴管喉部面積;采用多個獨立的燃燒室;調(diào)節(jié)脈沖參數(shù)(脈寬、重復頻率等)。針栓噴注器是一種通過可移動部件調(diào)節(jié)噴注面積的推力調(diào)節(jié)裝置,具有結(jié)構(gòu)簡單、可靠性高、操作安全[5]、燃燒穩(wěn)定[6]、可實現(xiàn)深度節(jié)流和面關機、可擴展性好[7]等優(yōu)良特性,可實現(xiàn)發(fā)動機大范圍(4∶1以上)推力調(diào)節(jié)。美國TRW(Thompson-Ramo-Wooldridge Inc)的MIRA5000發(fā)動機,在保證穩(wěn)定工作的前提下最大可實現(xiàn)35∶1的推力變比。采用針栓噴注器的發(fā)動機飛行成功率較高。針栓式發(fā)動機的推力適用范圍非常寬(達130 000∶1),最大推力的針栓式發(fā)動機為TRW的TR106發(fā)動機,采用LH2/LOX推進劑,推力為2 900 kN;最小推力的針栓式發(fā)動機為TRW用于“光亮鵝卵石計劃”的一款發(fā)動機,采用N2O4/肼推進劑,推力為22 N。該噴注器的推進劑適用性較好,迄今為止已成功試驗過20余種不同的推進劑組合,其中包括煤粉和空氣組合。此外,該噴注器還有燃燒穩(wěn)定的先天優(yōu)勢,在迄今所有的針栓式發(fā)動機地面試驗及飛行試驗中,很少發(fā)生實質(zhì)性的燃燒不穩(wěn)定[8]。由此可見,針栓噴注器及以其為基礎的推力調(diào)節(jié)技術(shù)具有極大的應用潛力。本文在總結(jié)回顧針栓推力調(diào)節(jié)技術(shù)發(fā)展歷程和應用實例的基礎上,重點從針噴注器霧化特性及針栓式發(fā)動機燃燒流動問題兩個方面進行了總結(jié)分析,提出了針栓噴注器及相應的發(fā)動機技術(shù)研究中需要重點關注的問題。

    1 針栓噴注器技術(shù)的發(fā)展

    1.1 針栓噴注器的起源與典型結(jié)構(gòu)

    針栓噴注器起源于19世紀50年代美國JPL(Jet Propulsion Laboratory)用于研究推進劑反應特性的裝置,此裝置就是針栓噴注器的雛形,如圖1所示。1960年JPL相關人員轉(zhuǎn)至TRW繼續(xù)進行針栓噴注器研制工作,隨著研究的不斷深入,針栓噴注器的構(gòu)型不斷改進,實現(xiàn)了推力可調(diào),可快速脈沖工作及可實現(xiàn)面關機等功能[8]。針栓噴注器的典型結(jié)構(gòu)[9]如圖2所示,其中心徑向流道(綠色部分)是一個帶針栓帽的管道或頂端開孔的中空針栓,環(huán)縫型軸向流道(紫色部分)由調(diào)節(jié)套筒與底座形成,燃料和氧化劑分別從兩個流道流出并發(fā)生碰撞霧化,通過改變調(diào)節(jié)套筒的位置可以在保證所需混合比的情況下實現(xiàn)推進劑流量的大范圍同步調(diào)節(jié),進而調(diào)節(jié)發(fā)動機推力。針栓噴注器的工作狀態(tài)如圖3所示[7],圖3(a)為推進劑僅從軸向流道流出的情形,圖3(b)為推進劑僅從徑向流道流出的情形,圖3(c)為徑向流與軸向流的合成流。

    圖1 JPL用于研究推進劑反應特性的裝置

    圖2 針栓噴注器典型結(jié)構(gòu)

    圖3 針栓噴注器工作狀態(tài)

    1.2 針栓式發(fā)動機的發(fā)展與應用

    美國在針栓式發(fā)動機研究領域處于領先地位,有很多典型的應用。國內(nèi)在該領域的研究起步較晚,目前唯一的實際應用是嫦娥三號的7 500 N下降級發(fā)動機。此外,開展相關研究的國家還包括日本、韓國、印度、德國等。

    1.2.1 國外的發(fā)展與應用

    美國在針栓式發(fā)動機研究領域起步較早。1961年TRW首次將針栓噴注器應用于MIRA500發(fā)動機,該發(fā)動機的推力可在111.2~222 4 N范圍內(nèi)調(diào)節(jié)[10]。1963年TRW開始為JPL研制作為“勘探者發(fā)動機”備份方案的MIRA150A發(fā)動機,推力為133~667 N[1]。1963到1972年間TRW負責研制了LMDE(Lunar Excursion Module Descent Engine)[6,11],該發(fā)動機也使用了針栓噴注器,最大推力為44.52 kN,具有10∶1的推力調(diào)節(jié)能力[4]。LMDE在10次任務中表現(xiàn)出了優(yōu)異的性能,由其衍生而來的固定推力發(fā)動機TR201應用于德爾塔運載火箭的二級發(fā)動機,飛行成功率為100%[8]。從1962年開始,TRW還進行了大量研究,試圖將針栓噴注器應用于大型助推發(fā)動機[12],1995年美國測試了迄今為止世界上最大的針栓式發(fā)動機TR106,該發(fā)動機推力達到了290 t級[3,13],此項研究工作直到現(xiàn)在還在繼續(xù)。

    在執(zhí)行LMDE計劃的同時,TRW還在研制小推力針栓式發(fā)動機。1965年開始研制推力為53~800 N的“月球跳躍者發(fā)動機”。同年開始研制C-1發(fā)動機,推力為440 N,該發(fā)動機可進行脈沖工作。1966年開發(fā)的URSA系列發(fā)動機,具有140、440、890 N 3種不同的固定推力,脈沖寬度小到0.02 s。20世紀80到90年代初生產(chǎn)了迄今為止最小的針栓噴注器,推力為22 N,僅重135 g[8]。

    根據(jù)“碰撞-殺傷”型導彈攔截器的需求,TRW利用面關機特性研制了36.5 kN的“哨兵”導彈俯仰偏航發(fā)動機,該發(fā)動機可在線性度為±8%條件下實現(xiàn)19∶1的推力變比,并可在任一推力水平產(chǎn)生小到8 ms的脈沖[14]。TRW還設計了一種具有面關機能力的針栓式發(fā)動機,該發(fā)動機用于空軍最早的“戰(zhàn)略防御主動動能武器”計劃之一,推進劑采用N2O4/MMH,脈沖響應為12 ms[15]。

    為了滿足約翰遜航天中心進行月球著陸器原型機的大氣飛行測試需求,普渡大學2012年設計測試了一個推力可調(diào)的液氧/液甲烷推力室,該推力室采用固定結(jié)構(gòu)的針栓噴注器,推力在5.85~18.9 kN范圍內(nèi)可調(diào)[16]。2017年以來,普渡大學的學生在該推力室的基礎上設計制造了一枚探空火箭,并基本完成發(fā)射前的準備工作[17]。美國SpaceX公司獵鷹9號火箭的梅林1D發(fā)動機也采用了針栓式噴注器,推力可在50%~100%范圍內(nèi)調(diào)節(jié),2015年12月該火箭一級部分首次實現(xiàn)了陸上垂直降落回收[18]。

    從國外發(fā)展與應用的情況來看,美國在針栓式發(fā)動機研究方面起步較早,目前處于領先地位,早期TRW公司在針栓噴注器和針栓式發(fā)動機方面處于主導地位,研制出了具備深度節(jié)流能力、推力涵蓋范圍大、具有面關機能力等特點的針栓式發(fā)動機,型號應用很多且比較成熟。目前SpaceX公司的多款發(fā)動機也采用了針栓噴注器,其發(fā)動機性能處于領先地位。

    1.2.2 國內(nèi)的發(fā)展與應用

    國內(nèi)在變推力火箭發(fā)動機領域的研究起步較晚。20世紀70年代國防科技大學率先開展針栓式發(fā)動機研究工作,并于1983年完成第一臺雙組元單調(diào)針栓式變推力發(fā)動機試車,推力調(diào)節(jié)范圍為245~1 221 N[19]。隨后研制了杠桿雙調(diào)針栓式變推力發(fā)動機SBF-03,推力為250~1 250 N[20-21]。上世紀90年代初,西安航天動力研究所成功研制了流量定位雙調(diào)針栓式變推力發(fā)動機,真空推力為250~1 250 N[22];本世紀初,西安航天動力研究所與西北工業(yè)大學等單位合作,成功研制了推力變比為10∶1的雙調(diào)變推力發(fā)動機[23],并在此基礎上發(fā)展了國內(nèi)目前唯一實際應用的針栓式變推力發(fā)動機,即嫦娥三號下降級7 500 N發(fā)動機,推力調(diào)節(jié)比為5∶1[24]。2018年,北京航空航天大學設計了一款基于機械定位雙調(diào)系統(tǒng)的氣氧/煤油針栓式變推力發(fā)動機,推力變比達15∶1,最高燃燒效率97.14%,但還沒有進行實際應用[10]。

    我國在針栓式發(fā)動機研究領域起步較晚,研制出的針栓式發(fā)動機推力較小,調(diào)節(jié)范圍有限,且大多還未進行實際的型號應用,目前我國唯一實際應用的針栓式變推力發(fā)動機為嫦娥三號下降級7 500 N發(fā)動機,其推力調(diào)節(jié)范圍為5∶1,與國外比還有較大差距。

    2 針栓噴注器霧化性能

    在液體火箭發(fā)動機工作過程中,推進劑經(jīng)過噴射、霧化、蒸發(fā)、混合反應、膨脹加速產(chǎn)生反推力來提供發(fā)動機的動力。噴注器的霧化性能對推進劑的蒸發(fā)、混合和燃燒有至關重要的影響,并直接影響發(fā)動機的性能。如圖4和圖5所示,針栓噴注器推力室流場結(jié)構(gòu)與平面噴注器推力室流場結(jié)構(gòu)有明顯的不同,平面噴注器推力室中推進劑沿推力室截面均勻分布,液滴軌跡近似沿軸線方向;而針栓式推力室中液滴初始軌跡與燃燒室軸線有很大的角度,推力室中存在兩個較大的回流區(qū),即頭部回流區(qū)和中心回流區(qū)?;亓鲄^(qū)的大小和位置受到噴射錐角等參數(shù)的影響,并影響推進劑的蒸發(fā)、混合及燃燒過程。平面噴注器推力室中液滴沿軸線運動且分布較均勻,液滴間發(fā)生相互碰撞的概率低,二次破碎的程度相對較低;而在針栓噴注器推力室中,由于液滴與室壁的撞擊作用及回流區(qū)的存在,液滴間發(fā)生相互碰撞的概率高,二次破碎程度相對較高。因此針栓噴注器和平面噴注器液流破碎霧化的動力學過程有所不同,需要進行深入研究。目前研究人員主要采用試驗、理論和數(shù)值仿真等方法開展相關的研究工作。

    圖4 平面噴注器發(fā)動機流場結(jié)構(gòu)示意圖

    圖5 針栓噴注器發(fā)動機流場結(jié)構(gòu)示意圖

    目前針栓噴注器按推進劑相態(tài)進行分類可分為兩類,即液/液和氣/液針栓噴注器,其中氣/液針栓噴注器為了實現(xiàn)較好的霧化效果,氣態(tài)推進劑一般都從外部的軸向通道噴出。下面從霧化過程、噴射錐角和霧化粒徑3個方面來討論這兩類噴注器的研究進展。

    2.1 霧化過程

    2.1.1 液/液噴注器

    實驗研究方面,Cheng等[25]對簡化噴注單元進行霧化實驗,結(jié)果表明液流霧化破碎模式與LMR(Local Momentum Ratio)及雷諾數(shù)密切相關,霧化場呈斗篷狀,如圖6所示。Marchi等[26]利用實驗研究了噴注器流道出口形狀對液膜穩(wěn)定性的影響,結(jié)果表明流道內(nèi)有臺階結(jié)構(gòu)時液膜更加穩(wěn)定。

    理論分析方面,Ninish等[27]對噴注器液膜進行理論分析,結(jié)果表明徑向射流的沖擊使得錐形液膜失穩(wěn),液膜表面形成擾動波(見圖7),擾動波增長導致液膜發(fā)生波動破碎。方昕昕等[28]基于小擾動假設,利用離心式噴注器錐形液膜破碎模型對針栓噴注器液膜破碎過程進行了線性不穩(wěn)定分析,結(jié)果表明正弦模式擾動波比曲張模式擾動波更不穩(wěn)定,在液膜破碎過程中占優(yōu);無量綱破碎時間和破碎長度隨壓降的增大均減小,隨噴注流道長寬比的增大近似線性增大,隨液膜半錐角的增大均減小,隨正弦模式加權(quán)因子的增大均減小。

    2.1.2 氣/液噴注器

    Son等[29]開發(fā)了二維軸對稱針栓噴注器尖端旋流場仿真方法,并通過徑向氣流速度分布作為標準與試驗結(jié)果進行對比驗證了該數(shù)值仿真方法的正確性,如圖8所示。

    圖8 徑向氣流速度分布的計算值與實驗值對比

    2.2 噴射錐角

    噴射錐角是針栓噴注器的關鍵性能參數(shù)之一,直接影響噴注器頭部和燃燒室壁面的熱環(huán)境,并且決定了回流區(qū)結(jié)構(gòu),影響推力室效率。

    2.2.1 液/液噴注器

    早期Santoro等[30]通過試驗研究了LOX/RP-1推進劑體系下噴射錐角的影響因素,發(fā)現(xiàn)其與TMR(Total Momentum Ratio)直接相關,而與阻塞因數(shù)幾乎無關。

    Freeberg等[31]以水為模擬介質(zhì)進行了液/液針栓噴注器冷流試驗,分析了噴射錐角、TMR(公式中簡記為TTMR)、LMR及針栓尖端噴孔幾何形狀之間的關系。結(jié)果表明,當建立噴射錐角和LMR之間的關系時,噴射錐角與噴孔的幾何形狀基本無關,這提供了一種預測噴射錐角的簡單方法。噴射錐角的理論研究開始的相對較晚,Ashgriz等[32]首先基于橫向射流噴射錐角模型提出了針栓噴注器的噴射錐角模型

    (1)

    Cheng等[33]在無黏不可壓條件下導出了液/液針栓噴注器的噴射錐角表達式

    θ=arccos[1/(1+TTMR)]

    (2)

    通過與數(shù)值仿真及試驗研究結(jié)果的對比驗證了式(2)的有效性,如圖9所示。

    圖9 噴射錐角隨總動量比的變化關系

    王凱等[34]基于三相CLSVOF方法(Coupled Level Set and VOF Method)對貼壁液膜/自由液膜撞擊的噴射錐角進行了數(shù)值模擬,并與自由液膜/自由液膜撞擊的噴射錐角進行了對比,揭示了入口TMR預測理論不適用于貼壁液膜/自由液膜撞擊的根源在于兩者撞擊形成的高壓區(qū)分布顯著不同導致垂直于壁面方向的動量不守恒。

    2.2.2 氣/液噴注器

    Son等[35]對氣/液針栓噴注器進行霧化試驗,噴注器結(jié)構(gòu)如圖10所示。試驗結(jié)果表明噴射錐角隨針栓開度的增加而增加,隨氣液速度比的增加而減??;在小流量工況下,可通過改變針栓錐度(θpt)而不能通過改變液體路針栓開度(Lopen)來改變噴射錐角,并利用無量綱數(shù)給出了該噴注器噴射錐角α的估算公式

    (3)

    其中

    s=1.15+1.35ξ

    p=1.30+0.90ξ

    Son等[36]研究了TMR和韋伯數(shù)We對氣/液針栓噴注器噴射錐角的影響,結(jié)果表明噴射半角隨著TMR/We的增加指數(shù)遞減,如圖11所示。Fang等[37]對針栓噴注器進行霧化試驗,針栓式發(fā)動機結(jié)構(gòu)如圖12所示。試驗結(jié)果表明:噴射錐角隨氣液流量比、氧化劑噴注縫寬度h0、液氧噴注半角α0的增加而減小,如圖13和圖14所示。Lee等[38]用液氮和氮氣模擬液氧和甲烷,對氣/液針栓噴注器進行低溫噴霧試驗,用陰影成像技術(shù)測量噴射錐角,結(jié)果表明噴射錐角隨TMR的增大而增大,當TMR一定時,針栓錐度(圖10中的θpt)增大,噴射錐角減小。

    圖11 噴射半角與TMR/We的關系

    圖12 針栓式發(fā)動機結(jié)構(gòu)示意圖

    圖13 氣液流量比與針栓錐度對噴射錐角的影響

    圖14 氣液流量比與噴注環(huán)縫寬度對噴射錐角的影響

    Radhakrishnan等[39]開展了噴注器二維數(shù)值仿真研究,湍流模型采用k-ε模型,一次和二次破碎分別使用單噴射模型和波動破碎模型,分析了噴射錐角及混合質(zhì)量隨參數(shù)NK=1/(1+TMR)的變化關系,發(fā)現(xiàn)隨著NK的增加,噴射錐角減小,混合質(zhì)量增加,且當NK大于0.83后,變化趨于穩(wěn)定,如圖15所示。

    圖15 噴射錐角和混合質(zhì)量隨NK的變化規(guī)律

    2.3 霧化粒徑

    噴霧場的液滴尺寸分布范圍較廣,常用分布函數(shù)來描述液滴尺寸。為方便起見,許多關于噴霧場的研究都采用液滴的平均直徑。噴霧場中液滴平均直徑的表示方法有很多,其中索泰爾直徑SMD(Sauter Mean Diameter)可以用來反映液霧燃燒性能的好壞,是燃燒流場中常用的霧化性能評價指標。

    2.3.1 液/液噴注器

    張紫豪等[9]采用馬爾文粒度儀對液/液針栓噴注器進行冷態(tài)噴霧粒徑測量,不同孔型測點分布如圖16所示。結(jié)果表明,徑向孔射流直接與軸向液膜作用時(2#孔型)噴霧的SMD要明顯小于兩孔射流交匯后再與液膜作用(1#孔型)的SMD,如圖17所示,分析認為撞擊所產(chǎn)生的波幅不同是影響SMD的主要原因。Ninish[27]等對液/液針栓噴注器的霧化試驗研究表明:軸向與徑向TMR越大,錐形液膜厚度和擾動波長越大,形成的液滴尺寸越大。

    圖16 不同孔型的測點分布圖

    圖17 不同孔型的SMD隨測點位置h的變化規(guī)律

    2.3.2 氣/液噴注器

    Son等[35]對氣/液針栓噴注器進行霧化試驗(噴注器結(jié)構(gòu)見圖10),試驗結(jié)果表明:SMD隨液體路針栓開度的增加而增加,隨對數(shù)韋伯數(shù)的增加而減小。參數(shù)敏感性分析表明,SMD對軸向流道間隙的變化最敏感,并利用無量綱數(shù)給出了該噴注器的SMD的估算公式

    SSMD=103(Lopen)(ξ-1)exp[4.0-q(We0.1)]

    (4)

    方昕昕等[40]研究了部分噴注器結(jié)構(gòu)參數(shù)對SMD和粒徑分布均勻度指數(shù)的影響,結(jié)果表明:SMD沿噴注器軸向變化不明顯,沿徑向有所增加;隨氣液流量比的增大,SMD減小,粒徑分布均勻度指數(shù)先減小而后增大;隨狹縫寬度的增大,SMD增大,粒徑分布均勻度指數(shù)減?。浑S液膜半錐角的增大,SMD減小,粒徑分布均勻度指數(shù)增大。

    3 針栓式發(fā)動機的燃燒流動問題

    在液體火箭發(fā)動機中,燃燒流動過程將推進劑的化學能轉(zhuǎn)化為燃氣的內(nèi)能,進而轉(zhuǎn)化為發(fā)動機的動能,在發(fā)動機工作過程中尤為重要。針栓式發(fā)動機具有特殊的流場結(jié)構(gòu),首先需要研究發(fā)動機、噴注器結(jié)構(gòu)參數(shù)和流動參數(shù)對燃燒性能的影響,評定燃燒效率;其次需要對發(fā)動機和噴注器的熱環(huán)境進行分析,避免發(fā)動機工作過程中各部件過熱損壞;最后還需要對發(fā)動機中的燃燒不穩(wěn)定問題進行分析,確保發(fā)動機安全可靠工作。

    對上述3方面問題進行研究的方法主要有實驗和數(shù)值仿真2種。關于針栓式發(fā)動機燃燒流動過程及整機性能的試驗研究中,嘗試了液氧/液甲烷[16]、液氧/煤油[18]、過氧化氫/甲醇[41]、液氧/酒精[42-43]、MMH/NTO[44]等不同推進劑體系。表1總結(jié)了針栓噴注器試驗過的所有推進劑組合。針栓式發(fā)動機推力室噴霧燃燒過程的數(shù)值仿真中,湍流模型一般采用k-ε模型,兩相反應流動中氣相方程在歐拉坐標系下建立,液相方程在拉格朗日坐標系下建立,兩相之間的耦合用氣液相互作用的源項考慮。

    表1 針栓噴注器試驗過的推進劑組合

    3.1 發(fā)動機及噴注器局部結(jié)構(gòu)對燃燒性能的影響

    3.1.1 局部結(jié)構(gòu)參數(shù)對燃燒性能的影響

    Austin等[42]對燃燒室特征長度L*、燃燒室直徑與噴注器直徑之比、TMR、噴注器長度與直徑之比等參數(shù)對發(fā)動機燃燒性能的影響研究(發(fā)動機結(jié)構(gòu)如圖18所示)表明:穩(wěn)態(tài)工作條件下燃燒性能對噴孔尺寸的變化不敏感,減小噴注器長度會提高燃燒效率;脈沖工作條件下L*略微增大可以提高燃燒效率,TMR對燃燒效率的變化無顯著影響。

    圖18 試驗發(fā)動機結(jié)構(gòu)

    周康等[45]的研究表明:對于氣氧/氣甲烷針栓式發(fā)動機,跳躍距離減小時,燃燒效率增大,火焰錐角不變;環(huán)縫寬度增加時,燃燒效率降低,火焰錐角增大;噴孔直徑增大時,燃燒效率基本不變,火焰錐角減小;火焰錐角在35°~68°范圍內(nèi),燃燒效率隨火焰錐角的增大而減小,如圖19所示。

    圖19 火焰錐角對燃燒效率的影響

    3.1.2 部分流動參數(shù)對燃燒性能的影響

    流動參數(shù)主要包括質(zhì)量流率、TMR、噴注壓降等。Sakaki等[43]對針栓式發(fā)動機燃燒室霧化及火焰結(jié)構(gòu)進行了光學觀測,在推進劑撞擊點觀測到了強烈的CH化學光(見圖20),并發(fā)現(xiàn)燃料中心方案比氧化劑中心方案的燃燒效率高。

    圖20 觀測到的瞬時火焰圖像(左)和CH化學光(右)

    李進賢等[46]的針栓式發(fā)動機內(nèi)流場仿真結(jié)果表明:質(zhì)量流率較大時,發(fā)動機燃燒性能越好;針栓深入燃燒室的長度增加,燃燒效率提高。

    Son等[47]對氣氧/甲烷針栓式發(fā)動機數(shù)值仿真結(jié)果表明:低質(zhì)量流率工況下火焰錐角減小會導致燃燒效率降低,此時減小針閥開度是提高燃燒效率的有效方法。俞南嘉等[3]對液氧/煤油針栓式發(fā)動機數(shù)值仿真結(jié)果表明,當TMR為1時,霧化混合效果最好,燃燒效率最高,如圖21所示。

    圖21 TMR對燃燒效率的影響

    張連博等[44]考慮流場與結(jié)構(gòu)之間的耦合傳熱效應(耦合計算模型如圖22所示),仿真結(jié)果表明:噴注器壓降對發(fā)動機的壓力、溫度和比沖的影響要比燃料和氧化劑的TMR對其影響顯著;對于燃料中心的針栓噴注器而言,改善燃料路的霧化效果更有利于推力室性能的提高。

    圖22 考慮結(jié)構(gòu)與流場之間耦合傳熱效應的針栓式發(fā)動機計算模型

    3.2 發(fā)動機熱環(huán)境分析

    Dmitry等[48]對液氧中心的液氧/液甲烷針栓式發(fā)動機熱試,研究表明:在針栓尖端用液氧進行主動冷卻不會降低燃燒效率,且能改善噴注器的熱環(huán)境;在噴注器軸向流道加導流板使得甲烷遠離針栓尖端區(qū)域,避免推進劑撞擊點發(fā)生在燃燒室中心區(qū)域,能顯著提高噴注器的壽命,但燃燒效率會降低。3種構(gòu)型的噴注器如圖23所示。

    圖23 3種構(gòu)型的針栓噴注器

    楊振寧等[49]對有面關機能力的針栓式發(fā)動機進行熱環(huán)境分析,結(jié)果表明:發(fā)動機在工作過程中頭部集液腔內(nèi)推進劑溫升很小,關機后推進劑有不超過200 ℃的溫升。

    3.3 發(fā)動機燃燒不穩(wěn)定問題

    Sakaki等[50]通過試驗研究了平面和軸對稱針栓噴注器的燃燒特性(平面和軸對稱針栓噴注器發(fā)動機結(jié)構(gòu)如圖24所示)。在軸對稱針栓噴注器發(fā)動機試驗中,當液氧噴注壓強小于燃燒室壓強時,觀察到了聲振,振幅超過40%,如圖25所示,當液氧噴注壓強達到燃燒室壓強的1.8倍時,振幅降到了20%。

    圖24 平面和軸對稱針栓噴注器發(fā)動機結(jié)構(gòu)

    圖25 試驗中觀察到的不穩(wěn)定燃燒現(xiàn)象

    Son等[47]對氣氧/氣甲烷針栓式發(fā)動機中的燃燒進行了數(shù)值仿真,研究了流動參數(shù)和幾何參數(shù)對燃燒特性的影響。研究表明當回流區(qū)太深時,會產(chǎn)生一個過熱點,這會引起冷卻問題和燃燒不穩(wěn)定。Kazuki等[51]對針栓式發(fā)動機燃燒室的縱向燃燒不穩(wěn)定性進行了試驗研究,將燃燒振蕩的時間尺度與燃燒室內(nèi)的各種現(xiàn)象進行了比較,以確定各燃燒振蕩模式的主導現(xiàn)象,結(jié)果表明TMR<1.14時,出現(xiàn)熱聲耦合振蕩,TMR>2時出現(xiàn)低頻振蕩。

    4 總結(jié)與展望

    針栓噴注器為大流量液體噴射問題提供了一個獨特的解決方案。由于針栓噴注器具有結(jié)構(gòu)簡單、可靠性高、燃燒穩(wěn)定、可深度節(jié)流等優(yōu)點,這類噴注器現(xiàn)已在液體火箭發(fā)動機領域獲得了廣泛應用,并表現(xiàn)出良好的應用潛力。然而,公開文獻中有關針栓噴注器的基礎研究(如噴注器霧化機理、燃燒穩(wěn)定性機理等)較少,其設計工作很大程度上依賴于工程經(jīng)驗(如將過熱部分材料替換為耐熱材料,依靠試驗確定使得霧化性能良好的噴注器關鍵參數(shù)等),因此需要對針栓噴注器進行深入研究,以提高其工作性能。

    4.1 關于針栓噴注器霧化性能

    針栓噴注器霧化性能方面的研究大致可分為霧化過程、噴射錐角、霧化粒徑3個方面。霧化過程研究方面,研究人員分析了液/液針栓噴注器的液膜破碎機理為徑向射流的沖擊在液膜表面形成不穩(wěn)定擾動波,擾動波增長導致液膜破碎,并認為正弦擾動波在液膜破碎過程中占優(yōu);對于氣/液針栓噴注器,液膜破碎理論分析較少。噴射錐角研究方面,研究人員分析了諸多參數(shù)對噴射錐角的影響,但改變這些參數(shù)本質(zhì)上是改變了噴注器的動量比,對于不同類型的針栓噴注器,動量比都是影響噴射錐角的主導因素。霧化粒徑研究方面,研究人員分析了部分參數(shù)對于液/液和氣/液針栓噴注器SMD的影響,分析認為,影響這2種噴注器SMD的主導因素也是動量比。

    針栓噴注器霧化性能研究集中于噴射錐角和霧化粒徑影響因素方面的研究,研究結(jié)果多為定性分析,缺乏定量研究。此外,對于噴注器液膜破碎霧化機理、霧化過程的研究較少。應建立針栓噴注器液膜不穩(wěn)定模型和霧化模型,研究液膜破碎霧化機理,并對霧化過程進行仿真模擬。

    4.2 關于針栓式發(fā)動機的燃燒流動

    針栓式發(fā)動機燃燒流動問題方面的研究可分為3個方面:發(fā)動機及噴注器局部參數(shù)對燃燒性能的影響、發(fā)動機熱環(huán)境分析和發(fā)動機燃燒不穩(wěn)定問題。

    發(fā)動機及噴注器局部參數(shù)對燃燒性能的影響方面,研究人員分析了部分結(jié)構(gòu)參數(shù)和流動參數(shù)對發(fā)動機燃燒效率的影響,但多為定性分析,沒有詳細分析影響機理。

    熱環(huán)境分析和燃燒不穩(wěn)定方面的研究較少。其中發(fā)動機熱環(huán)境分析方面,多為試驗測試,仿真模擬及理論分析較少,沒有深入分析熱環(huán)境發(fā)生變化的機理;燃燒不穩(wěn)定研究多為壓強振蕩模式的辨識,對燃燒不穩(wěn)定的機理性研究較少。

    4.3 目前需要重點關注的問題

    目前在針栓噴注器及針栓式發(fā)動機研發(fā)和應用過程中需要重點關注的主要問題包括:

    1)噴注器霧化問題。結(jié)構(gòu)參數(shù)和流動參數(shù)對噴注器霧化效果和發(fā)動機燃燒性能的影響基本只進行了定性研究,應進行定量分析和機理性研究,建立針栓噴注器霧化模型,同時研究霧化效果和燃燒效率對各參數(shù)的敏感性。

    2)發(fā)動機燃燒流動問題。發(fā)動機燃燒流動問題的研究集中于穩(wěn)態(tài)過程研究。發(fā)動機啟動過程和推力調(diào)節(jié)過程屬于非穩(wěn)態(tài)過程,對發(fā)動機的安全可靠工作極其重要,因此應研究非穩(wěn)態(tài)過程中發(fā)動機的工作特性。

    3)過熱問題。一方面是針栓尖端過熱問題。針栓尖端位于中心回流區(qū)的折轉(zhuǎn)位置,熱流密度高。另一方面是推進劑撞擊點附近過熱問題。對于自燃推進劑,發(fā)生在徑向射流根部小回流區(qū)的局部燃燒可能會導致針栓結(jié)構(gòu)的破壞。應對針栓尖端和撞擊點附近的熱環(huán)境進行分析研究,在保證燃燒效率的前提下通過合理的設計緩解針栓尖端和撞擊點附近的熱裕度壓力。

    4.4 發(fā)展趨勢分析

    針栓噴注器憑借其獨特優(yōu)勢在變推力發(fā)動機上具有極大的應用潛力,目前可回收利用是運載火箭的一個發(fā)展趨勢,將針栓噴注器用在可回收利用的運載火箭上是一個很好的選擇,因此針栓噴注器將來需要向大型化發(fā)展,另外,隨著環(huán)保意識的增強,針栓噴注器還將向環(huán)境友好型方向發(fā)展。

    猜你喜歡
    發(fā)動機研究
    FMS與YBT相關性的實證研究
    2020年國內(nèi)翻譯研究述評
    遼代千人邑研究述論
    2015款寶馬525Li行駛中發(fā)動機熄火
    2012年奔馳S600發(fā)動機故障燈偶爾點亮
    視錯覺在平面設計中的應用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    新版C-NCAP側(cè)面碰撞假人損傷研究
    新一代MTU2000發(fā)動機系列
    發(fā)動機的怠速停止技術(shù)i-stop
    亚洲精品久久成人aⅴ小说| 久久毛片免费看一区二区三区| 亚洲色图 男人天堂 中文字幕| 岛国毛片在线播放| 欧美在线黄色| 女同久久另类99精品国产91| 精品少妇内射三级| 淫妇啪啪啪对白视频| 亚洲国产av新网站| 亚洲精品自拍成人| videosex国产| 一个人免费看片子| 久久久精品区二区三区| 757午夜福利合集在线观看| 黑人猛操日本美女一级片| 精品国产一区二区久久| 丝袜人妻中文字幕| 1024视频免费在线观看| 极品教师在线免费播放| 午夜福利在线免费观看网站| 亚洲熟妇熟女久久| 男女床上黄色一级片免费看| 热99久久久久精品小说推荐| 涩涩av久久男人的天堂| 欧美日韩精品网址| 日韩免费高清中文字幕av| 亚洲色图综合在线观看| 国产精品 欧美亚洲| 一进一出抽搐动态| 又大又爽又粗| 欧美精品人与动牲交sv欧美| 精品国产一区二区三区久久久樱花| 一进一出好大好爽视频| 欧美日韩亚洲高清精品| 亚洲欧美激情在线| 女性被躁到高潮视频| 精品少妇内射三级| 丁香六月欧美| 人人澡人人妻人| 搡老乐熟女国产| 丰满人妻熟妇乱又伦精品不卡| 91麻豆av在线| 久久久精品国产亚洲av高清涩受| 欧美乱妇无乱码| 久9热在线精品视频| 大香蕉久久成人网| 久久久久久久久久久久大奶| 亚洲av第一区精品v没综合| 亚洲av欧美aⅴ国产| 夜夜爽天天搞| 日韩人妻精品一区2区三区| 一级片免费观看大全| 日韩欧美免费精品| 国产一区二区三区在线臀色熟女 | 91国产中文字幕| 又大又爽又粗| tocl精华| 18在线观看网站| a级毛片在线看网站| 久久精品亚洲av国产电影网| 亚洲一码二码三码区别大吗| 夜夜爽天天搞| 香蕉丝袜av| 天天添夜夜摸| 久久人妻福利社区极品人妻图片| 一夜夜www| 男女下面插进去视频免费观看| 男女边摸边吃奶| 视频在线观看一区二区三区| 最近最新免费中文字幕在线| 日韩人妻精品一区2区三区| 亚洲精华国产精华精| 久久亚洲真实| 最新的欧美精品一区二区| 国产成人影院久久av| 国产精品 国内视频| 精品免费久久久久久久清纯 | 高清毛片免费观看视频网站 | 国产激情久久老熟女| 色综合欧美亚洲国产小说| 正在播放国产对白刺激| 淫妇啪啪啪对白视频| av不卡在线播放| 色老头精品视频在线观看| 亚洲熟女毛片儿| 99久久人妻综合| 午夜福利视频在线观看免费| 国产日韩欧美在线精品| 一进一出好大好爽视频| 国产欧美日韩精品亚洲av| 亚洲精品美女久久久久99蜜臀| 国产在线视频一区二区| 国产精品.久久久| 亚洲av日韩精品久久久久久密| 日韩视频一区二区在线观看| 国产成人精品久久二区二区91| 欧美精品一区二区大全| 亚洲黑人精品在线| 又大又爽又粗| 香蕉国产在线看| 日韩有码中文字幕| 两人在一起打扑克的视频| 高清av免费在线| svipshipincom国产片| 久久国产亚洲av麻豆专区| av电影中文网址| a级毛片黄视频| 久久久久久免费高清国产稀缺| 伊人久久大香线蕉亚洲五| 少妇 在线观看| 超碰成人久久| av国产精品久久久久影院| 国产日韩一区二区三区精品不卡| 热99re8久久精品国产| 美女主播在线视频| 久久精品亚洲av国产电影网| 99riav亚洲国产免费| 满18在线观看网站| 国产有黄有色有爽视频| 国产三级黄色录像| 国产av精品麻豆| 露出奶头的视频| 十八禁网站网址无遮挡| 亚洲男人天堂网一区| 国产成人精品在线电影| 丝袜美腿诱惑在线| 国产亚洲精品久久久久5区| 欧美日本中文国产一区发布| 国产91精品成人一区二区三区 | 国产精品偷伦视频观看了| 亚洲七黄色美女视频| 久久中文字幕一级| 黑丝袜美女国产一区| 免费高清在线观看日韩| 国产欧美亚洲国产| 亚洲精品国产精品久久久不卡| 成人亚洲精品一区在线观看| 国产在线免费精品| 国产欧美亚洲国产| 少妇裸体淫交视频免费看高清 | 国产黄频视频在线观看| 国产av精品麻豆| 男女床上黄色一级片免费看| 老汉色av国产亚洲站长工具| 一个人免费在线观看的高清视频| 亚洲一码二码三码区别大吗| 交换朋友夫妻互换小说| 又大又爽又粗| 大香蕉久久成人网| 桃红色精品国产亚洲av| 久久性视频一级片| 久久久久久久精品吃奶| 久久久久久久久久久久大奶| 欧美一级毛片孕妇| 亚洲免费av在线视频| 一区二区三区国产精品乱码| 在线av久久热| 国产亚洲午夜精品一区二区久久| 大码成人一级视频| 美女主播在线视频| 国产精品成人在线| 一级,二级,三级黄色视频| 久久精品亚洲av国产电影网| 亚洲av国产av综合av卡| 丝袜人妻中文字幕| 国产精品秋霞免费鲁丝片| 欧美精品一区二区免费开放| 久久久久国产一级毛片高清牌| 18禁裸乳无遮挡动漫免费视频| 久久久久久亚洲精品国产蜜桃av| 日本wwww免费看| 美女主播在线视频| 黑人猛操日本美女一级片| 国产高清videossex| 久久影院123| 天堂8中文在线网| 在线观看免费高清a一片| 成人永久免费在线观看视频 | 欧美精品亚洲一区二区| 三上悠亚av全集在线观看| 国产日韩欧美亚洲二区| 欧美黑人精品巨大| 蜜桃在线观看..| 久久精品国产亚洲av高清一级| 亚洲国产欧美一区二区综合| 免费看十八禁软件| 欧美人与性动交α欧美精品济南到| 亚洲精品国产一区二区精华液| 50天的宝宝边吃奶边哭怎么回事| 在线观看舔阴道视频| 中文字幕制服av| 天天躁狠狠躁夜夜躁狠狠躁| 99re在线观看精品视频| 国产在线视频一区二区| 免费看十八禁软件| 99精品久久久久人妻精品| 免费不卡黄色视频| 高清在线国产一区| 一本大道久久a久久精品| 亚洲欧美日韩另类电影网站| 国产一区二区三区视频了| 亚洲国产欧美日韩在线播放| 精品一品国产午夜福利视频| 久久中文看片网| 久久中文字幕一级| 99精品久久久久人妻精品| 欧美成人免费av一区二区三区 | 国产欧美亚洲国产| 丝袜美腿诱惑在线| 欧美日韩福利视频一区二区| e午夜精品久久久久久久| 国产片内射在线| 高潮久久久久久久久久久不卡| 成人永久免费在线观看视频 | 精品亚洲成a人片在线观看| 亚洲精品久久成人aⅴ小说| 在线观看人妻少妇| 搡老岳熟女国产| 亚洲国产毛片av蜜桃av| 国产在线一区二区三区精| tube8黄色片| 1024香蕉在线观看| 久久久水蜜桃国产精品网| 免费人妻精品一区二区三区视频| tocl精华| 熟女少妇亚洲综合色aaa.| 国产在视频线精品| 亚洲精品国产区一区二| 亚洲精品国产区一区二| 在线观看人妻少妇| 19禁男女啪啪无遮挡网站| 狠狠婷婷综合久久久久久88av| 欧美精品人与动牲交sv欧美| 亚洲男人天堂网一区| 91成人精品电影| 在线av久久热| 久9热在线精品视频| 黄色视频在线播放观看不卡| 一边摸一边抽搐一进一出视频| 嫩草影视91久久| 久久精品国产a三级三级三级| 99国产极品粉嫩在线观看| 精品视频人人做人人爽| 久久精品国产亚洲av香蕉五月 | 老熟女久久久| 午夜福利一区二区在线看| 国内毛片毛片毛片毛片毛片| 这个男人来自地球电影免费观看| 黄频高清免费视频| 制服诱惑二区| 女人精品久久久久毛片| 亚洲伊人久久精品综合| 不卡一级毛片| 国产又爽黄色视频| avwww免费| 国产亚洲精品久久久久5区| 日韩人妻精品一区2区三区| 亚洲国产av新网站| 黄色丝袜av网址大全| 免费一级毛片在线播放高清视频 | 亚洲人成电影观看| 视频区图区小说| 99久久精品国产亚洲精品| 久久久国产一区二区| 亚洲精品av麻豆狂野| 99九九在线精品视频| 免费一级毛片在线播放高清视频 | 黑人巨大精品欧美一区二区mp4| 国产精品久久久久久精品古装| 欧美在线黄色| 国产精品.久久久| 99国产综合亚洲精品| 波多野结衣一区麻豆| 制服人妻中文乱码| 色婷婷久久久亚洲欧美| 久久久精品区二区三区| 久久精品aⅴ一区二区三区四区| 久久久精品免费免费高清| 国产免费福利视频在线观看| 国产精品偷伦视频观看了| 婷婷丁香在线五月| 国产精品一区二区免费欧美| 日韩中文字幕欧美一区二区| 午夜两性在线视频| 精品亚洲成国产av| 国产一卡二卡三卡精品| 国产精品一区二区在线观看99| e午夜精品久久久久久久| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 亚洲自偷自拍图片 自拍| 国产xxxxx性猛交| videosex国产| 大型av网站在线播放| 国产伦人伦偷精品视频| 99香蕉大伊视频| 大片电影免费在线观看免费| 亚洲欧美一区二区三区久久| 国产精品亚洲一级av第二区| 欧美激情 高清一区二区三区| 大陆偷拍与自拍| 久久天躁狠狠躁夜夜2o2o| 后天国语完整版免费观看| 欧美日韩成人在线一区二区| 国产精品免费大片| 法律面前人人平等表现在哪些方面| 色综合婷婷激情| 国产在线一区二区三区精| 人人妻人人爽人人添夜夜欢视频| 国产精品av久久久久免费| 国产色视频综合| av国产精品久久久久影院| 日韩中文字幕欧美一区二区| 巨乳人妻的诱惑在线观看| 国产精品偷伦视频观看了| 亚洲avbb在线观看| 成人永久免费在线观看视频 | 夜夜爽天天搞| 菩萨蛮人人尽说江南好唐韦庄| 我要看黄色一级片免费的| 97人妻天天添夜夜摸| 久久久久久免费高清国产稀缺| 国产高清激情床上av| 男女免费视频国产| 久久精品国产99精品国产亚洲性色 | 成人国产av品久久久| 99热网站在线观看| 久久午夜综合久久蜜桃| 欧美性长视频在线观看| 精品国产一区二区三区四区第35| 岛国毛片在线播放| 69精品国产乱码久久久| 少妇猛男粗大的猛烈进出视频| 精品国产亚洲在线| 精品视频人人做人人爽| 国内毛片毛片毛片毛片毛片| 成人三级做爰电影| 国产老妇伦熟女老妇高清| 999精品在线视频| 成人免费观看视频高清| 香蕉国产在线看| 国产成人免费无遮挡视频| 亚洲伊人色综图| av有码第一页| 老司机午夜十八禁免费视频| 国产精品久久久久久精品电影小说| 中文字幕制服av| 美女高潮到喷水免费观看| 中文字幕人妻丝袜制服| 老司机午夜十八禁免费视频| 欧美精品av麻豆av| 老司机亚洲免费影院| 欧美在线黄色| 黄频高清免费视频| 在线观看免费午夜福利视频| 国产福利在线免费观看视频| 国产成人av教育| 亚洲av日韩精品久久久久久密| 最新在线观看一区二区三区| 欧美国产精品一级二级三级| 麻豆乱淫一区二区| 狠狠精品人妻久久久久久综合| 国产一卡二卡三卡精品| 亚洲av第一区精品v没综合| 伦理电影免费视频| 黄色视频在线播放观看不卡| 午夜福利乱码中文字幕| 亚洲五月婷婷丁香| 久久人妻av系列| 欧美成人午夜精品| 日韩 欧美 亚洲 中文字幕| 99riav亚洲国产免费| 国产精品久久电影中文字幕 | 欧美精品一区二区大全| 少妇的丰满在线观看| 美女国产高潮福利片在线看| 不卡一级毛片| 亚洲全国av大片| 露出奶头的视频| e午夜精品久久久久久久| √禁漫天堂资源中文www| 另类精品久久| 国产欧美日韩一区二区三| 亚洲黑人精品在线| 一区二区av电影网| 9色porny在线观看| 水蜜桃什么品种好| 热99久久久久精品小说推荐| 中文欧美无线码| 午夜激情久久久久久久| 亚洲欧美激情在线| 日本精品一区二区三区蜜桃| 日韩大片免费观看网站| 亚洲专区字幕在线| 亚洲国产看品久久| 久久久国产成人免费| 91老司机精品| 亚洲全国av大片| 最近最新中文字幕大全免费视频| 国产精品美女特级片免费视频播放器 | 精品国产国语对白av| 色在线成人网| 丝袜人妻中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 亚洲伊人久久精品综合| 亚洲国产欧美网| 国产精品免费一区二区三区在线 | 91成人精品电影| 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区 | 叶爱在线成人免费视频播放| 久久ye,这里只有精品| 久久久久久免费高清国产稀缺| 蜜桃在线观看..| 国产精品亚洲av一区麻豆| 两性午夜刺激爽爽歪歪视频在线观看 | 久久九九热精品免费| 人人澡人人妻人| 人人妻,人人澡人人爽秒播| av超薄肉色丝袜交足视频| 多毛熟女@视频| 精品国产亚洲在线| 99在线人妻在线中文字幕 | 三级毛片av免费| 一本综合久久免费| 国产真人三级小视频在线观看| 国产成人一区二区三区免费视频网站| 久久精品亚洲精品国产色婷小说| av线在线观看网站| 免费观看人在逋| 人妻久久中文字幕网| 美女国产高潮福利片在线看| 欧美日本中文国产一区发布| 99久久国产精品久久久| 成人18禁高潮啪啪吃奶动态图| 欧美av亚洲av综合av国产av| 嫁个100分男人电影在线观看| 性少妇av在线| 一区二区三区激情视频| 人人妻人人澡人人爽人人夜夜| 亚洲成国产人片在线观看| 激情视频va一区二区三区| 久久av网站| 国产成人一区二区三区免费视频网站| 国产精品亚洲一级av第二区| 一个人免费看片子| 1024视频免费在线观看| 国产成+人综合+亚洲专区| 涩涩av久久男人的天堂| 国产成人影院久久av| 热99国产精品久久久久久7| 美女午夜性视频免费| 免费人妻精品一区二区三区视频| 少妇 在线观看| 超碰97精品在线观看| 美女午夜性视频免费| 亚洲性夜色夜夜综合| 国产成人影院久久av| 又黄又粗又硬又大视频| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 欧美黄色淫秽网站| 亚洲午夜理论影院| 亚洲成人免费av在线播放| 亚洲精品一二三| tocl精华| 国产在视频线精品| videosex国产| 欧美精品一区二区免费开放| 日日爽夜夜爽网站| 欧美黄色片欧美黄色片| 色婷婷av一区二区三区视频| 精品国产国语对白av| 国产免费av片在线观看野外av| 视频区图区小说| 国产成人欧美| 少妇裸体淫交视频免费看高清 | 成在线人永久免费视频| 18禁裸乳无遮挡动漫免费视频| 亚洲精品成人av观看孕妇| 亚洲欧美一区二区三区黑人| 在线播放国产精品三级| 午夜激情av网站| 在线永久观看黄色视频| 青草久久国产| 国产成人av激情在线播放| 久久毛片免费看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美一区二区三区久久| 免费看a级黄色片| 久久性视频一级片| 中文字幕精品免费在线观看视频| 91成年电影在线观看| 熟女少妇亚洲综合色aaa.| 亚洲一区中文字幕在线| 国产男女超爽视频在线观看| 久久九九热精品免费| 老司机亚洲免费影院| 99精品欧美一区二区三区四区| 2018国产大陆天天弄谢| 999精品在线视频| 亚洲av第一区精品v没综合| xxxhd国产人妻xxx| 精品视频人人做人人爽| 最新在线观看一区二区三区| 成人特级黄色片久久久久久久 | 露出奶头的视频| 日日摸夜夜添夜夜添小说| 国产淫语在线视频| 在线观看免费午夜福利视频| 丝袜在线中文字幕| 亚洲av日韩在线播放| 午夜福利视频精品| 少妇被粗大的猛进出69影院| 免费人妻精品一区二区三区视频| 国产日韩一区二区三区精品不卡| 成年动漫av网址| 日本av免费视频播放| 女人爽到高潮嗷嗷叫在线视频| 久久久水蜜桃国产精品网| 午夜久久久在线观看| 纯流量卡能插随身wifi吗| 后天国语完整版免费观看| 国产精品久久久久成人av| kizo精华| 高清在线国产一区| 老熟女久久久| 日本一区二区免费在线视频| 一区二区三区乱码不卡18| 日韩欧美一区二区三区在线观看 | a级片在线免费高清观看视频| 国产黄色免费在线视频| 下体分泌物呈黄色| 色94色欧美一区二区| 欧美国产精品一级二级三级| 午夜福利欧美成人| 91国产中文字幕| 美国免费a级毛片| 午夜福利在线免费观看网站| 国产精品欧美亚洲77777| 在线亚洲精品国产二区图片欧美| 欧美日韩视频精品一区| 久热爱精品视频在线9| 成在线人永久免费视频| 最近最新中文字幕大全电影3 | av不卡在线播放| 91精品三级在线观看| 久久午夜综合久久蜜桃| 在线天堂中文资源库| 91大片在线观看| 午夜福利乱码中文字幕| e午夜精品久久久久久久| 国产精品99久久99久久久不卡| 男女免费视频国产| 亚洲自偷自拍图片 自拍| 亚洲黑人精品在线| 精品久久久精品久久久| 一个人免费在线观看的高清视频| 在线观看免费高清a一片| 99在线人妻在线中文字幕 | 高清av免费在线| 高清视频免费观看一区二区| 亚洲情色 制服丝袜| 女性生殖器流出的白浆| 一区二区三区国产精品乱码| 性色av乱码一区二区三区2| 热99国产精品久久久久久7| 又大又爽又粗| 建设人人有责人人尽责人人享有的| 电影成人av| 97人妻天天添夜夜摸| 国产免费福利视频在线观看| 啦啦啦在线免费观看视频4| 日韩熟女老妇一区二区性免费视频| 精品高清国产在线一区| 成人免费观看视频高清| 免费日韩欧美在线观看| 亚洲一码二码三码区别大吗| 激情视频va一区二区三区| 1024视频免费在线观看| 精品久久久精品久久久| 国产精品av久久久久免费| 热re99久久国产66热| 免费观看av网站的网址| 一本大道久久a久久精品| 最新的欧美精品一区二区| 丰满饥渴人妻一区二区三| 精品久久久久久电影网| av线在线观看网站| 女人精品久久久久毛片| 国产精品久久久av美女十八| 久久影院123| 人人妻,人人澡人人爽秒播| 精品一区二区三区四区五区乱码| 999精品在线视频| 后天国语完整版免费观看| 欧美精品啪啪一区二区三区| 人成视频在线观看免费观看| 咕卡用的链子| 午夜福利乱码中文字幕| 脱女人内裤的视频| 黄片小视频在线播放| 在线 av 中文字幕| 日韩三级视频一区二区三区| 国产aⅴ精品一区二区三区波| 国产精品亚洲av一区麻豆| 两个人看的免费小视频| 老汉色∧v一级毛片| 亚洲av国产av综合av卡| 下体分泌物呈黄色| 成人手机av| 妹子高潮喷水视频| av福利片在线| 12—13女人毛片做爰片一| 午夜福利在线观看吧|