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

    邊界層轉(zhuǎn)捩與壓縮拐角分離流動(dòng)的非定常作用

    2017-03-15 07:06:00解少飛楊學(xué)軍
    關(guān)鍵詞:拐角層流邊界層

    解少飛, 宮 建,*, 高 波, 楊學(xué)軍, 沈 清

    (1. 中國(guó)航天空氣動(dòng)力技術(shù)研究院, 北京 100074; 2. 北京宇航系統(tǒng)工程研究所, 北京 100076)

    邊界層轉(zhuǎn)捩與壓縮拐角分離流動(dòng)的非定常作用

    解少飛1, 宮 建1,*, 高 波2, 楊學(xué)軍2, 沈 清1

    (1. 中國(guó)航天空氣動(dòng)力技術(shù)研究院, 北京 100074; 2. 北京宇航系統(tǒng)工程研究所, 北京 100076)

    為了認(rèn)識(shí)激波與轉(zhuǎn)捩邊界層之間的相互作用,選取壓縮拐角模型為研究對(duì)象,在高超聲速風(fēng)洞中開(kāi)展了激波/轉(zhuǎn)捩邊界層干擾試驗(yàn)研究。試驗(yàn)在FD-20炮風(fēng)洞中進(jìn)行,試驗(yàn)馬赫數(shù)為8,雷諾數(shù)0.5×107~2×107/m。試驗(yàn)采用了薄膜電阻溫度計(jì)和常規(guī)壓力傳感器,分別測(cè)量壓縮拐角周圍的熱流和壓力分布。根據(jù)干擾區(qū)上下游的邊界層流態(tài),將試驗(yàn)分為兩部分:層流/湍流干擾和轉(zhuǎn)捩/湍流干擾。對(duì)比分析了邊界層轉(zhuǎn)捩發(fā)生在干擾區(qū)時(shí)的熱流分布、壓力分布以及脈動(dòng)熱流的非定常特性。研究結(jié)果表明,激波/轉(zhuǎn)捩邊界層干擾的熱流和壓力分布特征,不同于常規(guī)的層流干擾和湍流干擾,其介于二者之間。層流/湍流干擾的熱流和壓力分布特征類似于層流干擾;轉(zhuǎn)捩/湍流干擾的熱流和壓力分布特征類似于湍流干擾?;ルp譜分析結(jié)果表明,當(dāng)邊界層轉(zhuǎn)捩發(fā)生在分離區(qū)時(shí),轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)同時(shí)出現(xiàn)并增長(zhǎng)。當(dāng)兩者幅值足夠大時(shí),轉(zhuǎn)捩脈動(dòng)會(huì)與分離泡脈動(dòng)發(fā)生非線性耦合作用。這種耦合作用會(huì)誘導(dǎo)出新的頻率特征的脈動(dòng)結(jié)構(gòu),從而使分離區(qū)內(nèi)的脈動(dòng)量顯著增加。

    邊界層轉(zhuǎn)捩;壓縮拐角;轉(zhuǎn)捩脈動(dòng);分離泡脈動(dòng);非定常作用

    0 引 言

    激波/邊界層干擾是高超聲速飛行器上普遍存在的流動(dòng)現(xiàn)象。在逆壓梯度下誘導(dǎo)邊界層分離,分離區(qū)域出現(xiàn)分離激波、再附激波、激波振蕩等復(fù)雜流動(dòng)現(xiàn)象,飛行器的氣動(dòng)力/熱性能受到嚴(yán)重影響。高超聲速飛行器在飛行過(guò)程中,其激波分離區(qū)附近的邊界層或從層流過(guò)渡到湍流(如再入階段和巡航加速階段),或是從湍流過(guò)渡到層流(如爬升階段),轉(zhuǎn)捩邊界層不可避免的會(huì)與激波發(fā)生相互作用。多年來(lái),研究者對(duì)層流和湍流條件下的激波/邊界層干擾開(kāi)展了大量的實(shí)驗(yàn)和計(jì)算研究[1-4],然而激波/轉(zhuǎn)捩邊界層干擾問(wèn)題較少研究。盡管如此,Dolling[5-6]指出這類激波-邊界層干擾很可能發(fā)生在進(jìn)氣道入口和飛行器其他外部位置。因此,開(kāi)展這項(xiàng)研究在理論上和應(yīng)用上都具有很大的意義。

    低超聲速方面,Chapman等[7]實(shí)驗(yàn)研究了低馬赫數(shù)(Ma=0.4~3.6)下的激波/邊界層干擾現(xiàn)象的流動(dòng)規(guī)律,涉及壓縮拐角、前臺(tái)階、入射激波等多種幾何外形。研究發(fā)現(xiàn),相同來(lái)流條件下,層流和轉(zhuǎn)捩干擾時(shí)分離區(qū)比湍流情況大很多;雷諾數(shù)增加轉(zhuǎn)捩干擾的分離區(qū)長(zhǎng)度減小,但層流和湍流情況下變化不大。

    高超聲速方面,Heffner[8]在馬赫數(shù)5條件下研究了雷諾數(shù)對(duì)激波/轉(zhuǎn)捩邊界層干擾的影響,并指出雷諾數(shù)增加導(dǎo)致分離區(qū)減小。Kaufman等[9]在研究鈍翼干擾時(shí)發(fā)現(xiàn)來(lái)流邊界層為轉(zhuǎn)捩時(shí),對(duì)稱面分離區(qū)尺度隨雷諾數(shù)增加而減小。然而,Pate[10]則發(fā)現(xiàn)馬赫數(shù)5條件下的激波/轉(zhuǎn)捩邊界層干擾,雷諾數(shù)對(duì)分離區(qū)長(zhǎng)度影響不大。Schuelein[11]在馬赫數(shù)6條件下實(shí)驗(yàn)測(cè)量了入射激波中的轉(zhuǎn)捩干擾。研究發(fā)現(xiàn),相同雷諾數(shù)下轉(zhuǎn)捩干擾的峰值熱流高于湍流干擾,激波入射位置的來(lái)流邊界層湍流度越高,分離泡越小。Benay等[12]在軸對(duì)稱外形(中空?qǐng)A柱與裙的組合)上也發(fā)現(xiàn)了轉(zhuǎn)捩干擾下熱流高于湍流干擾的現(xiàn)象。Schrijer和Scarano[13]對(duì)這一現(xiàn)象進(jìn)行了解釋,認(rèn)為這是轉(zhuǎn)捩出現(xiàn)在再附點(diǎn)R附近導(dǎo)致的。上述研究雖涉及了激波/邊界層干擾中的轉(zhuǎn)捩問(wèn)題,但僅局限于雷諾數(shù)效應(yīng)、熱流分布等現(xiàn)象研究,對(duì)其非定常特性的研究尚屬空白。

    國(guó)內(nèi)李素循研究團(tuán)隊(duì)[14-15]也對(duì)激波/邊界層干擾現(xiàn)象開(kāi)展了系統(tǒng)的實(shí)驗(yàn)和計(jì)算研究,包括不同的邊界層流態(tài)(層流和湍流),不同外形(鈍舵,壓縮拐角和噴流等)以及不同的實(shí)驗(yàn)狀態(tài)(馬赫數(shù)和雷諾數(shù)等),但對(duì)激波/邊界層干擾中的轉(zhuǎn)捩問(wèn)題,尚未涉及。

    本文在FD-20風(fēng)洞中,以楔角15°的壓縮拐角為研究對(duì)象,利用薄膜電阻溫度計(jì)和壓力傳感器,開(kāi)展激波/邊界層干擾與邊界層轉(zhuǎn)捩的研究,取得邊界層轉(zhuǎn)捩發(fā)生在分離區(qū)周圍時(shí)的流場(chǎng)結(jié)構(gòu)特征,探索邊界層轉(zhuǎn)捩與壓縮拐角分離流動(dòng)之間的非定常作用。

    1 實(shí)驗(yàn)設(shè)備、模型與狀態(tài)

    1.1 風(fēng) 洞

    實(shí)驗(yàn)在中國(guó)航天空氣動(dòng)力技術(shù)研究院的FD-20炮風(fēng)洞中進(jìn)行,如圖 1所示。這是一座高超聲速脈沖型風(fēng)洞。驅(qū)動(dòng)段長(zhǎng)10 m,內(nèi)徑160 mm,被驅(qū)動(dòng)段長(zhǎng)15 m,內(nèi)徑130 mm,噴管出口直徑為480 mm,試驗(yàn)段為直徑1.6 m×2 m。該風(fēng)洞可提供馬赫數(shù)為5~15的均勻流場(chǎng),相應(yīng)的Reynolds數(shù)模擬范圍為Re=1.0×106~6.0×107(1/m)。試驗(yàn)氣流總溫可達(dá)1200 K,根據(jù)不同高低壓匹配運(yùn)行方式,炮風(fēng)洞運(yùn)行時(shí)間為15~60ms,激波風(fēng)洞運(yùn)行時(shí)間為1~5 ms。

    1.2 實(shí)驗(yàn)?zāi)P团c測(cè)點(diǎn)位置

    圖2給出了模型上測(cè)點(diǎn)位置示意圖。平板長(zhǎng)500 mm,寬400 mm;楔板長(zhǎng)200 mm,寬400 mm,楔角15°。楔板前緣(即拐角位置)到平板前緣的距離為300 mm。模型表面沿流向安裝兩列測(cè)點(diǎn),分別為L(zhǎng)ine-1和Line-2。Line-1上安裝薄膜電阻溫度計(jì),用于表面熱流測(cè)量;Line-2上安裝壓力傳感器,用于表面壓力測(cè)量。為了盡可能精確的捕捉激波分離點(diǎn)S和再附點(diǎn)R的位置,在Line-1上使用了整體式傳感器。整體式傳感器可以使測(cè)點(diǎn)的間距更小,本次實(shí)驗(yàn)采用的間距為5 mm。其余位置相鄰測(cè)點(diǎn)間距為10 mm。

    1.3 實(shí)驗(yàn)狀態(tài)

    流場(chǎng)實(shí)際校測(cè)自由來(lái)流馬赫數(shù)為Ma=7.97。驅(qū)動(dòng)段壓力為6 MPa~30 MPa,實(shí)現(xiàn)來(lái)流單位雷諾數(shù)Re=5×106~1.8 ×107/m。風(fēng)洞有效運(yùn)行時(shí)間約為20 ms。詳細(xì)實(shí)驗(yàn)狀態(tài)如表1所示。

    表1 實(shí)驗(yàn)狀態(tài)Table 1 Experimental condition

    2 實(shí)驗(yàn)結(jié)果

    激波/轉(zhuǎn)捩邊界層干擾,可能存在兩種情況:1) 初始邊界層為層流(即分離點(diǎn)S上游為層流),轉(zhuǎn)捩發(fā)生在分離區(qū)內(nèi);2) 分離點(diǎn)S上游為轉(zhuǎn)捩邊界層。因此,按照分離區(qū)上下游的邊界層流態(tài),將實(shí)驗(yàn)中的激波/邊界層干擾分為兩類:層流/湍流干擾和轉(zhuǎn)捩/湍流干擾,并對(duì)這兩種激波/邊界層干擾分別討論分析。

    2.1 熱流和壓力分布

    兩類干擾的流場(chǎng)結(jié)構(gòu)無(wú)較大差異,均由分離激波、再附激波和分離泡等組成,如圖3所示。主要的差異在于熱流和壓力分布。

    圖4顯示了層流/湍流干擾情況下的熱流和壓力分布規(guī)律。Re=5.12×106/m時(shí),壓力在分離點(diǎn)S(x=-0.075 m)前緩慢上升,分離點(diǎn)S下游出現(xiàn)一個(gè)較長(zhǎng)的壓力平臺(tái),平臺(tái)壓力大約為未擾動(dòng)靜壓的2倍。在拐角上游,壓力再次上升,在x=0.05 m(P點(diǎn))附近達(dá)到峰值,之后逐漸下降為靜壓9.2倍時(shí)趨于常數(shù)。此時(shí)與二維無(wú)黏流斜激波理論計(jì)算得到的激波前后壓力比p2/p1=9.25基本吻合。熱流在分離前緩慢下降,與基于參考溫度法[16-17]計(jì)算的平板層流理論解吻合。分離點(diǎn)S下游,受到分離泡影響,熱流迅速降低,在分離區(qū)熱流達(dá)到最低點(diǎn)后開(kāi)始上升,形成明顯的“凹”型分布。之后拐角附近熱流迅速上升并達(dá)到峰值,熱流峰值位置與壓力峰值基本一致。雷諾數(shù)增加到Re=6.56×106/m時(shí),參數(shù)分布形態(tài)基本不變,但分離區(qū)長(zhǎng)度減小。

    圖5顯示了轉(zhuǎn)捩/湍流干擾情況下的熱流和壓力分布規(guī)律。Re=1.40×107/m時(shí),分離區(qū)內(nèi)的熱流凹陷減小,這是轉(zhuǎn)捩或湍流剪切層中脈動(dòng)結(jié)構(gòu)導(dǎo)致的熱流上升和分離泡導(dǎo)致的熱流下降相互作用的結(jié)果,此時(shí)壓力平臺(tái)也減?。籖e=1.75×107/m時(shí),轉(zhuǎn)捩在分離點(diǎn)S下游附近迅速完成,使得湍流剪切層在分離區(qū)內(nèi)起主導(dǎo)作用,湍流產(chǎn)生的高熱流導(dǎo)致熱流凹陷幾乎消失。

    總之,層流/湍流干擾時(shí),分離區(qū)內(nèi)的剪切層流態(tài)以層流和轉(zhuǎn)捩為主,熱流和壓力受到分離泡的影響較大,熱流和壓力分布類似于純層流干擾的情況,呈現(xiàn)出了熱流凹陷和壓力平臺(tái)的形態(tài);轉(zhuǎn)捩/湍流干擾時(shí),分離點(diǎn)S下游附近轉(zhuǎn)捩迅速完成,湍流是分離區(qū)及其下游的主要流態(tài),此時(shí)分離泡較小,熱流凹陷和壓力平臺(tái)受湍流結(jié)構(gòu)影響逐漸消失。

    2.2 非定常作用分析

    在高超聲速激波/邊界層干擾過(guò)程中,存在三種擾動(dòng)結(jié)構(gòu):1) 邊界層內(nèi)的擾動(dòng),如第二模態(tài)不穩(wěn)定波和G?rtler渦,其頻率最高大約為幾十千到幾百千赫茲量級(jí);2) 分離區(qū)內(nèi)的擾動(dòng),主要由不穩(wěn)定的分離泡產(chǎn)生,其頻率大約為幾千到十幾千赫茲量級(jí);3) 激波振蕩誘導(dǎo)的擾動(dòng),一般出現(xiàn)在湍流干擾的分離點(diǎn)S和再附點(diǎn)R附近,其頻率大約為幾百到一千赫茲量級(jí)。由于三種擾動(dòng)的頻率跨度較大,很難用一種傳感器同時(shí)捕捉,因此,本部分僅考慮25 kHz以下的擾動(dòng)量變化,分析這些擾動(dòng)的流向分布、主導(dǎo)頻率以及相互作用。傳感器采用薄膜電阻溫度計(jì),選用50 kHz采樣頻率。需要指出的是,由于激波振蕩擾動(dòng)出現(xiàn)的區(qū)域很小,大約僅為幾倍的邊界層厚度量級(jí),本文實(shí)驗(yàn)未捕捉到該現(xiàn)象。

    · 熱流脈動(dòng)量分布

    熱流的均方根值是反映熱流脈動(dòng)程度最直觀的參數(shù)之一。本文利用熱流脈動(dòng)與平均熱流值的比值的均方根值分析不同類型干擾下的熱流脈動(dòng)量分布。

    圖6顯示了不同類型干擾下熱流脈動(dòng)量分布與未擾動(dòng)平板熱流脈動(dòng)的比較結(jié)果。圖中,B表示無(wú)擾動(dòng)平板邊界層的轉(zhuǎn)捩起始位置,E表示轉(zhuǎn)捩終止位置,x=0為拐角位置。在無(wú)擾動(dòng)平板上,熱流脈動(dòng)的變化僅受邊界層流態(tài)的影響。層流邊界層中,熱流脈動(dòng)較低,并隨擾動(dòng)波增長(zhǎng),熱流脈動(dòng)小幅增加;轉(zhuǎn)捩邊界層中,擾動(dòng)波快速增加,熱流脈動(dòng)相應(yīng)增加;隨著擾動(dòng)波的破碎,脈動(dòng)熱流降低,并逐漸進(jìn)入湍流邊界層。圖中可以看出,轉(zhuǎn)捩/湍流干擾時(shí),轉(zhuǎn)捩在分離點(diǎn)S上游發(fā)生,在分離區(qū)內(nèi)結(jié)束。在分離區(qū)內(nèi),轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)發(fā)生耦合,誘導(dǎo)出更高幅值的熱流脈動(dòng),因此,熱流脈動(dòng)量峰值稍高于未擾動(dòng)平板的脈動(dòng)峰值。層流/湍流干擾時(shí),轉(zhuǎn)捩完全發(fā)生在分離區(qū)內(nèi)。轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)的耦合現(xiàn)象更加明顯,熱流脈動(dòng)量顯著增加,熱流脈動(dòng)量的增長(zhǎng)速度也高于未擾動(dòng)平板邊界層??傮w來(lái)說(shuō),相對(duì)于轉(zhuǎn)捩/湍流干擾,層流/湍流干擾會(huì)引起更高的脈動(dòng)熱流峰值。

    ·分離區(qū)內(nèi)脈動(dòng)結(jié)構(gòu)的非線性作用

    當(dāng)轉(zhuǎn)捩在分離區(qū)內(nèi)發(fā)生時(shí),轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)在分離區(qū)內(nèi)同時(shí)出現(xiàn)。它們是否會(huì)發(fā)生非線性相互作用,以哪種形式相互作用。本部分將利用互雙譜分析方法[18]對(duì)這一問(wèn)題展開(kāi)討論。

    互雙譜分析方法,首先利用連續(xù)小波變換計(jì)算時(shí)間序列信號(hào)中不同頻率下的小波系數(shù),然后通過(guò)計(jì)算小波系數(shù)之間的相關(guān)系數(shù),最終得到不同頻率之間的耦合關(guān)系。其基本公式為:

    Cbio=

    其中,W(f,t)、W(f1,t)和W(f2,t)分別表示時(shí)間序列信號(hào)中不同尺度對(duì)應(yīng)的小波系數(shù)。Cbio表示不同頻率對(duì)應(yīng)的小波系數(shù)的相關(guān)系數(shù)。如果頻率為f、f1和f2的特征信號(hào)相互間是獨(dú)立的,那么相關(guān)系數(shù)近似為0;若相關(guān)系數(shù)較高(一般認(rèn)為大于0.5),則說(shuō)明這三個(gè)信號(hào)的耦合程度較高;若f1和f2完美耦合,其相關(guān)系數(shù)會(huì)是1,但這僅僅是理想的情況。因此,互雙譜系數(shù)值在0到1之間。

    圖7和圖8顯示了層流/湍流干擾(Re=5.12×106/m)時(shí),不同位置擾動(dòng)波的非線性相互作用云圖。層流邊界層在激波的干擾下,在x=-0.075 m處發(fā)生分離,但分離區(qū)內(nèi)流動(dòng)仍保持層流。在分離泡擾動(dòng)的影響下,x=-0.032 m~-0.022 m附近熱流開(kāi)始增長(zhǎng),轉(zhuǎn)捩開(kāi)始。轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)的非線性作用首先出現(xiàn)在x=-0.022 m處。擾動(dòng)沿流向發(fā)展,非線性作用逐漸增長(zhǎng),在x=-0.017 m~-0.012 m附近達(dá)到最強(qiáng)。之后逐漸減弱,在x=-0.007m下游消失。在整個(gè)非線性作用過(guò)程中,分離泡脈動(dòng)的頻率在0~5 kHz之間,轉(zhuǎn)捩脈動(dòng)則在5 kHz~20 kHz左右。

    圖9和圖10顯示了轉(zhuǎn)捩/湍流干擾(Re=1.75×107/m)時(shí),不同位置擾動(dòng)波的非線性相互作用云圖。x=-0.065 m位于分離點(diǎn)S上游的轉(zhuǎn)捩邊界層中,未出現(xiàn)擾動(dòng)的非線性作用。在分離點(diǎn)S附近x=-0.032 m處,轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)發(fā)生非線性作用,主要位于f1=3~8 kHz,f2=15~23 kHz范圍內(nèi),屬于和頻作用。由于激波/邊界層干擾發(fā)生在邊界層轉(zhuǎn)捩后期,邊界層內(nèi)的擾動(dòng)已經(jīng)開(kāi)始衰減,非線性作用也沿流向逐漸降低。在x=-0.022 m處下游消失。需要指出的是,x=-0.027 m處的低頻區(qū)域(f1=-2~2 kHz,f2=2~4 kHz)出現(xiàn)的非線性作用是不同頻段的分離泡脈動(dòng)相互耦合造成的。

    3 結(jié) 論

    本文利用試驗(yàn)手段對(duì)壓縮拐角誘導(dǎo)的激波/邊界層干擾開(kāi)展研究,對(duì)比了層流/湍流干擾和轉(zhuǎn)捩/湍流干擾的現(xiàn)象差異,發(fā)現(xiàn)上游邊界層流態(tài)會(huì)影響熱流和壓力分布特征。對(duì)比兩種干擾與未擾動(dòng)平板的脈動(dòng)熱流分布,發(fā)現(xiàn)層流/湍流干擾會(huì)導(dǎo)致比轉(zhuǎn)捩/湍流干擾更高的脈動(dòng)峰值。利用互雙譜分析方法,從頻域角度分析了轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)的非耦合作用。研究發(fā)現(xiàn),當(dāng)轉(zhuǎn)捩發(fā)生在分離區(qū)時(shí),轉(zhuǎn)捩脈動(dòng)與分離泡脈動(dòng)的相互作用,會(huì)誘導(dǎo)出更高的脈動(dòng)能量。這會(huì)嚴(yán)重影響飛行器舵面、進(jìn)氣道唇口等部位的結(jié)構(gòu)強(qiáng)度。

    本文雖對(duì)激波/邊界層干擾問(wèn)題得到了新的認(rèn)識(shí),但仍留有遺憾。高超聲速邊界層中主要的脈動(dòng)結(jié)構(gòu)頻率很高(約100 kHz以上),而本文僅對(duì)25 kHz以下的脈動(dòng)結(jié)構(gòu)開(kāi)展了研究。另外,傳感器間距較大,未捕捉到激波震蕩。上述兩點(diǎn),將是后期研究的重點(diǎn)。

    [1]Holden M S. Historical review of experimental studies and prediction methods to describe laminar and turbulent shock wave/boundary layer interactions in hypersonic flows[R]. AIAA 2006-494.

    [2]Holden M S, Wadhams T P. Code validation study of laminar shock-boundary layer and shock-shock interactions in hypersonic flow, part a:experimental measurement[R]. AIAA 2001-1031.

    [3]Clemens N T, Narayanaswamy V. Shock/turbulent boundary layer interactions: review of recent work on sources of unsteadiness[R]. AIAA 2009-3710[4]Rizzetta D P, Visbal M R. Large Eddy Simulation of supersonic compression ramp flows[R]. AIAA 2001-2858.

    [5]Dolling D S. 50 years of shock wave/boundary layer interaction research-What next?[R] AIAA 2000-2596.

    [6]Erengil M E, Dolling D S. Unsteady wave structure near separation in a mach 5 compression ramp interaction[J]. AIAA Journal, 1991, 29: 728-735.

    [7]Chapman D R, Kuehn D M, Larson H K. Investigation of separated flows in supersonic and subsonic streams with emphasis on the effect of transition[R]. NACA Report 1356, 1958.

    [8]Heffner K S, Chpoun A, Lengrand J C. Experimental study of transitional axisymmetric shock-boundary layer interactions at Mach 5[R]. AIAA 1993-3131.

    [9]Kaufman L, Korkegi R, Morton L. Shock impingement caused by boundary layer separation ahead of blunt fins[R]. Technical Report ARL 72-0118, Aerospace Research Labs, Wright-Patterson Air Force Base, 1972.

    [10]Pate S R. Investigation of flow separation on a two-dimensional flat plate having a variable-span trailing-edge flap atMa=3 and 5[R]. AEDC-TDR-64-14, AD0432831, 1964.

    [11]Schuelein E. Effects of laminar-turbulent transition on the shock-wave/ boundary-layer interaction[R]. AIAA 2014-3332.

    [12]Benay R, Chanetz B, Mangin B, et al. Shock wave/transitional boundary-layer interactions in hypersonic flow[J]. AIAA Journal, 2006, 44(6): 1243-1254.

    [13]Schrijer F F J, Scarano F. Experiments on hypersonic boundary layer separation and reattachment on a blunt cone-flare using quantitative infrared thermography[R]. AIAA 2003-6967.

    [14]Li S X. Complex flow reduced by shock and boundary layer[M]. Beijing: Science Press, 2007. (in Chinese)李素循. 激波與邊界層主導(dǎo)的復(fù)雜流動(dòng)[M]. 北京: 科學(xué)出版社, 2007.

    [15]Li S X, Ni Z Y. Investigation of Laminar interactive flowfield in hypersonic flow[J]. Journal of Astronautics, 2003, 24(6): 547-551. (in Chinese)李素循, 倪招勇. 高超聲速層流干擾流場(chǎng)研究[J]. 宇航學(xué)報(bào), 2003, 24(6): 547-551.

    [16]Eckert E K G.Engineering relation for friction and heat transfer to surface in high velocity flow[J]. Journal of Aeronautical Sciences, 1955, 22(4): 231-238.

    [17]Eckert E K G.Engineering relation for friction and heat transfer to surface in high velocity laminar and turbulent boundary layer flow over surfaces with constant pressure and temperature[J]. American Society of Mechanical Engineers, 1956, 78: 1273-1283.

    [18]Chokani N. Nonlinear evolution of Mack modes in a hypersonic boundary layer[J]. Physics of fluids, 2005,17: 014102.

    Unsteady interaction between transitional boundary layer and flow separation in compression corner

    Xie Shaofei1, Gong Jian1,*, Gao Bo2, Yang Xuejun2, Shen Qing1

    (1.ChinaAcademyofAerospaceAerodynamics,Beijing100074,China; 2.ChinaAstronautResearchandTrainingCenter,Beijing100076,China)

    The experiments of shock wave/transitional boundary layer interaction were performed in FD-20 wind tunnel to study the interactions between shock wave and transitional boundary layer. The experimental model was compression corner. The Mach number is 8, and the Reynolds number 0.5×107~2×107/m. Film resistance thermometers and pressure sensors were used in the experiments to measure the distributions of the heat flux and the pressure. According to the state of boundary layer in the upstream and downstream of the separation region, the experiment was divided into laminar/turbulence interference part and transition/turbulence interference part. The distributions of parameters and unsteady characteristics of heat flux fluctuations were compared in separation region where transition occurred. The experimental results show that the characteristics of the distributions of transitional interaction are different from those of the turbulence interaction and the laminar interaction. Cross-bispectrum analysis shows that the transitional fluctuations and the separation bubble fluctuations appear and grow simultaneously with the boundary layer transition occurring in separation region. The transitional fluctuations may interact with the separation bubble fluctuations nonlinearly, when the amplitudes of the both fluctuations are large enough. The fluctuation structures with new frequency characteristics are induced, leading to remarkable increment of the amplitudes of the fluctuations.

    boundary layer transition; compression corner; transitional fluctuations; separation bubble fluctuations; unsteady interaction

    0258-1825(2017)01-0129-07

    2015-07-21;

    2015-11-19

    國(guó)家自然科學(xué)基金項(xiàng)目(基金編號(hào):11372296);武器裝備預(yù)研基金項(xiàng)目(9140A13040413HT71001)

    解少飛(1986-),男,河北,博士,工程師,研究方向:高超聲速?gòu)?fù)雜流動(dòng)實(shí)驗(yàn)及實(shí)驗(yàn)技術(shù). E-mail: flying_1125@126.com

    宮建*(1976-),男,遼寧,碩士,研究員,研究方向:高超聲速風(fēng)洞試驗(yàn)技術(shù). E-mail: 13641040905@139.com

    解少飛, 宮建, 高波, 等. 邊界層轉(zhuǎn)捩與壓縮拐角分離流動(dòng)的非定常作用[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2017, 35(1): 129-135.

    10.7638/kqdlxxb-2015.0104 Xie S F, Gong J, Gao B, et al.Unsteady interaction between transitional boundary layer and flow separation in compression corner[J]. Acta Aerodynamica Sinica, 2017, 35(1): 129-135.

    V211.3

    A doi: 10.7638/kqdlxxb-2015.0104

    猜你喜歡
    拐角層流邊界層
    拐 角
    層流輥道電機(jī)IP56防護(hù)等級(jí)結(jié)構(gòu)設(shè)計(jì)
    摻氫對(duì)二甲醚層流燃燒特性的影響
    Where Is My Home?
    層流切應(yīng)力誘導(dǎo)microRNA-101下調(diào)EZH2抑制血管新生
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    走過(guò)那一個(gè)拐角
    美文(2017年4期)2017-02-23 14:26:12
    拐角遇到奇跡
    一類具有邊界層性質(zhì)的二次奇攝動(dòng)邊值問(wèn)題
    非特征邊界的MHD方程的邊界層
    亚洲国产欧美在线一区| av黄色大香蕉| 综合色丁香网| 老女人水多毛片| www日本在线高清视频| 麻豆精品久久久久久蜜桃| 91成人精品电影| 国产高清三级在线| 国产又爽黄色视频| 纯流量卡能插随身wifi吗| 99久久人妻综合| 国产一级毛片在线| 国产精品蜜桃在线观看| 国产一级毛片在线| 黄网站色视频无遮挡免费观看| 久久久久国产精品人妻一区二区| 在现免费观看毛片| 中文字幕av电影在线播放| 中文字幕最新亚洲高清| 制服人妻中文乱码| 爱豆传媒免费全集在线观看| 最新的欧美精品一区二区| 欧美精品一区二区大全| 亚洲av电影在线进入| 国产免费一级a男人的天堂| 免费女性裸体啪啪无遮挡网站| 熟女av电影| 国产成人一区二区在线| 亚洲国产毛片av蜜桃av| 黄网站色视频无遮挡免费观看| 在线免费观看不下载黄p国产| 久久国内精品自在自线图片| av不卡在线播放| 丰满少妇做爰视频| 日韩欧美精品免费久久| 国产男人的电影天堂91| 亚洲内射少妇av| 26uuu在线亚洲综合色| 国产欧美日韩综合在线一区二区| 亚洲av男天堂| 国产在线免费精品| 蜜臀久久99精品久久宅男| 精品久久蜜臀av无| 青青草视频在线视频观看| 大香蕉97超碰在线| 超色免费av| 2022亚洲国产成人精品| 蜜桃在线观看..| 亚洲精品456在线播放app| 97人妻天天添夜夜摸| 热re99久久国产66热| 久久热在线av| 9色porny在线观看| 波多野结衣一区麻豆| 18在线观看网站| 久久99精品国语久久久| 永久免费av网站大全| 欧美+日韩+精品| 麻豆精品久久久久久蜜桃| 久久人人爽人人片av| 婷婷色综合www| 久久综合国产亚洲精品| av片东京热男人的天堂| 国国产精品蜜臀av免费| 国产国拍精品亚洲av在线观看| av福利片在线| 91国产中文字幕| 国产成人精品久久久久久| 人成视频在线观看免费观看| 少妇高潮的动态图| 亚洲成av片中文字幕在线观看 | 欧美另类一区| 亚洲天堂av无毛| 少妇高潮的动态图| 亚洲久久久国产精品| 亚洲精华国产精华液的使用体验| 国产欧美亚洲国产| 精品亚洲成国产av| 精品人妻在线不人妻| 国产深夜福利视频在线观看| 超色免费av| 满18在线观看网站| 大香蕉久久网| 97超碰精品成人国产| 午夜视频国产福利| 欧美日本中文国产一区发布| 18禁国产床啪视频网站| 亚洲,欧美,日韩| 日本-黄色视频高清免费观看| 人人妻人人澡人人爽人人夜夜| 边亲边吃奶的免费视频| 国产成人精品无人区| 99视频精品全部免费 在线| 中文字幕另类日韩欧美亚洲嫩草| 在线天堂最新版资源| 亚洲欧美精品自产自拍| 18在线观看网站| 欧美日韩视频精品一区| 各种免费的搞黄视频| av播播在线观看一区| 天堂中文最新版在线下载| 18+在线观看网站| 国产亚洲一区二区精品| 欧美97在线视频| 美女中出高潮动态图| 亚洲成色77777| 欧美精品一区二区大全| 大香蕉久久成人网| www.色视频.com| 在线观看免费高清a一片| 伊人久久国产一区二区| 亚洲av国产av综合av卡| 午夜福利视频精品| 青春草视频在线免费观看| 伊人亚洲综合成人网| 91久久精品国产一区二区三区| a级毛色黄片| 亚洲国产日韩一区二区| av又黄又爽大尺度在线免费看| 精品一区二区三区视频在线| 国产一区二区三区综合在线观看 | 乱码一卡2卡4卡精品| 寂寞人妻少妇视频99o| 97在线视频观看| 国产亚洲一区二区精品| 一级毛片黄色毛片免费观看视频| 成人影院久久| 老熟女久久久| 婷婷色麻豆天堂久久| 黄片播放在线免费| 91成人精品电影| 天天操日日干夜夜撸| 丰满迷人的少妇在线观看| 欧美日韩亚洲高清精品| 寂寞人妻少妇视频99o| 男女国产视频网站| 久久精品国产鲁丝片午夜精品| av.在线天堂| 激情五月婷婷亚洲| 天天躁夜夜躁狠狠躁躁| 国产精品 国内视频| 制服人妻中文乱码| 新久久久久国产一级毛片| 欧美精品一区二区大全| 国产成人av激情在线播放| 精品99又大又爽又粗少妇毛片| 精品国产一区二区三区久久久樱花| 久久久国产一区二区| 欧美成人午夜精品| 成人亚洲欧美一区二区av| 少妇的逼好多水| 建设人人有责人人尽责人人享有的| 岛国毛片在线播放| 夫妻午夜视频| 中文字幕制服av| 全区人妻精品视频| 看免费av毛片| 国产日韩欧美在线精品| av福利片在线| 久久精品夜色国产| 在线天堂中文资源库| 18禁裸乳无遮挡动漫免费视频| 午夜激情av网站| 免费看光身美女| 女人精品久久久久毛片| 欧美3d第一页| www.av在线官网国产| 十分钟在线观看高清视频www| 亚洲国产精品成人久久小说| av有码第一页| 制服人妻中文乱码| 国产日韩欧美在线精品| 日韩精品有码人妻一区| 五月天丁香电影| 中国美白少妇内射xxxbb| 精品久久久久久电影网| 97精品久久久久久久久久精品| 精品少妇黑人巨大在线播放| 精品久久国产蜜桃| av黄色大香蕉| 日本午夜av视频| 99久久精品国产国产毛片| 久热这里只有精品99| 免费在线观看黄色视频的| 国产精品嫩草影院av在线观看| 日韩一区二区视频免费看| 亚洲熟女精品中文字幕| 9热在线视频观看99| 欧美人与善性xxx| 亚洲一级一片aⅴ在线观看| 久久人人爽人人片av| 老司机影院毛片| 男女边吃奶边做爰视频| videos熟女内射| 亚洲激情五月婷婷啪啪| 99久久精品国产国产毛片| 亚洲,欧美精品.| 精品一区二区三区四区五区乱码 | 狠狠婷婷综合久久久久久88av| 欧美精品av麻豆av| av在线老鸭窝| 精品国产乱码久久久久久小说| www日本在线高清视频| 国产探花极品一区二区| 中国三级夫妇交换| 熟妇人妻不卡中文字幕| 蜜臀久久99精品久久宅男| 成人无遮挡网站| 另类亚洲欧美激情| 色吧在线观看| 国产精品久久久久久久久免| 国产黄色免费在线视频| 老司机影院成人| 在线精品无人区一区二区三| 免费av中文字幕在线| 久久精品久久久久久久性| 人人妻人人澡人人爽人人夜夜| 熟女电影av网| 天堂8中文在线网| 久久久国产欧美日韩av| 最黄视频免费看| 美女国产高潮福利片在线看| 夫妻性生交免费视频一级片| 国产成人av激情在线播放| 天天影视国产精品| 亚洲综合色惰| 成年av动漫网址| 美女大奶头黄色视频| 日本爱情动作片www.在线观看| 一边摸一边做爽爽视频免费| 插逼视频在线观看| 国产日韩一区二区三区精品不卡| 国产色婷婷99| 亚洲综合色惰| 国产一级毛片在线| 欧美亚洲 丝袜 人妻 在线| 午夜91福利影院| 亚洲精华国产精华液的使用体验| 日韩三级伦理在线观看| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 丰满少妇做爰视频| 日韩成人伦理影院| 国产精品欧美亚洲77777| 亚洲美女搞黄在线观看| 成人毛片60女人毛片免费| 卡戴珊不雅视频在线播放| 极品少妇高潮喷水抽搐| 18在线观看网站| 国产精品久久久久成人av| 国产成人一区二区在线| 国产亚洲最大av| av在线app专区| 在线观看免费日韩欧美大片| 欧美 亚洲 国产 日韩一| 黄色配什么色好看| 国产视频首页在线观看| 巨乳人妻的诱惑在线观看| 最新中文字幕久久久久| 飞空精品影院首页| 国产精品偷伦视频观看了| 日韩欧美精品免费久久| 看免费成人av毛片| 观看美女的网站| 久久精品aⅴ一区二区三区四区 | 国产日韩一区二区三区精品不卡| 久久久精品免费免费高清| 国产成人免费观看mmmm| 黑人巨大精品欧美一区二区蜜桃 | 黑人巨大精品欧美一区二区蜜桃 | 亚洲av在线观看美女高潮| 成人毛片a级毛片在线播放| 国产成人免费无遮挡视频| 欧美变态另类bdsm刘玥| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 免费在线观看黄色视频的| 人妻系列 视频| 成人毛片60女人毛片免费| 国产免费一区二区三区四区乱码| 九草在线视频观看| 青青草视频在线视频观看| h视频一区二区三区| 我的女老师完整版在线观看| 十八禁高潮呻吟视频| 亚洲天堂av无毛| 一级黄片播放器| 永久网站在线| 亚洲精品国产av成人精品| 精品少妇黑人巨大在线播放| 国产亚洲最大av| 亚洲三级黄色毛片| av又黄又爽大尺度在线免费看| 国产深夜福利视频在线观看| 国产国语露脸激情在线看| 日韩一区二区三区影片| 欧美精品一区二区大全| 国产老妇伦熟女老妇高清| 亚洲成色77777| 日本爱情动作片www.在线观看| 欧美变态另类bdsm刘玥| www.熟女人妻精品国产 | 国产片内射在线| 欧美97在线视频| 精品久久久久久电影网| 久久精品国产鲁丝片午夜精品| 观看av在线不卡| 韩国av在线不卡| 边亲边吃奶的免费视频| a级毛片黄视频| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| 国产在线视频一区二区| 91在线精品国自产拍蜜月| 久久久国产欧美日韩av| 色网站视频免费| 国产亚洲av片在线观看秒播厂| 久久久久久久久久人人人人人人| 国产成人精品婷婷| 黄色 视频免费看| 国产av精品麻豆| 久久久久人妻精品一区果冻| 国产色爽女视频免费观看| 纵有疾风起免费观看全集完整版| 国产片内射在线| 伊人久久国产一区二区| 久久99蜜桃精品久久| 午夜老司机福利剧场| 涩涩av久久男人的天堂| 边亲边吃奶的免费视频| 亚洲情色 制服丝袜| 在线亚洲精品国产二区图片欧美| 免费观看在线日韩| 亚洲情色 制服丝袜| 亚洲精品国产av蜜桃| 国产精品.久久久| 日本猛色少妇xxxxx猛交久久| 在线免费观看不下载黄p国产| 国产成人一区二区在线| 视频区图区小说| 性色avwww在线观看| av视频免费观看在线观看| 国产一区二区在线观看日韩| 国产亚洲精品久久久com| 免费看不卡的av| 精品一区二区免费观看| 黄色视频在线播放观看不卡| 久久精品国产亚洲av涩爱| 亚洲成国产人片在线观看| 亚洲精品第二区| 久久久久久久久久人人人人人人| 亚洲精品国产色婷婷电影| 亚洲伊人色综图| 国产黄色免费在线视频| 国语对白做爰xxxⅹ性视频网站| 日韩在线高清观看一区二区三区| 日日撸夜夜添| 欧美3d第一页| 人成视频在线观看免费观看| 91精品国产国语对白视频| 国产爽快片一区二区三区| 国产一区二区三区av在线| 另类精品久久| 精品国产乱码久久久久久小说| 女人被躁到高潮嗷嗷叫费观| 人人妻人人澡人人看| 高清毛片免费看| 一区二区av电影网| 国产在线一区二区三区精| 婷婷色综合大香蕉| 男人爽女人下面视频在线观看| 日韩视频在线欧美| av线在线观看网站| 国产高清不卡午夜福利| 国产激情久久老熟女| 亚洲av免费高清在线观看| 成年av动漫网址| av在线观看视频网站免费| 欧美国产精品va在线观看不卡| 成年人免费黄色播放视频| 国产又色又爽无遮挡免| 亚洲国产精品999| 各种免费的搞黄视频| 三上悠亚av全集在线观看| 母亲3免费完整高清在线观看 | 国产欧美日韩综合在线一区二区| 亚洲人与动物交配视频| 一级毛片黄色毛片免费观看视频| av国产精品久久久久影院| 亚洲精品久久成人aⅴ小说| 99久久中文字幕三级久久日本| a 毛片基地| 日韩,欧美,国产一区二区三区| 亚洲国产av新网站| 一边摸一边做爽爽视频免费| 国产有黄有色有爽视频| 黄片播放在线免费| 精品国产一区二区三区四区第35| 国产欧美日韩综合在线一区二区| 侵犯人妻中文字幕一二三四区| 人体艺术视频欧美日本| 人人妻人人添人人爽欧美一区卜| 美女中出高潮动态图| 草草在线视频免费看| 天堂俺去俺来也www色官网| 久久人妻熟女aⅴ| 久久韩国三级中文字幕| 国产又色又爽无遮挡免| 精品午夜福利在线看| 国产欧美日韩一区二区三区在线| 免费高清在线观看视频在线观看| www日本在线高清视频| 国产精品.久久久| 最近中文字幕高清免费大全6| 国产在线免费精品| 国产在视频线精品| 美女大奶头黄色视频| 久久精品夜色国产| 亚洲天堂av无毛| 亚洲美女搞黄在线观看| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 久久午夜综合久久蜜桃| 精品酒店卫生间| 午夜91福利影院| 永久网站在线| 人妻系列 视频| 观看av在线不卡| 亚洲伊人色综图| 五月玫瑰六月丁香| 另类亚洲欧美激情| 久久97久久精品| 99视频精品全部免费 在线| 最后的刺客免费高清国语| 成年人免费黄色播放视频| av天堂久久9| 亚洲美女搞黄在线观看| 一级,二级,三级黄色视频| 欧美日韩国产mv在线观看视频| 男女午夜视频在线观看 | 大片免费播放器 马上看| 亚洲国产毛片av蜜桃av| 国产精品国产三级专区第一集| 校园人妻丝袜中文字幕| 欧美精品高潮呻吟av久久| 男女啪啪激烈高潮av片| 国产精品女同一区二区软件| 亚洲av在线观看美女高潮| 国产免费一级a男人的天堂| 免费大片18禁| 女人被躁到高潮嗷嗷叫费观| 国产永久视频网站| 夜夜骑夜夜射夜夜干| 欧美最新免费一区二区三区| 久久青草综合色| 亚洲成av片中文字幕在线观看 | 婷婷色综合大香蕉| 久久99热6这里只有精品| 日韩一本色道免费dvd| 国产老妇伦熟女老妇高清| 日本免费在线观看一区| 丝袜美足系列| 国产在线免费精品| 日韩成人av中文字幕在线观看| 日韩在线高清观看一区二区三区| 男女边吃奶边做爰视频| 久久狼人影院| 精品午夜福利在线看| 免费黄频网站在线观看国产| 日韩大片免费观看网站| 两个人免费观看高清视频| 久久热在线av| 成年人午夜在线观看视频| 亚洲一区二区三区欧美精品| 如何舔出高潮| 国产精品国产三级国产av玫瑰| 久久久精品免费免费高清| 日本av手机在线免费观看| 精品国产一区二区三区久久久樱花| 人妻少妇偷人精品九色| 美女主播在线视频| 免费黄频网站在线观看国产| 国产精品一区二区在线不卡| 国产成人精品福利久久| 久久久久网色| av免费观看日本| 少妇精品久久久久久久| 欧美精品国产亚洲| 久久精品国产亚洲av涩爱| 一区二区三区乱码不卡18| 欧美xxⅹ黑人| 国产69精品久久久久777片| 亚洲欧洲国产日韩| 日韩一区二区视频免费看| xxxhd国产人妻xxx| 永久网站在线| 国产精品久久久久久精品古装| 菩萨蛮人人尽说江南好唐韦庄| 青春草视频在线免费观看| 人人妻人人添人人爽欧美一区卜| 日本免费在线观看一区| 久久久精品区二区三区| 我的女老师完整版在线观看| 欧美老熟妇乱子伦牲交| videossex国产| 亚洲激情五月婷婷啪啪| 亚洲美女搞黄在线观看| 色哟哟·www| 岛国毛片在线播放| 丝袜脚勾引网站| 亚洲熟女精品中文字幕| 亚洲第一av免费看| 日本av手机在线免费观看| 免费女性裸体啪啪无遮挡网站| 麻豆精品久久久久久蜜桃| 亚洲国产精品专区欧美| 97超碰精品成人国产| 精品卡一卡二卡四卡免费| 日韩人妻精品一区2区三区| 国产精品久久久久久精品古装| 亚洲精品乱久久久久久| 日韩成人伦理影院| av电影中文网址| 欧美日韩成人在线一区二区| videos熟女内射| 在线观看国产h片| 最近中文字幕2019免费版| 免费黄色在线免费观看| 日本vs欧美在线观看视频| 国产精品人妻久久久久久| av国产久精品久网站免费入址| 母亲3免费完整高清在线观看 | 精品国产一区二区三区四区第35| 久久精品久久精品一区二区三区| 最近最新中文字幕大全免费视频 | 18禁裸乳无遮挡动漫免费视频| 尾随美女入室| 久久国内精品自在自线图片| 国产精品一国产av| 巨乳人妻的诱惑在线观看| 看免费av毛片| av网站免费在线观看视频| 熟妇人妻不卡中文字幕| 精品一区在线观看国产| 99精国产麻豆久久婷婷| 国产69精品久久久久777片| 免费不卡的大黄色大毛片视频在线观看| 亚洲精华国产精华液的使用体验| 人妻人人澡人人爽人人| 久久99精品国语久久久| 一个人免费看片子| 日本爱情动作片www.在线观看| 亚洲一码二码三码区别大吗| 久久久久久久国产电影| 男女高潮啪啪啪动态图| 99国产综合亚洲精品| 热re99久久精品国产66热6| av又黄又爽大尺度在线免费看| 日本黄色日本黄色录像| 国产片特级美女逼逼视频| 久久精品国产自在天天线| 青春草视频在线免费观看| 国产免费视频播放在线视频| 亚洲一区二区三区欧美精品| 涩涩av久久男人的天堂| 综合色丁香网| 国产黄色视频一区二区在线观看| 亚洲av成人精品一二三区| 国产又色又爽无遮挡免| 国产在视频线精品| 国产熟女午夜一区二区三区| 久久精品aⅴ一区二区三区四区 | 最近的中文字幕免费完整| 亚洲国产av新网站| 午夜福利,免费看| 在线 av 中文字幕| 精品国产乱码久久久久久小说| 午夜视频国产福利| 一本久久精品| 婷婷色麻豆天堂久久| 久久久a久久爽久久v久久| 看免费av毛片| 男女啪啪激烈高潮av片| av黄色大香蕉| 久久99热这里只频精品6学生| 母亲3免费完整高清在线观看 | 国产成人午夜福利电影在线观看| 婷婷色综合www| 中文字幕人妻熟女乱码| 久久精品aⅴ一区二区三区四区 | 少妇熟女欧美另类| 日韩中文字幕视频在线看片| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区国产| 女人精品久久久久毛片| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区国产| av免费在线看不卡| 精品少妇黑人巨大在线播放| 最近2019中文字幕mv第一页| av片东京热男人的天堂| 欧美人与善性xxx| 国产成人aa在线观看| 成人国产av品久久久| 国产亚洲av片在线观看秒播厂| 久久久久久久精品精品| 亚洲激情五月婷婷啪啪| 又黄又爽又刺激的免费视频.| 91在线精品国自产拍蜜月| 亚洲精品日韩在线中文字幕| 国产精品久久久av美女十八| av在线播放精品| 2022亚洲国产成人精品| 成人综合一区亚洲| 亚洲国产最新在线播放|