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

    基于導(dǎo)流器噴水的導(dǎo)彈燃?xì)馍淞鹘禍匦Ч芯?/h1>
    2022-03-20 09:07:32張京力張曼曼
    關(guān)鍵詞:液態(tài)水導(dǎo)流壁面

    張京力,姜 毅,張曼曼,楊 瑩

    (北京理工大學(xué)宇航學(xué)院,北京 100081)

    0 引言

    車載導(dǎo)彈發(fā)射過程中,發(fā)動(dòng)機(jī)產(chǎn)生的高溫燃?xì)鈺?huì)對發(fā)射車和底部導(dǎo)流器等造成明顯的高溫?zé)g作用,這種熱沖擊效應(yīng)會(huì)給導(dǎo)彈安全發(fā)射帶來許多問題[1-2]。為了提高導(dǎo)彈發(fā)射安全性,延長發(fā)射裝置使用壽命,需要采取有效措施降低高溫燃?xì)鈱Πl(fā)射裝置的燒蝕作用。

    對于燃?xì)馍淞鲊娝禍匾?guī)律,國內(nèi)外學(xué)者已進(jìn)行了充分的理論與試驗(yàn)研究。美國的Geery等[3]最早開始了燃?xì)馍淞鲊娝禍貦C(jī)理的研究,通過改變?nèi)細(xì)饬髁亢蛧娝康谋壤菇禍匦Ч_(dá)到最佳,為燃?xì)馍淞鲊娝禍氐难芯康於嘶A(chǔ)。Giordan等[4]為火箭發(fā)射平臺設(shè)計(jì)了噴水降溫系統(tǒng),并利用數(shù)值仿真軟件對注水后的燃?xì)饬鲌鲞M(jìn)行了仿真計(jì)算,指出對燃?xì)饬鲌鲎⑺梢杂行p弱燃?xì)馍淞鞯臒釠_擊作用。Zhou等[5-7]對火箭燃?xì)饬鲌鲞M(jìn)行了數(shù)值模擬,并研究了不同噴水角度對流場降溫效果的影響,為發(fā)射場注水冷卻系統(tǒng)的設(shè)計(jì)提供理論參考。目前,燃?xì)馍淞鲊娝禍胤椒ㄒ呀?jīng)在美國肯尼迪航天發(fā)射中心得到應(yīng)用[8]。周帆、劉伯偉等[9-10]利用Mixture多相流模型,對多種發(fā)射方式下燃?xì)馍淞鲊娝鲌鲞M(jìn)行了數(shù)值仿真,得出了最優(yōu)化設(shè)計(jì)方法。紀(jì)添源[11]建立了運(yùn)載火箭縮比模型,引入離散項(xiàng)模型,對燃?xì)鈬娝鲌鲞M(jìn)行了數(shù)值仿真,研究了噴水強(qiáng)度、水滴直徑等參數(shù)對流場降溫效果的影響。上述研究均借用外部噴嘴、水管等設(shè)備進(jìn)行噴水降溫,占用空間較大,噴水系統(tǒng)較為復(fù)雜;同時(shí),將噴嘴、水管等設(shè)備直接安裝在發(fā)射裝置表面,可導(dǎo)致導(dǎo)彈發(fā)射系統(tǒng)的可靠性和穩(wěn)定性降低。

    基于以上分析,以單面導(dǎo)流器為例,文中提出一種噴水降溫方案,即在導(dǎo)流器內(nèi)布設(shè)豎直向上的噴水管道,利用液態(tài)水的汽化吸熱原理以及兩相流的沖擊作用,降低燃?xì)馍淞鲗?dǎo)流器和發(fā)射車壁面的燒蝕作用;同時(shí),在導(dǎo)流器表面直接開孔進(jìn)行噴水,不再使用外部噴嘴裝置,這大大降低了噴水系統(tǒng)的復(fù)雜程度,進(jìn)而提高導(dǎo)彈發(fā)射系統(tǒng)的可靠性。通過CFD仿真計(jì)算,對噴水裝置的降溫效果進(jìn)行研究,分析不同噴水速度與降溫效果之間的變化規(guī)律,得到最佳噴水方案,有效減弱高溫燃?xì)鈱Πl(fā)射裝置的燒蝕作用,進(jìn)一步提高發(fā)射設(shè)備的安全性。

    1 數(shù)值仿真理論基礎(chǔ)

    1.1 Mixture多相流模型

    Mixture多相流模型由Manninen等[12]在1996年提出,通過求解連續(xù)性方程、動(dòng)量方程、能量方程來模擬混合物中的各相,得到混合物各相的速度、溫度等參數(shù),通過計(jì)算體積分?jǐn)?shù)得到各相的分布,各相在計(jì)算單元內(nèi)的體積分?jǐn)?shù)之和為1。Mixture多相流模型能夠較好地處理高溫、高壓、可壓縮性氣液多相流動(dòng)等復(fù)雜問題[13]。因此,該模型能夠更加精確地模擬燃?xì)馍淞髋c液態(tài)水之間的相互作用。其中,連續(xù)性方程、動(dòng)量方程和能量方程如式(1)~式(3)所示[12]。

    連續(xù)性方程為:

    (1)

    式中:vm為質(zhì)量平均速度,ρm為混合物密度。

    混合物動(dòng)量方程可以由各個(gè)相的動(dòng)量方程疊加獲得:

    (2)

    式中:k表示第k相;F為體積力;μm為混合物粘性;vdr,k為次要相k的漂移速度。

    能量方程可以表示為:

    (3)

    式中:keff表示有效傳導(dǎo)率;SE表示體積熱源。

    式(3)中,對于可壓縮相有:

    (4)

    對于不可壓縮相有:

    Ek=hk

    (5)

    其中:Ek為第k相的能量;hk表示第k相的熱焓;vk為第k相的速度大小。

    1.2 液態(tài)水的物態(tài)變化方程

    高溫燃?xì)馀c液態(tài)水接觸時(shí),會(huì)發(fā)生相互作用。從發(fā)動(dòng)機(jī)噴出的燃?xì)鉁囟瓤梢赃_(dá)到3 000 K左右,因此在燃?xì)馀c液態(tài)水相互作用的區(qū)域會(huì)產(chǎn)生強(qiáng)烈的汽化現(xiàn)象。為了對液態(tài)水的汽化冷凝過程進(jìn)行數(shù)值模擬,將液態(tài)水的物態(tài)變化方程耦合到Mixture多相流模型中進(jìn)行計(jì)算。

    在實(shí)際計(jì)算中,通過計(jì)算每個(gè)單元格內(nèi)水的飽和溫度,即可得到水的汽化率。當(dāng)水的溫度高于飽和溫度時(shí),液態(tài)水吸收熱量并轉(zhuǎn)化為水蒸汽;當(dāng)水的溫度低于飽和溫度時(shí),水蒸汽放熱凝結(jié)為液態(tài)水。采用的液態(tài)水物態(tài)變化方程如式(6)、式(7)所示[14]。

    液態(tài)水汽化方程為:

    (6)

    水蒸汽凝結(jié)方程為:

    (7)

    2 計(jì)算模型

    2.1 物理模型

    研究的物理模型如圖1所示,主要由導(dǎo)彈、導(dǎo)彈發(fā)動(dòng)機(jī)噴管、發(fā)射箱、發(fā)射車壁、注水管道、導(dǎo)流器等組成。在導(dǎo)彈發(fā)射過程中,導(dǎo)彈發(fā)動(dòng)機(jī)噴管噴出高溫高壓燃?xì)?,?huì)對底部導(dǎo)流器和發(fā)射車壁等造成燒蝕作用?,F(xiàn)在導(dǎo)流器內(nèi)布設(shè)豎直向上的噴水管道,通過噴水管道向燃?xì)饬鲌鰢娚湟欢ㄋ俾实囊簯B(tài)水,利用液態(tài)水的汽化吸熱原理降低燃?xì)馍淞鞯臏囟龋瑫r(shí)液態(tài)水與燃?xì)馍淞鬟\(yùn)動(dòng)方向相反,可以進(jìn)一步降低高溫高壓燃?xì)鈱?dǎo)流器的直接沖擊及燒蝕作用。

    圖1 導(dǎo)流器噴水方案物理模型

    2.2 邊界條件設(shè)置

    圖2 對稱面網(wǎng)格劃分和邊界條件示意圖

    整個(gè)計(jì)算域分為入口邊界、出口邊界、壁面邊界和對稱面邊界。

    1)入口邊界條件:發(fā)動(dòng)機(jī)噴管入口截面設(shè)置為壓力入口,溫度為3 400 K,壓力為6.3 MPa;注水管入口截面設(shè)置為速度入口,噴水速度為20 m/s。

    2)出口邊界條件:計(jì)算域外邊界設(shè)為壓力出口,環(huán)境壓強(qiáng)為101 325 Pa,溫度為300 K。

    3)壁面邊界條件:發(fā)動(dòng)機(jī)噴管壁面、彈體表面、發(fā)射箱壁面、地面、導(dǎo)流器表面和注水管壁面均設(shè)置為無滑移絕熱壁面邊界條件。

    4)對稱面邊界條件:計(jì)算模型1/2對稱面位置。

    由于主要研究燃?xì)馍淞鲗Πl(fā)射車壁面和底部導(dǎo)流器的沖擊和燒蝕作用,在發(fā)射車壁面和導(dǎo)流器上分別選取3個(gè)溫度監(jiān)測點(diǎn),如圖3所示,分別為Point1~Point6。此外,將發(fā)射車壁面在對稱面上的投影線設(shè)為Line1,在Line1上等間距選取100個(gè)數(shù)據(jù)點(diǎn),用于后續(xù)網(wǎng)格無關(guān)性驗(yàn)證工作。

    圖3 監(jiān)測點(diǎn)位置示意圖

    2.3 網(wǎng)格無關(guān)性驗(yàn)證

    數(shù)值計(jì)算過程中,網(wǎng)格數(shù)量會(huì)影響計(jì)算的精確度?;诰W(wǎng)格疏密原則,共劃分了70萬、170萬、300萬和500萬4套網(wǎng)格。在圖3中的Line1上,由上而下等間距選取100個(gè)數(shù)據(jù)點(diǎn),對不同網(wǎng)格模型下數(shù)據(jù)點(diǎn)的溫度變化進(jìn)行對比分析,并對計(jì)算結(jié)果進(jìn)行后處理,畫出的溫度曲線如圖4所示。

    (5)實(shí)施維修方案。鏜缸、磨軸、鑲套、泄漏試驗(yàn)、探傷試驗(yàn)、動(dòng)平衡試驗(yàn)等基礎(chǔ)修復(fù)、清潔零部件、機(jī)體(含表面處理),吹干并測量檢測相關(guān)部件、組裝(含相關(guān)運(yùn)動(dòng)部件裝配間隙檢測、稱重、平衡試驗(yàn)、扭矩、扭曲、彎曲、平面度、平行度、漏光度、彈性、密封實(shí)驗(yàn)等)、盤機(jī)。

    圖4 Line1溫度分布曲線圖

    從圖4中可以看出,對于上述四套不同網(wǎng)格模型,Line1上的溫度總體變化趨勢相同。其中,70萬網(wǎng)格模型下的溫度曲線與其他三套網(wǎng)格模型下的計(jì)算結(jié)果相差較大,最大溫差約為140 K,計(jì)算誤差約為8%。170萬、300萬和500萬網(wǎng)格模型下的溫度曲線幾乎重合,但170萬網(wǎng)格模型的計(jì)算時(shí)間最短,其計(jì)算工時(shí)僅為300萬網(wǎng)格模型的1/2、500萬網(wǎng)格模型的1/3。綜合考慮計(jì)算精度和計(jì)算成本,選取170萬網(wǎng)格模型最為合適。

    2.4 數(shù)值模型試驗(yàn)驗(yàn)證

    將燃?xì)庾⑺鲌龇抡娼Y(jié)果與相關(guān)文獻(xiàn)實(shí)驗(yàn)結(jié)果進(jìn)行對比,確保數(shù)據(jù)仿真中所采用的邊界條件與實(shí)驗(yàn)中各參數(shù)保持一致,以驗(yàn)證所采用的數(shù)值模型的準(zhǔn)確性。圖5為文獻(xiàn)[15]所開展的噴水試驗(yàn),其與文中所研究的內(nèi)容具有很大的相似性,故以此實(shí)驗(yàn)作為參考,其中試驗(yàn)系統(tǒng)由縮比發(fā)動(dòng)機(jī)、固定臺架及噴水裝置等組成。發(fā)動(dòng)機(jī)點(diǎn)火前打開噴水管,液態(tài)水進(jìn)行相互交匯,點(diǎn)火后高溫燃?xì)馍淞鳑_入低溫液態(tài)水。整個(gè)過程由高速攝像機(jī)進(jìn)行記錄。

    圖5 噴水試驗(yàn)過程圖片

    圖6為噴水試驗(yàn)高速攝影圖與數(shù)值計(jì)算結(jié)果的溫度分布云圖之間的對比,可以發(fā)現(xiàn)兩者流場形態(tài)非常接近,且波節(jié)位置幾乎一致,即仿真計(jì)算得到的結(jié)果與試驗(yàn)現(xiàn)象十分吻合。該圖驗(yàn)證了采用耦合液態(tài)水汽化方程的Mixture多相流模型計(jì)算該類問題是可行的,仿真結(jié)果具有較好的可信度。

    圖6 高速攝影圖與溫度云圖對比

    3 仿真結(jié)果分析

    3.1 未噴水工況與噴水工況結(jié)果對比分析

    圖7展示了未噴水工況與噴水工況對稱面上溫度分布云圖。圖7(a)中可以看到明顯的激波結(jié)構(gòu),燃?xì)馍淞鹘?jīng)過激波后,溫度由1 300 K升高至2 200 K左右;當(dāng)燃?xì)馍淞鞯竭_(dá)導(dǎo)流器表面后,導(dǎo)流器有效地將噴管噴出的高溫燃?xì)鈱?dǎo)向了左側(cè),但是仍有一部分高溫燃?xì)庋刂l(fā)射車壁面進(jìn)行擴(kuò)散,使得發(fā)射車壁面靠近底部區(qū)域處于較高溫度環(huán)境中;由于燃?xì)馍淞鳑_擊導(dǎo)流器表面后速度明顯下降,產(chǎn)生滯止現(xiàn)象,因此導(dǎo)流器表面處燃?xì)鉁囟壬仙_(dá)到總溫3 400 K左右,導(dǎo)流器受到燃?xì)馍淞鞯臒g作用明顯。圖7(b)中,經(jīng)過導(dǎo)流器上的豎直管道進(jìn)行噴水之后,對稱面上燃?xì)饬鲌龅母邷貐^(qū)域明顯縮小,燃?xì)馍淞髟诎l(fā)射車壁面處的擴(kuò)散減弱,發(fā)射車壁面區(qū)域已基本沒有高溫流動(dòng);同時(shí)導(dǎo)流器表面最高溫度下降到2 100 K左右。結(jié)合圖8可以得出,向上噴射的液態(tài)水與燃?xì)馍淞靼l(fā)生碰撞后,在導(dǎo)流器表面和發(fā)射車壁面形成了一層水膜,水膜有效阻擋了燃?xì)馍淞鲗?dǎo)流器和發(fā)射車壁的直接燒蝕作用,使得發(fā)射車壁面和導(dǎo)流器大部分區(qū)域基本沒有高溫流動(dòng)。

    圖7 未噴水與噴水工況對稱面溫度分布云圖

    圖8 對稱面水體積分?jǐn)?shù)分布云圖

    圖9為未噴水工況與噴水工況三維流場溫度分布云圖。從圖9(a)中可以看出,由于高溫燃?xì)獾臒g作用,導(dǎo)流器表面溫度整體偏高,并且越靠近中心區(qū)域溫度越高。此外,高溫燃?xì)庠诎l(fā)射車壁面區(qū)域發(fā)生擴(kuò)散,使得發(fā)射車壁面底部區(qū)域的溫度達(dá)到了1 500 K左右。相比于圖9(a),圖9(b)中三維流場的高溫區(qū)域明顯縮小,主要集中在導(dǎo)流器表面中部區(qū)域,并且最高溫度明顯降低。同時(shí),燃?xì)馍淞鲗Πl(fā)射車壁面的燒蝕作用大幅減弱,發(fā)射車壁面溫度基本保持在700 K以下。

    圖9 未噴水工況與噴水工況三維流場溫度分布云圖

    圖10為未噴水工況與噴水工況下發(fā)射車壁面上監(jiān)測點(diǎn)Point1、Point2和Point3溫度變化曲線圖。從圖10中可以看出,未噴水工況下監(jiān)測點(diǎn)Point1、Point2和Point3的溫度上升時(shí)間有明顯差異,監(jiān)測點(diǎn)Point3在0.01 s左右迅速上升至1 800 K以上,這是由于燃?xì)馍淞髟诎l(fā)射車壁面底部的擴(kuò)散較慢,監(jiān)測點(diǎn)Point2升溫時(shí)間滯后于監(jiān)測點(diǎn)Point3,監(jiān)測點(diǎn)Point1次之;當(dāng)流場穩(wěn)定后,監(jiān)測點(diǎn)Point1、監(jiān)測點(diǎn)Point2和監(jiān)測點(diǎn)Point3的溫度分別穩(wěn)定在900 K、1 100 K和1 500 K左右。經(jīng)過導(dǎo)流器內(nèi)布設(shè)的豎直管道噴水后,發(fā)射車壁面溫度變化無明顯振蕩波動(dòng)現(xiàn)象,3個(gè)監(jiān)測點(diǎn)的溫度均穩(wěn)定在400 K以下。

    圖10 未噴水工況與噴水工況監(jiān)測點(diǎn)Point1、Point2和Point3溫度變化曲線圖

    圖11為未噴水工況與噴水工況下導(dǎo)流器表面監(jiān)測點(diǎn)Point4、Point5和Point6溫度變化曲線圖。從圖11中可以看出,未噴水工況下,監(jiān)測點(diǎn)Point4、Point5和Point6的溫度在0.01 s內(nèi)迅速上升至2 000 K以上,在0.15 s左右燃?xì)饬鲌鲞_(dá)到穩(wěn)定,各監(jiān)測點(diǎn)溫度保持不變,其中監(jiān)測點(diǎn)Point5位置的溫度最高,高達(dá)3 100 K。對高溫燃?xì)饬鲌鲞M(jìn)行噴水降溫之后,導(dǎo)流器表面上各監(jiān)測點(diǎn)溫度下降明顯,其中監(jiān)測點(diǎn)Point4、Point6溫度下降了2 000 K以上,監(jiān)測點(diǎn)Point5溫度大約降低了900 K。

    圖11 未噴水工況與噴水工況監(jiān)測點(diǎn)Point4、Point5和Point6溫度變化曲線圖

    表1為未噴水和噴水兩種工況下6個(gè)監(jiān)測點(diǎn)最高溫度對比。經(jīng)過導(dǎo)流器內(nèi)的豎直管道噴水之后,燃?xì)馍淞鲗?dǎo)流器表面和發(fā)射車壁面的燒蝕作用明顯減弱,監(jiān)測點(diǎn)Point1~Point4,Point6的最高溫度均下降了1 000 K以上,而監(jiān)測點(diǎn)Point5溫度下降程度較小。結(jié)合圖8可看出,監(jiān)測點(diǎn)Point5所在區(qū)域沒有形成水膜,導(dǎo)致高溫燃?xì)庵苯記_擊該區(qū)域,故溫度明顯高于其他區(qū)域溫度;由于液態(tài)水的汽化吸熱作用,監(jiān)測點(diǎn)Point5的溫度峰值下降了大約500 K。

    表1 未噴水與噴水工況監(jiān)測點(diǎn)最高溫度對比 單位:K

    通過以上分析可知,利用導(dǎo)流器內(nèi)布設(shè)的豎直向上管道對高溫燃?xì)饬鲌鲞M(jìn)行噴水,可對導(dǎo)流器和發(fā)射車壁等發(fā)射設(shè)備起到有效保護(hù)作用。

    3.2 噴水降溫優(yōu)化分析

    液態(tài)水豎直向上噴射的過程中,會(huì)與向下噴射的燃?xì)馍淞靼l(fā)生碰撞,兩者發(fā)生動(dòng)量交換。在燃?xì)馍淞鲃?dòng)量不變的情況下,可通過調(diào)節(jié)噴水速度,改變液態(tài)水的動(dòng)量,進(jìn)而改變?nèi)細(xì)馍淞鞯牧鲃?dòng)范圍。同時(shí)通過改變噴水速度,可以實(shí)現(xiàn)液態(tài)水汽化吸熱速率的變化,提高降溫效果?,F(xiàn)選取噴水速度10 m/s、20 m/s、40 m/s、50 m/s、55 m/s、60 m/s六種典型工況進(jìn)行對比分析,研究不同噴水速度與降溫效果之間的變化規(guī)律。

    對計(jì)算數(shù)據(jù)進(jìn)行后處理得到圖12、圖13的結(jié)果。圖12為噴水工況下對稱面液態(tài)水流線圖,圖13為噴水工況下對稱面燃?xì)饬骶€圖。

    圖12 不同噴水工況下對稱面液態(tài)水流線圖

    圖13 不同噴水工況下對稱面燃?xì)饬骶€圖

    從圖12中可以看出,隨著噴水速度逐漸增加,導(dǎo)流器表面和發(fā)射車壁面處的水膜厚度逐漸增加,這是由于噴水速度增大,使液態(tài)水的沖量變大,導(dǎo)致導(dǎo)流器表面上液態(tài)水與燃?xì)庾饔妹娴奈恢貌粩嗌弦?;同時(shí),噴水速度增大后,相同時(shí)間內(nèi)注入流場的液態(tài)水的量變多,在燃?xì)獾臎_擊作用下,發(fā)射車壁面形成了更厚的水膜。從圖13中可以看出,水膜對高溫燃?xì)馄鸬搅溯^好的隔離作用,并且水膜越厚,隔離效果越好。然而,噴水速度達(dá)到55 m/s時(shí),繼續(xù)增大噴水速度,導(dǎo)流器與發(fā)射車壁面的水膜厚度變化不明顯。

    圖14為未噴水工況與6種噴水工況下監(jiān)測點(diǎn)溫度變化曲線圖。可以看出,相較于未噴水工況,噴水后,監(jiān)測點(diǎn)Point1~Point4,Point6的溫度明顯降低,待流場穩(wěn)定后,5個(gè)監(jiān)測點(diǎn)溫度均穩(wěn)定在600 K以下,監(jiān)測點(diǎn)Point5溫度下降了800 K以上,這說明該噴水方案有效降低了高溫燃?xì)鈱?dǎo)流器和發(fā)射車壁的燒蝕作用。此外,隨著噴水速度的增加,監(jiān)測點(diǎn)Point1、Point2和Point3的溫度變化較小,基本保持在400 K以下;監(jiān)測點(diǎn)Point4和Point6的溫度隨著噴水速度的增加而減小;流場穩(wěn)定后,隨著噴水速度增加,Point5的溫度整體上呈下降趨勢,這說明燃?xì)饬鲌龅慕禍匦ЧS噴水速度的增加整體呈現(xiàn)變好趨勢。

    圖14 監(jiān)測點(diǎn)溫度變化曲線圖

    對于監(jiān)測點(diǎn)Point5,從圖14(e)中可以看出,噴水后溫度曲線出現(xiàn)波峰,之后達(dá)到穩(wěn)定,這是由于Point5所在區(qū)域?yàn)槿細(xì)庵苯記_擊區(qū)域,燃?xì)鈴膰姽車姵龊?,可在很短時(shí)間內(nèi)到達(dá)Point5所在區(qū)域,未能與液態(tài)水發(fā)生充分作用,且燃?xì)馇熬墯饬鳒囟容^高,使得監(jiān)測點(diǎn)Point5溫度出現(xiàn)波峰;隨后,氣液兩相流的耦合作用達(dá)到動(dòng)態(tài)平衡,監(jiān)測點(diǎn)5溫度達(dá)到穩(wěn)定狀態(tài)。隨著噴水速度的增加,監(jiān)測點(diǎn)Point5溫度峰值出現(xiàn)時(shí)間發(fā)生滯后,這是由于液態(tài)水的沖量隨噴水速度的改變而改變,液態(tài)水速度越大,對高溫燃?xì)獾臎_擊作用越明顯,高溫燃?xì)獾竭_(dá)Point5所在區(qū)域的時(shí)間延長。表2為在不同噴水速度工況下監(jiān)測點(diǎn)Point5處液態(tài)水的汽化速率。由表2可看出,隨著噴水速度的增大,監(jiān)測點(diǎn)Point5處的汽化速率先減小后增大;其中,噴水速度55 m/s工況下Point5處汽化速率最大,相同時(shí)間內(nèi)吸熱量最多,這與圖14(e)監(jiān)測點(diǎn)Point5溫度變化曲線中噴水速度55 m/s工況時(shí)溫度最低現(xiàn)象一致。

    表2 不同工況Point5汽化速率對比

    通過以上分析可知,高溫燃?xì)獾膰娝禍匦ЧQ于燃?xì)馀c液態(tài)水的動(dòng)量交換和能量交換。一方面,高溫燃?xì)馀c液態(tài)水的動(dòng)量交換會(huì)限制燃?xì)獾牧鲃?dòng)范圍,避免燃?xì)庵苯記_擊導(dǎo)流器表面和發(fā)射車壁面;另一方面,通過液態(tài)水的汽化相變吸熱,降低高溫燃?xì)獾臏囟?,從而減弱燃?xì)鈱Πl(fā)射裝置的燒蝕作用。通過圖14中6個(gè)監(jiān)測點(diǎn)的溫度變化曲線可以看出,在噴水速度55 m/s工況下,降溫效果最佳。

    4 結(jié)論

    為實(shí)現(xiàn)發(fā)射裝置熱防護(hù),提出了一種在導(dǎo)流器內(nèi)布設(shè)豎直向上噴水管道的新型方案,解決導(dǎo)彈發(fā)射過程中發(fā)射裝置的燒蝕問題;由于不再使用外部噴嘴裝置,進(jìn)而降低了噴水系統(tǒng)復(fù)雜程度?;贛ixture多相流模型,耦合液態(tài)水的汽化方程與組分輸運(yùn)模型,對不同噴水速度工況下的燃?xì)饬鲌鲞M(jìn)行CFD仿真計(jì)算,分析該方案的降溫機(jī)理及效果。主要結(jié)論如下:

    1)在導(dǎo)流器內(nèi)布設(shè)豎直向上的噴水管道方案可對導(dǎo)流器和發(fā)射車壁等發(fā)射設(shè)備起到有效降溫作用。導(dǎo)流器內(nèi)的豎直管道噴水后,導(dǎo)流器表面和發(fā)射車壁面形成水膜,對高溫燃?xì)馄鸬搅烁綦x作用,導(dǎo)流器表面和發(fā)射車壁面等發(fā)射裝置的溫度明顯降低。

    2)對比未噴水流場,對燃?xì)饬鲌鲞M(jìn)行噴水后,發(fā)射車壁面和導(dǎo)流器大部分區(qū)域已基本無高溫流動(dòng)。噴水后的燃?xì)饬鲌鲋?,高溫燃?xì)饬鲃?dòng)主要集中在導(dǎo)流器表面中部區(qū)域,且最高溫度從3 100 K下降到2 100 K左右,降溫幅度達(dá)到30%以上;在發(fā)射車壁面處,燃?xì)馍淞鞯臄U(kuò)散明顯減弱,發(fā)射車壁面處基本沒有高溫流動(dòng),其溫度基本保持在700 K以下,降溫效果明顯。

    3)隨著噴水速度的增大,流場降溫效果呈現(xiàn)出先變好再變差的趨勢,其中噴水速度55 m/s工況降溫效果最佳。噴水速度越大,導(dǎo)流器表面和發(fā)射車壁面形成的水膜厚度越厚,對燃?xì)獾母綦x作用不斷提高;當(dāng)噴水速度達(dá)到55 m/s時(shí),燃?xì)馍淞髦苯記_擊區(qū)域的液態(tài)水汽化速率最大,相同時(shí)間內(nèi)吸熱量最多,該區(qū)域的溫度達(dá)到最低,故噴水速度55 m/s對流場的降溫效果最佳。

    猜你喜歡
    液態(tài)水導(dǎo)流壁面
    二維有限長度柔性壁面上T-S波演化的數(shù)值研究
    基于微波輻射計(jì)的張掖地區(qū)水汽、液態(tài)水變化特征分析
    Ka/Ku雙波段毫米波雷達(dá)功率譜數(shù)據(jù)反演液態(tài)水含量方法研究
    新型分支相間導(dǎo)流排
    零下溫度的液態(tài)水
    某水利樞紐工程施工導(dǎo)流及水流控制研究
    PEMFC氣體擴(kuò)散層中液態(tài)水傳輸實(shí)驗(yàn)研究綜述
    導(dǎo)流堤在小型引水工程中的應(yīng)用
    壁面溫度對微型內(nèi)燃機(jī)燃燒特性的影響
    大崗山水電站導(dǎo)流隧洞下閘后滲漏處理

    麻豆乱淫一区二区| 国产成人免费观看mmmm| 永久免费av网站大全| 99国产精品一区二区蜜桃av | 麻豆av在线久日| 制服诱惑二区| 亚洲国产av新网站| 爱豆传媒免费全集在线观看| 免费不卡黄色视频| 国产97色在线日韩免费| 一边亲一边摸免费视频| av国产精品久久久久影院| 一边摸一边做爽爽视频免费| 精品亚洲成国产av| 免费高清在线观看日韩| 国产国语露脸激情在线看| 在线观看国产h片| a 毛片基地| 人人妻,人人澡人人爽秒播 | 在线观看国产h片| 欧美精品亚洲一区二区| 日韩欧美一区视频在线观看| 亚洲国产欧美一区二区综合| 激情五月婷婷亚洲| 秋霞在线观看毛片| av视频免费观看在线观看| 欧美在线一区亚洲| 国产精品一区二区在线不卡| 建设人人有责人人尽责人人享有的| 成人午夜精彩视频在线观看| 老鸭窝网址在线观看| 高清av免费在线| 久久久久网色| 日韩 欧美 亚洲 中文字幕| 精品国产乱码久久久久久男人| 少妇人妻 视频| 久久久久久久久免费视频了| 日韩一区二区三区影片| 一本久久精品| 十八禁高潮呻吟视频| 99久久99久久久精品蜜桃| 在现免费观看毛片| 精品少妇一区二区三区视频日本电影| 亚洲精品在线美女| av在线老鸭窝| 欧美黑人欧美精品刺激| 亚洲av欧美aⅴ国产| 一区二区av电影网| av又黄又爽大尺度在线免费看| xxx大片免费视频| 精品欧美一区二区三区在线| 国产高清videossex| 久久 成人 亚洲| 老司机深夜福利视频在线观看 | 亚洲欧美清纯卡通| 精品少妇一区二区三区视频日本电影| 亚洲av成人不卡在线观看播放网 | 国产成人一区二区在线| 黑丝袜美女国产一区| 国产成人影院久久av| 男女免费视频国产| 亚洲国产精品成人久久小说| 手机成人av网站| 国产无遮挡羞羞视频在线观看| 国产亚洲精品第一综合不卡| 欧美激情 高清一区二区三区| 久久国产精品大桥未久av| 又大又爽又粗| 日韩一卡2卡3卡4卡2021年| 多毛熟女@视频| 男男h啪啪无遮挡| 老司机在亚洲福利影院| 亚洲精品国产av蜜桃| 80岁老熟妇乱子伦牲交| 夫妻性生交免费视频一级片| 一本色道久久久久久精品综合| 高清不卡的av网站| 97在线人人人人妻| 天天操日日干夜夜撸| 久久久国产一区二区| a 毛片基地| 午夜免费男女啪啪视频观看| 69精品国产乱码久久久| 人妻人人澡人人爽人人| 国产精品久久久人人做人人爽| 1024视频免费在线观看| 日日夜夜操网爽| 狠狠精品人妻久久久久久综合| 国产在线免费精品| 亚洲精品乱久久久久久| 国产亚洲av高清不卡| 久久亚洲国产成人精品v| 一区二区三区精品91| 欧美精品av麻豆av| 超碰97精品在线观看| 在线天堂中文资源库| 99热国产这里只有精品6| 午夜福利视频精品| 国产淫语在线视频| 国产精品成人在线| 黄片播放在线免费| 女人高潮潮喷娇喘18禁视频| 黄频高清免费视频| 国产日韩一区二区三区精品不卡| 午夜视频精品福利| 看十八女毛片水多多多| 成人国产av品久久久| 69精品国产乱码久久久| 午夜福利视频精品| 国产在线观看jvid| cao死你这个sao货| 国产高清videossex| 久久久久久免费高清国产稀缺| 中文字幕人妻丝袜一区二区| 18禁裸乳无遮挡动漫免费视频| 中文乱码字字幕精品一区二区三区| 欧美黑人精品巨大| bbb黄色大片| 天天躁日日躁夜夜躁夜夜| 国产xxxxx性猛交| 日韩一本色道免费dvd| 国产日韩一区二区三区精品不卡| 午夜精品国产一区二区电影| 夫妻性生交免费视频一级片| 亚洲精品国产一区二区精华液| 一边摸一边做爽爽视频免费| 99精国产麻豆久久婷婷| 一级毛片电影观看| 亚洲av在线观看美女高潮| 国产免费福利视频在线观看| 黑人猛操日本美女一级片| 一区二区三区激情视频| 欧美日韩亚洲综合一区二区三区_| 国产男人的电影天堂91| 国产精品久久久久久精品古装| 视频区欧美日本亚洲| 亚洲成国产人片在线观看| 9热在线视频观看99| 亚洲成人国产一区在线观看 | 欧美av亚洲av综合av国产av| 免费人妻精品一区二区三区视频| 欧美乱码精品一区二区三区| 免费在线观看完整版高清| 18禁黄网站禁片午夜丰满| 人妻人人澡人人爽人人| 亚洲,欧美,日韩| 少妇人妻 视频| 色婷婷久久久亚洲欧美| 久久鲁丝午夜福利片| 日韩一卡2卡3卡4卡2021年| 男人爽女人下面视频在线观看| 宅男免费午夜| 国产成人av教育| 老司机影院成人| 久久天堂一区二区三区四区| 色婷婷av一区二区三区视频| 日韩一区二区三区影片| 精品久久久久久久毛片微露脸 | 天天操日日干夜夜撸| 99香蕉大伊视频| av国产久精品久网站免费入址| 亚洲av综合色区一区| 免费在线观看完整版高清| 又大又黄又爽视频免费| 一级a爱视频在线免费观看| 高清视频免费观看一区二区| 首页视频小说图片口味搜索 | 欧美日韩一级在线毛片| a 毛片基地| 国产在线视频一区二区| 国产精品免费大片| 老熟女久久久| 嫁个100分男人电影在线观看 | 伊人亚洲综合成人网| 精品少妇久久久久久888优播| 午夜免费成人在线视频| 又大又爽又粗| 国产精品熟女久久久久浪| 国产欧美日韩一区二区三 | 久久久国产一区二区| 一边摸一边做爽爽视频免费| 美女中出高潮动态图| 色网站视频免费| 亚洲五月婷婷丁香| www.熟女人妻精品国产| 成人国产一区最新在线观看 | 亚洲精品在线美女| 国产熟女欧美一区二区| 日本av手机在线免费观看| av在线app专区| 性色av乱码一区二区三区2| 丰满迷人的少妇在线观看| 久久精品国产a三级三级三级| cao死你这个sao货| 91精品国产国语对白视频| 国产免费又黄又爽又色| 日本色播在线视频| 国产精品99久久99久久久不卡| 亚洲精品美女久久久久99蜜臀 | 91麻豆精品激情在线观看国产 | 免费在线观看视频国产中文字幕亚洲 | 久久人人97超碰香蕉20202| 在线看a的网站| 久久久久精品国产欧美久久久 | 精品一区二区三卡| av在线播放精品| 国产亚洲av片在线观看秒播厂| 欧美日韩成人在线一区二区| 老司机深夜福利视频在线观看 | 中文字幕人妻熟女乱码| 高清欧美精品videossex| 宅男免费午夜| 免费久久久久久久精品成人欧美视频| 亚洲精品一区蜜桃| 日本色播在线视频| 精品福利观看| 肉色欧美久久久久久久蜜桃| 黄片播放在线免费| av视频免费观看在线观看| 久久这里只有精品19| 国产精品99久久99久久久不卡| 日韩熟女老妇一区二区性免费视频| 欧美中文综合在线视频| 精品人妻1区二区| 国产黄色免费在线视频| 欧美乱码精品一区二区三区| 亚洲国产日韩一区二区| 青春草亚洲视频在线观看| 国产在线一区二区三区精| 国产黄色免费在线视频| av不卡在线播放| 深夜精品福利| www.精华液| 成年av动漫网址| 日韩av在线免费看完整版不卡| 精品少妇内射三级| 91麻豆精品激情在线观看国产 | 曰老女人黄片| 你懂的网址亚洲精品在线观看| 两个人看的免费小视频| 国产成人系列免费观看| 黄色片一级片一级黄色片| 亚洲精品一卡2卡三卡4卡5卡 | 婷婷色麻豆天堂久久| 精品卡一卡二卡四卡免费| 亚洲精品一二三| 成人国语在线视频| 成人国产一区最新在线观看 | 亚洲精品一区蜜桃| 精品国产乱码久久久久久小说| 午夜激情av网站| 亚洲av电影在线进入| 亚洲国产精品一区三区| 久久久精品免费免费高清| 日本wwww免费看| 丰满饥渴人妻一区二区三| 性色av乱码一区二区三区2| 一区二区三区四区激情视频| 黄片小视频在线播放| 久久久久国产一级毛片高清牌| 中文字幕制服av| 视频在线观看一区二区三区| 国产片特级美女逼逼视频| 国产精品一区二区在线观看99| 中文字幕最新亚洲高清| 精品视频人人做人人爽| 中文字幕色久视频| av天堂在线播放| 两性夫妻黄色片| 黄色视频在线播放观看不卡| 高清黄色对白视频在线免费看| 国产1区2区3区精品| 亚洲人成网站在线观看播放| 中文精品一卡2卡3卡4更新| 狠狠精品人妻久久久久久综合| 日韩大码丰满熟妇| 国产精品免费视频内射| av网站免费在线观看视频| 一级毛片我不卡| 国产一区二区三区综合在线观看| 黄片播放在线免费| 欧美激情 高清一区二区三区| 午夜激情久久久久久久| 国产成人系列免费观看| 18禁裸乳无遮挡动漫免费视频| 999久久久国产精品视频| 国产精品国产三级专区第一集| 老司机亚洲免费影院| 色94色欧美一区二区| 黄色a级毛片大全视频| 午夜免费观看性视频| 操美女的视频在线观看| 国产一级毛片在线| 国产色视频综合| 一边摸一边抽搐一进一出视频| 夫妻性生交免费视频一级片| 亚洲中文日韩欧美视频| 岛国毛片在线播放| 免费在线观看日本一区| www.999成人在线观看| www.精华液| 国产成人系列免费观看| 中文字幕色久视频| 国产精品久久久av美女十八| 女人精品久久久久毛片| 精品亚洲成国产av| 欧美黄色淫秽网站| 女人高潮潮喷娇喘18禁视频| 99久久精品国产亚洲精品| 蜜桃国产av成人99| 中文欧美无线码| 水蜜桃什么品种好| 一级,二级,三级黄色视频| 免费日韩欧美在线观看| 狠狠婷婷综合久久久久久88av| 真人做人爱边吃奶动态| 看十八女毛片水多多多| 亚洲中文日韩欧美视频| 久久女婷五月综合色啪小说| 99国产精品99久久久久| 五月开心婷婷网| 高潮久久久久久久久久久不卡| 亚洲成人国产一区在线观看 | 国产一区二区激情短视频 | 欧美日韩亚洲高清精品| 成人黄色视频免费在线看| 日韩一区二区三区影片| 免费人妻精品一区二区三区视频| 少妇 在线观看| 男女免费视频国产| 美女国产高潮福利片在线看| av网站免费在线观看视频| 国产精品熟女久久久久浪| 国产熟女欧美一区二区| 欧美人与性动交α欧美软件| 国产男女超爽视频在线观看| 久久综合国产亚洲精品| kizo精华| 婷婷丁香在线五月| 建设人人有责人人尽责人人享有的| 久久精品国产亚洲av涩爱| 亚洲美女黄色视频免费看| 欧美另类一区| 精品少妇一区二区三区视频日本电影| 性色av乱码一区二区三区2| 精品卡一卡二卡四卡免费| 19禁男女啪啪无遮挡网站| 国产日韩一区二区三区精品不卡| kizo精华| 亚洲av电影在线观看一区二区三区| 成年动漫av网址| 日韩中文字幕视频在线看片| tube8黄色片| 免费不卡黄色视频| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 无遮挡黄片免费观看| 精品一区二区三区四区五区乱码 | 国产成人精品久久二区二区91| 亚洲国产欧美网| 亚洲av日韩精品久久久久久密 | 欧美97在线视频| 亚洲成国产人片在线观看| 人人妻,人人澡人人爽秒播 | 啦啦啦啦在线视频资源| 一区二区三区精品91| 国产精品免费视频内射| 久久精品成人免费网站| 一级a爱视频在线免费观看| cao死你这个sao货| 777久久人妻少妇嫩草av网站| 国产福利在线免费观看视频| 免费日韩欧美在线观看| 午夜福利一区二区在线看| 国产亚洲欧美精品永久| 一区福利在线观看| 一边亲一边摸免费视频| 日韩精品免费视频一区二区三区| 国产男女内射视频| 久久精品熟女亚洲av麻豆精品| 国产精品亚洲av一区麻豆| av天堂在线播放| 久久这里只有精品19| 亚洲 国产 在线| 婷婷色综合大香蕉| 好男人电影高清在线观看| 别揉我奶头~嗯~啊~动态视频 | 精品人妻在线不人妻| 亚洲综合色网址| tube8黄色片| 各种免费的搞黄视频| 极品少妇高潮喷水抽搐| 欧美黄色淫秽网站| 国产免费又黄又爽又色| 国产精品三级大全| 久久毛片免费看一区二区三区| 亚洲欧洲日产国产| 国产一卡二卡三卡精品| 青青草视频在线视频观看| 深夜精品福利| 日日夜夜操网爽| 九草在线视频观看| 99国产精品99久久久久| 色网站视频免费| 天天躁日日躁夜夜躁夜夜| 黑丝袜美女国产一区| 咕卡用的链子| 免费女性裸体啪啪无遮挡网站| 一本久久精品| 国产成人精品久久二区二区91| 777久久人妻少妇嫩草av网站| 麻豆av在线久日| 丰满饥渴人妻一区二区三| 麻豆国产av国片精品| 18禁观看日本| 国产97色在线日韩免费| 久久精品熟女亚洲av麻豆精品| 欧美另类一区| 精品人妻熟女毛片av久久网站| 亚洲 欧美一区二区三区| 久久精品人人爽人人爽视色| 欧美人与性动交α欧美软件| 久久免费观看电影| 国产成人欧美| 亚洲免费av在线视频| 在线观看一区二区三区激情| √禁漫天堂资源中文www| 久久天堂一区二区三区四区| 国产av国产精品国产| 老司机影院毛片| 亚洲欧美清纯卡通| 亚洲激情五月婷婷啪啪| 亚洲国产欧美日韩在线播放| e午夜精品久久久久久久| 国产精品麻豆人妻色哟哟久久| 日本猛色少妇xxxxx猛交久久| 18禁国产床啪视频网站| 国产1区2区3区精品| 亚洲九九香蕉| a 毛片基地| 国产精品成人在线| 亚洲av在线观看美女高潮| av又黄又爽大尺度在线免费看| 80岁老熟妇乱子伦牲交| 国产99久久九九免费精品| 国产高清不卡午夜福利| 欧美日韩av久久| 国产黄频视频在线观看| 狠狠精品人妻久久久久久综合| 亚洲精品国产区一区二| 久久性视频一级片| 亚洲av在线观看美女高潮| 久久久久国产一级毛片高清牌| 欧美变态另类bdsm刘玥| 久久99热这里只频精品6学生| 波野结衣二区三区在线| 下体分泌物呈黄色| 波多野结衣一区麻豆| 精品人妻1区二区| 日韩av免费高清视频| 亚洲av美国av| 亚洲国产最新在线播放| 日韩av不卡免费在线播放| 婷婷成人精品国产| 亚洲av男天堂| 亚洲人成网站在线观看播放| 高清不卡的av网站| 久久精品国产亚洲av高清一级| 亚洲综合色网址| 热99久久久久精品小说推荐| 国产免费视频播放在线视频| 高清黄色对白视频在线免费看| 午夜两性在线视频| 人人妻人人澡人人看| 一个人免费看片子| 国产高清国产精品国产三级| 成人亚洲精品一区在线观看| 9热在线视频观看99| 国产精品一国产av| 波多野结衣av一区二区av| 欧美中文综合在线视频| 黄色视频不卡| 国产老妇伦熟女老妇高清| 国产黄频视频在线观看| 涩涩av久久男人的天堂| 成人影院久久| 老司机影院毛片| 97人妻天天添夜夜摸| 久9热在线精品视频| 波多野结衣av一区二区av| 成人影院久久| 老司机影院毛片| 水蜜桃什么品种好| 99热国产这里只有精品6| 久久国产精品男人的天堂亚洲| 国产高清国产精品国产三级| 免费高清在线观看日韩| 777久久人妻少妇嫩草av网站| 国产熟女午夜一区二区三区| 亚洲国产毛片av蜜桃av| 亚洲av电影在线观看一区二区三区| 人人妻人人澡人人爽人人夜夜| 欧美精品啪啪一区二区三区 | 欧美黄色淫秽网站| 久久亚洲国产成人精品v| 精品国产国语对白av| a级毛片在线看网站| 欧美成人精品欧美一级黄| 亚洲免费av在线视频| 国产不卡av网站在线观看| 亚洲中文日韩欧美视频| 久久免费观看电影| 韩国精品一区二区三区| 欧美少妇被猛烈插入视频| 女性生殖器流出的白浆| 国产精品99久久99久久久不卡| 亚洲人成电影免费在线| 成年动漫av网址| 一区二区三区四区激情视频| 男女下面插进去视频免费观看| 老司机靠b影院| 国产精品亚洲av一区麻豆| 亚洲,欧美精品.| 天天躁日日躁夜夜躁夜夜| 这个男人来自地球电影免费观看| 亚洲精品一区蜜桃| 国产精品一区二区在线不卡| 国产av精品麻豆| 日韩一区二区三区影片| 成人国产av品久久久| 久久久久久久国产电影| 欧美日韩亚洲国产一区二区在线观看 | 制服诱惑二区| 国产日韩欧美在线精品| 亚洲欧美日韩高清在线视频 | 成人亚洲欧美一区二区av| 久久人人爽人人片av| 色婷婷av一区二区三区视频| 欧美激情 高清一区二区三区| 美国免费a级毛片| 国产黄频视频在线观看| 国产精品久久久av美女十八| 免费在线观看日本一区| 精品福利永久在线观看| www.自偷自拍.com| 亚洲精品国产色婷婷电影| 色网站视频免费| 一区二区三区激情视频| 国产一区亚洲一区在线观看| 精品人妻熟女毛片av久久网站| www日本在线高清视频| 亚洲熟女精品中文字幕| 国产精品香港三级国产av潘金莲 | 欧美精品一区二区大全| 韩国精品一区二区三区| 国产精品一区二区在线不卡| 少妇裸体淫交视频免费看高清 | 国产xxxxx性猛交| 脱女人内裤的视频| 母亲3免费完整高清在线观看| 亚洲av成人不卡在线观看播放网 | 亚洲精品国产区一区二| 高清av免费在线| 成人亚洲精品一区在线观看| 99热全是精品| 各种免费的搞黄视频| 亚洲国产日韩一区二区| 亚洲精品一卡2卡三卡4卡5卡 | 国产视频首页在线观看| 香蕉丝袜av| 免费在线观看影片大全网站 | 欧美激情高清一区二区三区| 国产精品久久久久久精品古装| 午夜激情久久久久久久| 看十八女毛片水多多多| 国产日韩欧美视频二区| 日韩大片免费观看网站| 日韩,欧美,国产一区二区三区| a级毛片黄视频| 久久久精品国产亚洲av高清涩受| 宅男免费午夜| 色婷婷久久久亚洲欧美| 少妇人妻 视频| 亚洲第一青青草原| 欧美 亚洲 国产 日韩一| 在线观看免费午夜福利视频| 国产又色又爽无遮挡免| 免费在线观看视频国产中文字幕亚洲 | 女人精品久久久久毛片| 一二三四社区在线视频社区8| 精品福利观看| 国产97色在线日韩免费| 国产一区二区三区综合在线观看| 国产精品一区二区精品视频观看| 国产精品秋霞免费鲁丝片| 亚洲视频免费观看视频| 热99国产精品久久久久久7| 国产熟女午夜一区二区三区| 午夜福利一区二区在线看| 亚洲欧美清纯卡通| 99热国产这里只有精品6| 1024香蕉在线观看| 欧美精品啪啪一区二区三区 | 欧美国产精品一级二级三级| 国产男女超爽视频在线观看| 国产精品久久久久久精品古装| √禁漫天堂资源中文www| 久久中文字幕一级| 中文字幕亚洲精品专区| 久久午夜综合久久蜜桃| 国产成人精品久久二区二区91| 国产在线一区二区三区精| 国产熟女午夜一区二区三区| 国产av一区二区精品久久| 高清欧美精品videossex|