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

    流向彎曲壁超聲速湍流邊界層研究進(jìn)展

    2020-05-20 02:44:08孫明波王前程
    關(guān)鍵詞:壓力梯度邊界層流向

    孫明波, 王前程, 王 旭, 劉 源

    (國(guó)防科技大學(xué) 高超聲速?zèng)_壓發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 長(zhǎng)沙 410073)

    0 引 言

    自普朗特1904年首先提出邊界層理論以來(lái)以來(lái),邊界層問(wèn)題始終是流體力學(xué)研究的核心問(wèn)題之一。一百多年來(lái),針對(duì)邊界層開展的研究極大地豐富了我們的物理認(rèn)知,顯著地推動(dòng)了航空工業(yè)的發(fā)展。在實(shí)際應(yīng)用中,飛行器表面通常是彎曲的,彎曲壁面的邊界層與平面邊界層有明顯差異,對(duì)航空飛行器部件性能有很大影響。近年來(lái)國(guó)際上高超聲速飛行器成為關(guān)注的焦點(diǎn),其氣動(dòng)布局以及吸氣式推進(jìn)系統(tǒng)設(shè)計(jì)中采用了包括流向形式在內(nèi)的各種曲率,如圖1所示。流向曲率影響下超聲速邊界層流動(dòng)行為成為高超聲速飛行器設(shè)計(jì)關(guān)注的重要內(nèi)容之一。

    (a) 美國(guó)波音公司公布的高超聲速飛機(jī)方案

    (b) 美國(guó)空軍公布的高超聲速打擊武器方案圖1 典型高超聲速飛行器方案:飛行器表面存在包括 流向形式在內(nèi)的各種曲率(圖片來(lái)自互聯(lián)網(wǎng))

    Fig.1 Two Typical hypersonic plans, showing streamwise curvature existence for both external and internal flows (Pictures are from the internet)

    彎曲邊界層流動(dòng)的穩(wěn)定性不同于平板邊界層。基于無(wú)黏分析,Reyleigh最早給出了不同型面彎曲形式下速度剖面與流動(dòng)穩(wěn)定性之間的關(guān)系[1],指出彎曲流動(dòng)產(chǎn)生離心不穩(wěn)定的主要原因是d(Γ2)/dr<0,其中Γ=rU為旋度,U=U(r)為流向速度,r是到圓心的距離??紤]到邊界層內(nèi)的速度分布,凸曲壁對(duì)邊界層具有穩(wěn)定效果,而凹曲壁則會(huì)促進(jìn)邊界層失穩(wěn),如圖2所示。G?rtler[2]通過(guò)線性穩(wěn)定性分析,指出凹曲壁會(huì)在邊界層中產(chǎn)生沿流向?qū)D(zhuǎn)的G?rtler渦[3-9],如圖3所示。G?rtler渦的形成與二次失穩(wěn)破碎主導(dǎo)凹邊界層的失穩(wěn)和轉(zhuǎn)捩過(guò)程[10-11]。在流向曲率的影響下,邊界層內(nèi)流體摩阻系數(shù)、速度型、湍流度、形狀因子、對(duì)流換熱系數(shù)以及雷諾應(yīng)力均會(huì)受到影響,凹曲壁和凸曲壁對(duì)時(shí)均參數(shù)具有相反的作用效果[12-17]。

    亞聲速條件下的研究顯示,凹曲壁能夠增強(qiáng)湍流邊界層的湍流度和壁面摩阻,而凸曲壁則具有抑制湍流脈動(dòng)并減小壁面摩阻的效果。相較于亞聲速流動(dòng)邊界層,超聲速條件下曲壁邊界層的影響因素更為復(fù)雜,邊界層除受主應(yīng)變率?U/?y的影響外,還受到流線彎曲?V/?x、法向壓力梯度?p/?y、流向壓力梯度?p/?x和體積膨脹·U等“額外應(yīng)變率”的綜合影響,這些因素在超聲速條件下都會(huì)有所不同[18]。由于影響因素較多,為深入分析超聲速曲壁邊界層,需要分別探究壓力梯度、曲率本身和可壓縮性對(duì)邊界層的影響規(guī)律,這也使得現(xiàn)有在低速條件下得到的研究結(jié)論難以直接用于超聲速流動(dòng)。此外,超聲速條件下流動(dòng)各向異性特點(diǎn)相較于亞聲速流動(dòng)更加顯著,使得研究單一因素對(duì)超聲速湍流的影響效果變得更為困難[19-20]。激波的存在更會(huì)使湍流流動(dòng)與激波本身產(chǎn)生相互作用,從而導(dǎo)致激波非定常運(yùn)動(dòng)甚至激波面扭曲,穿過(guò)激波的湍流由于熵變也會(huì)發(fā)生相應(yīng)的變化,特別是不同尺度湍流結(jié)構(gòu)對(duì)于激波的響應(yīng)也存在差異。當(dāng)激波強(qiáng)度足夠高時(shí),流動(dòng)還會(huì)發(fā)生分離,進(jìn)一步增加了流動(dòng)的非線性特性。在超聲速條件下,受擾動(dòng)超聲速湍流中不同尺度湍流結(jié)構(gòu)在壓力梯度、額外應(yīng)變率作用下的演化規(guī)律依然不十分清晰,而這些尺度的演化規(guī)律與流動(dòng)的宏觀力學(xué)、熱力學(xué)特性緊密相關(guān),直接影響著湍流建模的準(zhǔn)確性。

    圖2 邊界層的速度型對(duì)應(yīng)于不同曲壁下的穩(wěn)定性分析Fig.2 Stability of boundary layer under the impact of different curvatures

    圖3 凹曲壁邊界層中的G?rtler渦示意圖[10]Fig.3 A schematic of G?rtler vortices generated in the concave boundary layer[10]

    Floryan[21]和Saric[1]曾就彎曲邊界層的穩(wěn)定性和轉(zhuǎn)捩問(wèn)題進(jìn)行過(guò)綜述,但關(guān)于充分發(fā)展的湍流邊界層受流向彎曲影響的問(wèn)題未有提及。早期Spina等[22]在關(guān)于超聲速湍流邊界層結(jié)構(gòu)的綜述中,簡(jiǎn)要敘述了壓力梯度、流向曲率等額外應(yīng)變率的對(duì)邊界層時(shí)均和統(tǒng)計(jì)特性的影響。本文在這些經(jīng)典綜述文章的基礎(chǔ)上,針對(duì)超聲速湍流邊界層對(duì)彎曲壁的響應(yīng)特性進(jìn)行梳理,總結(jié)凸曲壁和凹曲壁兩種流向曲率形式下的超聲速湍流邊界層的時(shí)均和統(tǒng)計(jì)特性,并重點(diǎn)對(duì)近年來(lái)彎曲壁超聲速湍流邊界層的結(jié)構(gòu)演化特征和湍流受擾機(jī)制方面進(jìn)行綜述,提煉發(fā)展趨勢(shì),為后續(xù)研究及其應(yīng)用提供參考。

    1 凸曲壁超聲速湍流邊界層

    超聲速條件下,邊界層受流向凸曲壁的影響,其湍流統(tǒng)計(jì)和結(jié)構(gòu)特性均會(huì)發(fā)生變化。本節(jié)首先將對(duì)凸曲壁超聲速邊界層的典型特征進(jìn)行綜述,然后將對(duì)誘導(dǎo)凸曲壁邊界層湍流衰減的成因——體積膨脹和流向順壓梯度等因素的影響——進(jìn)行總結(jié)。

    1.1 凸曲壁的影響機(jī)制

    與平板邊界層相比,凸曲壁具有衰減超聲速湍流的效果[13, 23-24]。在凸曲率對(duì)超聲速湍流邊界層各類湍流參數(shù)影響的研究中,Bowersox和Buter[25]通過(guò)馬赫數(shù)2.9條件下的實(shí)驗(yàn)發(fā)現(xiàn),凸曲壁對(duì)雷諾剪應(yīng)力的大小和分布均有較大的影響。整個(gè)邊界層內(nèi),剪切應(yīng)力存在幅度達(dá)50%~100%的降低,類似的實(shí)驗(yàn)[26]也報(bào)告了邊界層剪應(yīng)力的大幅下降。Arnette等[27]使用流動(dòng)可視化和激光多普勒測(cè)速(LDV)相結(jié)合的方法,在馬赫數(shù)3的條件下對(duì)比測(cè)試了4個(gè)構(gòu)型的流場(chǎng)特性。Luker等[28]使用LDV對(duì)馬赫數(shù)2.9的凸曲壁邊界層的平均和湍流特性進(jìn)行了詳細(xì)的實(shí)驗(yàn)測(cè)量,并留意到了湍流逆轉(zhuǎn)捩現(xiàn)象。Ekoto等[29]通過(guò)平面粒子圖像測(cè)速技術(shù)(PIV)研究了馬赫數(shù)為2.9時(shí)的兩個(gè)不同曲率的曲面,同時(shí)考慮了表面粗糙度的影響。這些研究表明,受凸曲壁影響的超聲速湍流邊界層的軸向湍流強(qiáng)度可以降低到70%~90%,且減弱程度取決于流向曲率的大小。Luker等[28]觀察到近壁區(qū)域中雷諾切應(yīng)力與平板邊界層相比減小了約25%(如圖 4所示),然而整個(gè)邊界層內(nèi)雷諾切應(yīng)力對(duì)壁面彎曲響應(yīng)卻是不同的。在外部區(qū)域(如y/δ>0.5,其中δ是邊界層厚度),雷諾剪應(yīng)力通常為負(fù),而主應(yīng)變率則保持為正。Arnette[27]和Tichenor等[30]觀察到了類似的趨勢(shì)。Luker等[28]的應(yīng)變率和相關(guān)的額外量測(cè)量表明在邊界層的外側(cè)(y/δ>0.5)的總湍流產(chǎn)生是負(fù)的,這意味著部分能量從波動(dòng)的流場(chǎng)反饋回了平均流。

    圖4 凸曲率和零壓力梯度下湍流特性對(duì)比曲線[28]Fig.4 Comparison of turbulence fluctuations between the convex and zero-pressure-gradient boundary layers[28]

    在超聲速凸曲壁邊界層結(jié)構(gòu)特征的研究方面,Humble等[31]通過(guò)實(shí)驗(yàn)研究了馬赫數(shù)為4.9的湍流邊界層對(duì)凸曲率壁面的響應(yīng)規(guī)律,發(fā)現(xiàn)邊界層外無(wú)旋流動(dòng)進(jìn)入邊界層的頻率明顯降低,湍流/非湍流界面的不規(guī)則性隨曲率的增大而減小, 如圖5所示。Luker[28]等觀察了凸曲率影響下湍流間歇性的增加,同時(shí)間歇流開始朝向邊界層的邊緣移動(dòng),引起負(fù)湍流生成。Wang等[32]則通過(guò)激光平面納米粒子散射技術(shù)發(fā)現(xiàn)超聲速湍流邊界層在進(jìn)入凸曲壁段后,湍流結(jié)構(gòu)的尺度明顯放大(如圖6所示),同時(shí)湍流結(jié)構(gòu)與壁面之間的傾角沿流向增加。Gillis等[33]和Goldfeld等[34]則指出,超聲速湍流在受凸曲壁影響的膨脹和恢復(fù)過(guò)程中,邊界層外層的大尺度結(jié)構(gòu)對(duì)順壓梯度的響應(yīng)較為緩慢,而靠近壁面處的流動(dòng)結(jié)構(gòu)尺寸在順壓梯度的作用下明顯減小,上游大尺度結(jié)構(gòu)的保留被認(rèn)為是湍流生成主要局限于具有較大平均梯度的近壁區(qū)的結(jié)果。

    圖5 流向凸曲率對(duì)超聲速湍流邊界層的影響[31]Fig.5 Impact of the convex wall on the supersonic turbulent boundary layer[31]

    圖6 流向凸曲率影響的超聲速湍流邊界層觀測(cè)[32]Fig.6 Visualization of the supersonic convex turbulent boundary layer[32]

    關(guān)于引起凸曲壁超聲速邊界層湍流統(tǒng)計(jì)特性和結(jié)構(gòu)特征變化的原因一直是研究人員最為關(guān)注的問(wèn)題之一。由于超聲速條件下,彎曲邊界層問(wèn)題的影響因素復(fù)雜,除顯而易見的流向曲率本身的影響之外,壓力梯度、體積膨脹效應(yīng)等均是影響邊界層特性的重要因素[23],近年來(lái)的研究也揭示了體積膨脹和流向順壓梯度對(duì)于超聲速凸曲壁邊界層湍流衰減的貢獻(xiàn)[19, 35]。

    1.2 體積膨脹和順壓梯度的影響

    在以往的許多研究中,研究人員發(fā)現(xiàn)凸曲壁與膨脹拐角表面形成的邊界層存在諸多相似之處,如對(duì)衰減湍流、擴(kuò)展速度剖面中的對(duì)數(shù)區(qū)、較于零壓梯度邊界層對(duì)數(shù)律上移、尾跡區(qū)強(qiáng)度降低等[22,26,36-37]。體積膨脹和流向順壓梯度在誘導(dǎo)邊界層湍流特性發(fā)生以上變化的過(guò)程中發(fā)揮了重要作用。這兩種效應(yīng)在對(duì)于邊界層外層湍流具有相似的衰減效果。但在近壁區(qū),二者對(duì)邊界層的影響則存在較為顯著的區(qū)別:其中體積膨脹主要在膨脹斜坡的拐角處對(duì)近壁湍流存在顯著的衰減作用,而在隨后的發(fā)展中,近壁湍流特征則逐漸恢復(fù)[35];而與體積膨脹的影響不同,流向順壓梯度對(duì)湍流的衰減作用則主要表現(xiàn)在外層,其在整個(gè)壓力梯度影響區(qū)域內(nèi)對(duì)近壁湍流的衰減則較為微弱[19]。

    Knight等[38]指出,湍流衰減是超聲速氣流通過(guò)膨脹拐角的一個(gè)重要特征。Arnette等[39]發(fā)現(xiàn),在膨脹過(guò)程中,湍流水平顯著降低,在靠近壁面處會(huì)更加顯著。在強(qiáng)膨脹作用下,甚至?xí)霈F(xiàn)逆轉(zhuǎn)捩現(xiàn)象,表現(xiàn)為湍動(dòng)能水平的急劇下降和負(fù)雷諾應(yīng)力值的出現(xiàn)。Wang等[32]利用基于納米粒子的平面激光散射方法,在膨脹拐角后的湍流邊界層近壁區(qū)域同樣觀察到這種再層流化的現(xiàn)象,如圖7所示。Bradshaw[18]指出,體積膨脹是影響湍流剪切層的主要因素,尤其是在降低雷諾應(yīng)力方面。Humble等[31]證實(shí),在膨脹過(guò)程中,小尺度結(jié)構(gòu)因體積膨脹而湮滅。Wang等[32]也通過(guò)實(shí)驗(yàn)指出,體積膨脹對(duì)邊界層湍流衰減的貢獻(xiàn)要強(qiáng)于離心力效應(yīng)。Dussauge和Gaviglio[40]用解析和實(shí)驗(yàn)方法研究了體積膨脹對(duì)超聲速湍流邊界層湍流脈動(dòng)的影響。他們得出結(jié)論,體積膨脹是導(dǎo)致雷諾應(yīng)力變化的主要原因。Teramoto等[35]通過(guò)大渦模擬發(fā)現(xiàn)了湍流邊界層的再層流化(如圖8所示,x=0到5之間渦結(jié)構(gòu)明顯減小),同樣指出由體積膨脹所引入的負(fù)湍流生成項(xiàng)導(dǎo)致了該現(xiàn)象的發(fā)生。

    圖7 膨脹拐角表面的超聲速湍流邊界層結(jié)構(gòu)[32]Fig.7 Visualizations of the instantaneous supersonic turbulent boundary layer after expansions

    圖8 膨脹拐角影響下的湍流渦結(jié)構(gòu)變化[35]Fig.8 Turbulence structures under the impact of an expansion ramp[35]

    近期,F(xiàn)ang等[41]利用直接數(shù)值模擬研究了超聲速湍流邊界層在串聯(lián)膨脹-壓縮拐角中的響應(yīng)規(guī)律。超聲速湍流邊界層經(jīng)過(guò)膨脹拐角后出現(xiàn)了明顯的雙層化特點(diǎn),外層湍流結(jié)構(gòu)被明顯抑制(如圖9所示)。Tong等[42]利用直接數(shù)值模擬研究了馬赫數(shù)2.9條件下膨脹拐角對(duì)激波邊界層干擾問(wèn)題的影響。研究發(fā)現(xiàn),膨脹拐角能夠明顯減小激波所致流動(dòng)分離區(qū)的大小及由流動(dòng)分離帶來(lái)的三維效應(yīng),同時(shí)有利于湍流恢復(fù)。該研究同樣觀察到膨脹拐角下游的流場(chǎng)中,邊界層外層湍流水平被明顯抑制,而內(nèi)層湍流仍持續(xù)生成的現(xiàn)象。Sun等[43]通過(guò)DNS(直接數(shù)值模擬)研究了湍流邊界層通過(guò)膨脹拐角后的恢復(fù)機(jī)制,研究表明:與平板邊界層相比,膨脹拐角增加了近壁面區(qū)條帶結(jié)構(gòu)間距,在邊界層內(nèi)部湍動(dòng)能分布和生成呈現(xiàn)出典型的雙層結(jié)構(gòu),近壁區(qū)域內(nèi)湍流在內(nèi)層比在外層擁有更快的恢復(fù)速度。

    在受順壓梯度效應(yīng)的影響的超聲速湍流邊界層中也存在分層現(xiàn)象。揭示流向順壓梯度影響的前提是解耦流向曲率的影響,并避免早期研究中采用彎曲壁面引入壓力梯度的方法,因此需要將曲率效應(yīng)與壓力梯度及體積膨脹效應(yīng)解耦開來(lái)。Wang等[15,20]采用的通過(guò)外源反射膨脹波調(diào)節(jié)壁面附近壓力分布的等效壓力梯度方法為此提供了一個(gè)可靠的解決方案(如圖10)。在此基礎(chǔ)上Wang等[19]通過(guò)對(duì)順壓梯度作用下的平板超聲速湍流邊界層進(jìn)行直接數(shù)值模擬研究發(fā)現(xiàn),與膨脹拐角流動(dòng)中內(nèi)層湍流受到快速抑制不同,當(dāng)湍流邊界層只受到順壓梯度影響時(shí),內(nèi)層湍流衰減的速度要明顯慢于外層,如圖11所示,同為分層,但體積膨脹和流向順壓梯度所誘導(dǎo)形成的分層方式截然不同。

    (a) 瞬態(tài)密度紋影圖

    (b) 瞬態(tài)渦量等值面圖圖9 超聲速串聯(lián)膨脹-壓縮拐角湍流結(jié)構(gòu)[41]Fig.9 Turbulent structures over an expansion- compression corner[41]

    圖10 等效壓力梯度方法產(chǎn)生順壓梯度示意圖[19]flat-plate boundary layer[19]

    最近,Wu等[44]在比較充分發(fā)展湍流邊界層與轉(zhuǎn)捩邊界層的時(shí)空湍流結(jié)構(gòu)時(shí),發(fā)現(xiàn)在零壓梯度充分發(fā)展湍流邊界層中可能存在湍流型湍流渦斑(turbulent-turbulent spots, TUTs)。 Wang等[19]在零壓力梯度的平板邊界層中同樣觀察到類似的TUTs結(jié)構(gòu),并且TUTs結(jié)構(gòu)在順壓梯度作用下,會(huì)受到明顯抑制,如圖12所示。從結(jié)構(gòu)特征上來(lái)看,TUTs結(jié)構(gòu)的減弱也是超聲速凸曲壁邊界層湍流特性減弱的直觀體現(xiàn)。

    圖11 順壓梯度邊界層對(duì)外層湍流的衰減作用[19]

    Fig.11 The supression of turbulence in the outer layer of a supersonic turbulent boundary layer subjected to favorable pressure gradient[19]

    (a) 零壓力梯度邊界層

    (b) 順壓梯度邊界層圖12 順壓梯度下湍流邊界層中的TUTs示意圖[19]Fig.12 TUTs in a supersonic flat-plate turbulent boundary layer with streamwise zero and favorable pressure gradients[19]

    流向凸曲壁能夠誘導(dǎo)邊界層厚度增加、湍流度降低、摩阻系數(shù)減小、對(duì)流換熱系數(shù)減小等。體積膨脹和流向順壓梯度是誘導(dǎo)超聲速湍流邊界層的時(shí)均和統(tǒng)計(jì)特性變化的重要因素。在這兩個(gè)因素的影響下,湍流邊界層的分層響應(yīng)機(jī)制依然是目前研究人員關(guān)心的問(wèn)題,也是進(jìn)一步深入研究超聲速凸曲壁邊界層中湍流結(jié)構(gòu)和流動(dòng)特性的變化規(guī)律的突破點(diǎn)之一。

    2 凹曲壁超聲速湍流邊界層

    和流向凸曲壁邊界層類似,超聲速條件下的流向凹曲壁邊界層中也存在流向逆壓梯度、流向曲率等多重因素影響的問(wèn)題。近年來(lái)的研究逐步揭示了超聲速湍流邊界層對(duì)流向彎曲的響應(yīng)機(jī)制,流向壓力梯度對(duì)湍流結(jié)構(gòu)和統(tǒng)計(jì)特性改變的貢獻(xiàn)也逐漸清晰。本節(jié)將首先對(duì)流向逆壓梯度的影響及其作用機(jī)制進(jìn)行梳理,然后就超聲速湍流邊界層對(duì)凹曲壁的響應(yīng)進(jìn)行總結(jié)。

    2.1 逆壓梯度對(duì)超聲速湍流邊界層的影響

    在早期的研究中,流向壓力梯度主要由凹曲率壁面及壓縮拐角等具有收縮特性的幾何構(gòu)型產(chǎn)生,湍流特性的變化受到流線彎曲的耦合作用。Smith和Smits[45]以及Fernando[46]等首先嘗試分離相互間的干擾,來(lái)對(duì)比研究流向逆壓力梯度和流向彎曲的影響,結(jié)果顯示壓力梯度和純流向曲率對(duì)邊界層時(shí)均和湍流特性的影響均非常顯著,但凹曲率導(dǎo)致的湍流脈動(dòng)則比單純的流向逆壓梯度增強(qiáng)近50%,如圖13所示。

    壓力梯度通常用Clauser壓力梯度參數(shù)β=(dp/dx)(τw/δ*)(其中δ*是位移厚度;τw是壁面剪應(yīng)力;dp/dx是流向壓力梯度)表征[47]。正值β表示逆壓梯度,負(fù)值表示順壓提速。當(dāng)β=0時(shí),表示零壓梯度。由于角動(dòng)量守恒,體積壓縮-·U和橫向擴(kuò)散都可以減小流體單元的展向截面面積,從而增加展向渦量。前者由“額外應(yīng)變率”引起,會(huì)減小邊界層的厚度;后者即使在很小的水平上也可能對(duì)邊界層流動(dòng)產(chǎn)生較大的影響[18, 48-49]。與零壓梯度湍流邊界層相比,逆壓梯度作用下的不可壓縮壁面流動(dòng)沿流動(dòng)方向具有呈降低趨勢(shì)的表面摩擦系數(shù)[47, 50-53],而超聲速條件下,表面摩擦系數(shù)則會(huì)沿流動(dòng)方向不斷增大[54-55]。除了表面摩擦系數(shù)會(huì)在逆壓梯度作用下發(fā)生顯著的改變,平均速度沿邊界層厚度方向的分布相較于零壓梯度邊界層通常也會(huì)發(fā)生顯著的改變。然而由于現(xiàn)有數(shù)據(jù)多來(lái)源于實(shí)驗(yàn)研究,受限于速度測(cè)量的精度,不同實(shí)驗(yàn)得到的速度剖面數(shù)據(jù)往往具有一定出入。根據(jù)文獻(xiàn)[51, 56-57]的報(bào)道,逆壓梯度可能導(dǎo)致不可壓縮流動(dòng)速度剖面的對(duì)數(shù)律區(qū)低于經(jīng)典的對(duì)數(shù)律分布,而文獻(xiàn)[58]則得出即使存在逆壓梯度,經(jīng)典對(duì)數(shù)律依然存在的結(jié)論。對(duì)于超聲速邊界層,文獻(xiàn)[59]觀測(cè)到van-Driest變換速度剖面在逆壓梯度的作用下會(huì)發(fā)生相對(duì)于零壓梯度對(duì)數(shù)律區(qū)的上移。文獻(xiàn)[46]則表示對(duì)數(shù)律區(qū)基本維持不變。

    (a) 超聲速逆壓梯度平板邊界層在不同流向位置的湍流度剖面

    (b) 超聲速凹曲壁邊界層在不同流向位置的湍流度剖面圖13 逆壓梯度和凹曲壁條件下超聲速邊界層的在 不同流向位置湍流度分布對(duì)比[45]

    Fig.13 A comparison of the turbulence intensity between the adverse - pressure - gradient and concave boundary layers[45]

    Wang等[15, 20]通過(guò)采用等效壓力梯度方法來(lái)分離流向壓力梯度和流向曲率的影響,并開展的實(shí)驗(yàn)(如圖 14所示)[20]和直接數(shù)值模擬研究[60],得出了相似的結(jié)論。除了對(duì)數(shù)律區(qū)的可能偏移外,文獻(xiàn)[46, 61]還發(fā)現(xiàn)受逆壓梯度的影響,對(duì)數(shù)律區(qū)的斜率也會(huì)發(fā)生變化,這可能是由于邊界層內(nèi)湍流運(yùn)動(dòng)的長(zhǎng)度尺度增加導(dǎo)致的。在逆壓梯度的影響下,速度剖面的尾跡區(qū)相較于零壓梯度邊界層明顯上移。逆壓梯度能顯著地促進(jìn)不穩(wěn)定性的增長(zhǎng)并加速流動(dòng)轉(zhuǎn)捩,當(dāng)逆壓梯度效應(yīng)足夠強(qiáng)時(shí),邊界層還可能發(fā)生流動(dòng)分離[36, 62-63]。與零壓梯度邊界層流動(dòng)相比,逆壓梯度作用下的高階統(tǒng)計(jì)特征和流動(dòng)結(jié)構(gòu)會(huì)發(fā)生明顯變化。湍動(dòng)能和雷諾應(yīng)力通常會(huì)被顯著放大,特別是在外層,當(dāng)用壁面單位無(wú)量綱化湍流脈動(dòng)時(shí),邊界層外層區(qū)域可以觀察到第二峰值[50-51, 64]。

    圖14 等效逆壓梯度超聲速湍流邊界層實(shí)驗(yàn)設(shè)計(jì)示意圖[20]

    Fig.14 A schematic of experimental setup to generate streamwise adverse pressure gradient for a supersonic flat-plate boundary layer[20]

    湍流結(jié)構(gòu)特征在逆壓梯度的作用下也與零壓梯度邊界層不同。具體表現(xiàn)為大尺度結(jié)構(gòu)(large-scale motions)在邊界層內(nèi)更為活躍,渦結(jié)構(gòu)能量增益顯著[52, 53, 65-68]。由大尺度渦結(jié)構(gòu)產(chǎn)生的上噴、下掃和展向事件被認(rèn)為是導(dǎo)致平均流動(dòng)特性和湍流統(tǒng)計(jì)量改變的主要誘因。Bobke等[69]研究了在恒定和可變逆壓梯度下,不同摩擦雷諾數(shù)湍流邊界層的歷史效應(yīng)。結(jié)果表明,平均速度和雷諾應(yīng)力分布依賴于邊界層流動(dòng)的下游發(fā)展。由于存在歷史效應(yīng),逆壓梯度對(duì)湍流邊界層有累積效應(yīng)。此外,針對(duì)凹曲壁超聲速湍流邊界層和逆壓梯度作用下的平板超聲速湍流邊界層的數(shù)值模擬[14, 60, 70]均發(fā)現(xiàn)了類似于湍流型湍流渦斑的流場(chǎng)結(jié)構(gòu)(如圖15所示),這種結(jié)構(gòu)在逆壓梯度的作用下呈現(xiàn)出增強(qiáng)的趨勢(shì),其與邊界層外層大尺度結(jié)構(gòu)的增強(qiáng)之間的關(guān)聯(lián)仍有待進(jìn)一步研究。

    圖15 (a)逆壓梯度作用下湍流邊界層中的類湍流型湍流渦斑(TUTs)結(jié)構(gòu);(b)零壓梯度作用下湍流邊界層中的類湍流型湍流渦斑結(jié)構(gòu)[60]

    Fig.15 (a) TUTs in a supersonic adverse-pressure-gradient boundary layer; (b) TUTs in a supersonic zero-pressure-gradient boundary layer[60]

    流向逆壓梯度對(duì)邊界層時(shí)均和湍流特性的顯著影響源于其對(duì)邊界層中法向輸運(yùn)的增強(qiáng)。在凹曲壁超聲速邊界層流動(dòng)中,邊界層除受到壓力梯度的影響外,由流向彎曲帶來(lái)的離心不穩(wěn)定的作用會(huì)進(jìn)一步加強(qiáng)邊界層的法向輸運(yùn),并改變邊界層特性。

    2.2 凹曲壁超聲速湍流邊界層

    近些年來(lái),隨著實(shí)驗(yàn)測(cè)試和數(shù)值仿真手段的改善,關(guān)于凹曲壁超聲速湍流邊界層的認(rèn)識(shí)在湍流統(tǒng)計(jì)特性和湍流結(jié)構(gòu)的演化機(jī)制方面均逐步深入。

    實(shí)驗(yàn)研究方面,Jayaram等[61]在自由來(lái)流馬赫數(shù)為2.87的條件下測(cè)試了流向凹曲壁和壓縮拐角上形成的邊界層,他們發(fā)現(xiàn)這兩種構(gòu)型均能夠顯著增強(qiáng)邊界層的湍流特性,致使其流向速度型在對(duì)數(shù)律層存在一個(gè)凹陷,低于正常的平板邊界層的對(duì)數(shù)律層,他們將其歸因?yàn)榘记谶吔鐚又型牧鬟\(yùn)動(dòng)的長(zhǎng)度和時(shí)間尺度的增加所誘導(dǎo)的速度分布的變化。此外,Hoffmann等[71]也在實(shí)驗(yàn)中發(fā)現(xiàn)了類似的現(xiàn)象。Donovan等[55]在馬赫數(shù)2.86條件下發(fā)現(xiàn),與平板邊界層相比,凹曲壁邊界層壁面切應(yīng)力增長(zhǎng)了約125%。不同的是,他們?cè)趯?shí)驗(yàn)的彎曲段并沒(méi)有發(fā)現(xiàn)邊界層速度剖面中存在的凹陷現(xiàn)象,但卻在彎曲段下游的恢復(fù)段中發(fā)現(xiàn)了速度剖面中的凹陷。Wang等[16]在實(shí)驗(yàn)中發(fā)現(xiàn)超聲速平板邊界層進(jìn)入彎曲段后,原本的較大尺度結(jié)構(gòu)迅速轉(zhuǎn)變?yōu)榇罅啃〕叨冉Y(jié)構(gòu),如圖 16所示。

    圖16 實(shí)驗(yàn)觀測(cè)到的超聲速凹曲壁邊界層中的 瞬時(shí)流場(chǎng)結(jié)構(gòu)[16]Fig.16 Instantaneous turbulent structures visualized in experiment[16]

    近期,F(xiàn)laherty和Austin[72]對(duì)最高馬赫數(shù)達(dá)到7.45的高超聲速曲壁邊界層開展了研究。紋影結(jié)果顯示,相較于平板邊界層,凹曲壁段的邊界層厚度可以看到明顯變薄的現(xiàn)象,如圖17所示。此外,他們還對(duì)不同曲率彎曲壁面的近壁熱流進(jìn)行了詳細(xì)測(cè)量,結(jié)果表面,相較于平板邊界層,凹曲壁能夠顯著提高壁面的換熱系數(shù)。并且他們發(fā)現(xiàn),對(duì)于不同曲率的邊界層,當(dāng)以流向偏轉(zhuǎn)角為橫坐標(biāo)時(shí),不同曲壁邊界層的壁面熱流數(shù)據(jù)分布具有自相似性,如圖18所示,這對(duì)于工程設(shè)計(jì)具有重要意義,意味著工程中可以基于壁面偏轉(zhuǎn)角來(lái)對(duì)彎曲壁面的表面熱流進(jìn)行估算。

    (a) 平板邊界層

    (b) 凹曲壁邊界層圖17 平板邊界層和凹曲壁邊界層的紋影圖[72]Fig.17 Schlieren of flat-plate boundary layer and concave boundary layer[72]

    邊界層時(shí)均和湍流特性的變化主要來(lái)源于流場(chǎng)結(jié)構(gòu)的變化。對(duì)于凹曲壁超聲速層流邊界層,G?rtler渦的形成及其二次失穩(wěn)是邊界層轉(zhuǎn)捩的重要因素,符松課題組[6-7, 10]采用穩(wěn)定性分析研究了超聲速層流邊界層在曲率影響下其中G?rtler渦的非線性增長(zhǎng)過(guò)程,如圖19所示,揭示了不同模態(tài)波的影響和作用機(jī)制。

    對(duì)于湍流邊界層,由于流場(chǎng)結(jié)構(gòu)的復(fù)雜性,其中的G?rtler不穩(wěn)定難以直接辨識(shí),但是邊界層在進(jìn)入彎曲段后,其中低速和高速流動(dòng)沿展向的交錯(cuò)與周期性分布表明,即使在湍流邊界層中,G?rtler不穩(wěn)定對(duì)流場(chǎng)結(jié)構(gòu)的影響依然顯著,如圖20所示。Hoffmann等[71]以及Barlow和Johnston[73-74]均通過(guò)實(shí)驗(yàn)中發(fā)現(xiàn)G?rtler不穩(wěn)定能夠在邊界層中引入大尺度的流向渦旋結(jié)構(gòu),很多研究中也稱其為類G?rtler渦,這些結(jié)構(gòu)能夠誘導(dǎo)邊界層中形成明顯強(qiáng)于一般平板邊界層的內(nèi)外層之間的法向?qū)α髋c動(dòng)量交換運(yùn)動(dòng)。但是不同于層流邊界層中產(chǎn)生的G?rtler渦,湍流邊界層中這些流向渦旋結(jié)構(gòu)并沒(méi)有對(duì)應(yīng)的渦核和集中渦量。

    圖19 馬赫數(shù)4.5邊界層中G?rtler渦及低速與 高度條帶的形成過(guò)程[10]Fig.19 Formation of G?rtler vortices and the streaks within a Mach 4.5 boundary layer[10]

    (a) 瞬時(shí)渦結(jié)構(gòu)

    (b) 時(shí)均摩阻分布圖20 凹曲壁邊界層中瞬時(shí)渦結(jié)構(gòu)和時(shí)均摩阻分布云圖[75]

    Fig.20 Instantaneous vortices and time-averaged wall friction for a supersonic concave boundary layer[75]

    Tong等[75]在馬赫數(shù)2.9條件下的直接數(shù)值模擬研究也發(fā)現(xiàn)了平板邊界層在進(jìn)入到彎曲段之后,邊界層中會(huì)形成明顯大尺度的類G?rtler渦結(jié)構(gòu),同時(shí)這種大尺度結(jié)構(gòu)會(huì)誘導(dǎo)邊界層的時(shí)均摩阻分布形成展向周期性分布(如圖20所示),這也表明類G?rtler渦結(jié)構(gòu)的展向位置是總體穩(wěn)定的。但是,需要強(qiáng)調(diào)的是,從渦的任何定義來(lái)看,這些大尺度的類G?rtler渦結(jié)構(gòu)都難以稱之為單一的渦。此外,Tong等[75]模擬結(jié)果還顯示流向曲率對(duì)邊界層對(duì)數(shù)率層的影響較小,這也表明,在關(guān)于流向曲率對(duì)邊界層對(duì)數(shù)率層的影響方面依然存在爭(zhēng)論。此外,他們的數(shù)值模擬結(jié)果同樣顯示了邊界層中湍流脈動(dòng)的顯著增強(qiáng),其中在邊界層外層的增長(zhǎng)尤為明顯。

    近年來(lái),Wang等[14,16-17]和Sun等[70]通過(guò)系統(tǒng)的實(shí)驗(yàn)和直接數(shù)值模擬研究,揭示了G?rtler不穩(wěn)定誘導(dǎo)邊界層中形成大尺度運(yùn)動(dòng)的機(jī)制。

    Wang等[14]指出,邊界層中大尺度結(jié)構(gòu)的變化是導(dǎo)致其湍流統(tǒng)計(jì)特性變化的直接原因。湍流邊界層超級(jí)結(jié)構(gòu)(VLSM)中存在豐富的發(fā)卡渦結(jié)構(gòu),這些結(jié)構(gòu)在G?rtler不穩(wěn)定的作用下其強(qiáng)度會(huì)明顯增強(qiáng),并同時(shí)誘導(dǎo)形成更多的流向結(jié)構(gòu)。這些結(jié)構(gòu)的綜合作用顯著增強(qiáng)了邊界層中的動(dòng)量交換,放大平板邊界層中VLSM結(jié)構(gòu)并導(dǎo)致了展向不均勻,形成了Barlow和Johnston[73-74]在凹曲壁邊界層中觀測(cè)到的大尺度的類G?rtler渦結(jié)構(gòu)。凹曲壁邊界層中大尺度法向運(yùn)動(dòng)的形成,進(jìn)一步促使邊界層中上噴(Q2)和下洗(Q4)事件明顯增強(qiáng),由于Q2和Q4事件是邊界層中湍流生成的最主要貢獻(xiàn)者,邊界層的湍流脈動(dòng)也隨之增強(qiáng),并在邊界層外層形成了湍流脈動(dòng)的第二峰值[14],如圖 21所示。

    Sun等[70]指出,凹曲壁上超聲速湍流邊界層內(nèi)部存在的大尺度類G?rtler結(jié)構(gòu)會(huì)扭曲當(dāng)?shù)孛芏忍荻龋c密度和壓力梯度有關(guān)的斜壓渦量生成項(xiàng)相比逆壓力梯度邊界層會(huì)顯著增強(qiáng),該項(xiàng)對(duì)小尺度渦的形成具有顯著的貢獻(xiàn),這一機(jī)制同樣誘導(dǎo)了流向逆壓梯度邊界層中大量小尺度湍流結(jié)構(gòu)的形成,如圖22所示。關(guān)于G?rtler不穩(wěn)定在超聲速凹曲壁邊界層中增強(qiáng)并誘導(dǎo)形成更多流向結(jié)構(gòu)的發(fā)現(xiàn)也進(jìn)一步解釋了Wang等[16]在實(shí)驗(yàn)中發(fā)現(xiàn)的超聲速平板邊界層進(jìn)入彎曲段后出現(xiàn)大量小尺度結(jié)構(gòu)(如圖22所示)的原因。

    (a) 流向凹曲壁邊界層中的瞬時(shí)大尺度結(jié)構(gòu)

    (b) 流向湍流脈動(dòng)在不同流向位置的剖面圖21 流向凹曲壁邊界層中的瞬時(shí)大尺度結(jié)構(gòu)及其 誘導(dǎo)形成的外層湍流脈動(dòng)第二峰值[14]

    Fig.21 Instantaneous structures of a supersonic concave boundary layer and the secondary bump induced by the large scale motions[14]

    (a) G?rtler不穩(wěn)定引起當(dāng)?shù)靥荻茸兓?/p>

    (b) 斜壓效應(yīng)

    (c) 類G?rtler渦促進(jìn)了外層的交換圖22 超聲速凹曲壁邊界層中的G?rtler不穩(wěn)定及由斜壓渦量 生成項(xiàng)誘導(dǎo)形成類G?rtler渦結(jié)構(gòu)的原理簡(jiǎn)圖[70]Fig.22 Schematic of G?rtler instabilities and the small vortices formed by the large baroclinic production[70]

    總的來(lái)看,流向凹曲壁作用于超聲速湍流邊界層會(huì)增大摩阻系數(shù)、減小形狀因子、增強(qiáng)對(duì)流換熱系數(shù)等[45, 54-55, 71-72, 76-78]。由于凹曲壁誘導(dǎo)的湍流結(jié)構(gòu)特性的改變,尤其是大尺度類G?rtler結(jié)構(gòu)的形成是導(dǎo)致時(shí)均和湍流統(tǒng)計(jì)特性改變和關(guān)鍵,這方面的研究目前正獲得越來(lái)越多的關(guān)注。

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

    多年來(lái)國(guó)內(nèi)外研究人員關(guān)于超聲速邊界層中流向曲率和壓力梯度影響方面的研究取得了豐碩的成果。在早期發(fā)現(xiàn)的流向凸曲壁對(duì)湍流衰減作用和流向凹曲壁對(duì)湍流增強(qiáng)作用的基礎(chǔ)上,近年來(lái)的研究在流向彎曲壁影響超聲速湍流邊界層的物理機(jī)制方面取得了顯著進(jìn)展:揭示了體積膨脹效應(yīng)在超聲速湍流衰減中的重要作用,其在凸曲壁邊界層中的影響強(qiáng)于離心效應(yīng)和流向順壓梯度;發(fā)現(xiàn)了在凹曲壁邊界層中導(dǎo)致湍流特性增強(qiáng)的多重影響因素,包括G?rtler不穩(wěn)定對(duì)流向結(jié)構(gòu)的放大效應(yīng)、斜壓渦量生成以及流向逆壓梯度作用等,揭示了超聲速湍流邊界層對(duì)流向凹曲壁的響應(yīng)機(jī)制。這些研究成果對(duì)于湍流?;桶l(fā)展基于型面優(yōu)化的流動(dòng)控制技術(shù)具有指導(dǎo)意義。

    但是總體來(lái)看,目前不論是從理論認(rèn)識(shí)還是從指際工程實(shí)踐方面來(lái)看,依然存在諸多問(wèn)題需要進(jìn)一步深入研究:

    1)近壁湍流對(duì)流向曲率和壓力梯度響應(yīng)規(guī)律。目前的研究多關(guān)注于外層湍流結(jié)構(gòu)演化特性,但考慮到湍流邊界層的維持依賴于近壁湍流的自持循環(huán),揭示近壁湍流在流向曲率和壓力梯度影響下的響應(yīng)規(guī)律,對(duì)于深刻認(rèn)識(shí)流向彎曲壁超聲速湍流邊界層并對(duì)其演化進(jìn)行預(yù)測(cè)具有重要意義。

    2)高馬赫數(shù)、高雷諾數(shù)、不同壁面溫度條件下流向彎曲壁對(duì)超聲速湍流邊界層的影響。目前研究結(jié)論主要源于馬赫數(shù)3以下的實(shí)驗(yàn)和仿真,不論是從理論研究還是從指導(dǎo)工程實(shí)踐角度來(lái)看,關(guān)于邊界層湍流特性、壁面熱流、壁面摩阻等隨馬赫數(shù)、雷諾數(shù)、壁面溫度等參數(shù)變化規(guī)律的研究都亟待開展。

    3)超聲速三維流場(chǎng)中流向彎曲壁的影響機(jī)制。在工程問(wèn)題中,流向彎曲的影響往往同激波/邊界層干擾、展向彎曲等共同存在,三維性強(qiáng),相互作用機(jī)制復(fù)雜。進(jìn)一步需要針對(duì)激波與彎曲邊界層相互作用、多彎曲擾動(dòng)條件下的超聲速湍流邊界層發(fā)展特性等問(wèn)題開展研究,這對(duì)于高超聲速飛行器內(nèi)外流道的設(shè)計(jì)具有重要意義。

    致謝:感謝清華大學(xué)符松教授、中科院力學(xué)所李新亮研究員以及國(guó)防科技大學(xué)趙玉新教授在撰寫本文過(guò)程中的探討以及提供的資料。

    猜你喜歡
    壓力梯度邊界層流向
    小溪??!流向遠(yuǎn)方
    井岡教育(2020年6期)2020-12-14 03:04:42
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    壓力梯度在油田開發(fā)中的應(yīng)用探討
    十大漲幅、換手、振副、資金流向
    流向逆轉(zhuǎn)的啟示
    疊加原理不能求解含啟動(dòng)壓力梯度滲流方程
    一類具有邊界層性質(zhì)的二次奇攝動(dòng)邊值問(wèn)題
    非特征邊界的MHD方程的邊界層
    致密砂巖啟動(dòng)壓力梯度數(shù)值的影響因素
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    鄭州市春季邊界層風(fēng)氣候變化研究
    河南科技(2014年23期)2014-02-27 14:19:08
    免费人妻精品一区二区三区视频| 亚洲精品国产色婷婷电影| 中文字幕人妻熟人妻熟丝袜美| 大片免费播放器 马上看| 日韩av免费高清视频| 日本午夜av视频| 日韩一本色道免费dvd| 欧美精品一区二区大全| 亚洲成色77777| 亚洲人成77777在线视频| 国产精品不卡视频一区二区| 国产乱来视频区| 久久久久国产精品人妻一区二区| 99久国产av精品国产电影| 亚洲精品av麻豆狂野| 国产日韩一区二区三区精品不卡 | 久久久国产一区二区| 欧美人与性动交α欧美精品济南到 | 考比视频在线观看| 精品久久久噜噜| 老女人水多毛片| 自线自在国产av| 久久久久久久久久久丰满| 麻豆精品久久久久久蜜桃| 亚洲欧美一区二区三区国产| 国产老妇伦熟女老妇高清| 午夜91福利影院| 午夜福利视频精品| 午夜福利视频精品| 热re99久久国产66热| 久久 成人 亚洲| 国产成人免费无遮挡视频| 精品一区在线观看国产| 777米奇影视久久| av免费在线看不卡| 草草在线视频免费看| 97在线人人人人妻| 三级国产精品片| 国产精品久久久久久精品古装| 狂野欧美激情性xxxx在线观看| 99热国产这里只有精品6| 国产精品一区二区在线观看99| 五月开心婷婷网| 大香蕉久久网| www.av在线官网国产| videosex国产| 国产av码专区亚洲av| 日本-黄色视频高清免费观看| 亚洲成色77777| 一级毛片黄色毛片免费观看视频| 亚洲第一区二区三区不卡| 永久网站在线| 亚洲精品视频女| 母亲3免费完整高清在线观看 | 中文精品一卡2卡3卡4更新| 啦啦啦在线观看免费高清www| 亚洲内射少妇av| 久久影院123| xxxhd国产人妻xxx| av国产久精品久网站免费入址| 99九九在线精品视频| 能在线免费看毛片的网站| 亚洲国产色片| 51国产日韩欧美| 国产精品一二三区在线看| 香蕉精品网在线| 久久久久久久久久成人| 日韩成人伦理影院| 免费大片黄手机在线观看| 久久精品人人爽人人爽视色| 午夜福利在线观看免费完整高清在| 亚洲欧美一区二区三区黑人 | 99热全是精品| 下体分泌物呈黄色| 国产无遮挡羞羞视频在线观看| 内地一区二区视频在线| 午夜福利视频在线观看免费| 97在线人人人人妻| 国产淫语在线视频| 亚洲无线观看免费| 不卡视频在线观看欧美| 国产高清不卡午夜福利| 欧美日韩一区二区视频在线观看视频在线| 国产高清不卡午夜福利| 免费观看在线日韩| 美女大奶头黄色视频| 日本与韩国留学比较| 日韩伦理黄色片| 亚洲一级一片aⅴ在线观看| av又黄又爽大尺度在线免费看| 亚洲欧洲日产国产| av在线播放精品| 国产精品99久久99久久久不卡 | 亚洲国产av影院在线观看| 成人国产av品久久久| 丰满饥渴人妻一区二区三| 麻豆精品久久久久久蜜桃| 欧美成人午夜免费资源| 国产一区有黄有色的免费视频| 国产亚洲av片在线观看秒播厂| 亚洲少妇的诱惑av| 久久久国产欧美日韩av| 一级爰片在线观看| 国产精品不卡视频一区二区| 免费看光身美女| 人人妻人人澡人人爽人人夜夜| 久久99一区二区三区| 亚洲精品亚洲一区二区| 久久热精品热| 亚洲精品自拍成人| 哪个播放器可以免费观看大片| 亚洲欧美成人综合另类久久久| 日韩大片免费观看网站| 男女国产视频网站| 好男人视频免费观看在线| 男人爽女人下面视频在线观看| 国产精品国产三级国产专区5o| 中文字幕制服av| 亚洲第一av免费看| 考比视频在线观看| 女人久久www免费人成看片| 久久人妻熟女aⅴ| 你懂的网址亚洲精品在线观看| 欧美xxⅹ黑人| 欧美亚洲 丝袜 人妻 在线| 又大又黄又爽视频免费| 国产精品免费大片| 在线观看美女被高潮喷水网站| 少妇猛男粗大的猛烈进出视频| 成年av动漫网址| 日本欧美视频一区| 国精品久久久久久国模美| 国产免费现黄频在线看| 满18在线观看网站| 国产精品国产三级专区第一集| 亚洲伊人久久精品综合| 精品国产一区二区久久| xxx大片免费视频| 人人妻人人添人人爽欧美一区卜| 一区二区三区乱码不卡18| 97精品久久久久久久久久精品| 亚洲精品国产av蜜桃| 黄色欧美视频在线观看| kizo精华| 熟女电影av网| 亚洲av国产av综合av卡| 少妇被粗大猛烈的视频| 五月玫瑰六月丁香| 国产精品一区二区在线不卡| 国产精品久久久久久av不卡| 一级毛片我不卡| 久久人人爽人人片av| 超碰97精品在线观看| 亚洲国产精品成人久久小说| 9色porny在线观看| 一级毛片 在线播放| 免费人妻精品一区二区三区视频| 国产在线视频一区二区| 国产色爽女视频免费观看| 三级国产精品片| 国产精品熟女久久久久浪| 久久青草综合色| 多毛熟女@视频| 久久久久久伊人网av| 91aial.com中文字幕在线观看| 天堂中文最新版在线下载| 人人妻人人爽人人添夜夜欢视频| 国产精品人妻久久久久久| 日韩视频在线欧美| 欧美成人午夜免费资源| 91精品国产九色| 国产精品一区www在线观看| 中文字幕亚洲精品专区| 尾随美女入室| 熟女电影av网| 亚洲熟女精品中文字幕| 黑人高潮一二区| 大片免费播放器 马上看| 一级黄片播放器| 久久久久久久亚洲中文字幕| 女性生殖器流出的白浆| 久久久久久久久久成人| 国产永久视频网站| 色婷婷av一区二区三区视频| 黑人欧美特级aaaaaa片| 国产片特级美女逼逼视频| 久久久久久人妻| 午夜福利,免费看| 日韩人妻高清精品专区| 美女国产视频在线观看| 观看美女的网站| 欧美xxxx性猛交bbbb| 极品少妇高潮喷水抽搐| 亚洲性久久影院| 国产成人精品久久久久久| 搡女人真爽免费视频火全软件| 国产高清国产精品国产三级| 中文字幕精品免费在线观看视频 | 国产精品久久久久成人av| 99热6这里只有精品| 日本91视频免费播放| 精品人妻熟女毛片av久久网站| 永久免费av网站大全| 精品国产一区二区久久| 免费大片18禁| a级毛片黄视频| 国产日韩欧美亚洲二区| 日本黄色日本黄色录像| 在线看a的网站| 久久精品国产自在天天线| 久久综合国产亚洲精品| 我要看黄色一级片免费的| 久久精品国产a三级三级三级| 人人妻人人爽人人添夜夜欢视频| 我的老师免费观看完整版| 国产一区二区在线观看av| 国产片内射在线| 伦理电影免费视频| 久久久国产一区二区| 亚州av有码| 又粗又硬又长又爽又黄的视频| 99久久中文字幕三级久久日本| 欧美日韩精品成人综合77777| 欧美亚洲 丝袜 人妻 在线| 18禁在线播放成人免费| 91久久精品电影网| 欧美日韩综合久久久久久| 久久久久久久久久成人| 欧美老熟妇乱子伦牲交| 高清在线视频一区二区三区| 自线自在国产av| 18禁动态无遮挡网站| 亚洲精品456在线播放app| 成人二区视频| 狂野欧美白嫩少妇大欣赏| 欧美xxxx性猛交bbbb| 亚洲精品第二区| 久久精品夜色国产| 午夜影院在线不卡| .国产精品久久| 国产精品无大码| a级毛片免费高清观看在线播放| 亚洲精品日本国产第一区| 老熟女久久久| 亚洲国产日韩一区二区| 日韩av在线免费看完整版不卡| 精品久久蜜臀av无| 少妇猛男粗大的猛烈进出视频| 国产精品蜜桃在线观看| 国产精品人妻久久久久久| 91午夜精品亚洲一区二区三区| 丰满乱子伦码专区| 国产老妇伦熟女老妇高清| 亚洲国产成人一精品久久久| 亚洲人成网站在线观看播放| 大话2 男鬼变身卡| 亚洲美女黄色视频免费看| 9色porny在线观看| 韩国av在线不卡| 赤兔流量卡办理| 国产免费福利视频在线观看| 亚洲综合色惰| 丰满少妇做爰视频| 亚洲精品一区蜜桃| 国产成人精品福利久久| 免费观看a级毛片全部| 黄片播放在线免费| 人妻系列 视频| 久热这里只有精品99| 人妻夜夜爽99麻豆av| 亚洲av.av天堂| 亚洲色图综合在线观看| 婷婷成人精品国产| a级片在线免费高清观看视频| 中文天堂在线官网| 97超视频在线观看视频| 99九九在线精品视频| 欧美成人午夜免费资源| videosex国产| 国产精品秋霞免费鲁丝片| 国产精品久久久久久av不卡| 亚洲美女搞黄在线观看| videosex国产| 精品一区二区免费观看| 精品少妇久久久久久888优播| 丰满饥渴人妻一区二区三| 国产成人午夜福利电影在线观看| 成人毛片a级毛片在线播放| 亚洲一级一片aⅴ在线观看| av卡一久久| 999精品在线视频| 精品人妻熟女毛片av久久网站| 国产午夜精品久久久久久一区二区三区| 日韩强制内射视频| 99热全是精品| 在线精品无人区一区二区三| 国产欧美日韩综合在线一区二区| 啦啦啦视频在线资源免费观看| 婷婷色av中文字幕| kizo精华| 国产黄色免费在线视频| 成人18禁高潮啪啪吃奶动态图 | 在线观看人妻少妇| 肉色欧美久久久久久久蜜桃| 亚洲一级一片aⅴ在线观看| 精品卡一卡二卡四卡免费| 亚洲美女搞黄在线观看| 麻豆成人av视频| 午夜av观看不卡| 久久人人爽人人爽人人片va| 久久影院123| 久久国产亚洲av麻豆专区| 热99国产精品久久久久久7| 久久99热6这里只有精品| 久久精品国产a三级三级三级| 成人综合一区亚洲| 日韩欧美一区视频在线观看| 亚洲在久久综合| 国产有黄有色有爽视频| 日韩人妻高清精品专区| 3wmmmm亚洲av在线观看| 国产精品人妻久久久影院| 成人影院久久| 美女福利国产在线| 亚洲国产精品成人久久小说| 一区二区三区四区激情视频| 亚洲精品,欧美精品| 赤兔流量卡办理| 看十八女毛片水多多多| 久久精品国产a三级三级三级| 欧美日韩国产mv在线观看视频| 狠狠婷婷综合久久久久久88av| 精品久久久精品久久久| 老司机亚洲免费影院| 国产亚洲精品第一综合不卡 | 岛国毛片在线播放| 99久久中文字幕三级久久日本| 亚洲美女搞黄在线观看| 久久久久国产网址| 成人国产av品久久久| 视频在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 看免费成人av毛片| 日本vs欧美在线观看视频| 2021少妇久久久久久久久久久| 成人午夜精彩视频在线观看| 日本色播在线视频| 最近手机中文字幕大全| 久久狼人影院| 成人黄色视频免费在线看| 亚洲婷婷狠狠爱综合网| 日韩精品免费视频一区二区三区 | 99国产精品免费福利视频| 少妇人妻 视频| 精品人妻熟女av久视频| 精品人妻在线不人妻| 免费观看a级毛片全部| 一级毛片我不卡| av专区在线播放| 亚洲精品日韩av片在线观看| av电影中文网址| 免费看不卡的av| www.色视频.com| 亚洲精品日本国产第一区| 亚洲久久久国产精品| 久久久久视频综合| 国产精品一二三区在线看| 性高湖久久久久久久久免费观看| 啦啦啦在线观看免费高清www| 妹子高潮喷水视频| 精品国产露脸久久av麻豆| 激情五月婷婷亚洲| 亚洲人成网站在线播| 永久网站在线| 国产综合精华液| 91成人精品电影| 亚洲精品日韩在线中文字幕| 丝袜美足系列| 精品人妻在线不人妻| kizo精华| 国产欧美日韩综合在线一区二区| 日韩亚洲欧美综合| 午夜福利,免费看| 各种免费的搞黄视频| 91精品国产九色| 亚洲第一区二区三区不卡| 91精品国产九色| 永久网站在线| 久久97久久精品| 成人综合一区亚洲| av免费观看日本| 狂野欧美激情性bbbbbb| 亚洲中文av在线| xxx大片免费视频| 人妻 亚洲 视频| av在线观看视频网站免费| 免费观看的影片在线观看| 欧美激情国产日韩精品一区| 成人国语在线视频| 国产精品国产三级国产专区5o| 人人妻人人澡人人爽人人夜夜| 欧美国产精品一级二级三级| 国产深夜福利视频在线观看| 成人毛片60女人毛片免费| 成人18禁高潮啪啪吃奶动态图 | 看十八女毛片水多多多| av又黄又爽大尺度在线免费看| 免费看av在线观看网站| 免费观看av网站的网址| 久久ye,这里只有精品| 精品亚洲乱码少妇综合久久| 久久久亚洲精品成人影院| 91精品伊人久久大香线蕉| 午夜视频国产福利| 男女边摸边吃奶| 高清不卡的av网站| 成年女人在线观看亚洲视频| av免费观看日本| 狂野欧美激情性bbbbbb| 免费高清在线观看日韩| 综合色丁香网| 久久精品熟女亚洲av麻豆精品| 免费播放大片免费观看视频在线观看| 在现免费观看毛片| 成年人免费黄色播放视频| 亚洲精品国产色婷婷电影| 午夜福利视频在线观看免费| 婷婷色麻豆天堂久久| 国精品久久久久久国模美| 一区二区日韩欧美中文字幕 | 国产欧美亚洲国产| 日韩三级伦理在线观看| 91精品国产九色| 22中文网久久字幕| 一级毛片我不卡| 精品少妇黑人巨大在线播放| 国产男女内射视频| 精品亚洲乱码少妇综合久久| 国产精品无大码| 人妻少妇偷人精品九色| videosex国产| 久久影院123| 青春草视频在线免费观看| 精品久久久久久电影网| 日本黄色片子视频| 国产免费一级a男人的天堂| av电影中文网址| 不卡视频在线观看欧美| 视频区图区小说| 国产乱来视频区| 我的老师免费观看完整版| 亚洲国产最新在线播放| 欧美激情 高清一区二区三区| 日韩中文字幕视频在线看片| 99久久人妻综合| 国产 一区精品| 国产老妇伦熟女老妇高清| 久久久久久久久久久免费av| 观看美女的网站| 精品人妻在线不人妻| 国产成人91sexporn| 99热这里只有是精品在线观看| 日本黄大片高清| 精品国产国语对白av| 亚洲,欧美,日韩| 亚洲精品国产色婷婷电影| 久久婷婷青草| av又黄又爽大尺度在线免费看| 国产精品久久久久久精品古装| 久久久亚洲精品成人影院| 国产在视频线精品| 人人妻人人澡人人爽人人夜夜| 亚洲国产欧美在线一区| 十八禁高潮呻吟视频| 少妇丰满av| 妹子高潮喷水视频| 美女大奶头黄色视频| 精品国产国语对白av| 国产日韩一区二区三区精品不卡 | 十八禁高潮呻吟视频| 看非洲黑人一级黄片| a级毛色黄片| 一级,二级,三级黄色视频| 国产老妇伦熟女老妇高清| 国产精品久久久久久久电影| 亚洲精品久久午夜乱码| 午夜视频国产福利| 中文精品一卡2卡3卡4更新| 在线观看人妻少妇| 人人妻人人澡人人看| 国产69精品久久久久777片| 男女啪啪激烈高潮av片| 国产精品熟女久久久久浪| 女人久久www免费人成看片| 制服人妻中文乱码| 一区二区日韩欧美中文字幕 | freevideosex欧美| 欧美日韩视频高清一区二区三区二| 国产乱人偷精品视频| 欧美亚洲 丝袜 人妻 在线| 日韩 亚洲 欧美在线| 边亲边吃奶的免费视频| 插阴视频在线观看视频| 亚洲av二区三区四区| 亚洲,一卡二卡三卡| 亚洲无线观看免费| 亚洲情色 制服丝袜| 日韩一区二区视频免费看| av又黄又爽大尺度在线免费看| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人精品一二三区| 丝瓜视频免费看黄片| 免费播放大片免费观看视频在线观看| 美女视频免费永久观看网站| 中文欧美无线码| a级毛色黄片| 亚洲精品亚洲一区二区| 黄色视频在线播放观看不卡| 永久免费av网站大全| 美女内射精品一级片tv| 国产精品一区二区三区四区免费观看| 精品少妇久久久久久888优播| 亚洲,欧美,日韩| 免费大片18禁| 精品人妻熟女毛片av久久网站| 一级二级三级毛片免费看| 亚洲精品久久成人aⅴ小说 | 成年人午夜在线观看视频| 男的添女的下面高潮视频| 老司机亚洲免费影院| 三级国产精品欧美在线观看| 亚洲国产最新在线播放| 老司机影院毛片| 最后的刺客免费高清国语| 插逼视频在线观看| 精品久久久精品久久久| 欧美日韩综合久久久久久| 美女内射精品一级片tv| 久久久国产精品麻豆| 一二三四中文在线观看免费高清| 国产欧美亚洲国产| 日日摸夜夜添夜夜爱| www.av在线官网国产| 啦啦啦中文免费视频观看日本| 精品少妇黑人巨大在线播放| 成人18禁高潮啪啪吃奶动态图 | 日产精品乱码卡一卡2卡三| 韩国高清视频一区二区三区| 亚洲精品日韩在线中文字幕| 少妇 在线观看| 亚洲av中文av极速乱| 一边摸一边做爽爽视频免费| 18在线观看网站| 午夜激情福利司机影院| 久久久亚洲精品成人影院| 少妇熟女欧美另类| 在线观看免费日韩欧美大片 | 欧美亚洲 丝袜 人妻 在线| 啦啦啦啦在线视频资源| 午夜av观看不卡| 青春草国产在线视频| 男女啪啪激烈高潮av片| 一级黄片播放器| 99精国产麻豆久久婷婷| 高清不卡的av网站| 国产白丝娇喘喷水9色精品| 国模一区二区三区四区视频| 国产av精品麻豆| videossex国产| 久久人妻熟女aⅴ| 丝袜脚勾引网站| 日产精品乱码卡一卡2卡三| a级毛色黄片| 五月伊人婷婷丁香| 高清毛片免费看| 99久久中文字幕三级久久日本| 黄色欧美视频在线观看| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 一级二级三级毛片免费看| 国产欧美亚洲国产| 热re99久久国产66热| 成人亚洲精品一区在线观看| 18禁在线无遮挡免费观看视频| 在线天堂最新版资源| av不卡在线播放| 久久免费观看电影| 久久久精品区二区三区| 日韩伦理黄色片| 国产亚洲av片在线观看秒播厂| 午夜激情久久久久久久| 国产午夜精品一二区理论片| 日韩成人av中文字幕在线观看| 欧美激情极品国产一区二区三区 | 热re99久久国产66热| 99国产综合亚洲精品| 国产片特级美女逼逼视频| 蜜臀久久99精品久久宅男| 晚上一个人看的免费电影| 亚洲国产精品专区欧美| 欧美97在线视频| 制服人妻中文乱码| 大香蕉久久网| 性色avwww在线观看| 亚洲国产精品一区二区三区在线| 国产精品99久久99久久久不卡 | 男女高潮啪啪啪动态图| 中国美白少妇内射xxxbb| 亚洲,一卡二卡三卡| 一级毛片黄色毛片免费观看视频| 简卡轻食公司| 如日韩欧美国产精品一区二区三区 | 51国产日韩欧美| 男人操女人黄网站| 伊人久久国产一区二区|