• <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方程的邊界層
    日韩欧美国产一区二区入口| 听说在线观看完整版免费高清| 国产av在哪里看| 精品久久久久久久人妻蜜臀av| 九色国产91popny在线| 欧美极品一区二区三区四区| 美女cb高潮喷水在线观看 | 国产日本99.免费观看| 日本 欧美在线| 国产av不卡久久| 日本精品一区二区三区蜜桃| 国产真人三级小视频在线观看| 国产精品一区二区精品视频观看| 国产aⅴ精品一区二区三区波| 18美女黄网站色大片免费观看| 色综合亚洲欧美另类图片| 婷婷精品国产亚洲av| 一本精品99久久精品77| 1024手机看黄色片| 国产精品自产拍在线观看55亚洲| 99re在线观看精品视频| 女同久久另类99精品国产91| 亚洲欧洲精品一区二区精品久久久| 男人的好看免费观看在线视频| 99国产综合亚洲精品| 男女下面进入的视频免费午夜| 国产三级中文精品| 老熟妇乱子伦视频在线观看| 精品一区二区三区av网在线观看| 一本一本综合久久| 淫秽高清视频在线观看| 精品人妻1区二区| 国产精品一区二区三区四区免费观看 | 成在线人永久免费视频| 女人高潮潮喷娇喘18禁视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲va日本ⅴa欧美va伊人久久| 日本黄色片子视频| 校园春色视频在线观看| 亚洲国产欧美网| 身体一侧抽搐| 精品国产美女av久久久久小说| 亚洲熟妇中文字幕五十中出| а√天堂www在线а√下载| 无人区码免费观看不卡| 午夜福利成人在线免费观看| 真人做人爱边吃奶动态| 动漫黄色视频在线观看| 九九久久精品国产亚洲av麻豆 | 久久这里只有精品19| 两个人看的免费小视频| 国产精品乱码一区二三区的特点| 国产av在哪里看| 十八禁网站免费在线| 一夜夜www| 美女cb高潮喷水在线观看 | 国产精品香港三级国产av潘金莲| 嫩草影院精品99| 亚洲精华国产精华精| 99久久国产精品久久久| 日韩大尺度精品在线看网址| 免费av不卡在线播放| 嫩草影视91久久| 色综合亚洲欧美另类图片| 午夜亚洲福利在线播放| 午夜激情福利司机影院| 很黄的视频免费| 老司机福利观看| 日日摸夜夜添夜夜添小说| 免费无遮挡裸体视频| 久久香蕉精品热| 一级a爱片免费观看的视频| 亚洲精品粉嫩美女一区| 亚洲av日韩精品久久久久久密| 成人无遮挡网站| 成人三级做爰电影| 在线播放国产精品三级| 高潮久久久久久久久久久不卡| 黄色日韩在线| 97碰自拍视频| 一a级毛片在线观看| 亚洲无线在线观看| 日本一二三区视频观看| 午夜福利成人在线免费观看| 韩国av一区二区三区四区| 老司机深夜福利视频在线观看| www日本在线高清视频| 国产精品野战在线观看| 中文字幕久久专区| 制服人妻中文乱码| 小蜜桃在线观看免费完整版高清| 亚洲专区国产一区二区| 日韩免费av在线播放| 丁香欧美五月| АⅤ资源中文在线天堂| 最近最新免费中文字幕在线| 一级作爱视频免费观看| 熟女人妻精品中文字幕| 一进一出抽搐动态| 午夜精品在线福利| 级片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 老司机深夜福利视频在线观看| 韩国av一区二区三区四区| 国产精品电影一区二区三区| 免费在线观看视频国产中文字幕亚洲| 麻豆一二三区av精品| 99精品在免费线老司机午夜| 久久久久久久久免费视频了| 成人18禁在线播放| 丰满人妻熟妇乱又伦精品不卡| 色综合婷婷激情| 久久精品91蜜桃| 男人舔女人下体高潮全视频| 99在线人妻在线中文字幕| 香蕉国产在线看| 狂野欧美激情性xxxx| 一级毛片高清免费大全| 精品电影一区二区在线| 丰满的人妻完整版| 久久天堂一区二区三区四区| 久久久成人免费电影| 在线永久观看黄色视频| 精品国产美女av久久久久小说| 久久久久久大精品| 午夜福利视频1000在线观看| 国产v大片淫在线免费观看| 久久国产精品人妻蜜桃| 一本精品99久久精品77| 日本在线视频免费播放| 国产av一区在线观看免费| 色综合欧美亚洲国产小说| 国产精品永久免费网站| 日本在线视频免费播放| 日韩国内少妇激情av| 免费av不卡在线播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩精品网址| 欧美激情久久久久久爽电影| 国产精品一区二区三区四区免费观看 | 国产精品香港三级国产av潘金莲| 国产不卡一卡二| 巨乳人妻的诱惑在线观看| www日本黄色视频网| 99精品欧美一区二区三区四区| 色综合站精品国产| 国产精品av久久久久免费| 精品久久久久久久毛片微露脸| 19禁男女啪啪无遮挡网站| 美女黄网站色视频| 五月玫瑰六月丁香| 久久久久久国产a免费观看| 搡老妇女老女人老熟妇| 一边摸一边抽搐一进一小说| 在线观看日韩欧美| 三级毛片av免费| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看视频国产中文字幕亚洲| 久久国产乱子伦精品免费另类| 午夜福利免费观看在线| 午夜精品在线福利| 最新美女视频免费是黄的| 国产成人aa在线观看| 国产成人啪精品午夜网站| 国产单亲对白刺激| 国产成人欧美在线观看| 成人亚洲精品av一区二区| 99热这里只有精品一区 | 成人精品一区二区免费| 真实男女啪啪啪动态图| 国产av麻豆久久久久久久| 一区二区三区国产精品乱码| 日韩欧美在线乱码| 国产黄片美女视频| 国产三级在线视频| 久久久成人免费电影| 久久久水蜜桃国产精品网| 亚洲av日韩精品久久久久久密| 99精品欧美一区二区三区四区| 国产精品爽爽va在线观看网站| 亚洲九九香蕉| 亚洲国产精品成人综合色| 色精品久久人妻99蜜桃| 嫩草影院入口| 男人和女人高潮做爰伦理| 欧美色视频一区免费| 热99re8久久精品国产| 国产成人啪精品午夜网站| 国产 一区 欧美 日韩| 不卡av一区二区三区| 97人妻精品一区二区三区麻豆| а√天堂www在线а√下载| 熟女人妻精品中文字幕| 小蜜桃在线观看免费完整版高清| 老熟妇仑乱视频hdxx| 伊人久久大香线蕉亚洲五| 国产主播在线观看一区二区| 国产精品自产拍在线观看55亚洲| 最近最新中文字幕大全免费视频| 欧美另类亚洲清纯唯美| 国产黄片美女视频| 欧美午夜高清在线| 一级作爱视频免费观看| 亚洲成av人片在线播放无| 国产午夜福利久久久久久| 男女下面进入的视频免费午夜| 国产成人精品无人区| 亚洲成a人片在线一区二区| 午夜精品一区二区三区免费看| 欧美三级亚洲精品| 99久久精品热视频| 亚洲av日韩精品久久久久久密| 日本一二三区视频观看| 国产久久久一区二区三区| 国产免费男女视频| 中文亚洲av片在线观看爽| 国内精品美女久久久久久| 久久久国产欧美日韩av| 国产精品免费一区二区三区在线| 嫩草影院精品99| a级毛片在线看网站| 国产亚洲av嫩草精品影院| 女同久久另类99精品国产91| 高潮久久久久久久久久久不卡| 脱女人内裤的视频| 香蕉久久夜色| 在线观看66精品国产| 伊人久久大香线蕉亚洲五| 最新中文字幕久久久久 | 两性午夜刺激爽爽歪歪视频在线观看| 麻豆一二三区av精品| 不卡av一区二区三区| 国产精品久久久久久久电影 | 又紧又爽又黄一区二区| 99久久久亚洲精品蜜臀av| 美女大奶头视频| 五月玫瑰六月丁香| 岛国在线免费视频观看| 日本一二三区视频观看| 狂野欧美激情性xxxx| 脱女人内裤的视频| 国产精品国产高清国产av| 欧美一区二区精品小视频在线| 大型黄色视频在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 天堂动漫精品| 性色av乱码一区二区三区2| 久久午夜亚洲精品久久| 国内精品久久久久久久电影| 亚洲国产精品久久男人天堂| 久久精品国产综合久久久| 亚洲专区国产一区二区| 国产v大片淫在线免费观看| 小蜜桃在线观看免费完整版高清| 熟女电影av网| 国产成人系列免费观看| 国产一区二区激情短视频| 99在线视频只有这里精品首页| 真人做人爱边吃奶动态| 久久久久亚洲av毛片大全| 国产伦人伦偷精品视频| 丰满人妻熟妇乱又伦精品不卡| 人妻夜夜爽99麻豆av| 十八禁网站免费在线| 99久久精品国产亚洲精品| 欧美日本视频| av国产免费在线观看| 两个人视频免费观看高清| 精品日产1卡2卡| 国产激情偷乱视频一区二区| 国产精品99久久99久久久不卡| 欧美日本亚洲视频在线播放| 噜噜噜噜噜久久久久久91| 男女视频在线观看网站免费| 成年免费大片在线观看| 国内精品久久久久精免费| 亚洲av五月六月丁香网| 午夜免费成人在线视频| 亚洲熟妇中文字幕五十中出| 日本 欧美在线| 美女 人体艺术 gogo| 级片在线观看| 国产免费男女视频| 又粗又爽又猛毛片免费看| 日韩欧美一区二区三区在线观看| 我要搜黄色片| 三级毛片av免费| 午夜精品一区二区三区免费看| 黄色日韩在线| 他把我摸到了高潮在线观看| 欧美激情在线99| 夜夜夜夜夜久久久久| 美女被艹到高潮喷水动态| 欧美三级亚洲精品| 少妇的丰满在线观看| 亚洲 欧美一区二区三区| 我的老师免费观看完整版| 亚洲性夜色夜夜综合| 国产69精品久久久久777片 | 亚洲国产色片| 亚洲男人的天堂狠狠| 欧美不卡视频在线免费观看| 国产精品久久久久久亚洲av鲁大| 亚洲激情在线av| 岛国在线免费视频观看| 桃色一区二区三区在线观看| 成人一区二区视频在线观看| 国产一级毛片七仙女欲春2| 国内揄拍国产精品人妻在线| 免费看美女性在线毛片视频| 国产真实乱freesex| 国产久久久一区二区三区| 国产伦一二天堂av在线观看| 黄色 视频免费看| 午夜a级毛片| 最近最新中文字幕大全免费视频| 国产69精品久久久久777片 | 久久99热这里只有精品18| 日日干狠狠操夜夜爽| 久久这里只有精品19| 99热精品在线国产| 国产乱人视频| 亚洲av日韩精品久久久久久密| av视频在线观看入口| 99久久国产精品久久久| 亚洲人成网站高清观看| 老熟妇乱子伦视频在线观看| 日韩精品青青久久久久久| 成人特级av手机在线观看| 国产欧美日韩一区二区精品| 中文在线观看免费www的网站| 久久国产精品影院| 男人舔女人下体高潮全视频| а√天堂www在线а√下载| 桃色一区二区三区在线观看| 成人一区二区视频在线观看| 亚洲精品在线美女| 日韩欧美三级三区| 欧美在线一区亚洲| 欧美一级毛片孕妇| 国产精品99久久99久久久不卡| 日韩 欧美 亚洲 中文字幕| 日本免费a在线| 国产精品亚洲av一区麻豆| 一区二区三区国产精品乱码| 日韩人妻高清精品专区| 午夜福利免费观看在线| 狂野欧美白嫩少妇大欣赏| 一区二区三区国产精品乱码| 久久久久久人人人人人| 免费高清视频大片| 美女 人体艺术 gogo| 激情在线观看视频在线高清| 久久人人精品亚洲av| 一个人免费在线观看的高清视频| 三级毛片av免费| 在线国产一区二区在线| 好男人电影高清在线观看| 黄频高清免费视频| 精品99又大又爽又粗少妇毛片 | 亚洲av成人不卡在线观看播放网| 久久亚洲精品不卡| 嫩草影院入口| 精品久久久久久久久久久久久| h日本视频在线播放| 国产免费av片在线观看野外av| 免费看十八禁软件| 怎么达到女性高潮| 亚洲中文字幕一区二区三区有码在线看 | 一级a爱片免费观看的视频| 91在线精品国自产拍蜜月 | 一卡2卡三卡四卡精品乱码亚洲| 全区人妻精品视频| 亚洲av电影在线进入| av黄色大香蕉| 午夜影院日韩av| 69av精品久久久久久| 亚洲午夜理论影院| 国产精品免费一区二区三区在线| 久久久国产成人精品二区| 熟妇人妻久久中文字幕3abv| 久久亚洲精品不卡| 欧美xxxx黑人xx丫x性爽| 两个人看的免费小视频| 亚洲中文字幕一区二区三区有码在线看 | 琪琪午夜伦伦电影理论片6080| 很黄的视频免费| 黄色成人免费大全| 女生性感内裤真人,穿戴方法视频| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 日韩欧美免费精品| 亚洲国产精品久久男人天堂| 久久香蕉精品热| 国产一区二区在线av高清观看| bbb黄色大片| 国产亚洲精品综合一区在线观看| 午夜福利18| cao死你这个sao货| 一级毛片女人18水好多| 国产三级在线视频| 国产高清有码在线观看视频| 国产精品 欧美亚洲| 成人特级黄色片久久久久久久| 十八禁网站免费在线| 午夜激情福利司机影院| 欧美日韩国产亚洲二区| 美女午夜性视频免费| 在线观看日韩欧美| 黄色日韩在线| 国产精品精品国产色婷婷| 欧美一级毛片孕妇| 国产精品亚洲一级av第二区| www.999成人在线观看| 精品欧美国产一区二区三| 成人精品一区二区免费| av欧美777| 禁无遮挡网站| 麻豆一二三区av精品| 搡老岳熟女国产| 午夜福利视频1000在线观看| 一级毛片高清免费大全| 五月伊人婷婷丁香| 日韩高清综合在线| 亚洲av成人一区二区三| 少妇人妻一区二区三区视频| 午夜福利在线观看吧| 欧美日本视频| 亚洲精品美女久久av网站| 国产精品精品国产色婷婷| 成人亚洲精品av一区二区| 成年女人看的毛片在线观看| a级毛片在线看网站| 无限看片的www在线观看| 亚洲av美国av| 99精品久久久久人妻精品| 国产免费男女视频| 久久精品国产清高在天天线| 国产麻豆成人av免费视频| 国内精品一区二区在线观看| 国产精品日韩av在线免费观看| 欧美国产日韩亚洲一区| 亚洲成人久久性| 亚洲av电影在线进入| 亚洲精品456在线播放app | 亚洲午夜理论影院| av欧美777| 国产淫片久久久久久久久 | 天堂动漫精品| 午夜久久久久精精品| 巨乳人妻的诱惑在线观看| x7x7x7水蜜桃| 久久精品人妻少妇| 亚洲欧美日韩高清专用| 国产精品亚洲一级av第二区| 综合色av麻豆| 成人一区二区视频在线观看| 丰满的人妻完整版| 国产久久久一区二区三区| 男人舔女人下体高潮全视频| 老司机午夜福利在线观看视频| 国内少妇人妻偷人精品xxx网站 | 久久精品91蜜桃| 热99re8久久精品国产| 久久国产精品影院| 日本 欧美在线| 亚洲午夜精品一区,二区,三区| 亚洲av第一区精品v没综合| 亚洲无线观看免费| 欧美黑人欧美精品刺激| 国产 一区 欧美 日韩| 亚洲片人在线观看| 婷婷丁香在线五月| 波多野结衣高清无吗| 在线观看日韩欧美| 女同久久另类99精品国产91| 两个人视频免费观看高清| 人人妻人人看人人澡| 99热只有精品国产| 桃色一区二区三区在线观看| 国产视频一区二区在线看| 噜噜噜噜噜久久久久久91| 国产 一区 欧美 日韩| 日韩成人在线观看一区二区三区| 俺也久久电影网| 免费在线观看视频国产中文字幕亚洲| 亚洲精品粉嫩美女一区| 黑人欧美特级aaaaaa片| 国产一区二区在线观看日韩 | 在线观看美女被高潮喷水网站 | av中文乱码字幕在线| 欧美成狂野欧美在线观看| 国产精品99久久久久久久久| 18禁国产床啪视频网站| 美女高潮的动态| av在线蜜桃| 精品久久久久久久末码| 欧美乱码精品一区二区三区| 国产私拍福利视频在线观看| 性欧美人与动物交配| ponron亚洲| 男人和女人高潮做爰伦理| 欧美激情久久久久久爽电影| 亚洲狠狠婷婷综合久久图片| 国产精品一区二区免费欧美| 精华霜和精华液先用哪个| 三级毛片av免费| 精品国产亚洲在线| 久久久久久久精品吃奶| 亚洲成av人片在线播放无| 亚洲,欧美精品.| 法律面前人人平等表现在哪些方面| 亚洲成人中文字幕在线播放| 在线观看日韩欧美| 国产亚洲av高清不卡| 最近最新中文字幕大全免费视频| 看黄色毛片网站| 中文字幕人妻丝袜一区二区| 国产久久久一区二区三区| 精品久久久久久,| 99久久精品国产亚洲精品| 国内少妇人妻偷人精品xxx网站 | 黄色日韩在线| 欧美黑人巨大hd| 成人国产综合亚洲| 性色av乱码一区二区三区2| 人人妻人人看人人澡| 人妻夜夜爽99麻豆av| 九九热线精品视视频播放| 国产欧美日韩精品一区二区| 亚洲av成人av| 亚洲在线观看片| 97超视频在线观看视频| 日韩欧美国产在线观看| 婷婷六月久久综合丁香| 国产精品久久电影中文字幕| 美女午夜性视频免费| 桃色一区二区三区在线观看| 亚洲中文字幕一区二区三区有码在线看 | 男人的好看免费观看在线视频| 在线观看66精品国产| a级毛片在线看网站| 天天一区二区日本电影三级| 中文字幕熟女人妻在线| 日日摸夜夜添夜夜添小说| 欧美黄色片欧美黄色片| 欧洲精品卡2卡3卡4卡5卡区| 美女免费视频网站| 久久精品国产99精品国产亚洲性色| 男人舔奶头视频| 熟女人妻精品中文字幕| 两个人的视频大全免费| 禁无遮挡网站| 女同久久另类99精品国产91| 老鸭窝网址在线观看| 黄色日韩在线| 亚洲五月婷婷丁香| 亚洲成人久久性| 可以在线观看的亚洲视频| 国产精品久久久久久亚洲av鲁大| 国产精品美女特级片免费视频播放器 | 亚洲真实伦在线观看| 亚洲成人久久爱视频| 999久久久精品免费观看国产| 亚洲精品久久国产高清桃花| 亚洲中文字幕日韩| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 亚洲av熟女| 亚洲国产精品合色在线| 少妇裸体淫交视频免费看高清| 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 成人欧美大片| 天堂√8在线中文| 脱女人内裤的视频| 免费高清视频大片| 两个人看的免费小视频| 亚洲激情在线av| 高清在线国产一区| 男女做爰动态图高潮gif福利片| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 亚洲欧美激情综合另类| 久久久久久人人人人人| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 黄片大片在线免费观看| 国产视频一区二区在线看| 在线观看午夜福利视频| 国产精品98久久久久久宅男小说| 两个人视频免费观看高清| 国产精品乱码一区二三区的特点| 日韩欧美免费精品| 国产伦在线观看视频一区| 国产精品av久久久久免费| 香蕉久久夜色| 91av网站免费观看| 国产伦在线观看视频一区| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 舔av片在线| 无限看片的www在线观看| 亚洲精华国产精华精| 久久精品91蜜桃| 美女 人体艺术 gogo| 亚洲国产精品999在线| 黄色丝袜av网址大全| 国产黄a三级三级三级人| 欧美日本视频| cao死你这个sao货| 精品一区二区三区视频在线观看免费| 狂野欧美白嫩少妇大欣赏| 国产精品久久久av美女十八| 黄色丝袜av网址大全|