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

    迎風(fēng)凹腔及其組合體減阻防熱技術(shù)研究進展

    2021-09-18 06:19:42張杰肖鋒黃偉顏力孟玉珊
    航空兵器 2021年4期

    張杰 肖鋒 黃偉 顏力 孟玉珊

    摘 要: 熱防護技術(shù)是高超聲速飛行器面向工程實際應(yīng)用時必須要解決的關(guān)鍵技術(shù)之一。 為實現(xiàn)更有效的減阻防熱, 學(xué)者們提出了多種熱防護方案, 其中迎風(fēng)凹腔及其組合方式是最有效、 最有發(fā)展前景的主動熱防護方案之一。 本文介紹了單一迎風(fēng)凹腔及其與逆向射流、 能量沉積、 發(fā)散冷卻等組合方案的作用機理并總結(jié)歸納了其研究現(xiàn)狀。 最后提出了對未來的展望, 尤其是跟進并創(chuàng)新已有組合方案, 并進行更多地面試驗佐證數(shù)值模擬結(jié)果, 以期早日應(yīng)用于工程實際。

    關(guān)鍵詞:減阻; 防熱; 迎風(fēng)凹腔; 逆向射流; 能量沉積; 高超聲速飛行器; 熱防護技術(shù)

    中圖分類號:V211.3?? 文獻標(biāo)識碼:??? A?? 文章編號:1673-5048(2021)04-0016-08

    0 引? 言

    近年來, 高超聲速飛行器以其強突防能力、 快速反應(yīng)能力與高飛行速度等特點, 日益成為世界各航空航天大國竟相發(fā)展的目標(biāo)。 不論在軍事領(lǐng)域還是民用范疇, 高超聲速飛行器都極具發(fā)展?jié)摿蛻?yīng)用前景。 因其長時間高馬赫數(shù)飛行, 嚴酷的氣動熱環(huán)境成為高超聲速飛行器設(shè)計過程中一個無法回避的問題, 熱防護技術(shù)直接關(guān)系到飛行器的飛行安全, 也是決定高超聲速飛行器能否應(yīng)用于工程實際的關(guān)鍵[1-2]。

    飛行器在高超聲速狀態(tài)下作業(yè)時, 除受摩擦阻力外, 還會受到弓形激波壓縮產(chǎn)生的阻力, 可占全部阻力的50%以上, 而且該比重隨來流馬赫數(shù)的擴大持續(xù)增長, 這將嚴重影響飛行器的氣動性能, 加大燃料消耗, 減少有效載荷, 因此減小波阻的重要性不言而喻[3-4]。 由于氣體粘性阻滯作用, 飛行器壁面熱流激增使得壁面溫度升高, 尤其以頭部最為顯著, 高溫高壓的氣體通過機體外殼傳熱至內(nèi)部, 影響飛行器內(nèi)各類電子設(shè)備工作性能甚至導(dǎo)致其失效[5]。 Huang[6]給出了三種典型再入飛行器的熱流對比, 如圖1所示, 可以直觀看出峰值熱流約為500 kW/m2, 且大多數(shù)時間熱流強度不低于100 kW/m2。

    面對如此強大的氣動阻力以及嚴峻的氣動加熱現(xiàn)象, 設(shè)計出能有效保護飛行器表面的減阻防熱系統(tǒng)至關(guān)重要。 近年來, 國內(nèi)外學(xué)者提出許多主動冷卻方法如逆向射流[7]、 加裝減阻桿[8]、 迎風(fēng)凹腔[9]、 能量沉積[10]及許多組合方式來實現(xiàn)減阻防熱功能。 這類主動冷卻方法主要通過噴射冷卻工質(zhì)或者設(shè)計特定機械結(jié)構(gòu), 加裝減阻設(shè)備等方式改變鼻錐繞流流場以實現(xiàn)減阻防熱目的。 其中, Knight[11]對能量沉積方案減阻防熱的研究進展進行了全面總結(jié), 但此方法仍停留在理論階段, 工程應(yīng)用存在困難。 Ahmed等[12]對高超聲速飛行器加裝減阻桿實現(xiàn)減阻防熱方案的研究進展進行了系統(tǒng)總結(jié), 指出了減阻桿會破壞氣動外形等一些問題, 并相應(yīng)提出研究展望。 進而, Huang等[13]也分別對減阻桿及其組合體、 逆向射流及其組合體[6]在超/高超聲速氣流中的減阻防熱研究進展進行了比較系統(tǒng)的總結(jié)。

    迎風(fēng)凹腔構(gòu)型自提出以來, 學(xué)者們做了很多研究與試驗。 本文對迎風(fēng)凹腔及其組合體構(gòu)型的研究現(xiàn)狀進行總結(jié)歸納, 主要分為迎風(fēng)凹腔單一構(gòu)型和組合體兩部分, 組合體部分包括逆向射流、 能量沉積、 發(fā)散冷卻與迎風(fēng)凹腔的組合方式。 本文介紹了單一構(gòu)型和組合體的作用機理和國內(nèi)外研究進展, 包括數(shù)值模擬和地面試驗研究。 最后對主要組合方式的優(yōu)勢與研究中存在的問題進行了總結(jié), 并對未來的發(fā)展進行了展望。

    1 單一迎風(fēng)凹腔

    1.1 結(jié)構(gòu)來源及基本原理

    1922年, Hartmann[14]發(fā)現(xiàn)將迎風(fēng)凹腔置于高超聲速流場中, 會發(fā)出高強度、 不連續(xù)的聲音, 可以將其作為聲源裝置, 其目的并非應(yīng)用于高超聲速鈍頭體頭部的減阻防熱。 隨后, 將其應(yīng)用于熱防護方案的構(gòu)想在20世紀(jì)50年代由Burbank等[15]提出, 并成為最受關(guān)注的減阻防熱方案之一, 其結(jié)構(gòu)示意圖見圖2。 圖中直觀展示了在飛行器頭部設(shè)置的軸對稱迎風(fēng)凹腔, 與高超聲速來流正對, 與鼻錐同軸。 就腔體而言, 整體構(gòu)成比較簡單, 僅有一圓形底面及一圈側(cè)壁面, 側(cè)壁前端與鼻錐外壁面相接處為唇口。 因此, 凹腔底壁直徑(D)和側(cè)壁長度(L)是值得關(guān)注并開展研究的結(jié)構(gòu)參數(shù)變量。

    高超聲速來流在凹腔中產(chǎn)生震蕩, 來流的能量通過震蕩耗散從而降低飛行器的表面熱流。 陸海波等人曾具體分析其作用機理, 簡要總結(jié)如下: 由于凹腔的振蕩特性, 激波處于凹腔唇口上游的某個位置, 并不斷震動。 此時, 就會存在兩種情況: 一是激波靠近腔體, 二是激波遠離腔體。 第一種情況下, 弓形激波愈發(fā)靠近鼻錐, 激波壓縮加劇, 波后流體馬赫數(shù)降低, 溫度上升, 真實氣體定壓比熱上升, 在總焓值不變的條件下, 波后流體的滯止溫度下降, 其攜帶的能量也隨之減少, 從而減小了頭部的氣動加熱。 與之相反的第二種情況下, 雖然滯止溫度上升, 但由于激波遠離腔體, 下游氣體會補充到兩者之間, 起到兩個作用: 第一, 該部分氣體總溫比波后流體低, 起到一定的冷卻降溫作用, 減小唇口燒蝕程度; 第二, 導(dǎo)致流動分離形成回流區(qū)。 這兩點保護了鼻錐表面, 降低了鼻錐滯止區(qū)域的氣動加熱[17]。

    Johnson[18]研究了迎風(fēng)凹腔的流場分布, 發(fā)現(xiàn)該流場具有顯著的非定常性, 這一結(jié)論也在后人的試驗[19]中得到了驗證。 圖3所示的激波紋影圖像, 可以明顯觀察到激波的前后往復(fù)運動。

    Hartmann[14]通過試驗方法研究了迎風(fēng)凹腔流動的振蕩特性, 之后很多學(xué)者也通過數(shù)值方法對此特性開展了相關(guān)性研究, 探究其影響因素。 對于淺腔, 引發(fā)流動震蕩的是自由來流的噪聲, 且整體流動對來流擾動非常敏感, 但是對于高超聲速流場中的深腔, 流動自發(fā)振蕩。 其不但增加了流場預(yù)測的難度, 而且也會引發(fā)結(jié)構(gòu)震動從而導(dǎo)致控制困難, 這也是凹腔結(jié)構(gòu)的一個顯著缺點, 高速飛行下的操控靈敏性降低可能會埋下隱患。 這一減阻防熱方案會在凹腔唇口產(chǎn)生嚴重的燒蝕, 但這種燒蝕也遠低于原始駐點處的燒蝕程度, 其結(jié)構(gòu)簡單、 效果顯著, 引起了許多學(xué)者的關(guān)注與研究。 采用迎風(fēng)凹腔結(jié)構(gòu)可以獲得一定防熱效果, 但不一定能夠獲得減阻效果, 且迎風(fēng)凹腔的冷卻效果也并沒有達到單一結(jié)構(gòu)應(yīng)對極端加熱條件的程度, 因此需要研究更多的組合構(gòu)型。

    1.2 研究進展

    自迎風(fēng)凹腔構(gòu)型被提出以來, 國內(nèi)外學(xué)者對這一構(gòu)型開展了廣泛的試驗與數(shù)值研究, 尤其是20世紀(jì)末, 其試驗愈發(fā)增多。 Yuceil等[20]通過紅外攝像儀對比發(fā)現(xiàn), 相較于無腔體結(jié)構(gòu)的球形鼻錐, 有凹腔結(jié)構(gòu)的鼻錐在唇口的溫度更低, 并稱其為“冷環(huán)”。 而后又對球頭圓柱進行風(fēng)洞試驗[21-22], 將 L/D=0.4作為淺腔與深腔的分界線, L/D=0.7作為深腔與更深腔的分界線, 試驗發(fā)現(xiàn)當(dāng)長徑比值為2時整體流場穩(wěn)定且防熱效果最突出。 對于淺腔而言, 鼻區(qū)附近Δφ=5°~6°范圍內(nèi)存在“冷環(huán)”結(jié)構(gòu), 而且非常穩(wěn)定, 隨著長徑比增大至大于0.5時, “冷環(huán)”開始消失。 Engblom等[23-24]提出對尖唇口鈍化處理可消除原先唇口處的回流, 減少產(chǎn)熱,? 并且當(dāng)凹腔長徑比大于2時, 由于深腔的自振蕩效應(yīng), “冷環(huán)”現(xiàn)象重現(xiàn)且更為明顯。

    Saravanan等[19]對導(dǎo)彈外形試驗件展開了在來流條件為馬赫數(shù)8情況下的減阻防熱性能探究試驗, 試驗裝置示意如圖4所示。 研究中, 采用一個無量綱參數(shù)——斯坦頓數(shù)St沿壁面的分布來表征壁面熱流分布。 通過試驗對比, 迎風(fēng)凹腔構(gòu)型的壁面最大傳熱率可降低35%, 甚至更多。 在減阻方面, 正對來流情況下, 當(dāng)長徑比為4時, 阻力下降5.12%, 但當(dāng)該值為2時, 阻力卻增加了, 這證明迎風(fēng)凹腔可以有效防熱, 但減阻效果并不穩(wěn)定, 甚至?xí)幼琛?同時, 在試驗中還發(fā)現(xiàn)唇口壓強存在著類似正弦分布的震蕩現(xiàn)象, 這也驗證了流場的非定常性。 進一步的, 在后續(xù)的工作過程中, 使用非結(jié)構(gòu)網(wǎng)格在穩(wěn)態(tài)情況下對凹腔直徑6 mm和12 mm的模型進行數(shù)值模擬。 根據(jù)仿真結(jié)果顯示, 表面熱流分布和試驗吻合得較好, 但是壓強與阻力系數(shù)的結(jié)果差強人意。 文中還指出, 結(jié)合有限體積方法和有限差分方法, 可以更快速地仿真出全流場的準(zhǔn)確信息。

    陸海波等[25-26]根據(jù)前人的研究工作, 采用數(shù)值模擬方法繼續(xù)深入研究此迎風(fēng)凹腔構(gòu)型, 這也為之后研究逆向射流與迎風(fēng)凹腔組合體減阻防熱方案打下了堅實的基礎(chǔ)。 使用三維結(jié)構(gòu)網(wǎng)格, 引入k-ε湍流模型, 進行穩(wěn)態(tài)假設(shè), 對比了無凹腔、 凹腔半徑6 mm及凹腔半徑12 mm三種情況下的馬赫數(shù)云圖以及溫度云圖(見圖5), 從圖中也可以觀察到唇口的冷環(huán)現(xiàn)象。

    Ladoon等[27]繼續(xù)完善前人對振蕩特性的研究, 設(shè)計了可變深度的凹腔構(gòu)型, 借助激光對比研究淺腔與深腔的不同特性, 并給出穩(wěn)定振蕩下長徑比的取值。 該方法在后來學(xué)者的研究中得以復(fù)現(xiàn)[28], 證明長徑比在一定程度上可以表征振蕩阻尼, 是影響振蕩衰減率的重要參量。

    Silton等[29-30]用冰制作了大量直徑不同、? 長徑比不同的模型, 用以研究前緣凹腔減緩燒蝕的效應(yīng), 指出加長凹腔深度或者鈍化唇口半徑都可以延緩鼻區(qū)燒蝕, 并給出了一個最佳的構(gòu)型以供參考。 與此同時, 對唇口鈍化的模型進行了數(shù)值模擬, 結(jié)果吻合較好, 指出利用前緣鈍化方法, 可以達到材料的均勻燒蝕。

    很長一段時間內(nèi), 學(xué)者們專注于對工作參數(shù)的研究, 近些年有學(xué)者開始將目光轉(zhuǎn)向凹腔的結(jié)構(gòu)參數(shù), 對于傳統(tǒng)的迎風(fēng)凹腔構(gòu)型提出了創(chuàng)新。 Yadav等[31]在球頭圓柱的頭部加入了一個很淺的拋物線形迎風(fēng)凹腔, 如圖6所示。 數(shù)值模擬結(jié)果表明, 此拋物形凹腔越深, 整體減阻防熱效果越好, 但這個結(jié)果是建立在唇口加熱加劇的基礎(chǔ)上, 因此, 就設(shè)計更有效減阻防熱構(gòu)型的角度而言, 該嘗試不算成功, 但為學(xué)者們提供了一種新的思路, 極具參考價值。 之后, 孫喜萬[32]提出了一種基于近似最大推力噴管的新迎風(fēng)凹腔構(gòu)型, 并對其減阻防熱性能開展了數(shù)值模擬, 證明相比原先構(gòu)型, 新的迎風(fēng)凹腔構(gòu)型在減阻防熱方面的性能確有提升。

    總體而言, 迎風(fēng)凹腔是結(jié)構(gòu)最為簡單的主動熱防護方案, 不必額外加裝裝置, 不需要儲存冷卻工質(zhì), 也不需要設(shè)計噴流裝置或者發(fā)汗流道等, 只需改變前緣的結(jié)構(gòu)。 但凹腔流場的非定常性以及腔內(nèi)的震蕩是很明顯的缺點, 在應(yīng)用時必須考慮好通過何種方式穩(wěn)定前緣受力, 不能放大凹腔對操作靈敏性的影響。 同時, 對于凹腔的研究可以繼續(xù)考慮如何控制或者利用這種振蕩的特性, 以及如何處理在有攻角情況下, 凹腔對氣動力的影響。

    2 組合構(gòu)型

    目前的減阻防熱方案可以根據(jù)其是否可控分為兩大類, 即被動熱防護方法和主動減阻防熱方法[33]。 被動熱防護方法是通過耐熱材料或燒蝕材料達到防熱目的, 其發(fā)展主要依托新材料的研究[34]。 主動減阻防熱方法則是向飛行器前緣流場注入冷卻工質(zhì)或設(shè)計特定的結(jié)構(gòu), 以改變飛行器頭部繞流流場, 降低飛行器頭部壁面熱流達到熱防護目的, 具體可劃分為三類: 加注工質(zhì)的流體方法, 包括逆向射流、 發(fā)散冷卻等; 加裝裝置或改造頭部構(gòu)型的結(jié)構(gòu)方法, 如減阻桿、 迎風(fēng)凹腔等; 來流高能的能量沉積方法。 盡管單一原理的方案有巨大的應(yīng)用前景, 但是將多種機理結(jié)合起來, 會起到揚長避短的效果, 因此, 取主動減阻防熱方法中的兩類方法與迎風(fēng)凹腔結(jié)構(gòu)組合, 即可形成如下三種組合構(gòu)型。

    2.1 迎風(fēng)凹腔和逆向射流組合體

    逆向射流這種減阻防熱方案的提出旨在解決一些熱防護系統(tǒng)不能重復(fù)利用、 無法為高超聲速飛行器提供長期有效保護的問題[35]。 使用鈍頭體前端逆向噴流的方法可以顯著減小飛行器頭部的熱載荷和阻力。 圖7所示為迎風(fēng)凹腔與逆向射流組合構(gòu)型[36]。 就理論分析而言, 逆向射流可避免凹腔內(nèi)部流場振蕩, 凹腔反之也能起到降低流場穩(wěn)定臨界射流總壓, 提高射流效能的作用[16], 兩者可以很好互補。 組合構(gòu)型在實際運行時可分階段工作, 低馬赫時不開啟逆向射流, 僅使用凹腔結(jié)構(gòu), 而高馬赫時逆向射流也加入工作, 進一步提升減阻防熱的功效, 從而有效保護飛行器[37]。

    陸海波等[38]在研究不同參數(shù)對此構(gòu)型減阻防熱效能的影響之后,提出了迎風(fēng)凹腔和逆向射流組合減阻防熱方案, 圖8為數(shù)值模擬下的部分流場圖。? 當(dāng)凹腔長度足夠大時, 整個腔體類似一根加速噴管, 基于此腔體內(nèi)部

    構(gòu)型進行設(shè)計優(yōu)化, 起到減小激波損失加速射流, 提高效能的作用。 作者也進行了與之匹配的試驗, 結(jié)合對比試驗的紋影圖, 發(fā)現(xiàn)激波位置吻合較好, 整體結(jié)構(gòu)完整。

    Bibi等[39]研究了自由來流為馬赫數(shù)3.98條件下, 一種帶張角前置凹腔的幾何尺寸對半球體逆向射流減阻效果的影響, 發(fā)現(xiàn)使用帶張角前置凹腔的逆向射流比簡單的逆向射流可以獲得更好的減阻效果, 其研究的凹腔模型如圖9所示。

    Sun等[40]對圓柱形迎風(fēng)凹腔進行了改進, 提出了一種拋物線形凹腔, 研究在自由來流為馬赫數(shù)7.96條件下, 拋物線形凹腔與逆向射流結(jié)合的方法對鼻錐體的降熱效果, 得到采用拋物線形凹腔比圓柱形凹腔降熱效果更好的結(jié)論。 Sudarshan等[41]利用圖10所示的對比模型, 在相同試驗條件下也證實了這個結(jié)論的正確性。

    Zhang等[42]研究了自由來流為馬赫數(shù)3.98的條件下, 迎風(fēng)凹腔結(jié)合周期射流對鼻錐表面減阻防熱效果的瞬態(tài)特征, 得出結(jié)論: 迎風(fēng)凹腔的尺寸L/D=1.33時, 周期性射流可獲得最好減阻效果; 當(dāng)L/D=1.25 時, 則可獲得最好的降熱效果。

    總體而言, 該組合體可以顯著減小鼻錐體的阻力和表面熱流, 提高飛行器的控制性能, 同時使用小的噴流總壓有利于節(jié)省噴流能量, 提高噴流效率。 此外, 對于凹腔的結(jié)構(gòu)參數(shù)(如凹腔長度、 出口直徑、 初始膨脹半徑、 唇口鈍化半徑、 整體構(gòu)型等)與逆向射流的工作參數(shù)(如周期、 波形、 振幅等)改變, 一定程度上也會影響組合構(gòu)型的減阻防熱水平, 都可以進行進一步的數(shù)值模擬以及試驗研究。

    2.2 迎風(fēng)凹腔和能量沉積組合體

    相較于迎風(fēng)凹腔與逆向射流的組合構(gòu)型, 對能量沉積的研究就少很多, 其中一個原因是能量沉積方式在工程實現(xiàn)確實有難度, 但技術(shù)的實現(xiàn)需要契機, 這里簡單介紹下能量沉積[43]。

    如圖11所示, 點爆炸波造成的負壓區(qū)使弓形激波變形, 并引起滯止點附近橫向流動, 將滯止位置的空氣引向周圍, 形成低密度通道。 具體而言, 在飛行器上游點爆炸瞬間形成的高壓, 極熱氣體會突然膨脹, 從而通過爆炸波將周圍的介質(zhì)推開。 在介質(zhì)中這種局部“爆炸狀”行為的清除作用會形成一個非常熱、 低壓和低密度的氣泡區(qū)域, 當(dāng)被下游對流時, 氣泡會繼續(xù)膨脹[44]。 當(dāng)此波與高超聲飛行器之前的現(xiàn)有弓形激波相互作用時, 激波下游的高壓空氣在低密度氣泡內(nèi)膨脹。 另外, 氣泡內(nèi)部的高溫還具有將局部馬赫數(shù)降低到亞聲速值的作用。 最終, 駐點區(qū)周圍的壓力場達到較低穩(wěn)定的水平。 自上游到下游形成錐形或拋物形激波、 超聲速流、 弓形激波、 亞聲速流、 再附激波的結(jié)構(gòu)。 能量沉積的脈沖重復(fù)產(chǎn)生了更復(fù)雜的流場結(jié)構(gòu), 該結(jié)構(gòu)取決于能量水平、 重復(fù)頻率、 脈沖持續(xù)時間和沉積位置[45]。

    Bazyma等[46]研究了在自由來流為馬赫數(shù)3的條件下, 改變能量源的工作參數(shù), 包括能量密度及能量形狀等, 研究其對組合體方案的減阻防熱性能的影響。 根據(jù)結(jié)果顯示, 使用高頻的能量源甚至能夠抹去凹腔內(nèi)振蕩的變化, 將類似壓強的環(huán)境參數(shù)控制為一個定值, 較之于單一構(gòu)型, 組合方案似乎更為受控。 除此之外, 能量源的形狀也對減阻效能有影響, 其中球形的表現(xiàn)最佳, 如圖12所示。

    2.3 迎風(fēng)凹腔和發(fā)散冷卻組合體

    發(fā)散冷卻[47]是一種極具潛力的主動熱防護方式, 有很多優(yōu)點, 如冷卻能力強、 冷卻工質(zhì)用量小、 穩(wěn)定性高等, 但學(xué)者們通過數(shù)值模擬和試驗研究, 發(fā)現(xiàn)鼻錐模型表面極不均勻的壓力和溫度分布, 使得冷卻劑不易從駐點處流出, 從而達不到駐點處冷卻的目的[48]。 為提高駐點冷卻效率, 已有的直接優(yōu)化方法提升并不明顯, 因此, 在發(fā)散冷卻的基礎(chǔ)上尋求其他冷卻結(jié)構(gòu)的組合勢在必行。

    除常見組合構(gòu)型, 近些年也出現(xiàn)一些較為新穎的組合結(jié)構(gòu), 比如氣膜-發(fā)散冷卻[49]和逆噴-發(fā)散冷卻[50]等, 能明顯改善局部冷卻效果, 已成為近幾年的研究熱點, 同樣迎風(fēng)凹腔作為一種針對飛行器頭部駐點區(qū)域防熱卓有成效的方案, 自然也和發(fā)散冷卻存在著很強的互補性。 凹腔結(jié)構(gòu)可以有效降低駐點處的熱流, 解決發(fā)散冷卻中存在的頑疾問題; 而發(fā)散冷卻可以顯著擴大熱防護的范圍, 使凹腔防熱不僅僅局限于頭部駐點區(qū)域。

    欒蕓等[51]提出一種新型冷卻結(jié)構(gòu): 凹腔-發(fā)散組合冷卻, 結(jié)合了迎風(fēng)凹腔結(jié)構(gòu)與發(fā)散冷卻, 在上述機理的作用下, 取長補短。 圖13是楔形鼻錐凹腔-發(fā)散冷卻結(jié)構(gòu)示意圖, 在鼻錐的駐點區(qū)域切割出不可滲透的迎風(fēng)凹腔結(jié)構(gòu), 而多孔的發(fā)散面就是在迎風(fēng)凹腔的側(cè)壁面, 冷卻劑通入多孔材料中, 完成組合達到防熱的目的。

    文獻[51]還將單獨發(fā)散冷卻、 單獨迎風(fēng)凹腔以及組合構(gòu)型分別與純鼻腔進行冷卻效果對比, 對比結(jié)果如圖14所示, 可以看出, 迎風(fēng)凹腔發(fā)散冷卻組合構(gòu)型結(jié)合了兩者的優(yōu)點, 不僅可以提高駐點處的冷卻效果, 還能使鼻錐處降溫, 是一種高效的新型熱防護方式。 相較于另外兩種組合體的研究歷程, 發(fā)散冷卻與迎風(fēng)凹腔結(jié)合的方法無疑是最短的, 但試驗與數(shù)值模擬結(jié)果都不豐富, 其組合的合理性需要進一步深入探究。

    3 構(gòu)型評價及未來展望

    高超聲速飛行器熱防護技術(shù)直接關(guān)系到飛行器的飛行安全, 高效的減阻防熱方案亟待解決。 新材料的研發(fā)能直接提升飛行器的防熱性能, 此外, 新型熱防護技術(shù)的研究與組合方案的提出愈顯重要。

    基于上述對迎風(fēng)凹腔系列方案的介紹, 簡要總結(jié)這些思路的優(yōu)缺點。 迎風(fēng)凹腔這一方案的提出極大簡化了減阻防熱方案的復(fù)雜度與可實現(xiàn)難度, 相較于其他方案, 不需要預(yù)留通道, 不需要攜帶工質(zhì), 也不會導(dǎo)致前緣尖端存在極端熱環(huán)境, 當(dāng)然凹腔流動振蕩導(dǎo)致前緣受力的不穩(wěn)定以及相對較弱的冷卻能力是其不足之處。

    而三種組合方案都克服了凹腔本身的一些缺點, 有其各自的長處。 與逆向射流的組合方案, 凹腔提高了射流效能, 能獲得更好的防熱效果, 同時射流避免凹腔內(nèi)部流場振蕩。 總體而言, 兩種方案相輔相成, 有很強的可行性與實用性。 與能量沉積的組合方案, 高頻的能量源能夠抹去凹腔內(nèi)的振蕩變化, 使其更為易控, 但工程上實現(xiàn)難度較大。 與發(fā)散冷卻的結(jié)合方案中, 解決了冷卻效果偏弱的問題, 但在攻角存在下的可控性依然有待深入研究。

    結(jié)合減阻防熱研究的發(fā)展趨勢, 對于迎風(fēng)凹腔及其組合體的數(shù)值模擬與試驗, 得出以下三點展望:

    其一, 參考逆向射流系列方案, 學(xué)者們大多研究其工作參數(shù), 近些年才關(guān)注其結(jié)構(gòu)參數(shù), 如噴嘴形狀等。 對于迎風(fēng)凹腔系列方案, 同樣不能局限于單一長徑比的研究, 唇口鈍化半徑, 凹腔整體構(gòu)型也是值得深挖的結(jié)構(gòu)參數(shù)方向; 工作參數(shù)方面, 由于凹腔構(gòu)型的非定常性, 來流攻角的影響也極具參考價值。 對于總體設(shè)計而言, 不同參數(shù)的設(shè)計優(yōu)化可能存在矛盾, 解決此類問題, 多目標(biāo)設(shè)計優(yōu)化方法有其獨特的優(yōu)越性。

    其二, 前人的數(shù)值模擬過程中, 大多采用絕熱壁假設(shè), 但真實的熱環(huán)境下, 外部流場和內(nèi)部結(jié)構(gòu)也會相互作用, 考慮氣動熱效應(yīng), 耦合壁面與絕熱壁面所模擬出的結(jié)果會存在著很大差異, 因此數(shù)值模擬中采用流-熱耦合分析很有必要。 與此同時, 考慮到工程實際, 對于凹腔系列方案的研究必須在原先定常流動假設(shè)的基礎(chǔ)上開展非定常數(shù)值研究。

    其三, 所有的熱防護技術(shù)最終都要面向工程應(yīng)用, 因此與之匹配的試驗都需要開展, 一方面可以給數(shù)值模擬方法提供數(shù)據(jù)支持, 驗證模擬結(jié)果的正確性, 一方面也能檢驗熱防護方案的可行性, 加速該方法投入工程實際的進程。 高超聲速下的氣動加熱試驗本身就存在一定難度, 而迎風(fēng)凹腔又是一種不易試驗的構(gòu)型, 因此對此類方案試驗的開展提出了很高的要求。

    4 結(jié)? 論

    本文綜述了以迎風(fēng)凹腔為基礎(chǔ)的減阻防熱方案作用機理和研究現(xiàn)狀, 主要包括單一迎風(fēng)凹腔及其與逆向射流、 能量沉積、 發(fā)散冷卻等組合方案, 得出以下結(jié)論:

    (1) 迎風(fēng)凹腔是一個優(yōu)點與缺點同樣明顯的熱防護方案, 在應(yīng)用中要注意其非定常性, 通過組合方案來弱化其劣勢, 類似能量沉積這種應(yīng)用于工程實際比較難的組合方案, 應(yīng)盡力采用數(shù)值模擬的方法探究其特性, 也需要嘗試發(fā)掘更多新興的組合方式。

    (2) 迎風(fēng)凹腔與逆向射流的組合方案是凹腔組合方案中可行性最強的一種, 逆向射流可以避免凹腔內(nèi)部流場振蕩, 凹腔可以降低流場穩(wěn)定的臨界射流總壓, 提高射流效能, 兩者結(jié)合可達到揚長避短的效果。 下一步需要結(jié)合數(shù)值仿真和風(fēng)洞試驗方法對此類組合方案進行更深入的研究, 如改變傳統(tǒng)凹腔構(gòu)型、 引入減阻性能更優(yōu)秀的脈沖射流等, 尋找效能更高的組合構(gòu)型, 揭示其中蘊含的精細渦系結(jié)構(gòu)作用, 以此佐證地面試驗測試中觀察到的現(xiàn)象。

    (3) 高超聲速飛行器的減阻與防熱經(jīng)常是不能同時實現(xiàn)的, 高升阻比與飛行器前緣構(gòu)型的矛盾需要設(shè)計者很好的權(quán)衡, 利用基于代理模型的多目標(biāo)設(shè)計優(yōu)化算法來優(yōu)化迎風(fēng)凹腔的組合構(gòu)型, 了解設(shè)計參數(shù)和目標(biāo)功能之間的相互作用關(guān)系非常有價值, 可以對比選擇更有效的幾何模型。 除此之外, 還需要系統(tǒng)地研究工作參數(shù)和結(jié)構(gòu)參數(shù)對減阻防熱效果的影響, 這些信息對于熱防護系統(tǒng)的設(shè)計至關(guān)重要。

    (4) 隨著時代進步和技術(shù)發(fā)展, 新一代高超聲速飛行器必然會對熱防護系統(tǒng)提出更高要求, 這需要研究出更為有效的減阻防熱方案。 可以盡可能多地對迎風(fēng)凹腔組合方案進行風(fēng)洞試驗, 加速工程應(yīng)用的進程, 同時開展數(shù)值模擬時可以更多地采用流-熱耦合分析的方法, 更準(zhǔn)確地把握熱防護系統(tǒng)的特性。 除此之外, 需要拓寬思路, 除與現(xiàn)有的各種主動熱防護方案結(jié)合外, 另辟蹊徑, 探尋更為有效合理的熱防護方法。

    參考文獻:

    [1] Yang Y Z, Yang J L, Fang D N. Research Progress on Thermal Protection Materials and Structures of Hypersonic Vehicles[J]. Applied Mathematics and Mechanics, 2008, 29(1): 51-60.

    [2] 張志剛. 高超聲速飛行器熱防護系統(tǒng)設(shè)計方法[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2014.

    Zhang Zhigang. Design Method of Thermal Protection System on Hypersonic Vehicle[D]. Harbin: Harbin Institute of Technology, 2014. (in Chinese)

    [3] Schmisseur J D. Hypersonics into the 21st Century: A Perspective on AFOSR-Sponsored Research in Aerothermodynamics[J]. Progress in Aerospace Sciences, 2015, 72: 3-16.

    [4] 洪延姬, 李倩, 王殿愷, 等. 超聲速飛行器的激光空氣錐減阻方法[M]. 北京: 科學(xué)出版社, 2016.

    Hong Yanji, Li Qian, Wang Diankai, et al. Laser Air Cone Drag Reduction Method of Supersonic Vehicle[M]. Beijing: Science Press, 2016.(in Chinese)

    [5] Viviani A, Pezzella G. Aerodynamic and Aerothermodynamic Analy-sis of Space Mission Vehicles[M]. Cham: Springer International Publishing, 2015: 457-570.

    [6] Huang W. A Survey of Drag and Heat Reduction in Supersonic Flows by a Counterflowing Jet and Its Combinations[J]. Journal of Zhejiang University-Science A (Applied Physics & Engineering), 2015, 16(7): 551-561.

    [7] Shen B X, Liu W Q. Thermal Protection Performance of Opposing Jet Generating with Solid Fuel[J]. Acta Astronautica, 2018, 144: 90-96.

    [8] Qin Q H, Xu J L, Guo S. Fluid-Thermal Analysis of Aerodynamic Heating over Spiked Blunt Body Configurations[J]. Acta Astronautica, 2017, 132: 230-242.

    [9] Yadav R, Guven U. Aerodynamic Heating of a Hypersonic Projectile with Forward-Facing Ellipsoid Cavity at Nose[J]. Journal of Spacecraft and Rockets, 2014, 52(1): 157-165.

    [10] Ashwin Ganesh M, John B. Concentrated Energy Addition for Active Drag Reduction in Hypersonic Flow Regime[J]. Acta Astronautica, 2018, 142: 221-231.

    [11] Knight D. Survey of Aerodynamic Drag Reduction at High Speed by Energy Deposition[J]. Journal of Propulsion and Power, 2008, 24(6): 1153-1167.

    [12] Ahmed M Y M, Qin N. Recent Advances in the Aerothermodynamics of Spiked Hypersonic Vehicles[J]. Progress in Aerospace Sciences, 2011, 47(6): 425-449.

    [13] Huang W, Chen Z, Yan L, et al. Drag and Heat Flux Reduction Mechanism Induced by the Spike and Its Combinations in Supersonic Flows: A Review[J]. Progress in Aerospace Sciences, 2019, 105: 31-39.

    [14] Hartmann J. On a New Method for the Generation of Sound-Waves[J]. Physical Review, 1922, 20(6): 719-727.

    [15] Burbank P B, Stallings R L, Jr. Heat-Transfer and Pressure Measurements on a Flat Nose Cylinder at a Mach Number Range of 2.49 to 4.44 [R]. NASA-TX-X-19, 1959.

    [16] 張瑞瑞.逆向脈沖射流及其組合體在空天飛行器中的減阻防熱機理研究[D].長沙:國防科技大學(xué), 2018.

    Zhang Ruirui. Investigation of Drag and Heat Flux Reduction Mechanism of the Counterflowing Pulsed Jet and Its Combinations on Aerospace Vehicles[D].Changsha: National University of Defense Technology, 2018. (in Chinese)

    [17] 陸海波, 田世英. 迎風(fēng)凹腔: 一種有效的高超聲速飛行器熱防護選擇[J]. 飛航導(dǎo)彈, 2015(6): 11-15.

    Lu Haibo, Tian Shiying. Windward Cavity: An Effective Thermal Protection Option for Hypersonic Aircraft [J]. Aerodynamic Missile Journal, 2015(6): 11-15.(in Chinese)

    [18] Johnson R H. Instability in Hypersonic Flow about Blunt Bodies[J]. Physics of Fluids, 1959, 2(5): 526.

    [19] Saravanan S, Jagadeesh G, Reddy K P J. Investigation of Missile-Shaped Body with Forward-Facing Cavity at Mach 8[J]. Journal of Spacecraft and Rockets, 2009, 46(3): 577-591.

    [20] Yuceil B, Dolling D, Wilson D. A Preliminary Investigation of the Helmholtz Resonator Concept for Heat Flux Reduction[C]∥28th Thermophysics Conference,? 1993.

    [21] Yuceil K, Dolling D. IR Imaging and Shock Visualization of Flow over a Blunt Body with a Nose Cavity[C]∥34th Aerospace Sciences Meeting and Exhibit, 1996.

    [22] Yuceil K B, Dolling D S. Nose Cavity Effects on Blunt Body Pressure and Temperatures at Mach 5[J]. Journal of Thermophysics and Heat Transfer, 1995, 9(4): 612-619.

    [23] Engblom W A, Goldstein D B, Ladoon D, et al. Fluid Dynamics of Hypersonic Forward-Facing Cavity Flow[J]. Journal of Spacecraft and Rockets, 1997, 34(4): 437-444.

    [24] Engblom W A, Yuceil B, Goldstein D B, et al. Experimental and Numerical Study of Hypersonic Forward-Facing Cavity Flow[J]. Journal of Spacecraft and Rockets, 1996, 33(3): 353-359.

    [25] Lu H B, Liu W Q. Nummerical Simulation in Influence of Forward-Facing Cavity on Aerodynamic Heating of Hypersonic Vehicle[J]. Procedia Engineering, 2012, 29: 4096-4100.

    [26] 陸海波, 劉偉強. 高超聲速飛行器鼻錐迎風(fēng)凹腔結(jié)構(gòu)防熱效能研究[J]. 宇航學(xué)報, 2012, 33(8): 1013-1018.

    Lu Haibo, Liu Weiqiang. Investigation on Thermal Protection Efficiency of Hypersonic Vehicle Nose with Forward-Facing Cavity[J]. Journal of Astronautics, 2012, 33(8): 1013-1018.(in Chinese)

    [27] Ladoon D W, Schneider S P, Schmisseur J D. Physics of Resonance in a Supersonic Forward-Facing Cavity[J]. Journal of Spacecraft and Rockets, 1998, 35(5): 626-632.

    [28] Chou A, Schneider S P. Measurements of Resonance in a Forward-Facing Cavity at Mach Six[J]. Journal of Spacecraft and Rockets, 2015, 52(5): 1486-1494.

    [29] Silton S Z, Goldstein D B. Optimization of an Axial Nose-Tip Ca-vity for Delaying Ablation Onset in Hypersonic Flow[C]∥41st Aerospace Sciences Meeting and Exhibit, 2003.

    [30] Silton S I, Goldstein D B. Ablation Onset in Unsteady Hypersonic Flow about Nose Tip with Cavity[J]. Journal of Thermophysics and Heat Transfer, 2000, 14(3): 421-434.

    [31] Yadav R, Guven U. Aerothermodynamics of a Hypersonic Vehicle with a Forward-Facing Parabolic Cavity at Nose[J]. Journal of Aerospace Engineering, 2014, 228(10): 1863-1874.

    [32] 孫喜萬. 高超聲速再入飛行器頭部減阻防熱方案設(shè)計與優(yōu)化[D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 2016.

    Sun Xiwan. Design and Optimization of Drag and Heat Reduction Scheme of Hypersonic Re-Entry Vehicle Nose-Tip[D]. Changsha: National University of Defense Technology, 2016. (in Chinese)

    [33] Zheng Y Y, Ahmed N A. Thermal Protection Systems in Spacecraft Re-Entry:A Brief Overview[J]. JP Journal of Heat and Mass Transfer, 2013, 8(1): 99-118.

    [34] 王浚, 王佩廣. 高超聲速飛行器一體化防熱與熱控設(shè)計方法[J]. 北京航空航天大學(xué)學(xué)報, 2006, 32(10): 1129-1134.

    Wang Jun, Wang Peiguang. Integrated Thermal Protection and Control System Design Methodology for Hypersonic Vehicles[J]. Journal of Beijing University of Aeronautics and Astronautics, 2006, 32(10): 1129-1134.(in Chinese)

    [35] 鄧帆, 謝峰, 黃偉, 等. 逆向噴流技術(shù)在高超聲速飛行器上的應(yīng)用[J]. 空氣動力學(xué)學(xué)報, 2017, 35(4): 485-495.

    Deng Fan, Xie Feng, Huang Wei, et al. Applications of Counterflowing Jet Technology in Hypersonic Vehicle[J]. Acta Aerodynamica Sinica, 2017, 35(4): 485-495.(in Chinese)

    [36] Lu H B, Liu W Q. Research on Thermal Protection Mechanism of Forward-Facing Cavity and Opposing Jet Combinatorial Thermal Protection System[J]. Heat and Mass Transfer, 2014, 50(4): 449-456.

    [37] Lu H B, Liu W Q. Investigation of Thermal Protection System by Forward-Facing Cavity and Opposing Jet Combinatorial Configuration[J]. Chinese Journal of Aeronautics, 2013, 26(2): 287-293.

    [38] 陸海波. 迎風(fēng)凹腔與逆向噴流組合強化防熱結(jié)構(gòu)復(fù)雜流場和傳熱特性研究[D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 2012.

    Lu Haibo. Research on Complicated Flow Field and Heat Transfer Characteristic of Forward-Facing Cavity Combined with Opposing Jet Fortified Thermal Protection Configuration[D]. Changsha: National University of Defense Technology, 2012. (in Chinese)

    [39] Bibi A, Maqsood A, Sherbaz S, et al. Drag Reduction of Supersonic Blunt Bodies Using Opposing Jet and Nozzle Geometric Varia-tions[J]. Aerospace Science and Technology, 2017, 69: 244-256.

    [40] Sun X W, Guo Z Y, Huang W, et al. Drag and Heat Reduction Mechanism Induced by a Combinational Novel Cavity and Counterflowing Jet Concept in Hypersonic Flows[J]. Acta Astronautica, 2016, 126: 109-119.

    [41] Sudarshan B, Rao S M V, Jagadeesh G, et al. Effect of the Axial Cavity with an Opposing High-Pressure Jet Combination in a Mach 6 Flow Condition[J]. Acta Astronautica, 2021, 178: 335-348.

    [42] Zhang R R, Dong M Z, Huang W, et al. Drag and Heat Flux Reduction Mechanism Induced by the Combinational Forward-Facing Cavity and Pulsed Counterflowing Jet Configuration in Supersonic Flows[J]. Acta Astronautica, 2019, 160: 62-75.

    [43] Ahmed M Y M, Qin N. Forebody Shock Control Devices for Drag and Aero-Heating Reduction: A Comprehensive Survey with a Practical Perspective[J]. Progress in Aerospace Sciences, 2020, 112: 100585.

    [44] Adelgren R G, Yan H, Elliott G S, et al. Control of Edney IV Interaction by Pulsed Laser Energy Deposition[J]. AIAA Journal, 2005, 43(2): 256-269.

    [45] Borzov V Y, Mikhailov V M, Rybka I V, et al. Experimental Investigation of Supersonic Flow about an Obstacle with Power Deli-very into the Undisturbed Flow[J]. Journal of Engineering Physics and Thermophysics, 1994, 66(5): 449-454.

    [46] Bazyma L A, Rashkovan V M. Stabilization of Blunt Nose Cavity Flows by Using Energy Deposition[J]. Journal of Spacecraft and Rockets, 2005, 42(5): 790-794.

    [47] Langener T, von Wolfersdorf J, Selzer M, et al. Experimental Investigations of Transpiration Cooling Applied to C/C Material[J]. International Journal of Thermal Sciences, 2012, 54: 70-81.

    [48] Dittert C, Selzer M, Bhrk H. Flowfield and Pressure Decay Analysis of Porous Cones[J]. AIAA Journal, 2017, 55(3): 874-882.

    [49] Jiang P X, Huang G, Zhu Y H, et al. Experimental Investigation of Combined Transpiration and Film Cooling for Sintered Metal Porous Struts[J]. International Journal of Heat and Mass Transfer, 2017, 108: 232-243.

    [50] Shen B X, Liu W Q. Insulating and Absorbing Heat of Transpiration in a Combinational Opposing Jet and Platelet Transpiration Blunt Body for Hypersonic Vehicle[J]. International Journal of Heat and Mass Transfer, 2019, 138: 314-325.

    [51] 欒蕓, 賀菲, 王建華. 飛行器鼻錐凹腔-發(fā)散組合冷卻數(shù)值模擬[J]. 航空學(xué)報, 2021, 42(2): 623937.

    Luan Yun, He Fei, Wang Jianhua. Transpiration Cooling of Nose-Cone with Forward-Facing Cavity: Numerical Simulation[J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(2): 623937. (in Chinese)

    A Survey of Drag and Heat Reduction Induced by

    Forward-Facing Cavity and Its Combinations

    Zhang Jie, Xiao Feng, Huang Wei , Yan Li, Meng Yushan

    (Science and Technology on Scramjet Laboratory, College of Aerospace Science and Engineering,

    National University of Defense Technology, Changsha 410073, China)

    Abstract: Thermal protection technology is a key technology in the design of? hypersonic aircraft. In order to achieve more effective drag reduction and heat protection, scholars have proposed a variety of thermal protection schemes. The forward-facing cavity and its combinations is one of the most effective and promising active thermal protection schemes. This paper introduces the mechanism of single forward-facing cavity and its combination with counterflowing jet, energy deposition and divergent coding,? and summarizes their research status.

    It puts forward the outlook for the future, especially to follow up and innovate the existing combination schemes, and conduct more ground tests to support the numerical simulation results, with a view to early application in the engineering practice.

    Key words: ?drag reduction; thermal protection; forward-facing cavity; counterflowing jet; energy deposition; hypersonic vehicle; thermal protection technology

    91老司机精品| 久久久精品大字幕| 欧美乱码精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 老司机在亚洲福利影院| 日韩欧美精品v在线| 久久这里只有精品中国| 久久中文字幕人妻熟女| 非洲黑人性xxxx精品又粗又长| 级片在线观看| 成年版毛片免费区| 啦啦啦免费观看视频1| 天堂√8在线中文| 欧美在线一区亚洲| 精品午夜福利视频在线观看一区| 欧美乱码精品一区二区三区| 婷婷精品国产亚洲av在线| 久久婷婷人人爽人人干人人爱| 亚洲熟妇熟女久久| 国产精品久久久久久亚洲av鲁大| 99久久精品国产亚洲精品| 老鸭窝网址在线观看| 岛国在线免费视频观看| 一二三四社区在线视频社区8| av在线播放免费不卡| 少妇粗大呻吟视频| 亚洲专区中文字幕在线| 久久精品人妻少妇| 国产三级中文精品| 久久久久亚洲av毛片大全| 2021天堂中文幕一二区在线观| 亚洲精品一区av在线观看| 国产av在哪里看| 在线视频色国产色| 99re在线观看精品视频| 宅男免费午夜| 亚洲性夜色夜夜综合| 国产高清有码在线观看视频 | 免费在线观看日本一区| 不卡一级毛片| 在线免费观看的www视频| 久久久久久久精品吃奶| 蜜桃久久精品国产亚洲av| 久久久久免费精品人妻一区二区| 国产三级在线视频| 日韩三级视频一区二区三区| 1024手机看黄色片| 欧美一级a爱片免费观看看 | 精品无人区乱码1区二区| 怎么达到女性高潮| 国产一区二区三区视频了| 黄色女人牲交| 99riav亚洲国产免费| 欧美极品一区二区三区四区| 国产精品一区二区三区四区免费观看 | 中文字幕最新亚洲高清| 一边摸一边做爽爽视频免费| 免费av毛片视频| 久久精品国产亚洲av高清一级| 国产视频内射| 香蕉久久夜色| 97人妻精品一区二区三区麻豆| 亚洲无线在线观看| 99久久精品国产亚洲精品| 久久 成人 亚洲| 色哟哟哟哟哟哟| 香蕉av资源在线| cao死你这个sao货| 首页视频小说图片口味搜索| 日韩欧美在线二视频| av在线播放免费不卡| 亚洲成人国产一区在线观看| 国产亚洲欧美在线一区二区| 日本a在线网址| 青草久久国产| 午夜激情福利司机影院| 熟女少妇亚洲综合色aaa.| 色综合欧美亚洲国产小说| 欧美大码av| 国产一区二区三区视频了| av在线播放免费不卡| 51午夜福利影视在线观看| xxxwww97欧美| 国产精品乱码一区二三区的特点| 国产主播在线观看一区二区| 日日干狠狠操夜夜爽| 桃色一区二区三区在线观看| 久久久久免费精品人妻一区二区| 啦啦啦免费观看视频1| 久久人妻福利社区极品人妻图片| 国产麻豆成人av免费视频| 国产探花在线观看一区二区| 91字幕亚洲| 日韩三级视频一区二区三区| 亚洲人与动物交配视频| 韩国av一区二区三区四区| 国产精品久久久久久亚洲av鲁大| 国产黄色小视频在线观看| 亚洲免费av在线视频| 91在线观看av| 午夜福利欧美成人| 亚洲国产欧美一区二区综合| 丰满人妻熟妇乱又伦精品不卡| 成年版毛片免费区| 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| 亚洲人成网站高清观看| 老熟妇仑乱视频hdxx| 一二三四社区在线视频社区8| 欧美成狂野欧美在线观看| 精品一区二区三区av网在线观看| 99久久久亚洲精品蜜臀av| 亚洲18禁久久av| 毛片女人毛片| 国内精品久久久久精免费| 999精品在线视频| 亚洲熟妇中文字幕五十中出| 日韩欧美免费精品| 日本 av在线| 成人一区二区视频在线观看| 在线观看舔阴道视频| 日韩中文字幕欧美一区二区| 精品电影一区二区在线| 97人妻精品一区二区三区麻豆| 变态另类丝袜制服| 一边摸一边抽搐一进一小说| 国内精品久久久久久久电影| 国产精品亚洲av一区麻豆| 成人国语在线视频| 精品一区二区三区av网在线观看| 精品少妇一区二区三区视频日本电影| 欧美av亚洲av综合av国产av| 床上黄色一级片| 亚洲av第一区精品v没综合| 丝袜人妻中文字幕| 两人在一起打扑克的视频| www.精华液| 国产精品电影一区二区三区| 亚洲最大成人中文| 久久婷婷成人综合色麻豆| 五月伊人婷婷丁香| 19禁男女啪啪无遮挡网站| 男女视频在线观看网站免费 | 久久精品亚洲精品国产色婷小说| 熟女电影av网| 日韩欧美精品v在线| 亚洲精品中文字幕一二三四区| 久久久久国产一级毛片高清牌| 757午夜福利合集在线观看| avwww免费| 人成视频在线观看免费观看| 亚洲第一电影网av| 国产爱豆传媒在线观看 | 手机成人av网站| 99re在线观看精品视频| 国产高清激情床上av| 欧美一级a爱片免费观看看 | 日本 av在线| 丰满人妻一区二区三区视频av | 久久伊人香网站| 久久精品亚洲精品国产色婷小说| 久久婷婷人人爽人人干人人爱| 嫁个100分男人电影在线观看| 久久久精品大字幕| 久久国产乱子伦精品免费另类| 一夜夜www| www.精华液| 亚洲成人免费电影在线观看| 天天添夜夜摸| 亚洲第一电影网av| 国产精品影院久久| 免费搜索国产男女视频| 成人18禁高潮啪啪吃奶动态图| 国产一级毛片七仙女欲春2| 国产亚洲欧美在线一区二区| 午夜福利成人在线免费观看| 身体一侧抽搐| 在线观看美女被高潮喷水网站 | 久久久国产欧美日韩av| 欧美黑人巨大hd| 亚洲一区高清亚洲精品| 亚洲av成人精品一区久久| 亚洲国产欧美网| 欧美日韩精品网址| 啦啦啦免费观看视频1| 亚洲aⅴ乱码一区二区在线播放 | 动漫黄色视频在线观看| 国产蜜桃级精品一区二区三区| 舔av片在线| 色播亚洲综合网| 成人国语在线视频| 亚洲 欧美 日韩 在线 免费| 国产黄a三级三级三级人| 免费在线观看成人毛片| 国产精品98久久久久久宅男小说| 99在线人妻在线中文字幕| 精品一区二区三区四区五区乱码| 777久久人妻少妇嫩草av网站| av中文乱码字幕在线| 欧美一区二区国产精品久久精品 | 老汉色∧v一级毛片| 国产精品爽爽va在线观看网站| 久久草成人影院| av福利片在线观看| 欧美性猛交黑人性爽| 狂野欧美激情性xxxx| 小说图片视频综合网站| 床上黄色一级片| 首页视频小说图片口味搜索| 国产99白浆流出| 亚洲人成77777在线视频| 成人三级做爰电影| 12—13女人毛片做爰片一| 国产精品av久久久久免费| 精品欧美一区二区三区在线| 黄色 视频免费看| 国产精品自产拍在线观看55亚洲| 伦理电影免费视频| 好男人电影高清在线观看| 免费观看人在逋| 亚洲国产精品久久男人天堂| 久久欧美精品欧美久久欧美| 69av精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 男人舔女人的私密视频| 色精品久久人妻99蜜桃| 脱女人内裤的视频| 亚洲片人在线观看| 精品一区二区三区四区五区乱码| 久久香蕉精品热| 久久久久久国产a免费观看| 国产单亲对白刺激| 老司机靠b影院| 在线观看美女被高潮喷水网站 | 99re在线观看精品视频| 女警被强在线播放| 国产视频一区二区在线看| 久久久国产精品麻豆| 亚洲av电影不卡..在线观看| 97碰自拍视频| 国产精华一区二区三区| 日本熟妇午夜| 男女床上黄色一级片免费看| 成人三级做爰电影| 国产精品综合久久久久久久免费| 亚洲精品中文字幕一二三四区| 黑人欧美特级aaaaaa片| 日本五十路高清| 日本 欧美在线| 一级毛片精品| 亚洲激情在线av| 亚洲人成伊人成综合网2020| 亚洲专区国产一区二区| 欧美日韩中文字幕国产精品一区二区三区| 正在播放国产对白刺激| 久久人人精品亚洲av| 久久天躁狠狠躁夜夜2o2o| 可以在线观看毛片的网站| av国产免费在线观看| 国产亚洲av嫩草精品影院| 中文资源天堂在线| 99在线人妻在线中文字幕| 99久久综合精品五月天人人| 人妻夜夜爽99麻豆av| 丝袜美腿诱惑在线| 成人欧美大片| 91麻豆精品激情在线观看国产| 三级毛片av免费| 午夜激情av网站| 欧美一区二区国产精品久久精品 | 精品久久久久久,| 黄片大片在线免费观看| 又黄又爽又免费观看的视频| 亚洲无线在线观看| 午夜精品一区二区三区免费看| 搡老岳熟女国产| 久久久国产成人免费| 超碰成人久久| 一个人免费在线观看的高清视频| 亚洲在线自拍视频| 久久婷婷人人爽人人干人人爱| 国产成人aa在线观看| 777久久人妻少妇嫩草av网站| 久久久国产精品麻豆| 在线播放国产精品三级| 国产三级中文精品| 久久午夜亚洲精品久久| 亚洲专区国产一区二区| 欧美性猛交╳xxx乱大交人| 久久午夜综合久久蜜桃| 色噜噜av男人的天堂激情| 久99久视频精品免费| 国产黄片美女视频| 中文字幕熟女人妻在线| 中文字幕久久专区| 香蕉av资源在线| 国产69精品久久久久777片 | 国产精品爽爽va在线观看网站| 午夜久久久久精精品| 精品久久久久久,| 真人一进一出gif抽搐免费| 岛国在线观看网站| 欧美日韩亚洲综合一区二区三区_| 制服人妻中文乱码| 国产高清激情床上av| 欧美色欧美亚洲另类二区| 精品国产超薄肉色丝袜足j| а√天堂www在线а√下载| 亚洲精品av麻豆狂野| 国产三级在线视频| 亚洲精品中文字幕一二三四区| 99久久精品国产亚洲精品| 成人精品一区二区免费| 国产精品99久久99久久久不卡| 日韩欧美一区二区三区在线观看| 成年人黄色毛片网站| 国产片内射在线| 欧美成人午夜精品| a级毛片a级免费在线| 观看免费一级毛片| 国产又色又爽无遮挡免费看| 国产区一区二久久| 好男人电影高清在线观看| 一本久久中文字幕| 久久久久国产精品人妻aⅴ院| 老汉色av国产亚洲站长工具| 最近最新中文字幕大全电影3| 国产真人三级小视频在线观看| 热99re8久久精品国产| 精品不卡国产一区二区三区| 99久久无色码亚洲精品果冻| 日韩国内少妇激情av| 亚洲中文av在线| 亚洲天堂国产精品一区在线| 婷婷六月久久综合丁香| 可以免费在线观看a视频的电影网站| 午夜福利高清视频| 两性夫妻黄色片| 777久久人妻少妇嫩草av网站| 久久中文字幕一级| 一区二区三区国产精品乱码| 久久精品影院6| 国产av一区在线观看免费| 久久国产乱子伦精品免费另类| 日韩有码中文字幕| 每晚都被弄得嗷嗷叫到高潮| 午夜精品在线福利| 午夜两性在线视频| 日本 欧美在线| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美3d第一页| 精品人妻1区二区| 麻豆国产av国片精品| 小说图片视频综合网站| 色精品久久人妻99蜜桃| 久久久久亚洲av毛片大全| 99riav亚洲国产免费| 婷婷亚洲欧美| 日日干狠狠操夜夜爽| 成人欧美大片| 18禁黄网站禁片免费观看直播| 国产成人精品无人区| 色尼玛亚洲综合影院| 午夜福利高清视频| 五月伊人婷婷丁香| 不卡一级毛片| 国产av不卡久久| 精品无人区乱码1区二区| 免费一级毛片在线播放高清视频| 好男人电影高清在线观看| 久久中文字幕人妻熟女| √禁漫天堂资源中文www| 两个人的视频大全免费| 亚洲性夜色夜夜综合| 国产精品 欧美亚洲| 国产精品综合久久久久久久免费| 丝袜美腿诱惑在线| 日本在线视频免费播放| 久久久久久九九精品二区国产 | 中文字幕最新亚洲高清| 999久久久国产精品视频| 我要搜黄色片| svipshipincom国产片| 国内少妇人妻偷人精品xxx网站 | 三级毛片av免费| 老司机靠b影院| 久久性视频一级片| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 久热爱精品视频在线9| 国产成人精品无人区| 老司机深夜福利视频在线观看| 欧美性猛交╳xxx乱大交人| 精品高清国产在线一区| 国产三级黄色录像| 一级毛片高清免费大全| 久久久久久久午夜电影| 精品高清国产在线一区| 淫秽高清视频在线观看| 亚洲国产欧美网| 久久精品综合一区二区三区| 中文字幕高清在线视频| 久久久国产精品麻豆| 黄色毛片三级朝国网站| 亚洲中文av在线| 成人欧美大片| 色综合站精品国产| 人人妻人人看人人澡| 琪琪午夜伦伦电影理论片6080| 日日爽夜夜爽网站| av有码第一页| 久久久久久久久中文| 床上黄色一级片| 精品国产乱子伦一区二区三区| 国产成年人精品一区二区| 欧美久久黑人一区二区| 美女高潮喷水抽搐中文字幕| 最新在线观看一区二区三区| 欧美精品亚洲一区二区| 亚洲欧洲精品一区二区精品久久久| 免费无遮挡裸体视频| 午夜免费观看网址| av视频在线观看入口| 色尼玛亚洲综合影院| 日本免费一区二区三区高清不卡| 亚洲aⅴ乱码一区二区在线播放 | 搡老妇女老女人老熟妇| 欧美精品啪啪一区二区三区| 免费人成视频x8x8入口观看| 三级男女做爰猛烈吃奶摸视频| a级毛片在线看网站| 免费无遮挡裸体视频| 夜夜爽天天搞| 国产精品亚洲一级av第二区| 国产亚洲av嫩草精品影院| 在线观看www视频免费| 久久久久九九精品影院| 国产99白浆流出| 后天国语完整版免费观看| 久久久久久免费高清国产稀缺| 91国产中文字幕| 动漫黄色视频在线观看| 一级作爱视频免费观看| 男人舔奶头视频| 精品一区二区三区四区五区乱码| 午夜福利18| 国产亚洲欧美98| 国产三级在线视频| 91麻豆av在线| 国产精品av视频在线免费观看| 日本三级黄在线观看| 亚洲午夜精品一区,二区,三区| 欧美日韩亚洲综合一区二区三区_| 国产伦在线观看视频一区| 日韩欧美免费精品| 免费在线观看视频国产中文字幕亚洲| 久久精品国产亚洲av高清一级| 曰老女人黄片| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久久久毛片| 91国产中文字幕| 两个人的视频大全免费| 精品无人区乱码1区二区| 免费看a级黄色片| 女人被狂操c到高潮| 国产一区在线观看成人免费| 久久伊人香网站| 宅男免费午夜| 国产激情欧美一区二区| 精品久久久久久久毛片微露脸| 搡老熟女国产l中国老女人| 免费观看人在逋| 少妇熟女aⅴ在线视频| 国产成人精品无人区| 国产精品免费视频内射| 给我免费播放毛片高清在线观看| 久久久久久免费高清国产稀缺| 国产精品 国内视频| 亚洲天堂国产精品一区在线| 青草久久国产| 国产亚洲精品一区二区www| 女生性感内裤真人,穿戴方法视频| 欧美国产日韩亚洲一区| 午夜福利欧美成人| 给我免费播放毛片高清在线观看| 一级黄色大片毛片| a级毛片在线看网站| 日韩国内少妇激情av| 国产精品精品国产色婷婷| 黄色毛片三级朝国网站| 中文字幕人妻丝袜一区二区| 黄色丝袜av网址大全| 亚洲成人中文字幕在线播放| 好男人在线观看高清免费视频| 在线观看免费视频日本深夜| 亚洲国产精品成人综合色| 欧美乱码精品一区二区三区| 国产黄a三级三级三级人| 久久人妻av系列| 国产一区二区在线观看日韩 | 亚洲午夜精品一区,二区,三区| 欧美日韩乱码在线| 久久久久久久久中文| 欧美日韩中文字幕国产精品一区二区三区| 观看免费一级毛片| 一进一出好大好爽视频| 法律面前人人平等表现在哪些方面| 婷婷精品国产亚洲av在线| 国产精华一区二区三区| 午夜激情av网站| 国产精品av久久久久免费| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 久久草成人影院| 老熟妇仑乱视频hdxx| 午夜福利在线观看吧| 搡老熟女国产l中国老女人| 99精品欧美一区二区三区四区| 亚洲男人的天堂狠狠| 亚洲av成人不卡在线观看播放网| 成人手机av| 一区二区三区激情视频| 国产成人影院久久av| 手机成人av网站| 免费搜索国产男女视频| 欧美日韩国产亚洲二区| 制服丝袜大香蕉在线| 国产精品亚洲美女久久久| 国内毛片毛片毛片毛片毛片| 色综合婷婷激情| 成在线人永久免费视频| 国产亚洲精品综合一区在线观看 | 观看免费一级毛片| 草草在线视频免费看| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 亚洲熟妇熟女久久| 欧美日韩亚洲国产一区二区在线观看| 国产野战对白在线观看| 99热6这里只有精品| 亚洲全国av大片| 欧美一级毛片孕妇| 国产精品一及| 久久 成人 亚洲| 国产黄色小视频在线观看| 亚洲欧美精品综合一区二区三区| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| 天堂影院成人在线观看| 亚洲中文av在线| 又黄又爽又免费观看的视频| 香蕉国产在线看| 亚洲专区中文字幕在线| 国产精品美女特级片免费视频播放器 | 国产片内射在线| 亚洲片人在线观看| 最新在线观看一区二区三区| 国内精品久久久久精免费| 亚洲国产中文字幕在线视频| 日韩欧美精品v在线| 精品欧美一区二区三区在线| 日韩精品免费视频一区二区三区| 午夜成年电影在线免费观看| 午夜免费观看网址| 50天的宝宝边吃奶边哭怎么回事| 天堂动漫精品| 国产激情欧美一区二区| 午夜福利免费观看在线| 女人爽到高潮嗷嗷叫在线视频| 制服丝袜大香蕉在线| www日本在线高清视频| 国内毛片毛片毛片毛片毛片| 搞女人的毛片| АⅤ资源中文在线天堂| 国产1区2区3区精品| 99国产综合亚洲精品| 老司机深夜福利视频在线观看| 日日干狠狠操夜夜爽| 成人高潮视频无遮挡免费网站| а√天堂www在线а√下载| 久久精品国产综合久久久| 亚洲av电影在线进入| 久久久久免费精品人妻一区二区| 搡老妇女老女人老熟妇| 亚洲最大成人中文| 欧美日韩乱码在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人欧美在线观看| 亚洲精品色激情综合| 亚洲一区二区三区不卡视频| a在线观看视频网站| 亚洲精品在线观看二区| 级片在线观看| 成在线人永久免费视频| 国产视频一区二区在线看| 在线看三级毛片| 少妇被粗大的猛进出69影院| 久久久久久久久中文| 久久精品91无色码中文字幕| 午夜福利高清视频| 国产成人啪精品午夜网站| 国产蜜桃级精品一区二区三区| 中文在线观看免费www的网站 | 成人av一区二区三区在线看| 成人手机av| 国产一级毛片七仙女欲春2| 在线观看日韩欧美| av中文乱码字幕在线| 婷婷丁香在线五月| 亚洲成人国产一区在线观看| 丁香六月欧美| 嫩草影院精品99| 黑人操中国人逼视频| 91成年电影在线观看| 中文字幕久久专区|