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

    高超聲速雙楔激波干擾定常射流控制試驗研究

    2024-05-09 10:16:40謝瑋羅振兵周巖劉強吳建軍董昊
    航空學報 2024年7期
    關鍵詞:弓形來流氣源

    謝瑋,羅振兵,*,周巖,劉強,吳建軍,董昊

    1.國防科技大學 空天科學學院,長沙 410073

    2.南京航空航天大學 航空學院,南京 210016

    飛行器在超聲速/高超聲速飛行時,機體及各氣動部件(鈍前緣、進氣道、鈍舵、尾翼等)均會出現(xiàn)激波,激波干擾即是指由激波引發(fā)的干擾現(xiàn)象[1-3],一般包括激波/邊界層干擾與激波/激波干擾。激波/邊界層干擾是指激波與邊界層之間的干擾現(xiàn)象,對于實際飛行器,激波干擾更加復雜,往往還會產生多道激波之間的干擾,即激波/激波干擾。本文主要關注激波/激波干擾問題,但激波/激波干擾中也不可避免存在著激波/邊界層干擾[4-5],因此后文以激波干擾進行描述。

    從幾何空間上看,常見的激波干擾有兩種產生形式:一種是由飛行器氣動部件自身誘導產生的多道激波之間的相互干擾,例如雙楔構型中的兩道同族激波相互干擾、鈍舵構型中脫體激波與分離激波的相互干擾等;另一種是各氣動部件之間的激波相互干擾,例如飛行器的頭部激波與下游進氣道唇口處弓形激波以及尾翼/襟翼/鈍舵前緣激波之間的相互干擾等[4-7]。

    Edney[8]早期以斜激波入射弓形激波為入手點對激波干擾問題進行了系統(tǒng)的分析,研究發(fā)現(xiàn)激波干擾流場中會產生激波、膨脹波、剪切層和超聲速射流等復雜流動結構,導致壁面壓力和熱流急劇上升;同時,他們根據(jù)激波入射位置以及產生的流場結構的不同,將激波干擾現(xiàn)象分成6類,其中第Ⅳ類干擾現(xiàn)象最為嚴重,所產生的壓力和熱流可達駐點數(shù)倍。Olejniczak等[9]在高超聲速雙楔流動中也發(fā)現(xiàn)了4 種激波干擾類型,其中第Ⅳ、Ⅴ、Ⅵ類與Edney 的分類一致,此外還新發(fā)現(xiàn)了一種第Ⅳr 類干擾。Edney[8]和Olejniczak等[9]針對斜激波入射鈍頭和雙楔的兩種分類奠定了激波干擾分析和研究的基礎,這兩類構型也成為了激波干擾問題研究的典型構型[10-14]。

    國際上以美國為代表的眾多國家和單位,包括美國國家航空航天局(National Aeronautics and Space Administration,NASA)蘭利研究中心[15]、卡爾斯本大學巴法羅研究中心[16]、伊利諾伊大學[17]、羅格斯-新澤西州立大學[18]等對激波干擾致力致熱特性開展了持續(xù)性研究。例如,NASA 蘭利研究中心的Keyes和Hains 等的試驗研究表明,第Ⅲ類和第Ⅳ類斜激波入射型激波干擾會造成嚴重的氣動載荷,壓力和熱流極值最大約為無干擾時的7.5和17.0倍[15]。國內中國科學院力學研究所[19]、中國科學技術大學[20-21]、中國空氣動力研究與發(fā)展中心[22]以及國防科技大學[23]等單位也對激波干擾問題開展了研究。近期,中國科學院力學研究所彭俊[19]在JF-12 高焓高超聲速風洞中開展了大尺寸激波干擾試驗;結果表明激波干擾作用下壓力極值增大了5倍,熱流極值增大了12倍。由此可見采取適當?shù)牧鲃涌刂品绞綄げǜ蓴_區(qū)域實施控制具有重要意義。

    然而,目前對于激波干擾問題的主動流動控制研究主要針對激波/邊界層干擾控制[24-27],對于激波干擾問題的主動流動控制研究還較少。Albertson和Venkat[28]初步開展了第Ⅲ、Ⅳ類激波干擾逆向射流控制試驗研究;結果表明逆向射流可以顯著改變第Ⅲ、Ⅳ類激波干擾結構,避免其直接撞擊唇口前緣,但遺憾的是他們未開展定量的試驗測量。吳文堂[29]、王殿愷[30]等開展了高重頻激光能量沉積控制第Ⅳ類激波干擾的數(shù)值模擬研究,結果表明典型工況中頻率150 kHz、單脈沖能量5 mJ 的高重頻激光可實現(xiàn)鈍頭體表面峰值駐點壓力、熱流和波阻分別降低40%、33%和23%。謝瑋等[31]開展了等離子體合成射流高焓高超聲速流場雙楔第Ⅵ類激波干擾控制試驗,結果表明等離子體合成射流可以顯著消除二楔斜激波和激波干擾點。唐孟瀟[32]、孔亞康[33]和張傳標[34]等基于陣列表面電弧放電激勵器也開展了雙楔第Ⅵ類激波干擾控制,取得了類似的控制效果。但由于放電過程中強電磁干擾的影響,上述試驗中均未能測得有效定量試驗數(shù)據(jù)。綜上所述,已開展的激波干擾主動流動控制研究結果缺少試驗定量驗證,且控制效果也有待進一步提升。因此,本文采用紋影和熱流測量技術,針對雙楔第Ⅴ類、第Ⅳ類激波干擾開展了定常射流控制試驗研究,以驗證定常射流對激波干擾流場的控制作用及降熱效果。射流產生方式主要為高壓氣源,此外,基于課題組高速來流能量利用思想[3],提出了雙楔激波干擾自持射流控制方法。

    1 試驗設置

    1.1 風洞及試驗模型

    試驗在南京航空航天大學?500 mm 高超聲速風洞(Nanjing University of Aeronautics and Astronautics ?500 mm Hypersonic Wind Tunnel,NHW)中開展,NHW 是一座高壓下吹-真空吸氣暫沖式高超聲速風洞。風洞試驗來流馬赫數(shù)Ma∞=8,來流總壓P0=0.94 MPa,來流總溫T0=518 K,單位雷諾數(shù)Reunit=3.9×106m-1,有效試驗時間約為7 s。風洞試驗共采用了6 個風洞試驗模型,即30-60 基準、30-60氣源、30-60自持、30-75基準、30-75氣源和30-75自持試驗模型?!?0”是指一楔角度為30°,“60/75”是指二楔角度為60°/75°,“基準”是指無射流控制的雙楔試驗模型,“氣源”是指通過高壓氣源產生氣源射流的雙楔試驗模型,“自持”是指通過集氣通道收集高速來流產生自持射流的雙楔試驗模型。

    以30-60 氣源模型為例展示試驗設置。如圖1(a)和圖1(b)所示(L 表示長度,?表示直徑),雙楔寬60 mm,一楔長L1=169.9 mm,二楔長L2=89.7 mm。在一楔上設置40 mm×2 mm長方形縫狀射流出口,射流出口距離一楔、二楔交界處24.5 mm,射流噴射方向與一楔表面成60°夾角。當前設置模型的展向寬度偏窄,若模型加寬,由于射流展向控制范圍有限,因此射流出口展向寬度也需要相應加寬,而射流從氣源流向射流出口,壓力損失大,且射流出口寬度越大,損失越大,易導致形成的射流偏弱。此外,若試驗模型過寬,風洞堵塞度可能過大。射流出口位置則是在滿足熱電偶安裝的情況下盡可能靠近受控激波干擾區(qū),以期取得更好的控制效果。射流出口下游模型中線上共布置15 個?2 mm 熱電偶安裝孔用于安裝熱電偶,其中一楔共布置3個,間距5 mm;二楔共布置12個,間距6 mm,熱電偶用于測量射流控制前后激波干擾區(qū)熱流變化。射流入口為直徑7 mm 的圓柱型通道,射流入口插入內徑4 mm、外徑6 mm 圓截面的高壓氣管,隨后射流經過寬為40 mm、厚為4 mm 的長方形截面過渡水平通道,再流入豎直方向寬為40 mm、厚為2 mm 的射流通道,并從射流出口噴出。過渡水平通道的設計即是為了使射流在流入豎直射流通道前基本轉變成縫隙狀,從而使縫狀射流更加均勻??紤]到安裝問題,在距離射流出口約32 mm 的位置設置了壓力傳感器安裝孔,用于安裝NS-3 壓力傳感器測量射流出口靜壓。

    圖1 試驗設置Fig.1 Experimental setup

    30-60 自持模型將射流通道改為集氣通道,如圖1(c)和圖1(d)所示。自持射流的基本原理類似于進氣道,即充分利用激波增壓和來流動能增壓,對于本文試驗模型,來流首先通過雙楔模型下方斜激波進行了初步增壓,此時來流速度仍然為超聲速,仍然具有較高的動能;隨后來流進入集氣通道,由于集氣通道為收縮型結構,來流將進一步減速,由能量守恒方程可知,高速來流的動能進一步轉化為壓力能;但此時,來流速度可能仍然較高,因此在過渡通道可能會出現(xiàn)滯止激波,進一步減速增壓,并經過射流通道在射流出口形成自持射流。

    此外,30-60 基準模型相比30-60 氣源模型取消了射流通道和射流孔,30-75 系列3 個模型將二楔角度改為75°,其他設置與對應的30-60 系列模型保持一致。以30-60 自持模型為例展示了原點和坐標軸的設定,如圖1(c)和圖1(d)所示,原點為一楔和二楔交線的中點,X 方向與一楔平行并指向流向下游,Y 方向垂直于一楔,因此在計算二楔上熱電偶X 坐標時,30-60 模型和30-75 模型需要分別按30°和45°進行投影。

    風洞試驗工況如表1 所示,其中工況1-1 采用的是30-60 基準模型;工況1-2~1-8 采用的是30-60 氣源模型,通過高壓氣源(常溫)產生不同壓比的氣源定常射流進行流動控制,射流壓比(Pressure Ratio,PR)定義為射流出口靜壓與一楔斜激波波后靜壓(約2.9 kPa)的比值;工況1-9為30-60 自持模型,通過集氣通道集氣產生自持射流進行流動控制;工況2-1 為30-75 基準模型,工況2-2~2-11 為30-75 氣源模型,工況2-12 為30-75 自持模型。

    表1 試驗工況Table 1 Experimental cases

    1.2 試驗測量裝置

    熱電偶采用的是中國科學院力學研究所研制的E 型同軸熱電偶,如圖2 所示,該熱電偶是基于表面溫度測量熱傳導理論研制的熱流傳感器,通過測定溫度變化歷程反算得到熱流。熱電偶熱流信號響應時間小于5 μs,溫度測量范圍為-100~500 ℃,熱流測量范圍為0~50 MW/m2,熱流測量精度<7%。

    圖2 E 型同軸熱電偶Fig.2 E-type coaxial thermocouple

    傳感器測量熱流時,首先測得溫度隨時間的變化曲線,隨后換算成熱流值,熱流計算公式為

    式中:T 為溫度;t 為時刻;τ 為時間變量。

    離散后的熱流計算公式為

    式中:n 為采集點序列數(shù);i 為循環(huán)數(shù);ti為時刻;T(ti)為ti時刻熱電偶輸出溫度值,其余下標含義以此類推;ρck 為熱電偶平均物性參數(shù),試驗采用的熱電偶(ρck)1/2=8 160 J/(s1/2·m2·K)。

    采用上海天沐公司的NS-3 壓力傳感器測量試驗模型射流出口處靜壓,傳感器量程為0~400 kPa,綜合測試精度為±0.25%F.S.(F.S.即Full Scale,滿量程)。

    風洞試驗段兩側設置有兩個?350 mm 光學玻璃窗口,試驗中通過該窗口進行了紋影拍攝。紋影系統(tǒng)包括光源、兩個凹面鏡、兩個平面鏡、刀口以及高速相機,采用典型的Z 型布置,相機采用的是ACS-3 高速相機,拍攝幀頻為1 000 fps,曝光時間為1.1 μs。

    數(shù)據(jù)采集系統(tǒng)采用的是東華DH8302 動態(tài)信號測試分析系統(tǒng),該數(shù)據(jù)采集儀量程范圍5~10 000 mV、采樣頻率最高可達1 MHz、噪聲<3 μVrms,適用于低信號值(63 μV/K)熱電偶的數(shù)據(jù)采集。此外DH8302 配備有熱電偶采集模塊,便于試驗測量,試驗中共使用了15 個采集通道。

    2 無控激波干擾流場與熱流分析

    30-60 基準模型和30-75 基準模型紋影、流場結構示意圖和熱流測量結果如圖3、圖4 所示,熱流測量結果的橫坐標為熱電偶測點的幾何相對位置X/L1(X 為各熱電偶測點的橫坐標,L1為一楔長度)。對比文獻[9]基于無黏流場的分類可知,30-60 基準模型中流場結構為典型的第Ⅴ類激波干擾。一楔產生的斜激波IP 與二楔產生的弓形激波PW 相互干擾,產生透射激波PT、馬赫桿TQ,三波點P和T 發(fā)出的兩道滑移線C1和C2形成射流邊界,包裹著超聲速射流[4,9]。此外,一楔和二楔之間產生了一個分離區(qū),分離區(qū)上游產生了分離激波SO,且分離激波末端打在透射激波PT上,分離區(qū)下游產生再附激波BQ,與馬赫桿TQ 相連接。需要指出的是,文獻[4,9]中,三波點Q 發(fā)出的滑移線包圍著超聲速氣流向下游移動,QR 是反射激波,下游可能還存在相應的激波結構,但在試驗中上述流動現(xiàn)象未能通過紋影觀測到。從熱流沿流向的變化結果可以看出,激波干擾區(qū)的熱流顯著高于非激波干擾區(qū),熱流極值位于激波干擾區(qū)與二楔壁面的接觸區(qū)域,即X/L1=0.08(測點6,再附點R和反射激波點B 之間),約91.9 kW/m2。

    圖3 30-60 基準和30-75 基準模型紋影測量結果及流場結構示意圖Fig.3 Schlieren measurement results and schematic diagrams of flow structures for 30-60 reference and 30-75 reference models

    圖4 30-60 基準和30-75 基準模型熱流測量結果Fig.4 Heat flux measurement results for 30-60 reference and 30-75 reference models

    由于本文模型展向寬度偏窄,可能存在著三維效應的影響,為了確認流場波系結構,利用激波極線[4]方法對30-60 基準流場進行了理論分析。如圖5 所示,圖中θ 為氣流偏轉角,P/P∞為流場靜壓與來流靜壓的比值,激波極線圖中某一特定的點即對應于圖3 流場示意圖相應區(qū)域的流場參數(shù)。來流馬赫數(shù)Ma∞=8,原點即點對應來流參數(shù)。一楔后馬赫數(shù)Ma1=2.58,流場參數(shù)由跳至①點。分離激波波后馬赫數(shù)Ma2=2.41,流場參數(shù)由①點跳至②點。透射激波PT 波后馬赫數(shù)Ma3=1.98,流場參數(shù)由弱解的②點跳至強解的③④點,其中④點是弓形激波PW 后的流場參數(shù),透射激波PT 的作用即是對②區(qū)的流場進一步增壓,以實現(xiàn)③區(qū)和④區(qū)的壓力匹配。斜激波QB 波后馬赫數(shù)Ma7=1.25,流場參數(shù)由②點跳至⑦點。激波極線中①⑦點位置由一楔角30°、二楔角60°確定,②點由紋影圖中分離激波角(55.8°)確定。③④點位于②區(qū)分離激波極線與區(qū)來流激波極線交點處,與之類似,⑤⑥點位于③區(qū)激波極線與②區(qū)分離激波極線交點處,⑧⑨點位于②區(qū)分離激波極線與⑦區(qū)激波極線交點處。從激波極線的分析及與文獻[9]無黏結果的對比可以確認,圖3給出的第Ⅴ類激波干擾流場示意圖是合理的。后文中第Ⅳ類激波干擾基準流場及射流控制流場參考30-60 基準模型給出了流場結構示意圖,但由于黏性影響大且彎曲激波較多,僅在射流噴出等影響較小的區(qū)域利用激波極線進行了一定的輔助分析。

    圖5 30-60 基準流場對應激波極線Fig.5 Shock polar diagram for 30-60 reference flow field

    30-75 基準模型中流場結構則為典型的第Ⅳ類激波干擾。此時,斜激波IP和弓形激波PW 干擾產生的超聲速射流結構直接沖擊二楔壁面,并在沖擊處產生高熱流,即X/L1=0.12(測點8,滑移線C1和C2與壁面交點之間),熱流極值約為111.1 kW/m2,高于30-60 基準模型中第Ⅴ類激波干擾產生的熱流極值。此外,除30-60 基準模型激波干擾區(qū)對應的幾個測點,30-75 基準模型各測點對應的熱流均要更高。由于二楔角度增大,分離區(qū)明顯變大,此時分離激波SO 與一楔斜激波IP 相互干擾,相交于點S1并導致斜激波的后半段S1P 角度增大,而分離激波SO 的后半段角度則顯著變小,且略為彎曲。需要指出的是,紋影圖中馬赫桿TQ直接與分離區(qū)相連接,且未觀察到明顯的再附激波,再附激波可能與馬赫桿TQ 融合,由此可見分離區(qū)對激波干擾區(qū)域影響較大。試驗中第Ⅴ、Ⅳ類激波干擾流場中均觀察到了分離激波的非定常運動,但激波干擾區(qū)域流場結構及熱流測量結果隨時間振蕩變化較小,因此分析認為分離激波的非定常運動對本文研究影響較小。

    3 雙楔激波干擾射流控制

    3.1 第Ⅴ類激波干擾控制

    接著開展了雙楔激波干擾射流控制研究,首先是30-60 模型中第Ⅴ類激波干擾的控制,紋影及示意圖如圖6、圖7 所示。隨著氣源射流的噴出,射流誘導產生一道弓形射流激波JP2,射流和射流激波共同對下游激波干擾區(qū)產生控制作用,使激波干擾流場發(fā)生顯著變化,不再呈現(xiàn)典型的第Ⅴ類激波干擾結構。弓形射流激波JP2與一楔斜激波IP1相交于P1,并相互干擾產生膨脹波系P1L[9,32]。而弓形 射流激波JP2與二楔 弓形激 波P2W 相互干擾,誘導產生透射激波P2N和馬赫桿P2Q。射流與來流之間形成剪切層JM,在一定程度上隔絕了激波干擾區(qū)與壁面之間的干擾。強射流使得分離區(qū)及分離激波SS1向上游移動,并與一楔斜激波相交于S1點。由于射流寬度沒有展向全覆蓋,因此紋影流場中可能存在三維效應的影響,因此采取激波極線方法進行輔助分析驗證是必要的。然而,對于射流控制下雙楔復雜流場,黏性的影響難以消除,且流場中包含較多的彎曲激波,因此綜合考慮,針對可能出現(xiàn)三維效應的射流出口附近激波干擾區(qū)進行理論分析。以PR=1.2 流場為例,如圖6(a)、圖7(a)、圖8,來流馬赫數(shù)Ma∞=8,原點即點對應來流參數(shù)。一楔后馬赫數(shù)Ma1=2.58,流場參數(shù)由點跳至①點。分離激波角相比基準流場略微增大,約為56.6°,波后馬赫數(shù)Ma2=2.37,流場參數(shù)由①點跳至②點。將弓形射流激波JP2靠近射流出口部分近似看做直線斜激波,測得激波角度約為70°,其波后馬赫數(shù)Ma3=1.88,流場參數(shù)由②點跳至③點。從圖8 激波極線可以看出,弓形激波波后激波極線未與來流激波極線相交,即③區(qū)靜壓高于⑤區(qū)靜壓,因此,為了實現(xiàn)壓力匹配,產生了膨脹波P1L。而P1K 理論上應為滑移線,紋影圖中對應位置可能是紋影圖三維積分效應的影響。而對于P2點附近的激波干擾,主要由射流弓形激波JP2誘導,因此,分析認為三維效應影響較小。由于射流壓比越大,弓形射流激波角越大,因此③區(qū)與⑤區(qū)的壓力差越大,膨脹波越強,但其對應射流噴出區(qū)域的流場結構均一致。

    圖6 第Ⅴ類激波干擾氣源/自持射流控制紋影測量結果Fig.6 Schlieren measurement results of Type-Ⅴ shock interaction control by steady/self-sustaining jets

    圖7 第Ⅴ類激波干擾射流控制示意圖Fig.7 Schematic diagrams of Type-Ⅴ shock interaction control by jets

    圖8 第Ⅴ類激波干擾射流控制流場(PR=1.2)激波極線Fig.8 Shock polar diagram of Type-Ⅴ shock interaction control flow field by jets(PR=1.2)

    隨著射流壓比的變化,P2點附近流場結構有一定的變化,并對降熱效果有顯著的影響。當氣源射流壓比較小時(PR≤4.8,圖6(a)~圖6(c)、圖7(a)),弓形射流激波JP2相對較弱,對馬赫桿P2Q 的強度影響較小,馬赫桿P2Q 仍與壁面有著較強的相互干擾,馬赫桿P2Q 所在位置附近也是壁面高熱流所在,如圖9 所示。但可以發(fā)現(xiàn),隨著射流壓比增大,馬赫桿P2Q 逐漸被射流推向下游,雙楔表面熱流值也相應改變。如圖9 所示,基準工況中熱流極值出現(xiàn)在X/L1=0.08(測點6),壓比1.2和3.0 時熱流極值均出現(xiàn)在X/L1=0.17(測點9),而壓比4.8 時熱流極值則出現(xiàn)在X/L1=0.20(測點10)。當射流壓比較大時(PR>4.8,圖6(d)~圖6(g)、圖7(b)),弓形射流激波JP2與二楔弓形激波P2W 相互干擾,在馬赫桿P2Q 上游產生透射激波P2Q1并與剪切層JM 相交于Q1點。在透射激波P2Q1和強射流作用下,馬赫桿P2Q 被顯著削弱,其原理類似于將強激波分散成多個弱激波系的控制方法。相應的,熱流也顯著降低,不再存在明顯的極值,如圖9 所示,熱流較高區(qū)域均出現(xiàn)在X/L1=0.27、0.30(測點12、13)??傮w而言,熱流極值隨著射流壓比的增大而不斷降低,相比基準工況時的91.9 kW/m2,射流壓比12.0 時的熱流極值僅為17.3 kW/m2,降低了約81.2%。

    圖9 第Ⅴ類激波干擾氣源/自持射流控制熱流測量結果Fig.9 Heat flux measurement results of Type-Ⅴ shock interaction control by steady/self-sustaining jets

    結合不同射流壓比下的紋影、示意圖及熱流測量結果,分析認為氣源射流控制下雙楔表面熱流降低主要有兩方面原因:一是隔絕作用,即氣源射流噴射到激波干擾流場中后,覆蓋了雙楔壁面并與來流相互作用產生剪切層,在一定程度上隔絕了激波干擾區(qū)域與壁面之間的熱量傳導;二是結構改變作用,即氣源射流顯著改變了激波干擾結構,避免了局部強激波干擾區(qū)的產生。不同射流壓比下兩種因素起的作用不同:當射流壓比較小時,激波干擾結構雖然有所改變,但激波干擾作用仍然較強,隔絕作用占主導;而當射流壓比較大時,兩種因素共同起重要作用,激波與壁面的干擾顯著減弱,因此熱流顯著降低。

    圖6(h)展示了第Ⅴ類激波干擾自持射流控制結果。高速來流經過試驗模型下方的激波系壓縮,并在集氣通道內滯止,形成高壓區(qū),從而可經過射流通道在射流出口形成自持射流,經測量,自持射流壓比約為1.7。在自持射流控制下,激波干擾區(qū)同樣被推向下游,此時熱流極值出現(xiàn)在X/L1=0.14(測點8),約為73.2 kW/m2,相比無控狀態(tài)降低了約20%,如圖9 所示。然而,對比發(fā)現(xiàn),壓比1.7 的自持射流的降熱效果低于壓比為1.2 時的氣源射流(降熱約37.9%)。這是由于氣源射流由高壓氣源直接產生,總溫為室溫,而自持射流的總溫則為高速來流滯止后的總溫,接近風洞來流總溫。由此可見射流溫度對于激波干擾控制降熱效果也具有較為顯著的影響。

    3.2 第Ⅳ類激波干擾控制

    30-75 模型中第Ⅳ類激波干擾控制的紋影結果及示意圖如圖10、圖11 所示。與第Ⅴ類激波干擾控制類似,隨著氣源射流的噴出,射流誘導產生一道弓形射流激波,射流和射流激波共同對下游激波干擾區(qū)產生控制作用,且射流壓比對流場結構及降熱效果影響同樣顯著。當射流壓比較小時(PR≤4.8,圖10(a)~圖10(c)、圖11(a)),第Ⅳ類激波干擾流場結構未發(fā)生明顯變化,此時,弓形射流激波JP 與一楔斜激波IP 以及二楔弓形激波PW 相交于P點,并相互干擾產生超聲速射流和膨脹波系PL。強射流使得分離區(qū)及分離激波SS1向上游移動,并與一楔斜激波相交于S1點。但是射流及剪切層JM 對超聲速射流與壁面的干擾有一定的隔絕作用,相應的,激波干擾區(qū)熱流極值也有一定的降低。射流壓比4.8時降為約87.4 kW/m2,相比基準狀態(tài)降低了約21.3%,熱流極值位置均在X/L1=0.12(測點8)。當射流壓比較大時(PR>4.8,圖10(d)~圖10(j)、圖11(b)),激波干擾流場開始出現(xiàn)明顯的變化,弓形射流激波JP1與一楔斜激波IP3以及二楔弓形激波P2W 通過透射激波P3P2和P3P1相連接,且X/L1=0.04~0.12(測點5~8)附近逐漸出現(xiàn)一道滯止激波P2Q。當射流壓比大于10.2 時(圖10(g)~圖10(j)),滯止激 波P2Q 逐漸與二楔弓形激波P2W 直接相連,且隨著射流壓比越來越大,滯止激波P2Q 也逐漸增強。滯止激波P2Q 的出現(xiàn)使得無控時激波干擾區(qū)產生的超聲速射流完全消失,超聲速流動經過滯止激波P2W 轉變?yōu)閬喡曀倭鲃?,從而避免了超聲速射流直接沖擊二楔壁面。相應的,雙楔壁面熱流也顯著降低,且熱流極值均出現(xiàn)在X/L1=0.07~0.12(測點6~8)之間,這也與滯止激波所在的位置密切相關。當射流壓比為17.4時,熱流極值僅為22.7 kW/m2,降低了79.6%,如圖12 所示。可見對于第Ⅳ類激波干擾,當射流壓比較大時,同樣是隔絕作用和結構改變作用共同導致了雙楔壁面熱流降低。

    圖10 第Ⅳ類激波干擾氣源/自持射流控制紋影測量結果Fig.10 Schlieren measurement results of Type-Ⅳ shock interaction control by steady/self-sustaining jets

    圖11 第Ⅳ類激波干擾射流控制示意圖Fig.11 Schematic diagrams of Tpe-Ⅳ shock interaction control by jets

    圖12 第Ⅳ類激波干擾氣源/自持射流控制熱流測量結果Fig.12 Heat flux measurement results of Type-Ⅳshock interaction control by steady/selfsustaining jets

    圖10(k)展示了第Ⅳ類激波干擾自持射流控制結果。與第Ⅴ類激波干擾自持射流控制一致,高速來流在集氣通道內滯止,并在射流出口形成壓比約為1.7 的自持射流。由于自持射流壓比較低,激波干擾流場未發(fā)生顯著改變。此時熱流極值仍出現(xiàn)在X/L1=0.12(測點8),約為106.1 kW/m2,僅降低了約4.5%。相同自持射流控制下降熱效果顯著低于第Ⅴ類激波干擾控制,可見第Ⅳ類激波干擾控制對于射流壓比的要求更高。在后續(xù)的研究中,可通過優(yōu)化進氣口設計減少壓力損失以及將進氣口布置于飛行器駐點等高壓處提高自持射流壓比,進而提升控制效果。

    4 結論

    開展了基于氣源/自持定常射流的高超聲速雙楔流場第Ⅴ類、第Ⅳ類激波干擾主動流動控制試驗研究,進行了紋影及熱流測量,得到的主要結論如下:

    1)雙楔激波干擾流場控制降熱機理主要體現(xiàn)在兩方面:一是隔絕作用,射流與來流干擾形成剪切層,減弱了激波干擾區(qū)與壁面的相互作用;二是結構改變作用,射流顯著改變了激波干擾結構,使其不再是傳統(tǒng)的第Ⅴ類、第Ⅳ類激波干擾,避免了局部強激波干擾區(qū)的產生。射流壓比較小時,隔絕作用占主導;射流壓比較大時,兩種因素共同起重要作用。

    2)射流壓比對激波干擾流場控制降熱效果具有顯著影響,射流壓比越大,隔絕作用及結構改變作用均越強,熱流極值降低比例也越大,第Ⅴ類、第Ⅳ類激波干擾的熱流極值最高降低分別約為81.2%和79.6%。

    3)提出了基于高速來流能量利用的雙楔激波干擾自持射流主動流動控制方式,在自持射流(壓比為1.7)控制下,第Ⅴ類、第Ⅳ類激波干擾熱流極值分別降低約20%和4.5%。提升集氣壓力是提高自持射流壓比,進而提升激波干擾控制降熱效果的關鍵。

    致謝

    衷心感謝南京航空航天大學NHW 風洞實驗室的同學和實驗員在試驗完成過程中提供的幫助,衷心感謝中國科學技術大學王軍在論文完成過程中提供的幫助。

    猜你喜歡
    弓形來流氣源
    兩種典型來流條件下風力機尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    降低弓形換位線圈換位尺寸偏差研究
    防爆電機(2020年4期)2020-12-14 03:11:16
    不同來流條件對溢洪道過流能力的影響
    飛機地面氣源機組設計及關鍵技術
    專用汽車(2018年10期)2018-11-02 05:32:30
    一次由“弓形回波”所引發(fā)的冰雹天氣過程雷達特征分析
    彈發(fā)匹配驗證試驗系統(tǒng)來流快速啟動技術研究
    豬弓形體病的診治
    大型飛機氣源起動裝備供氣流程設計與計算
    氣源凈化技術問答
    壓縮機技術(2014年4期)2014-03-20 15:55:40
    MOCVD氣源輸運技術研究
    機械與電子(2014年3期)2014-02-28 02:07:50
    真人做人爱边吃奶动态| 午夜激情av网站| 男女下面插进去视频免费观看| 国产欧美日韩一区二区三| 丰满饥渴人妻一区二区三| 亚洲自偷自拍图片 自拍| 国产亚洲午夜精品一区二区久久| 久久ye,这里只有精品| 久久人妻av系列| 大码成人一级视频| 精品久久蜜臀av无| 一区二区三区精品91| 老鸭窝网址在线观看| 国产福利在线免费观看视频| 青草久久国产| 日韩有码中文字幕| 久久天躁狠狠躁夜夜2o2o| 久久午夜亚洲精品久久| 99香蕉大伊视频| 美女扒开内裤让男人捅视频| 午夜免费鲁丝| 精品高清国产在线一区| 欧美日韩av久久| netflix在线观看网站| 亚洲午夜精品一区,二区,三区| 国产精品久久久久久精品电影小说| 久久精品国产综合久久久| 精品久久久精品久久久| 日日爽夜夜爽网站| 少妇 在线观看| 天天躁夜夜躁狠狠躁躁| 老司机在亚洲福利影院| 热99国产精品久久久久久7| 国产亚洲av高清不卡| 久久久精品区二区三区| 欧美+亚洲+日韩+国产| 啦啦啦免费观看视频1| 丁香六月天网| 纯流量卡能插随身wifi吗| 757午夜福利合集在线观看| 热re99久久国产66热| 在线观看免费视频日本深夜| 法律面前人人平等表现在哪些方面| 桃红色精品国产亚洲av| 亚洲人成伊人成综合网2020| 自拍欧美九色日韩亚洲蝌蚪91| 国产aⅴ精品一区二区三区波| 大码成人一级视频| 麻豆成人av在线观看| 精品一区二区三卡| 中国美女看黄片| 多毛熟女@视频| 变态另类成人亚洲欧美熟女 | 亚洲精品中文字幕一二三四区 | 国产免费av片在线观看野外av| 69av精品久久久久久 | 91精品三级在线观看| 欧美日韩成人在线一区二区| 国精品久久久久久国模美| 少妇粗大呻吟视频| 国产一区二区三区在线臀色熟女 | 久久九九热精品免费| 最新美女视频免费是黄的| 女同久久另类99精品国产91| 国产精品二区激情视频| 国产精品秋霞免费鲁丝片| 成人精品一区二区免费| 日本欧美视频一区| 视频在线观看一区二区三区| 国产av国产精品国产| 午夜精品国产一区二区电影| 亚洲九九香蕉| 国产精品久久电影中文字幕 | 丝袜在线中文字幕| 国产精品一区二区在线观看99| 精品一区二区三区四区五区乱码| 欧美精品啪啪一区二区三区| 国产一卡二卡三卡精品| 99香蕉大伊视频| 热re99久久国产66热| 国产男靠女视频免费网站| 欧美激情高清一区二区三区| √禁漫天堂资源中文www| 午夜激情av网站| 国产精品国产av在线观看| 免费看十八禁软件| 国精品久久久久久国模美| 成人免费观看视频高清| 91国产中文字幕| 色视频在线一区二区三区| 超碰97精品在线观看| 国产欧美日韩精品亚洲av| 操出白浆在线播放| 亚洲一区中文字幕在线| 宅男免费午夜| 国产国语露脸激情在线看| 免费在线观看日本一区| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 在线播放国产精品三级| 动漫黄色视频在线观看| 久久久国产精品麻豆| 亚洲国产av影院在线观看| av电影中文网址| 国产男靠女视频免费网站| 一夜夜www| 高清av免费在线| 97人妻天天添夜夜摸| 蜜桃在线观看..| 欧美精品一区二区免费开放| 咕卡用的链子| 两个人免费观看高清视频| 亚洲av国产av综合av卡| 久久 成人 亚洲| 久久精品国产a三级三级三级| 亚洲性夜色夜夜综合| 色综合欧美亚洲国产小说| 黄色视频在线播放观看不卡| 男女床上黄色一级片免费看| 女人精品久久久久毛片| 在线观看免费视频日本深夜| 国产一区二区三区在线臀色熟女 | 热99re8久久精品国产| 亚洲午夜理论影院| 精品人妻在线不人妻| 免费人妻精品一区二区三区视频| 精品卡一卡二卡四卡免费| 亚洲美女黄片视频| 最近最新免费中文字幕在线| 老司机午夜福利在线观看视频 | 好男人电影高清在线观看| av免费在线观看网站| 午夜免费成人在线视频| 亚洲男人天堂网一区| 69av精品久久久久久 | 欧美+亚洲+日韩+国产| www.999成人在线观看| 亚洲av日韩在线播放| 高清av免费在线| 曰老女人黄片| 在线观看一区二区三区激情| 国产av精品麻豆| 国产97色在线日韩免费| 男男h啪啪无遮挡| 日韩免费av在线播放| 国产午夜精品久久久久久| 久9热在线精品视频| 两个人看的免费小视频| 一级黄色大片毛片| 久热这里只有精品99| 国产免费福利视频在线观看| 日韩一卡2卡3卡4卡2021年| 国产99久久九九免费精品| 久久久精品区二区三区| 欧美变态另类bdsm刘玥| 日韩人妻精品一区2区三区| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 手机成人av网站| 涩涩av久久男人的天堂| 国产午夜精品久久久久久| 亚洲五月婷婷丁香| 正在播放国产对白刺激| 国产精品国产高清国产av | 国产精品香港三级国产av潘金莲| 成人黄色视频免费在线看| 老司机福利观看| 精品福利观看| 91成年电影在线观看| 国产亚洲欧美精品永久| 国产成人欧美| 久久亚洲真实| 亚洲av成人一区二区三| 亚洲av美国av| 欧美黄色淫秽网站| 不卡一级毛片| 久久99一区二区三区| 9色porny在线观看| 99re在线观看精品视频| 日本一区二区免费在线视频| 成年版毛片免费区| 91成年电影在线观看| www.精华液| 三上悠亚av全集在线观看| 高清黄色对白视频在线免费看| 国产精品一区二区免费欧美| 啦啦啦免费观看视频1| 久久精品亚洲av国产电影网| 欧美精品啪啪一区二区三区| videosex国产| 9色porny在线观看| 国产精品香港三级国产av潘金莲| 日本撒尿小便嘘嘘汇集6| 久久国产精品影院| 黄频高清免费视频| 青青草视频在线视频观看| 国产精品影院久久| 国产色视频综合| 久久精品国产综合久久久| 国产欧美日韩一区二区三| 久久这里只有精品19| 黄色成人免费大全| 久久国产精品大桥未久av| 国产高清videossex| av线在线观看网站| 久久九九热精品免费| 夜夜骑夜夜射夜夜干| 国产伦理片在线播放av一区| 精品久久蜜臀av无| 亚洲精品在线观看二区| 亚洲熟女精品中文字幕| 亚洲精品成人av观看孕妇| 亚洲av电影在线进入| 色尼玛亚洲综合影院| 丝瓜视频免费看黄片| 天天操日日干夜夜撸| 免费在线观看日本一区| 久久久欧美国产精品| 老司机午夜十八禁免费视频| 999久久久精品免费观看国产| 一夜夜www| 免费在线观看日本一区| 18禁国产床啪视频网站| 99久久精品国产亚洲精品| 国产亚洲精品一区二区www | a在线观看视频网站| 大片免费播放器 马上看| 一区二区日韩欧美中文字幕| 少妇的丰满在线观看| 蜜桃国产av成人99| 国产成人影院久久av| 成人永久免费在线观看视频 | 亚洲欧美一区二区三区久久| 久久久久精品人妻al黑| 男女高潮啪啪啪动态图| 男女午夜视频在线观看| 免费在线观看黄色视频的| 水蜜桃什么品种好| 欧美日韩中文字幕国产精品一区二区三区 | 欧美乱码精品一区二区三区| 久久久久久久国产电影| 十八禁网站网址无遮挡| 国产欧美日韩一区二区精品| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久精品电影小说| av欧美777| 极品教师在线免费播放| 国产国语露脸激情在线看| 久久精品国产99精品国产亚洲性色 | 热99re8久久精品国产| 91老司机精品| 国产av一区二区精品久久| 精品国产一区二区久久| 蜜桃国产av成人99| 麻豆av在线久日| 久久精品国产综合久久久| 欧美黄色片欧美黄色片| 久久久精品94久久精品| 亚洲国产看品久久| 日韩大片免费观看网站| 欧美+亚洲+日韩+国产| 国产成人影院久久av| av不卡在线播放| 国产午夜精品久久久久久| 久久久久精品国产欧美久久久| 精品亚洲成国产av| 精品国产一区二区久久| 日韩一区二区三区影片| 国产在线观看jvid| 久久久久久久精品吃奶| 亚洲精品国产一区二区精华液| 日韩免费高清中文字幕av| 亚洲国产毛片av蜜桃av| 久久亚洲真实| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品一二三| 免费人妻精品一区二区三区视频| 丁香六月欧美| 精品亚洲成国产av| av超薄肉色丝袜交足视频| 国产精品一区二区免费欧美| 老熟妇仑乱视频hdxx| 久久中文字幕人妻熟女| 免费久久久久久久精品成人欧美视频| 亚洲精品一卡2卡三卡4卡5卡| 美女高潮喷水抽搐中文字幕| 视频区图区小说| 久久久久网色| 亚洲成av片中文字幕在线观看| 色婷婷av一区二区三区视频| 成年人午夜在线观看视频| 99在线人妻在线中文字幕 | 黄色成人免费大全| 日韩人妻精品一区2区三区| 一本大道久久a久久精品| 女人久久www免费人成看片| 91大片在线观看| 在线观看免费视频网站a站| 成人国语在线视频| 新久久久久国产一级毛片| av福利片在线| 亚洲国产成人一精品久久久| 后天国语完整版免费观看| av免费在线观看网站| 男女免费视频国产| 超碰成人久久| 亚洲人成77777在线视频| 看免费av毛片| 啦啦啦 在线观看视频| 久久ye,这里只有精品| 精品视频人人做人人爽| 18禁裸乳无遮挡动漫免费视频| 亚洲av成人不卡在线观看播放网| 亚洲熟女精品中文字幕| 亚洲欧美激情在线| 日韩大码丰满熟妇| 国产精品免费视频内射| 久久精品91无色码中文字幕| 在线 av 中文字幕| 满18在线观看网站| 亚洲午夜精品一区,二区,三区| 肉色欧美久久久久久久蜜桃| 中文字幕高清在线视频| 日韩欧美三级三区| 国产国语露脸激情在线看| 俄罗斯特黄特色一大片| 亚洲国产av影院在线观看| 成人永久免费在线观看视频 | 久久精品亚洲av国产电影网| av不卡在线播放| 黄色视频不卡| 日韩大片免费观看网站| 国产一区二区在线观看av| 多毛熟女@视频| 999精品在线视频| 黄频高清免费视频| 日日夜夜操网爽| 午夜福利乱码中文字幕| 欧美日韩福利视频一区二区| 精品福利观看| 久久av网站| 黄色片一级片一级黄色片| 下体分泌物呈黄色| 国产精品二区激情视频| 欧美亚洲 丝袜 人妻 在线| 久久精品亚洲熟妇少妇任你| 十八禁网站免费在线| 亚洲免费av在线视频| 在线观看66精品国产| 免费在线观看完整版高清| 亚洲欧美激情在线| 亚洲精品中文字幕一二三四区 | 99re在线观看精品视频| 国产欧美日韩综合在线一区二区| 精品福利观看| 亚洲中文av在线| 极品少妇高潮喷水抽搐| 久久亚洲真实| 男男h啪啪无遮挡| 欧美午夜高清在线| 欧美人与性动交α欧美精品济南到| 国产精品一区二区在线不卡| 欧美日韩黄片免| 亚洲精品国产精品久久久不卡| 在线av久久热| 亚洲熟女毛片儿| 国产男女超爽视频在线观看| 黑人猛操日本美女一级片| 欧美+亚洲+日韩+国产| 纵有疾风起免费观看全集完整版| 看免费av毛片| 成年动漫av网址| 热99re8久久精品国产| 一二三四在线观看免费中文在| 久久久国产精品麻豆| 成人影院久久| 国产亚洲精品第一综合不卡| 久久精品亚洲av国产电影网| 亚洲自偷自拍图片 自拍| 国产精品二区激情视频| 高清av免费在线| 国产精品影院久久| 欧美精品一区二区免费开放| 久久精品亚洲熟妇少妇任你| 日本黄色日本黄色录像| 日日夜夜操网爽| 性高湖久久久久久久久免费观看| 精品久久久精品久久久| 满18在线观看网站| 又紧又爽又黄一区二区| 久久精品亚洲av国产电影网| 三上悠亚av全集在线观看| 国产成人一区二区三区免费视频网站| 在线观看人妻少妇| www.精华液| 蜜桃在线观看..| 久久人人爽av亚洲精品天堂| 亚洲中文av在线| 午夜日韩欧美国产| 亚洲精品中文字幕一二三四区 | 叶爱在线成人免费视频播放| 不卡av一区二区三区| cao死你这个sao货| 麻豆乱淫一区二区| 亚洲天堂av无毛| 久久久久久人人人人人| 老司机午夜十八禁免费视频| 国产xxxxx性猛交| √禁漫天堂资源中文www| 欧美午夜高清在线| 99精品在免费线老司机午夜| 99精国产麻豆久久婷婷| 免费日韩欧美在线观看| 色精品久久人妻99蜜桃| 亚洲精品在线美女| 99热国产这里只有精品6| 黑人巨大精品欧美一区二区蜜桃| 精品久久蜜臀av无| 亚洲视频免费观看视频| xxxhd国产人妻xxx| 91成年电影在线观看| 夜夜骑夜夜射夜夜干| 69av精品久久久久久 | 人人妻人人澡人人看| 女人高潮潮喷娇喘18禁视频| 高清av免费在线| 热99久久久久精品小说推荐| 亚洲色图综合在线观看| 欧美日韩亚洲综合一区二区三区_| www日本在线高清视频| 最近最新免费中文字幕在线| 美女视频免费永久观看网站| aaaaa片日本免费| 天天躁狠狠躁夜夜躁狠狠躁| √禁漫天堂资源中文www| 精品乱码久久久久久99久播| 中文字幕高清在线视频| 高清黄色对白视频在线免费看| 一级a爱视频在线免费观看| 国产色视频综合| 亚洲全国av大片| 色婷婷av一区二区三区视频| 午夜福利乱码中文字幕| 人人妻人人添人人爽欧美一区卜| 黄色视频,在线免费观看| 99国产精品一区二区三区| 中亚洲国语对白在线视频| 一区二区三区激情视频| 精品国产一区二区久久| 国产精品香港三级国产av潘金莲| 久久久久久久久免费视频了| 久久精品人人爽人人爽视色| 老司机亚洲免费影院| 汤姆久久久久久久影院中文字幕| 91字幕亚洲| 亚洲第一欧美日韩一区二区三区 | 在线av久久热| 久久天堂一区二区三区四区| 亚洲色图综合在线观看| 久久人妻熟女aⅴ| 国产一区二区 视频在线| av欧美777| 青青草视频在线视频观看| 91成年电影在线观看| 757午夜福利合集在线观看| 这个男人来自地球电影免费观看| 亚洲专区中文字幕在线| 色综合欧美亚洲国产小说| 宅男免费午夜| 色播在线永久视频| 久久久久久久久久久久大奶| 老司机靠b影院| 成人国产av品久久久| 久久久久视频综合| 一进一出好大好爽视频| 中文字幕高清在线视频| 亚洲黑人精品在线| 12—13女人毛片做爰片一| 国产区一区二久久| 91麻豆精品激情在线观看国产 | 丰满迷人的少妇在线观看| 女人久久www免费人成看片| 欧美性长视频在线观看| 一区二区三区精品91| av线在线观看网站| 免费av中文字幕在线| 91精品三级在线观看| 美国免费a级毛片| 中国美女看黄片| 18禁观看日本| 亚洲综合色网址| 精品乱码久久久久久99久播| 在线观看舔阴道视频| 国产精品久久久久久人妻精品电影 | 十八禁网站网址无遮挡| 欧美人与性动交α欧美软件| 这个男人来自地球电影免费观看| 国产深夜福利视频在线观看| 国产男女内射视频| 岛国毛片在线播放| 老熟妇乱子伦视频在线观看| 亚洲熟女精品中文字幕| 在线观看一区二区三区激情| 黄色a级毛片大全视频| 国产精品自产拍在线观看55亚洲 | 他把我摸到了高潮在线观看 | 亚洲黑人精品在线| 午夜视频精品福利| 国产亚洲一区二区精品| 亚洲国产av新网站| 考比视频在线观看| 精品少妇内射三级| 亚洲一区二区三区欧美精品| 亚洲av成人不卡在线观看播放网| 色尼玛亚洲综合影院| 女性生殖器流出的白浆| 国产精品98久久久久久宅男小说| 国产精品美女特级片免费视频播放器 | 99热国产这里只有精品6| 美女午夜性视频免费| 国产精品一区二区免费欧美| 狠狠婷婷综合久久久久久88av| 日本黄色视频三级网站网址 | 我要看黄色一级片免费的| 超色免费av| 国产野战对白在线观看| 最新美女视频免费是黄的| 亚洲av电影在线进入| 变态另类成人亚洲欧美熟女 | 啦啦啦免费观看视频1| 久久国产精品人妻蜜桃| 久久中文字幕一级| 在线观看www视频免费| 两人在一起打扑克的视频| 另类亚洲欧美激情| 亚洲成人手机| 久久天躁狠狠躁夜夜2o2o| 老熟女久久久| 欧美大码av| 国产有黄有色有爽视频| 欧美日韩黄片免| 大码成人一级视频| 啦啦啦 在线观看视频| 国产成+人综合+亚洲专区| 高清毛片免费观看视频网站 | 不卡av一区二区三区| tocl精华| 日韩精品免费视频一区二区三区| 少妇猛男粗大的猛烈进出视频| a级毛片在线看网站| 99久久人妻综合| 国产精品一区二区在线不卡| 91av网站免费观看| 免费在线观看影片大全网站| 免费久久久久久久精品成人欧美视频| 性高湖久久久久久久久免费观看| 少妇 在线观看| 日本vs欧美在线观看视频| 国产精品久久久人人做人人爽| 精品卡一卡二卡四卡免费| 99riav亚洲国产免费| 一个人免费看片子| 亚洲 欧美一区二区三区| 久久久久久久久久久久大奶| 亚洲精品在线观看二区| 9色porny在线观看| 免费女性裸体啪啪无遮挡网站| 妹子高潮喷水视频| 精品欧美一区二区三区在线| 啦啦啦 在线观看视频| 国产精品99久久99久久久不卡| 精品久久久久久电影网| 一本一本久久a久久精品综合妖精| 免费在线观看视频国产中文字幕亚洲| 美女主播在线视频| 国产av一区二区精品久久| 亚洲欧美一区二区三区久久| 成年人黄色毛片网站| 男女无遮挡免费网站观看| 亚洲成人手机| 成人手机av| 国产淫语在线视频| 9191精品国产免费久久| 两性夫妻黄色片| 久久精品熟女亚洲av麻豆精品| 丝袜美足系列| 久久性视频一级片| 精品视频人人做人人爽| netflix在线观看网站| 如日韩欧美国产精品一区二区三区| 啦啦啦免费观看视频1| 中文字幕高清在线视频| 在线观看免费日韩欧美大片| 精品国产一区二区三区四区第35| 色播在线永久视频| 欧美精品一区二区免费开放| h视频一区二区三区| 精品国产国语对白av| 如日韩欧美国产精品一区二区三区| 亚洲美女黄片视频| 青青草视频在线视频观看| 国产精品久久久久久人妻精品电影 | 欧美亚洲 丝袜 人妻 在线| 在线观看一区二区三区激情| 精品高清国产在线一区| 久久精品亚洲av国产电影网| 国产淫语在线视频| 在线观看舔阴道视频| 亚洲一卡2卡3卡4卡5卡精品中文| 他把我摸到了高潮在线观看 | 国产精品成人在线| 美女午夜性视频免费| 精品国产国语对白av| 国产欧美日韩一区二区精品| 久久毛片免费看一区二区三区|