亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        磁加速等離子體激勵(lì)器激勵(lì)特性

        2021-03-01 09:49:20陽(yáng)鵬宇賴慶仁
        氣體物理 2021年2期
        關(guān)鍵詞:電弧射流等離子體

        陽(yáng)鵬宇, 唐 坤, 張 鑫, 賴慶仁

        (1. 中國(guó)空氣動(dòng)力研究與發(fā)展中心空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 四川綿陽(yáng) 621000; 2. 中國(guó)空氣動(dòng)力研究與發(fā)展中心低速空氣動(dòng)力研究所, 四川綿陽(yáng) 621000)

        引 言

        等離子體主動(dòng)流動(dòng)控制技術(shù)以其響應(yīng)時(shí)間短、 激勵(lì)頻帶寬等技術(shù)優(yōu)點(diǎn), 在改善飛行器空氣動(dòng)力性能方面具有潛在應(yīng)用前景[1-5]. 分離流等離子體氣動(dòng)激勵(lì)控制研究是等離子體流動(dòng)控制研究領(lǐng)域中的重要分支. 大量研究表明等離子體在低風(fēng)速下可有效抑制翼型/機(jī)翼的前緣失速分離、 推遲失速迎角、 提高升阻比[6-11]. 但是在較高風(fēng)速下,等離子體的控制效果嚴(yán)重下降. 一個(gè)重要原因是高風(fēng)速下, 等離子體激勵(lì)器產(chǎn)生的擾動(dòng)強(qiáng)度降低. 如何增強(qiáng)激勵(lì)器的擾動(dòng)強(qiáng)度, 是等離子體激勵(lì)器在高風(fēng)速下實(shí)現(xiàn)分離流控制的關(guān)鍵.

        已有研究從等離子體激勵(lì)器的材料、 結(jié)構(gòu)、 布局和放電波形等方面進(jìn)行優(yōu)化, 以期望提高氣動(dòng)擾動(dòng)強(qiáng)度[12-18]. 另一種思路是引入磁場(chǎng), 改變等離子體激勵(lì)器工作時(shí)的能量分配比例. 具體講, 等離子體激勵(lì)器在工作時(shí), 放電電源輸出的電能主要轉(zhuǎn)換為氣體電離能, 氣體Joule熱和氣體動(dòng)能. 而通過(guò)磁場(chǎng)對(duì)等離子體施加Lorentz力, 可定向加速等離子體. 加速后的等離子體與中性粒子碰撞, 將更多動(dòng)能傳遞給氣體, 從而提高放電電能轉(zhuǎn)換為氣體動(dòng)能的比例, 實(shí)現(xiàn)等離子體氣動(dòng)激勵(lì)強(qiáng)度的提升.

        采用以上思路設(shè)計(jì)的激勵(lì)器結(jié)構(gòu)一般如圖1所示, 在壁面布置一對(duì)平行的正負(fù)電極, 沿垂直于壁面的方向施加磁場(chǎng), 磁場(chǎng)覆蓋電極所在區(qū)域. 當(dāng)正負(fù)電極間施加高電壓并出現(xiàn)放電時(shí), 產(chǎn)生的等離子體受到正交電磁場(chǎng)作用, 沿著電極軌道定向運(yùn)動(dòng), 并將動(dòng)能傳遞給周圍氣體, 形成氣動(dòng)激勵(lì)效果.

        利用磁加速等離子體改善航空/航天推進(jìn)器性能的研究已有很多報(bào)道, 例如俄羅斯提出的AJAX項(xiàng)目中采用磁流體加速提高超燃沖壓發(fā)動(dòng)機(jī)推進(jìn)效率[19], 航天領(lǐng)域中采用磁加速等離子體方式為航天器提供推力[20]. 對(duì)于提高等離子體氣動(dòng)激勵(lì)強(qiáng)度的研究, Pafford等利用磁場(chǎng)加速等離子體使放電電弧滑動(dòng)速度達(dá)到約100 m/s, 在來(lái)流16 m/s流場(chǎng)下誘導(dǎo)壁面射流達(dá)到32 m/s[21]. Zaidi等通過(guò)磁場(chǎng)加速方式將超聲速來(lái)流附面層中放電電弧的速度由350 m/s加速到2 000 m/s[22].

        已有研究結(jié)果表明通過(guò)引入磁場(chǎng)提高等離子體擾動(dòng)強(qiáng)度具有較大潛力. 然而, 目前, 此方面的研究還處于初步探索階段, 尤其是對(duì)磁加速等離子體激勵(lì)器誘導(dǎo)流場(chǎng)的研究鮮有報(bào)道. 本文針對(duì)不同磁場(chǎng)強(qiáng)度下放電等離子體的電學(xué)特性和誘導(dǎo)流場(chǎng)的結(jié)構(gòu)形態(tài)開展實(shí)驗(yàn)研究, 摸清磁場(chǎng)強(qiáng)度對(duì)激勵(lì)器電學(xué)特性與誘導(dǎo)流場(chǎng)特性的影響規(guī)律, 為提高等離子體在較高風(fēng)速下的控制效果提供支撐.

        1 實(shí)驗(yàn)系統(tǒng)

        實(shí)驗(yàn)系統(tǒng)由激勵(lì)器、 放電電源、 磁場(chǎng)產(chǎn)生設(shè)備和測(cè)量系統(tǒng)組成. 如圖2所示, 等離子體激勵(lì)器由陶瓷板和電極組成. 電極的內(nèi)邊緣呈現(xiàn)喇叭形狀, 以提高激勵(lì)器放電穩(wěn)定性, 內(nèi)邊緣長(zhǎng)度長(zhǎng)約 30 mm, 擴(kuò)張角約為6°, 最小間距dmin有 1 mm 和2.5 mm兩種. 放電電源輸出電壓為準(zhǔn)DC, 波形是半正弦波, 頻率為6 kHz, 峰值在6 kV以內(nèi)連續(xù)可調(diào). 磁場(chǎng)產(chǎn)生設(shè)備采用鐵釹硼磁鐵和電磁鐵, 磁場(chǎng)覆蓋激勵(lì)器的整個(gè)放電區(qū)域, 方向垂直于激勵(lì)器陶瓷板, 強(qiáng)度B在2 000 mT內(nèi)可調(diào). 測(cè)量系統(tǒng)采集獲得激勵(lì)器的電壓/電流波形、 磁場(chǎng)強(qiáng)度、 等離子體圖像和近壁面的誘導(dǎo)流場(chǎng). 等離子體圖像采用高速相機(jī)拍攝, 拍攝頻率為10 kfps. 近壁面誘導(dǎo)流場(chǎng)采用煙流流場(chǎng)顯示技術(shù)和PIV技術(shù)獲得, 示蹤粒子為橄欖油發(fā)煙顆粒, 拍攝頻率為0.05 kfps, 流場(chǎng)拍攝區(qū)域有橫向截面1和縱向截面2, 3(如圖2所示).

        2 激勵(lì)器放電特性

        2.1 放電形態(tài)

        圖3給出了激勵(lì)器工作時(shí)的典型放電圖像, 放電覆蓋的區(qū)域從電極喉道(電極間距最小處)到電極末端, 放電產(chǎn)生的等離子體貼于壁面.

        圖4給出了高速相機(jī)拍攝到的一系列放電圖像. 相鄰兩張圖像的時(shí)間間隔為0.1 ms. 由圖可知, 當(dāng)激勵(lì)器工作時(shí), 電極間實(shí)際產(chǎn)生了一條電弧. 電弧受到空間正交電磁場(chǎng)力的作用, 沿著電極從喉道滑動(dòng)到電極末端. 電弧在末端熄滅后重新在喉道產(chǎn)生, 并向下運(yùn)動(dòng), 由此重復(fù)進(jìn)行.

        圖3 激勵(lì)器工作時(shí)放電圖像Fig. 3 Discharge image of the plasma actuator

        圖4 B=270 mT, dmin=1 mm不同時(shí)刻電弧空間位置Fig. 4 Discharge arc position with B=270 mT, dmin=1 mm at different time instants

        2.2 等離子體定向運(yùn)動(dòng)速度

        根據(jù)高速相機(jī)拍攝的電弧圖像, 由不同時(shí)刻電弧的位置可計(jì)算電弧移動(dòng)速度, 即放電等離子體定向運(yùn)動(dòng)速度. 圖5展示了dmin=1 mm電極在不同磁場(chǎng)強(qiáng)度B下電弧的平均移動(dòng)速度v, 兩者大致成線性關(guān)系. 同樣dmin=2.5 mm時(shí)得到一組電弧速度與磁場(chǎng)強(qiáng)度的數(shù)據(jù)見表1. 與dmin=1 mm電極相比, 電弧移動(dòng)速度明顯提高, 例如在B=360 mT時(shí)v由3.92 m/s增加到4.98 m/s, 但對(duì)比同一磁場(chǎng)條件下測(cè)試的3組數(shù)據(jù)發(fā)現(xiàn)其速度波動(dòng)較大. 實(shí)驗(yàn)得到的最大電弧速度約6.76 m/s. 進(jìn)一步提高電弧移動(dòng)速度, 須增加磁場(chǎng)強(qiáng)度或極間電場(chǎng)強(qiáng)度. 實(shí)驗(yàn)發(fā)現(xiàn)提高電場(chǎng)強(qiáng)度, 電極容易燒蝕, 使放電穩(wěn)定性降低; 當(dāng)磁場(chǎng)強(qiáng)度大于1 000 mT時(shí), 電弧在電極上移動(dòng)出現(xiàn)不連續(xù)情況, 難以通過(guò)電弧的位置推算等離子體移動(dòng)速度.

        圖5 電弧移動(dòng)速度隨磁場(chǎng)強(qiáng)度的變化曲線(dmin=1 mm)Fig. 5 Curve of arc moving speed with magnetic field intensity(dmin=1 mm)

        表1 不同磁場(chǎng)強(qiáng)度電弧移動(dòng)速度(dmin=2.5 mm)

        2.3 電壓電流波形分析

        2.3.1 放電電壓電流特性

        圖6為典型條件下的放電電壓/電流波形圖, 可看到波形具有明顯周期性, 頻率為6 kHz, 這與放電電源輸出電壓頻率一致. 電壓/電流波形的一個(gè)周期從電流脈沖峰值開始, 此時(shí)電極間氣體剛被擊穿, 并產(chǎn)生電弧. 電弧使放電回路導(dǎo)通, 電源輸出電壓幾乎全部加載到限流電阻兩端, 電極間電壓從擊穿前的最大值直接下掉到接近零. 隨后電極間電壓逐漸升高, 回路電流則保持在幾十毫安量值, 對(duì)應(yīng)圖6中第1階段. 此階段電流波形穩(wěn)定, 量值較小, 可認(rèn)為放電處于正常輝光放電模式, 放電通道大致分為陰極區(qū)和正柱區(qū). 陰極區(qū)電壓在小電流狀態(tài)時(shí)由電極材質(zhì)決定, 在幾百伏量值, 不隨電流變化. 在6 kHz低頻條件下, 正柱區(qū)可等效為一純電阻, 阻值與放電等離子體密度、 遷移率等參數(shù)有關(guān). 因?yàn)殡姌O間電壓逐漸上升, 而電流基本保持不變, 表明正柱區(qū)電阻在不斷增大, 推知等離子體密度在逐漸減小. 當(dāng)?shù)入x子體密度減小到一定程度后, 放電通道變得不穩(wěn)定, 電壓波形出現(xiàn)鋸齒形波動(dòng), 而電流逐漸降低, 隨后降為零, 電弧熄滅, 放電

        進(jìn)入如圖6所示第2個(gè)階段. 在第2個(gè)階段, 電極間電壓逐漸升高, 當(dāng)達(dá)到擊穿電壓時(shí), 氣體再次被擊穿, 產(chǎn)生電弧, 放電進(jìn)入下一個(gè)周期.

        圖6 典型條件下激勵(lì)器放電電壓電流波形Fig. 6 Discharge voltage and current waveform of plasma actuator under typical conditions

        以上周期性的放電特征在電弧亮度上也有體現(xiàn), 圖7給出了時(shí)間間隔為0.05 ms的一組電弧圖像. 從圖上可以看出, 電弧亮度呈現(xiàn)強(qiáng)弱交替變化. 由于圖像采集頻率為20 kHz, 與放電電壓頻率6 kHz 并不是成倍數(shù), 所以圖片展示的電弧亮度變化存在相移, 周期性不明顯, 但基本反映了放電的周期性特征. 圖7中第7張照片上沒有電弧圖像, 該拍攝時(shí)刻對(duì)應(yīng)放電波形中第2階段, 此時(shí)電流為零, 電弧熄滅. 值得注意的是, 雖然電弧熄滅, 但該區(qū)域依然存在放電產(chǎn)生的殘余等離子體. 當(dāng)下一個(gè)放電周期的輸出電壓達(dá)到一定幅值時(shí), 殘余等離子體會(huì)促使電弧在當(dāng)?shù)禺a(chǎn)生, 這使得電弧在周期性產(chǎn)生和熄滅的同時(shí)依然沿著電極向下移動(dòng), 而不至于駐留在電極的喉道.

        圖7 間隔0.5 ms電弧發(fā)光強(qiáng)度變化Fig. 7 Change of arc luminous intensity attime interval of 0.5 ms

        2.3.2 不同參數(shù)下放電波形

        圖8為Un=4 000 V,B=400 mT條件得到的多個(gè)放電周期波形(Un為電源輸出有效電壓值,B為施加的磁場(chǎng)強(qiáng)度). 由圖可知, 波形的包絡(luò)也存在周期性. 每個(gè)周期內(nèi)包絡(luò)的電壓逐漸增大, 包絡(luò)的電流逐漸減小. 這一周期性的變化對(duì)應(yīng)圖4中電弧沿電極移動(dòng)的過(guò)程. 電弧向電極末端移動(dòng)時(shí)當(dāng)?shù)胤烹婇g距離增大, 用于維持放電的電壓也越大. 當(dāng)Un降到2 700 V時(shí), 由圖9可知, 維持電壓的包絡(luò)周期變短, 此時(shí)Un難以維持較大的極間距離放電, 使得電弧未移動(dòng)到電極的末端就返回電極喉道開始下一個(gè)周期的運(yùn)動(dòng).

        圖8 Un=4 000 V, B=400 mT放電電壓電流波形Fig. 8 Discharge voltage and current waveform at Un=4 000 V,B=400 mT

        圖9 Un=2 700 V, B=400 mT放電電壓電流波形Fig. 9 Discharge voltage and current waveform at Un=2 700 V,B=400 mT

        當(dāng)磁場(chǎng)強(qiáng)度由400 mT增加到1 000 mT時(shí), 由圖10, 11可知電壓波形包絡(luò)的原有周期性基本消失, 此時(shí)電弧不再周期性地沿著電極移動(dòng).

        圖10 Un=3 000 V, B=1 000 mT放電電壓電流波形Fig. 10 Discharge voltage and current waveform at Un=3 000 V,B=1 000 mT

        圖11 Un=4 000 V, B=1 000 mT放電電壓電流波形Fig. 11 Discharge voltage and current waveform at Un=4 000 V,B=1 000 mT

        根據(jù)電壓電流波形數(shù)據(jù)可由公式計(jì)算激勵(lì)器的能耗, 表2展示了4種工況下放電功耗大小, 初步規(guī)律是隨著Un和B的增大功耗會(huì)隨之增大. 消耗的能量一部分轉(zhuǎn)化為Joule熱, 一部分用于氣體電離所需, 一部分能量轉(zhuǎn)化為等離子體以及周圍氣體的定向運(yùn)動(dòng).

        表2 不同激勵(lì)條件下激勵(lì)器功耗

        3 激勵(lì)器誘導(dǎo)流場(chǎng)特性

        3.1 誘導(dǎo)流場(chǎng)形態(tài)

        圖12給出了流場(chǎng)拍攝截面1的煙流圖像. 從圖中可以看到激勵(lì)器誘導(dǎo)的流場(chǎng)包含了不同尺度的旋渦. 這些旋渦在電極周圍產(chǎn)生, 并貼著壁面運(yùn)動(dòng), 通過(guò)卷入上層氣體逐漸變大, 最后脫離壁面向上發(fā)展. 另外, 圖13和14展示了截面2和3的煙流圖像, 同樣可以看到有明顯的旋渦產(chǎn)生, 但旋渦僅僅分布在靠近陰極的一側(cè). 對(duì)比截面2和3的旋渦區(qū)域可知, 該方向的旋渦同時(shí)平行于電極向外發(fā)展. 結(jié)合以上3個(gè)截面的煙流圖像可知, 激勵(lì)器誘導(dǎo)的流場(chǎng)區(qū)域形狀如圖15陰影部分所示.

        該激勵(lì)器誘導(dǎo)渦產(chǎn)生的基本原理是通過(guò)等離子體定向運(yùn)動(dòng)誘導(dǎo)近壁面射流, 射流與周圍靜止氣體相互剪切形成旋渦. 根據(jù)橫向和縱向截面的煙流圖像, 可將激勵(lì)器誘導(dǎo)的射流細(xì)分為兩種, 一是電弧沿電極導(dǎo)軌滑動(dòng), 帶動(dòng)附近氣體運(yùn)動(dòng)形成沿縱向射流; 二是電弧形成過(guò)程中, 氣體擊穿產(chǎn)生的大量正離子從陽(yáng)極遷移到陰極, 帶動(dòng)附近氣體運(yùn)動(dòng)形成沿橫向射流. 兩種激勵(lì)方式同時(shí)作用, 使誘導(dǎo)的流場(chǎng)具有較強(qiáng)的三維效應(yīng). 另外, 由于激勵(lì)器工作時(shí)等離子體的產(chǎn)生和運(yùn)動(dòng)具有明顯的周期性, 射流具有脈沖特點(diǎn), 因此誘導(dǎo)的流場(chǎng)呈現(xiàn)強(qiáng)非定常性.

        圖12 典型條件下在截面1上的誘導(dǎo)煙流圖Fig. 12 Induced flow smoke visualization image of actuator on plane 1 under typical conditions

        圖13 典型激勵(lì)條件下在截面2上的誘導(dǎo)煙流圖Fig. 13 Induced flow smoke visualization image of conditionson plane 2 under typical condition

        圖14 典型激勵(lì)條件下在截面3上的誘導(dǎo)煙流圖Fig. 14 Induced flow smoke visualization image of actuator on plane 3 under typical conditions

        圖15 激勵(lì)器誘導(dǎo)流場(chǎng)區(qū)域形狀示意圖Fig. 15 Shape of flow field induced by actuator

        3.2 誘導(dǎo)流場(chǎng)強(qiáng)度

        采用PIV技術(shù)獲得了激勵(lì)器誘導(dǎo)流場(chǎng)的速度分布, 以此分析等離子體激勵(lì)強(qiáng)度. 圖16為截面1某個(gè)時(shí)刻的流場(chǎng)速度云圖和流線, 可以看到射流從左向右發(fā)展. 由于射流區(qū)域靠近壁面且速度變化大, 使得這一區(qū)域的示蹤粒子濃度非常低, 難以準(zhǔn)確計(jì)算出射流速度大小. 就當(dāng)前測(cè)試結(jié)果, 最大誘導(dǎo)氣流速度出現(xiàn)在圖中紅色區(qū)域, 約為 0.15 m/s. 圖17為截面3得到的速度場(chǎng), 得到了類似結(jié)果, 最大誘導(dǎo)氣流速度約0.2 m/s, 顯然測(cè)得的最大氣流速度并不能反映射流速度大小. 為進(jìn)一步確定激勵(lì)產(chǎn)生的射流速度大小, 實(shí)驗(yàn)采用了粒子跡線法進(jìn)行估算. 圖18為給定曝光時(shí)間典型激勵(lì)條件下的示蹤粒子移動(dòng)軌跡. 通過(guò)軌跡長(zhǎng)度(依據(jù)壁板尺寸得到)除以曝光時(shí)間可大致估算近壁面射流速度大小, 統(tǒng)計(jì)結(jié)果表明射流速度在1~2 m/s, 與等離子體定向運(yùn)動(dòng)速度(3~7 m/s)在一個(gè)量級(jí).

        圖16 典型激勵(lì)條件下在截面1上的誘導(dǎo)流場(chǎng)速度場(chǎng)圖Fig. 16 Velocity field on plane 1 under typical conditions

        圖17 典型激勵(lì)條件下在截面3上的誘導(dǎo)流場(chǎng)速度場(chǎng)圖Fig. 17 Velocity field on plane 3 under typical conditions

        圖18 截面1上的粒子跡線圖Fig. 18 Particle trace on plane 1

        鑒于該激勵(lì)器誘導(dǎo)的流場(chǎng)伴隨大量旋渦, 為高度非定常流, 因此僅通過(guò)射流大小評(píng)估激勵(lì)強(qiáng)度顯然不夠, 對(duì)此須引入時(shí)間量, 通過(guò)誘導(dǎo)流場(chǎng)變化的劇烈程度(擾動(dòng)能力)來(lái)評(píng)估激勵(lì)強(qiáng)度將更為合理, 此項(xiàng)工作要求流場(chǎng)采樣頻率和精度更高, 有待改進(jìn)測(cè)量手段后進(jìn)一步開展.

        4 結(jié)論

        當(dāng)激勵(lì)器以6 kHz頻率工作時(shí), 放電電弧與產(chǎn)生的等離子體在電磁場(chǎng)作用下沿電極定向運(yùn)動(dòng), 運(yùn)動(dòng)速度約5 m/s. 激勵(lì)器誘導(dǎo)流場(chǎng)存在兩種作用方式, 一是放電等離子體產(chǎn)生時(shí)誘導(dǎo)近壁面橫向射流, 二是等離子體沿電極運(yùn)動(dòng)時(shí)誘導(dǎo)近壁面縱向射流. 誘導(dǎo)射流速度約2 m/s. 兩種激勵(lì)方式同時(shí)作用, 誘導(dǎo)的流場(chǎng)具有明顯三維特征. 兩種激勵(lì)方式均為周期脈沖, 使得誘導(dǎo)的流場(chǎng)具有顯著非定常性.

        增大電源輸出電壓和增大磁場(chǎng)強(qiáng)度可以提高等離子體的激勵(lì)強(qiáng)度, 所產(chǎn)生的能耗也會(huì)隨之增大. 之后將通過(guò)改變激勵(lì)器材料、 優(yōu)化激勵(lì)器布局、 提高電場(chǎng)強(qiáng)度及磁場(chǎng)強(qiáng)度的方式, 進(jìn)一步提高磁加速等離子體激勵(lì)器的擾動(dòng)強(qiáng)度. 此外, 通過(guò)層析PIV對(duì)激勵(lì)器誘導(dǎo)三維流場(chǎng)進(jìn)行進(jìn)一步刻畫.

        猜你喜歡
        電弧射流等離子體
        深海逃逸艙射流注水均壓過(guò)程仿真分析
        低壓天然氣泄漏射流擴(kuò)散特性研究
        煤氣與熱力(2022年4期)2022-05-23 12:45:00
        故障電弧探測(cè)器與故障電弧保護(hù)裝置在工程中的應(yīng)用分析
        連續(xù)磁活動(dòng)對(duì)等離子體層演化的影響
        基于低溫等離子體修飾的PET/PVC浮選分離
        2219鋁合金激光電弧復(fù)合焊接及其溫度場(chǎng)的模擬
        等離子體種子處理技術(shù)介紹
        航空電氣系統(tǒng)中故障電弧的分析
        電子制作(2017年22期)2017-02-02 07:10:15
        射流齒形噴嘴射流流場(chǎng)與氣動(dòng)聲學(xué)分析
        地鐵站臺(tái)活塞風(fēng)附壁射流起始段的實(shí)測(cè)和實(shí)驗(yàn)驗(yàn)證
        黑人老外3p爽粗大免费看视频| 亚洲欧美在线视频| 中文字幕国内一区二区| 二区视频在线免费观看| 亚洲 中文 欧美 日韩 在线| 少妇高潮惨叫喷水在线观看| 无码精品人妻一区二区三区98| 日本精品一级二区三级| 中文字字幕人妻中文| 玩两个丰满老熟女| 日韩欧美亚洲国产一区二区三区| 精品福利一区二区三区| 国产精品videossex国产高清| 国产精品久久久久久人妻精品| 国产一线视频在线观看高清 | 人妻少妇猛烈井进入中文字幕| 成年女人黄小视频| 久久香蕉国产线看观看网| 久久熟女乱一区二区三区四区| 国产三级不卡一区不卡二区在线| 精品免费久久久久久久| 精品伊人久久香线蕉| 在线一区二区三区免费视频观看| 国产精品蝌蚪九色av综合网| 羞羞视频在线观看| 亚洲AV无码专区国产H小说| 亚州中文字幕乱码中文字幕| 看久久久久久a级毛片| 韩国19禁主播深夜福利视频| 国产精品一级av一区二区| 国产区女主播一区在线| 激情第一区仑乱| 久久久国产精品免费无卡顿| 免费人成在线观看播放视频| 国产偷国产偷精品高清尤物| 久久久国产精品ⅤA麻豆| 国产91熟女高潮一曲区| 国产在线无码精品无码| 久久99精品国产99久久6男男| 无码一区二区丝袜| 国产精品美女久久久网站三级|