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

    深空探測用低溫推進劑貯箱熱力排氣系統(tǒng)研究

    2022-01-20 07:16:32楊云帆厲彥忠
    宇航學報 2021年11期
    關鍵詞:貯箱節(jié)流箱體

    劉 展,楊云帆,陳 虹,厲彥忠,3

    (1.中國礦業(yè)大學力學與土木工程學院,徐州 221116;2.航天低溫推進劑技術國家重點實驗室,北京 100028;3.西安交通大學能源與動力工程學院,西安 710049)

    0 引 言

    在未來載人小行星、載人登月和載人探火任務中,為增加有效載荷,有必要改進推進系統(tǒng)的整體性能。低溫推進劑[1-2],如液氫(LH2)、液態(tài)甲烷(LCH4)與液氧(LO2),憑借優(yōu)異性能,已被大規(guī)模應用于航天領域。然而,由于低溫推進劑貯存溫度較低,受熱易蒸發(fā),進而造成貯箱壓力升高,給低溫推進劑長期在軌安全貯存帶來極大隱患。因此,需采取合理有效的方式來減少推進劑的相變蒸發(fā),實現(xiàn)推進劑貯箱壓力的有效控制[2-3]。

    目前主要的低溫貯箱壓力控制方式[4]包括:直接排氣降壓、混合噴射降壓、熱力排氣系統(tǒng)(Thermodynamic vent system,TVS)控壓以及主動制冷控壓。對于直接排氣方式,在外部漏熱作用下,低溫貯箱壓力會持續(xù)升高,一旦升高到所設壓力上限,就通過排出部分氣相燃料來降低箱體壓力。待箱體壓力降到所設壓力下限時,排氣閥關閉。該方式雖然操作簡單,但所造成的排氣損失卻很大?;旌蠂娚浣祲悍绞绞窃谙潴w底部設置一噴射裝置,將箱體底部冷流體噴射到氣相區(qū),通過冷卻氣相的方式達到降低箱體壓力的目的。由于沒有冷量輸入,實際漏入箱內的熱量并沒有排出,該方式僅是把氣相區(qū)積聚的熱量轉移到過冷液相中,并不能從根本上消除外部漏熱,因此只適用于低溫推進劑短期貯存。對于TVS控壓以及帶有制冷機的主動制冷控壓方式,兩者中前者通過節(jié)流小部分流體產(chǎn)生冷量,并將冷量帶入貯箱,用以冷卻氣液相流體;后者則通過低溫制冷機產(chǎn)生制冷量,通過冷卻循環(huán)流體來達到降低貯箱壓力的目的。兩控壓方式都能從根本上消除外部漏熱,實現(xiàn)貯箱壓力控制。然而,基于目前科技水平,用于大型低溫貯箱主動制冷的低溫制冷機仍有待開發(fā),為獲得20.0 K的溫度,現(xiàn)有的逆布雷頓制冷機制冷量僅在幾十到幾百瓦,遠不能有效消除外部環(huán)境向低溫貯箱的漏熱;并且低溫制冷機體型、重量均較大,占地面積較多,這些都大大增加了發(fā)射成本,目前仍不適合低溫貯箱長期在軌壓力控制。然而,TVS通過節(jié)流小部分流體,利用該部分流體的汽化潛熱來冷卻箱內流體,僅以犧牲小部分流體為代價,就實現(xiàn)了貯箱壓力控制。只要初期設計中考慮到該部分排氣損失,在發(fā)射前,通過適當增加推進劑質量,就可以滿足低溫貯箱長期在軌運行需求。總的來說,采用TVS對低溫推進劑貯箱進行壓力控制具有較大的應用前景及開發(fā)潛力。

    本文詳細對比了國內外不同機構所開展的有關TVS方面的研究進展,梳理了TVS運行關鍵技術與主要影響因素,指出了中國在低溫推進劑貯箱TVS壓力調控方面的發(fā)展規(guī)劃。

    1 TVS工作原理

    圖1(a)展示了貯箱TVS壓力控制示意圖[5]??梢钥闯?TVS主要包括循環(huán)泵、Joule-Thompson(J-T)節(jié)流裝置、同心套管式換熱器、噴射棒及附屬設施。循環(huán)泵是整個系統(tǒng)運行的動力裝置;節(jié)流裝置是TVS產(chǎn)生冷量的關鍵部件;而同心套管式換熱器是將排氣潛熱合理利用,實現(xiàn)兩相流與單相流體換熱的核心部件。

    圖1 TVS結構示意圖與工作p-h圖Fig.1 Diagram of TVS and operation p-h diagram

    為實現(xiàn)箱體壓力控制,TVS主要有兩種工作模式。第一為混合降壓模式,此時箱內流體需具有一定過冷度。其工作原理為:當液相溫度未達到設定溫度而氣枕壓力達到TVS開啟壓力時,循環(huán)泵將低溫流體從貯箱底部抽出,穿過套管式換熱器(此時套管換熱器不工作),再經(jīng)噴射棒噴口呈放射狀噴入貯箱內部,通過抑制流體熱分層來降低貯箱壓力。該工作模式僅是實現(xiàn)箱內熱量從高溫氣相向低溫液相的轉移,并沒有從根本上消除外部環(huán)境漏熱。因此,隨著時間的持續(xù),箱內液體溫度將逐漸升高,當液體溫度升高到所設箱體壓力下限對應的飽和溫度時,再通過循環(huán)噴射將會造成液相的大量蒸發(fā),噴射近飽和液體已不能起到有效控制箱體壓力的目的。為此需開啟J-T節(jié)流制冷環(huán)路,進入TVS第二種工作模式,即節(jié)流制冷模式。該模式下,經(jīng)循環(huán)泵的流體分為兩部分,其中大股流體進入套管換熱器內管,小股流體通過J-T閥等焓膨脹,變成低溫低壓的氣液兩相流后進入套管換熱器環(huán)管,與大股流體進行換熱。在套管換熱器中,小股氣液兩相流在換熱器中相變蒸發(fā),自身被加熱成飽和或過熱氣體后排出箱體。而大股流體則被小股流體相變產(chǎn)生的冷量所冷卻,被冷卻的大股流體,在循環(huán)泵驅動下通過噴射棒送入箱體對箱內流體進行冷卻降溫。TVS通過小部分流體的汽化潛熱來消除漏熱,達到降低箱體壓力的目的。

    圖1(b)展示了TVS工作過程所對應p-h圖。這里僅對TVS第二種工作模式進行分析,認為箱內液相溫度已達到箱體壓力下限pmin對應飽和溫度時的狀態(tài)點1(p1),此時TVS節(jié)流制冷模式開啟。在外部漏熱下,箱體壓力從p1增加到壓力上限pmax(p2),循環(huán)泵與J-T閥同時打開。在泵的提升作用下,箱體壓力從p2增加到p3。之后,從箱體底部抽吸的流體被分為兩股,其中大股流體進入內管,小股流體通過節(jié)流閥,壓力從p3減小為p4,然后進入環(huán)管。兩股流體在換熱器中進行熱量交換。在內管中,大股流體從T3被冷卻到T5,然后經(jīng)噴射棒噴射到箱體,通過冷卻氣液相溫度來降低箱體壓力;而在環(huán)管中,小股流體從氣液兩相流變?yōu)轱柡蜌怏w或過熱氣體而排出箱體。至此TVS一個節(jié)流制冷工作周期結束。TVS持續(xù)工作,箱體壓力一直降低,當其降低到所設置的壓力下限時,TVS停止工作。在外部漏熱下,貯箱壓力將再次升高,當箱體壓力升高到箱體所設定壓力上限時,TVS再次啟動,并進入周期性啟停工作模式。有關TVS工作過程的熱力分析詳見文獻[6],此處不再贅述。

    2 TVS國內外研究現(xiàn)狀

    TVS概念是以航天為背景,最早由美國國家航空航天局NASA于20世紀90年代初提出,旨在為滿足低溫推進劑長期貯存而設計。

    DiStefano等[7]介紹了微重力條件下,采用節(jié)流換熱裝置冷凝貯箱氣枕以實現(xiàn)液體獲取的目的,此即為TVS的雛形,亦可稱為傳統(tǒng)型TVS。由于傳統(tǒng)型TVS自身結構的限制,其控壓性能與消除流體熱分層的能力均具有局限性。Flachbart等[5]對傳統(tǒng)型TVS進行了改進,逐步形成了噴射棒TVS(圖1a所示)。與傳統(tǒng)型TVS相比,噴射棒TVS采用低溫循環(huán)泵取代混合器、采用套管式換熱器替代了緊湊式換熱器。由于循環(huán)泵、J-T閥均布置在貯箱外部,噴射棒TVS不僅維修方便,而且設備產(chǎn)熱也未帶入箱內。再者,節(jié)流后的低溫流體通過貫穿于氣枕與液相之間的噴射棒徑向噴射,所以無論氣液相界面位置如何,噴射棒TVS均能有效地消除貯箱內部流體熱分層,實現(xiàn)箱體壓力控制。

    鑒于噴射棒TVS具有良好的控壓性能,本節(jié)將對國內外所開展有關噴射棒TVS的研究工作進行詳細梳理,以掌握相關技術的研究進展與現(xiàn)狀。

    2.1 國外研究現(xiàn)狀

    為實現(xiàn)低溫推進劑長期在軌壓力控制,目前國際上主要有美國、法國、日本以及韓國的科研機構開展了TVS的實驗及模擬研究。

    2.1.1美國

    以美國NASA為主的研究人員,在TVS實驗研究以及數(shù)值計算方面開展了大量研究。1994年,Nguyen[8]開發(fā)了一種能夠描述低溫貯箱TVS運行性能的集總參數(shù)模型。結果表明,在低溫貯箱自增壓階段,計算模型預測值偏高;而在TVS控壓階段,模型預測值與試驗結果吻合較好。

    20世紀90年代初,美國NASA搭建了TVS性能測試的多目標氫試驗平臺(Multipurpose hydrogen test bed,MHTB)[9](如圖2所示),并開展了不同工況的實驗研究。為模擬在軌工況,MHTB測試罐外表面采用發(fā)泡+多層隔熱材料處理,整個試驗罐放置于真空腔內。試驗過程中,真空腔內最低真空度達1.33×10-6Pa,真空腔溫度控制在80~300 K。

    圖2 NASA多目標氫測試平臺[9]Fig.2 Multipurpose hydrogen test bed built by NASA[9]

    1996~1998年,Hedayat等[10-11]在MHTB上開展了液氫LH2貯箱TVS控壓性能試驗。研究人員常采用充注液體體積與罐體總體積之比,即充注率來表示充注流體的質量。表1羅列了研究人員所開展的不同測試工況設置。結果表明:在18.7~54.1 W的模擬漏熱下,對不同流體充注率工況,噴射棒TVS(如圖3a所示)均可有效破壞箱內流體熱分層。當采用全GH2增壓時,TVS將貯箱壓力控制在131~138 kPa內。而當采用GHe增壓時,TVS可將貯箱壓力控制在207~214 kPa內。由于氦氣會阻礙LH2的蒸發(fā)和GH2的冷凝,從而抑制箱內降溫降壓過程,因此當采用GHe增壓時,J-T裝置要比全GH2增壓模式運行更長的時間才能將外部漏熱帶走。

    表1 LH2貯箱TVS不同運行機制[10-11]Table 1 Different operation codes of TVS in LH2 tank[10-11]

    為研究微重力下增壓GHe對LH2貯箱TVS控壓性能的影響,美國馬歇爾空間飛行中心[12]于2005年在MHTB上開展了為期兩周的試驗研究,所開展測試工況見表2。研究人員對比了單獨混合噴射以及泵+J-T節(jié)流制冷兩運行模式的控壓性能。結果表明:當漏熱為70 W時,在不同GHe增壓工況,TVS均能將箱體壓力波動控制在±3.45 kPa內。與僅依靠GH2增壓相比,采用GHe增壓,TVS控壓周期將延長37~68%,這與1998年所開展GHe增壓LH2貯箱實驗結果一致[11]。實驗結果還表明,對于低充注率工況,在氣枕熱分層明顯的情況下,僅通過開啟J-T閥就可實現(xiàn)對LH2貯箱的壓力控制。

    表2 GHe增壓LH2儲罐TVS測試工況[12]Table 2 TVS Test case in LH2 tank pressurized with GHe[12]

    2004年,Flachbart等[13]通過更換大流量循環(huán)泵以及節(jié)流裝置,在MHTB上開展了液氮LN2貯箱TVS控壓性能測試。試驗裝置平行安裝了一大一小兩個J-T閥用于節(jié)流制冷,具體如圖3b所示。所開展的測試工況見表3。通過對四種充注率工況進行試驗測試,結果表明:改進的TVS可有效控制箱體壓力,將箱內流體溫度變化由開始的1.7~7.5 K降低至0.6~0.8 K。當節(jié)流制冷模式開啟時,氣枕壓力線性下降,直到氣液溫度大致相等為止。另外,試驗還發(fā)現(xiàn):當流體充注率降低時,TVS控壓周期也降低。通過對兩J-T閥流量對比得出,小流量J-T閥使控壓周期減短,但泵的運行時間將變長,輸入功率增加;而大流量J-T閥則恰好相反。

    表3 LN2貯箱TVS測試工況[13]Table 3 Test conditions of TVS in LN2 tank[13]

    圖3 兩種噴射棒型TVS示意圖Fig.3 Schematic of two kinds of spray bar TVS

    當采用GHe增壓LN2貯箱時,研究人員[13]對液體充注率為50%與25%兩工況進行了TVS控壓性能測試。試驗中貯箱漏熱量為450 W,TVS控壓范圍為158.8~165.3 kPa。測試結果表明:無論GHe是否存在于貯箱氣枕區(qū),TVS均能有效的將貯箱壓力控制在設定范圍內。另外,Hedayat等[14]還采用一維模型對GHe增壓LN2貯箱TVS控壓過程進行了數(shù)值預測。對比發(fā)現(xiàn),采用一維模型預測的大部分流體飽和壓力與試驗結果存在較大偏差。

    2006年11月,Flachbart等[15]對裝有液態(tài)甲烷LCH4的低溫貯箱進行了13天TVS控壓實驗研究。試驗用TVS示意圖如圖3b所示。測試結果表明,在充注率為37%、貯箱總漏熱為715 W時,開啟TVS可使液相飽和壓力降低至44.2 kPa。向LCH4貯箱注入GHe可使箱體壓力升高至165.5 kPa,達到TVS控壓上限。由于試驗用TVS是針對LH2貯箱設計,其最大漏熱消除能力為55 W,而LCH4貯箱最大熱負荷高達715 W。因此對于LCH4箱體,TVS能夠有效控制氣枕壓力在設定范圍內變化,但不能有效抑制液相溫度的升高。

    2014年,VanDresar[16]針對GHe增壓液氧LO2貯箱,實驗研究了TVS的控壓性能。實驗工況見表4,主要對比了軸向噴射、頂部噴射以及軸向+頂部噴射相結合等三種模式對箱體壓力控制的效果。結果表明:組合模式可以將箱內氣液相溫度維持在合理的溫度波動范圍;而軸向噴射模式僅能夠將液體溫度保持相對穩(wěn)定,不能抑制氣相溫度的升高。頂部噴射模式在一定程度上能夠控制液體溫度變化。當不考慮GHe增壓LO2貯箱時,TVS可將箱體壓力維持在所設控壓帶范圍內。而當考慮GHe增壓LO2貯箱,TVS將箱體壓力控制在合理范圍內時,液相溫度將逐漸升高;如果TVS將流體溫度控制在合理波動范圍,箱體壓力將逐漸降低。

    表4 LO2貯箱TVS測試工況[16]Table 4 Test conditions of TVS in LO2 tank[16]

    截止2014年,美國NASA在MHTB上已對LH2、LN2、LCH4以及LO2進行了不同控壓范圍、不同漏熱水平、不同充注率、不同控壓模式以及不同增壓氣體等方面的試驗研究,積累了大量的控壓試驗實測數(shù)據(jù)與操作經(jīng)驗。

    2014~2016年間,美國NASA開展了部分TVS控壓方面的數(shù)值模擬研究[17-18]?;贛HTB上不同TVS控壓試驗,NASA構建了用于預測箱體自增壓過程以及TVS噴射降壓過程的歐拉—拉格朗日數(shù)值模型,著重考慮了氣液相間熱質交換、噴射流體粒子運動及換熱模型、流體流動湍流因子等因素對TVS控壓過程的影響。

    2016年,Majumdar等[19]基于通用流體系統(tǒng)模擬程序對MHTB上TVS工作過程建立了數(shù)值模型,模型中考慮了噴射棒與套管換熱器模型、不同換熱方式、箱內氣液相間以及與壁面間的換熱,預測了低溫推進劑貯箱TVS壓力控制過程。

    2.1.2法國

    2014年,法國國家空間研究中心(Centre national d’etudes spatiales,CNRS)以及格勒諾布爾阿爾卑斯大學的Thibault等[20]搭建了圖4所示的地面TVS實驗平臺。與噴射棒TVS不同的是,該實驗裝置中大股流體在被節(jié)流流體冷卻后,通過設置在氣相區(qū)的噴頭將冷卻流體噴射到箱體內部。實驗以氟化酮(NOVEC1230)為模擬工質,外部模擬漏熱為0~360 W,箱內氣液相溫度均為60 ℃。實驗結果表明,開啟TVS可有效降低箱體壓力,較好地實現(xiàn)外部漏熱的移除與箱體壓力的控制。

    圖4 法國TVS示意圖[20-23]Fig.4 Diagram of TVS in France [20-23]

    2016年,Mer等[21]提出了均相熱力模型,對帶有噴射回路的TVS控壓系統(tǒng)進行了優(yōu)化設計。在模擬漏熱為26.0 W、充注率為66%的設置下,采用NOVEC1230為工質,預測了TVS的控壓性能。結果表明:計算模型較好地反映了罐體冷卻降壓過程。

    同年,Mer等[22]提出采用主動絕熱技術控制箱體壁面漏熱,仍以NOVEC1230為測試工質,研究了箱體自增壓以及TVS主動制冷過程。實驗結果表明,對于給定的噴射流體溫度,箱體最終所達到的平衡狀態(tài)依賴于外部漏熱大小以及噴射流體質量流量;而對于冷卻降壓動態(tài)過程,TVS性能主要由罐體充注率以及噴射流體質量流量決定。

    2018年,Mer等[23]實驗研究了不凝性GN2對NOVEC1230罐體TVS控壓性能的影響。他們發(fā)現(xiàn),GN2的存在使得自增壓過程出現(xiàn)流體逆熱分層現(xiàn)象。而在TVS工作階段,噴射流體使得部分GN2溶入NOVEC1230,造成了罐體內部強烈的熱力不平衡現(xiàn)象。隨著噴射過程的持續(xù),罐內逐步趨于平衡狀態(tài)。

    不同于美國NASA所以開展的實驗,法國研究人員只研究了TVS單次節(jié)流制冷后箱體的降壓過程,并且在TVS工作一個周期后,氣相溫度出現(xiàn)近似指數(shù)變化的降溫曲線。這與美國NASA在MHTB上測試出的溫度周期性波動變化曲線是不同的。因此,法國所開展的試驗研究,并不能充分體現(xiàn)TVS對低溫推進劑長期在軌儲存具有良好的控壓效果。

    2.1.3日本

    為響應空間探測的需求,日本室蘭工業(yè)大學(Muroran Institute of Technology,MIT)開展了低溫貯箱TVS控壓的實驗研究。2018~2019年間,Imai等[24-26]搭建了如圖5所示的低溫TVS控壓測試裝置。不同于噴射棒TVS,Imai主要測試了低溫流體節(jié)流后,采用噴嘴的軸向噴射降壓以及采用噴射棒的徑向噴射降壓過程。實驗測試罐內徑為210.7 mm,高度為800 mm,測試工質為LN2,增壓氣體為GHe。通過開展多工況測試,結果表明:當循環(huán)流量為0.5~0.6 kg/min時,采用噴嘴軸向噴射可以消除流體熱分層。在較高的質量流量和較低的入口溫度下,TVS控壓效果顯著。另外,當進口溫度較低時,熱浮力較小,混合射流不易到達自由界面,TVS的降壓功能變弱。

    圖5 TVS裝置(a)與實驗流程圖(b)[24-26]Fig.5 Diagram of TVS device (a)and flow chart (b)[24-26]

    不同于美國NASA以及法國CNRS所開展的TVS實驗研究,Imai在實驗過程中,節(jié)流過冷的流體從LN2罐供給,以至于測試罐內流體液位逐漸升高。由于液位高度以及混合射流條件的不同,軸向混合噴射的流體無法到達罐體內部自由界面附近,不能較好實現(xiàn)箱體壓力控制。再者,實驗研究中射流速度流速和溫度不能獨立設定,且所開展的實驗工況較少,因此Imai所開展的實驗僅僅驗證了過冷流體噴射消除流體熱分層的能力,并不能充分驗證TVS的控壓性能。

    2.1.4韓國

    2018年,以韓國高等科學技術學院(Korea advanced Iinstitute of Science and Technology,KAIST)為主的研究單位開展了低溫貯箱TVS控壓性能實驗研究[27]。所搭建的測試裝置如圖6所示。整個裝置由測試罐、盤管換熱器、小型LN2罐以及真空泵組成。測試罐為置于真空腔中的銅制柱狀罐體,其內徑198 mm、高度216 mm;盤管換熱器平均直徑90 mm、長度3000 mm、換熱管內徑6.35 mm。實驗中采用LN2為工質,控壓范圍為120~140 kPa,所開展的實驗工況見表5。研究人員發(fā)現(xiàn):不同實驗工況均驗證了TVS的控壓效果。當TVS開啟時,貯箱壓力以及流體溫度均呈現(xiàn)出明顯的降低趨勢。這是因為本實驗中盤管換熱器置于貯箱氣枕區(qū),節(jié)流孔口浸沒在液體中。在TVS運行時,液體經(jīng)節(jié)流孔口節(jié)流制冷,兩相流進入盤管換熱器對氣枕進行充分冷卻,以致于貯箱壓力以及流體溫度均出現(xiàn)明顯降低。同時研究人員還構建了用于預測該過程的熱力預測模型,預測了盤管換熱器內壓力、溫度和流體出口溫度等參數(shù)變化。通過與實驗結果對比發(fā)現(xiàn),熱力模型預測結果與實驗數(shù)據(jù)存在一定偏差。

    表5 LN2貯箱被動TVS運行參數(shù)[27]Table 5 Operation parameters of passive TVS in LN2 tank[27]

    圖6 被動TVS裝置[27]Fig.6 Passive TVS device[27]

    通過與美國、法國以及日本等機構所開展的TVS實驗研究對比可知,韓國KAIST所搭建的實驗裝置最為簡單,僅通過液體節(jié)流冷卻氣枕來實現(xiàn)箱體壓力控制,該裝置僅能稱為被動型TVS,并不屬于噴射棒TVS。另外,縱觀四國所開展的低溫貯箱控壓測試,僅有法國以及韓國所開展實驗中出現(xiàn)了氣相冷凝與氣相溫度的明顯降低。分析論證發(fā)現(xiàn),這主要與實驗初始設置有關。兩實驗測試罐體內部初始氣液相溫度均設置在箱體壓力對應的飽和溫度附近,因此一旦有液體節(jié)流制冷量帶入氣枕,必然會引起氣相冷凝以及罐體壓力降低。而其他國家所開展的實驗工況中,測試罐體氣相區(qū)一般處于過熱狀態(tài),因此在TVS運行期間,氣相冷凝量較少,氣相溫度整體上呈相對穩(wěn)定的波動變化。

    2.2 國內研究現(xiàn)狀

    中國在低溫推進劑長期在軌貯存與管理方面的研究基礎相對薄弱,在載人登月以及深空探測相關背景的牽引下,低溫推進劑長期在軌貯箱壓力控制的工作才逐步開展。2016年之前,國內僅有蘭州物理研究所[28]、中國運載火箭研究院[29]、西安交通大學[30]、上海交通大學[31]以及北京控制工程研究所[32]等單位開展了TVS文獻調研與方案論證工作,有關TVS的實驗研究在2016年以后才逐步展開,主要的研究機構包括西安交通大學、上海交通大學、航天低溫推進劑國家重點實驗室以及中國運載火箭技術研究院。

    1)西安交通大學

    2016年,西安交通大學厲彥忠團隊[6]對TVS工作過程進行了熱力分析。同時,針對TVS同心套管換熱器,詳細考慮了內管單相流、環(huán)管兩相流以及管外自然對流等不同傳熱方式的耦合,提出了一種準穩(wěn)態(tài)換熱模型來預測TVS換熱器的傳熱性能。

    之后,厲彥忠團隊與中國運載火箭技術研究院聯(lián)合采用R123開展了TVS地面原理實驗研究[33-35]。實驗流程圖如圖7所示。實驗臺主要由測試罐、噴射棒換熱器、循環(huán)泵、節(jié)流閥、冷凝罐、緩沖罐等組成。測試罐體直徑與柱段高度均為800 mm;上下封頭高度為218 mm。罐體金屬壁為4.0 mm厚的304不銹鋼。套管換熱器長度為1.2 m,換熱器內管尺寸φ15.88×1.0 mm;外管尺寸φ22.2×1.2 mm;噴射棒尺寸φ9.52×1.0 mm。套管換熱器共設置四根噴射棒,每根噴射棒開有直徑為2.0 mm的12個單面噴口。基于該實驗裝置,兩單位主要開展了6種不同實驗工況(詳見表6)研究。通過改變循環(huán)流量、控壓帶以及外部熱負荷研究了TVS的運行性能,實驗驗證了TVS在實現(xiàn)箱體壓力控制以及消除外部漏熱的有效性和可行性。他們發(fā)現(xiàn):流體過冷度對推進劑長期貯存具有重要意義;在混合噴射降壓階段,單獨開啟循環(huán)泵噴射過冷流體即可達到箱體降壓的目的,并且在該階段,增大循環(huán)流量對實現(xiàn)箱體壓力控制以及氣枕冷卻具有積極意義;在TVS節(jié)流制冷階段,循環(huán)流量與節(jié)流比需合理設計。節(jié)流比較小時,不能對液體進行較好冷卻,以至于箱體升壓過快,TVS頻繁啟停。節(jié)流比較大時,換熱不充分使得排氣帶液,造成推進劑的浪費。與直接排氣相比,在相同初始設置下,TVS可節(jié)省30%的排氣損失。

    表6 采用R123所開展的TVS測試工況[33-35]Table 6 TVS test conditions of TVS conducted with R123[33-35]

    圖7 實驗流程圖[33-35]Fig.7 Experimental flow chart[33-35]

    2020年,厲彥忠團隊[36]針對TVS節(jié)流過程建立了熱力模型,深入揭示了焦湯節(jié)流制冷效應,分析了背壓對節(jié)流過程的影響,優(yōu)化了TVS運行參數(shù)。

    2)上海交通大學

    2016年,上海交通大學黃永華團隊[37-40]采用R141b開展了地面小型TVS實驗研究。測試罐直徑為450 mm、高度為785 mm、總容積約113L。實驗中采用纏繞紫銅絲的套管換熱器,換熱器直徑φ14.0×0.5 mm、外管尺寸φ19.0×0.5 mm、噴射棒直徑φ5.0×0.5 mm,每個噴射棒上布置了9個直徑φ2.0 mm的噴口。貯箱外部包裹4個加熱瓦,可用于模擬不同加熱方式?;趫D8所示實驗裝置,黃永華根據(jù)箱體壓力以及液相溫度設計了三種不同壓力控制策略,所開展的實驗工況見表7,對比研究了不同因素對TVS控壓性能的影響。結果表明:TVS可實現(xiàn)對R141b儲罐的壓力控壓。在TVS運行階段,大部分漏熱均可有效消除,并且J-T制冷階段的箱體壓增速率小于混合降壓階段的壓增速率。對于高充注率工況,TVS作用后,氣枕升壓速率變快,同時TVS單次工作時間變短,啟停切換更加頻繁。隨著貯箱控壓帶的增加,TVS排氣損失先降低后增加。

    表7 采用R141b開展的TVS工況設置[37-40]Table 7 Test condition of TVS with R141b[37-40]

    圖8 以R141b為工質的TVS測試流程[37-40]Fig.8 TVS Test flow chart with R141b[37-40]

    除實驗研究外,黃永華團隊還詳細考慮了氣液相變以及相間換熱,構建了TVS性能預測數(shù)學模型[41],并與美國NASA所開展的實驗結果進行對比,一定程度上驗證了TVS控壓模型的有效性。

    3)航天低溫推進劑國家重點實驗室

    2017年,上海交通大學黃永華教授在航天低溫推進劑國家重點實驗室所搭建的高效低溫流體貯存測試平臺(Efficient cryogenic fluid storage test platform,ECFSTP)上開展了LN2貯箱TVS控壓實驗研究[42]。實驗臺如圖9所示,其主要包括測試罐、冷屏、真空腔、MLI絕熱層、真空泵、換熱器、J-T閥等裝置。測試罐內徑與高度為1.0 m與2.0 m;罐內容積為1.36 m3。真空腔體真空度約0.01 Pa。在充注率為50%,外部漏熱為9.31 W/m2的設置下,對比了220~240 kPa、220~260 kPa、220~280 kPa、400~440 kPa與640~680 kPa五種不同控壓帶對LN2貯箱控壓性能的影響。實驗共設置了混合噴射、節(jié)流制冷兩種運行模式及兩種控壓策略。結果表明,當開啟TVS節(jié)流制冷模式時,過熱氣枕得到了充分冷卻,箱體壓力快速降低。相比于單純混合噴射模式,TVS節(jié)流制冷控壓效果更佳,TVS停止運行后箱體壓增速率明顯降低。對于不同控壓帶設置,TVS均表現(xiàn)出良好的控壓性能。隨著控壓帶范圍的增大,TVS的綜合性能得到改善。同時,黃永華等[43]還開展了控制策略對TVS性能的影響研究。他們發(fā)現(xiàn):與單純節(jié)流液相相比,TVS對氣液相節(jié)流可獲得更快的貯箱降低速率。TVS啟停判據(jù)應綜合考慮氣枕壓力以及液體飽和溫度的影響,否則將造成較大的排氣損失。與GN2增壓相比,增壓GHe增加了TVS的工作時間。然而,即使不凝性氣體存在,TVS仍具有良好的控壓效果。

    圖9 高效低溫流體貯存測試平臺[42-43]Fig.9 Efficient cryogenic fluid storage test platform[42-43]

    4)中國運載火箭技術研究院

    2018年,中國運載火箭技術研究院周振君[44]也在ECFSTP上開展了LN2貯箱增壓以及TVS控壓性能的研究。測試了直接排氣、節(jié)流制冷及并行模式下的貯箱內部熱力參數(shù)變化。實驗發(fā)現(xiàn):與直接排氣對比,采用節(jié)流制冷的并行模式可獲得更快的貯箱壓降,并且可以節(jié)省97%的排氣損失。

    2019年,周振君等[45]建立了低溫貯箱直接排氣控壓與TVS控壓計算模型,以排氣損失為研究目標,對比分析了在特定條件下LH2貯箱與LO2貯箱的最佳儲存方式。結果表明:相比于低溫LH2,LO2更容易長時間在軌貯存。對于LH2貯箱,可根據(jù)在軌貯存時間的長短,選擇TVS無損貯存、TVS排氣以及直接排氣等不同控壓模式。

    不同于ECFSTP被動TVS裝置,中國運載火箭技術研究院于2020年搭建了LN2噴射棒主動TVS控壓裝置[46],具體如圖10所示。測試裝置主要包括LN2貯箱、噴射棒TVS、數(shù)據(jù)采集系統(tǒng)、真空腔、LN2加注系統(tǒng)與真空系統(tǒng)。實驗中真空度可達0.01 Pa,共設置了三個控壓帶,研究了噴射混合模式與節(jié)流制冷模式對低溫貯箱控壓性能的影響。不同工況運行設置見表8。實驗結果表明:當采用單純噴射降壓模式時,貯箱升壓與降壓時間隨著控壓周期的增多而減小,并且液相溫度隨著時間的持續(xù)而波動升高。當采用TVS節(jié)流制冷模式時,貯箱壓力得到良好控制;在背壓不變時,制冷量隨著節(jié)流前壓力的增大而增大。另外,研究人員還發(fā)現(xiàn):當貯箱壓力降低到所設壓力下限時,TVS停止工作,此時箱體壓力仍繼續(xù)降低。也就是說TVS控壓周期的增壓階段具有一定的滯后。

    表8 不同工況運行參數(shù)[46]Table 8 Operation parameter of different conditions[46]

    圖10 液氮TVS測試平臺[46]Fig.10 Test device of TVS with

    2020年,上海宇航系統(tǒng)工程研究所也曾聯(lián)合上海交通大學開展地面液氮TVS控壓性能的實驗研究,但相關研究成果至今仍未報導。

    綜上可知,自2016年之后,國內有關TVS的實驗研究與理論分析開始走上快車道。西安交通大學與上海交通大學通過開展常溫流體噴射棒TVS性能實驗,總結并掌握了TVS運行過程中的關鍵技術、主要影響因素與調控策略,積累了一定的實驗操作經(jīng)驗。航天低溫推進劑國家重點實驗室搭建了ECFSTP平臺,并采用LN2開展了貯箱控壓測試,但因TVS換熱器浸沒在液體中,同時在氣相以及液相區(qū)設置了節(jié)流裝置,這與典型噴射棒TVS有所差別,因此仍屬于被動TVS。近期,中國運載火箭研究院所開展的低溫LN2貯箱壓力控制實驗才真正驗證了噴射棒TVS的控壓效果與消除外部漏熱的能力。

    3 對比分析

    圖11所示給出了不同國家在TVS控壓研究方面的關鍵節(jié)點與研究內容綜合對比。從中可以看出,從1990年提出TVS概念開始,美國已開展了長達30年TVS控壓方面的研究,研究內容涉及理論建模、試驗測試以及數(shù)值模擬。至今,美國NASA已開展了LH2、LN2、LCH4、LO2等不同流體TVS控壓性能試驗,同時研究了GHe對不同低溫燃料貯箱TVS運行性能的影響。通過20多年的試驗研究,美國已深入掌握了地面TVS運行的主要關鍵技術與影響規(guī)律,將TVS的技術成熟度指標(Technology readiness level,TRL)提高到5。

    圖11 國內外TVS控壓研究綜合對比Fig.11 Comprehensive comparison of TVS research at home and abroad

    法國采用常溫流體開展了TVS試驗研究,并且僅研究了TVS工作一次后的降壓效果,并沒有完全證實其長期運行的有效性。日本雖然采用LN2開展了低溫控壓試驗,但試驗裝置與常規(guī)噴射棒TVS不同,節(jié)流流體由LN2罐補給,造成測試罐體液位升高。雖然試驗驗證了噴射流體對熱分層的削弱能力,但沒有較好反映TVS的控壓性能。韓國搭建了被動TVS裝置,采用LN2測試了其控壓性能。以上三國所開展的TVS研究均與典型噴射棒TVS有所不同,并不能從真正意義上驗證TVS的控壓效果,因此其TRL在2~3內。

    中國在低溫貯箱壓力控制方面的研究起步較晚。在2016年以前,只有幾家研究機構開展了TVS的方案論證。而在2016年之后,國內主要的四家研究機構開始投入大量的精力從事TVS的理論與試驗研究。經(jīng)過近5年的努力,中國已搭建了不同功能的TVS裝置,測試工質也從常溫R123、R141b推廣到低溫LN2。同時開展了部分GHe增壓LN2貯箱的TVS控壓試驗,在一定程度上掌握了TVS裝置的設計方法與試驗關鍵技術,TRL值在3~4內。

    4 總結與展望

    4.1 主要結論

    通過對國內外在低溫貯箱用TVS所開展的30多年的研究成果整理分析,所獲主要結論如下:

    1)國內外不同試驗結果均表明:噴射棒TVS可實現(xiàn)對低溫貯箱壓力的良好控制以及外部漏熱的有效消除。被動型TVS在短時間內能夠起到良好的控壓效果,但不能抑制低溫流體溫度的波動升高,因此不適于低溫推進劑長期貯存。

    2)低溫燃料初始過冷度對其長時間在軌貯存具有重要意義。在低溫燃料處于過冷狀態(tài)時,僅通過開啟循環(huán)泵進行混合噴射就可以達到降低箱體壓力的效果;而當箱體內部大部分流體已接近飽和狀態(tài)時,此時需通過開啟TVS節(jié)流制冷模式進行降溫降壓。

    3)在TVS節(jié)流制冷模式中,當貯箱壓力降低到所設置壓力下限時,TVS停止工作,此時貯箱壓力仍繼續(xù)降低,貯箱升壓存在一定的滯后。這主要是由于過冷噴射流體與過熱氣枕來不及充分換熱所致。

    4)由于He分子量比H2大,其阻礙了LH2的蒸發(fā)和GH2的冷凝。因此,與GH2增壓LH2貯箱相比,采用GHe增壓LH2時,TVS節(jié)流制冷模式需要運行更長的時間才能達到良好的降壓效果。然而,在采用GHe增壓其他低溫流體時,對TVS工作效果影響不大。

    4.2 展望

    雖然中國在TVS研究方面取得了一定成績,但仍存在很多不足之處。例如,目前中國已搭建的低溫測試罐體積均較小,所測數(shù)據(jù)不能較好反映罐體內部熱物理場分布;低溫貯箱測量系統(tǒng)精度偏低,智能化程度不高;中國至今還沒有開展LH2、LCH4、LO2等低溫流體TVS控壓性能試驗。因此,為實現(xiàn)低溫推進劑長期在軌安全貯存,中國在TVS方面的研究仍有很長的路要走。

    盡管國內外針對TVS工作過程開展了深入的理論分析與實驗研究,掌握了TVS實驗的部分關鍵技術,但至今仍未開發(fā)出完善的TVS數(shù)值預測模型。另外,現(xiàn)有TVS實驗均為地面實驗,如何開展空間低重力工況TVS測試仍是國內外面臨的棘手難題。為促進地月轉移、載人小行星、深空探測等空間科學項目,實現(xiàn)低溫推進劑長期在軌貯存,以上研究仍是中國亟需開展的工作重點與努力方向。

    猜你喜歡
    貯箱節(jié)流箱體
    天然氣井井下節(jié)流器研究現(xiàn)狀及應用前景
    運載火箭貯箱補償器結構剛度的試驗研究
    超高壓氣井井下節(jié)流技術應用和設計方法
    基于Surface Evolver的推進劑貯箱氣液界面分析
    貯箱爆炸碎片初始速度及影響因素
    貯箱輕量化設計幾何參數(shù)優(yōu)化方法
    PR方程模擬節(jié)流效應的數(shù)值研究
    山東化工(2019年12期)2019-07-05 08:44:26
    高牌號灰鐵前端箱體質量提升
    “節(jié)流”是核心和重點
    超大型冷剪箱體加工難點分析
    中国美女看黄片| 国产亚洲精品综合一区在线观看| 韩国av在线不卡| 国产精品一二三区在线看| 麻豆精品久久久久久蜜桃| 夜夜夜夜夜久久久久| 性欧美人与动物交配| 18禁在线无遮挡免费观看视频 | 99热这里只有精品一区| 麻豆成人午夜福利视频| 1024手机看黄色片| 1024手机看黄色片| 一本一本综合久久| 国产精品久久视频播放| 亚洲精品国产av成人精品 | 国产精品国产高清国产av| 精品熟女少妇av免费看| 男人舔奶头视频| 日韩人妻高清精品专区| 婷婷精品国产亚洲av在线| 国内精品美女久久久久久| 精品不卡国产一区二区三区| 不卡一级毛片| 午夜福利高清视频| 99久久久亚洲精品蜜臀av| 欧美最新免费一区二区三区| 婷婷精品国产亚洲av| 成人亚洲欧美一区二区av| 精品不卡国产一区二区三区| 欧美不卡视频在线免费观看| 看免费成人av毛片| 国产精品国产高清国产av| 欧美日韩一区二区视频在线观看视频在线 | 夜夜夜夜夜久久久久| 精品无人区乱码1区二区| 成年女人毛片免费观看观看9| 国产精品伦人一区二区| 桃色一区二区三区在线观看| 成人av在线播放网站| 亚洲四区av| 麻豆成人午夜福利视频| 美女高潮的动态| 久久精品国产亚洲av天美| 亚洲av电影不卡..在线观看| 日韩欧美国产在线观看| 亚洲第一电影网av| 久久人人爽人人片av| 一级毛片我不卡| videossex国产| 美女xxoo啪啪120秒动态图| 亚洲最大成人av| 久久久精品欧美日韩精品| 嫩草影院精品99| 国产精品三级大全| 色综合站精品国产| 国产免费一级a男人的天堂| 国产男靠女视频免费网站| 国产毛片a区久久久久| 网址你懂的国产日韩在线| 又黄又爽又免费观看的视频| 亚洲av五月六月丁香网| 亚洲经典国产精华液单| 久久久久国产精品人妻aⅴ院| 少妇的逼水好多| 波多野结衣巨乳人妻| 国产精品人妻久久久久久| 欧美高清成人免费视频www| 亚洲欧美日韩高清在线视频| 欧美潮喷喷水| 欧美另类亚洲清纯唯美| 一进一出抽搐动态| 韩国av在线不卡| 成人二区视频| 亚洲国产精品国产精品| 中文在线观看免费www的网站| 精品久久久久久久久久免费视频| 精品国内亚洲2022精品成人| 中文在线观看免费www的网站| 午夜精品国产一区二区电影 | 国产黄片美女视频| 黄色配什么色好看| 国内精品一区二区在线观看| 黄色配什么色好看| 波野结衣二区三区在线| 国产亚洲精品综合一区在线观看| 看黄色毛片网站| 国产伦精品一区二区三区四那| 国内精品一区二区在线观看| 晚上一个人看的免费电影| 国产视频内射| 午夜爱爱视频在线播放| 男女视频在线观看网站免费| 少妇的逼好多水| 99在线人妻在线中文字幕| 成人漫画全彩无遮挡| 综合色av麻豆| 久久久精品欧美日韩精品| 亚洲七黄色美女视频| videossex国产| 18禁在线无遮挡免费观看视频 | 日本一本二区三区精品| 男女那种视频在线观看| 国内精品久久久久精免费| 真实男女啪啪啪动态图| 亚洲三级黄色毛片| 亚洲无线观看免费| 欧美最黄视频在线播放免费| 日韩欧美在线乱码| 看免费成人av毛片| 身体一侧抽搐| 亚洲精品在线观看二区| 日韩精品有码人妻一区| 国产精品av视频在线免费观看| 高清日韩中文字幕在线| 久久婷婷人人爽人人干人人爱| 美女黄网站色视频| 亚洲国产欧美人成| 日本精品一区二区三区蜜桃| 舔av片在线| 中国国产av一级| 国产黄色小视频在线观看| 91在线精品国自产拍蜜月| 午夜福利视频1000在线观看| 免费观看的影片在线观看| 99在线视频只有这里精品首页| av女优亚洲男人天堂| 最新在线观看一区二区三区| 欧美高清性xxxxhd video| 午夜久久久久精精品| 小蜜桃在线观看免费完整版高清| 亚洲,欧美,日韩| 欧美人与善性xxx| 国产午夜精品论理片| 亚洲专区国产一区二区| 国产午夜福利久久久久久| 成人三级黄色视频| 日韩欧美三级三区| 听说在线观看完整版免费高清| 久久6这里有精品| 男女之事视频高清在线观看| 97热精品久久久久久| 日韩亚洲欧美综合| 色视频www国产| 在线播放国产精品三级| 欧美潮喷喷水| 日韩精品中文字幕看吧| 悠悠久久av| 嫩草影视91久久| 一a级毛片在线观看| 少妇的逼好多水| 美女内射精品一级片tv| 亚洲成人久久爱视频| 日本撒尿小便嘘嘘汇集6| 尤物成人国产欧美一区二区三区| 亚洲图色成人| 成人精品一区二区免费| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品亚洲一级av第二区| 99国产极品粉嫩在线观看| 久久这里只有精品中国| 久久人人爽人人爽人人片va| 综合色丁香网| 一本久久中文字幕| 国产精品一二三区在线看| 国产 一区 欧美 日韩| 久久久久久久久久成人| 精品熟女少妇av免费看| 少妇人妻一区二区三区视频| 色综合亚洲欧美另类图片| 男女边吃奶边做爰视频| 日韩av在线大香蕉| 中文字幕精品亚洲无线码一区| 亚洲av电影不卡..在线观看| 亚洲人成网站在线观看播放| 久久久久免费精品人妻一区二区| 在线观看免费视频日本深夜| 日韩欧美精品v在线| 亚洲精品粉嫩美女一区| 乱人视频在线观看| 午夜老司机福利剧场| 久久精品国产亚洲网站| 国产aⅴ精品一区二区三区波| 伦精品一区二区三区| 国产精品久久久久久久电影| 无遮挡黄片免费观看| 亚洲av免费在线观看| 国产精品精品国产色婷婷| 亚洲国产色片| av黄色大香蕉| av专区在线播放| 人妻制服诱惑在线中文字幕| 身体一侧抽搐| 亚洲美女视频黄频| 国产av不卡久久| 最好的美女福利视频网| 精品一区二区三区人妻视频| 亚洲国产精品sss在线观看| 婷婷六月久久综合丁香| 亚洲无线在线观看| 亚洲四区av| 高清毛片免费看| 91在线观看av| 悠悠久久av| 亚洲一级一片aⅴ在线观看| 午夜亚洲福利在线播放| 欧美色欧美亚洲另类二区| 成人一区二区视频在线观看| 日韩精品中文字幕看吧| 日产精品乱码卡一卡2卡三| 亚洲综合色惰| 成人午夜高清在线视频| 国产黄色视频一区二区在线观看 | 亚洲av中文av极速乱| 亚洲精品一区av在线观看| 日韩欧美一区二区三区在线观看| 偷拍熟女少妇极品色| 99在线视频只有这里精品首页| 免费看av在线观看网站| 国产在视频线在精品| 欧美绝顶高潮抽搐喷水| 国产午夜福利久久久久久| 国产精品av视频在线免费观看| 黄色视频,在线免费观看| 日韩亚洲欧美综合| 成人鲁丝片一二三区免费| 99精品在免费线老司机午夜| 少妇的逼水好多| 大香蕉久久网| 国产av在哪里看| 九色成人免费人妻av| 久久国内精品自在自线图片| 蜜桃久久精品国产亚洲av| 97超碰精品成人国产| 国产极品精品免费视频能看的| 一a级毛片在线观看| 国产高清有码在线观看视频| 久久精品国产亚洲网站| 联通29元200g的流量卡| 国产老妇女一区| 国产一区二区在线观看日韩| 狂野欧美激情性xxxx在线观看| 最后的刺客免费高清国语| 亚洲人成网站在线播| 成人高潮视频无遮挡免费网站| 午夜福利在线在线| 精品乱码久久久久久99久播| 久久天躁狠狠躁夜夜2o2o| 夜夜夜夜夜久久久久| 欧美高清成人免费视频www| 亚洲欧美日韩无卡精品| 性插视频无遮挡在线免费观看| 男人狂女人下面高潮的视频| 在线播放无遮挡| 午夜福利18| 欧美日韩乱码在线| 欧美色欧美亚洲另类二区| 亚洲国产日韩欧美精品在线观看| 99久久久亚洲精品蜜臀av| 精品久久久久久久久av| 99在线视频只有这里精品首页| 亚洲在线自拍视频| 国产一级毛片七仙女欲春2| 国产成年人精品一区二区| 男插女下体视频免费在线播放| 中文字幕av在线有码专区| 3wmmmm亚洲av在线观看| 久久午夜亚洲精品久久| 国产成人一区二区在线| 久久久久久九九精品二区国产| 国产精品久久久久久精品电影| 少妇裸体淫交视频免费看高清| 啦啦啦啦在线视频资源| 日日摸夜夜添夜夜添av毛片| 日韩欧美精品v在线| 亚洲av熟女| 一级av片app| 不卡视频在线观看欧美| 久久国产乱子免费精品| 两个人视频免费观看高清| 成熟少妇高潮喷水视频| 黄色一级大片看看| 国产精品免费一区二区三区在线| 卡戴珊不雅视频在线播放| 欧美成人免费av一区二区三区| 国产亚洲精品综合一区在线观看| 亚洲精品日韩av片在线观看| 国产黄色小视频在线观看| 成年女人毛片免费观看观看9| 亚洲,欧美,日韩| 日本与韩国留学比较| 一级毛片久久久久久久久女| 三级毛片av免费| 欧美一级a爱片免费观看看| 午夜福利高清视频| 久久久久久久午夜电影| 国产美女午夜福利| 中国美白少妇内射xxxbb| 久久精品夜夜夜夜夜久久蜜豆| 欧美一区二区亚洲| 日本三级黄在线观看| 国产熟女欧美一区二区| 亚洲美女黄片视频| 欧美激情在线99| 国产乱人偷精品视频| 你懂的网址亚洲精品在线观看 | 不卡视频在线观看欧美| 午夜福利视频1000在线观看| 男女视频在线观看网站免费| 色在线成人网| 国产 一区精品| 欧美激情久久久久久爽电影| 丝袜喷水一区| 如何舔出高潮| 国产亚洲av嫩草精品影院| 国产一区二区三区在线臀色熟女| 深夜精品福利| 亚洲av美国av| 日本精品一区二区三区蜜桃| 亚洲自偷自拍三级| 国产一区二区三区在线臀色熟女| 久久午夜亚洲精品久久| 久久草成人影院| 可以在线观看毛片的网站| 日韩,欧美,国产一区二区三区 | 不卡一级毛片| 舔av片在线| 黄色欧美视频在线观看| 久久精品人妻少妇| 亚洲四区av| 天堂影院成人在线观看| 大型黄色视频在线免费观看| 亚洲四区av| 欧美人与善性xxx| 国产一区二区在线av高清观看| 成人二区视频| 99国产极品粉嫩在线观看| 97热精品久久久久久| 人妻制服诱惑在线中文字幕| 在线观看一区二区三区| 观看美女的网站| 一夜夜www| 免费av不卡在线播放| 亚洲欧美精品自产自拍| 日韩欧美精品v在线| 2021天堂中文幕一二区在线观| 国产精品一区二区三区四区久久| 国内精品宾馆在线| 国产激情偷乱视频一区二区| 看片在线看免费视频| 亚洲色图av天堂| 色综合亚洲欧美另类图片| 91久久精品国产一区二区成人| 少妇人妻精品综合一区二区 | 激情 狠狠 欧美| 最近2019中文字幕mv第一页| 久久精品国产亚洲av涩爱 | 一个人看的www免费观看视频| 少妇被粗大猛烈的视频| а√天堂www在线а√下载| 国产v大片淫在线免费观看| 亚洲av二区三区四区| 国产一级毛片七仙女欲春2| 五月伊人婷婷丁香| 色综合站精品国产| 最近的中文字幕免费完整| 别揉我奶头 嗯啊视频| 久久久久久久久久成人| 国产精品野战在线观看| 老司机福利观看| 亚洲美女搞黄在线观看 | 一个人看视频在线观看www免费| 网址你懂的国产日韩在线| 久久欧美精品欧美久久欧美| 欧美性猛交黑人性爽| 亚洲四区av| 免费黄网站久久成人精品| 欧美日韩国产亚洲二区| 国产精品1区2区在线观看.| 亚洲精品日韩av片在线观看| 在线看三级毛片| 久久久久国产精品人妻aⅴ院| 国产黄色小视频在线观看| 亚洲在线自拍视频| 亚州av有码| 国产69精品久久久久777片| 国产伦在线观看视频一区| 免费av毛片视频| 亚洲一区高清亚洲精品| 国产乱人视频| 美女xxoo啪啪120秒动态图| 老师上课跳d突然被开到最大视频| 亚洲美女视频黄频| 内地一区二区视频在线| 国产男靠女视频免费网站| 非洲黑人性xxxx精品又粗又长| 三级经典国产精品| 国产av不卡久久| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 欧美成人精品欧美一级黄| 亚洲欧美日韩无卡精品| 精品人妻熟女av久视频| 嫩草影视91久久| 毛片一级片免费看久久久久| av福利片在线观看| 国产成人a∨麻豆精品| 高清毛片免费观看视频网站| 夜夜看夜夜爽夜夜摸| 99riav亚洲国产免费| 精品不卡国产一区二区三区| 免费人成视频x8x8入口观看| 插逼视频在线观看| 日韩强制内射视频| 午夜视频国产福利| 国产国拍精品亚洲av在线观看| 国产久久久一区二区三区| 亚洲成人中文字幕在线播放| 99热这里只有是精品在线观看| 欧美日本亚洲视频在线播放| 国产亚洲精品av在线| 美女黄网站色视频| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 久久人人爽人人片av| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 少妇熟女aⅴ在线视频| 色噜噜av男人的天堂激情| 国产成人影院久久av| 国产亚洲精品综合一区在线观看| 中文亚洲av片在线观看爽| 色噜噜av男人的天堂激情| 嫩草影院新地址| 99久久成人亚洲精品观看| 亚洲精品亚洲一区二区| 久久人人爽人人片av| 国产高潮美女av| 国产麻豆成人av免费视频| 在线观看66精品国产| 性插视频无遮挡在线免费观看| av专区在线播放| 国产黄a三级三级三级人| 欧洲精品卡2卡3卡4卡5卡区| 香蕉av资源在线| 久久九九热精品免费| 97超碰精品成人国产| 欧美潮喷喷水| 精品一区二区三区av网在线观看| 蜜臀久久99精品久久宅男| 禁无遮挡网站| 天堂影院成人在线观看| 人妻夜夜爽99麻豆av| 一区二区三区四区激情视频 | 深夜精品福利| 在线看三级毛片| 成人高潮视频无遮挡免费网站| 精品乱码久久久久久99久播| 日韩欧美精品免费久久| 亚洲成a人片在线一区二区| 成人高潮视频无遮挡免费网站| 日本精品一区二区三区蜜桃| 赤兔流量卡办理| 久久久精品94久久精品| 精品少妇黑人巨大在线播放 | 亚洲成人av在线免费| 高清毛片免费观看视频网站| 国产精品综合久久久久久久免费| 亚洲国产精品成人久久小说 | 99久久中文字幕三级久久日本| 精品久久久久久久久久免费视频| 99热网站在线观看| 成人综合一区亚洲| 国产成年人精品一区二区| 国产精品野战在线观看| 国产精品爽爽va在线观看网站| 欧美绝顶高潮抽搐喷水| 日本爱情动作片www.在线观看 | 亚洲欧美清纯卡通| 亚洲专区国产一区二区| 日韩 亚洲 欧美在线| 色在线成人网| 国产精品伦人一区二区| 精品乱码久久久久久99久播| 亚洲人成网站在线观看播放| 久久久久久大精品| 看十八女毛片水多多多| 久久九九热精品免费| 久久精品91蜜桃| 一进一出抽搐动态| 自拍偷自拍亚洲精品老妇| 女人被狂操c到高潮| 香蕉av资源在线| 国产久久久一区二区三区| 日韩制服骚丝袜av| 大又大粗又爽又黄少妇毛片口| 国产精品人妻久久久影院| 尾随美女入室| 淫秽高清视频在线观看| 色综合色国产| a级毛片a级免费在线| 午夜视频国产福利| 99久久无色码亚洲精品果冻| 日韩,欧美,国产一区二区三区 | 久久久久久伊人网av| 免费观看精品视频网站| av免费在线看不卡| av专区在线播放| 天美传媒精品一区二区| 久久久久国产精品人妻aⅴ院| 国内精品宾馆在线| 女同久久另类99精品国产91| 春色校园在线视频观看| 亚洲乱码一区二区免费版| 精品人妻熟女av久视频| 午夜a级毛片| 亚洲欧美日韩高清在线视频| 哪里可以看免费的av片| 内射极品少妇av片p| 我要看日韩黄色一级片| av天堂中文字幕网| 日韩亚洲欧美综合| 国内少妇人妻偷人精品xxx网站| 国产单亲对白刺激| 村上凉子中文字幕在线| 成年av动漫网址| 免费av毛片视频| 亚洲在线观看片| 露出奶头的视频| 中文字幕av在线有码专区| 最近中文字幕高清免费大全6| 内地一区二区视频在线| 国产男人的电影天堂91| 熟妇人妻久久中文字幕3abv| 国产一区二区亚洲精品在线观看| 亚洲精华国产精华液的使用体验 | 欧美xxxx黑人xx丫x性爽| 国产女主播在线喷水免费视频网站 | 久久欧美精品欧美久久欧美| 狂野欧美白嫩少妇大欣赏| 国产精品av视频在线免费观看| 国产精品久久久久久久久免| 色吧在线观看| 成人精品一区二区免费| 亚洲精品国产成人久久av| 变态另类成人亚洲欧美熟女| 成年av动漫网址| 高清日韩中文字幕在线| 老司机午夜福利在线观看视频| 精品一区二区三区视频在线观看免费| 亚洲aⅴ乱码一区二区在线播放| 国产av麻豆久久久久久久| avwww免费| 男人舔女人下体高潮全视频| 禁无遮挡网站| 国产大屁股一区二区在线视频| av中文乱码字幕在线| 日韩欧美免费精品| av在线老鸭窝| 国产蜜桃级精品一区二区三区| aaaaa片日本免费| 日产精品乱码卡一卡2卡三| 国产人妻一区二区三区在| 一进一出抽搐gif免费好疼| 91久久精品国产一区二区成人| 色综合色国产| 久久午夜福利片| 成年版毛片免费区| videossex国产| 国产精品一区二区三区四区免费观看 | 久久99热这里只有精品18| 国产欧美日韩一区二区精品| 我的老师免费观看完整版| 最近中文字幕高清免费大全6| 一个人看视频在线观看www免费| 亚洲内射少妇av| 亚洲三级黄色毛片| 少妇被粗大猛烈的视频| 色尼玛亚洲综合影院| 高清日韩中文字幕在线| 在线播放国产精品三级| 免费电影在线观看免费观看| av女优亚洲男人天堂| 免费一级毛片在线播放高清视频| 国产一区二区亚洲精品在线观看| 国产91av在线免费观看| 中文字幕熟女人妻在线| 亚洲中文字幕日韩| 午夜免费激情av| 国产三级在线视频| 大又大粗又爽又黄少妇毛片口| 久久久久久久久久黄片| 日韩成人伦理影院| a级毛片a级免费在线| 久久精品91蜜桃| 如何舔出高潮| 精品久久久久久久久亚洲| 日日摸夜夜添夜夜添小说| 最近视频中文字幕2019在线8| 国产成人freesex在线 | 又爽又黄无遮挡网站| 精品一区二区免费观看| 97超碰精品成人国产| 黄片wwwwww| 免费人成在线观看视频色| 色播亚洲综合网| 亚洲人成网站在线观看播放| 女人十人毛片免费观看3o分钟| 亚洲精品在线观看二区| 国产v大片淫在线免费观看| 亚洲无线在线观看| 欧美区成人在线视频| 国产成人影院久久av| 99国产极品粉嫩在线观看| 青春草视频在线免费观看| 丰满的人妻完整版|