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

    低雷諾數(shù)條件下湍流度及湍流梯度對(duì)翼型氣動(dòng)特性的影響

    2021-09-07 02:01:20張陽(yáng)周洲李旭
    關(guān)鍵詞:迎角雷諾數(shù)升力

    張陽(yáng), 周洲, 李旭

    (西北工業(yè)大學(xué) 航空學(xué)院, 陜西 西安 710072)

    近年來(lái)對(duì)垂直起降(vertical take-off and landing,VTOL)固定翼無(wú)人機(jī)的研究備受關(guān)注,極具代表性的有XV-24“雷擊”無(wú)人機(jī)、NASA的GL-10驗(yàn)證機(jī)、XC-142A驗(yàn)證機(jī)等多種型號(hào)[1]。該類無(wú)人機(jī)普遍采用分布式動(dòng)力來(lái)提供起降的升力,由于分布式動(dòng)力系統(tǒng)對(duì)空氣的驅(qū)動(dòng),其后方的流場(chǎng)變得極其復(fù)雜,不僅被加速,并且伴有湍流度變大、氣流紊亂等特點(diǎn)。機(jī)翼處于剛啟動(dòng)的分布式動(dòng)力系統(tǒng)的噴流中,特征雷諾數(shù)極低,流場(chǎng)中的流動(dòng)情況復(fù)雜多變,可能存在層流邊界層、層流分離、轉(zhuǎn)捩、湍流邊界層和湍流分離等現(xiàn)象[2]。另外,在分布式動(dòng)力系統(tǒng)噴流的作用下,湍流度的變化也較為劇烈,其大小也同樣影響著機(jī)翼的氣動(dòng)特性[3]。一般認(rèn)為湍流強(qiáng)度小于或等于1%的湍流為低強(qiáng)度湍流,大于10%為高強(qiáng)度湍流,而機(jī)翼在實(shí)際應(yīng)用中多處于大于5%湍流強(qiáng)度的氣流中。剛啟動(dòng)的小型風(fēng)機(jī),其后方流動(dòng)在受到大氣湍流和風(fēng)機(jī)尾流湍流影響時(shí),湍流度通??梢源笥?0%[4];航空發(fā)動(dòng)機(jī)由于燃燒室中大量摻混射流和火焰筒內(nèi)復(fù)雜冷卻結(jié)構(gòu)的存在,使得其出口的湍流度達(dá)15%~20%[5]。

    VTOL無(wú)人機(jī)分布式動(dòng)力系統(tǒng)從啟動(dòng)到正常運(yùn)轉(zhuǎn)的過(guò)程中,噴流的多變性影響著處于動(dòng)力系統(tǒng)后機(jī)翼的特征雷諾數(shù)以及流場(chǎng)的湍流度變化。機(jī)翼的特征雷諾數(shù)可以由幾千變化到幾十萬(wàn),經(jīng)歷的湍流度可以大于20%,同時(shí)機(jī)翼在其弦長(zhǎng)范圍內(nèi)亦會(huì)經(jīng)歷湍流度的急劇變化。因此,面對(duì)VTOL無(wú)人機(jī)分布式動(dòng)力噴流后方的復(fù)雜流場(chǎng),研究低雷諾數(shù)(Re<5×105)條件下,受湍流度及湍流度梯度影響的翼型氣動(dòng)特性顯得尤為重要。

    國(guó)內(nèi)外在低雷諾數(shù)條件下湍流度對(duì)翼型影響的研究主要針對(duì)低來(lái)流速度下的低速翼型。研究人員[6-8]通過(guò)實(shí)驗(yàn)現(xiàn)象分析了低雷諾數(shù)下湍流度的增加對(duì)機(jī)翼空氣動(dòng)力學(xué)性能提升的機(jī)理以及湍流度與雷諾數(shù)對(duì)翼型氣動(dòng)性能影響的相似性,Lasse[9]、王庶等[10]通過(guò)數(shù)值仿真對(duì)低雷諾數(shù)下機(jī)翼表面由于不同湍流度而產(chǎn)生不同程度的流動(dòng)分離等流場(chǎng)細(xì)節(jié)進(jìn)行了進(jìn)一步的研究。已有研究主要是對(duì)特定雷諾數(shù)下的來(lái)流湍流度影響進(jìn)行分析,本文針對(duì)這一背景,剝離分布式動(dòng)力噴流的三維效應(yīng)、旋流等因素的影響,重點(diǎn)對(duì)低雷諾數(shù)、湍流度、湍流梯度等因素對(duì)翼型的氣動(dòng)影響開(kāi)展更深入的研究。

    本文以VTOL飛行器起飛時(shí)分布式動(dòng)力系統(tǒng)變工況及復(fù)雜流場(chǎng)變化為研究背景,基于計(jì)算流體力學(xué)(computational fluid dynamics,CFD)方法使用商業(yè)軟件Fluent對(duì)處于不同湍流度、雷諾數(shù)以及湍流度梯度條件下的翼型氣動(dòng)性能進(jìn)行了研究分析。通過(guò)NACA0012翼型氣動(dòng)力與實(shí)驗(yàn)值的對(duì)比,驗(yàn)證了本文采用的數(shù)值計(jì)算方法。主要研究了翼型在低雷諾數(shù)條件下湍流度及雷諾數(shù)對(duì)翼型氣動(dòng)特性及流場(chǎng)特征的影響,不同湍流度梯度變化規(guī)律影響下流場(chǎng)特征的差異;重點(diǎn)分析了翼型受不同湍流度影響時(shí)氣動(dòng)力發(fā)生變化的原因,湍流度/雷諾數(shù)/湍流梯度變化對(duì)翼型繞流轉(zhuǎn)捩產(chǎn)生影響的機(jī)理,分離泡的產(chǎn)生與演化等。

    1 數(shù)值研究方法

    1.1 湍流度估算方法

    湍流強(qiáng)度Tu(turbulence intensity)是衡量流體微團(tuán)無(wú)規(guī)則運(yùn)動(dòng)程度的物理量。不同于層流流動(dòng)中流體微團(tuán)保持互相平行的層狀運(yùn)動(dòng),湍流使得流體微團(tuán)處于無(wú)規(guī)則的旋渦式運(yùn)動(dòng)中,而不同尺寸旋渦團(tuán)間相互摻混,引起強(qiáng)烈的能量和熱量的交換和傳輸,并引起機(jī)械能量的迅速耗散。湍流強(qiáng)度Tu定義為脈動(dòng)速度u′的均方根與平均速度U的比值

    (1)

    式中:

    湍流能量主要集中在大尺度渦結(jié)構(gòu)中,湍流尺度l(turbulence length scale)就是與攜帶湍流能量的大尺度渦結(jié)構(gòu)有關(guān)的物理量。湍流的特征長(zhǎng)度取決于湍流發(fā)展具有決定性影響的幾何長(zhǎng)度。自由流湍流效應(yīng)隨長(zhǎng)度尺度的增大而減小,至少部分歸因于壁面法向分量速度漲落的減小;自由剪切層則不是這樣[3,6]。入口處的流動(dòng)為受到壁面限制并且?guī)в型牧鬟吔鐚拥牧鲃?dòng),可以利用邊界層厚度δ99%,通過(guò)公式l=0.4δ99%來(lái)計(jì)算湍流長(zhǎng)度尺度。對(duì)于入口處收格柵影響的流動(dòng),湍流尺度與格柵孔徑d近似相等,即l≈d。

    湍流黏性比μt/μ(turbulence and viscosity ratio)表示湍流黏性與層流黏性的比值。湍流黏性的大小直接與湍流雷諾數(shù)成比例。在高雷諾數(shù)邊界層,剪切層和充分發(fā)展的管流中湍流黏性比較大,大約為100~1 000量級(jí)。但在大多數(shù)外流的自由來(lái)流邊界層中湍流黏性比相當(dāng)小,通常為1<μt/μ<10。

    在相對(duì)風(fēng)洞靜止的參考系中觀察,風(fēng)洞格柵后的流動(dòng)是定常湍流流動(dòng),其強(qiáng)度從實(shí)驗(yàn)段入口處的Tuinlet開(kāi)始沿流向不斷衰減,根據(jù)文獻(xiàn)[11]的推導(dǎo),可以得到從實(shí)驗(yàn)段入口處開(kāi)始風(fēng)洞湍流的湍流度隨流向距離的增大而衰減的關(guān)系式Tu=F(Rex)。

    (2)

    式中:β,β*為湍流模型中的經(jīng)驗(yàn)常數(shù),其值分別取β=0.09;β*=0.082 8。

    (3)

    (4)

    與SST流動(dòng)模型耦合方程為

    (5)

    其他參數(shù)具體定義可參考文獻(xiàn)[12-13]。

    2 數(shù)值模擬方法驗(yàn)證

    參考文獻(xiàn)[10],對(duì)NACA0012翼型在超低雷諾數(shù)(Re=5 300)條件下進(jìn)行了數(shù)值模擬。翼型迎角α變化范圍為0°~25°,弦長(zhǎng)c=0.002 5 m,空氣入口邊界來(lái)流速度為V=30.96 m/s。計(jì)算域半徑為弦長(zhǎng)的50倍,采用C型結(jié)構(gòu)網(wǎng)格進(jìn)行劃分,網(wǎng)格及邊界條件的定義見(jiàn)圖1。網(wǎng)格劃分時(shí)給定第一層網(wǎng)格高度d1=6.53×10-4c,y+=0.25,網(wǎng)格單元量為125 000。

    在Tu∞=6%的條件下,計(jì)算得到圖2所示的數(shù)值計(jì)算值與實(shí)驗(yàn)值的對(duì)比。從圖中可以看出,升力系數(shù)在10°迎角以內(nèi)擬合良好,在12°迎角處誤差稍大,阻力特性曲線的變化趨勢(shì)與實(shí)驗(yàn)值保持一致。

    圖1 網(wǎng)格劃分及邊界條件 圖2 數(shù)值計(jì)算與實(shí)驗(yàn)結(jié)果對(duì)比

    從以上的對(duì)比結(jié)果可以看出,本文采用的計(jì)算方法適用于低雷諾數(shù)條件下不同湍流度變化對(duì)翼型影響的復(fù)雜流場(chǎng)計(jì)算,并且具有較高精度。下文繼續(xù)選用較為常用的對(duì)稱翼型NACA0012作為基礎(chǔ)翼型,來(lái)研究處于不同湍流度/雷諾數(shù)/湍流度梯度流場(chǎng)的翼型的氣動(dòng)特性。

    3 不同湍流度/雷諾數(shù)對(duì)翼型的影響

    選取Re1=26 500,Re2=53 000的工況對(duì)翼型繞流進(jìn)行數(shù)值模擬,計(jì)算過(guò)程中保證翼型前緣處的湍流度分別為T(mén)u=6%,Tu=3.2%及Tu=0.6%。圖3給出數(shù)值計(jì)算得到的升阻力特性對(duì)比(其中,Cd,p為翼型壓差阻力因數(shù),Cd,f為翼型摩擦阻力因數(shù))。

    通過(guò)升力特性對(duì)比可以看到:在Tu=6%時(shí),2種雷諾數(shù)下的翼型在中小迎角范圍內(nèi)(0°~8°)升力系數(shù)隨著迎角的增加呈線性增長(zhǎng),失速迎角均為12°。在較小湍流度Tu=0.6%時(shí),翼型升力系數(shù)的非線性較為明顯,同時(shí)失速迎角變小,在Re1=26 500時(shí)升力特性變差,失速迎角僅為8°。

    類似的,湍流度的變化對(duì)翼型升力系數(shù)也有較大影響。在Re1=26 500時(shí),翼型相同迎角的升力系數(shù)隨著湍流度的增大而增大;而在Re2=53 000時(shí),翼型小迎角范圍內(nèi)(2°~5°)的升力系數(shù)在較小湍流度時(shí)反而更大,這很可能是湍流度較小時(shí)翼型上表面產(chǎn)生層流分離泡造成的結(jié)果。

    圖3 不同雷諾數(shù)下的氣動(dòng)特性對(duì)比

    另外,翼型的阻力特性也隨著湍流度/雷諾數(shù)的降低而變差。翼型在10°~12°迎角附近出現(xiàn)了阻力激增現(xiàn)象,圖3c)反映了在該迎角范圍內(nèi)翼型的壓差阻力突然變大,這是由于隨著迎角的增加,翼型上表面的層流分離泡向后緣發(fā)展成為后緣分離泡,逐步演化為大迎角流動(dòng)分離,翼型失速的結(jié)果。

    各個(gè)湍流度/雷諾數(shù)條件下的翼型升阻特性表現(xiàn)出較大的差異,一方面說(shuō)明在低雷諾數(shù)條件下湍流度/雷諾數(shù)的變化對(duì)翼型氣動(dòng)特性影響至關(guān)重要,另一方面也反映了湍流度/雷諾數(shù)的增大對(duì)翼型的氣動(dòng)影響具有相似的作用:隨著湍流度/雷諾數(shù)的增大,翼型升力特性及失速特性均有改善。下文通過(guò)翼型繞流流動(dòng)特性及流場(chǎng)形態(tài)進(jìn)一步研究二者對(duì)翼型氣動(dòng)特性影響的機(jī)理。

    3.1 湍流度對(duì)翼型的影響

    在雷諾數(shù)Re1=26 500條件下,選取翼型繞流具有明顯差異的5°,8°迎角狀態(tài)進(jìn)行分析,圖4給出了處于不同湍流度的翼型表面壓力因數(shù)分布對(duì)比,圖5為受不同湍流度影響的翼型表面(弦向)摩擦阻力因數(shù)的變化。圖中,Cp為壓力因數(shù),Cf,X為弦向摩擦阻力因數(shù)。

    圖4 受不同湍流度影響的翼型壓力因數(shù)分布 圖5 受不同湍流度影響的翼型摩阻因數(shù)

    對(duì)比不同湍流度下的壓力分布、摩阻系數(shù)曲線,結(jié)合對(duì)應(yīng)狀態(tài)的流場(chǎng)細(xì)節(jié)(見(jiàn)圖6至7),可以看到:

    當(dāng)翼型迎角為5°時(shí),處于Tu=6%的翼型表面壓力分布曲線光滑,流線附著良好,說(shuō)明翼型表面流動(dòng)穩(wěn)定(如圖4a)及圖6a)所示)。處于Tu=3.2%的翼型壓力分布在20%c~70%c附近出現(xiàn)輕微鼓包,其Cf,X從21%c開(kāi)始變?yōu)樨?fù)值,直至67%c處回到零值,說(shuō)明在該范圍內(nèi)出現(xiàn)了流動(dòng)分離、轉(zhuǎn)捩、再附的過(guò)程,形成了典型的層流分離泡(如圖5a)及圖6b)所示)。而處于Tu=0.6%的翼型壓力分布在其上表面約30%c~60%c附近形成壓力平臺(tái),上表面的Cf,X從18%c開(kāi)始變?yōu)樨?fù)值,此時(shí)流動(dòng)發(fā)生分離,上表面流線包裹細(xì)長(zhǎng)的類似分離泡的結(jié)構(gòu)(如圖6c)所示),但是流動(dòng)并未重新附著于翼面上。

    圖6 α=5°時(shí)不同湍流度的流場(chǎng)結(jié)構(gòu)

    當(dāng)翼型迎角為8°時(shí),處于Tu=6%的翼型在其上表面5%c~25%c形成短分離泡(如圖5b)及圖7a)所示),轉(zhuǎn)捩位置發(fā)生在上表面13%c附近(Cf,X谷值的位置反映了流動(dòng)轉(zhuǎn)捩的開(kāi)始)。處于Tu=3.2%的翼型前緣分離泡長(zhǎng)度增長(zhǎng),為4%c~33%c,轉(zhuǎn)捩位置23%c有所推遲。處于Tu=0.6%的翼型其上表面的流動(dòng)分離區(qū)域更加明顯(如圖7b)所示),流動(dòng)在4%c附近發(fā)生分離,并一直延伸至后緣,盡管未能形成再附,但是通過(guò)圖5b)可以看到其轉(zhuǎn)捩位置更加推遲,為33%c??偟膩?lái)說(shuō),隨著湍流度減小,翼型上表面流動(dòng)分離點(diǎn)前移,轉(zhuǎn)捩位置推遲,分離泡長(zhǎng)度增長(zhǎng),流動(dòng)分離區(qū)域更大,翼型失速迎角變小。

    圖7 α=8°時(shí)不同湍流度的流場(chǎng)結(jié)構(gòu)

    通過(guò)翼型在Tu=0.6%時(shí)的流場(chǎng)變化可以看出在低湍流度下翼型上表面流動(dòng)分離的演化過(guò)程:隨著迎角的增加,翼型后緣附近邊界層的逆壓梯度不斷增大,并引發(fā)邊界層的分離,隨后分離點(diǎn)逐漸向前移動(dòng),直到大迎角時(shí)翼型上表面完全分離,最后失速。而當(dāng)翼型處于Tu=6%時(shí),上述流動(dòng)分離受到很大程度的抑制,使得翼型大迎角特性得到了改善。這說(shuō)明,邊界層外較大的湍流度使得流動(dòng)穩(wěn)定性增強(qiáng),為繞翼型的上游層流注入能量,增強(qiáng)了抵抗附面層內(nèi)逆壓梯度的能力,延緩了層流分離的發(fā)生。另外湍流度的增大也使得流動(dòng)轉(zhuǎn)捩提前發(fā)生,為流動(dòng)在壁面的再附提供了可能。

    3.2 雷諾數(shù)對(duì)翼型的影響

    雷諾數(shù)對(duì)翼型的影響在文獻(xiàn)中已經(jīng)有很多討論,本節(jié)主要通過(guò)雷諾數(shù)由26 500增大53 000時(shí)翼型氣動(dòng)特性的變化來(lái)反映雷諾數(shù)對(duì)以下幾方面的影響:失速迎角的提高、轉(zhuǎn)捩的提前及分離泡的產(chǎn)生。圖8為Re2=53 000,翼型8°迎角工況下的翼型壓力分布與表面摩擦阻力變化。

    圖8 α=8°受不同湍流度影響的翼型壓力分布及摩阻系數(shù)

    由8圖可以看出,翼型在不同湍流度的影響下在前緣附近的流動(dòng)都發(fā)生了轉(zhuǎn)捩并產(chǎn)生短分離泡。對(duì)比Re1=26 500的翼型特性(見(jiàn)圖4b)及圖5b)),翼型在Re2=53 000,Tu=0.6%時(shí),隨著雷諾數(shù)的增大,由翼型后緣開(kāi)始的上表面大迎角分離轉(zhuǎn)變?yōu)楹缶壐浇男》秶蛛x,流動(dòng)分離的程度受到了抑制。這說(shuō)明雷諾數(shù)的增大,翼型表面剪切效應(yīng)增強(qiáng),動(dòng)能更充沛,抗逆壓梯度能力增強(qiáng),使得失速迎角提高。

    同時(shí),結(jié)合圖5b)及圖8可以看到,雷諾數(shù)增大使得翼型上表面的流動(dòng)分離提前發(fā)生。在Tu=3.2%時(shí),翼型上表面分離點(diǎn)由4%c前移到3%c,并且流動(dòng)在15%c重新附著在壁面上,層流分離泡的長(zhǎng)度由4%c~33%c縮短為3%c~15%c。通過(guò)Cf,X谷值的位置可以看到,隨著雷諾數(shù)的增大,轉(zhuǎn)捩位置由23%c前移到11%c。這個(gè)結(jié)果與文獻(xiàn)[14-15]描述一致:雷諾數(shù)的增大會(huì)導(dǎo)致流動(dòng)的提前轉(zhuǎn)捩及再附,也因此會(huì)縮短分離泡的長(zhǎng)度。分離泡的產(chǎn)生改變了翼型的有效形狀,增大了翼型上表面曲率,使得翼型在小湍流度下的升力系數(shù)反而高于高湍流度,導(dǎo)致了如圖3a)所示的小迎角下升力系數(shù)的非線性現(xiàn)象。然而當(dāng)分離泡尺度較小或產(chǎn)生于翼型后緣附近時(shí),對(duì)翼型的升力特性影響較小。

    對(duì)不同湍流度/雷諾數(shù)條件下翼型氣動(dòng)特性的對(duì)比表明,低雷諾數(shù)條件下翼型周圍的流場(chǎng)十分不穩(wěn)定,翼型上表面的流動(dòng)現(xiàn)象豐富,湍流度及雷諾數(shù)對(duì)流場(chǎng)特征影響相似,主要表現(xiàn)在翼型升力的提高,失速迎角的增大,但是引發(fā)轉(zhuǎn)捩提前的機(jī)理較為不同。湍流度增加導(dǎo)致的轉(zhuǎn)捩提前是因?yàn)榱鲃?dòng)較高的速度脈動(dòng)注入了邊界層。而雷諾數(shù)增加引起的轉(zhuǎn)捩提前則是因?yàn)榧羟袑拥暮穸仍诹鲃?dòng)分離點(diǎn)附近減小,從而改變了邊界層的穩(wěn)定性。Simoni等[16]的研究表明雷諾數(shù)的增加顯著影響了分離點(diǎn)處邊界層的結(jié)構(gòu),導(dǎo)致渦脫落現(xiàn)象較高的成長(zhǎng)率和不同的動(dòng)態(tài)屬性。

    4 湍流度梯度變化規(guī)律對(duì)翼型的影響

    由于流場(chǎng)的湍流度變化會(huì)隨著氣流的擾動(dòng)產(chǎn)生較大波動(dòng),機(jī)翼在其弦長(zhǎng)范圍內(nèi)亦會(huì)經(jīng)歷不同的湍流度梯度而產(chǎn)生不同的氣動(dòng)特性。本文將湍流度梯度定義為

    (6)

    式中:ΔTu為翼型前緣到后緣湍流度的變化量;c為翼型弦長(zhǎng),以此值的大小來(lái)表示湍流度梯度對(duì)翼型在其弦長(zhǎng)范圍內(nèi)的氣動(dòng)影響。下文將研究湍流度梯度變化對(duì)翼型氣動(dòng)性能的影響。

    4.1 雷諾數(shù)對(duì)湍流度梯度的影響

    計(jì)算雷諾數(shù)分別選取Re0=5 300,Re1=26 500,參考自由來(lái)流、風(fēng)機(jī)尾流、航空發(fā)動(dòng)機(jī)尾流等不同湍流度狀態(tài)的流動(dòng)情況,計(jì)算邊界流動(dòng)入口處的湍流度分別設(shè)置為T(mén)u,inlet=9%,15%,25%。保證翼型(弦長(zhǎng)0.1 m)前緣處湍流強(qiáng)度均為T(mén)u=6%,而翼型后緣由于湍流度梯度變化規(guī)律的不同而產(chǎn)生差異,流場(chǎng)湍流度梯度變化分別為?Tu/?x=0.040,0.035,0.022。

    計(jì)算得到不同雷諾數(shù)下3種湍流度梯度影響的翼型升力曲線,如圖9所示??梢钥吹?在Re0=5 300時(shí),翼型升力系數(shù)的差異主要體現(xiàn)在迎角α=5°~15°范圍內(nèi),翼型升力系數(shù)隨著湍流度梯度的增大而增大,?Tu/?x=0.040狀態(tài)下的翼型產(chǎn)生了較大的升力。而在Re1=26 500時(shí),湍流度梯度對(duì)翼型升力特性影響有限,升力系數(shù)差別較小,升力曲線趨近一致。

    圖9 湍流度梯度對(duì)翼型升力特性的影響

    下面對(duì)翼型迎角為10°時(shí),不同湍流度梯度影響下的翼型表面摩擦阻力系數(shù)(見(jiàn)圖10)進(jìn)行分析。

    圖10 湍流度梯度對(duì)翼型在α=10°時(shí)摩阻系數(shù)的影響

    由可以看到,在Re0=5 300時(shí),在較大湍流度梯度(?Tu/?x=0.040)的影響下,翼型上表面的流動(dòng)分離最早發(fā)生,但是由于雷諾數(shù)過(guò)低,流動(dòng)分離后無(wú)法再附著到壁面上,而流動(dòng)轉(zhuǎn)捩依然存在,翼型表面摩擦阻力系數(shù)反映了在較大湍流度梯度時(shí)流動(dòng)轉(zhuǎn)捩也更早。在Re1=26 500時(shí),翼型表面摩擦阻力因數(shù)曲線在不同湍流度梯度的影響下幾乎重合,但是依然可以發(fā)現(xiàn),翼型上表面流動(dòng)分離后的再附位置略有差異:較大的湍流度梯度對(duì)應(yīng)著較早的流動(dòng)再附。

    圖11為在Re0=5 300,翼型迎角10°時(shí)不同湍流度梯度對(duì)應(yīng)的流場(chǎng)結(jié)構(gòu),可以看到:翼型在不同的湍流度梯度下,翼型流線形態(tài)的差異主要體現(xiàn)在翼型上表面的后半部分:隨著湍流度梯度的增大,流動(dòng)分離的現(xiàn)象受到一定程度抑制,分離渦中心位置略微前移,分離區(qū)域后半部分逐漸被壓縮而更靠近壁面。另外,湍動(dòng)能集中的位置隨著湍流度梯度的變大而發(fā)生了很小的前移和下移,湍動(dòng)能能量集中位置的變化說(shuō)明流動(dòng)轉(zhuǎn)捩提前發(fā)生,分離渦中心位置略微前移,這使得在轉(zhuǎn)捩點(diǎn)之后的湍流剪切層被外部自由流體提前注入能量,流動(dòng)的分離得到了更早的抑制作用,從而使得之后的流線向壁面靠近。這樣的分離渦改變了翼型的有效形狀,流線經(jīng)歷了更大的流動(dòng)曲率的變化,從而使得處于較高湍流梯度的翼型產(chǎn)生了更大的升力。

    圖11 在Re0=5,300,α=10°時(shí)不同湍流度梯度的流場(chǎng)結(jié)構(gòu)

    4.2 湍流度對(duì)湍流度梯度的影響

    以雷諾數(shù)Re0=5 300,計(jì)算邊界流動(dòng)入口湍流度為T(mén)u,inlet=25%,15%的工況,研究幾種不同前緣湍流度Tu=12.0%,9.4%,6.0%,2.6%對(duì)湍流度梯度的影響。翼型10°迎角時(shí)的摩阻系數(shù)以及幾種不同湍流度對(duì)應(yīng)的湍流度梯度如圖12所示。

    圖12 不同湍流度及湍流度梯度的摩阻系數(shù)

    可以看到,隨前緣位置湍流度增大,不同入口湍流度Tu,inlet(15%,25%)下對(duì)應(yīng)的湍流度梯度也增大(?Tu/?x由0.007/0.011增大到0.171/0.203),同時(shí)其對(duì)翼型上表面流動(dòng)的影響越來(lái)越大。在前緣湍流度從2.6%增大到12.0%的過(guò)程中,流動(dòng)的轉(zhuǎn)捩位置(圖中以圓○標(biāo)明)由64%c提前到38%c。而湍流度梯度的增大使得在同一湍流度狀態(tài)下的流動(dòng)分離位置(圖中以菱形◇標(biāo)明)有所提前,同時(shí)分離點(diǎn)距離相差越來(lái)越大。

    選取湍流度梯度對(duì)翼型10°迎角時(shí)影響較大的狀態(tài)(Tu=12.0%,α=10°)進(jìn)行進(jìn)一步的流場(chǎng)分析,如圖13所示。

    圖13 不同湍流度梯度的流場(chǎng)結(jié)構(gòu)

    可以看到,湍流度梯度的增大使得流場(chǎng)中的湍動(dòng)能衰減更加迅速:在?Tu/?x=0.171時(shí),流場(chǎng)中湍動(dòng)能為20的量值可以延伸至接近尾緣的位置,而在?Tu/?x=0.203時(shí),在接近弦長(zhǎng)40%的位置湍動(dòng)能便衰減到20以下。翼型上表面的分離泡演化與圖11(Tu=6.0%)一致,同樣是隨著湍流度梯度的增大,分離泡后半部分受到了一定程度的壓縮,而不同的是分離泡整體厚度變得更小,后半部分被壓縮得更加明顯。

    在其他較小前緣湍流度的工況下,盡管使湍流度梯度增大可以使得流動(dòng)分離、轉(zhuǎn)捩、再附等現(xiàn)象提前,但是湍流度梯度對(duì)翼型的氣動(dòng)影響及流動(dòng)特征的改變變得更加有限,說(shuō)明湍流度的大小制約著湍流度梯度對(duì)流場(chǎng)的影響能力。

    5 結(jié) 論

    本文以VTOL飛行器起降時(shí)分布式動(dòng)力噴流產(chǎn)生的復(fù)雜流動(dòng)為背景,研究了不同湍流度/雷諾數(shù)、不同湍流度梯度對(duì)翼型氣動(dòng)特性的影響,得到以下結(jié)論:

    1) 湍流度/雷諾數(shù)的增大會(huì)導(dǎo)致流動(dòng)轉(zhuǎn)捩提前,翼型邊界層抗逆壓梯度能力增強(qiáng),進(jìn)而使得氣動(dòng)特性得到改善。但是湍流度/雷諾數(shù)導(dǎo)致流動(dòng)轉(zhuǎn)捩提前的機(jī)理有所不同:湍流度的增大會(huì)使得具有較高速度脈動(dòng)的流體微團(tuán)進(jìn)入翼型邊界層;而雷諾數(shù)的增大會(huì)減小剪切層厚度進(jìn)而改變流動(dòng)分離點(diǎn)后方邊界層的穩(wěn)定性。

    2) 湍流度/雷諾數(shù)的增大會(huì)使得翼型上表面層流分離泡提前出現(xiàn),但是分離泡的長(zhǎng)度會(huì)隨之縮短。分離泡的產(chǎn)生會(huì)改變翼型的有效形狀,增大翼型上表面曲率,從而導(dǎo)致小迎角下升力系數(shù)的非性現(xiàn)象。

    3) 湍流度梯度的增大在一定程度上會(huì)使得翼型表面流動(dòng)分離位置、轉(zhuǎn)捩位置及流動(dòng)再附位置提前,但是這種影響受到雷諾數(shù)和湍流度數(shù)值的制約,會(huì)隨著雷諾數(shù)的增大而衰弱,隨著湍流度的增大而加強(qiáng)。

    猜你喜歡
    迎角雷諾數(shù)升力
    高速列車車頂–升力翼組合體氣動(dòng)特性
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    無(wú)人機(jī)升力測(cè)試裝置設(shè)計(jì)及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    升力式再入飛行器體襟翼姿態(tài)控制方法
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計(jì)研究
    民機(jī)高速風(fēng)洞試驗(yàn)的阻力雷諾數(shù)效應(yīng)修正
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    久久99一区二区三区| 国产永久视频网站| 国产精品国产三级国产av玫瑰| 国产综合精华液| 亚洲欧美中文字幕日韩二区| 国产日韩欧美亚洲二区| 亚洲经典国产精华液单| 国产 一区精品| 国产亚洲最大av| 久久久久网色| 欧美成人午夜免费资源| 国产毛片在线视频| 日产精品乱码卡一卡2卡三| 亚洲av福利一区| 六月丁香七月| 亚洲精品国产av蜜桃| 日韩成人伦理影院| 2022亚洲国产成人精品| 丰满乱子伦码专区| 国产精品99久久99久久久不卡 | 秋霞伦理黄片| 国产男女超爽视频在线观看| 少妇人妻 视频| 偷拍熟女少妇极品色| 亚洲欧美日韩另类电影网站| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 美女主播在线视频| 偷拍熟女少妇极品色| 我要看黄色一级片免费的| 伊人久久国产一区二区| 国产白丝娇喘喷水9色精品| 一区在线观看完整版| 青青草视频在线视频观看| 99热全是精品| 99热这里只有是精品50| 国产精品久久久久久久电影| 日韩大片免费观看网站| 国语对白做爰xxxⅹ性视频网站| 久久久久久久大尺度免费视频| a级一级毛片免费在线观看| 两个人的视频大全免费| 国国产精品蜜臀av免费| 亚洲人与动物交配视频| 美女中出高潮动态图| 亚洲av二区三区四区| 黄片无遮挡物在线观看| 午夜免费鲁丝| 免费看av在线观看网站| 中国美白少妇内射xxxbb| 国产精品久久久久久久电影| 777米奇影视久久| 2021少妇久久久久久久久久久| 精品亚洲成国产av| 欧美性感艳星| 一级黄片播放器| 免费看不卡的av| 男人舔奶头视频| 免费av中文字幕在线| 国产 精品1| 国产淫语在线视频| 男的添女的下面高潮视频| 免费久久久久久久精品成人欧美视频 | 午夜视频国产福利| 精品人妻偷拍中文字幕| 久久午夜综合久久蜜桃| 婷婷色av中文字幕| 亚洲国产日韩一区二区| 久久精品国产自在天天线| 久久精品久久精品一区二区三区| 午夜精品国产一区二区电影| 免费av不卡在线播放| 亚洲伊人久久精品综合| 国产在视频线精品| 国产色婷婷99| 美女大奶头黄色视频| 国产精品欧美亚洲77777| 欧美精品亚洲一区二区| 亚洲精品一二三| 蜜臀久久99精品久久宅男| 在线观看国产h片| 99热这里只有是精品在线观看| 日韩在线高清观看一区二区三区| 日日啪夜夜撸| 人人妻人人添人人爽欧美一区卜| 中文字幕亚洲精品专区| 欧美日韩综合久久久久久| 国产伦精品一区二区三区视频9| 精品国产国语对白av| 高清在线视频一区二区三区| 日本-黄色视频高清免费观看| 亚洲自偷自拍三级| 狂野欧美白嫩少妇大欣赏| 国产精品人妻久久久久久| 免费高清在线观看视频在线观看| 一区二区三区四区激情视频| 久久久久久久久久久丰满| 大陆偷拍与自拍| 各种免费的搞黄视频| 亚洲四区av| 久久国产精品男人的天堂亚洲 | 99视频精品全部免费 在线| 一本色道久久久久久精品综合| 久久国内精品自在自线图片| 两个人免费观看高清视频 | 乱人伦中国视频| 在线天堂最新版资源| xxx大片免费视频| 国产精品一区www在线观看| 国产精品人妻久久久久久| 91aial.com中文字幕在线观看| 免费观看性生交大片5| 美女视频免费永久观看网站| av福利片在线观看| 日韩中文字幕视频在线看片| 国产中年淑女户外野战色| 国产精品99久久久久久久久| 日本wwww免费看| 春色校园在线视频观看| 国产在线男女| 五月伊人婷婷丁香| 一本大道久久a久久精品| 久久婷婷青草| 国产午夜精品久久久久久一区二区三区| 久久99热这里只频精品6学生| 中文字幕免费在线视频6| .国产精品久久| 亚洲精品色激情综合| 欧美精品亚洲一区二区| 国产在线视频一区二区| 我的女老师完整版在线观看| 欧美精品亚洲一区二区| 亚洲,欧美,日韩| 亚洲怡红院男人天堂| 大香蕉97超碰在线| 国产精品不卡视频一区二区| 免费播放大片免费观看视频在线观看| 中文字幕精品免费在线观看视频 | 观看免费一级毛片| 男女啪啪激烈高潮av片| 一边亲一边摸免费视频| 日韩免费高清中文字幕av| 如日韩欧美国产精品一区二区三区 | 男女国产视频网站| 午夜免费鲁丝| 精品国产国语对白av| 国产精品国产三级国产专区5o| 中文在线观看免费www的网站| 国产伦精品一区二区三区四那| av福利片在线观看| 亚洲真实伦在线观看| 一本—道久久a久久精品蜜桃钙片| 国产精品女同一区二区软件| 老司机亚洲免费影院| 丝瓜视频免费看黄片| 伦精品一区二区三区| 午夜日本视频在线| 一级a做视频免费观看| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区三区四区免费观看| 最新的欧美精品一区二区| 久久久久精品久久久久真实原创| 久久青草综合色| 国产精品伦人一区二区| 欧美97在线视频| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 国产成人a∨麻豆精品| 秋霞在线观看毛片| 免费看不卡的av| 欧美xxⅹ黑人| 国产精品一区二区性色av| 精品久久久久久久久亚洲| 夫妻午夜视频| 18禁在线无遮挡免费观看视频| 国产精品国产三级国产专区5o| 精品亚洲乱码少妇综合久久| 色哟哟·www| 一级二级三级毛片免费看| 国产精品一二三区在线看| 高清在线视频一区二区三区| 国产探花极品一区二区| 男女啪啪激烈高潮av片| 亚洲伊人久久精品综合| 欧美国产精品一级二级三级 | 免费黄频网站在线观看国产| 免费大片黄手机在线观看| 亚洲av免费高清在线观看| 亚洲国产欧美日韩在线播放 | 最新中文字幕久久久久| 国产成人免费无遮挡视频| 18禁动态无遮挡网站| 国产成人一区二区在线| 亚洲一区二区三区欧美精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产色片| 男人添女人高潮全过程视频| 亚洲一区二区三区欧美精品| 天堂8中文在线网| 纯流量卡能插随身wifi吗| 精品久久久噜噜| 久久国产精品大桥未久av | 日韩成人av中文字幕在线观看| 又大又黄又爽视频免费| 激情五月婷婷亚洲| 性色avwww在线观看| 中文字幕制服av| 国产亚洲av片在线观看秒播厂| 日韩成人伦理影院| 日韩欧美 国产精品| 免费播放大片免费观看视频在线观看| 嫩草影院入口| 啦啦啦啦在线视频资源| 男人添女人高潮全过程视频| 久久久久久久久久久免费av| 精品久久久久久久久亚洲| 久久ye,这里只有精品| 欧美丝袜亚洲另类| a级一级毛片免费在线观看| 婷婷色麻豆天堂久久| 国产综合精华液| 极品教师在线视频| 欧美日本中文国产一区发布| 3wmmmm亚洲av在线观看| 亚洲情色 制服丝袜| 99热全是精品| 大香蕉久久网| 观看av在线不卡| 欧美精品国产亚洲| 午夜老司机福利剧场| 中文字幕制服av| 国产伦在线观看视频一区| 寂寞人妻少妇视频99o| 搡女人真爽免费视频火全软件| 国产综合精华液| 老司机亚洲免费影院| 秋霞伦理黄片| 亚洲第一区二区三区不卡| 日韩大片免费观看网站| 大片免费播放器 马上看| 国产爽快片一区二区三区| 欧美三级亚洲精品| 亚洲内射少妇av| 丝袜脚勾引网站| 免费人成在线观看视频色| 插逼视频在线观看| 极品人妻少妇av视频| 高清不卡的av网站| 国产 精品1| 老女人水多毛片| 国产黄片美女视频| 免费不卡的大黄色大毛片视频在线观看| 国产免费视频播放在线视频| 在线亚洲精品国产二区图片欧美 | 亚洲国产精品成人久久小说| 中文在线观看免费www的网站| 成人毛片a级毛片在线播放| 爱豆传媒免费全集在线观看| 国产综合精华液| 国产精品久久久久成人av| 亚洲av在线观看美女高潮| 日韩在线高清观看一区二区三区| 交换朋友夫妻互换小说| 国产国拍精品亚洲av在线观看| 男人和女人高潮做爰伦理| 久久精品国产亚洲av天美| 国产探花极品一区二区| 国产真实伦视频高清在线观看| 欧美国产精品一级二级三级 | 亚洲第一区二区三区不卡| 丰满人妻一区二区三区视频av| 亚洲av免费高清在线观看| 日韩视频在线欧美| 日韩欧美 国产精品| 男女国产视频网站| 国产成人aa在线观看| 中文字幕久久专区| 大话2 男鬼变身卡| 午夜久久久在线观看| 在现免费观看毛片| 国产精品国产av在线观看| 国内少妇人妻偷人精品xxx网站| 99热网站在线观看| 久久青草综合色| 国产日韩欧美亚洲二区| 亚洲国产精品一区二区三区在线| 亚洲精品乱久久久久久| 国产在视频线精品| 久久精品国产鲁丝片午夜精品| 十八禁网站网址无遮挡 | 欧美xxxx性猛交bbbb| 久热这里只有精品99| 成人毛片60女人毛片免费| 美女国产视频在线观看| 女性被躁到高潮视频| a级一级毛片免费在线观看| 91久久精品电影网| 伦理电影免费视频| 国产乱来视频区| 成人二区视频| 尾随美女入室| 交换朋友夫妻互换小说| 免费av中文字幕在线| 亚洲四区av| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 成人国产麻豆网| 一级二级三级毛片免费看| 少妇 在线观看| 久久久a久久爽久久v久久| 一级毛片黄色毛片免费观看视频| 亚州av有码| 18禁在线无遮挡免费观看视频| 亚洲精品乱码久久久v下载方式| 中文天堂在线官网| 我要看日韩黄色一级片| 亚洲国产色片| 日韩欧美一区视频在线观看 | 午夜老司机福利剧场| 少妇精品久久久久久久| 亚洲欧美精品自产自拍| 国产亚洲91精品色在线| 精品人妻熟女av久视频| 99国产精品免费福利视频| 亚洲av国产av综合av卡| 久久久久人妻精品一区果冻| 一级黄片播放器| 五月开心婷婷网| 有码 亚洲区| 极品教师在线视频| 国产av精品麻豆| 我的女老师完整版在线观看| 51国产日韩欧美| 高清黄色对白视频在线免费看 | 一级,二级,三级黄色视频| 国内精品宾馆在线| 在线看a的网站| 久久国内精品自在自线图片| 国产av码专区亚洲av| 一区二区三区免费毛片| 九九在线视频观看精品| 午夜福利影视在线免费观看| 欧美最新免费一区二区三区| 最后的刺客免费高清国语| 欧美最新免费一区二区三区| 王馨瑶露胸无遮挡在线观看| 十八禁网站网址无遮挡 | 王馨瑶露胸无遮挡在线观看| 卡戴珊不雅视频在线播放| 久热久热在线精品观看| 国产熟女欧美一区二区| 亚洲国产av新网站| av女优亚洲男人天堂| 欧美+日韩+精品| 99热国产这里只有精品6| 久久6这里有精品| 国产日韩欧美在线精品| 午夜福利在线观看免费完整高清在| 国产成人免费观看mmmm| 高清欧美精品videossex| 高清毛片免费看| 亚洲av日韩在线播放| av天堂久久9| 中文字幕人妻熟人妻熟丝袜美| 在线观看av片永久免费下载| 六月丁香七月| 少妇的逼水好多| 亚洲伊人久久精品综合| 麻豆乱淫一区二区| 日韩成人伦理影院| 久久久久视频综合| 少妇熟女欧美另类| 欧美日韩综合久久久久久| 青春草亚洲视频在线观看| 十八禁高潮呻吟视频 | 能在线免费看毛片的网站| 中文资源天堂在线| 久久久久视频综合| 午夜精品国产一区二区电影| 亚洲美女视频黄频| 在线观看人妻少妇| 亚洲av福利一区| 夜夜爽夜夜爽视频| 亚洲精品国产成人久久av| 97超碰精品成人国产| 亚洲欧洲精品一区二区精品久久久 | 国产精品久久久久久精品古装| 在线观看免费日韩欧美大片 | 国产探花极品一区二区| 综合色丁香网| 大香蕉久久网| 王馨瑶露胸无遮挡在线观看| 精品国产露脸久久av麻豆| 美女cb高潮喷水在线观看| 亚洲电影在线观看av| 国产精品久久久久成人av| 亚洲国产精品一区二区三区在线| 国产黄色视频一区二区在线观看| 国产亚洲5aaaaa淫片| 国产精品久久久久久精品古装| 精华霜和精华液先用哪个| 亚洲精品456在线播放app| 欧美日韩精品成人综合77777| 久久久久久人妻| 色婷婷av一区二区三区视频| 国产色婷婷99| 久久久久人妻精品一区果冻| 日韩,欧美,国产一区二区三区| 免费观看无遮挡的男女| 欧美日韩综合久久久久久| 免费观看的影片在线观看| 国产成人精品无人区| 狂野欧美激情性bbbbbb| 午夜av观看不卡| 在线亚洲精品国产二区图片欧美 | 一本久久精品| 22中文网久久字幕| 麻豆精品久久久久久蜜桃| 午夜av观看不卡| 久久久久久人妻| 99热网站在线观看| 久久婷婷青草| 欧美日韩视频高清一区二区三区二| 国产欧美亚洲国产| 一本久久精品| 久久免费观看电影| 成人二区视频| 99热网站在线观看| 国产欧美日韩精品一区二区| 免费大片18禁| 国产男女超爽视频在线观看| 国产淫片久久久久久久久| 午夜日本视频在线| 亚洲国产成人一精品久久久| 色吧在线观看| 成人午夜精彩视频在线观看| 嫩草影院新地址| 亚洲成色77777| 高清欧美精品videossex| 边亲边吃奶的免费视频| 欧美丝袜亚洲另类| 五月伊人婷婷丁香| 国产色爽女视频免费观看| 一级a做视频免费观看| 日韩成人伦理影院| 六月丁香七月| 中文天堂在线官网| 韩国高清视频一区二区三区| 美女主播在线视频| 日本黄色片子视频| a级毛片在线看网站| 国产亚洲精品久久久com| 亚洲国产精品999| 尾随美女入室| av福利片在线观看| 久久久久人妻精品一区果冻| 久久久国产欧美日韩av| 精品一区二区免费观看| 久久久久久久国产电影| 日韩中字成人| 国产乱来视频区| 男女免费视频国产| 久久97久久精品| 日本与韩国留学比较| av天堂久久9| 国产淫语在线视频| 日韩av免费高清视频| 大又大粗又爽又黄少妇毛片口| 99热网站在线观看| 内地一区二区视频在线| 免费观看a级毛片全部| 极品教师在线视频| 亚洲熟女精品中文字幕| 熟女人妻精品中文字幕| 中文字幕人妻丝袜制服| 国产亚洲精品久久久com| 亚洲欧美一区二区三区黑人 | 国产免费福利视频在线观看| 亚洲国产精品一区三区| 大又大粗又爽又黄少妇毛片口| 亚州av有码| 日韩视频在线欧美| 欧美bdsm另类| a 毛片基地| 欧美xxⅹ黑人| 91久久精品国产一区二区成人| 18+在线观看网站| 国产一区亚洲一区在线观看| 欧美 日韩 精品 国产| 亚洲精品乱码久久久久久按摩| 亚洲va在线va天堂va国产| 久久久久久久久久成人| 有码 亚洲区| 精品亚洲乱码少妇综合久久| 久久精品久久久久久噜噜老黄| 亚洲精品第二区| 天堂中文最新版在线下载| 一区二区三区乱码不卡18| 中文字幕人妻丝袜制服| 国产成人a∨麻豆精品| 国产伦在线观看视频一区| 日本猛色少妇xxxxx猛交久久| 十八禁高潮呻吟视频 | 晚上一个人看的免费电影| 男女啪啪激烈高潮av片| 欧美 亚洲 国产 日韩一| 人妻制服诱惑在线中文字幕| 免费av中文字幕在线| 97在线人人人人妻| 亚洲欧洲国产日韩| 热re99久久精品国产66热6| 一级,二级,三级黄色视频| 最后的刺客免费高清国语| 国产成人freesex在线| √禁漫天堂资源中文www| 中文在线观看免费www的网站| 最黄视频免费看| 免费在线观看成人毛片| 日本-黄色视频高清免费观看| 美女福利国产在线| 人人妻人人爽人人添夜夜欢视频 | 最近手机中文字幕大全| av有码第一页| 色视频在线一区二区三区| 欧美老熟妇乱子伦牲交| 欧美精品国产亚洲| 丝袜喷水一区| 国产一区二区在线观看av| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 少妇 在线观看| 美女大奶头黄色视频| 成年人午夜在线观看视频| 熟女电影av网| 亚洲国产最新在线播放| 免费大片黄手机在线观看| 国产日韩一区二区三区精品不卡 | 久久久久精品性色| 国产精品久久久久久久电影| 91精品伊人久久大香线蕉| 欧美激情国产日韩精品一区| 国产高清国产精品国产三级| 日韩电影二区| 国产精品人妻久久久影院| 亚洲图色成人| 我的女老师完整版在线观看| 女的被弄到高潮叫床怎么办| 一区在线观看完整版| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 极品教师在线视频| 一区在线观看完整版| 韩国高清视频一区二区三区| 欧美精品一区二区免费开放| 免费黄网站久久成人精品| 亚洲av成人精品一二三区| 欧美成人精品欧美一级黄| 久久久久精品久久久久真实原创| 新久久久久国产一级毛片| 99热国产这里只有精品6| 亚洲成人一二三区av| 亚洲真实伦在线观看| 在线观看www视频免费| av在线观看视频网站免费| 水蜜桃什么品种好| 国产精品一区二区性色av| 热99国产精品久久久久久7| 99精国产麻豆久久婷婷| 亚洲婷婷狠狠爱综合网| 99视频精品全部免费 在线| 一二三四中文在线观看免费高清| 国产欧美亚洲国产| 最近中文字幕高清免费大全6| 十八禁网站网址无遮挡 | 日韩,欧美,国产一区二区三区| 成人漫画全彩无遮挡| 日日摸夜夜添夜夜添av毛片| 免费av中文字幕在线| 另类亚洲欧美激情| 大片电影免费在线观看免费| 亚洲av成人精品一区久久| 在线观看三级黄色| 黄色欧美视频在线观看| 啦啦啦在线观看免费高清www| 欧美一级a爱片免费观看看| 大码成人一级视频| 日韩精品免费视频一区二区三区 | 国产黄色视频一区二区在线观看| 高清黄色对白视频在线免费看 | 熟女电影av网| 大香蕉97超碰在线| 国产成人精品久久久久久| 色婷婷久久久亚洲欧美| 亚洲av欧美aⅴ国产| 五月天丁香电影| 国产精品一区二区性色av| 精品一区二区三区视频在线| 久久久久久久久久久久大奶| 久久人人爽人人片av| 日韩av在线免费看完整版不卡| 国产黄色视频一区二区在线观看| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 亚洲精品乱码久久久久久按摩| 日韩一区二区三区影片| 少妇人妻一区二区三区视频| 美女视频免费永久观看网站| 五月伊人婷婷丁香| 久久久国产欧美日韩av| 久久久精品94久久精品| 麻豆成人av视频| 一区二区三区免费毛片| 日韩人妻高清精品专区| 亚洲欧美日韩另类电影网站| 日本色播在线视频|