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

    低溫推進(jìn)劑在軌零蒸發(fā)貯存研究進(jìn)展

    2017-04-11 01:05:45王博杰
    載人航天 2017年2期
    關(guān)鍵詞:液氫貯箱制冷機(jī)

    王博杰,耑 銳,張 亮,齊 超,王 文

    (1.上海交通大學(xué)制冷與低溫工程研究所,上海200240;2.上海宇航系統(tǒng)工程研究所,上海201108)

    低溫推進(jìn)劑在軌零蒸發(fā)貯存研究進(jìn)展

    王博杰1,耑 銳2,張 亮2,齊 超1,王 文1

    (1.上海交通大學(xué)制冷與低溫工程研究所,上海200240;2.上海宇航系統(tǒng)工程研究所,上海201108)

    采用低溫制冷機(jī)與低溫貯箱耦合的方式,可以實(shí)現(xiàn)低溫推進(jìn)劑在太空環(huán)境下的零蒸發(fā)貯存。從系統(tǒng)分析、機(jī)理研究及試驗(yàn)研究三個(gè)方面對(duì)低溫推進(jìn)劑零蒸發(fā)貯存技術(shù)的現(xiàn)狀進(jìn)行了綜述,闡述了現(xiàn)有的系統(tǒng)分析方法、理論分析模型及零蒸發(fā)試驗(yàn)方案,可為低溫推進(jìn)劑空間貯存的發(fā)展提供參考。

    低溫推進(jìn)劑;零蒸發(fā);在軌貯存;壓力控制;熱分層

    1 引言

    低溫推進(jìn)劑可以應(yīng)用于上面級(jí)、探月、深空探測(cè)及空間燃料站等多個(gè)航天領(lǐng)域。但低溫推進(jìn)劑沸點(diǎn)低、極易蒸發(fā)、難以存儲(chǔ)等特點(diǎn)限制了其長(zhǎng)時(shí)間在軌使用。同時(shí),對(duì)于低溫推進(jìn)劑長(zhǎng)時(shí)間貯存的研究需要考慮空間微重力、復(fù)雜外熱流、各種輻射和粒子的影響并兼顧功耗低、重量輕的嚴(yán)格要求,使得其具有極大的技術(shù)挑戰(zhàn)性[1]。

    低溫推進(jìn)劑的共同特點(diǎn)是具有較低的飽和溫度,盡管絕熱技術(shù)的發(fā)展能使低溫貯箱的絕熱性能達(dá)到很高的水平,但還是因?yàn)槁釙?huì)導(dǎo)致低溫推進(jìn)劑升溫或蒸發(fā),致使箱內(nèi)壓力上升[2]。為確保低溫貯箱的安全,當(dāng)箱內(nèi)壓力上升到一定值時(shí),需要進(jìn)行排氣,這種控壓方式使得低溫液體在長(zhǎng)時(shí)間的在軌運(yùn)行中會(huì)大量損失[3]。另外,在太空環(huán)境中,低溫液體的排放有很多限制。首先,微重力環(huán)境使氣液界面不確定,若無(wú)氣液分離器,排氣過(guò)程中將有大量液體被排出;其次,排出的氣體由于微重力的影響會(huì)彌漫在航天器周圍,如果是易燃易爆的低溫推進(jìn)劑,將嚴(yán)重危害航天器以及宇航員的安全。另外,定期的排放使低溫液體的攜帶量遠(yuǎn)遠(yuǎn)大于額定需液量,這就直接影響到貯箱的質(zhì)量,增加了運(yùn)行負(fù)荷和發(fā)射費(fèi)用[4]。

    低溫推進(jìn)劑零蒸發(fā)貯存(Zero Boil?off,ZBO)的概念是隨著低溫制冷機(jī)與絕熱技術(shù)的發(fā)展被提出來(lái)的[5]。低溫推進(jìn)劑零蒸發(fā)貯存技術(shù)在航天器長(zhǎng)期運(yùn)行中較常規(guī)的排氣方式有很大的優(yōu)越性,其主要體現(xiàn)在對(duì)低溫推進(jìn)劑攜帶量與發(fā)射成本的收益,Plachta等[6]認(rèn)為在短期的貯存中,ZBO系統(tǒng)因?yàn)樾枰~外增加主動(dòng)冷卻裝置,較傳統(tǒng)被動(dòng)系統(tǒng)優(yōu)勢(shì)不大,但在貯存時(shí)間超過(guò)一定天數(shù)時(shí)(液氧是1周,液態(tài)甲烷是2周,液氫是2個(gè)月),ZBO系統(tǒng)在低溫推進(jìn)劑攜帶量與發(fā)射成本的收益會(huì)越來(lái)越大。

    NASA從二十世紀(jì)九十年代開始對(duì)ZBO進(jìn)行了理論與實(shí)驗(yàn)的研究,積累了很多經(jīng)驗(yàn)[7]。而我國(guó)在這方面尚處于起步階段,中國(guó)科學(xué)院及蘭州空間技術(shù)物理研究所等開展了相關(guān)的理論研究[8?10]。低溫推進(jìn)劑零蒸發(fā)貯存技術(shù)的研究對(duì)我國(guó)未來(lái)探月及深空探測(cè)意義重大。因此本文就低溫推進(jìn)劑零蒸發(fā)貯存技術(shù)的現(xiàn)狀從系統(tǒng)組成、理論機(jī)理以及試驗(yàn)研究等方面進(jìn)行了詳細(xì)地闡述與分析,指出了低溫推進(jìn)劑在軌零蒸發(fā)貯存所面臨的技術(shù)難題以及我國(guó)在該領(lǐng)域急需開展的關(guān)鍵技術(shù)與研究。

    2 低溫推進(jìn)劑零蒸發(fā)貯存概念

    低溫流體在長(zhǎng)期貯存過(guò)程中不可避免會(huì)受到各種漏熱帶來(lái)的困擾。低溫推進(jìn)劑零蒸發(fā)貯存的實(shí)現(xiàn)是將被動(dòng)絕熱技術(shù)與主動(dòng)制冷技術(shù)有機(jī)結(jié)合,通過(guò)低溫制冷機(jī)將貯箱系統(tǒng)漏熱完全補(bǔ)償,以實(shí)現(xiàn)貯存過(guò)程中低溫貯箱的壓力控制和低溫液體的零蒸發(fā)[3]。圖1是NASA描述的一個(gè)典型的零蒸發(fā)系統(tǒng)示意圖[7],系統(tǒng)中主要包括低溫制冷機(jī)、換熱器、攪拌風(fēng)扇以及多層絕熱結(jié)構(gòu)。低溫制冷機(jī)為系統(tǒng)帶走熱量;多層絕熱結(jié)構(gòu)用來(lái)保證系統(tǒng)的漏熱在制冷機(jī)最大制冷量之內(nèi);而攪拌風(fēng)扇一方面用于消除貯箱內(nèi)的熱分層,控制壓力水平,另一方面保證氣泡在零重力環(huán)境下也處于換熱器的位置,便于跟制冷機(jī)更有效地進(jìn)行換熱;另外在實(shí)際應(yīng)用中,因?yàn)橐WC制冷機(jī)的長(zhǎng)期運(yùn)行,需要增加太陽(yáng)能采集系統(tǒng)。在實(shí)際運(yùn)行中,當(dāng)貯箱內(nèi)壓力上升到一個(gè)閾值時(shí),攪拌風(fēng)扇開啟,用來(lái)打破熱分層,從而控制壓力,但攪拌混流的方式并不能快速降低壓力,當(dāng)僅靠攪拌風(fēng)扇無(wú)法再控制壓力時(shí),制冷機(jī)開啟,對(duì)貯箱內(nèi)低溫推進(jìn)劑進(jìn)行冷卻,從而達(dá)到對(duì)貯箱內(nèi)壓力真正的控制。在ZBO技術(shù)的發(fā)展中,箱體內(nèi)的混流裝置逐漸出現(xiàn)了循環(huán)泵、噴嘴或是噴桿等結(jié)構(gòu)[11?13]。而低溫冷卻屏技術(shù)(Broad Area Cooling,BAC)的提出,使箱內(nèi)冷卻與混流結(jié)構(gòu)可以得到一定程度的簡(jiǎn)化,Plachta等[14]在實(shí)驗(yàn)中發(fā)現(xiàn)僅通過(guò)高效的BAC系統(tǒng)能夠完全實(shí)現(xiàn)液氧的零蒸發(fā)貯存,以及液氫的低蒸發(fā)貯存。

    3 低溫推進(jìn)劑零蒸發(fā)貯存系統(tǒng)分析

    低溫推進(jìn)劑零蒸發(fā)貯存的系統(tǒng)分析是零蒸發(fā)系統(tǒng)設(shè)計(jì)與優(yōu)化重要的一部分,它包括針對(duì)各種不同任務(wù)要求,對(duì)ZBO系統(tǒng)總質(zhì)量、系統(tǒng)功耗等因素的分析。在文獻(xiàn)[3]中,一種簡(jiǎn)單的分析模型被用來(lái)進(jìn)行這方面的初步分析。模型給出了一系列ZBO系統(tǒng)中重要部分的質(zhì)量與能耗評(píng)估,包括低溫貯箱的體積與質(zhì)量的關(guān)系、多層絕熱結(jié)構(gòu)的性能與質(zhì)量的關(guān)系、蒸汽冷卻屏的質(zhì)量、以及低溫制冷機(jī)系統(tǒng)的冷卻性能與質(zhì)量及能耗的關(guān)系,如表1所示。

    Haberbusch等[19]通過(guò)上述模型對(duì)低溫推進(jìn)劑ZBO貯存進(jìn)行了系統(tǒng)分析。分析對(duì)象為球形液氫貯箱,貯箱外布置多層絕熱與冷卻屏,如圖2所示。作者分三種情況:(液氫存儲(chǔ)量250 kg和4000 kg,無(wú)外部冷卻屏;液氫存儲(chǔ)量250 kg,80 K的外部冷卻屏;液氫存儲(chǔ)量4000 kg,80 K的外部冷卻屏)討論了多層絕熱的層數(shù)以及內(nèi)部冷卻屏溫度對(duì)目標(biāo)ZBO系統(tǒng)質(zhì)量與能耗的影響。結(jié)果表明:1)對(duì)于同時(shí)使用內(nèi)部冷卻屏和外部冷卻屏的情況,外部多層絕熱達(dá)到75層時(shí)系統(tǒng)的質(zhì)量最??;無(wú)論是大儲(chǔ)液量(4000 kg)還是小儲(chǔ)液量(250 kg),系統(tǒng)能耗都會(huì)隨著多層絕熱層數(shù)的增加而降低,但當(dāng)多層絕熱超過(guò)100層時(shí),系統(tǒng)能耗趨于不變,如圖3所示;2)對(duì)于只使用內(nèi)部冷卻屏的情況,兩種儲(chǔ)液量在多層絕熱層數(shù)為150層時(shí)系統(tǒng)質(zhì)量達(dá)到最小;系統(tǒng)的能耗在大小儲(chǔ)液量下有很大的區(qū)別,對(duì)于小儲(chǔ)液量情況,系統(tǒng)能耗在多層絕熱層數(shù)150層時(shí)達(dá)到最小,而對(duì)于大儲(chǔ)液量情況,這個(gè)值要遠(yuǎn)大于150層,如圖4所示。這些結(jié)論為ZBO系統(tǒng)的優(yōu)化設(shè)計(jì)提供了一定的理論指導(dǎo)。

    表1 ZBO系統(tǒng)各部分質(zhì)量與能耗關(guān)系式Table 1 Mass and power relationships for ZBO system components

    4 零蒸發(fā)系統(tǒng)機(jī)理研究

    優(yōu)化ZBO系統(tǒng)使其性能達(dá)到最好需要對(duì)低溫貯箱內(nèi)的傳熱與運(yùn)輸機(jī)理的理解,ZBO系統(tǒng)的機(jī)理研究主要包括貯箱內(nèi)自增壓過(guò)程以及熱分層研究。

    4.1自增壓過(guò)程

    4.1.1 模型

    對(duì)低溫箱體內(nèi)增壓過(guò)程的研究始于20世紀(jì)50年代,最早采用的計(jì)算模型叫做均勻熱力學(xué)模型,模型認(rèn)為氣液兩相都處于飽和狀態(tài),且溫度一致[20]。這個(gè)假設(shè)過(guò)高估計(jì)了液相對(duì)平衡外部漏熱的貢獻(xiàn),所以均勻熱力學(xué)模型的壓力預(yù)測(cè)值往往低于實(shí)際值。之后Aydelott[21]提出了表面蒸發(fā)模型,認(rèn)為外部漏熱全部用來(lái)蒸發(fā),這顯然忽略了氣液相顯熱在整個(gè)能量平衡中的影響,所以表面蒸發(fā)模型的壓力預(yù)測(cè)值往往高于實(shí)際值。

    顯然,前面兩種模型在對(duì)箱內(nèi)增壓過(guò)程研究中過(guò)于理想,與實(shí)驗(yàn)數(shù)據(jù)的偏差也比較大[22]。之后的學(xué)者將質(zhì)量與能量的傳遞考慮到模型當(dāng)中,提出一系列的模型,統(tǒng)稱為多區(qū)域模型(mutizone model)。其中包括Riemer[23]提出的兩區(qū)域模型,他將氣液兩相分成兩個(gè)區(qū)域分別進(jìn)行研究;Estey等[24]在兩區(qū)域模型基礎(chǔ)上增加了獨(dú)立的氣液相界面區(qū)域,提出了三區(qū)域模型;Epstein和 Geor?gius[25]更是將氣相、液相以及箱壁分為數(shù)個(gè)軸向的區(qū)域進(jìn)行研究,提出了復(fù)雜的多區(qū)域模型;Schallhorn等[26]提出了邊界層模型,即在三區(qū)域模型基礎(chǔ)上增加了邊界層區(qū);Barsi等[27]在三區(qū)域模型基礎(chǔ)上增加了壁面區(qū)域。這一系列多區(qū)域模型大大豐富了箱內(nèi)增壓過(guò)程的研究,準(zhǔn)確性也較之前的模型有了很大的提升。

    但是上述的模型過(guò)于依賴區(qū)域之間換熱關(guān)聯(lián)式,不同的換熱關(guān)聯(lián)式得到的結(jié)果差異將會(huì)很大,所以選擇合適的換熱關(guān)聯(lián)式對(duì)模型的計(jì)算準(zhǔn)確性至關(guān)重要。

    4.1.2 CFD數(shù)值模擬

    考慮到之前的壓力模型過(guò)于依賴區(qū)域之間的換熱關(guān)聯(lián)式,同時(shí)也為了更好地理解低溫箱體內(nèi)增壓過(guò)程,更多的學(xué)者是采用CFD數(shù)值模擬的方式對(duì)箱內(nèi)增壓過(guò)程進(jìn)行研究。

    Panzarella和Kassemi[28]提出了lump?vapor模型,將其用于研究常重力下低溫貯箱的自增壓過(guò)程。這種模型將氣泡區(qū)集總為性質(zhì)一致的區(qū)域,而液體區(qū)通過(guò)數(shù)值模擬的方式聯(lián)立能量、動(dòng)量以及連續(xù)性方程來(lái)求解。模型考慮到了液體的熱分層可能帶來(lái)的影響,討論了三種不同加熱方式下貯箱內(nèi)的壓力升高過(guò)程,即對(duì)氣泡加熱、對(duì)液體加熱以及對(duì)整個(gè)貯箱加熱。與集總模型計(jì)算結(jié)果相比較,兩者具有相同的壓力增長(zhǎng)速度,但是壓力值有所偏差,這是因?yàn)樵谝后w區(qū)存在不同程度的熱分層對(duì)氣液界面的傳熱過(guò)程造成了一定影響。

    Panzarella等[29]還采用lump?vapor模型與數(shù)值模擬結(jié)合的方法對(duì)微重力下貯箱內(nèi)增壓過(guò)程進(jìn)行了研究。研究發(fā)現(xiàn)即使在微重力條件下,浮升力還是能很快將氣泡移動(dòng)至箱體頂部;另外,自然對(duì)流依舊可以在液體區(qū)造成相當(dāng)明顯的熱分層,兩者在微重力條件下依舊不容忽視。在對(duì)箱內(nèi)增壓過(guò)程預(yù)測(cè)中,在經(jīng)過(guò)初始的一段時(shí)間之后,箱內(nèi)壓力的變化速度趨于穩(wěn)定,與均勻熱力學(xué)模型預(yù)測(cè)結(jié)果接近,但因?yàn)槌跏茧A段的影響,最終的壓力預(yù)測(cè)結(jié)果與均勻熱力學(xué)模型還是有較大的差異。

    Venkat和Sherif[30]采用數(shù)值模擬的方法對(duì)低溫推進(jìn)劑儲(chǔ)存系統(tǒng)在常重力和微重力下的自增壓過(guò)程進(jìn)行了研究,采用Arnett和Voth提出的改進(jìn)模型[31]分析了不同的漏熱量、充液率及重力水平等因素對(duì)貯箱內(nèi)壓力增長(zhǎng)速度、蒸發(fā)速度以及熱分層速度的影響。研究結(jié)果表明,壓力增長(zhǎng)速度隨著漏熱量的增大而增大;在低充液率下的壓力增長(zhǎng)速度相比高充液率要快;另外重力加速度越小,貯箱內(nèi)壓力增長(zhǎng)速度與熱分層速度也越慢,這是因?yàn)橹亓铀俣仍叫?,?dǎo)致熱分層的自然對(duì)流就越弱。

    Barsi和 Kassemi[32?33]通過(guò)采用兩相的 CFD模型來(lái)研究常重力與微重力下液氫貯箱的自增壓過(guò)程。模擬結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比,認(rèn)為兩相的CFD模型能很好地預(yù)測(cè)不同充液率下的低溫貯箱自增壓過(guò)程。模擬結(jié)果在中等充液率水平下與試驗(yàn)結(jié)果出現(xiàn)了偏差,作者認(rèn)為這可能是不均勻的加熱負(fù)荷引起的。另外,均勻熱力學(xué)模型在壓力增長(zhǎng)速度方面有很好地預(yù)測(cè)性,而對(duì)于壓力值的計(jì)算上還有所偏差。

    4.2熱分層

    低溫箱體中緊貼箱壁的低溫流體在受到外熱流加熱后,密度變小,在自然對(duì)流的作用下,熱流體被帶到氣液界面處,之后向箱體中心擴(kuò)散并產(chǎn)生湍流擾動(dòng),擾動(dòng)產(chǎn)生的高溫流體渦團(tuán)會(huì)向下滲透,形成一定厚度的熱層,即所謂的熱分層現(xiàn)象[34]。熱分層會(huì)使局部蒸發(fā)加劇,導(dǎo)致貯箱內(nèi)壓力的過(guò)快增長(zhǎng)。對(duì)熱分層現(xiàn)象的研究始于二十世紀(jì)六七十年代,一系列的理論模型被提出用來(lái)研究熱分層現(xiàn)象的機(jī)理。

    Bailey等[35]將箱體內(nèi)部分為氣相區(qū)、液相主流區(qū)以及邊界層區(qū)來(lái)研究低溫流體的熱分層現(xiàn)象,并通過(guò)求解能量方程獲得了熱分層的溫度分布;Schmidt等[36]采用誤差函數(shù)來(lái)預(yù)測(cè)液氫箱體內(nèi)部界面溫度分布,并通過(guò)地面試驗(yàn)進(jìn)行了驗(yàn)證;Barnett等[37]采用雙曲線函數(shù)預(yù)測(cè)箱內(nèi)界面溫度分布,并通過(guò)開發(fā)新的模型來(lái)預(yù)測(cè)箱內(nèi)分層壓增等參數(shù)的變化;Neff等[38]采用自然對(duì)流半經(jīng)驗(yàn)公式,對(duì)低溫箱體在外部熱源作用下內(nèi)部的熱分層現(xiàn)象進(jìn)行了研究,指出熱分層的曲線可由無(wú)量綱容積與無(wú)量綱溫度擬合的多項(xiàng)式表示;Robbins等[39]基于湍流邊界層理論開發(fā)了一種新的熱分層模型,模型考慮了邊界層的脫離、擾動(dòng)造成的平均溫度降低等因素;Daigle等[40]針對(duì)液氫貯箱,提出了熱分層的簡(jiǎn)化熱力學(xué)模型,該模型將氣液相區(qū)分布劃分計(jì)算網(wǎng)格,并考慮了邊界層對(duì)熱分層的影響,采用集總參數(shù)法計(jì)算,并考慮氣液界面的熱質(zhì)交換,此模型在研究低溫箱體熱分層和壓力變化時(shí)具有良好的通用性,也是應(yīng)用比較廣泛的熱分層模型。

    熱分層模型的準(zhǔn)確性往往受到箱體尺寸、熱流大小等因素的影響,不同的模型適用范圍也各異,熱分層研究中涉及的瞬態(tài)溫度計(jì)算本身就是極其復(fù)雜的過(guò)程,所以需要通過(guò)更多的試驗(yàn)進(jìn)行驗(yàn)證,以及提出更完善、更具通用性的理論模型。

    5 主動(dòng)壓力控制過(guò)程研究

    主動(dòng)壓力控制過(guò)程控壓方式各異,控壓過(guò)程復(fù)雜,尚未有合適的模型對(duì)其進(jìn)行研究,一般都是采用CFD數(shù)值模擬對(duì)低溫箱體內(nèi)的壓力控制過(guò)程進(jìn)行研究?,F(xiàn)有的主動(dòng)壓力控制過(guò)程主要包括熱管方式、噴頭方式、噴桿方式等。

    5.1熱管主動(dòng)控壓方式

    通過(guò)熱管與低溫制冷機(jī)結(jié)合來(lái)實(shí)現(xiàn)低溫貯箱的壓力控制是一種較有效的方式。Ho等[41]采用CFD模擬對(duì)液氫貯箱進(jìn)行了瞬態(tài)熱分析,研究了液氫在主動(dòng)冷卻未開啟狀態(tài)與開啟狀態(tài)貯箱內(nèi)的溫度分布,如圖5所示。作者分三個(gè)階段對(duì)其進(jìn)行研究,第一階段,主動(dòng)冷卻未開啟時(shí),貯箱內(nèi)溫度最高的點(diǎn)位于貯箱角落;第二階段,主動(dòng)冷卻開啟時(shí),貯箱內(nèi)最高溫度點(diǎn)轉(zhuǎn)移到貯箱底部中間位置;第三階段,當(dāng)主動(dòng)冷卻剛關(guān)閉,尚存在微弱流體流動(dòng)時(shí),貯箱內(nèi)的最高溫度點(diǎn)又轉(zhuǎn)移到貯箱中間高度的箱壁上。這些預(yù)測(cè)結(jié)果為貯箱內(nèi)溫度監(jiān)測(cè)點(diǎn)的布置提供了參考,而這些溫度監(jiān)測(cè)點(diǎn)將是低溫貯箱間歇式運(yùn)行的依據(jù)。之后Ho等做了更細(xì)致的研究,研究發(fā)現(xiàn),當(dāng)噴管出流速度較大時(shí),箱內(nèi)流體的熱分層明顯消除,貯箱壓力得到有效控制,而處于下方的泵噴管較處于側(cè)面的泵噴管對(duì)貯箱內(nèi)液體混流和冷卻效果更好。

    5.2圓盤噴頭主動(dòng)控壓方式

    Ho等[12]在2006年提出了低溫制冷機(jī)與徑向圓盤噴頭結(jié)合的強(qiáng)制循環(huán)對(duì)流冷卻方式,并對(duì)其做了數(shù)值模擬研究,如圖6所示。這種結(jié)構(gòu)是由一根圓管一端連接換熱器,一端連接一個(gè)帶有多個(gè)噴嘴的圓盤,這些噴嘴以中心線為軸組成不同直徑的同心圓。當(dāng)貯箱內(nèi)壓力升高到一定值時(shí),一部分流體從貯箱頂部被引出,與制冷機(jī)冷端進(jìn)行換熱,再經(jīng)噴淋裝置回到貯箱中,這種方式為貯箱內(nèi)部帶來(lái)冷量同時(shí),還破壞了貯箱內(nèi)的熱分層。作者探討了入口管徑、噴頭深度、噴嘴間距等對(duì)系統(tǒng)性能的影響。模擬結(jié)果表明:入口管徑越大,表明循環(huán)流量也越大,系統(tǒng)的冷卻性能也越好,不過(guò)泵的能耗也相應(yīng)增大;圓盤噴頭位于貯箱中部時(shí),冷卻與控壓性能最好。

    5.3噴桿主動(dòng)控壓方式

    Kartuzova和Kassemi[42]在2014年針對(duì)NASA多功能氫試驗(yàn)床對(duì)噴桿方式的主動(dòng)控壓過(guò)程進(jìn)行了CFD模擬研究湍流模型采用的是k?ω模型,兩相模型采用的是VOF模型,對(duì)于液滴噴射過(guò)程采用的是拉格朗日噴射模型。噴射過(guò)程中考慮了有質(zhì)量交換和無(wú)質(zhì)量交換兩種情況,并對(duì)兩者進(jìn)行了比較,模擬的結(jié)果如圖7所示。從圖中可以發(fā)現(xiàn),考慮質(zhì)量交換的模擬結(jié)果更接近實(shí)驗(yàn)數(shù)據(jù),不考慮質(zhì)量交換的壓力降低速度較實(shí)驗(yàn)數(shù)據(jù)小很多,相應(yīng)地,有質(zhì)量交換的蒸發(fā)量較無(wú)質(zhì)量交換的小。這結(jié)果表明,在對(duì)噴桿主動(dòng)冷卻方式的模擬中,需要考慮到液滴與氣枕的質(zhì)量交換。

    6 零蒸發(fā)貯存試驗(yàn)研究

    從零蒸發(fā)貯存的概念提出以來(lái),NASA就做了一系列的地面試驗(yàn),主要針對(duì)的是低溫推進(jìn)劑中液氫與液氧的儲(chǔ)存,由于液氮與液氧性質(zhì)比較接近,故在很多液氧零蒸發(fā)試驗(yàn)中,都采用液氮作為替代工質(zhì)。

    6.1液氫試驗(yàn)

    6.1.1 低溫制冷機(jī)冷凝器與銅葉片混合冷卻

    Plachta等[2]在1999年對(duì)真空腔內(nèi)的球形液氫貯箱進(jìn)行了ZBO試驗(yàn)。如圖8所示,液氫貯箱直徑1.39 m,其中充注有90%的液氫。制冷機(jī)為兩級(jí)的商用GM制冷機(jī),其參數(shù)為第一級(jí)20 W/35 K,第二級(jí)17.5 W/18 K,布置于箱體頂部。第一級(jí)冷頭連接銅葉片,通過(guò)冷卻銅葉片來(lái)阻擋熱量的進(jìn)入;第二級(jí)冷頭連接一個(gè)冷凝器直接與箱內(nèi)氣枕接觸,通過(guò)直接冷卻氣枕的方式來(lái)降低箱內(nèi)壓力。試驗(yàn)結(jié)果表明,無(wú)論是混合使用冷凝器和銅葉片,還是單獨(dú)使用冷凝器,都可以實(shí)現(xiàn)液氫的ZBO儲(chǔ)存,而前者的壓力控制效率更高。

    6.1.2 低溫制冷機(jī)與噴桿結(jié)合制冷

    2002年,NASA的MSFC中心、ARC中心以及GRC中心共同合作在多功能氫試驗(yàn)床(MHTB)上,采用一臺(tái)30 W/20 K的GB37型低溫制冷機(jī),對(duì)容積為18 m3的液氫貯箱,采用噴桿強(qiáng)制循環(huán)對(duì)流的方式,進(jìn)行了ZBO試驗(yàn)[13]。如圖9所示,試驗(yàn)中低溫循環(huán)泵將液氫從貯箱中引出,與低溫制冷機(jī)換熱,經(jīng)過(guò)冷卻后的冷流體再經(jīng)過(guò)噴桿噴射回貯箱中,如此循環(huán),實(shí)現(xiàn)低溫貯箱內(nèi)流體的冷卻與壓力控制。在不同的液氫充液率工況(95%、50%、25%)進(jìn)行了試驗(yàn)研究,結(jié)果表明這種低溫制冷機(jī)與噴桿結(jié)合的冷卻方式可以實(shí)現(xiàn)液氫的零蒸發(fā)儲(chǔ)存。

    6.1.3 低溫冷卻屏系統(tǒng)

    低溫冷卻屏技術(shù)(Broad Area Cooling,BAC)首先由NASA格林研究中心在2005年提出,他們認(rèn)為在儲(chǔ)箱外包一個(gè)氣體冷卻屏,能大大減少儲(chǔ)箱內(nèi)低溫流體的蒸發(fā)量[43]。2014年,格林研究中心針對(duì)氫箱,僅采用BAC系統(tǒng)進(jìn)行了ZBO試驗(yàn)嘗試[44]。如圖10所示,試驗(yàn)采用15 W/90 K的逆布雷頓低溫制冷機(jī),BAC管路布置于絕熱結(jié)構(gòu)之間的冷卻屏上,采用傳統(tǒng)的多層絕熱與自支撐的多層絕熱結(jié)構(gòu)(Self?support MLI,SS?MLI)分別進(jìn)行了試驗(yàn)研究。因?yàn)榈蜏刂评錂C(jī)性能的制約,無(wú)法達(dá)到液氫的零蒸發(fā),但兩者分別有48%和45%蒸發(fā)量的減小。另外,試驗(yàn)結(jié)果表明SS?MLI較傳統(tǒng)MLI絕熱性能強(qiáng)18%左右。

    6.2液氧試驗(yàn)

    6.2.1 低溫制冷機(jī)與熱管結(jié)合冷卻

    在2004年NASA格林中心進(jìn)行的低溫推進(jìn)劑貯存ZBO實(shí)驗(yàn)中,箱體內(nèi)采用的就是熱管冷卻[45],如圖11所示。試驗(yàn)采用氮為工質(zhì),制冷機(jī)為NASA自主研發(fā)的高效能制冷機(jī)(HEC)在95 K有10 W制冷量?;旌媳迷?0 W功率下?lián)P程1.5 m,流速33 L/min。熱管蒸發(fā)段部分布置有大量的翅片,使其與箱體內(nèi)的流體充分換熱。并且在此套系統(tǒng)中,還布置了一個(gè)浸沒式泵進(jìn)行混流,用來(lái)破壞熱分層。試驗(yàn)結(jié)果表明,冷端性能良好,熱管表現(xiàn)出來(lái)的優(yōu)良性能使得制冷機(jī)安裝可以更遠(yuǎn)離貯箱;混合器增加的熱量小于預(yù)期值;低溫制冷機(jī)在散熱溫度311 K、制冷溫度75 K、功率消耗150 W條件下成功轉(zhuǎn)移出了6.8 W的熱量。

    6.2.2 低溫冷卻屏系統(tǒng)

    2014年NASA格林研究中心進(jìn)行了液氧貯箱的ZBO試驗(yàn)研究[44],采用的就是已經(jīng)發(fā)展完善的低溫冷卻屏技術(shù)。貯箱采用75層的傳統(tǒng)型多層絕熱材料(MLI)包裹,密度為24層/cm,BAC系統(tǒng)采用“Tube?on?Tank”布置,使用10根管徑6.4 mm的換熱管,BAC系統(tǒng)的制冷機(jī)采用15 W/90 K的逆布雷頓循環(huán)制冷機(jī),圖12是試驗(yàn)裝置。試驗(yàn)結(jié)果表明當(dāng)制冷機(jī)制冷量在8.5 W時(shí),功耗145 W,就能實(shí)現(xiàn)液氧的零蒸發(fā)儲(chǔ)存;如果進(jìn)一步提高輸入功率為212 W,實(shí)現(xiàn)零蒸發(fā)的同時(shí),它的壓力降低速率將大幅提高。這表明通過(guò)被動(dòng)絕熱與BAC系統(tǒng)已經(jīng)足以實(shí)現(xiàn)液氧的ZBO儲(chǔ)存與壓力控制。

    7 結(jié)論

    本文對(duì)ZBO技術(shù)的現(xiàn)狀進(jìn)行了概括與總結(jié),可以得出以下結(jié)論:

    1)對(duì)低溫推進(jìn)劑兩相過(guò)程的研究主要包括自增壓過(guò)程以及熱分層研究,影響它們的因素主要包括充液率、熱流邊界、工作壓力以及重力加速度等。研究者為此建立了一系列的壓力與熱分層模型,具有一定的準(zhǔn)確性,而在近幾年中研究者更多地采用兩相的CFD模擬來(lái)預(yù)測(cè)貯箱內(nèi)低溫推進(jìn)劑壓力變化及熱分層現(xiàn)象。目前主流的壓力控制方式有熱管方式、噴頭方式以及噴桿方式,主要通過(guò)CFD數(shù)值模擬對(duì)各種控壓方式進(jìn)行研究。

    2)現(xiàn)有ZBO技術(shù)主要包括低溫制冷機(jī)冷凝器與銅葉片混合冷卻技術(shù)、低溫制冷機(jī)與熱管結(jié)合冷卻技術(shù)、低溫制冷機(jī)與箱內(nèi)循環(huán)強(qiáng)制對(duì)流結(jié)合冷卻技術(shù)以及低溫冷卻屏技術(shù)。最新提出的低溫冷卻屏技術(shù)僅通過(guò)箱外冷卻就實(shí)現(xiàn)了液氧的零蒸發(fā)貯存。

    3)目前的ZBO試驗(yàn)當(dāng)中主要以液氫、液氧為主,液氫的零蒸發(fā)貯存較于液氧難度更高。從NASA目前的報(bào)道來(lái)講,由于20 K溫區(qū)的空間低溫制冷機(jī)技術(shù)成熟度不高,液氫的在軌零蒸發(fā)貯存至今無(wú)法實(shí)現(xiàn)。

    鑒于低溫推進(jìn)劑零蒸發(fā)貯存在未來(lái)深空探測(cè)及探月上的重要性,有必要對(duì)低溫推進(jìn)劑零蒸發(fā)貯存的原理與實(shí)現(xiàn)方式進(jìn)行更深入地研究。我國(guó)對(duì)于該領(lǐng)域尚處起步階段,為了更好地推進(jìn)我國(guó)航天事業(yè)的發(fā)展,除了機(jī)理研究之外,對(duì)于多層絕熱結(jié)構(gòu)、低溫制冷機(jī)等關(guān)鍵單機(jī)技術(shù)的研究同樣刻不容緩。

    (References)

    [1] 李鵬,孫培杰,包軼穎,等.低溫推進(jìn)劑長(zhǎng)期在軌儲(chǔ)存技術(shù)研究概述[J].載人航天,2012,18(1):30?36. Li P,Sun P J,Bao Y Y,et al.An overview on long?term storage of cryogenic propellant in?orbit[J].Manned Space?flight, 2012,18(1):30?36.(in Chinese)

    [2] Plachta D W.Hybrid thermal control testing of a cryogenic propellant tank[J].Advances in cryogenic engineering,2000,45:465?472.

    [3] Kittel P,Plachta D.Propellant preservation for mars missions[J].Advances in cryogenic engineering,2000,(45):443?450.

    [4] 李廣武,安剛,李娜,等.低溫液體無(wú)損儲(chǔ)存技術(shù)的發(fā)展與應(yīng)用[J].真空與低溫,2008,14(3):172?176. Li G W,An G,Li N,et al.Development and application of cryogenic liquid zero boil off storage technology[J].Vacuum&Cryogenics,2008,14(3):172?176.(in Chinese)

    [5] Plachta D W,Christie R J,Jurns J M,et al.Passive ZBO storage of liquid hydrogen and liquid oxygen applied to space science mission concepts[J].Cryogenics,2006,46(2):89?97.

    [6] Plachta D W,Kittel P.An updated zero boil?off cryogenic propellant storage analysis applied to upper stages or depots in a LEO environment[C]//38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference&Exhibit,Indianapolis,Indian?a,2002:3589.

    [7] Hastings L J,Plachta D W,Salerno L,et al.An overview of NASA efforts on zero boil off storage of cryogenics propellants[J].Cryogenics,2002,41:833?839.

    [8] 程進(jìn)杰,朱建炳,冶文蓮,等.基于ZBO存儲(chǔ)的低溫儲(chǔ)箱漏熱分析[J].低溫工程,2012,(6):46?50. Cheng J J,Zhu J B,Ye W L,et al.Heat leakage analysis of cryogenic tank based on ZBO storage[J].Cryogenics,2012,(6):46?50.(in Chinese)

    [9] 冶文蓮,王麗紅,王田剛,等.帶攪拌器的ZBO低溫貯箱三維模擬[J].低溫工程,2013,(6):5?9. Ye W L,Wang L H,Wang T G,et al.Three?dimensional a?nalysis of zero boil?off cryogenic storage tank used mixer[J]. Cryogenics,2013,(6):5?9.(in Chinese)

    [10] 程進(jìn)杰,孫郁,楊建斌,等.ZBO存儲(chǔ)低溫儲(chǔ)箱內(nèi)的壓力變化模擬分析[J].低溫與超導(dǎo),2014,42(1):17?20,34. Cheng J J,Sun Y,Yang J B,et al.Modeling the pressure variation law of cryogenic tank in the process of ZBO storage[J].Cryogenics&Superconductivity,2014,42(1):17?20,34.(in Chinese)

    [11] Plachta D.Results of an advanced development zero boil?off cryogenic propellant storage test[C]//40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit,F(xiàn)ort Lauderdale,F(xiàn)lorida,2004:3837.

    [12] Ho S H,Rahman M M.Zero boil?off cryogenic liquid hydro?gen storage tank with axial cold?spray system[C]//Proceed?ings of 2006 ASME International Mechanical Engineering Con?gress and Exposition,IMECE 2006?15341,Chicago,Illinois,2006:569?576.

    [13] Hedayat A,Hastings L J,Bryant C,et al.Large scale dem?onstration of liquid hydrogen storage with zero boiloff[J].Ad?vances in Cryogenic Engineering,2002,47:1276?1283.

    [14] Plachta D W,Christie R J,F(xiàn)eller J R,et al.Cryogenic boil?off reduction system testing[C]//50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference,Cleveland,OH,2014:3579.

    [15] Martin J J,Hastings L.Large?scale liquid hydrogen testing of a variable density multilayer insulation with a foam substrate[R].NASA Marshall Space Flight Center,NASA/TM,2001.

    [16] Ter Brake H J M.State of the art review on low power cryo?coolers[C]//Proceedings of the 19th International Cryogenic Engineering Conference.Alpha Science Int'l Ltd Grenoble,F(xiàn)rance,2003:307.

    [17] Humble R,Henry G,Larson W.Space Propulsion System A?nalysis and Design[M].New York McGraw?Hill,1995.

    [18] Moran M.Conceptual study of on orbit production of cryogenic propellants by water electrolysis[C]//27th Joint Propulsion Conference.Sacramento,1991:1844.

    [19] Haberbusch M S,Stochl R J,Culler A J.Thermally opti?mized zero boil?off densified cryogen storage system for space[J].Cryogenics,2004,44(6):485?491.

    [20] Aydelott J C.Self Pressurization of Liquid Hydrogen Tankage[M].New York:Cornell University,1967.

    [21] Aydelott J C.Normal gravity self?pressurization of 9?inch(23 cm)diameter spherical liquid hydrogen tankage[R].NASA/TN D?4171,1967.

    [22] Barsi S.Ventless Pressure Control of Cryogenic Storage Tanks[D].Cleveland Case Western Reserve University,2011:251.

    [23] Riemer D H.Cryogenic Tank Stratification a Simpler Ap?proach[M].New York:Advances in Cryogenic Engineering. Springer US,1986:957?962.

    [24] Estey P N,Lewis D H,Connor M.Prediction of a propellant tank pressure history using state space methods[J].Space?craft and Rockets,1983,20(1):49?54.

    [25] Epstein M,Georgius H,Anderson R E.A generalized pro?pellant tank?pressurization analysis[J].Int.Advances in Cryogenic Engineering,Plenum,1965:2:290?302.

    [26] Schallhorn P,Campbell D,Chase S,et al.Upper stage tank thermodynamic modeling using SINDA/FLUINT[C]//42nd AIAA/ASME/SAE/ASEE Joint Propulsion Conference&Ex?hibit,Sacramento,2008:5051.

    [27] Barsi S,Kassemi M.Investigation of tank pressurization and pressure control?part II:numerical modeling[J].Journal of Thermal Science and Engineering Applications,2013,5(4):041006.

    [28] Panzarella C H,Kassemi M.On the validity of purely thermo?dynamic descriptions of two?phase cryogenic fluid storage[J]. Journal of Fluid Mechanics,2003,484(484):41?68.

    [29] Panzarella C H,Plachta D,Kassemi M.Pressure control of large cryogenic tanks in microgravity[J].Cryogenics,2004,44(6):475?483.

    [30] Venkat S,Sherif S A.Self?pressurization and thermal stratifi?cation in a liquid hydrogen tank under varying gravity condi?tions[C]//Proceedings of the 42nd AIAA Aerospace Sciences Meeting and Exhibit,Reno,Nevada,2004:10844?10854.

    [31] Arnett R W,Voth R O.A computer program for calculation of thermal stratification and self?pressurization in a liquid hydro?gen tank[R].NASA CR?2026,1972.

    [32] Barsi S,Kassemi M.Numerical and experimental comparisons of the self?pressurization behavior of an LH2 tank in normal gravity[J].Cryogenics,2008,48(3?4):122?129.

    [33] Barsi S,Kassemi M.Validation of tank self?pressurization models in normal gravity[C]//45th AIAA Aerospace Sciences Meeting and Exhibit,Reno,Nevada,2007:952.

    [34] Clark J A.A review of pressurization,stratification and inter?facial phenomena[J].Int.Advances in Cryogenic Engineer?ing,1964,10:259?283.

    [35] Bailey T,Jefferson T,Skartvedt G,et al.Cryogenic propel?lant stratification analysis and test data correlation[J].AIAA Journal,1963,1(7):1657?1659.

    [36] Schmidt A F,Purcell J R,Wilson W A,et al.An experi?mental study concerning the pressurization and stratification of liquid hydrogen[J].Advances in Cryogenic Engineering, 1960,65C(2):487?497.

    [37] Barnett D O,Winstead T W,McReynolds L S.Investigation of liquid?hydrogen stratification in large cylindrical tank of Saturn configuration[J].Adv.Cryog.Eng,1964,10:314?323.

    [38] Neff B D,Chiang C W.Free Convection in a Container of Cryogenic Fluid[M].New York:Advances in Cryogenic En?gineering.Springer US,1967,12:112?124.

    [39] Robbins J H,ROGERS J.An analysis on predicting thermal stratification in liquid hydrogen[J].Journal of Spacecraft and Rockets,1966,3(1):40?45.

    [40] Daigle M J,Smelyanskiy V N,Boschee J,et al.Temperature stratification in a cryogenic fuel tank[J].Journal of Thermo?physics and Heat Transfer,2012,27(1):116?126.

    [41] Ho S H,Rahman M M.Transient analysis of cryogenic liquid hydrogen storage tank with intermittent forced circulation[J]. Journal of Thermophysics and Heat Transfer,2010,24(2):374?380.

    [42] Kartuzova O,Kassemi M,Agui J.Self?pressurization and spray cooling simulations of the Multipurpose Hydrogen Test Bed(MHTB)ground?based experiment[C]//50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference,Cleveland,OH,2014:3578.

    [43] Christie R J,Tomsik T M,Elchert J P,et al.Broad area cooler concepts for cryogenic propellant tanks[C]//Thermal&Fluids Analysis Workshop,Hampton,VA,2011:23.

    [44] Plachta D W,Johnson W L,F(xiàn)eller J R.Cryogenic boil?off re?duction system testing[C]//50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference,Cleveland,OH,2014:3579.

    (責(zé)任編輯:龐迎春)

    Research Progress of Zero Boil off Storage of Cryogenic Propellant On?orbit

    WANG Bojie1,ZHUAN Rui2,ZHANG Liang2,QI Chao1,WANG Wen1

    (1.Institute of Refrigeration and Cryogenics,Shanghai Jiao Tong University,Shanghai 200240,China;2.Shanghai Institute of Aerospace Systems Engineering,Shanghai 201108,China)

    A cryocooler integrated with the cryogenic tank can achieve zero boil?off storage of cryo?genic propellant in space.This paper summarized the theoretical and technical development of zero boil off storage in orbit,and listed several analysis methods,theoretical models and technical pro?grams,and analyzed their advantages and disadvantages.It may provide reference for the future de?velopment of cryogenic propellant storage in orbit.

    cryogenic propellant;zero boil?off;storage in orbit;pressure control;thermal stratification

    V51

    :A

    :1674?5825(2017)02?0236?09

    2016?06?14;

    2017?03?01

    王博杰,男,博士研究生,研究方向?yàn)榈蜏匾后w的儲(chǔ)存與管理。E?mail:bojiew@126.com

    猜你喜歡
    液氫貯箱制冷機(jī)
    運(yùn)載火箭貯箱補(bǔ)償器結(jié)構(gòu)剛度的試驗(yàn)研究
    甘肅隴西液氫生產(chǎn)及碳減排示范基地項(xiàng)目開工
    四川化工(2022年1期)2022-03-12 04:26:56
    3項(xiàng)液氫國(guó)家標(biāo)準(zhǔn)正式實(shí)施
    國(guó)家標(biāo)準(zhǔn)委批準(zhǔn)發(fā)布3項(xiàng)液氫國(guó)家標(biāo)準(zhǔn)
    全國(guó)液氫產(chǎn)業(yè)正式進(jìn)入快車道
    基于Surface Evolver的推進(jìn)劑貯箱氣液界面分析
    貯箱爆炸碎片初始速度及影響因素
    貯箱輕量化設(shè)計(jì)幾何參數(shù)優(yōu)化方法
    制冷機(jī)的制造與改造廣州圣嘉機(jī)電設(shè)備有限公司
    雙溫區(qū)雙冷指斯特林制冷機(jī)連管的設(shè)計(jì)
    真空與低溫(2015年4期)2015-06-18 10:47:26
    午夜福利在线观看免费完整高清在| 内地一区二区视频在线| 黄色欧美视频在线观看| 国产精品.久久久| 亚洲av国产av综合av卡| 欧美另类一区| 日韩亚洲欧美综合| 高清在线视频一区二区三区| 少妇被粗大的猛进出69影院 | 晚上一个人看的免费电影| 啦啦啦中文免费视频观看日本| 国产 精品1| 亚洲av在线观看美女高潮| 99国产精品免费福利视频| 伊人亚洲综合成人网| 亚洲av男天堂| 人妻人人澡人人爽人人| 丁香六月天网| 国产淫语在线视频| 欧美成人精品欧美一级黄| 观看av在线不卡| 日本免费在线观看一区| 久久久久国产网址| 十八禁网站网址无遮挡| 日韩欧美一区视频在线观看| 少妇精品久久久久久久| av在线老鸭窝| 精品亚洲乱码少妇综合久久| 久久久久久久大尺度免费视频| 欧美日韩亚洲高清精品| 18禁动态无遮挡网站| 亚洲精品乱码久久久v下载方式| 色5月婷婷丁香| a级毛片在线看网站| 爱豆传媒免费全集在线观看| 中文字幕精品免费在线观看视频 | 国产成人精品无人区| 亚洲一区二区三区欧美精品| 97超视频在线观看视频| 国产在视频线精品| 日本91视频免费播放| 成人免费观看视频高清| 日本wwww免费看| 久久精品久久精品一区二区三区| 精品99又大又爽又粗少妇毛片| 色5月婷婷丁香| 最新的欧美精品一区二区| 最近最新中文字幕免费大全7| 男人操女人黄网站| 亚洲精品456在线播放app| 国产精品99久久99久久久不卡 | 亚洲av成人精品一二三区| 日韩伦理黄色片| 亚洲欧美成人综合另类久久久| 搡女人真爽免费视频火全软件| 欧美日韩精品成人综合77777| 久久精品国产自在天天线| 69精品国产乱码久久久| 日本与韩国留学比较| 哪个播放器可以免费观看大片| 日产精品乱码卡一卡2卡三| 亚洲精品亚洲一区二区| 国产一区二区三区av在线| 久久ye,这里只有精品| 国产精品久久久久成人av| 亚洲av在线观看美女高潮| 亚洲,欧美,日韩| 18在线观看网站| 丰满乱子伦码专区| 狂野欧美白嫩少妇大欣赏| 亚洲第一av免费看| 亚洲精品aⅴ在线观看| 十分钟在线观看高清视频www| 精品卡一卡二卡四卡免费| 自线自在国产av| av有码第一页| 99久久中文字幕三级久久日本| 亚洲高清免费不卡视频| 精品99又大又爽又粗少妇毛片| 亚洲中文av在线| 中文欧美无线码| 欧美日韩精品成人综合77777| 99久久精品一区二区三区| 国产成人freesex在线| 人体艺术视频欧美日本| 伊人亚洲综合成人网| 久久99蜜桃精品久久| 欧美日韩成人在线一区二区| av在线观看视频网站免费| 免费高清在线观看视频在线观看| 日本av免费视频播放| 亚洲美女视频黄频| 成人国产av品久久久| 九草在线视频观看| 亚洲中文av在线| 边亲边吃奶的免费视频| 国产亚洲午夜精品一区二区久久| 国产视频内射| a级片在线免费高清观看视频| 这个男人来自地球电影免费观看 | 日韩成人伦理影院| 少妇被粗大猛烈的视频| 一级毛片我不卡| 老女人水多毛片| 欧美三级亚洲精品| 国产高清国产精品国产三级| 丰满饥渴人妻一区二区三| 久久久a久久爽久久v久久| 国产亚洲精品第一综合不卡 | 亚洲精品日韩在线中文字幕| 成人黄色视频免费在线看| 26uuu在线亚洲综合色| av在线播放精品| 亚洲精品国产av蜜桃| 国产免费现黄频在线看| 日韩av免费高清视频| 黄色视频在线播放观看不卡| 99国产精品免费福利视频| 只有这里有精品99| 国产一区亚洲一区在线观看| 夜夜看夜夜爽夜夜摸| 热re99久久国产66热| 午夜福利视频精品| 国模一区二区三区四区视频| 2021少妇久久久久久久久久久| 欧美xxxx性猛交bbbb| 欧美亚洲日本最大视频资源| 亚洲精品乱码久久久v下载方式| 亚洲三级黄色毛片| 国产精品99久久久久久久久| 国产成人精品福利久久| 欧美最新免费一区二区三区| 22中文网久久字幕| 日韩中字成人| 在线观看国产h片| 国产69精品久久久久777片| 不卡视频在线观看欧美| 日本wwww免费看| av不卡在线播放| 少妇猛男粗大的猛烈进出视频| 亚洲欧美日韩另类电影网站| 国产精品.久久久| 日韩精品有码人妻一区| av国产久精品久网站免费入址| 丁香六月天网| 国产精品偷伦视频观看了| 久久综合国产亚洲精品| 人人澡人人妻人| 欧美激情 高清一区二区三区| av视频免费观看在线观看| 成年人免费黄色播放视频| 欧美激情 高清一区二区三区| 亚洲成人手机| 欧美日韩国产mv在线观看视频| 午夜老司机福利剧场| 丝袜脚勾引网站| 91精品国产九色| 午夜福利视频在线观看免费| 18+在线观看网站| 国产精品99久久久久久久久| 一级毛片黄色毛片免费观看视频| 亚洲美女黄色视频免费看| 亚洲av.av天堂| 欧美性感艳星| 亚洲欧美日韩卡通动漫| 国产精品嫩草影院av在线观看| 高清毛片免费看| 制服丝袜香蕉在线| 色视频在线一区二区三区| 久久久久久久久久久丰满| 日本黄大片高清| 一级毛片电影观看| 亚洲人成网站在线播| 欧美三级亚洲精品| 国产亚洲最大av| 成人二区视频| 交换朋友夫妻互换小说| 国产成人a∨麻豆精品| 亚洲av成人精品一区久久| 日韩人妻高清精品专区| 毛片一级片免费看久久久久| 两个人免费观看高清视频| 国产午夜精品久久久久久一区二区三区| 黄色毛片三级朝国网站| 少妇被粗大的猛进出69影院 | 永久免费av网站大全| 中文字幕最新亚洲高清| 久久精品国产a三级三级三级| 欧美三级亚洲精品| av国产精品久久久久影院| 高清在线视频一区二区三区| 成人手机av| 丝袜美足系列| 日韩成人伦理影院| 国产男人的电影天堂91| 自线自在国产av| 天天操日日干夜夜撸| 美女xxoo啪啪120秒动态图| 国产欧美日韩一区二区三区在线 | 一级黄片播放器| 亚洲中文av在线| 欧美精品国产亚洲| 狠狠精品人妻久久久久久综合| 国产有黄有色有爽视频| 天堂8中文在线网| 一级片'在线观看视频| 日韩成人av中文字幕在线观看| 国产日韩欧美亚洲二区| 黄片无遮挡物在线观看| 欧美日韩av久久| 精品少妇黑人巨大在线播放| 中国三级夫妇交换| 久久久久久久大尺度免费视频| 制服丝袜香蕉在线| 亚洲精品av麻豆狂野| 国产精品麻豆人妻色哟哟久久| 日韩av不卡免费在线播放| 亚洲色图综合在线观看| 国产av一区二区精品久久| 欧美日韩视频精品一区| 曰老女人黄片| 九九在线视频观看精品| 18禁在线播放成人免费| 黄色欧美视频在线观看| 精品久久久久久久久亚洲| 国产亚洲精品久久久com| 亚洲性久久影院| 亚洲综合色惰| 女的被弄到高潮叫床怎么办| 大话2 男鬼变身卡| 亚洲色图综合在线观看| 亚洲精华国产精华液的使用体验| 欧美人与性动交α欧美精品济南到 | 97超碰精品成人国产| 成人综合一区亚洲| 成人免费观看视频高清| 亚洲精品日韩在线中文字幕| 国产男人的电影天堂91| 免费观看的影片在线观看| 老女人水多毛片| 精品熟女少妇av免费看| 色吧在线观看| 亚洲熟女精品中文字幕| 久久久久久久久久久免费av| 亚洲国产精品一区二区三区在线| 欧美精品一区二区大全| 精品久久久久久电影网| 日日摸夜夜添夜夜爱| 国产精品久久久久久精品古装| 久久婷婷青草| 我的老师免费观看完整版| 日日撸夜夜添| 日韩免费高清中文字幕av| 亚洲人成网站在线播| 大香蕉久久网| 成人亚洲精品一区在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品视频女| videossex国产| 欧美激情 高清一区二区三区| 久久久久视频综合| 国产精品久久久久久久久免| 纵有疾风起免费观看全集完整版| 国产成人免费无遮挡视频| av视频免费观看在线观看| 色婷婷av一区二区三区视频| 亚洲精品国产av蜜桃| 久久免费观看电影| 国产亚洲精品久久久com| 人妻制服诱惑在线中文字幕| 中文字幕最新亚洲高清| 亚洲精品久久成人aⅴ小说 | 999精品在线视频| 亚洲欧美精品自产自拍| videossex国产| 在线精品无人区一区二区三| 欧美另类一区| 在线观看免费视频网站a站| 亚洲国产av新网站| 中文字幕亚洲精品专区| 在线观看免费日韩欧美大片 | 18在线观看网站| 久久国内精品自在自线图片| 亚洲无线观看免费| 亚洲欧洲精品一区二区精品久久久 | 中国美白少妇内射xxxbb| 18禁裸乳无遮挡动漫免费视频| 999精品在线视频| 久久久久久久久久久久大奶| 五月天丁香电影| 久久国产精品大桥未久av| 国产日韩一区二区三区精品不卡 | 搡女人真爽免费视频火全软件| 人妻系列 视频| 午夜福利在线观看免费完整高清在| 91精品国产国语对白视频| 欧美激情 高清一区二区三区| 国产亚洲欧美精品永久| 少妇精品久久久久久久| 精品国产乱码久久久久久小说| 一级黄片播放器| 亚洲欧美精品自产自拍| 一本大道久久a久久精品| 伦理电影大哥的女人| 制服人妻中文乱码| 欧美一级a爱片免费观看看| 国产精品秋霞免费鲁丝片| 麻豆精品久久久久久蜜桃| 亚洲精品aⅴ在线观看| 性色avwww在线观看| 免费大片黄手机在线观看| av一本久久久久| 一个人免费看片子| 日韩大片免费观看网站| 亚洲国产精品国产精品| 亚洲美女视频黄频| 国内精品宾馆在线| 亚洲av在线观看美女高潮| 亚洲精品第二区| 在线播放无遮挡| 午夜av观看不卡| 女性被躁到高潮视频| 久久国产精品大桥未久av| av国产精品久久久久影院| 欧美激情 高清一区二区三区| 国产精品免费大片| 日本91视频免费播放| 亚洲美女视频黄频| 国产深夜福利视频在线观看| 久久久久精品久久久久真实原创| 久热久热在线精品观看| 国产精品无大码| 免费黄色在线免费观看| 欧美精品人与动牲交sv欧美| 亚洲国产精品成人久久小说| 天天操日日干夜夜撸| 我要看黄色一级片免费的| 免费av不卡在线播放| 高清不卡的av网站| 成年人免费黄色播放视频| 国产视频首页在线观看| 亚洲精品国产av蜜桃| 少妇高潮的动态图| 亚洲av在线观看美女高潮| 高清黄色对白视频在线免费看| 777米奇影视久久| 少妇 在线观看| 岛国毛片在线播放| 国产成人免费观看mmmm| 国产精品99久久99久久久不卡 | 99国产精品免费福利视频| 国产欧美日韩综合在线一区二区| 久久精品久久精品一区二区三区| 欧美日韩亚洲高清精品| 国产高清不卡午夜福利| 免费播放大片免费观看视频在线观看| av国产久精品久网站免费入址| 丝袜喷水一区| 国产精品 国内视频| 久久青草综合色| 亚洲av成人精品一区久久| av播播在线观看一区| 日韩欧美精品免费久久| 国产又色又爽无遮挡免| www.色视频.com| 亚洲欧美一区二区三区国产| 黄色一级大片看看| 精品久久久久久久久av| 欧美 日韩 精品 国产| 满18在线观看网站| 天美传媒精品一区二区| 国产一区有黄有色的免费视频| 如日韩欧美国产精品一区二区三区 | 亚洲精品色激情综合| 亚洲第一区二区三区不卡| 成人影院久久| 伊人久久国产一区二区| 三级国产精品欧美在线观看| 午夜激情久久久久久久| 蜜桃国产av成人99| 日本-黄色视频高清免费观看| 免费av中文字幕在线| 99热这里只有精品一区| 三级国产精品欧美在线观看| 久久99热这里只频精品6学生| 亚洲av成人精品一区久久| 制服诱惑二区| 简卡轻食公司| 国产成人精品在线电影| 搡老乐熟女国产| 97在线人人人人妻| 免费人妻精品一区二区三区视频| 一本大道久久a久久精品| 亚洲久久久国产精品| xxx大片免费视频| 99热这里只有精品一区| www.色视频.com| 乱码一卡2卡4卡精品| 色婷婷av一区二区三区视频| 久久久久久久亚洲中文字幕| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 少妇的逼好多水| 久久狼人影院| 免费av中文字幕在线| 赤兔流量卡办理| 国产高清国产精品国产三级| 十八禁网站网址无遮挡| 99九九线精品视频在线观看视频| 97超碰精品成人国产| 欧美性感艳星| 黄色一级大片看看| 日韩欧美精品免费久久| 三级国产精品欧美在线观看| 人妻少妇偷人精品九色| 午夜日本视频在线| 婷婷色av中文字幕| 少妇丰满av| 亚洲不卡免费看| 亚洲伊人久久精品综合| tube8黄色片| 一级片'在线观看视频| 国产免费又黄又爽又色| √禁漫天堂资源中文www| 国产成人午夜福利电影在线观看| 一级二级三级毛片免费看| 伦精品一区二区三区| 精品少妇久久久久久888优播| 岛国毛片在线播放| 黄片播放在线免费| 成年av动漫网址| 国产爽快片一区二区三区| 免费高清在线观看日韩| 日本猛色少妇xxxxx猛交久久| 啦啦啦啦在线视频资源| 蜜桃在线观看..| 精品国产乱码久久久久久小说| 亚洲av免费高清在线观看| 国产白丝娇喘喷水9色精品| 91精品国产国语对白视频| 国产免费现黄频在线看| 国产免费一区二区三区四区乱码| 男女免费视频国产| 久久精品夜色国产| 欧美人与性动交α欧美精品济南到 | 日本黄色日本黄色录像| 日韩电影二区| 亚洲欧洲日产国产| 国产毛片在线视频| 亚洲欧美日韩卡通动漫| 高清午夜精品一区二区三区| 熟女人妻精品中文字幕| 午夜av观看不卡| 亚洲不卡免费看| 97超视频在线观看视频| 日韩av在线免费看完整版不卡| 下体分泌物呈黄色| 国产黄色免费在线视频| 纯流量卡能插随身wifi吗| 久久久久久伊人网av| 大又大粗又爽又黄少妇毛片口| 黑人猛操日本美女一级片| 欧美+日韩+精品| 丝袜美足系列| 91久久精品电影网| 国产一级毛片在线| 日本vs欧美在线观看视频| 国产一区有黄有色的免费视频| 国产亚洲午夜精品一区二区久久| 欧美人与善性xxx| 国产精品久久久久久av不卡| 国产免费一级a男人的天堂| 日韩一区二区三区影片| 国产一区二区在线观看日韩| av在线老鸭窝| 久久午夜综合久久蜜桃| 久久久久久久久久久免费av| 中文乱码字字幕精品一区二区三区| 国产极品天堂在线| 中文字幕亚洲精品专区| 一本色道久久久久久精品综合| 亚洲丝袜综合中文字幕| 国产黄片视频在线免费观看| 九九在线视频观看精品| 久久精品夜色国产| 国产极品粉嫩免费观看在线 | 欧美精品一区二区免费开放| 高清在线视频一区二区三区| 男女边吃奶边做爰视频| 黄色配什么色好看| 免费观看的影片在线观看| 亚洲av国产av综合av卡| 亚洲精品乱码久久久v下载方式| 在线观看三级黄色| 欧美人与善性xxx| 亚洲国产最新在线播放| 美女大奶头黄色视频| 人妻少妇偷人精品九色| av网站免费在线观看视频| 国精品久久久久久国模美| 高清欧美精品videossex| 欧美激情极品国产一区二区三区 | 成人国产麻豆网| 秋霞在线观看毛片| 精品人妻偷拍中文字幕| 日本黄色片子视频| 美女国产视频在线观看| 热99久久久久精品小说推荐| 午夜福利视频精品| 91久久精品电影网| 久久久久久久大尺度免费视频| 各种免费的搞黄视频| 少妇的逼水好多| 国产免费视频播放在线视频| 交换朋友夫妻互换小说| 亚洲国产成人一精品久久久| 一级片'在线观看视频| 91精品三级在线观看| 亚洲精品国产av成人精品| 2022亚洲国产成人精品| 日日啪夜夜爽| 久久久久人妻精品一区果冻| 亚洲成人av在线免费| 简卡轻食公司| 精品国产乱码久久久久久小说| 精品久久久噜噜| 中文字幕最新亚洲高清| 看免费成人av毛片| 国产成人免费无遮挡视频| 五月开心婷婷网| 国产免费现黄频在线看| 国产精品国产三级国产av玫瑰| 久久久久久久久久久丰满| 精品国产一区二区久久| 亚洲人成网站在线播| 国产亚洲欧美精品永久| 黄色配什么色好看| 国产精品久久久久久精品电影小说| 日本色播在线视频| 青春草视频在线免费观看| 最近中文字幕高清免费大全6| 制服丝袜香蕉在线| 香蕉精品网在线| 黑人高潮一二区| 亚洲久久久国产精品| 99久久精品国产国产毛片| 91aial.com中文字幕在线观看| av在线老鸭窝| 日本av手机在线免费观看| 夜夜看夜夜爽夜夜摸| 欧美丝袜亚洲另类| 成人国语在线视频| 伦精品一区二区三区| 99久久精品国产国产毛片| 国产午夜精品久久久久久一区二区三区| 精品亚洲成a人片在线观看| 中国国产av一级| 免费高清在线观看日韩| 国产乱来视频区| 色婷婷av一区二区三区视频| 日韩成人伦理影院| 国产精品久久久久成人av| 欧美日韩视频精品一区| 久久狼人影院| 国产成人a∨麻豆精品| 亚洲精品乱码久久久久久按摩| 内地一区二区视频在线| 国产视频首页在线观看| 熟女电影av网| 婷婷色av中文字幕| 午夜激情久久久久久久| 日韩熟女老妇一区二区性免费视频| 下体分泌物呈黄色| 少妇人妻 视频| 久久精品熟女亚洲av麻豆精品| 美女福利国产在线| 2022亚洲国产成人精品| 亚洲国产毛片av蜜桃av| 日韩成人av中文字幕在线观看| 国产免费一级a男人的天堂| 18禁观看日本| 久久ye,这里只有精品| 夜夜看夜夜爽夜夜摸| 久热久热在线精品观看| 国产精品三级大全| 大陆偷拍与自拍| 妹子高潮喷水视频| 少妇的逼好多水| 欧美最新免费一区二区三区| 韩国高清视频一区二区三区| 我的女老师完整版在线观看| 精品少妇黑人巨大在线播放| 亚洲欧美中文字幕日韩二区| 亚洲欧美成人综合另类久久久| 国产精品不卡视频一区二区| 热99国产精品久久久久久7| 最黄视频免费看| 另类亚洲欧美激情| 高清欧美精品videossex| 精品久久久久久久久亚洲| av视频免费观看在线观看| 国产在视频线精品| 最黄视频免费看| 热99久久久久精品小说推荐| 免费观看的影片在线观看| 精品国产乱码久久久久久小说| 热99久久久久精品小说推荐| 免费观看的影片在线观看| 久久精品久久精品一区二区三区| 日本-黄色视频高清免费观看| 亚洲精品国产av蜜桃| 久久ye,这里只有精品| 熟女电影av网| 纯流量卡能插随身wifi吗| 成人影院久久|