• <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新网站| 免费观看性生交大片5| 国产黄色视频一区二区在线观看| 美女视频免费永久观看网站| 不卡视频在线观看欧美| 欧美日韩精品成人综合77777| 男女国产视频网站| 亚洲成人久久爱视频| 国产高清三级在线| 人体艺术视频欧美日本| 亚洲国产av新网站| 春色校园在线视频观看| 久久久久久久大尺度免费视频| 欧美高清成人免费视频www| 国产av码专区亚洲av| 视频区图区小说| 在线观看美女被高潮喷水网站| 狂野欧美激情性bbbbbb| 18+在线观看网站| 亚洲欧美精品自产自拍| 国产av不卡久久| 亚洲欧美清纯卡通| 精品一区二区三卡| 国产午夜福利久久久久久| 精品久久久久久电影网| 一本久久精品| 真实男女啪啪啪动态图| 国产精品一区二区性色av| 高清日韩中文字幕在线| 日韩欧美精品免费久久| 在线观看一区二区三区激情| 1000部很黄的大片| 777米奇影视久久| 国精品久久久久久国模美| 国产69精品久久久久777片| 赤兔流量卡办理| 两个人的视频大全免费| a级毛片免费高清观看在线播放| 亚洲精品日本国产第一区| 国产 一区 欧美 日韩| 国产男女超爽视频在线观看| 亚洲精品乱久久久久久| 精品一区二区三区视频在线| 又黄又爽又刺激的免费视频.| 国产老妇伦熟女老妇高清| 国产欧美另类精品又又久久亚洲欧美| 看十八女毛片水多多多| 亚洲国产精品999| 别揉我奶头 嗯啊视频| 成人二区视频| 国产爽快片一区二区三区| av黄色大香蕉| 午夜免费观看性视频| 免费看光身美女| 亚洲欧美精品自产自拍| 亚洲欧美精品专区久久| 免费av毛片视频| 国产视频首页在线观看| 嫩草影院入口| 99热这里只有是精品在线观看| 中文在线观看免费www的网站| 国产 一区 欧美 日韩| 人妻制服诱惑在线中文字幕| 高清日韩中文字幕在线| 国产淫片久久久久久久久| 国产色婷婷99| 91午夜精品亚洲一区二区三区| 亚洲欧洲日产国产| 美女高潮的动态| 午夜视频国产福利| 午夜激情福利司机影院| 夜夜爽夜夜爽视频| 美女视频免费永久观看网站| 国产成人精品一,二区| 又爽又黄a免费视频| 欧美成人精品欧美一级黄| 成年免费大片在线观看| 毛片女人毛片| 国产av国产精品国产| 成人黄色视频免费在线看| 少妇 在线观看| 亚洲av免费高清在线观看| 18禁动态无遮挡网站| 少妇人妻精品综合一区二区| 在线免费十八禁| 人妻系列 视频| 亚洲国产精品国产精品| 久久精品熟女亚洲av麻豆精品| 极品少妇高潮喷水抽搐| 内地一区二区视频在线| 亚洲丝袜综合中文字幕| 白带黄色成豆腐渣| 国产免费一级a男人的天堂| 日本午夜av视频| 午夜激情久久久久久久| 一区二区三区四区激情视频| 亚洲婷婷狠狠爱综合网| 韩国高清视频一区二区三区| 久久久久久久亚洲中文字幕| 不卡视频在线观看欧美| 男女啪啪激烈高潮av片| 国产亚洲午夜精品一区二区久久 | 久久久精品免费免费高清| 日韩av不卡免费在线播放| 精品一区在线观看国产| 久久99热6这里只有精品| 人妻一区二区av| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 黄色欧美视频在线观看| 午夜日本视频在线| 在线观看免费高清a一片| 日本三级黄在线观看| 亚洲电影在线观看av| 自拍欧美九色日韩亚洲蝌蚪91 | 日本一本二区三区精品| 欧美精品人与动牲交sv欧美| 视频中文字幕在线观看| 又黄又爽又刺激的免费视频.| 国内精品美女久久久久久| 欧美另类一区| 精品99又大又爽又粗少妇毛片| 大陆偷拍与自拍| 精品国产乱码久久久久久小说| 欧美成人a在线观看| 在线观看人妻少妇| 国产乱人视频| 欧美激情在线99| 亚洲av在线观看美女高潮| 插逼视频在线观看| 搡老乐熟女国产| 国产精品久久久久久精品古装| 99久久九九国产精品国产免费| 少妇人妻久久综合中文| 久久久久精品性色| 免费看光身美女| 久久久久久伊人网av| 男人和女人高潮做爰伦理| 人妻 亚洲 视频| 国产黄频视频在线观看| xxx大片免费视频| 免费播放大片免费观看视频在线观看| 99热6这里只有精品| 18禁在线播放成人免费| 一级二级三级毛片免费看| 女人十人毛片免费观看3o分钟| 国产精品嫩草影院av在线观看| 如何舔出高潮| 午夜免费男女啪啪视频观看| 亚洲第一区二区三区不卡| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩视频高清一区二区三区二| 国产精品一区www在线观看| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 深夜a级毛片| 我要看日韩黄色一级片| 国产片特级美女逼逼视频| 亚洲国产精品国产精品| 日韩电影二区| 午夜视频国产福利| 久久人人爽人人爽人人片va| 一级a做视频免费观看| 在线观看美女被高潮喷水网站| 久久午夜福利片| 97精品久久久久久久久久精品| 欧美潮喷喷水| 亚洲av二区三区四区| 久久久亚洲精品成人影院| 免费不卡的大黄色大毛片视频在线观看| 午夜福利视频精品| 韩国av在线不卡| 免费看不卡的av| 22中文网久久字幕| xxx大片免费视频| 午夜福利视频1000在线观看| 国产爽快片一区二区三区| 人妻制服诱惑在线中文字幕| 国产精品福利在线免费观看| 久久久精品94久久精品| 日本欧美国产在线视频| 高清在线视频一区二区三区| 九草在线视频观看| 日韩伦理黄色片| 91在线精品国自产拍蜜月| 国产亚洲av片在线观看秒播厂| 亚洲av国产av综合av卡| 插阴视频在线观看视频| 在线精品无人区一区二区三 | 免费大片18禁| tube8黄色片| 亚洲精品久久久久久婷婷小说| 欧美激情在线99| 国产免费又黄又爽又色| 久久精品国产自在天天线| 成年人午夜在线观看视频| 天天一区二区日本电影三级| 精品久久国产蜜桃| 亚洲美女搞黄在线观看| 91aial.com中文字幕在线观看| 国模一区二区三区四区视频| 中国美白少妇内射xxxbb| 久久久久久伊人网av| 日韩一区二区三区影片| 两个人的视频大全免费| 欧美日韩国产mv在线观看视频 | 午夜日本视频在线| 久久久久九九精品影院| 国产精品.久久久| 欧美国产精品一级二级三级 | 99热这里只有是精品在线观看| 色网站视频免费| 精品一区在线观看国产| 成人美女网站在线观看视频| 国产午夜精品久久久久久一区二区三区| 欧美日韩国产mv在线观看视频 | 深夜a级毛片| 日本欧美国产在线视频| 老司机影院毛片| 青春草亚洲视频在线观看| 国产乱人偷精品视频| 免费高清在线观看视频在线观看| 久久久久九九精品影院| 久久久久久久久久久免费av| 丰满乱子伦码专区| 亚洲av在线观看美女高潮| 久久久午夜欧美精品| 欧美成人一区二区免费高清观看| 欧美人与善性xxx| 黑人高潮一二区| 欧美 日韩 精品 国产| 国内精品美女久久久久久| 一区二区三区四区激情视频| 尾随美女入室| a级毛片免费高清观看在线播放| 尤物成人国产欧美一区二区三区| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产| 少妇高潮的动态图| 欧美一区二区亚洲| 欧美高清成人免费视频www| 亚洲一级一片aⅴ在线观看| 久久精品久久久久久久性| av.在线天堂| 一区二区三区四区激情视频| 人妻一区二区av| 国产成人精品婷婷| 晚上一个人看的免费电影| 永久免费av网站大全| 91午夜精品亚洲一区二区三区| av卡一久久| 91狼人影院| 五月玫瑰六月丁香| 免费看日本二区| 精品人妻偷拍中文字幕| 中文乱码字字幕精品一区二区三区| 一个人看视频在线观看www免费| 高清av免费在线| 直男gayav资源| 亚洲精品视频女| 亚洲,一卡二卡三卡| 亚洲天堂国产精品一区在线| 亚洲美女视频黄频| 成人二区视频| 街头女战士在线观看网站| 在线免费十八禁| 最近中文字幕2019免费版| 99九九线精品视频在线观看视频| 80岁老熟妇乱子伦牲交| 日本wwww免费看| 中文乱码字字幕精品一区二区三区| 一个人看的www免费观看视频| 深爱激情五月婷婷| 亚洲人成网站在线观看播放| 高清av免费在线| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美精品v在线| av女优亚洲男人天堂| 成人亚洲欧美一区二区av| 精品少妇久久久久久888优播| 赤兔流量卡办理| 国产精品人妻久久久影院| 丝瓜视频免费看黄片| av在线观看视频网站免费| 热99国产精品久久久久久7| 一个人看的www免费观看视频| 在线观看美女被高潮喷水网站| 久热这里只有精品99| 日日撸夜夜添| av在线app专区| av免费观看日本| 国产 精品1| 日韩,欧美,国产一区二区三区| 久久精品夜色国产| 老女人水多毛片| 日韩一区二区视频免费看| 不卡视频在线观看欧美| 中文精品一卡2卡3卡4更新| 成人二区视频| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 毛片一级片免费看久久久久| 日韩电影二区| 国产男女内射视频| 成人一区二区视频在线观看| 日韩制服骚丝袜av| 亚洲精品国产av成人精品| 九色成人免费人妻av| 97在线视频观看| 亚洲精品第二区| 人人妻人人爽人人添夜夜欢视频 | 六月丁香七月| 亚洲av成人精品一区久久| 亚洲精品乱码久久久久久按摩| 久久99热这里只频精品6学生| 国产在线一区二区三区精| 精品人妻偷拍中文字幕| 久久精品久久精品一区二区三区| 尤物成人国产欧美一区二区三区| 天天躁夜夜躁狠狠久久av| 最近中文字幕高清免费大全6| 老师上课跳d突然被开到最大视频| 精品人妻偷拍中文字幕| 久久韩国三级中文字幕| 成人一区二区视频在线观看| 免费观看a级毛片全部| 麻豆精品久久久久久蜜桃| 成人毛片a级毛片在线播放| 久久久久九九精品影院| 天堂中文最新版在线下载 | 伦精品一区二区三区| 欧美少妇被猛烈插入视频| 欧美潮喷喷水| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 大码成人一级视频| 成人美女网站在线观看视频| 久久人人爽人人爽人人片va| 亚洲丝袜综合中文字幕| 一级毛片电影观看| 在线播放无遮挡| 男女边吃奶边做爰视频| 日韩欧美一区视频在线观看 | 亚洲伊人久久精品综合| 亚洲国产精品国产精品| 丰满人妻一区二区三区视频av| 亚洲美女搞黄在线观看| 最新中文字幕久久久久| 久久精品国产a三级三级三级| 国内揄拍国产精品人妻在线| 高清毛片免费看| 欧美日韩精品成人综合77777| 国产精品伦人一区二区| 少妇丰满av| 国产欧美日韩一区二区三区在线 | 日韩视频在线欧美| 亚洲国产欧美在线一区| 成人欧美大片| 久久久精品欧美日韩精品| 免费av观看视频| 欧美少妇被猛烈插入视频| 午夜视频国产福利| 99九九线精品视频在线观看视频| 久久精品国产亚洲av涩爱| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放| 日日摸夜夜添夜夜爱| 精品酒店卫生间| 91久久精品国产一区二区成人| 亚洲国产日韩一区二区| 亚洲高清免费不卡视频| 国产老妇女一区| 一区二区三区乱码不卡18| 国产高清不卡午夜福利| 亚洲最大成人手机在线| 日韩人妻高清精品专区| 午夜激情福利司机影院| 91aial.com中文字幕在线观看| 免费看av在线观看网站| 又大又黄又爽视频免费| 久久久久久久久久久丰满| 日韩不卡一区二区三区视频在线| 亚洲av在线观看美女高潮| 国产成人91sexporn| 伊人久久国产一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本与韩国留学比较| 天堂俺去俺来也www色官网| 国产精品国产三级专区第一集| 国产精品一二三区在线看| 91aial.com中文字幕在线观看| 嘟嘟电影网在线观看| 日本色播在线视频| 亚洲人成网站在线播| 又大又黄又爽视频免费| 观看免费一级毛片| 久久人人爽人人爽人人片va| 免费av观看视频| 青青草视频在线视频观看| 高清视频免费观看一区二区| 国产成人免费无遮挡视频| 国产精品一区二区在线观看99| 好男人视频免费观看在线| 高清午夜精品一区二区三区| 亚洲欧美中文字幕日韩二区| 亚洲精品乱码久久久久久按摩| 国产乱来视频区| 少妇人妻一区二区三区视频| 久久精品国产亚洲av天美| 日韩制服骚丝袜av| 夜夜爽夜夜爽视频| 舔av片在线| 十八禁网站网址无遮挡 | 在线观看国产h片| 亚洲在久久综合| 一个人观看的视频www高清免费观看| 国产精品.久久久| 丝瓜视频免费看黄片| 草草在线视频免费看| 精品亚洲乱码少妇综合久久| 国产成人精品婷婷| 欧美日本视频| 国产av国产精品国产| 老司机影院毛片| 大香蕉97超碰在线| 在线看a的网站| 精品一区二区三区视频在线| 免费高清在线观看视频在线观看| 成人亚洲精品av一区二区| 国产成人91sexporn| 成人美女网站在线观看视频| 日韩强制内射视频| 色综合色国产| 亚洲精品国产av成人精品| 纵有疾风起免费观看全集完整版| 秋霞在线观看毛片| 少妇丰满av| 国产亚洲av片在线观看秒播厂| 国产一区亚洲一区在线观看| 亚洲精华国产精华液的使用体验| 亚洲国产欧美人成| 91精品一卡2卡3卡4卡| 久久精品久久精品一区二区三区| 欧美区成人在线视频| 国产精品熟女久久久久浪| 亚洲精品亚洲一区二区| 亚洲av不卡在线观看| 七月丁香在线播放| 欧美变态另类bdsm刘玥| 国产av国产精品国产| 99视频精品全部免费 在线| 国产精品久久久久久精品电影小说 | 成人欧美大片| 午夜免费鲁丝| 熟女电影av网| 亚洲精品成人av观看孕妇| 日韩欧美精品v在线| av线在线观看网站| 国产亚洲91精品色在线| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久久大av| 男女无遮挡免费网站观看| 狂野欧美激情性xxxx在线观看| 在线精品无人区一区二区三 | 亚洲天堂国产精品一区在线| 欧美国产精品一级二级三级 | 王馨瑶露胸无遮挡在线观看| 欧美xxxx黑人xx丫x性爽| 尾随美女入室| 各种免费的搞黄视频| 狠狠精品人妻久久久久久综合| 天天躁夜夜躁狠狠久久av| 黄色视频在线播放观看不卡| 欧美区成人在线视频| 黄色怎么调成土黄色| 99久久精品一区二区三区| 精品久久久久久久久av| 亚洲av在线观看美女高潮| 3wmmmm亚洲av在线观看| 一区二区av电影网| 春色校园在线视频观看| 特级一级黄色大片| 狂野欧美激情性xxxx在线观看| 在线观看一区二区三区激情| 国产伦精品一区二区三区视频9| 国产一级毛片在线| 国产真实伦视频高清在线观看| 国产69精品久久久久777片| 免费高清在线观看视频在线观看| 少妇人妻精品综合一区二区| 91午夜精品亚洲一区二区三区| 国产 精品1| 国产爱豆传媒在线观看| 国产精品久久久久久精品电影小说 | 婷婷色av中文字幕| 黄片无遮挡物在线观看| 国产亚洲一区二区精品| 国产精品成人在线| 啦啦啦啦在线视频资源| 亚洲精品日韩av片在线观看| 欧美成人精品欧美一级黄| 黄片wwwwww| 亚洲精品456在线播放app| 国产一区二区三区av在线| av在线app专区| 日韩强制内射视频| 国产高清不卡午夜福利| 免费黄色在线免费观看| 久久99热这里只有精品18| 久久久精品免费免费高清| 国精品久久久久久国模美| 三级国产精品欧美在线观看| 极品教师在线视频| 大片电影免费在线观看免费| 十八禁网站网址无遮挡 | 亚洲高清免费不卡视频| 天天躁夜夜躁狠狠久久av| 久久久久久久国产电影| 一级a做视频免费观看| 九草在线视频观看| 午夜爱爱视频在线播放| 看非洲黑人一级黄片| 久久精品熟女亚洲av麻豆精品| 色视频在线一区二区三区| av.在线天堂| 啦啦啦啦在线视频资源| 日韩三级伦理在线观看| 日韩国内少妇激情av| 久久久色成人| 久久国产乱子免费精品| 国产大屁股一区二区在线视频| 成人漫画全彩无遮挡| 国产熟女欧美一区二区| 超碰97精品在线观看| 亚洲自偷自拍三级| 日韩,欧美,国产一区二区三区| 日本欧美国产在线视频| 欧美日韩精品成人综合77777| 联通29元200g的流量卡| 久久韩国三级中文字幕| 亚洲精品自拍成人| 日本三级黄在线观看| 国产精品爽爽va在线观看网站| 成人亚洲欧美一区二区av| a级毛色黄片| 亚洲精品乱久久久久久| 伊人久久精品亚洲午夜| 下体分泌物呈黄色| 一区二区av电影网| 在线播放无遮挡| 美女cb高潮喷水在线观看| 久久99精品国语久久久| 日日啪夜夜爽| 听说在线观看完整版免费高清| 夜夜看夜夜爽夜夜摸| 免费看日本二区| 一边亲一边摸免费视频| 亚洲经典国产精华液单| 久久精品夜色国产| 国产成人91sexporn| 色5月婷婷丁香| 日本一二三区视频观看| 亚洲精品成人久久久久久| 亚洲三级黄色毛片| 亚洲aⅴ乱码一区二区在线播放| 中文天堂在线官网| 久久久久久久久久久丰满| 欧美性感艳星| 我的老师免费观看完整版| 免费黄网站久久成人精品| 狂野欧美激情性xxxx在线观看| 国产精品久久久久久久电影| 欧美3d第一页| 中文精品一卡2卡3卡4更新| 女人久久www免费人成看片| 亚洲成人久久爱视频| 国产成人91sexporn| 赤兔流量卡办理| av国产精品久久久久影院| .国产精品久久| 亚洲国产色片| 九九爱精品视频在线观看| 精华霜和精华液先用哪个| 新久久久久国产一级毛片| 啦啦啦在线观看免费高清www| 夜夜爽夜夜爽视频| 亚洲精品国产色婷婷电影| 欧美+日韩+精品| 免费大片18禁| 久久久久久久久久人人人人人人| 亚洲综合精品二区| 亚洲av欧美aⅴ国产| 国产美女午夜福利| 深夜a级毛片| 亚洲精品一区蜜桃| 成人国产av品久久久| 亚洲综合色惰| 午夜免费观看性视频| 欧美性感艳星| av天堂中文字幕网| 人妻少妇偷人精品九色| 精品人妻视频免费看| 97在线人人人人妻| 高清av免费在线| 国产探花极品一区二区| 国语对白做爰xxxⅹ性视频网站| 天天躁夜夜躁狠狠久久av| 女人被狂操c到高潮| 国产精品无大码| 香蕉精品网在线| 美女国产视频在线观看| 小蜜桃在线观看免费完整版高清| 午夜老司机福利剧场|