• <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
    偷拍熟女少妇极品色| 88av欧美| 88av欧美| 精品久久久久久久毛片微露脸| 19禁男女啪啪无遮挡网站| 久久香蕉国产精品| 哪里可以看免费的av片| 久久精品综合一区二区三区| 久久久久久大精品| 一夜夜www| 久久天躁狠狠躁夜夜2o2o| 免费看十八禁软件| 久久久久九九精品影院| 99在线视频只有这里精品首页| 中文字幕精品亚洲无线码一区| 香蕉国产在线看| 别揉我奶头~嗯~啊~动态视频| av女优亚洲男人天堂 | 日韩欧美 国产精品| 亚洲一区二区三区色噜噜| 老熟妇仑乱视频hdxx| 男人舔奶头视频| 国产精品av视频在线免费观看| 免费av毛片视频| 国产黄a三级三级三级人| 我要搜黄色片| 久久久久国内视频| 女警被强在线播放| 免费在线观看成人毛片| 老司机福利观看| 日韩欧美国产在线观看| 看免费av毛片| 一级毛片精品| 好看av亚洲va欧美ⅴa在| 国产成人精品久久二区二区免费| 欧美3d第一页| 亚洲avbb在线观看| bbb黄色大片| 亚洲欧美一区二区三区黑人| 日本 av在线| 日韩欧美一区二区三区在线观看| 亚洲精品美女久久av网站| av在线天堂中文字幕| 亚洲国产精品成人综合色| 久久午夜亚洲精品久久| 一级作爱视频免费观看| 国产单亲对白刺激| 人妻丰满熟妇av一区二区三区| 国产亚洲精品av在线| 白带黄色成豆腐渣| 两个人的视频大全免费| 成人国产一区最新在线观看| 丁香欧美五月| 香蕉av资源在线| 亚洲片人在线观看| 男女做爰动态图高潮gif福利片| 草草在线视频免费看| 久久天堂一区二区三区四区| 国产亚洲欧美在线一区二区| 欧美成人性av电影在线观看| а√天堂www在线а√下载| 桃红色精品国产亚洲av| 可以在线观看的亚洲视频| 免费在线观看亚洲国产| 日本撒尿小便嘘嘘汇集6| 美女扒开内裤让男人捅视频| 久久久久久久久免费视频了| 国产午夜精品论理片| 三级男女做爰猛烈吃奶摸视频| 一区二区三区高清视频在线| 美女午夜性视频免费| 在线免费观看的www视频| 欧美一区二区国产精品久久精品| 成人午夜高清在线视频| 精品一区二区三区视频在线 | av在线蜜桃| 日韩免费av在线播放| 亚洲国产精品999在线| 美女免费视频网站| 久久精品国产清高在天天线| 国产私拍福利视频在线观看| 长腿黑丝高跟| 国产黄色小视频在线观看| 亚洲性夜色夜夜综合| 91久久精品国产一区二区成人 | 国产精品九九99| 两个人视频免费观看高清| 男人和女人高潮做爰伦理| 特大巨黑吊av在线直播| 国产在线精品亚洲第一网站| 黑人操中国人逼视频| h日本视频在线播放| 人妻夜夜爽99麻豆av| 日韩欧美精品v在线| 美女黄网站色视频| 精品国产乱码久久久久久男人| 18禁裸乳无遮挡免费网站照片| 久99久视频精品免费| 日韩欧美在线二视频| 看片在线看免费视频| 欧美日韩精品网址| 亚洲成人久久性| 99re在线观看精品视频| 长腿黑丝高跟| 99久久无色码亚洲精品果冻| 午夜精品久久久久久毛片777| 欧美极品一区二区三区四区| 一级毛片女人18水好多| 曰老女人黄片| 午夜免费成人在线视频| 18禁国产床啪视频网站| 欧美日韩乱码在线| 亚洲五月婷婷丁香| 老司机午夜福利在线观看视频| 真实男女啪啪啪动态图| 亚洲av成人一区二区三| 国产亚洲精品久久久久久毛片| 午夜免费成人在线视频| 99热只有精品国产| 我要搜黄色片| 欧美色视频一区免费| 日韩欧美在线二视频| 色综合婷婷激情| 日本成人三级电影网站| 欧美一区二区精品小视频在线| 精品国产超薄肉色丝袜足j| 国产精品国产高清国产av| 中文字幕高清在线视频| 怎么达到女性高潮| 日本在线视频免费播放| 国产单亲对白刺激| 操出白浆在线播放| 久久久久国内视频| 日韩av在线大香蕉| 最近视频中文字幕2019在线8| 一二三四社区在线视频社区8| 色噜噜av男人的天堂激情| 999久久久国产精品视频| 国内精品美女久久久久久| 禁无遮挡网站| 午夜精品在线福利| 精品熟女少妇八av免费久了| 久久这里只有精品19| 精品午夜福利视频在线观看一区| 亚洲在线观看片| 亚洲精品一卡2卡三卡4卡5卡| 色在线成人网| АⅤ资源中文在线天堂| 欧美激情在线99| 国产日本99.免费观看| 国产高清三级在线| 午夜视频精品福利| 中文字幕人成人乱码亚洲影| 久久精品影院6| 久久中文字幕一级| 别揉我奶头~嗯~啊~动态视频| 国产人伦9x9x在线观看| 超碰成人久久| 午夜免费成人在线视频| 亚洲精品粉嫩美女一区| 色精品久久人妻99蜜桃| 最新在线观看一区二区三区| 欧美最黄视频在线播放免费| 老司机在亚洲福利影院| 午夜免费成人在线视频| 波多野结衣高清作品| 18美女黄网站色大片免费观看| 中出人妻视频一区二区| 久久久久免费精品人妻一区二区| 国产99白浆流出| 欧美午夜高清在线| 欧美zozozo另类| 免费在线观看成人毛片| 啦啦啦观看免费观看视频高清| 国产精品国产高清国产av| 女人高潮潮喷娇喘18禁视频| 美女黄网站色视频| 中国美女看黄片| 老司机福利观看| 国产成+人综合+亚洲专区| 午夜福利在线观看免费完整高清在 | 免费在线观看影片大全网站| 一区福利在线观看| 黄频高清免费视频| 非洲黑人性xxxx精品又粗又长| 啦啦啦韩国在线观看视频| 精品不卡国产一区二区三区| 亚洲自拍偷在线| 亚洲无线在线观看| 亚洲成av人片免费观看| 日本在线视频免费播放| 成年版毛片免费区| 午夜福利欧美成人| 免费大片18禁| 桃色一区二区三区在线观看| 国产aⅴ精品一区二区三区波| 美女午夜性视频免费| 国产不卡一卡二| 搡老妇女老女人老熟妇| 两人在一起打扑克的视频| 日日夜夜操网爽| 欧美三级亚洲精品| 色综合亚洲欧美另类图片| 国产精品久久久久久人妻精品电影| 国产日本99.免费观看| 少妇的逼水好多| 久久精品人妻少妇| 中文字幕人妻丝袜一区二区| 成人无遮挡网站| 日韩欧美在线二视频| 色老头精品视频在线观看| 精华霜和精华液先用哪个| 19禁男女啪啪无遮挡网站| 1024手机看黄色片| 九色国产91popny在线| 在线a可以看的网站| 国产精品久久久人人做人人爽| 91久久精品国产一区二区成人 | 久久午夜亚洲精品久久| 久久久精品大字幕| 国产精品久久久久久久电影 | 琪琪午夜伦伦电影理论片6080| 国产淫片久久久久久久久 | 国产精品综合久久久久久久免费| 丁香欧美五月| 久久精品夜夜夜夜夜久久蜜豆| 久久久国产精品麻豆| 午夜精品久久久久久毛片777| 级片在线观看| 亚洲av五月六月丁香网| 男人舔女人下体高潮全视频| 亚洲欧美日韩卡通动漫| 亚洲一区高清亚洲精品| 无限看片的www在线观看| 97人妻精品一区二区三区麻豆| 免费大片18禁| 精品无人区乱码1区二区| 91在线观看av| 亚洲成人精品中文字幕电影| 中文字幕熟女人妻在线| 狂野欧美白嫩少妇大欣赏| 色综合婷婷激情| 女人高潮潮喷娇喘18禁视频| av在线天堂中文字幕| 首页视频小说图片口味搜索| 久久久久国产精品人妻aⅴ院| 精华霜和精华液先用哪个| 免费人成视频x8x8入口观看| 成人亚洲精品av一区二区| 曰老女人黄片| 在线观看一区二区三区| 男女视频在线观看网站免费| 亚洲精华国产精华精| 国产精品乱码一区二三区的特点| 欧美日本亚洲视频在线播放| 极品教师在线免费播放| 色视频www国产| 国产精品国产高清国产av| 国产伦精品一区二区三区视频9 | 国产高清激情床上av| 悠悠久久av| 国内揄拍国产精品人妻在线| 夜夜爽天天搞| 亚洲avbb在线观看| 在线十欧美十亚洲十日本专区| 国产一区二区在线观看日韩 | 波多野结衣高清作品| 操出白浆在线播放| 中文字幕高清在线视频| 国产乱人视频| 亚洲性夜色夜夜综合| 色视频www国产| 日本熟妇午夜| 久久精品aⅴ一区二区三区四区| 日本 欧美在线| 成人18禁在线播放| 久久国产精品影院| 宅男免费午夜| 久久久国产成人免费| 久久香蕉国产精品| av在线天堂中文字幕| 欧美中文日本在线观看视频| 亚洲一区二区三区不卡视频| 国产极品精品免费视频能看的| 少妇人妻一区二区三区视频| 一个人免费在线观看电影 | 精品福利观看| 美女免费视频网站| 色综合站精品国产| 久99久视频精品免费| 真人做人爱边吃奶动态| 九九在线视频观看精品| 18禁国产床啪视频网站| 成人av在线播放网站| 亚洲在线自拍视频| 老司机在亚洲福利影院| 88av欧美| 亚洲五月天丁香| av黄色大香蕉| 三级国产精品欧美在线观看 | 欧美中文日本在线观看视频| 欧美日韩精品网址| 99久久精品热视频| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 禁无遮挡网站| 亚洲一区高清亚洲精品| 亚洲av五月六月丁香网| 黄色女人牲交| 久久久久久大精品| 欧美黄色淫秽网站| 欧美日韩综合久久久久久 | av黄色大香蕉| 午夜福利高清视频| 亚洲狠狠婷婷综合久久图片| 女人被狂操c到高潮| 免费观看精品视频网站| 伊人久久大香线蕉亚洲五| 国内揄拍国产精品人妻在线| 哪里可以看免费的av片| 夜夜爽天天搞| xxx96com| 久久精品人妻少妇| 亚洲av日韩精品久久久久久密| 老司机午夜十八禁免费视频| 看免费av毛片| 人人妻人人看人人澡| 久久久国产欧美日韩av| 日韩欧美免费精品| 色综合欧美亚洲国产小说| 久久久久国产精品人妻aⅴ院| 香蕉av资源在线| 欧美性猛交╳xxx乱大交人| 久久久久国产一级毛片高清牌| 91九色精品人成在线观看| 欧美日韩精品网址| 五月玫瑰六月丁香| 久久精品国产99精品国产亚洲性色| 日本三级黄在线观看| 国产一级毛片七仙女欲春2| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 又黄又粗又硬又大视频| 村上凉子中文字幕在线| av在线天堂中文字幕| 亚洲avbb在线观看| 免费在线观看日本一区| 午夜免费成人在线视频| 90打野战视频偷拍视频| 久久这里只有精品19| 不卡av一区二区三区| 国内揄拍国产精品人妻在线| 国产精品久久久av美女十八| 桃红色精品国产亚洲av| 色尼玛亚洲综合影院| 91字幕亚洲| 18禁黄网站禁片免费观看直播| 日本一二三区视频观看| 久久伊人香网站| 特级一级黄色大片| av国产免费在线观看| 极品教师在线免费播放| 亚洲人与动物交配视频| 欧美最黄视频在线播放免费| 久久人人精品亚洲av| 99久久精品热视频| 日韩有码中文字幕| 黑人欧美特级aaaaaa片| 人人妻人人澡欧美一区二区| 亚洲狠狠婷婷综合久久图片| 久久伊人香网站| 91老司机精品| 一级a爱片免费观看的视频| 国产精品久久久久久人妻精品电影| svipshipincom国产片| 日本黄大片高清| 极品教师在线免费播放| 99久久精品国产亚洲精品| 免费看光身美女| 美女免费视频网站| 97超视频在线观看视频| 操出白浆在线播放| 成人精品一区二区免费| 国产亚洲精品久久久久久毛片| 成人特级av手机在线观看| 日本在线视频免费播放| 男女之事视频高清在线观看| 日本 欧美在线| 国产1区2区3区精品| 国产又色又爽无遮挡免费看| 午夜亚洲福利在线播放| 久久国产精品人妻蜜桃| 成人欧美大片| 热99在线观看视频| 午夜精品久久久久久毛片777| 亚洲av片天天在线观看| 免费av毛片视频| 免费在线观看亚洲国产| 国产在线精品亚洲第一网站| 狂野欧美白嫩少妇大欣赏| 国产激情久久老熟女| 国产高潮美女av| 69av精品久久久久久| 中文字幕久久专区| 日本 av在线| 一进一出抽搐动态| h日本视频在线播放| 制服丝袜大香蕉在线| 性欧美人与动物交配| 1000部很黄的大片| 久久亚洲真实| 精品久久久久久久久久久久久| 两个人的视频大全免费| 欧美成人性av电影在线观看| 久久午夜综合久久蜜桃| 日韩成人在线观看一区二区三区| 国产亚洲精品久久久com| 欧美高清成人免费视频www| 免费av毛片视频| 97超级碰碰碰精品色视频在线观看| 在线a可以看的网站| 我要搜黄色片| 怎么达到女性高潮| 国产精品av视频在线免费观看| 18禁黄网站禁片午夜丰满| 老熟妇乱子伦视频在线观看| 久久久国产精品麻豆| 亚洲欧洲精品一区二区精品久久久| 日本五十路高清| 一本久久中文字幕| 久久久久久久久免费视频了| 国产野战对白在线观看| 午夜精品久久久久久毛片777| 一级黄色大片毛片| 舔av片在线| 国产单亲对白刺激| 在线a可以看的网站| 国产一区二区激情短视频| 欧美日韩亚洲国产一区二区在线观看| 变态另类成人亚洲欧美熟女| 亚洲狠狠婷婷综合久久图片| 色综合亚洲欧美另类图片| 亚洲人成伊人成综合网2020| 99精品久久久久人妻精品| 国产精品99久久99久久久不卡| 国产高潮美女av| 一夜夜www| 精品一区二区三区av网在线观看| 午夜福利在线在线| 搡老岳熟女国产| 国产成人一区二区三区免费视频网站| 亚洲一区二区三区色噜噜| 在线免费观看的www视频| 一卡2卡三卡四卡精品乱码亚洲| 国产三级在线视频| 视频区欧美日本亚洲| 日本一本二区三区精品| 日韩免费av在线播放| 香蕉国产在线看| 中文资源天堂在线| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 亚洲精品中文字幕一二三四区| 99热6这里只有精品| 国内精品久久久久精免费| 久久久成人免费电影| 又爽又黄无遮挡网站| 国产一区二区三区视频了| 国产av不卡久久| 夜夜躁狠狠躁天天躁| 婷婷精品国产亚洲av在线| 五月伊人婷婷丁香| 国产日本99.免费观看| av在线蜜桃| 久久亚洲真实| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 熟女少妇亚洲综合色aaa.| 99国产极品粉嫩在线观看| 国产成年人精品一区二区| 欧美国产日韩亚洲一区| ponron亚洲| a在线观看视频网站| 99热这里只有是精品50| 亚洲片人在线观看| 女人被狂操c到高潮| 日本五十路高清| 欧美激情在线99| 欧美日韩国产亚洲二区| 天天添夜夜摸| 99热这里只有精品一区 | 国产精品亚洲美女久久久| 噜噜噜噜噜久久久久久91| 日日干狠狠操夜夜爽| 欧美日韩国产亚洲二区| 亚洲精品美女久久久久99蜜臀| 99久久综合精品五月天人人| 国产精品国产高清国产av| 日韩欧美免费精品| 免费看光身美女| 国产男靠女视频免费网站| 无人区码免费观看不卡| 一进一出抽搐gif免费好疼| 在线观看午夜福利视频| 免费无遮挡裸体视频| 国产高清视频在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 大型黄色视频在线免费观看| 精品一区二区三区四区五区乱码| 99视频精品全部免费 在线 | 看黄色毛片网站| 亚洲熟妇中文字幕五十中出| 天堂影院成人在线观看| 在线观看一区二区三区| 亚洲人成电影免费在线| 999精品在线视频| 国产精品爽爽va在线观看网站| 国产亚洲精品综合一区在线观看| www.999成人在线观看| 欧美丝袜亚洲另类 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品999在线| 亚洲精品久久国产高清桃花| 日韩中文字幕欧美一区二区| 18禁美女被吸乳视频| 熟女少妇亚洲综合色aaa.| 美女免费视频网站| 午夜免费激情av| 欧美日韩亚洲国产一区二区在线观看| 日韩精品青青久久久久久| 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 国产精品女同一区二区软件 | 国产欧美日韩精品一区二区| 国产91精品成人一区二区三区| 久久伊人香网站| 男女做爰动态图高潮gif福利片| 午夜福利在线观看吧| 欧美高清成人免费视频www| a级毛片a级免费在线| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利欧美成人| 午夜成年电影在线免费观看| 一区福利在线观看| 亚洲天堂国产精品一区在线| 亚洲午夜理论影院| 久久中文字幕一级| 国内精品一区二区在线观看| 国产亚洲欧美在线一区二区| 国产高清三级在线| 亚洲,欧美精品.| 亚洲天堂国产精品一区在线| 91麻豆精品激情在线观看国产| 夜夜爽天天搞| 久久国产精品影院| 欧美乱妇无乱码| 国产一级毛片七仙女欲春2| 给我免费播放毛片高清在线观看| 69av精品久久久久久| 九九在线视频观看精品| 国产精华一区二区三区| 在线观看66精品国产| 欧美xxxx黑人xx丫x性爽| 好看av亚洲va欧美ⅴa在| 51午夜福利影视在线观看| 亚洲色图av天堂| 禁无遮挡网站| 亚洲一区二区三区不卡视频| 日本a在线网址| 精品福利观看| 老鸭窝网址在线观看| 欧美三级亚洲精品| 国产精品免费一区二区三区在线| 亚洲 欧美一区二区三区| 国产私拍福利视频在线观看| 国产激情久久老熟女| 97人妻精品一区二区三区麻豆| 99久久国产精品久久久| 午夜福利18| www.精华液| 最新中文字幕久久久久 | 国产精品一区二区精品视频观看| 亚洲国产精品成人综合色| 在线观看66精品国产| 19禁男女啪啪无遮挡网站| 国产欧美日韩精品一区二区| 免费大片18禁| 婷婷亚洲欧美| 午夜福利在线观看吧| 午夜福利在线观看免费完整高清在 | 国产精品av视频在线免费观看| tocl精华| 波多野结衣高清无吗| 日韩有码中文字幕| 三级国产精品欧美在线观看 | 亚洲av成人不卡在线观看播放网| 两人在一起打扑克的视频| 一级黄色大片毛片| 国产亚洲av嫩草精品影院| 人妻夜夜爽99麻豆av| 97超级碰碰碰精品色视频在线观看| 久久国产乱子伦精品免费另类| 亚洲成a人片在线一区二区| 最好的美女福利视频网| 欧美乱妇无乱码| 亚洲真实伦在线观看| 日韩人妻高清精品专区| 99久久综合精品五月天人人| 国产精品久久久久久精品电影| 国产探花在线观看一区二区| 舔av片在线| 一个人看的www免费观看视频| 两个人视频免费观看高清| 精品一区二区三区视频在线观看免费| 国产美女午夜福利|