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

    有關(guān)近空間高超聲速飛行器邊界層轉(zhuǎn)捩和湍流的兩個(gè)問(wèn)題

    2017-04-28 04:52:37張涵信
    關(guān)鍵詞:馬赫數(shù)粘性邊界層

    周 恒, 張涵信

    (1.天津大學(xué) 力學(xué)系, 天津 300072; 2. 中國(guó)空氣動(dòng)力研究與發(fā)展中心 國(guó)家計(jì)算流體力學(xué)實(shí)驗(yàn)室, 北京 100191)

    ?

    有關(guān)近空間高超聲速飛行器邊界層轉(zhuǎn)捩和湍流的兩個(gè)問(wèn)題

    周 恒1,*, 張涵信2

    (1.天津大學(xué) 力學(xué)系, 天津 300072; 2. 中國(guó)空氣動(dòng)力研究與發(fā)展中心 國(guó)家計(jì)算流體力學(xué)實(shí)驗(yàn)室, 北京 100191)

    和一般的飛行器一樣,在近空間飛行器的研制中,其邊界層的轉(zhuǎn)捩和湍流也是需要考慮的兩個(gè)重要問(wèn)題。但即使是對(duì)一般的飛行器,“轉(zhuǎn)捩”和“湍流”也還是兩個(gè)歷經(jīng)百年而仍未很好解決的問(wèn)題,而對(duì)近空間飛行器來(lái)說(shuō),空氣動(dòng)力學(xué)本身就還存在若干新的需要研究解決的基礎(chǔ)問(wèn)題,邊界層的轉(zhuǎn)捩和湍流就更是沒有很好解決的問(wèn)題。本文討論了兩個(gè)問(wèn)題:1) 為增強(qiáng)對(duì)高超聲速飛行器邊界層轉(zhuǎn)捩預(yù)測(cè)的能力,需要開展哪些方面的研究工作及其困難;2) 是否有可能當(dāng)飛行器飛行高度足夠大時(shí),其邊界層就不會(huì)再有湍流問(wèn)題?

    近空間飛行器;轉(zhuǎn)捩;湍流;預(yù)測(cè)能力

    1 轉(zhuǎn)捩預(yù)測(cè)問(wèn)題

    轉(zhuǎn)捩問(wèn)題歷經(jīng)百年的研究,在低速流方面已經(jīng)有了不小的進(jìn)展。而對(duì)高超聲速邊界層的轉(zhuǎn)捩問(wèn)題,近年來(lái)的研究雖有一些進(jìn)展,但仍存在嚴(yán)重的不足[1-2]。

    可以肯定,整個(gè)轉(zhuǎn)捩過(guò)程可分為三個(gè)階段或作為三個(gè)問(wèn)題來(lái)研究。1) 感受性問(wèn)題,即邊界層外的擾動(dòng)如何引發(fā)邊界層內(nèi)以不穩(wěn)定波形式出現(xiàn)的擾動(dòng)?2) 被引發(fā)的不穩(wěn)定波在邊界層內(nèi)的演化。如果開始時(shí)幅值很小,例如一般認(rèn)為,如果不穩(wěn)定波的速度的幅值小于邊界層外緣速度的1%,則可以用線性理論;而大于1%,則要用非線性理論。3) 擾動(dòng)演化到什么情況將觸發(fā)轉(zhuǎn)捩?即預(yù)測(cè)轉(zhuǎn)捩發(fā)生的判據(jù)。

    以上的第二個(gè)問(wèn)題研究的時(shí)間最長(zhǎng)也最充分。特別是現(xiàn)在有了計(jì)算機(jī)和各種計(jì)算方法,只要知道擾動(dòng)的初始條件,已經(jīng)沒有原則上的困難了。

    第三個(gè)問(wèn)題原來(lái)是一個(gè)難題。長(zhǎng)期以來(lái)人們提出多種邊界層內(nèi)擾動(dòng)演化的非線性理論,但始終無(wú)法和第三個(gè)問(wèn)題,即演化到什么程度將觸發(fā)轉(zhuǎn)捩,聯(lián)系起來(lái)。在文獻(xiàn)[2]中,作者通過(guò)多種轉(zhuǎn)捩的直接數(shù)值模擬,發(fā)現(xiàn)了轉(zhuǎn)捩過(guò)程的機(jī)理是:擾動(dòng)通過(guò)非線性作用對(duì)平均流進(jìn)行修正,使得平均流剖面的線性穩(wěn)定性發(fā)生變化。如果修正的結(jié)果是線性不穩(wěn)定的參數(shù)區(qū)(例如線性穩(wěn)定性理論的中性曲線所包含的區(qū)域)擴(kuò)大,則會(huì)形成不穩(wěn)定性的正反饋,從而觸發(fā)轉(zhuǎn)捩。他們并總結(jié)出,在技術(shù)人員喜歡用的半經(jīng)驗(yàn)轉(zhuǎn)捩預(yù)測(cè)方法eN法中,當(dāng)不穩(wěn)定波幅值達(dá)到1.5%時(shí)轉(zhuǎn)捩將發(fā)生。這可以認(rèn)為第三個(gè)問(wèn)題在一定程度上已經(jīng)解決。

    現(xiàn)今已有不少作者從事這方面的計(jì)算研究工作,例如Zhong等采用激波裝配法[4-9],Balakurma等則采用了激波捕捉法[10-11],然而他們對(duì)激波捕捉法的正確使用和激波后的擾動(dòng)的非定常分析均存在缺陷。有文給出了線性擾動(dòng)解,還用平行流假定作聲波、熵波和渦波分析,對(duì)高超聲速繞流這都是欠適當(dāng)?shù)?。?duì)有復(fù)雜激波系的問(wèn)題,用激波裝配法來(lái)研究擾動(dòng)通過(guò)激波經(jīng)歷的變化是不現(xiàn)實(shí)的,而用激波捕捉法是否能可靠地反映這一變化則有待于新的能證實(shí)其可靠性的研究結(jié)果。

    因此,為了解決轉(zhuǎn)捩預(yù)測(cè)問(wèn)題,當(dāng)前最需要做的基礎(chǔ)研究就是超聲速/高超聲速條件下的感受性問(wèn)題。

    前文已述,解決感受性問(wèn)題的一個(gè)攔路虎就是飛行環(huán)境下的擾動(dòng)是什么不清楚,為此就要做實(shí)地測(cè)量的飛行試驗(yàn)。

    從Zhong等已有的直接數(shù)值模擬結(jié)果來(lái)看,激波后的聲波在感受性問(wèn)題中起主要作用,而激波前的聲波又在產(chǎn)生激波后的聲波中起主要作用。但在流體中,聲波的產(chǎn)生只有三種來(lái)源:一是由流體中的外來(lái)物體的體積周期性變化所引起;二是由外來(lái)物體對(duì)流體的周期性作用力所引起;三是由流體本身的運(yùn)動(dòng)產(chǎn)生了周期性的雷諾應(yīng)力所引起。而在高空,實(shí)際上不存在外來(lái)物體(不算飛行器本身),因此可以認(rèn)為不存在前兩種聲源。但第三種聲源激發(fā)聲波的效率很低,可以忽略不計(jì)。所以,飛行試驗(yàn)的目標(biāo)就是測(cè)出大氣中的以速度和溫度擾動(dòng)出現(xiàn)的渦波和熵波。

    但這是非常困難的任務(wù)。首先,如果用超聲速/高超聲速飛行器作為儀器載體,則儀器只能感受飛行器頭激波后的擾動(dòng)。要還原成大氣中原來(lái)的擾動(dòng),是要解一個(gè)反問(wèn)題,而反問(wèn)題的解是不唯一的。其次,儀器的傳感器要有非常高的靈敏度。例如,一般估計(jì)對(duì)馬赫數(shù)為20的飛行器,速度背景擾動(dòng)的幅值為萬(wàn)分之一的量級(jí)。而如果馬赫數(shù)為6時(shí),同樣的擾動(dòng)對(duì)應(yīng)的無(wú)量綱幅值就是萬(wàn)分之三的量級(jí)。因此,測(cè)量時(shí)儀器就要有能分辨速度的萬(wàn)分之一或萬(wàn)分之三的變化的能力。而頻率范圍則要高達(dá)幾十萬(wàn)赫茲,對(duì)溫度傳感器的要求也類似。如果要用亞聲速的飛行器作為儀器載體,則對(duì)儀器傳感器的靈敏度和頻響范圍的要求低得多,但如何將該飛行器送上高空后能有一段亞聲速的平飛段是一難題,而且亞聲速的飛行器本身對(duì)大氣流動(dòng)的干擾會(huì)反映到傳感器上,如何排除?有人認(rèn)為在30 km以下似乎可以用探空氣球作為儀器載體。但探空氣球的存在改變了周圍相當(dāng)大的一片區(qū)域的流場(chǎng),儀器的傳感器能否伸出到離氣球足夠遠(yuǎn)的地方?

    通常,在飛行器的設(shè)計(jì)中,地面的風(fēng)洞實(shí)驗(yàn)起了很大的作用,但要將地面風(fēng)洞作為工程問(wèn)題的轉(zhuǎn)捩預(yù)測(cè)的手段,是不夠的。如前所述,人們?cè)缇驼J(rèn)識(shí)到對(duì)超聲速流而言,小擾動(dòng)可分解為三種擾動(dòng)。風(fēng)洞的實(shí)驗(yàn)段收縮比加大后,三種擾動(dòng)更小,于是就認(rèn)為風(fēng)洞氣流幾乎沒有聲、熵、渦擾動(dòng)。這在最初只是從減小背景擾動(dòng)幅值的結(jié)果來(lái)考慮的。但現(xiàn)在知道,從感受性的角度看,聲波擾動(dòng)、渦波擾動(dòng)及熵波擾動(dòng)的作用是不同的。同樣幅值的聲波和渦波導(dǎo)致的轉(zhuǎn)捩結(jié)果是不同的[3]。而上面的討論中已提到,高空大氣中本質(zhì)上不存在聲波。而靜風(fēng)洞的出現(xiàn)則在不自覺的條件下正好滿足了這一要求(靜風(fēng)洞的出現(xiàn)比文獻(xiàn)[3]中所紀(jì)錄的聲波擾動(dòng)和渦波擾動(dòng)導(dǎo)致的小鈍頭錐上的轉(zhuǎn)捩線定性上就不同的研究結(jié)果要早得多)。但如果要模擬真實(shí)的轉(zhuǎn)捩,則靜風(fēng)洞中的剩余擾動(dòng)必須和真實(shí)情況下的背景擾動(dòng)一致(無(wú)量綱化后一致),而我們實(shí)際上并不知道實(shí)際擾動(dòng)是什么(如上所述,需要由飛行試驗(yàn)測(cè)定)。而且即使知道了,要在風(fēng)洞中實(shí)現(xiàn)(至少是對(duì)轉(zhuǎn)捩起作用的擾動(dòng)頻率段的譜要一致)也是很大的難題。

    此外,還有其他問(wèn)題使得現(xiàn)有的風(fēng)洞不能完全滿足預(yù)測(cè)工程問(wèn)題中的轉(zhuǎn)捩的要求。這是因?yàn)椋?)N-S方程無(wú)量綱化后,和不可壓縮流的情況不同,方程的系數(shù)并不是常數(shù),而還有無(wú)量綱化后的粘性系數(shù)和普朗特?cái)?shù)。粘性系數(shù)和溫度有關(guān),而且對(duì)高超聲速流來(lái)說(shuō),實(shí)際飛行中邊界層的溫度還可能很高,從而粘性系數(shù)還要受氣體分子內(nèi)能被激發(fā)、離解、甚至化學(xué)反應(yīng)的影響。因此,要做到模型和實(shí)際飛行器的無(wú)量綱方程中的粘性系數(shù)一致,只有二者的溫度場(chǎng)完全一致才行,而這在現(xiàn)有的風(fēng)洞中是不可能做到的。2)模型和實(shí)際飛行器的表面溫度還要一致,因?yàn)檫吔鐚觾?nèi)的不穩(wěn)定波的演化受溫度邊界條件的影響是很大的,而這一項(xiàng)的要求顯然也不可能在實(shí)驗(yàn)中實(shí)現(xiàn)。因此,做到風(fēng)洞試驗(yàn)和工程實(shí)際的完全相似目前是做不到的。

    既然這樣,那風(fēng)洞實(shí)驗(yàn)還能起什么作用呢?我們認(rèn)為,目前它能起的最大作用是能提供有較完備參數(shù)的實(shí)驗(yàn)結(jié)果,可以作為檢驗(yàn)轉(zhuǎn)捩預(yù)測(cè)方法的依據(jù)。即將轉(zhuǎn)捩預(yù)測(cè)方法先直接面對(duì)風(fēng)洞和模型的參數(shù),而不是針對(duì)實(shí)際飛行器的參數(shù),以檢驗(yàn)其預(yù)測(cè)轉(zhuǎn)捩的結(jié)果是否和風(fēng)洞實(shí)驗(yàn)結(jié)果相符。但這樣做,就要求給出風(fēng)洞實(shí)驗(yàn)段的背景擾動(dòng)的頻率譜及模型壁面溫度(不是實(shí)驗(yàn)開始前的溫度,而是實(shí)驗(yàn)過(guò)程中的溫度)分布。

    2 湍流問(wèn)題

    迄今為止,在近空間以高超聲速飛行的飛行器是否必須考慮湍流問(wèn)題還沒有一個(gè)可靠的說(shuō)法。實(shí)際飛行實(shí)驗(yàn)的結(jié)果非常少。經(jīng)常被引用的只有幾篇文章,都還是20世紀(jì)七、八十年代的。文獻(xiàn)[12]是關(guān)于一個(gè)長(zhǎng)近4 m、半錐角為5°的小鈍頭錐,從約100 km高空返回地面的過(guò)程中做的試驗(yàn)。從100 km高空到約13 km高空,其飛行馬赫數(shù)都在20左右。結(jié)果發(fā)現(xiàn)一直降到30 km左右,才開始在尾部出現(xiàn)湍流。隨著高度不斷降低,轉(zhuǎn)捩線往前移,即湍流區(qū)不斷擴(kuò)大。文獻(xiàn)[13]是1974年發(fā)表的,研究對(duì)象為長(zhǎng)67.23 cm、鼻錐半徑16.09 cm、底部直徑48.69 cm的9°鼻錐。由拖曳火箭把它送上天,以5 km/s、迎角12.5°再入飛行,錐面有泰氟隆涂層,實(shí)驗(yàn)結(jié)果用粗糙度稍大的轉(zhuǎn)捩準(zhǔn)則處理。當(dāng)飛行高度是36.58 km以上,氣流全為層流;30.48 km時(shí),湍流區(qū)已擴(kuò)展至鼻端,即高度為36.58 km以下會(huì)有湍流的影響。文獻(xiàn)[14]也有類似的結(jié)果。后來(lái)美國(guó)根據(jù)自由飛和風(fēng)洞實(shí)驗(yàn)結(jié)果,提出了不同粗糙度轉(zhuǎn)捩準(zhǔn)則的轉(zhuǎn)捩公式[15]。我們要求近空間飛行器和航天飛行器一樣,基本上是光滑的,即實(shí)際粗糙度比上述實(shí)驗(yàn)的粗糙度小,則可以合理地認(rèn)為,開始發(fā)生轉(zhuǎn)捩的高度至少應(yīng)在40 km到45 km間。2015年又有幾篇有關(guān)美國(guó)HIFiRE-1計(jì)劃的分析試驗(yàn)報(bào)告,但其飛行高度值僅達(dá)到稍高于25 km。

    型號(hào)設(shè)計(jì)中總提出一個(gè)問(wèn)題:什么高度就可以不考慮湍流和轉(zhuǎn)捩影響,特別是近空間滑翔飛行器?從上面的分析來(lái)看,似乎可以認(rèn)為飛行高度在45km以上就可以不考慮轉(zhuǎn)捩和湍流問(wèn)題,但實(shí)際上問(wèn)題也許不是這么簡(jiǎn)單。

    首先,為了驗(yàn)證我們的結(jié)論,我們努力收集了國(guó)外飛行器的飛行結(jié)果,包括雙子星座、Apollo、美國(guó)各種航天飛機(jī)、蘇聯(lián)的飛船、英國(guó)的Hermes航天器、日本的小型飛船、美國(guó)的HTV-2以及中國(guó)的飛船,特別是哥倫比亞失事調(diào)查、復(fù)飛以及HTV-2前部的破壞等,都沒有發(fā)現(xiàn)是因?yàn)樵?0~50 km以上沒有考慮轉(zhuǎn)捩和湍流所導(dǎo)致的。最近美國(guó)正在改造世界最大的高超聲速風(fēng)洞(九號(hào)),其一個(gè)重要目的是能夠做邊界層自然轉(zhuǎn)捩的實(shí)驗(yàn),他們把對(duì)應(yīng)的飛行高度上限定為50 km。 以上是從已有的實(shí)際情況總結(jié)出來(lái)的結(jié)論,或許我們也可以從另一個(gè)角度對(duì)此問(wèn)題做一些分析,即從氣體稀薄程度來(lái)判斷是否會(huì)有湍流。

    本文作者曾在另一篇文章中[16]指出,對(duì)不同的問(wèn)題,判斷是否要考慮氣體的稀薄效應(yīng)要用不同的標(biāo)準(zhǔn),否則按早年的觀點(diǎn),高度在80 km以下的范圍,對(duì)飛行器的設(shè)計(jì)來(lái)說(shuō),都可以用連續(xù)介質(zhì)模型。2005年黃章峰等在他們的一篇文章[17]中,報(bào)告了他們所做的一個(gè)在30 km高空、以馬赫數(shù)4.5飛行的平板邊界層的湍流直接數(shù)值模擬。結(jié)果發(fā)現(xiàn),所得到的湍流和低速邊界層的湍流有幾乎完全相同的特征。例如,在壁面區(qū)存在低速條紋,其間距約為100個(gè)粘性長(zhǎng)度,此間距會(huì)隨離壁面的距離增加而逐漸加大。在湍流近壁區(qū),存在相干結(jié)構(gòu),其特征也和低速邊界層中的湍流相似。因此,我們可以假定在近壁區(qū)出現(xiàn)低速條紋,且條紋間距約為100粘性長(zhǎng)度,是湍流邊界層的普適現(xiàn)象。在判斷是否有湍流問(wèn)題時(shí),判斷氣體稀薄與否的參考長(zhǎng)度也許就應(yīng)取這一長(zhǎng)度。

    我們也可以先假定連續(xù)介質(zhì)模型可用,對(duì)在更高空處飛行的平板做類似的計(jì)算,但要計(jì)算多種工況則工作量較大,而采用以下的從層流解來(lái)估計(jì)則要方便得多。

    設(shè)在某一高度有一平板以0°迎角及某一速度飛行,則可按相似性解得到層流邊界層的解,從而得到相應(yīng)的粘性長(zhǎng)度。

    例如,對(duì)位于10 000 m和30 000 m高空的平板層流邊界層,采用相似性解,得到馬赫數(shù)為0.85和4時(shí),位于離前緣2 m處的粘性長(zhǎng)度l+見表1。

    表1 不同高空、不同馬赫數(shù)下的粘性長(zhǎng)度Table 1 Viscous length at different altitude and different Mach number

    由以上結(jié)果看,馬赫數(shù)和壁面的溫度條件對(duì)粘性長(zhǎng)度的影響不大,而高度則影響較大。高度從10 000 m增至30 000 m,l+就增加為原來(lái)的約4.8倍。

    黃章峰的DNS顯示,從層流轉(zhuǎn)捩為湍流時(shí),壁面摩擦系數(shù)增為原來(lái)的3~4倍。其他人對(duì)不可壓邊界層做的轉(zhuǎn)捩直接數(shù)值模擬,也得到類似的結(jié)果。如果取3,則對(duì)應(yīng)于湍流邊界層,粘性長(zhǎng)度單位約相當(dāng)于層流的(1/3)1/2≈0.58倍。因此,對(duì)湍流邊界層來(lái)說(shuō),粘性長(zhǎng)度和高度變化的規(guī)律也不變,在高度從10 000 m增加到30 000 m時(shí),也是增為原來(lái)的4.8倍,或高度每增加10 000 m,粘性長(zhǎng)度約增加約2.2倍。在比30 000 m更高的高空,這一增加的比例變化不會(huì)有量級(jí)的變化。

    但是從10 000 m高空開始,高度只要每增加10 000 m,大氣分子自由程就增為原來(lái)的約4.7倍。因此,高度每增加10 000 m,湍流邊界層近壁區(qū)的低速條紋間距和當(dāng)?shù)胤肿幼杂沙讨染鸵禐樵瓉?lái)的約2.2/4.7=0.47倍。當(dāng)這個(gè)比值降到一定的程度,低速條紋就無(wú)法維持了,即不會(huì)有湍流了。

    10 000 m高空分子自由程約為0.22μ,30 000 m高空時(shí)約為4.7μ。對(duì)300 K壁溫情況:馬赫數(shù)為4時(shí),在靠近壁面處,溫度升高很少,當(dāng)?shù)胤肿幼杂沙毯蛠?lái)流的分子自由程基本相同。而相應(yīng)地湍流邊界層低速條紋間距(0.58×100l+)則分別約為1.322 mm和6.322 mm,間距和分子自由程之比分別為6009和1345。

    如果馬赫數(shù)達(dá)20,則對(duì)定溫壁(300 K),在層流邊界層厚度的1/8處,溫度約為來(lái)流溫度的11倍[18](對(duì)湍流邊界層應(yīng)更高一些)。相應(yīng)地,對(duì)30 000 m高空,該處的分子自由程可達(dá)52μ,或0.052 mm。一個(gè)低速條紋間距對(duì)應(yīng)僅121個(gè)分子自由程。

    按一般湍流研究結(jié)果看,湍流邊界層近壁區(qū)的低速條紋是一種“大”結(jié)構(gòu),充分發(fā)展湍流中還應(yīng)有比它尺度至少小一個(gè)量級(jí)的小結(jié)構(gòu)。而要形成一個(gè)結(jié)構(gòu),其尺度又不能太小,必須包含足夠多的分子,或者說(shuō),其尺度應(yīng)該至少比分子自由程大一個(gè)量級(jí)。因此,低速條紋的間距至少應(yīng)該比分子自由程大兩個(gè)量級(jí),或約100倍。否則,至少連續(xù)介質(zhì)模型下的那種湍流是否能存在要存疑。但這也許還不能完全否定可能存在湍流。事實(shí)上,也許在一般的湍流和完全沒有湍流之間存在一個(gè)過(guò)渡區(qū),就像從層流到湍流存在一個(gè)轉(zhuǎn)捩過(guò)程那樣。在那個(gè)范圍內(nèi),可能存在既非層流又非充分發(fā)展的湍流。

    同樣地分析,到40 000 m高空時(shí),一個(gè)低速條紋僅對(duì)應(yīng)約56個(gè)分子自由程。從物理上講,應(yīng)已不能形成通常意義下的湍流。

    如果近空間飛行器采取的是耐高溫材料防熱,則壁面接近于絕熱情況。馬赫數(shù)為8時(shí),壁面溫度即可達(dá)來(lái)流溫度的約11倍(層流相似性解),即在40 km高空已不會(huì)有湍流。

    由于長(zhǎng)航程的高超聲速飛行器一般都要飛行于較高的高度,例如至少是40 km或更高,飛行馬赫數(shù)也較大,因此可能的確不需要考慮湍流和轉(zhuǎn)捩的問(wèn)題。

    但是,不能無(wú)條件地應(yīng)用這一結(jié)論。

    實(shí)際上,前面的分析都是針對(duì)0°迎角飛行的平板而做的,其隱含的假定就是邊界層內(nèi)的壓力和環(huán)境壓力相同。對(duì)有迎角平板,情況就不一樣了。

    例如,假定平板以速度v,迎角θ飛行,則按牛頓理論,平板迎風(fēng)面單位面積所受壓力約為:

    在海平面,空氣單位重量約為1.29kg/m3,大氣壓約為10 339kg/m2。在30 km高空,空氣單位重量和大氣壓力分別降為約0.0189kg/m3和119kg/m2。而在50 km高空,空氣單位重量約為0.000 89kg/m3(以上數(shù)據(jù)不是精確值,因?yàn)楦呖盏拇髿鉁囟?、密度和壓力都是不斷變?dòng)的)。如果馬赫數(shù)為15,則其速度約為4500 m/s,代入式(1),得

    如果要求和30 km高空處大氣壓相當(dāng),即約為119kg/m2,則Sin2θ=0.0647,或θ=14.7°。即迎角達(dá)到14.7°時(shí),其邊界層內(nèi)的壓力就和在30 km高空以同樣速度飛行但迎角為0°的平板一樣,也就有可能發(fā)生轉(zhuǎn)捩。

    對(duì)于像錐體那樣的飛行體,即使有迎角,在迎風(fēng)面的壓力升高顯然遠(yuǎn)小于平板的情況,因此在飛行試驗(yàn)中未發(fā)現(xiàn)高度超過(guò)30 km時(shí)發(fā)生轉(zhuǎn)捩也是可以理解的。但對(duì)有升力面的飛行器,則在升力面上,要根據(jù)實(shí)際邊界層中的壓力和溫度,看邊界層中的分子自由程和按連續(xù)介質(zhì)模型所得湍流邊界層中低速條紋間距(約100個(gè)粘性長(zhǎng)度)相比的值??梢钥紤]暫時(shí)將這一臨界值定為1%,大于它則多半不能維持湍流,也就不存在轉(zhuǎn)捩問(wèn)題了。

    3 結(jié)束語(yǔ)

    上面討論了研制近空間飛行器需要研究考慮的兩個(gè)空氣動(dòng)力學(xué)問(wèn)題。給出了它們的研究現(xiàn)狀,提出了為解決它們需要開展的研究及可能遇到的困難。由于問(wèn)題的復(fù)雜,也由于已有的研究成果很少,因此提出的看法和建議不能說(shuō)是完全正確的,也不能說(shuō)是很全面的。隨著研究的開展,可能還會(huì)發(fā)現(xiàn)新的問(wèn)題,或者發(fā)現(xiàn)文中的觀點(diǎn)不準(zhǔn)確甚至不正確,這都是正常的。如果我們的觀點(diǎn)或建議對(duì)今后的研究中能起到一定的促進(jìn)作用,我們就很滿意了。

    [1]Berin J J, Cummings R M. Critical hypersonic aerothermodunamic phenomena[J]. Annual Review of Fluid Mechanics, 2006, 38: 129-157.

    [2]周恒, 蘇彩虹, 張永明. 超聲速/高超聲速邊界層的轉(zhuǎn)捩機(jī)理及預(yù)測(cè)[M]. 北京: 科學(xué)出版社, 2015: 1-49, 88-96.

    [3]Pope A, Goin V L. High speed wind tunnel testing[M]. United Kingdom: Wiley, 1965.

    [4]Zhong X. Leading-edgereceptivity to free-stream disturbance waves for hypersonic flow over aparabola[J]. Journal of Fluid Mechanics, 2001, 441: 315-367.

    [5]Ma Y, Zhong X. Receptivity of asupersonic boundary layer overaflatplate. Part1. Waves tructures and interactions[J]. Journal of Fluid Mechanics, 2003, 488: 31-78.

    [6]Ma Y, Zhong X. Receptivity of asupersonic boundary layer overaflatplate. Part2. Recep-tivity to free stream sound[J]. Journal of Fluid Mechanics, 2003, 488: 79-121.

    [7]Ma Y, Zhong X. Receptivity of asupersonic boundary layer overaflatplate. Part3. Effects of different types of free-stream disturbances[J]. Journal of Fluid Mechanics, 2005, 532: 63-109.

    [8]Ma Y, Zhong X. Boundary-layer receptivity of Mach 7.99 flow overabluntcone to free-stream acoustic waves[J]. Journal of Fluid Mechanics, 2006, 556(1): 55-103.

    [9]Zhong X, Wang X. Directnumerical simulation on there ceptivity, instability and transition of hypersonic boundary layers[J]. Annual Review of Fluid Mechanics, 2012, 44(1): 527-561.

    [10]Balakumar P. Receptivity of asupersonic boundary layer to acoustic disturbances[J]. AIAA J., 2009, 47(5): 1069-1078.

    [11]Balakumar P, Kegerise M A. Receptivity of hypersonic boundary layers over straight and flared cones[J]. AIAA J., 2010, 53.

    [12]Wright R L, Zoby E V. Flight measuement of boundary layer transition of a 50 half angle cone at free stream Mach number of 20 (reentry)[R]. NASA TM X-2253, 1971.

    [13]Hayer D T, Herskovitz S B, Lennon J F, et al. An ablation technique for enhancing reentry antenna performance, flight test results[R]. FACRL-TR-74-0572, ADA012250, 1974.

    [14]Hayer D T, Herskovitz S B, Lennon J F, et al. Flight test data comparing electron attachment by ablation products and by liquid injection. AIAA-75-181[R]. Reston: AIAA, 1975.

    [15]Batt R Q, Legner H H. A review of roughness-induced nosetiptransition[J]. AIAA J., 1983, 21(1).

    [16]周恒, 張涵信. 空氣動(dòng)力學(xué)的新問(wèn)題[J]. 中國(guó)科學(xué)(中文版), 2015, 40(10): 104709.

    [17]Huang Zhangfeng, Zhou Heng, Luo Jisheng. Direct numerical simulation of a supersonic turbulent boundary layer on a flat plate and its analysis[J]. Science in China G, 2005, 48(5): 626-640.

    [18]Anderson J D Jr. Hypersonic and high-temperature gas dynamics[M]. Second Edition, AIAA Education Series, 2006, 293-296.

    Two problems in the transition and turbulence for near space
    hypersonic flying vehicles

    Zhou Heng1,*, Zhang Hanxin2

    (1.MechanicalSystemofTianjinUniversity,Tianjin300072,China;2.NationalLaboratoryforComputationalFluidDynamics,ChinaAerodynamicsResearchandDevelopmentCenter,Beijing100191,China)

    For the research and development of near space flying vehicles, also as the same for conventional flying vehiclessuch as airplanes, the problems of transition and turbulence of the boundary layers are two important issues must be taken into consideration. However, even for conventional flying vehicles, these two problems are still not fully resolved, even though the investigations have been lasted for more than 100 years already. For near space flying vehicles, not only the related aerodynamics has its own unsolved fundamental scientific problems, let along the problems of transition and turbulence. In this paper, two related problems are focused on: 1) In order to enhance our capability of predicting the transition of the boundary layer of a hypersonic flying vehicle, what kinds of research work should we do and what difficulties we might face? 2) Would it be possible that there would be no problem of turbulence for its boundary layer if the attitude of the flying vehicle is sufficiently high?

    near space flying vehicles; transition; turbulence; capability of predicting

    0258-1825(2017)02-0151-05

    2017-02-17;

    2017-03-23

    周恒*(1929-),男,中國(guó)科學(xué)院院士,研究方向:邊界層轉(zhuǎn)捩,流動(dòng)穩(wěn)定性非線性理論. E-mail: hzhou1@tju.edu.cn

    周恒, 張涵信. 有關(guān)近空間高超聲速飛行器邊界層轉(zhuǎn)捩和湍流的兩個(gè)問(wèn)題[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2017, 35(2): 151-155.

    10.7638/kqdlxxb-2017.0016. Zhou H, Zhang H X. Two problems in the transition and turbulence for near space hypersonic flying vehicles[J]. Acta Aerodynamica Sinica, 2017, 35(2): 151-155.

    V211.3

    A doi: 10.7638/kqdlxxb-2017.0016

    猜你喜歡
    馬赫數(shù)粘性邊界層
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    一類具有粘性項(xiàng)的擬線性拋物型方程組
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    載荷分布對(duì)可控?cái)U(kuò)散葉型性能的影響
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    帶粘性的波動(dòng)方程組解的逐點(diǎn)估計(jì)
    粘性非等熵流體方程平衡解的穩(wěn)定性
    一類具有邊界層性質(zhì)的二次奇攝動(dòng)邊值問(wèn)題
    家庭醫(yī)生增強(qiáng)基層首診粘性
    非特征邊界的MHD方程的邊界層
    免费观看人在逋| 亚洲成a人片在线一区二区| 欧美老熟妇乱子伦牲交| 成人特级黄色片久久久久久久| 欧美另类亚洲清纯唯美| 国产一区二区三区在线臀色熟女| 日本五十路高清| 中文字幕高清在线视频| 国产亚洲欧美98| 18禁裸乳无遮挡免费网站照片 | 搡老妇女老女人老熟妇| 亚洲 国产 在线| 91九色精品人成在线观看| 大码成人一级视频| 日韩高清综合在线| 黄色视频不卡| 91字幕亚洲| 亚洲最大成人中文| 久久久久精品国产欧美久久久| 国产av又大| 丝袜美腿诱惑在线| 国产精品二区激情视频| 亚洲av成人一区二区三| 亚洲人成电影免费在线| 成人手机av| 亚洲七黄色美女视频| 女人高潮潮喷娇喘18禁视频| 好看av亚洲va欧美ⅴa在| 91麻豆av在线| 女人被狂操c到高潮| 国产97色在线日韩免费| 久久人妻福利社区极品人妻图片| 男女做爰动态图高潮gif福利片 | 国产精品99久久99久久久不卡| 国产成人免费无遮挡视频| 精品国产一区二区久久| 亚洲aⅴ乱码一区二区在线播放 | www.精华液| 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 少妇熟女aⅴ在线视频| 亚洲在线自拍视频| 亚洲一卡2卡3卡4卡5卡精品中文| 狠狠狠狠99中文字幕| 99精品欧美一区二区三区四区| 亚洲少妇的诱惑av| 韩国av一区二区三区四区| 少妇的丰满在线观看| 亚洲人成电影免费在线| av视频在线观看入口| 美女高潮到喷水免费观看| 男女下面进入的视频免费午夜 | 伦理电影免费视频| 99久久精品国产亚洲精品| 视频在线观看一区二区三区| 淫秽高清视频在线观看| 午夜免费观看网址| 欧美av亚洲av综合av国产av| 精品人妻1区二区| 久久久久久免费高清国产稀缺| 老司机深夜福利视频在线观看| 色综合欧美亚洲国产小说| 亚洲色图 男人天堂 中文字幕| 亚洲人成电影免费在线| 欧美黄色淫秽网站| 日日干狠狠操夜夜爽| 1024香蕉在线观看| av天堂在线播放| 午夜免费成人在线视频| 好男人电影高清在线观看| 亚洲成人精品中文字幕电影| 777久久人妻少妇嫩草av网站| 欧美不卡视频在线免费观看 | 啦啦啦 在线观看视频| 国产精品一区二区三区四区久久 | 久久精品亚洲精品国产色婷小说| 黄色毛片三级朝国网站| 国内精品久久久久精免费| 级片在线观看| www国产在线视频色| 麻豆av在线久日| 伦理电影免费视频| xxx96com| 成年女人毛片免费观看观看9| 久久香蕉国产精品| 美女午夜性视频免费| av中文乱码字幕在线| 涩涩av久久男人的天堂| 一边摸一边做爽爽视频免费| 又黄又粗又硬又大视频| 黄色a级毛片大全视频| 校园春色视频在线观看| 在线播放国产精品三级| 欧美丝袜亚洲另类 | 午夜免费鲁丝| 中文字幕久久专区| 日韩有码中文字幕| 国产精品久久久久久人妻精品电影| 亚洲熟妇熟女久久| 国产高清激情床上av| 午夜福利一区二区在线看| 成人国产综合亚洲| 午夜日韩欧美国产| 成人永久免费在线观看视频| 亚洲国产精品久久男人天堂| 母亲3免费完整高清在线观看| 久久久精品欧美日韩精品| 欧美av亚洲av综合av国产av| netflix在线观看网站| 精品久久久精品久久久| 成年版毛片免费区| 亚洲五月色婷婷综合| 精品久久久久久久人妻蜜臀av | 在线观看舔阴道视频| 亚洲国产中文字幕在线视频| 免费不卡黄色视频| 如日韩欧美国产精品一区二区三区| 久久精品国产亚洲av高清一级| 亚洲成人国产一区在线观看| 精品久久久久久,| 国产成人系列免费观看| 男人舔女人的私密视频| 国产一级毛片七仙女欲春2 | 一本大道久久a久久精品| 亚洲色图av天堂| 精品国内亚洲2022精品成人| 国内精品久久久久精免费| 制服人妻中文乱码| 正在播放国产对白刺激| 成人国产一区最新在线观看| 自线自在国产av| 国产又色又爽无遮挡免费看| av天堂在线播放| 侵犯人妻中文字幕一二三四区| 成人国产综合亚洲| 欧美日韩瑟瑟在线播放| 天天一区二区日本电影三级 | 免费高清视频大片| 免费在线观看完整版高清| 99热只有精品国产| 久久精品国产清高在天天线| 一进一出好大好爽视频| 亚洲视频免费观看视频| 国产精品久久久人人做人人爽| 亚洲片人在线观看| 美女大奶头视频| 国产又爽黄色视频| 亚洲国产精品999在线| 国产野战对白在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品亚洲一级av第二区| 欧美成人午夜精品| 亚洲一区二区三区色噜噜| 亚洲av成人不卡在线观看播放网| 国产一级毛片七仙女欲春2 | 成人18禁高潮啪啪吃奶动态图| 国产成人av教育| 校园春色视频在线观看| 九色国产91popny在线| 欧美中文综合在线视频| 亚洲精品国产一区二区精华液| 91麻豆av在线| 如日韩欧美国产精品一区二区三区| 亚洲人成电影免费在线| 国产av又大| 亚洲成国产人片在线观看| 亚洲国产中文字幕在线视频| 黄片播放在线免费| 最新美女视频免费是黄的| 波多野结衣一区麻豆| 十八禁人妻一区二区| 久久午夜亚洲精品久久| 日本撒尿小便嘘嘘汇集6| 无人区码免费观看不卡| 亚洲精品久久国产高清桃花| a在线观看视频网站| 如日韩欧美国产精品一区二区三区| 久9热在线精品视频| tocl精华| 三级毛片av免费| 日韩一卡2卡3卡4卡2021年| 在线观看免费午夜福利视频| 非洲黑人性xxxx精品又粗又长| 精品无人区乱码1区二区| 非洲黑人性xxxx精品又粗又长| 亚洲成a人片在线一区二区| 国产区一区二久久| 成熟少妇高潮喷水视频| 国产熟女午夜一区二区三区| 电影成人av| 首页视频小说图片口味搜索| 色尼玛亚洲综合影院| 日日夜夜操网爽| 日韩欧美国产在线观看| 91大片在线观看| 国产精品一区二区免费欧美| 香蕉久久夜色| 身体一侧抽搐| 一进一出抽搐gif免费好疼| www日本在线高清视频| 老司机在亚洲福利影院| 精品福利观看| 国产精品一区二区精品视频观看| 亚洲人成网站在线播放欧美日韩| 成人av一区二区三区在线看| 99国产极品粉嫩在线观看| 99久久久亚洲精品蜜臀av| 大型黄色视频在线免费观看| 亚洲精品久久成人aⅴ小说| 国产免费av片在线观看野外av| 黄色女人牲交| 岛国视频午夜一区免费看| 黄色视频不卡| 欧美日韩亚洲综合一区二区三区_| av福利片在线| 午夜老司机福利片| 在线av久久热| 午夜老司机福利片| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美日韩高清在线视频| 欧美最黄视频在线播放免费| 成人18禁高潮啪啪吃奶动态图| 一区二区三区国产精品乱码| 老司机在亚洲福利影院| 日本黄色视频三级网站网址| 亚洲成人精品中文字幕电影| 日韩欧美一区视频在线观看| 欧美午夜高清在线| 欧美日韩一级在线毛片| 国产精品av久久久久免费| 久久精品成人免费网站| 禁无遮挡网站| 无人区码免费观看不卡| 婷婷丁香在线五月| 精品高清国产在线一区| 欧美不卡视频在线免费观看 | 老汉色av国产亚洲站长工具| 两个人视频免费观看高清| 在线观看免费日韩欧美大片| 久9热在线精品视频| 亚洲成国产人片在线观看| 国产精品一区二区免费欧美| 熟妇人妻久久中文字幕3abv| 丁香欧美五月| 午夜a级毛片| 久久久久久大精品| 一区二区三区精品91| 国产区一区二久久| 久久精品国产综合久久久| 国产乱人伦免费视频| 亚洲全国av大片| 久久精品国产清高在天天线| 老鸭窝网址在线观看| 国产成人欧美在线观看| 午夜日韩欧美国产| 999久久久精品免费观看国产| 在线永久观看黄色视频| 久久中文字幕人妻熟女| 国产亚洲av嫩草精品影院| 精品久久久久久成人av| 看片在线看免费视频| 成人永久免费在线观看视频| 亚洲美女黄片视频| 成人亚洲精品一区在线观看| 制服诱惑二区| 国产在线精品亚洲第一网站| 侵犯人妻中文字幕一二三四区| 国产一卡二卡三卡精品| 久久午夜亚洲精品久久| 精品国产乱子伦一区二区三区| 91老司机精品| 在线国产一区二区在线| 身体一侧抽搐| 久久精品影院6| 嫁个100分男人电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 老司机靠b影院| 亚洲第一青青草原| 人人妻人人澡人人看| 久久久久久久久免费视频了| 九色亚洲精品在线播放| 久久亚洲精品不卡| 青草久久国产| 午夜久久久在线观看| 日韩中文字幕欧美一区二区| 国产成年人精品一区二区| 国内精品久久久久久久电影| 曰老女人黄片| 免费高清视频大片| 国产成人欧美在线观看| 亚洲成人免费电影在线观看| 淫秽高清视频在线观看| 中文字幕色久视频| 九色亚洲精品在线播放| 精品久久久久久,| 美女高潮到喷水免费观看| 日韩高清综合在线| 精品欧美一区二区三区在线| 国产成人免费无遮挡视频| 老汉色∧v一级毛片| av视频免费观看在线观看| 成在线人永久免费视频| 真人一进一出gif抽搐免费| 中出人妻视频一区二区| 脱女人内裤的视频| 午夜两性在线视频| 亚洲少妇的诱惑av| 国产一区二区三区综合在线观看| 国产亚洲精品久久久久5区| 亚洲人成电影观看| 国产97色在线日韩免费| 欧美大码av| 免费一级毛片在线播放高清视频 | 国产免费av片在线观看野外av| 长腿黑丝高跟| 亚洲精品中文字幕一二三四区| 99久久精品国产亚洲精品| 高潮久久久久久久久久久不卡| av免费在线观看网站| 国产精品1区2区在线观看.| 少妇粗大呻吟视频| 亚洲激情在线av| 日韩av在线大香蕉| 国产97色在线日韩免费| 中文字幕av电影在线播放| 国产一区二区三区视频了| 国产精品久久久av美女十八| 亚洲av成人av| 丝袜人妻中文字幕| 国产精品 欧美亚洲| 亚洲欧美激情在线| 亚洲国产高清在线一区二区三 | 中文字幕高清在线视频| 亚洲三区欧美一区| 亚洲第一青青草原| 午夜福利在线观看吧| 成人亚洲精品av一区二区| 日日爽夜夜爽网站| 一进一出抽搐gif免费好疼| 一区二区日韩欧美中文字幕| 麻豆av在线久日| 国产人伦9x9x在线观看| 侵犯人妻中文字幕一二三四区| 午夜a级毛片| 满18在线观看网站| 免费不卡黄色视频| 亚洲精华国产精华精| 免费不卡黄色视频| 真人做人爱边吃奶动态| 亚洲精品在线美女| 黄片小视频在线播放| 国产色视频综合| 正在播放国产对白刺激| 变态另类成人亚洲欧美熟女 | 国产精品 欧美亚洲| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利18| 一夜夜www| 色综合欧美亚洲国产小说| 欧美绝顶高潮抽搐喷水| 久热这里只有精品99| 一级,二级,三级黄色视频| 欧美久久黑人一区二区| 九色亚洲精品在线播放| 亚洲激情在线av| 国产熟女午夜一区二区三区| 日日夜夜操网爽| 纯流量卡能插随身wifi吗| 女生性感内裤真人,穿戴方法视频| 国产精华一区二区三区| 国产成人av激情在线播放| 亚洲电影在线观看av| 99热只有精品国产| 久久九九热精品免费| 中文字幕最新亚洲高清| 亚洲五月色婷婷综合| 丰满的人妻完整版| 亚洲人成电影观看| 欧美久久黑人一区二区| 亚洲精品国产区一区二| 国产一区二区三区综合在线观看| 99在线视频只有这里精品首页| 中文字幕久久专区| 少妇粗大呻吟视频| 成人免费观看视频高清| 日韩有码中文字幕| 丝袜美足系列| 国产精品影院久久| 色播亚洲综合网| 欧美中文综合在线视频| 亚洲电影在线观看av| 十八禁网站免费在线| 成人手机av| 看免费av毛片| www.自偷自拍.com| 51午夜福利影视在线观看| 狂野欧美激情性xxxx| 老熟妇乱子伦视频在线观看| 亚洲片人在线观看| 久久精品成人免费网站| 国产在线精品亚洲第一网站| 老司机靠b影院| 免费看美女性在线毛片视频| 国产成人欧美在线观看| 美女国产高潮福利片在线看| 亚洲片人在线观看| 在线观看免费日韩欧美大片| 日日夜夜操网爽| 久久久久久人人人人人| 色综合站精品国产| 国产精品秋霞免费鲁丝片| 午夜福利影视在线免费观看| 国产xxxxx性猛交| 亚洲精华国产精华精| 欧美久久黑人一区二区| 乱人伦中国视频| 欧美激情高清一区二区三区| 欧美中文综合在线视频| 亚洲 欧美 日韩 在线 免费| 国产一区二区三区综合在线观看| 国产精品秋霞免费鲁丝片| 欧美绝顶高潮抽搐喷水| 久久亚洲精品不卡| 91国产中文字幕| 亚洲中文日韩欧美视频| 村上凉子中文字幕在线| av网站免费在线观看视频| 久久人妻av系列| 午夜福利高清视频| 免费看美女性在线毛片视频| 可以免费在线观看a视频的电影网站| 老熟妇乱子伦视频在线观看| 久久天堂一区二区三区四区| av福利片在线| 亚洲va日本ⅴa欧美va伊人久久| 国产av一区二区精品久久| 国产精品美女特级片免费视频播放器 | 午夜福利18| 亚洲精品久久国产高清桃花| 99国产精品99久久久久| 一区福利在线观看| 三级毛片av免费| 亚洲av五月六月丁香网| 男女床上黄色一级片免费看| 性少妇av在线| 一边摸一边做爽爽视频免费| 成在线人永久免费视频| 亚洲免费av在线视频| 日韩精品中文字幕看吧| 欧美黄色片欧美黄色片| 欧美丝袜亚洲另类 | 69精品国产乱码久久久| 女性被躁到高潮视频| 天堂√8在线中文| 首页视频小说图片口味搜索| 在线观看免费视频日本深夜| 亚洲中文av在线| 人妻久久中文字幕网| 一区二区三区高清视频在线| 亚洲片人在线观看| 久久中文看片网| 久久久国产成人精品二区| 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频 | 久久人人97超碰香蕉20202| 免费人成视频x8x8入口观看| 一级毛片精品| 啦啦啦免费观看视频1| 久久香蕉激情| 亚洲色图av天堂| 国产午夜精品久久久久久| 午夜激情av网站| 久久精品国产亚洲av香蕉五月| 欧美日韩一级在线毛片| 午夜激情av网站| 亚洲av电影在线进入| 久久影院123| 亚洲国产精品sss在线观看| 真人一进一出gif抽搐免费| 欧美精品亚洲一区二区| 国产精品av久久久久免费| 国内精品久久久久精免费| 精品国产美女av久久久久小说| 女人被狂操c到高潮| 少妇粗大呻吟视频| 在线观看午夜福利视频| av电影中文网址| 日本精品一区二区三区蜜桃| 99久久综合精品五月天人人| 亚洲专区中文字幕在线| 1024视频免费在线观看| 国产精品亚洲美女久久久| 中文字幕人妻丝袜一区二区| 午夜福利成人在线免费观看| 韩国精品一区二区三区| 国产亚洲av嫩草精品影院| 757午夜福利合集在线观看| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 精品人妻在线不人妻| 亚洲人成伊人成综合网2020| 丝袜在线中文字幕| 亚洲av熟女| 亚洲国产日韩欧美精品在线观看 | av福利片在线| 亚洲五月色婷婷综合| 中文字幕另类日韩欧美亚洲嫩草| 看片在线看免费视频| 黄色女人牲交| 欧美色欧美亚洲另类二区 | 国产精品久久久久久人妻精品电影| 亚洲中文日韩欧美视频| 男女床上黄色一级片免费看| 欧美中文综合在线视频| 国产三级在线视频| 美女 人体艺术 gogo| 在线观看www视频免费| 亚洲黑人精品在线| 好看av亚洲va欧美ⅴa在| 欧美日本中文国产一区发布| 亚洲欧洲精品一区二区精品久久久| 男女午夜视频在线观看| 日韩有码中文字幕| 色哟哟哟哟哟哟| 大香蕉久久成人网| 国产精品免费视频内射| 亚洲aⅴ乱码一区二区在线播放 | 国产伦人伦偷精品视频| 中文字幕av电影在线播放| 精品乱码久久久久久99久播| 成人18禁在线播放| 精品一品国产午夜福利视频| 国产区一区二久久| 国产av在哪里看| 女警被强在线播放| 成人av一区二区三区在线看| 午夜精品在线福利| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久人妻精品电影| 久久久久久久久中文| 日韩三级视频一区二区三区| 国产日韩一区二区三区精品不卡| 国产在线观看jvid| aaaaa片日本免费| 亚洲全国av大片| 日韩高清综合在线| 亚洲五月色婷婷综合| 一级片免费观看大全| 国产精品国产高清国产av| 国产成人精品在线电影| 一级毛片精品| 国产精品精品国产色婷婷| 无限看片的www在线观看| 亚洲人成77777在线视频| 精品卡一卡二卡四卡免费| 欧美日韩亚洲综合一区二区三区_| 淫妇啪啪啪对白视频| 99热只有精品国产| 国产人伦9x9x在线观看| 两个人免费观看高清视频| 成人三级做爰电影| 久久久国产欧美日韩av| 亚洲成a人片在线一区二区| 日本免费一区二区三区高清不卡 | 精品国产亚洲在线| 18美女黄网站色大片免费观看| 亚洲国产精品合色在线| 最近最新中文字幕大全免费视频| 咕卡用的链子| 大陆偷拍与自拍| 色综合婷婷激情| 国产一区二区三区在线臀色熟女| 黑人操中国人逼视频| 女生性感内裤真人,穿戴方法视频| 国产成人免费无遮挡视频| 国产欧美日韩精品亚洲av| 久久久久精品国产欧美久久久| 国产不卡一卡二| 大码成人一级视频| 成人三级黄色视频| 色av中文字幕| 欧美黑人精品巨大| 色尼玛亚洲综合影院| 国产免费av片在线观看野外av| 国产片内射在线| 久久精品国产清高在天天线| 国产人伦9x9x在线观看| 久久中文看片网| 国内精品久久久久精免费| 国产黄a三级三级三级人| 日日爽夜夜爽网站| 女性生殖器流出的白浆| 热re99久久国产66热| 色综合站精品国产| 欧美日本视频| 亚洲七黄色美女视频| 国产乱人伦免费视频| 18禁裸乳无遮挡免费网站照片 | 国产精品一区二区免费欧美| 亚洲av美国av| 日日摸夜夜添夜夜添小说| 嫩草影院精品99| 日韩精品中文字幕看吧| 一边摸一边抽搐一进一小说| 欧美成人午夜精品| 久久午夜亚洲精品久久| 韩国精品一区二区三区| www国产在线视频色| 亚洲精品久久成人aⅴ小说| 久久久久精品国产欧美久久久| 免费在线观看完整版高清| 成人欧美大片| 99国产综合亚洲精品| 亚洲最大成人中文| 啪啪无遮挡十八禁网站|