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

    連續(xù)式跨聲速風(fēng)洞大開角段整流裝置設(shè)計(jì)數(shù)值模擬

    2015-03-28 08:07:11李紅喆廖達(dá)雄叢成華
    關(guān)鍵詞:風(fēng)洞氣流阻尼

    李紅喆,廖達(dá)雄,叢成華

    (1.中國(guó)空氣動(dòng)力研究與發(fā)展中心設(shè)備設(shè)計(jì)及測(cè)試技術(shù)研究所,四川綿陽(yáng) 621000;2.中國(guó)空氣動(dòng)力研究與發(fā)展中心空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川綿陽(yáng) 621000)

    連續(xù)式跨聲速風(fēng)洞大開角段整流裝置設(shè)計(jì)數(shù)值模擬

    李紅喆1,*,廖達(dá)雄1,2,叢成華1,2

    (1.中國(guó)空氣動(dòng)力研究與發(fā)展中心設(shè)備設(shè)計(jì)及測(cè)試技術(shù)研究所,四川綿陽(yáng) 621000;2.中國(guó)空氣動(dòng)力研究與發(fā)展中心空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川綿陽(yáng) 621000)

    采用阻尼網(wǎng)對(duì)大開角段內(nèi)的氣流分離進(jìn)行控制,并合理設(shè)置其參數(shù),是工程上有效的方法之一。為驗(yàn)證阻尼網(wǎng)工程設(shè)計(jì)方法的可靠性,以0.6m連續(xù)式跨聲速風(fēng)洞為背景,通過數(shù)值模擬,對(duì)工程設(shè)計(jì)方法的初步結(jié)果進(jìn)行了驗(yàn)證,并在此基礎(chǔ)上結(jié)合大開角段布置環(huán)境對(duì)阻尼網(wǎng)參數(shù)進(jìn)行了優(yōu)化。由計(jì)算結(jié)果知,采用方案3-4(兩層阻尼網(wǎng)損失系數(shù)分別為1.6和1.0)時(shí),大開角段出口截面的速度均方根偏差值(RMS)為14.5%;考慮布置環(huán)境影響,調(diào)整兩層阻尼網(wǎng)損失系數(shù)至0.8和1.0時(shí),RMS值為16.2%。研究結(jié)果表明,阻尼網(wǎng)工程設(shè)計(jì)方法結(jié)合數(shù)值模擬可以有效地應(yīng)用于大開角段整流裝置的設(shè)計(jì),達(dá)到了抑制大開角段內(nèi)氣流分離,降低壓力損失,提高出口速度均勻性的設(shè)計(jì)目標(biāo)。

    連續(xù)式風(fēng)洞;跨聲速;大開角段;阻尼網(wǎng);數(shù)值模擬;流場(chǎng)品質(zhì)

    0 引 言

    我國(guó)國(guó)防現(xiàn)代化建設(shè)和航天航空事業(yè)的飛速發(fā)展,對(duì)風(fēng)洞試驗(yàn)設(shè)備提出了更高的要求??缏曀亠L(fēng)洞試驗(yàn)設(shè)備是高速空氣動(dòng)力學(xué)研究的主要手段,建設(shè)大型連續(xù)式跨聲速風(fēng)洞將大大提高我國(guó)跨聲速領(lǐng)域試驗(yàn)?zāi)M能力。連續(xù)式跨聲速風(fēng)洞通常采用較大的收縮比提高試驗(yàn)段流場(chǎng)品質(zhì),需要設(shè)計(jì)擴(kuò)散段實(shí)現(xiàn)較大收縮比并降低風(fēng)洞壓力損失。若采用常規(guī)擴(kuò)散段,其長(zhǎng)度將占到風(fēng)洞回路長(zhǎng)度的50%,而采用大開角段可以有效減少風(fēng)洞回路長(zhǎng)度,節(jié)約風(fēng)洞占地,減少造價(jià)及運(yùn)行成本。

    通常,擴(kuò)散段擴(kuò)散全角大于5°、面積比大于2就可能發(fā)生氣流分離[1],這將使出口氣流分布十分不均勻,造成較大的壓力損失,并嚴(yán)重影響試驗(yàn)段的氣流品質(zhì)。大開角段的擴(kuò)散角度大,氣流流動(dòng)需要克服較大的逆壓梯度,因此,壁面附近邊界層增長(zhǎng)迅速,容易出現(xiàn)邊界層分離。氣流分離將會(huì)產(chǎn)生較大分離壓損,造成出口氣流的不均勻分布,影響下游換熱器效率,甚至是試驗(yàn)段流場(chǎng)品質(zhì)。恰當(dāng)布置整流裝置,以減小氣流分離,提高出口氣流均勻性是大開角擴(kuò)散段(下文簡(jiǎn)稱為大開角段)設(shè)計(jì)的核心。

    阻尼網(wǎng)在連續(xù)式跨聲速風(fēng)洞大開角段整流中應(yīng)用最為廣泛且整流效果最好[2-3]。有較多學(xué)者對(duì)阻尼網(wǎng)設(shè)計(jì)方法進(jìn)行了研究,如Mehta、Farell、Smith等人,得出了一些結(jié)論,但并未形成完整的理論設(shè)計(jì)方法,大多仍需要依靠工程經(jīng)驗(yàn)完成設(shè)計(jì)[4-6]。在前期的研究中,基于零靜壓恢復(fù)思想總結(jié)形成了一套較為完整的阻尼網(wǎng)工程設(shè)計(jì)方法[7]。本文以0.6m連續(xù)式跨聲速風(fēng)洞大開角段為應(yīng)用背景,采用阻尼網(wǎng)工程設(shè)計(jì)方法設(shè)計(jì)大開角段內(nèi)阻尼網(wǎng)布置參數(shù),通過數(shù)值模擬驗(yàn)證并優(yōu)化阻尼網(wǎng)工程設(shè)計(jì)結(jié)果,同時(shí)研究前后部段對(duì)阻尼網(wǎng)設(shè)計(jì)的影響。

    1 計(jì)算模型及計(jì)算條件

    1.1 計(jì)算模型

    如圖1所示,大開角段采用圓變方過渡,入口尺寸Φ2 600mm,出口尺寸3 500mm×3 500mm,總長(zhǎng)2 500mm,當(dāng)量擴(kuò)散全角30.2°,面積比為2.31。選取大開角段入口截面中心點(diǎn)為坐標(biāo)原點(diǎn),沿流向?yàn)閤軸正方向。計(jì)算時(shí),入口前增加5 000mm的平直段,出口增加11 700mm的平直段。采用結(jié)構(gòu)網(wǎng)格,網(wǎng)格約140萬(wàn)。在研究大開角段前后部段的影響時(shí),計(jì)算模型包含了壓縮機(jī)尾罩、第二擴(kuò)散段、大開角段、換熱器段,網(wǎng)格約230萬(wàn)。阻尼網(wǎng)及換熱器區(qū)域網(wǎng)格進(jìn)行了加密處理,以更準(zhǔn)確地捕捉流場(chǎng)特征。

    圖1 計(jì)算模型氣動(dòng)布局示意圖(mm)Fig.1 Aerodynamic layout of simulation model(mm)

    大開角段內(nèi)阻尼網(wǎng)設(shè)計(jì)基于零壓力恢復(fù)思想,依據(jù)阻尼網(wǎng)工程設(shè)計(jì)方法進(jìn)行設(shè)計(jì)[7]。阻尼網(wǎng)網(wǎng)絲直徑為0.5mm,開孔率β由Wieghardt及Cornell公式計(jì)算,阻尼網(wǎng)層數(shù)、安裝位置、損失系數(shù)通過阻尼網(wǎng)工程設(shè)計(jì)方法確定[7-8]。

    1.2 計(jì)算條件

    風(fēng)洞采用干燥空氣作為試驗(yàn)介質(zhì),設(shè)計(jì)點(diǎn)試驗(yàn)段馬赫數(shù)Ma=0.9,穩(wěn)定段總壓p0=2.5×105Pa,當(dāng)?shù)卮髿鈮簆a=96 980Pa,大開角段入口截面氣流速度為13.2m/s,壓縮機(jī)尾罩段入口截面氣流速度為64.5m/s。

    計(jì)算采用不可壓N-S方程,分離隱式算法,邊界條件為速度入口和壓力出口,采用SSTk-ω湍流模型。使用多孔跳躍邊界條件模擬流體通過阻尼網(wǎng)的壓力損失[9]。換熱器段僅考慮壓損影響,不計(jì)熱傳導(dǎo),在換熱器段入口處將換熱器段等效為阻尼網(wǎng)進(jìn)行計(jì)算,換熱器的當(dāng)?shù)負(fù)p失系數(shù)為30。

    1.3 評(píng)價(jià)指標(biāo)

    進(jìn)行擴(kuò)散段性能分析時(shí),通常采用大開角段出口截面速度均方根偏差值RMS和大開角段總壓力恢復(fù)系數(shù)C*p衡量擴(kuò)散段性能優(yōu)劣[10-12]。取大開角段出口截面中心區(qū)域35×35個(gè)點(diǎn)計(jì)算氣流軸向速度的均方根偏差RMS值,中心區(qū)域占出口截面積的94.4%,排除了邊界層速度劇烈變化的影響。風(fēng)洞設(shè)計(jì)中,部段出口速度均方根偏差RMS一般小于15時(shí)即可視為均勻流場(chǎng)。采用大開角段入口截面動(dòng)壓計(jì)算總壓力恢復(fù)系數(shù)C*p。

    1.4 數(shù)值驗(yàn)證

    采用TWG風(fēng)洞縮比大開角段模型進(jìn)行了計(jì)算方法驗(yàn)證[13]。該模型入口截面直徑420mm,面積比為3.26,擴(kuò)散全角46.6°。圖2為該模型大開角段出口截面中線速度分布,本文計(jì)算結(jié)果與風(fēng)洞試驗(yàn)結(jié)果、文獻(xiàn)計(jì)算結(jié)果吻合程度較高。另外,還進(jìn)行了網(wǎng)格無(wú)關(guān)性驗(yàn)證,限于篇幅這里不再贅述。

    圖2 TWG風(fēng)洞縮比大開角段模型加網(wǎng)后計(jì)算及實(shí)驗(yàn)結(jié)果對(duì)比Fig.2 Comparison of computational result and experimental result about wide angle diffuser with screens of the scaled TWG wind tunnel

    2 阻尼網(wǎng)設(shè)計(jì)

    圖3為大開角段對(duì)角面速度分布,從入口開始出現(xiàn)氣流分離,分離區(qū)域較大,主要集中在大開角段截面由圓變方過渡的四個(gè)角區(qū)。出口截面軸向速度RMS值為111.9%,大開角段總壓力恢復(fù)系數(shù)C*p為0.1。由于大開角段為圓變方構(gòu)型,角區(qū)附近的邊界層為兩個(gè)壁面邊界層的疊加,因此角區(qū)附近的邊界層增長(zhǎng)快,較之圓形壁面更容易出現(xiàn)流動(dòng)分離。大開角段出口速度分布十分不均勻,需要設(shè)計(jì)阻尼網(wǎng)以消除內(nèi)部氣流分離,提高出口氣流均勻性。

    圖3 大開角段未布置阻尼網(wǎng)時(shí)對(duì)角面速度分布Fig.3 Velocity distribution of diagonal surface in wide angle diffuser without screens

    2.1 初步設(shè)計(jì)

    根據(jù)阻尼網(wǎng)工程設(shè)計(jì)方法可知,0.6m風(fēng)洞大開角段構(gòu)型布置兩層阻尼網(wǎng)較為合理,確定方案3-4(指第一層阻尼網(wǎng)安裝在與大開角段入口面積比為1.3處,第二層阻尼網(wǎng)安裝位置與第一層阻尼網(wǎng)安裝位置截面面積比為1.4)為初步設(shè)計(jì)結(jié)果,另選取了15種阻尼網(wǎng)布置方案作為設(shè)計(jì)對(duì)比,各方案參數(shù)見表1。Ks為阻尼網(wǎng)當(dāng)?shù)負(fù)p失系數(shù),β為開孔率,下標(biāo)1、2表示層數(shù)。

    圖4給出了16種方案的阻尼網(wǎng)設(shè)計(jì)參數(shù)。從出口截面軸向速度RMS看,加入兩層阻尼網(wǎng)后大開角段出口速度均勻性大幅提高,從第1組到第4組,氣流均勻性呈上升趨勢(shì),且每組中,出口速度RMS也逐漸降低。阻尼網(wǎng)位置后移,有助于提高出口速度的均勻程度,但尚未達(dá)到高品質(zhì)流場(chǎng)標(biāo)準(zhǔn)。由此可見,根據(jù)工程設(shè)計(jì)方法得到的阻尼網(wǎng)損失系數(shù)略小,應(yīng)適當(dāng)增加損失系數(shù)。從總壓力恢復(fù)系數(shù)C*p可以得到相似結(jié)論,大開角段較大的面積比導(dǎo)致大開角段前后流速變化大,當(dāng)?shù)貕毫p失相同的阻尼網(wǎng)后移壓損減小,同時(shí),隨著阻尼網(wǎng)位置后移,擴(kuò)散段內(nèi)因氣流分離所引起的靜壓損失減小,壓力恢復(fù)系數(shù)呈上升趨勢(shì)。

    圖5為第二層阻尼網(wǎng)與第一層位置面積比為1.4時(shí)對(duì)角面速度流線圖。加入阻尼網(wǎng)后,氣流分離得到了較好的抑制,分離區(qū)域明顯減小,僅在四個(gè)角區(qū)存在小范圍分離。由于阻尼網(wǎng)后,邊界層厚度減小,壁面附近靜壓降低,氣流流向壁面,使邊界層可以重新承受較大的逆壓梯度,有效防止氣流發(fā)生分離。計(jì)算結(jié)果顯示,在保持間距不變的情況下,阻尼網(wǎng)位置后移則壓力損失減小,出口速度均勻性提高。但從局部看,網(wǎng)前分離區(qū)域擴(kuò)大,出口分離區(qū)域減小,為此必須進(jìn)一步改進(jìn)阻尼網(wǎng)設(shè)計(jì)以消除分離。

    表1 16組阻尼網(wǎng)初步設(shè)計(jì)參數(shù)Table 1 Preliminary design coefficients of 16kinds of screen arrangements

    圖4 16組阻尼網(wǎng)計(jì)算結(jié)果Fig.4 Simulation results of 16kinds of screen arrangements

    圖5 布置阻尼網(wǎng)后大開角段對(duì)角面速度分布Fig.5 Velocity distribution of diagonal surface in wide angle diffuser with screens

    2.2 優(yōu)化設(shè)計(jì)

    增大第一層阻尼網(wǎng)損失系數(shù)可以有效減弱網(wǎng)前氣流分離,對(duì)抑制網(wǎng)間、網(wǎng)后分離也十分重要。但大開角段截面面積變化引起出、入口較大的速度變化,增加第二層阻尼網(wǎng)損失帶來的壓力損失相對(duì)較小,且對(duì)提高擴(kuò)散段出口速度均勻性有重要作用。因此,需合理分配兩層阻尼網(wǎng)損失系數(shù)。

    為消除網(wǎng)前分離,在第二層網(wǎng)損失系數(shù)不變的情況下增加第一層阻尼網(wǎng)損失系數(shù),分別將第一層阻尼網(wǎng)的損失系數(shù)增加至0.8、1.2、1.6、2.0、3.0、4.0(分別對(duì)應(yīng)優(yōu)化方案1~6)。

    從表2可以看到,隨著第一層阻尼網(wǎng)損失系數(shù)的增加,出口速度分布均勻性提高,但壓力恢復(fù)系數(shù)有所下降。

    表2 四種優(yōu)化方案的RMS及值Table 2RMSandvalue of 4optimized schemes

    表2 四種優(yōu)化方案的RMS及值Table 2RMSandvalue of 4optimized schemes

    方案編號(hào)優(yōu)化方案1優(yōu)化方案2RMS/%C*pRMS/%C*p1-4 24.3-0.17 21.9-0.50 2-4 21.3-0.03 18.9-0.31 3-4 19.7 0.09 17.3-0.14 4-4 19.5 0.16 16.7 0.07方案編號(hào)優(yōu)化方案3優(yōu)化方案4RMS/%C*pRMS/%C*p1-4 20.8-0.83 20.0-1.16 2-4 17.5-0.59 16.6-0.75 3-4 15.9-0.37 15.1-0.60 4-4 15.4-0.23 14.6-0.33方案編號(hào)優(yōu)化方案5優(yōu)化方案6RMS/%C*pRMS/%C*p1-4 15.6-1.95 14.4-2.77 2-4 11.0-1.54 9.5-2.23 3-4 9.2-1.17 7.9-1.74 4-4 8.7-0.93 7.8-1.43

    圖6為優(yōu)化阻尼網(wǎng)布置后大開角段內(nèi)最大分離速度,可以看到,增加第一層網(wǎng)損失系數(shù),分離速度減小。對(duì)比四種阻尼網(wǎng)組合,1-4、2-4逆向速度較小,3-4有所增加,4-4大幅增加。初始方案中,分離在-3m/s左右。增加第一層網(wǎng)損失系數(shù)后,優(yōu)化方案1、2的分離區(qū)域和分離速度顯著減小。繼續(xù)增加第一層網(wǎng)損失系數(shù)(優(yōu)化方案3、4、5、6)對(duì)分離區(qū)域和分離速度的改變很小,但壓力恢復(fù)系數(shù)急劇降低。

    阻尼網(wǎng)前移有助于減小網(wǎng)前分離,但壓力損失急劇增加,阻尼網(wǎng)后移則增加網(wǎng)前分離,并導(dǎo)致最大逆向速度跳躍增加。為此,確定優(yōu)化方案3、4為合理設(shè)計(jì)方案,能有效抑制網(wǎng)前分離。2-4、3-4方案為較為優(yōu)化的阻尼網(wǎng)組合,阻尼網(wǎng)位置如圖7所示。

    圖6 六組優(yōu)化方案的大開角段內(nèi)最大分離速度Fig.6 Max adverse velocity of 6optimized schemes in wide angle diffuser

    圖7 阻尼網(wǎng)布置位置示意圖Fig.7 Schematic of screens location

    基于方案3增加第二層阻尼網(wǎng)的損失系數(shù)形成方案7,Ks1和Ks2分別為1.6和1.0,基于方案4形成方案8,Ks1和Ks2分別為2.0和1.0。

    從圖8可以看出,增加第二層網(wǎng)損失系數(shù)后RMS有所下降,但壓力恢復(fù)系數(shù)降低。方案7與4相比,適當(dāng)減小第一層網(wǎng)損失系數(shù),同時(shí)增加第二層網(wǎng)損失系數(shù)后,在降低擴(kuò)散段出口截面RMS的同時(shí),提高了擴(kuò)散段的壓力恢復(fù)系數(shù),減小了擴(kuò)散段損失。

    圖8 四組優(yōu)化方案計(jì)算結(jié)果Fig.8 Simulation results of 4kinds of optimized schemes

    由此,方案7的3-4組合為針對(duì)本文所研究大開角段模型的合理優(yōu)化阻尼網(wǎng)布置方式。若需要進(jìn)一步減小網(wǎng)前分離,可以繼續(xù)增加第一層網(wǎng)的損失系數(shù),但改善幅度很小,會(huì)帶來較大的壓損。如需獲得RMS小于15%的出口速度分布,可通過增加第二層網(wǎng)損失系數(shù)的方式實(shí)現(xiàn)。

    3 一體化模型阻尼網(wǎng)設(shè)計(jì)

    大開角段內(nèi)部流場(chǎng)直接受到前后部段的影響,考慮前后部段進(jìn)行模擬,能夠獲得更真實(shí)的流場(chǎng)分布,更有針對(duì)性地對(duì)其內(nèi)部的整流裝置進(jìn)行設(shè)計(jì)。

    圖9給出了布置壓縮機(jī)尾罩和第二擴(kuò)散段后空大開角段對(duì)角面速度分布,從大開角段入口開始出現(xiàn)氣流分離,分離區(qū)域主要集中在四個(gè)角區(qū),并逐漸向后發(fā)展。與圖3比較可以看出,布置壓縮機(jī)尾罩和第二擴(kuò)散段后,分離區(qū)域明顯向下游擴(kuò)展。由于尾罩的存在,壓縮機(jī)尾罩段呈擴(kuò)散狀,和第二擴(kuò)散段均會(huì)產(chǎn)生逆壓梯度,使邊界層變厚,導(dǎo)致氣流分離,對(duì)流場(chǎng)品質(zhì)產(chǎn)生不良影響。大開角段出口截面RMS值為122.4,壓力恢復(fù)系數(shù)為-0.01,兩個(gè)參數(shù)均有相當(dāng)程度的降低。

    圖9 布置壓縮機(jī)尾罩及第二擴(kuò)散段后模型對(duì)角面速度分布Fig.9 Velocity distribution of diagonal surface in wide angle diffuser with compressor tail cone and 2nd diffuser without screens

    從圖10可以看出,布置換熱器能有效減小分離尺度,限制分離區(qū)域向后發(fā)展,使分離集中在大開角段內(nèi)。大開角段出口截面RMS為76.2%,擴(kuò)散段總壓力恢復(fù)系數(shù)C*p為0.33。與未加換熱器時(shí)相比,加入換熱器有利于均勻出口氣流和壓力恢復(fù)。也要看到,由于換熱器阻力過大,導(dǎo)致氣流在靠近壁面附近出現(xiàn)了過沖現(xiàn)象[14]。由于大開角段內(nèi)壁面附近邊界層剝離嚴(yán)重,產(chǎn)生較大的分離損失,因此分離后,壁面附近的氣流靜壓較中心區(qū)域低。換熱器損失系數(shù)較大時(shí),為保持通過流量不變,氣流在換熱器后的壁面附近出現(xiàn)了速度過沖。損失系數(shù)越大,則過沖現(xiàn)象越顯著。在阻尼網(wǎng)設(shè)計(jì)時(shí),如果阻力系數(shù)過大,也會(huì)出現(xiàn)這一問題。

    圖11為一體化模型對(duì)角面速度分布。一體化模型大開角段出口截面RMS值為81.5%,C*p為0.3,分離區(qū)域集中出現(xiàn)在大開角段內(nèi)。

    對(duì)一體化模型的阻尼網(wǎng)設(shè)計(jì),采用方案2-4、3-4及優(yōu)化方案7。從表3可以看出兩種方案均可有效抑制大開角段內(nèi)的大范圍氣流分離,獲得均勻程度高的出口截面速度分布。

    圖10 布置換熱器后對(duì)角面速度分布Fig.10 Velocity distribution of diagonal surface in wide angle diffuser with heat exchanger

    圖11 一體化模型對(duì)角面速度分布Fig.11 Velocity distribution of diagonal surface in integrative model

    表3 一體化模型布置阻尼網(wǎng)后計(jì)算結(jié)果Table 3 Results of the integrative model with screens

    為抑制氣流分離,盡量降低壓力損失,確定3-4方案為最終方案??紤]前后部段的影響,將兩層阻尼網(wǎng)損失系數(shù)減小為0.8、1.0后出口RMS為16.2%,有效提高了出口速度均勻性,基本滿足設(shè)計(jì)要求;總壓力恢復(fù)系數(shù)C*p為-0.03,有效降低了壓力損失,與零壓力恢復(fù)設(shè)計(jì)思想相吻合。圖12為采用最終方案計(jì)算得到的一體化模型流場(chǎng)流線顯示圖,氣流分布十分均勻,僅在大開角段入口至第一層網(wǎng)前四個(gè)角區(qū)存在微小的回流區(qū),達(dá)到了抑制大開角段內(nèi)大范圍氣流分離的目的。

    圖12 方案3-4一體化模型流線圖Fig.12 Streamlines of the integrative model with screens scheme 3-4

    4 結(jié) 論

    通過數(shù)值模擬對(duì)0.6m連續(xù)式跨聲速風(fēng)洞大開角段內(nèi)阻尼網(wǎng)的工程設(shè)計(jì)結(jié)果進(jìn)行了優(yōu)化,結(jié)論如下:

    (1)阻尼網(wǎng)工程設(shè)計(jì)方法是可行的。其設(shè)計(jì)基礎(chǔ)是前無(wú)其它過渡部段的圓形截面擴(kuò)散段,在實(shí)際設(shè)計(jì)時(shí),應(yīng)根據(jù)實(shí)際的大開角段構(gòu)型,對(duì)阻尼網(wǎng)的設(shè)計(jì)損失系數(shù)進(jìn)行調(diào)整。對(duì)于0.6m風(fēng)洞,3-4組合是較為合理的方案,兩層阻尼網(wǎng)損失系數(shù)為1.6和1.0時(shí)可有效抑制氣流分離。

    (2)阻尼網(wǎng)設(shè)計(jì)應(yīng)充分考慮大開角段布置環(huán)境對(duì)其性能的影響。換熱器段有利于抑制大開角段內(nèi)氣流分離沿軸向發(fā)展,第二擴(kuò)散段和壓縮機(jī)尾罩段則加劇了分離的產(chǎn)生。根據(jù)0.6m風(fēng)洞布局,確定了方案3-4,兩層阻尼網(wǎng)損失系數(shù)分別為0.8、1.0。該方案實(shí)現(xiàn)了高性能大開角段無(wú)明顯氣流分離,降低壓力損失,提高出口速度均勻性的設(shè)計(jì)目標(biāo)。

    研究中發(fā)現(xiàn),使用損失系數(shù)較大的阻尼網(wǎng)時(shí),氣流在洞壁附近出現(xiàn)了速度過沖現(xiàn)象,未來工作中將考慮洞壁邊界層效應(yīng),研究變目數(shù)阻尼網(wǎng)對(duì)整流效果的影響。下一步還將在0.6m連續(xù)式跨聲速風(fēng)洞形成試驗(yàn)?zāi)芰r(shí)獲得大開角擴(kuò)散段試驗(yàn)數(shù)據(jù),與計(jì)算結(jié)果進(jìn)行對(duì)比,完善設(shè)計(jì)方法。

    [1] Liu Zhengchong.Aerodynamic and structure design of high and low speed wind tunnel[M].Beijing:National defence industry press,2003.(in Chinese)劉政崇.高低速風(fēng)洞氣動(dòng)與結(jié)構(gòu)設(shè)計(jì)[M].北京:國(guó)防工業(yè)出版社,2003.

    [2] Mehta R D.The aerodynamic design of blower tunnels with wideangle diffusers[J].Progress in Aerospace Sciences,1979,18: 59-120.doi:10.1016/0376-0421(77)90003-3

    [3] Mehta R D,Bradshaw P.Design rules for small low speed wind tunnels[J].Aeronaut.Journal,1979,83.

    [4] Mehta R D.The design of wide-angle diffusers[R].I.C.Aero Report 76-03,1976.doi:10.2514/3.44542

    [5] Farell C,Xia L.A note on the design of screen-filled wide-angle diffusers[J].Journal of Wind Engineering and Industrial Aerodynamics,1990,33(3):479-486.doi:10.1016/0167-6105(90)90002-T

    [6] Ward-Smith A J,Lane D,Reynolds A J,et al.Flow regimes in wide-angle screened diffusers[J].International Journal of Mechanical Sciences,1991,33(1):41-54.doi:10.1016/0020-7403(91)90026-Y

    [7] Li Hongzhe,Liao Daxiong,Cong Chenghua.A review of design methods of continuous transonic wind tunnel’s wide-angle diffuser[J].Journal of Experiments in Fluid Mechanics,2012,26(6):59-65.(in Chinese)李紅喆,廖達(dá)雄,叢成華.連續(xù)式跨聲速風(fēng)洞大開角擴(kuò)散段設(shè)計(jì)方法研究概述[J].實(shí)驗(yàn)流體力學(xué),2012,26(6):59-65.doi:10.3969/j.issn.1672-9897.2012.06.013

    [8] Seltsam M M.Experimental and theoretical study of wide-angle diffuser flow with screens[J].AIAA Journal,1995,33(11):2092-2100.doi:10.2514/3.12952

    [9] Fluent Inc.FLUENT user′s guide,2003.

    [10]Sahin B,Ward-Smith A J,Lane D.The pressure drop and flow characteristics of wide-angle screened diffusers of large area ratio[J].Journal of Wind Engineering and Industrial Aerodynamics,1995,58(1-2):33-50.doi:10.1016/j.expthermflusci.2004.10.002.

    [11]Sahin B,Ward-Smith A J.The pressure distribution in and flow characteristics of wide angle diffusers using perforated plates for the flow control with application to electrostatic precipitators[J].Int.J.of Mech.Engg.Sci.,1993,35:117-127.doi:10.1016/0020-7403(93)90070-B.

    [12]Noui-Mehid M N,Wu J,Sutalo I D,et al.Velocity distribution downstream of an asymmetric wide-angle diffuser[J].Experimental Thermal and Fluid Science,2005,29(6):649-657.doi:10.1016/j.expthermflusci.2004.10.002

    [13]Seltsam M M.A new tool for screen integration in internal flow systems and its application to wind tunnel design[C].Aerodynamics of Wind Tunnel Circuits and their Components,F(xiàn)rance,1997:1-12.

    [14]Mehta R D.Turbulent boundary layer perturbed by a screen[J].AIAA Journal,1985,23(9):1335-1342.doi:10.2514/3.9089

    Numerical simulation of flow conditioning device design in wide angle diffuser of continuous transonic wind tunnels

    Li Hongzhe1,Liao Daxiong1,2,Cong Chenghua1,2
    (1.FacilityDesignandInstrumentationInstituteofChinaAerodynamicsResearchandDevelopmentCenter,Mianyang621000,China2.StateKeyLaboratoryofAerodynamics,ChinaAerodynamicsResearchandDevelopmentCenter,Mianyang621000,China)

    Large area expansion can be achieved effectively by wide angle diffusers in wind tunnels,while the airflow separation may occur at the same time.In order to ensure the flow field quality,the separation must be suppressed.Utilizing damp-screens with reasonable parameters is a valid method in engineering.With the background of 0.6mcontinuum transonic wind tunnel,so as to test the engineering design method of damp-screens,we verify the preliminary result by the method through the way of numerical simulation.From the calculation,by taking the scheme 3-4,loss coefficients for two screens are 1.6and 1.0respectively,theRMSvalue at the outlet of the wide angle diffuser will be 14.5%;considering the influence by the environment and adjusting the two coefficients to 0.8and 1.0,theRMSvalue will be 16.2%.The research indicates that,the engineering design method,in addition to numerical simulations,can be effectively applied in the design of wide angle diffusers.Suppressing the airflow separation,reducing the pressure loss,and improving the exit velocity uniformity can be achieved simultaneously.

    continuous wind tunnel;transonic;wide angle diffuser;screen;numerical simulation;flow quality

    V211.74

    :Adoi:10.7638/kqdlxxb-2013.0078

    0258-1825(2015)02-0198-06

    2013-07-20;

    :2013-10-30

    李紅喆*(1987-),女,陜西寶雞人,碩士,工程師,從事風(fēng)洞氣動(dòng)設(shè)計(jì)及數(shù)值模擬研究.E-mail:cardc.lhz@gmail.com

    李紅喆,廖達(dá)雄,叢成華.連續(xù)式跨聲速風(fēng)洞大開角段整流裝置設(shè)計(jì)數(shù)值模擬[J].空氣動(dòng)力學(xué)學(xué)報(bào),2015,33(2):198-203.

    10.7638/kqdlxxb-2013.0078 Li H Z,Liao D Xi,Cong C H.Numerical simulation of flow conditioning device design in wide angle diffuser of continuous transonic wind tunnels[J].Acta Aerodynamica Sinica,2015,33(2):198203.

    猜你喜歡
    風(fēng)洞氣流阻尼
    氣流的威力
    N維不可壓無(wú)阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    斑頭雁進(jìn)風(fēng)洞
    黃風(fēng)洞貂鼠精
    基于NI cRIO平臺(tái)的脈沖燃燒風(fēng)洞控制系統(tǒng)設(shè)計(jì)
    固體運(yùn)載火箭變軌發(fā)動(dòng)機(jī)噴管氣流分離研究
    具阻尼項(xiàng)的Boussinesq型方程的長(zhǎng)時(shí)間行為
    飛片下的空氣形成的“超強(qiáng)高速氣流刀”
    国产精品一国产av| 亚洲专区中文字幕在线| 黄色视频不卡| 亚洲av成人精品一二三区| 制服人妻中文乱码| 日韩 亚洲 欧美在线| 国产野战对白在线观看| 精品免费久久久久久久清纯 | 九色亚洲精品在线播放| a级毛片在线看网站| 国产精品人妻久久久影院| 热re99久久国产66热| 黄色 视频免费看| 亚洲av美国av| 69精品国产乱码久久久| 欧美大码av| a 毛片基地| 国产精品久久久av美女十八| 高清黄色对白视频在线免费看| 久久久久久亚洲精品国产蜜桃av| 2021少妇久久久久久久久久久| 乱人伦中国视频| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 在线观看国产h片| 水蜜桃什么品种好| 国产午夜精品一二区理论片| 亚洲精品久久久久久婷婷小说| 久久久国产一区二区| 一本综合久久免费| 精品一区二区三卡| 电影成人av| 国产一区二区在线观看av| 中文字幕av电影在线播放| 蜜桃国产av成人99| 看十八女毛片水多多多| 一级毛片 在线播放| netflix在线观看网站| 嫩草影视91久久| 啦啦啦 在线观看视频| 后天国语完整版免费观看| 在线观看国产h片| 国产日韩欧美视频二区| 亚洲成av片中文字幕在线观看| 嫁个100分男人电影在线观看 | 亚洲中文av在线| 观看av在线不卡| 一级毛片女人18水好多 | 一级毛片 在线播放| 91九色精品人成在线观看| 韩国精品一区二区三区| 国产伦人伦偷精品视频| 多毛熟女@视频| 亚洲精品一卡2卡三卡4卡5卡 | 在线天堂中文资源库| 亚洲天堂av无毛| xxxhd国产人妻xxx| 99久久人妻综合| 久久精品人人爽人人爽视色| 少妇人妻久久综合中文| 老司机亚洲免费影院| 少妇精品久久久久久久| 免费观看av网站的网址| svipshipincom国产片| 欧美日韩精品网址| 亚洲成av片中文字幕在线观看| 如日韩欧美国产精品一区二区三区| 亚洲国产精品一区三区| 亚洲精品中文字幕在线视频| 午夜福利在线免费观看网站| xxx大片免费视频| 2021少妇久久久久久久久久久| 在线观看一区二区三区激情| 在线看a的网站| 熟女少妇亚洲综合色aaa.| 国产片特级美女逼逼视频| 青春草视频在线免费观看| 国产成人精品久久二区二区91| 美女扒开内裤让男人捅视频| 另类精品久久| 黄色一级大片看看| 可以免费在线观看a视频的电影网站| 日本欧美国产在线视频| 国产精品久久久人人做人人爽| 国产女主播在线喷水免费视频网站| av网站在线播放免费| a级毛片黄视频| 日韩欧美一区视频在线观看| 国产精品一国产av| av国产久精品久网站免费入址| 丁香六月天网| 1024视频免费在线观看| 免费一级毛片在线播放高清视频 | 晚上一个人看的免费电影| 女人久久www免费人成看片| 18禁黄网站禁片午夜丰满| 啦啦啦视频在线资源免费观看| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品无人区| 免费久久久久久久精品成人欧美视频| 又紧又爽又黄一区二区| 国产精品九九99| 久久这里只有精品19| 男男h啪啪无遮挡| 精品亚洲乱码少妇综合久久| 天天躁夜夜躁狠狠久久av| 欧美大码av| 欧美 日韩 精品 国产| 久久精品人人爽人人爽视色| 亚洲欧洲日产国产| 午夜视频精品福利| 天堂中文最新版在线下载| 午夜激情久久久久久久| 免费看av在线观看网站| 国产免费现黄频在线看| www.精华液| 亚洲欧美日韩另类电影网站| 精品少妇黑人巨大在线播放| 首页视频小说图片口味搜索 | 欧美日韩福利视频一区二区| av国产久精品久网站免费入址| 在线观看一区二区三区激情| 亚洲精品久久久久久婷婷小说| xxxhd国产人妻xxx| 国产又爽黄色视频| 国产一区二区三区av在线| 亚洲,一卡二卡三卡| 人人妻,人人澡人人爽秒播 | 一本综合久久免费| 久久av网站| 国产欧美日韩综合在线一区二区| 欧美国产精品va在线观看不卡| 久久精品熟女亚洲av麻豆精品| 亚洲人成电影观看| 一级毛片我不卡| 美女中出高潮动态图| 国产黄色免费在线视频| 亚洲国产欧美网| 亚洲第一av免费看| 久9热在线精品视频| av国产久精品久网站免费入址| 国产日韩欧美在线精品| 波多野结衣一区麻豆| a 毛片基地| 久久毛片免费看一区二区三区| 久久九九热精品免费| 少妇猛男粗大的猛烈进出视频| 涩涩av久久男人的天堂| 国产色视频综合| 中文字幕最新亚洲高清| 晚上一个人看的免费电影| 香蕉丝袜av| 多毛熟女@视频| 黄色 视频免费看| 欧美日韩亚洲综合一区二区三区_| 中文精品一卡2卡3卡4更新| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜一区二区| 精品国产一区二区久久| 日韩 亚洲 欧美在线| 捣出白浆h1v1| 成人三级做爰电影| 亚洲三区欧美一区| 精品国产乱码久久久久久男人| 男女高潮啪啪啪动态图| 老熟女久久久| 1024视频免费在线观看| 欧美乱码精品一区二区三区| 亚洲视频免费观看视频| av片东京热男人的天堂| 97精品久久久久久久久久精品| 黄色视频不卡| 男女高潮啪啪啪动态图| 99久久人妻综合| 母亲3免费完整高清在线观看| 久久人人爽av亚洲精品天堂| 久久国产精品男人的天堂亚洲| 汤姆久久久久久久影院中文字幕| 又黄又粗又硬又大视频| 亚洲精品美女久久av网站| 首页视频小说图片口味搜索 | 国产麻豆69| 美女国产高潮福利片在线看| 看免费成人av毛片| 欧美日韩亚洲国产一区二区在线观看 | 91字幕亚洲| 亚洲av日韩在线播放| 午夜免费鲁丝| 波多野结衣av一区二区av| 亚洲精品日韩在线中文字幕| 少妇猛男粗大的猛烈进出视频| 亚洲国产看品久久| 少妇人妻 视频| 下体分泌物呈黄色| 免费高清在线观看日韩| 久久久久精品国产欧美久久久 | 男人添女人高潮全过程视频| 少妇人妻 视频| 欧美精品亚洲一区二区| 色网站视频免费| 久久ye,这里只有精品| 久久人人97超碰香蕉20202| 国产女主播在线喷水免费视频网站| 校园人妻丝袜中文字幕| 一区二区日韩欧美中文字幕| 一本—道久久a久久精品蜜桃钙片| 一边摸一边做爽爽视频免费| 99国产综合亚洲精品| 国产人伦9x9x在线观看| 波多野结衣一区麻豆| 超碰成人久久| 性高湖久久久久久久久免费观看| 国产日韩一区二区三区精品不卡| 国产又爽黄色视频| 中文字幕人妻熟女乱码| 久热爱精品视频在线9| 两个人免费观看高清视频| 婷婷色麻豆天堂久久| 国产片内射在线| 欧美+亚洲+日韩+国产| 一级片免费观看大全| 中国国产av一级| 新久久久久国产一级毛片| 免费高清在线观看视频在线观看| 一本大道久久a久久精品| 国产精品二区激情视频| 国产精品 欧美亚洲| 韩国精品一区二区三区| 免费观看a级毛片全部| 老司机亚洲免费影院| 黄色视频不卡| 国产精品一区二区免费欧美 | 国产福利在线免费观看视频| 亚洲国产最新在线播放| 亚洲,欧美精品.| 狂野欧美激情性xxxx| 老司机靠b影院| 电影成人av| 9191精品国产免费久久| 在线天堂中文资源库| 国产精品欧美亚洲77777| 黑人巨大精品欧美一区二区蜜桃| 免费观看a级毛片全部| a 毛片基地| 亚洲人成77777在线视频| 中国国产av一级| 国产一区有黄有色的免费视频| 麻豆国产av国片精品| 一边摸一边抽搐一进一出视频| 五月开心婷婷网| 男女国产视频网站| 亚洲精品久久成人aⅴ小说| 一区二区av电影网| 国产日韩欧美视频二区| 日韩av不卡免费在线播放| 一级毛片我不卡| 捣出白浆h1v1| 精品一区在线观看国产| 青春草视频在线免费观看| 亚洲欧洲国产日韩| 亚洲欧美激情在线| 亚洲av欧美aⅴ国产| 亚洲av片天天在线观看| 成人免费观看视频高清| 亚洲欧美日韩高清在线视频 | 久久人人97超碰香蕉20202| 51午夜福利影视在线观看| 91麻豆av在线| 国产亚洲一区二区精品| 考比视频在线观看| 久久99精品国语久久久| 啦啦啦在线免费观看视频4| 日本91视频免费播放| 如日韩欧美国产精品一区二区三区| 亚洲国产精品成人久久小说| 天堂中文最新版在线下载| 99国产精品99久久久久| 精品人妻熟女毛片av久久网站| 美女扒开内裤让男人捅视频| av在线app专区| 香蕉国产在线看| a 毛片基地| 狠狠婷婷综合久久久久久88av| 国精品久久久久久国模美| 欧美老熟妇乱子伦牲交| 久久久国产一区二区| av在线app专区| 别揉我奶头~嗯~啊~动态视频 | 熟女av电影| 91麻豆av在线| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 性色av一级| 日韩av不卡免费在线播放| 在现免费观看毛片| 久久精品国产亚洲av高清一级| 99久久综合免费| 亚洲久久久国产精品| 一区在线观看完整版| 女性被躁到高潮视频| 人体艺术视频欧美日本| 精品久久久精品久久久| 黑丝袜美女国产一区| 国精品久久久久久国模美| 亚洲第一青青草原| 久久热在线av| 高清视频免费观看一区二区| 久久国产精品人妻蜜桃| kizo精华| 在线观看国产h片| 午夜免费观看性视频| 国产精品人妻久久久影院| 亚洲精品一区蜜桃| 色婷婷av一区二区三区视频| 亚洲国产精品国产精品| 99热网站在线观看| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 热re99久久精品国产66热6| 亚洲国产欧美一区二区综合| 国产片特级美女逼逼视频| 免费观看av网站的网址| 国产精品熟女久久久久浪| 久久人妻福利社区极品人妻图片 | cao死你这个sao货| 精品卡一卡二卡四卡免费| 欧美激情极品国产一区二区三区| 日日夜夜操网爽| 高清视频免费观看一区二区| 青青草视频在线视频观看| kizo精华| 久久人人爽av亚洲精品天堂| 十八禁网站网址无遮挡| 亚洲 欧美一区二区三区| 欧美人与善性xxx| 久久影院123| 精品人妻在线不人妻| 777米奇影视久久| tube8黄色片| 久久久久久亚洲精品国产蜜桃av| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久亚洲精品国产蜜桃av| 亚洲国产毛片av蜜桃av| 中文字幕色久视频| 国产又爽黄色视频| www.精华液| 久久精品熟女亚洲av麻豆精品| 91麻豆av在线| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美精品永久| 久久精品aⅴ一区二区三区四区| 啦啦啦在线观看免费高清www| 在线观看人妻少妇| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人免费av在线播放| 精品一区二区三区四区五区乱码 | 成人手机av| 秋霞在线观看毛片| 只有这里有精品99| 蜜桃国产av成人99| 久久精品久久久久久噜噜老黄| 中文字幕亚洲精品专区| 日韩制服骚丝袜av| 97在线人人人人妻| 精品久久蜜臀av无| 日韩欧美一区视频在线观看| 欧美大码av| 亚洲,欧美精品.| 9色porny在线观看| 19禁男女啪啪无遮挡网站| 日日夜夜操网爽| 亚洲中文av在线| 美女福利国产在线| 欧美日韩综合久久久久久| 色播在线永久视频| 国产伦人伦偷精品视频| 老司机在亚洲福利影院| 亚洲人成电影观看| 亚洲国产欧美日韩在线播放| 国产国语露脸激情在线看| 自线自在国产av| 精品国产一区二区久久| 一个人免费看片子| 国产精品免费大片| 亚洲成人国产一区在线观看 | 美女脱内裤让男人舔精品视频| 一区二区三区乱码不卡18| 欧美激情高清一区二区三区| 在线 av 中文字幕| 亚洲精品日韩在线中文字幕| 国产无遮挡羞羞视频在线观看| 国产精品一国产av| 韩国精品一区二区三区| 一本色道久久久久久精品综合| 久久久国产欧美日韩av| 色综合欧美亚洲国产小说| 男人添女人高潮全过程视频| 欧美另类一区| 99国产精品免费福利视频| av在线播放精品| 看十八女毛片水多多多| 欧美人与性动交α欧美软件| 国产精品 国内视频| 中文乱码字字幕精品一区二区三区| 欧美精品人与动牲交sv欧美| 久久久精品区二区三区| 日本一区二区免费在线视频| 国产亚洲精品第一综合不卡| 欧美日韩亚洲国产一区二区在线观看 | 精品福利永久在线观看| 久久国产精品影院| 日本午夜av视频| 999久久久国产精品视频| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av片在线观看秒播厂| 777米奇影视久久| 久久久精品免费免费高清| 国产视频首页在线观看| 亚洲精品美女久久av网站| 91精品伊人久久大香线蕉| 夫妻性生交免费视频一级片| 色婷婷久久久亚洲欧美| 麻豆av在线久日| 性色av一级| 国产欧美日韩综合在线一区二区| 一个人免费看片子| 久久久精品94久久精品| 99re6热这里在线精品视频| 午夜福利视频精品| 国产免费福利视频在线观看| 久久国产精品影院| av福利片在线| 国产一区有黄有色的免费视频| 午夜久久久在线观看| 欧美激情 高清一区二区三区| 脱女人内裤的视频| 欧美日韩亚洲综合一区二区三区_| 国产老妇伦熟女老妇高清| 麻豆av在线久日| 亚洲精品国产区一区二| 少妇人妻久久综合中文| h视频一区二区三区| 亚洲情色 制服丝袜| 欧美日本中文国产一区发布| 大片免费播放器 马上看| 欧美xxⅹ黑人| 韩国高清视频一区二区三区| 黄色片一级片一级黄色片| 少妇精品久久久久久久| av天堂久久9| 日本黄色日本黄色录像| 91九色精品人成在线观看| 纵有疾风起免费观看全集完整版| 啦啦啦啦在线视频资源| 免费不卡黄色视频| 久久久精品免费免费高清| 我的亚洲天堂| 国产熟女欧美一区二区| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| 韩国精品一区二区三区| 国产不卡av网站在线观看| 久久热在线av| 国产男女超爽视频在线观看| 欧美精品亚洲一区二区| 中文字幕人妻熟女乱码| 搡老岳熟女国产| 亚洲精品国产色婷婷电影| 精品少妇黑人巨大在线播放| 国语对白做爰xxxⅹ性视频网站| 人成视频在线观看免费观看| 国产一级毛片在线| 在线av久久热| a级片在线免费高清观看视频| 丁香六月欧美| 久久精品国产亚洲av高清一级| 国产日韩欧美视频二区| 久久这里只有精品19| 青春草视频在线免费观看| 国产精品一区二区免费欧美 | 七月丁香在线播放| 久久久欧美国产精品| 高潮久久久久久久久久久不卡| 成在线人永久免费视频| 少妇 在线观看| 中文字幕人妻熟女乱码| 国产免费又黄又爽又色| 50天的宝宝边吃奶边哭怎么回事| 最近中文字幕2019免费版| 亚洲欧美精品综合一区二区三区| 中文字幕人妻熟女乱码| 欧美97在线视频| 久久热在线av| 久久中文字幕一级| 亚洲三区欧美一区| 七月丁香在线播放| 宅男免费午夜| 99久久99久久久精品蜜桃| 精品福利永久在线观看| 中文字幕亚洲精品专区| 大片电影免费在线观看免费| 99国产精品免费福利视频| 女人爽到高潮嗷嗷叫在线视频| 电影成人av| 国产欧美日韩一区二区三区在线| 国产亚洲欧美精品永久| 啦啦啦 在线观看视频| 国产精品熟女久久久久浪| 国产色视频综合| 国产野战对白在线观看| 亚洲av电影在线观看一区二区三区| 午夜福利在线免费观看网站| 久热爱精品视频在线9| 视频区图区小说| 日韩电影二区| 免费在线观看日本一区| 亚洲欧美精品自产自拍| 蜜桃国产av成人99| 91麻豆av在线| 在线天堂中文资源库| 五月天丁香电影| 日本色播在线视频| 欧美在线黄色| 免费在线观看日本一区| 纵有疾风起免费观看全集完整版| 中文字幕制服av| 美女扒开内裤让男人捅视频| 国产亚洲欧美在线一区二区| 九色亚洲精品在线播放| 国产精品秋霞免费鲁丝片| 80岁老熟妇乱子伦牲交| 大片免费播放器 马上看| 久久天躁狠狠躁夜夜2o2o | av网站在线播放免费| 国产在线视频一区二区| 一边摸一边抽搐一进一出视频| 最黄视频免费看| 国产真人三级小视频在线观看| 免费在线观看日本一区| 天天操日日干夜夜撸| 亚洲专区中文字幕在线| 欧美少妇被猛烈插入视频| 亚洲成色77777| 久久久久国产精品人妻一区二区| 久久这里只有精品19| 久久久欧美国产精品| 咕卡用的链子| 中文欧美无线码| 天天影视国产精品| av国产精品久久久久影院| 天天操日日干夜夜撸| 亚洲 国产 在线| 精品一区在线观看国产| 叶爱在线成人免费视频播放| 狂野欧美激情性bbbbbb| 两人在一起打扑克的视频| 亚洲av美国av| 亚洲精品日本国产第一区| 香蕉国产在线看| 国产一区有黄有色的免费视频| xxxhd国产人妻xxx| 久久久久久久久久久久大奶| 久久精品国产a三级三级三级| 制服人妻中文乱码| 欧美xxⅹ黑人| 18禁观看日本| 国产高清国产精品国产三级| 狠狠婷婷综合久久久久久88av| 999精品在线视频| 国产激情久久老熟女| 亚洲av电影在线进入| 国产精品一二三区在线看| 国产一区亚洲一区在线观看| 大香蕉久久网| 国产亚洲av高清不卡| 爱豆传媒免费全集在线观看| 黄色视频不卡| 亚洲精品成人av观看孕妇| 七月丁香在线播放| 国产一区二区三区综合在线观看| 亚洲欧美一区二区三区国产| 99久久精品国产亚洲精品| 中文字幕人妻丝袜一区二区| 欧美日韩av久久| 日本vs欧美在线观看视频| avwww免费| 精品卡一卡二卡四卡免费| 青草久久国产| 一区福利在线观看| 丰满少妇做爰视频| 波多野结衣一区麻豆| 精品人妻在线不人妻| 一本综合久久免费| 最黄视频免费看| kizo精华| 在线 av 中文字幕| 精品久久久久久电影网| 中文字幕人妻丝袜一区二区| 黄色视频在线播放观看不卡| 免费不卡黄色视频| 捣出白浆h1v1| 又紧又爽又黄一区二区| 国产精品人妻久久久影院| 日韩熟女老妇一区二区性免费视频| 亚洲中文字幕日韩| 高清av免费在线| 女性生殖器流出的白浆| 亚洲精品成人av观看孕妇| 91麻豆av在线| 夫妻性生交免费视频一级片| 蜜桃国产av成人99| www.熟女人妻精品国产| 精品免费久久久久久久清纯 | 欧美国产精品一级二级三级| 欧美激情高清一区二区三区|