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

    噴射結(jié)構(gòu)對(duì)充液圓管內(nèi)氣幕特性影響的數(shù)值分析

    2016-05-08 08:21:01周良梁余永剛曹永杰
    含能材料 2016年7期
    關(guān)鍵詞:噴孔側(cè)壁斜面

    周良梁, 余永剛, 曹永杰

    (1. 南京理工大學(xué)能源與動(dòng)力工程學(xué)院, 江蘇 南京 210094; 2. 西北機(jī)電工程研究所, 陜西 咸陽 7120991)

    1 引 言

    常規(guī)水下槍炮發(fā)射通常采用淹沒式發(fā)射或密封式發(fā)射。淹沒式直接將發(fā)射管淹沒在水中,發(fā)射結(jié)構(gòu)簡(jiǎn)單,但發(fā)射過程中為抑制過高膛壓的出現(xiàn),通常發(fā)射初速較低[1-2]。而密封式發(fā)射通過管口密封,將發(fā)射管內(nèi)部與液體工質(zhì)隔絕,使得發(fā)射時(shí)管內(nèi)處于氣體環(huán)境中,能夠獲得較高的發(fā)射初速[3],但該種發(fā)射系統(tǒng)結(jié)構(gòu)較為復(fù)雜。氣幕式發(fā)射作為一種新型的導(dǎo)彈水下發(fā)射方式,可以通過簡(jiǎn)單的內(nèi)部結(jié)構(gòu),利用發(fā)動(dòng)機(jī)產(chǎn)生的廢氣在發(fā)射管外部為導(dǎo)彈創(chuàng)造一條氣體通道,減小發(fā)射過程中附加質(zhì)量與阻力[4]。

    對(duì)于氣幕式發(fā)射過程,Yagla等人[5]通過實(shí)驗(yàn)和數(shù)值模擬獲得了水下氣幕形態(tài),并對(duì)燃?xì)馍淞髋c液體作用界面穩(wěn)定性特征進(jìn)行了分析。曹嘉怡[6]對(duì)氣幕式發(fā)射彈體出筒過程進(jìn)行模擬,獲得了燃?xì)膺吔缫约皟?nèi)部流場(chǎng)隨時(shí)間變化特性。王亞東[7]等人采用可壓縮氣體模型,對(duì)氣幕式發(fā)射過程中氣體彈性對(duì)彈體運(yùn)動(dòng)的影響進(jìn)行了數(shù)值分析。目前國內(nèi)外有關(guān)水下氣幕式發(fā)射的研究主要集中在導(dǎo)彈水下發(fā)射領(lǐng)域,通過實(shí)驗(yàn)研究與數(shù)值分析,獲得導(dǎo)彈水下氣幕式發(fā)射過程流場(chǎng)特性,但尚未見到有關(guān)火炮水下氣幕式發(fā)射的研究報(bào)道。

    火炮水下氣幕式發(fā)射通過多股燃?xì)馍淞魉聟R聚生成氣幕,有關(guān)氣體射流與水相互作用問題相關(guān)學(xué)者已進(jìn)行了一些研究。Dai[8]對(duì)水下超聲速氣流流動(dòng)特性進(jìn)行了實(shí)驗(yàn)觀測(cè)。施紅輝等[9]對(duì)超音速射流脹鼓和回?fù)籼匦赃M(jìn)行了實(shí)驗(yàn)研究,結(jié)果表明脹鼓頻率越大,回?fù)纛l率越大。齊麗婷等[10]通過實(shí)驗(yàn)與數(shù)值分析對(duì)燃?xì)馍淞髟谝后w中擴(kuò)展的二維問題進(jìn)行了研究。曹嘉怡等[11]研究分析了水下超聲速噴管燃?xì)馍淞鲃?dòng)力學(xué)特性。莽珊珊[12]通過實(shí)驗(yàn),觀測(cè)了單股高溫燃?xì)馍淞髟诔錆M液體的圓柱形及漸擴(kuò)型容器中擴(kuò)展形態(tài)。Weiland等[13]研究了圓形射流入水過程中氣液邊界動(dòng)力學(xué)特性與氣液作用穩(wěn)定性的關(guān)系。Yu[14]對(duì)雙股射流在二維平面觀察室中摻混過程進(jìn)行了探討。Xue[15]研究了雙股射流在圓柱形充液室中摻混過程,并分析了結(jié)構(gòu)參數(shù)變化對(duì)射流發(fā)展的影響。Graaf等[16]對(duì)4股射流入水后四股射流湍流匯聚過程動(dòng)力學(xué)特性進(jìn)行了分析。目前,對(duì)氣體射流與水相互作用的研究主要側(cè)重于流場(chǎng)特性及射流匯聚過程,尚未涉及多股氣體射流生成管內(nèi)氣幕的研究。

    基于氣幕式發(fā)射原理,本研究擬采用多股燃?xì)馍淞髟诠軆?nèi)生成氣幕,為水下火炮彈丸發(fā)射提供氣體通道。為了解多股燃?xì)馍淞髟谏砉軆?nèi)生成氣幕過程,設(shè)計(jì)了多股燃?xì)馍淞髋潘M試驗(yàn)裝置,針對(duì)實(shí)驗(yàn)結(jié)果采用數(shù)值模擬,分析噴射結(jié)構(gòu)對(duì)氣幕生成過程的影響。

    2 數(shù)學(xué)物理模型

    2.1 物理模型

    針對(duì)高溫高壓燃?xì)馍淞髟谝后w工質(zhì)中擴(kuò)展特性,對(duì)多股燃?xì)馍淞鲾U(kuò)展過程進(jìn)行如下簡(jiǎn)化處理: (1)燃?xì)鉃槔硐肟蓧簹怏w,燃?xì)馍淞鬟M(jìn)入觀察室后不再發(fā)生化學(xué)反應(yīng),噴射過程中燃?xì)饨M分保持不變; (2)采用k-ε模型模擬流場(chǎng)中氣液湍流作用; (3)考慮到高溫燃?xì)馀c液體作用時(shí)間短暫(小于10 ms),可忽略液體工質(zhì)的汽化的影響; (4)物性參數(shù)取平均值。

    2.2 數(shù)學(xué)模型

    根據(jù)物理模型,采用Navies-Stokes方程組[17]、氣體狀態(tài)方程[17]與k-ε湍流方程組[18],建立多股燃?xì)馍淞魉聰U(kuò)展的三維非穩(wěn)態(tài)計(jì)算模型。

    (1)連續(xù)性方程:

    (1)

    式中,ρ為總密度,kg·m-3;v為速度矢量,m·s-1;S為質(zhì)量源項(xiàng),kg·m-3·s-1,由于射流發(fā)展過程中不考慮化學(xué)反應(yīng)與液體蒸發(fā),因此S=0; 為捕捉氣液相互界面,采用VOF模型[17]計(jì)算氣液組分,燃?xì)饨M分連續(xù)性方程:

    (2)

    式中,α1為氣體的體積分?jǐn)?shù),作為計(jì)算主項(xiàng),水的體積分?jǐn)?shù)α2=1-α1。

    (2)動(dòng)量方程:

    (3)

    式中,p為壓力,Pa;g為重力矢量,m·s-2。

    (3)能量方程:

    (4)

    式中,Ei為第i相流體總能,J;keff為混合物有效傳熱系數(shù),W·m-2·K-1;keff=k1α1+k2α2,式中k1為氣相傳熱系數(shù),W·m-2·K-1;k2為液相傳熱系數(shù),W·m-2·K-1;T為溫度,K。

    (4)氣體狀態(tài)方程:

    p=ρ1RgT

    (5)

    式中,Rg為氣相組份所對(duì)應(yīng)的氣體常數(shù),J·mol-1·K-1,計(jì)算過程入口溫度取為燃?xì)鉁囟绕骄礣=2000 K。

    (5)k-ε湍流方程

    (6)

    (7)

    式中,μt為湍動(dòng)粘度系數(shù),Pa·s;Gk、Gb、YM分別為速度梯度、浮力以及脈動(dòng)引起的湍流動(dòng)能,單位均為J;C1ε=1.44、C2ε=1.92、C3ε=0.09為經(jīng)驗(yàn)常數(shù),σk=1.0、σε=1.3分別為湍動(dòng)能與耗散率對(duì)應(yīng)的Prandtl數(shù)。

    2.3 計(jì)算域及邊界條件

    實(shí)驗(yàn)裝置及噴頭結(jié)構(gòu)如圖1、圖2所示。該裝置主要由觀察室、噴頭與燃?xì)獍l(fā)生器三部分組成。觀察室由透明有機(jī)玻璃制成,實(shí)驗(yàn)過程中內(nèi)部裝滿液體工質(zhì),觀察室內(nèi)徑為2 cm,總長度為1 m。噴頭內(nèi)部有密封膜片,防止觀察室內(nèi)液體灌入。噴頭頂部設(shè)計(jì)9個(gè)對(duì)稱分布的噴孔,即: 1個(gè)中心噴孔,4個(gè)斜面噴孔與4個(gè)側(cè)壁矩形噴孔,噴孔分布如圖2所示,圖中A1-A1為噴頭對(duì)稱截面。觀察室、噴頭與壓力傳感器分別固定在燃?xì)獍l(fā)射器上。燃?xì)獍l(fā)生器內(nèi)裝有速燃火藥。該裝置的工作原理為:首先由脈沖電源放電點(diǎn)燃燃?xì)獍l(fā)生器內(nèi)速燃火藥,發(fā)生器內(nèi)壓力迅速上升,當(dāng)壓力到達(dá)密封膜片閾值時(shí),高溫高壓燃?xì)鉀_破噴頭內(nèi)密封膜片,進(jìn)入充液觀察室內(nèi)。在觀察室中,多股燃?xì)馍淞髋c液體工質(zhì)發(fā)生氣液摻混,生成管內(nèi)氣幕,推動(dòng)噴頭上方液體運(yùn)動(dòng)。

    圖1 實(shí)驗(yàn)裝置結(jié)構(gòu)

    1—觀察室, 2—噴頭, 3—壓力傳感器, 4—燃?xì)獍l(fā)生器

    Fig.1 Structure of experiment device

    1—observation chamber, 2—sprayer, 3—pressure transducer, 4—gas generator

    由圖2可以看出,噴頭表面噴孔對(duì)稱分布,可以認(rèn)為四個(gè)側(cè)壁噴孔與四個(gè)斜面噴孔分布具有相同的噴射條件,而中心射流沿中心軸線擴(kuò)展,因此選取觀察室1/4區(qū)域作為計(jì)算區(qū)域。計(jì)算區(qū)域及邊界示意圖如圖3所示。

    a.top projection b. section view

    圖2 噴頭結(jié)構(gòu)示意圖

    1—中心噴孔, 2—斜面噴孔, 3—側(cè)壁噴孔, 4—密封膜片

    Fig.2 Structure of sprayer

    1—center nozzle, 2—slant nozzle, 3—side wall nozzle,4—sealing film

    圖3 計(jì)算區(qū)域示意圖

    Fig.3 Diagram of computational domain

    壓力入口取為密封膜片處,入口壓力根據(jù)實(shí)驗(yàn)測(cè)量結(jié)果進(jìn)行擬合; 壓力出口取為與大氣相通的觀察室頂部,出口條件取為大氣條件。計(jì)算采用壓力基求解方式,密度、動(dòng)量、能量與湍流項(xiàng)采用一階迎風(fēng)差分格式,壓力項(xiàng)采用PRESTO!插值格式,壓力與速度通過PISO方法耦合。

    為了分析比較噴射結(jié)構(gòu)對(duì)射流擴(kuò)展過程的影響,針對(duì)兩種不同噴射結(jié)構(gòu)的噴頭進(jìn)行了數(shù)值模擬,噴頭結(jié)構(gòu)尺寸如表1所示。兩種噴頭側(cè)壁噴孔尺寸相同,B噴頭在A噴頭的基礎(chǔ)上將斜面噴孔增大到2 mm。

    表1 噴頭結(jié)構(gòu)尺寸

    Table 1 Structure dimensions of the sprayers

    typediameterofcenternozzleΦ1/mmdiameterofobliquenozzlesΦ2/mmsizeofsidewallnozzles/mmA2.01.53.0×1.0B2.02.03.0×1.0

    計(jì)算初始時(shí),觀察室內(nèi)充滿液體工質(zhì)(α1=0),而噴頭內(nèi)部充滿氣體(α1=1),噴頭內(nèi)氣體與觀察室中液體相通,而觀察室中液體與外部大氣相通,初始參數(shù)為環(huán)境參數(shù)。計(jì)算開始時(shí)刻取為膜片被沖破的瞬間,高溫高壓燃?xì)鈴膲毫θ肟谔幜魅搿?/p>

    3 數(shù)值模擬與實(shí)驗(yàn)結(jié)果對(duì)比

    針對(duì)A型噴頭的實(shí)驗(yàn)工況,模擬了多股燃?xì)馍淞髟诔湟菏抑邪l(fā)展過程。壓力入口條件根據(jù)實(shí)驗(yàn)結(jié)果近似得到: 壓力p=(0.45+2.45e-t/0.0085) MPa,式中t為時(shí)間。圖4a為實(shí)驗(yàn)中高速錄像拍攝的A型噴頭多股燃?xì)馍淞鲾U(kuò)展序列過程。圖4b為對(duì)應(yīng)時(shí)刻數(shù)值模擬結(jié)果組分分布投影圖。對(duì)比實(shí)驗(yàn)照片與模擬圖像,可以看出模擬獲得的射流形態(tài)與實(shí)驗(yàn)結(jié)果基本一致。

    a. experiments b. simulations

    圖4 A型噴頭多股射流擴(kuò)展序列圖(pmax=2.9 MPa)

    Fig.4 The expansion sequence of A sprayer multi jets(pmax=2.9 MPa)

    通過讀取圖4b中多股燃?xì)馍淞黜敹说臍庖航缑娴奈恢米鴺?biāo),可以得到多股燃?xì)馍淞鬏S向擴(kuò)展位移變化過程,將模擬結(jié)果與實(shí)驗(yàn)值進(jìn)行比較,如圖5所示。從圖5中可以看出在射流發(fā)展初始階段,相對(duì)誤差較大,這主要是因?yàn)楦咚贁z像記錄時(shí)間間隔為0.5 ms,初始時(shí)刻捕捉存在一定誤差,實(shí)驗(yàn)結(jié)果中0.5 ms時(shí)讀取的初始位移存在初始誤差。但在后期兩者逐漸趨近,在8 ms時(shí)相對(duì)誤差接近5%,計(jì)算結(jié)果與實(shí)測(cè)值吻合較好。

    圖5 多股射流軸向擴(kuò)展位移計(jì)算值與實(shí)測(cè)值的對(duì)比

    Fig.5 Comparison between simulated axial displacement of the multi gas jets and experiment result

    4 計(jì)算結(jié)果與分析

    在以上工作的基礎(chǔ)上,針對(duì)A型噴頭多股燃?xì)馍淞鲾U(kuò)展特性進(jìn)行數(shù)值分析。采用與A型噴頭相同壓力工況,針對(duì)B型進(jìn)行數(shù)值預(yù)測(cè),對(duì)比分析噴射結(jié)構(gòu)參數(shù)變化對(duì)氣幕生成特性的影響。

    4.1 多股燃?xì)馍淞鲾U(kuò)展特性分析4.1.1 組分分布

    針對(duì)多股燃?xì)馍淞魉诘膶?duì)稱面A1-A1進(jìn)行流場(chǎng)分析,得到A型噴頭多股燃?xì)馍淞鲄R聚演化過程,如圖6所示。

    圖6 射流在對(duì)稱截面上的氣液組份分布

    Fig.6 The phase distribution of the multi jets in the cross section

    從圖6射流的擴(kuò)展過程可以看出,1 ms之前中心與斜面射流Taylor空腔分別沿著噴孔方向獨(dú)立擴(kuò)展。在1.5 ms時(shí)斜面射流撞擊到壁面,Taylor空腔擴(kuò)展方向發(fā)生改變,向軸向偏轉(zhuǎn),同時(shí)斜面射流界面上氣液湍流摻混作用增強(qiáng),界面出現(xiàn)卷吸; 斜面射流徑向擴(kuò)展抑制了中心射流尾部徑向發(fā)展,中心射流尾部收縮。在 2 ms時(shí)燃?xì)忾_始從側(cè)壁矩形噴孔中噴出,推動(dòng)斜面射流向軸向擴(kuò)展; 中心射流頂部由于Taylor不穩(wěn)定性的作用,射流中心凹陷,出現(xiàn)分叉變形,隨著中心射流發(fā)展,頂部分叉現(xiàn)象進(jìn)一步加劇。在2.5 ms時(shí),側(cè)壁射流開始匯入斜面射流,匯聚界面氣液湍流摻混增強(qiáng),3 ms時(shí)側(cè)壁射流完全匯入斜面射流生成貼壁發(fā)展的側(cè)面射流。隨著側(cè)面射流徑向發(fā)展,側(cè)面射流逐漸向中心射流靠攏。在6 ms時(shí),側(cè)面射流與中心射流開始匯聚。7~8 ms,側(cè)面射流繼續(xù)向中心擴(kuò)展,多股射流在觀察室內(nèi)生成與觀察室橫截面積接近的管內(nèi)氣幕。

    4.1.2 壓力分布

    A型噴頭多股燃?xì)馍淞髟趯?duì)稱截面上壓力演化過程如圖7所示。

    圖7 射流在對(duì)稱截面上的壓力分布

    Fig.7 The pressure distribution of the multi jets in the cross section

    由圖7可以看出,射流發(fā)展過程中,高壓區(qū)主要集中在噴頭內(nèi)。燃?xì)馍淞鲝膰娍讎姵龊笤谝后w中形成Taylor空腔,燃?xì)庋杆倥蛎?,壓力沿程降低? ms時(shí),中心射流與斜面射流Taylor空腔頂端受到液體的阻擋,射流頂部產(chǎn)生壓縮,出現(xiàn)局部高壓。由于斜面射流擴(kuò)展受阻,2 ms時(shí)斜面射流Taylor空腔內(nèi)部壓力上升; 而中心射流尾部因?yàn)槭湛s,在噴孔上方產(chǎn)生局部高壓。3 ms時(shí)隨著斜面射流與側(cè)壁射流的匯聚,交匯區(qū)因氣液摻混壓力上升,在壁面x=5 mm附近出現(xiàn)一個(gè)貼壁局部高壓區(qū),同時(shí)側(cè)面射流頂部軸向擴(kuò)展受到頂部液體工質(zhì)的阻礙,在頂部產(chǎn)生貼壁壓縮區(qū),壓力上升。隨著射流沿軸向的發(fā)展,中心射流與斜面射流的摻混,5 ms時(shí)中心射流摻混區(qū)內(nèi)部壓力上升,高壓區(qū)隨摻混區(qū)的軸向發(fā)展移動(dòng)。中心射流頂部氣體能量逐漸耗散,中心射流頂部壓縮作用減弱,頂部高壓區(qū)消失。

    4.1.3 速度分布

    針對(duì)多股燃?xì)馍淞髟趯?duì)稱面上的流動(dòng)特性分析,得到A型噴頭射流在對(duì)稱面上的軸向速度與流線演化過程,如圖8,圖9所示。

    圖8 射流在對(duì)稱面上的軸向速度分布

    Fig.8 The velocity distribution of the multi jets in the cross section

    圖9 射流在對(duì)稱面上的流線演化圖

    Fig.9 The vortex evolution process in the cross section

    結(jié)合圖8、圖9可以看出: 多股燃?xì)馍淞髟诎l(fā)展過程中,高速運(yùn)動(dòng)區(qū)主要集中在噴孔附近,且噴孔處速度梯度較大。結(jié)合圖6可以發(fā)現(xiàn): 射流發(fā)展初期在t=1 ms時(shí),高溫燃?xì)饧铀倭鞒鰢娍?,受到液體阻滯,在噴孔附近產(chǎn)生回流漩渦。在t=2 ms時(shí),隨著側(cè)壁射流出現(xiàn),側(cè)壁噴孔上方區(qū)域速度上升,在側(cè)壁噴孔上方出現(xiàn)一個(gè)低速流動(dòng)區(qū)域; Taylor空腔頂部回流減弱,回流區(qū)消失。在t=3 ms時(shí),側(cè)壁射流與斜面射流匯聚生成側(cè)面射流; 與斜面射流相比,側(cè)面射流軸向擴(kuò)展能力增強(qiáng),側(cè)面射流頂部氣液作用加劇,側(cè)面射流Taylor空腔頂部出現(xiàn)回流區(qū),抑制了中心射流與側(cè)面射流之間殘余液體的軸向運(yùn)動(dòng),在兩股射流之間形成低速區(qū)。隨著射流的軸向擴(kuò)展,回流區(qū)增大,中心也隨之向下游移動(dòng)。隨著側(cè)面射流的徑向發(fā)展,t=8 ms時(shí)中心射流Taylor空腔與側(cè)面射流Taylor空腔匯聚,氣液界面處的湍流摻混在空腔內(nèi)部形成低速旋渦區(qū),射流間低速區(qū)域徑向分布范圍進(jìn)一步增大,斜面射流噴孔上方回流區(qū)重新出現(xiàn)。

    4.2 噴射結(jié)構(gòu)對(duì)氣幕生成特性的影響

    針對(duì)2 mm斜面噴孔的B型噴頭射流進(jìn)行流場(chǎng)分析,得到B型噴頭射流在對(duì)稱面A1-A1上的組分分布如圖10所示,圖11給出了對(duì)應(yīng)時(shí)刻壓力分布。

    對(duì)比A型噴頭(斜面孔1.5 mm)圖6與圖10中組分分布,可以發(fā)現(xiàn): 斜面噴孔直徑增大后,斜面射流擴(kuò)展尺寸增大,使側(cè)壁射流發(fā)展受阻滯后,直至3 ms時(shí)才出現(xiàn),而A型噴頭側(cè)壁射流在2 ms時(shí)就出現(xiàn)。另外,斜面噴孔直徑增大,導(dǎo)致斜面射流與側(cè)壁射流匯聚生成的側(cè)面射流軸向擴(kuò)展能力增強(qiáng),側(cè)面射流頂部軸向位移增大。同時(shí)中心射流頂部氣液湍流摻混過程中Taylor不穩(wěn)定性減弱,直至5 ms時(shí),中心射流頂部才出現(xiàn)分叉變形,而A型噴頭中心射流在2 ms時(shí)就出現(xiàn)分叉變形。

    圖10 射流在對(duì)稱截面上的氣液組份分布 (B型噴頭)

    Fig.10 The phase distribution of the multi jets in the cross section (B sprayer)

    圖11 射流在對(duì)稱截面上的壓力分布(B型噴頭)

    Fig.11 The pressure distribution of the multi jets in the cross section (B sprayer)

    由圖11,可以發(fā)現(xiàn),射流發(fā)展初期,B型噴頭因斜面射流Taylor空腔內(nèi)部氣體聚集量大,產(chǎn)生的高壓區(qū)范圍寬,導(dǎo)致側(cè)壁噴孔上方也形成高壓區(qū),阻滯了側(cè)壁射流的發(fā)展。隨著側(cè)壁射流與斜面射流的匯聚,進(jìn)一步強(qiáng)化了高壓區(qū)的存在,而A型噴頭相應(yīng)位置的高壓區(qū)范圍較小。在射流發(fā)展后期,B型噴頭形成的氣幕頂部壓縮作用較強(qiáng),頂部始終有一個(gè)高壓區(qū)存在,而A型噴頭形成的氣幕頂部在5 ms后高壓區(qū)就消失了。

    結(jié)合射流頂部位移可以獲得不同斜面噴孔下氣幕頂部軸向擴(kuò)展速度,如圖12所示。由圖可見,3 ms之前,A型與B型噴頭射流頂部軸向擴(kuò)展速度呈波動(dòng)下降趨勢(shì)。在3 ms之后,B型噴頭射流平穩(wěn)上升; 而A型噴頭射流在3.5 ms之后呈波動(dòng)上升趨勢(shì); 射流發(fā)展后期,B型噴頭射流頂部軸向速度超過A型噴頭射流。這主要是由于B型噴頭氣幕前端有一個(gè)穩(wěn)定的高壓區(qū),而A型噴頭氣幕前端高壓區(qū)存在時(shí)間短,從而易受到上方液體擾動(dòng),導(dǎo)致速度波動(dòng)較大。

    圖12 不同斜面噴孔下氣幕頂部軸向擴(kuò)展速度

    Fig.12 Axial expansion velocity of the gas-curtain with different slant nozzles

    圖13為B型噴頭在t=8 ms時(shí)流場(chǎng)運(yùn)動(dòng)特性分布。對(duì)比圖8、圖9中A型噴頭8 ms時(shí)速度與流線分布,可以看出: 增大斜面噴孔直徑,側(cè)面射流回流作用增強(qiáng),斜面射流回流區(qū)向上方移動(dòng),低速區(qū)核心區(qū)域出現(xiàn)在斜面射流回流區(qū)中,而A型噴頭對(duì)應(yīng)低速核心區(qū)出現(xiàn)在側(cè)面射流回流區(qū)。B型噴頭側(cè)面射流回流區(qū)中心軸向位置基本不變,但徑向位置向中心軸靠近; 而中心射流內(nèi)渦流區(qū)徑向位置不變,但軸向位置向下移動(dòng)之x=13 mm處,而A型噴頭中對(duì)應(yīng)位置為x=17 mm,同時(shí)中心射流渦流區(qū)范圍也隨之減小。

    圖14為x=200 mm處液體工質(zhì)運(yùn)動(dòng)特性曲線。從圖14中可以看出: B型噴頭氣幕上方液體速度始終大于A型噴頭上方液體。在2.5 ms至3.5 ms期間,兩個(gè)噴頭上方液體的加速度基本相同,其他時(shí)刻B型噴頭上液體的加速度大于A型噴頭。在0 ms時(shí),使用A型噴頭時(shí)液體瞬時(shí)加速度為1940 m/s2; 使用B型噴頭時(shí)液體瞬時(shí)加速度為2200 m/s2。在8 ms時(shí),A型噴頭上方液體運(yùn)動(dòng)速度為8.26 m·s-1,瞬時(shí)加速度為660 m/s2; 而B型噴頭上方液體瞬時(shí)加速度增大到800 m/s2,運(yùn)動(dòng)速度提高到9.4 m·s-1。增大斜面噴孔,提高了射流整體排水性能。

    圖13 B型噴頭流場(chǎng)運(yùn)動(dòng)特性(t=8 ms)

    Fig.13 Kinetic characteristic of B sprayer flow(t=8 ms)

    a.v-tcurveb.a-tcurve

    圖14 不同斜面噴孔下液體工質(zhì)運(yùn)動(dòng)特性(x=200 mm)

    Fig.14 Kinetic characteristic of liquid medium with different slant nozzles(x=200 mm)

    5 結(jié) 論

    (1)建立了多股燃?xì)馍淞魉聰U(kuò)展的三維非穩(wěn)態(tài)數(shù)理模型,計(jì)算了多股燃?xì)馍淞髟诔湟簣A管內(nèi)通過湍流摻混生成管氣幕的過程,對(duì)比模擬結(jié)果與實(shí)驗(yàn)中氣液摻混照片,可以發(fā)現(xiàn)兩者形態(tài)基本一致,且計(jì)算獲得的射流頂端軸向擴(kuò)展位移與實(shí)測(cè)結(jié)果基本吻合,說明計(jì)算模型是合理的。

    (2)多股燃?xì)馍淞鲝膰娍讎姵龊?,在液體中形成多個(gè)Taylor空腔,Taylor空腔頂部受到液體阻滯,射流頂部壓力升高; 斜面射流撞壁與側(cè)壁射流匯合后生成側(cè)面射流,在射流發(fā)展后期,隨著中心射流與側(cè)面射流的匯聚,氣液湍流摻混劇烈,在匯聚界面產(chǎn)生低速渦流區(qū)域,頂部高壓區(qū)消失,但中心射流與側(cè)面射流匯聚區(qū)域產(chǎn)生局部高壓; 中心射流頂部由于Taylor不穩(wěn)定性作用,出現(xiàn)分叉變形。

    (3)噴頭斜面噴孔由1.5 mm增大到2 mm時(shí),斜面射流與側(cè)壁射流匯聚區(qū)域壓力升高,側(cè)壁射流發(fā)展受到抑制,斜面射流回流作用增強(qiáng),低速核心區(qū)出現(xiàn)在斜面射流回流核心區(qū)域; 側(cè)面射流整體軸向擴(kuò)展能力增強(qiáng),射流發(fā)展后期多股射流整體排水性能提升,8 ms時(shí)B型噴頭氣幕上方液體運(yùn)動(dòng)速度增大了1.14 m·s-1。

    參考文獻(xiàn):

    [1] 孔德仁, 王昌明, 柳光遼, 等. 水下槍械內(nèi)彈道基本方程組及其設(shè)計(jì)[J]. 南京理工大學(xué)學(xué)報(bào), 1999, 23(3): 9-12.

    KONG De-ren, WANG Chang-ming, LIU Guang-liao, et al. Fundamental equations of underwater gun′s interior trajectory and it′s design[J].JournalofNanjingUniversityofScienceandTechnology, 1999, 23(3): 9-12.

    [2] 王昌明, 陳(舟疌), 柳光遼, 等. 相關(guān)分析在水下槍械內(nèi)彈道計(jì)算中的應(yīng)用[J]. 南京理工大學(xué)學(xué)報(bào), 2005, 29(3): 260-263.

    WANG Chang-ming, CHEN Jie, LIU Guang-liao, et al. Application of correlation analysis used in interior ballistics calculation for underwater gun[J].JournalofNanjingUniversityofScienceandTechnology, 2005, 29(3): 260-263.

    [3] 劉育平, 李金新, 楊臻, 等. 水下炮內(nèi)彈道分析與數(shù)值仿真[J]. 火炮發(fā)射與控制學(xué)報(bào), 2007(4): 30-33.

    LIU Yu-ping, LI Jin-xin, YANG Zhen, et al. Interior ballistics analysis and numerical simulation of underwater gun[J].JournalofGunLuanch&Control, 2007(4): 30-33.

    [4] Yagla J, John B, Koski S. Launch Dynamics Environment of a Water Piercing Missile Launcher[C]∥24th International Symposium on Ballistics, New Orleans, Louisiana, 2008: 1-17.

    [5] Yagla J, Deloach K, Weiland C. Water Piercing Missile Launcher[C]∥22th International Symposium on Ballistics, Adelaide, AU, 2004: 0419-0423.

    [6] 曹嘉怡, 魯傳敬, 陳鑫, 等. 導(dǎo)彈水下熱發(fā)射出筒過程流動(dòng)特性[J]. 固體火箭技術(shù), 2011, 34(3): 281-284.

    CAO Jia-yi, LU Chuan-jing, CHEN Xin, et al. Flow behavior during uncorking process of an underwater hot-launched missile[J].JournalofSolidRocketTechnology, 2011, 34(3): 281-284.

    [7] 王亞東, 袁緒龍, 張宇文, 等. 氣泡彈性對(duì)同心筒水下發(fā)射影響研究[J]. 計(jì)算力學(xué)學(xué)報(bào), 2013, 30(2): 313-318.

    WANG Ya-dong, YUAN Xu-long, ZHANG Yu-wen, et al. On the effect of bubble elasticity to sub-launch process using concentric canister launcher[J].ChineseJournalofComputationalMechanics, 2013, 30(2): 313-318.

    [ 8] DAI Z, WANG B, QI L, et al. Experimental study on hydrodynamic behaviors of high-speed gas jets in still water[J].ActaMechanicaSinica, 2006, 22(5): 443-448.

    [9] 施紅輝, 郭強(qiáng), 王超, 等. 水下超音速氣體射流脹鼓和回?fù)舻年P(guān)聯(lián)性研究[J]. 力學(xué)學(xué)報(bào), 2010, 42(6): 1206-1210.

    SHI Hong-hui, GUO Qiang, WANG Chao, et al. Experiments on the relationship between bulging and back-attack of submerged supersonic gas jets[J].ChineseJournalofTheoreticalandAppliedMechanics, 2010, 42(6): 1206-1210.

    [10] 齊麗婷, 余永剛, 彭志國, 等. 含能氣體射流在液體工質(zhì)中擴(kuò)展的兩維模型及數(shù)值模擬[J]. 含能材料, 2008, 16(2): 131-137.

    QI Li-ting, YU Yong-gang, PENG Zhi-guo, et al. A 2-D model of energetic gas jet expansion process in liquid and numerical simulation[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2008, 16(2): 131-137.

    [11] 曹嘉怡, 魯傳敬, 李杰, 等. 水下超聲速燃?xì)馍淞鲃?dòng)力學(xué)特性研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展A輯, 2009(05): 575-582.

    CAO Jia-yi, LU Chuan-jing, LI Jie, et al. Research on dynamic characteristics of underwater superasonic gas jets[J].JournalofHydrodynamicsser.A, 2009(05): 575-582.

    [12] 莽珊珊, 余永剛. 高壓燃?xì)馍淞髟谡b液體中擴(kuò)展過程的實(shí)驗(yàn)和數(shù)值模擬[J]. 爆炸與沖擊, 2011, 31(3): 300-305.

    MANG Shan-shan, YU Yong-gang. Experiment and numerical simulation for high pressure combustible gas jet expansion process in a bulk-loaded liquid[J].ExplosionandShockWaves, 2011, 31(3): 300-305.

    [13] Weiland C, Vlachos P. Round gas jets submerged in water[J].InternationalJournalofMultiphaseFlow, 2013, 48: 46-57.

    [14] YU Yong-gang, YAN Shao-heng, ZHAO Na, et al. Study on expansion process and interaction of high speed twin combustion-gas jets in liquid[J].JournalofAppliedMechanics, 2010, 77(5): 51404.

    [15] 薛曉春, 余永剛, 張琦. 雙束燃?xì)馍淞髋c整裝式液體裝藥相互作用的實(shí)驗(yàn)和數(shù)值模擬[J]. 兵工學(xué)報(bào), 2013, 34(6): 669-677.

    XUE Xiao-chun, YU Yong-gang, ZHANG Qi. Experiment and numerical simulation for interaction of twin gas jets and bulk-loaded liquid charge[J].ActaArmamentarii, 2013, 34(6): 669-677.

    [16] Graaf K, Brandner P, Penesis I. Bubble dynamics of a seismic air gun[J].ExperimentalThermalandFluidScience, 2014, 55: 228-238.

    [17] 陳鑫, 魯傳敬, 李杰, 等. VOF和Mixture多相流模型在空泡流模擬中的應(yīng)用[C]∥第九屆全國水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨第二十二屆全國水動(dòng)力學(xué)研討會(huì), 2009: 324-331.

    CHEN Xin, LU Chuan-jing, LI Jie, et al. Application in simulating cavitating flows by using VOF and mixture multi phase models[C]∥Proceedings of 9thNational Conference on Hydrodynamics and 22ndNational Congress on Hydrodynamics, 2009: 324-331.

    [18] 王福軍. 計(jì)算流體動(dòng)力學(xué)分析—CFD軟件原理與應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2004.

    猜你喜歡
    噴孔側(cè)壁斜面
    斜面之上探動(dòng)能
    高速公路隧道洞口淺埋段施工中的雙側(cè)壁導(dǎo)坑法
    巧用“相對(duì)”求解光滑斜面體問題
    巧用“相對(duì)”求解光滑斜面體問題
    柴油機(jī)噴油嘴變截面噴孔內(nèi)壁粗糙度影響研究
    柳鋼6號(hào)高爐應(yīng)對(duì)爐缸側(cè)壁溫度超高的護(hù)爐實(shí)踐
    四川冶金(2019年4期)2019-11-18 09:30:50
    越野車側(cè)壁加強(qiáng)筋板工藝改進(jìn)
    模具制造(2019年3期)2019-06-06 02:10:52
    基于CFD的噴嘴結(jié)構(gòu)參數(shù)對(duì)各孔內(nèi)部流動(dòng)特性影響研究
    一題多變 搞定斜面上的運(yùn)動(dòng)
    層合板上層建筑側(cè)壁抗空爆性能研究
    亚洲国产精品一区三区| 久久狼人影院| 精品亚洲成国产av| 欧美成人午夜免费资源| 欧美成人午夜免费资源| 精品一区二区三卡| av电影中文网址| 亚洲天堂av无毛| 熟女电影av网| 老司机影院毛片| 国产免费福利视频在线观看| 成人毛片60女人毛片免费| 国产成人一区二区在线| 国产精品麻豆人妻色哟哟久久| 亚洲不卡免费看| 高清毛片免费看| 国产欧美日韩一区二区三区在线 | 考比视频在线观看| 蜜桃在线观看..| 在线精品无人区一区二区三| 精品少妇黑人巨大在线播放| 99re6热这里在线精品视频| 亚洲av二区三区四区| 国产在线免费精品| 国产国语露脸激情在线看| 国产成人av激情在线播放 | 99视频精品全部免费 在线| 午夜av观看不卡| 免费av中文字幕在线| 2021少妇久久久久久久久久久| 高清毛片免费看| 免费大片黄手机在线观看| 免费高清在线观看视频在线观看| 精品一区在线观看国产| 亚洲精品第二区| 男人添女人高潮全过程视频| 永久免费av网站大全| 最近最新中文字幕免费大全7| 校园人妻丝袜中文字幕| 丝袜在线中文字幕| 色5月婷婷丁香| 亚洲激情五月婷婷啪啪| 中国国产av一级| 精品人妻熟女av久视频| 97超碰精品成人国产| 边亲边吃奶的免费视频| 少妇的逼好多水| 欧美 亚洲 国产 日韩一| 久久久久久久久久久丰满| 久久鲁丝午夜福利片| 亚洲精品久久久久久婷婷小说| 日本猛色少妇xxxxx猛交久久| 人妻人人澡人人爽人人| 欧美97在线视频| 精品熟女少妇av免费看| 欧美精品国产亚洲| 在现免费观看毛片| 午夜免费观看性视频| 人成视频在线观看免费观看| 大片电影免费在线观看免费| 熟女人妻精品中文字幕| 日韩一区二区视频免费看| 国产精品人妻久久久久久| 国产成人一区二区在线| 永久网站在线| 婷婷成人精品国产| 久久狼人影院| videosex国产| 国产白丝娇喘喷水9色精品| 国产免费现黄频在线看| 国产成人freesex在线| h视频一区二区三区| 熟妇人妻不卡中文字幕| 麻豆乱淫一区二区| 丝袜美足系列| 熟女人妻精品中文字幕| 欧美日韩国产mv在线观看视频| 秋霞在线观看毛片| 亚洲国产日韩一区二区| 久久久久久久久大av| 欧美日韩国产mv在线观看视频| 91成人精品电影| 熟女电影av网| 国产视频首页在线观看| 国产精品蜜桃在线观看| 欧美三级亚洲精品| 欧美国产精品一级二级三级| 999精品在线视频| 黄色配什么色好看| 国产精品秋霞免费鲁丝片| 中文字幕av电影在线播放| 99热这里只有精品一区| 99国产综合亚洲精品| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 亚洲欧美色中文字幕在线| 成年美女黄网站色视频大全免费 | 18+在线观看网站| 大话2 男鬼变身卡| 一边亲一边摸免费视频| 人妻系列 视频| 十八禁网站网址无遮挡| 日本vs欧美在线观看视频| 欧美日韩国产mv在线观看视频| 少妇精品久久久久久久| 亚洲av中文av极速乱| 国产高清国产精品国产三级| 天堂俺去俺来也www色官网| 欧美人与善性xxx| 观看av在线不卡| 一二三四中文在线观看免费高清| 少妇精品久久久久久久| 中文欧美无线码| 最近中文字幕高清免费大全6| 少妇 在线观看| 亚洲三级黄色毛片| 大又大粗又爽又黄少妇毛片口| 久久亚洲国产成人精品v| 成人亚洲精品一区在线观看| 亚洲av日韩在线播放| 视频区图区小说| 中文字幕最新亚洲高清| 欧美 日韩 精品 国产| 悠悠久久av| 久久狼人影院| 国产三级黄色录像| 啪啪无遮挡十八禁网站| 无人区码免费观看不卡 | 欧美黑人精品巨大| 日韩欧美三级三区| 免费观看人在逋| 久久性视频一级片| 怎么达到女性高潮| 天天影视国产精品| 亚洲色图 男人天堂 中文字幕| 一级毛片精品| 在线观看免费视频网站a站| 大香蕉久久成人网| 日韩免费高清中文字幕av| 日韩 欧美 亚洲 中文字幕| 国产黄色免费在线视频| 亚洲国产欧美网| 免费少妇av软件| 麻豆av在线久日| 黄色怎么调成土黄色| 日韩免费av在线播放| 黄色视频,在线免费观看| 日韩免费高清中文字幕av| 九色亚洲精品在线播放| 亚洲五月婷婷丁香| 免费不卡黄色视频| 窝窝影院91人妻| 成年人午夜在线观看视频| 国产欧美日韩一区二区三区在线| 在线永久观看黄色视频| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 大码成人一级视频| 一本大道久久a久久精品| av在线播放免费不卡| 香蕉丝袜av| 中文字幕人妻熟女乱码| 国产精品一区二区在线不卡| 精品少妇内射三级| 国产精品免费视频内射| 波多野结衣av一区二区av| 欧美老熟妇乱子伦牲交| 国内毛片毛片毛片毛片毛片| 最新的欧美精品一区二区| 成人国产av品久久久| 9热在线视频观看99| 日韩免费av在线播放| 人人妻人人爽人人添夜夜欢视频| 亚洲av日韩在线播放| 日韩一卡2卡3卡4卡2021年| 青青草视频在线视频观看| 成年版毛片免费区| 在线观看人妻少妇| 亚洲av美国av| 在线看a的网站| 国产精品秋霞免费鲁丝片| 99精品在免费线老司机午夜| 久久久精品区二区三区| 亚洲av电影在线进入| 91成年电影在线观看| 国产精品美女特级片免费视频播放器 | 成年人黄色毛片网站| 久久人妻av系列| av线在线观看网站| 一级,二级,三级黄色视频| 日日夜夜操网爽| 搡老岳熟女国产| 国产精品偷伦视频观看了| 午夜福利视频精品| 欧美老熟妇乱子伦牲交| 自拍欧美九色日韩亚洲蝌蚪91| 国产无遮挡羞羞视频在线观看| 国产av精品麻豆| www.熟女人妻精品国产| 在线观看人妻少妇| 日韩有码中文字幕| 亚洲五月色婷婷综合| 国产亚洲欧美在线一区二区| 免费av中文字幕在线| 自线自在国产av| 老鸭窝网址在线观看| 免费在线观看日本一区| av欧美777| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲色图av天堂| 亚洲精品国产精品久久久不卡| 自拍欧美九色日韩亚洲蝌蚪91| 欧美亚洲日本最大视频资源| 在线看a的网站| 精品国内亚洲2022精品成人 | 久久av网站| 操美女的视频在线观看| 日韩中文字幕欧美一区二区| 中文字幕人妻丝袜制服| 亚洲情色 制服丝袜| 女人高潮潮喷娇喘18禁视频| 国产熟女午夜一区二区三区| 国产精品98久久久久久宅男小说| 国产亚洲精品久久久久5区| 9色porny在线观看| 精品一品国产午夜福利视频| 免费不卡黄色视频| 18禁观看日本| 桃红色精品国产亚洲av| 欧美日韩av久久| 视频区欧美日本亚洲| 午夜成年电影在线免费观看| 国产在线观看jvid| 操美女的视频在线观看| 老熟妇仑乱视频hdxx| 国产精品久久久久久人妻精品电影 | 成年动漫av网址| 国产精品美女特级片免费视频播放器 | 国产有黄有色有爽视频| 欧美大码av| 免费黄频网站在线观看国产| 国产熟女午夜一区二区三区| 日本黄色视频三级网站网址 | 一级片'在线观看视频| 老熟妇仑乱视频hdxx| 久久免费观看电影| 香蕉久久夜色| 午夜福利在线免费观看网站| av国产精品久久久久影院| 国产xxxxx性猛交| 国产欧美日韩一区二区精品| 一级黄色大片毛片| 久久久国产精品麻豆| 久久影院123| 国产xxxxx性猛交| 免费在线观看黄色视频的| 变态另类成人亚洲欧美熟女 | 成人特级黄色片久久久久久久 | 9191精品国产免费久久| 国产欧美亚洲国产| 国产一区二区 视频在线| 免费少妇av软件| 99久久国产精品久久久| 另类精品久久| 99riav亚洲国产免费| 国产精品熟女久久久久浪| 国产三级黄色录像| 国产日韩一区二区三区精品不卡| 人妻 亚洲 视频| 欧美黄色淫秽网站| 久久精品国产a三级三级三级| 一区二区三区精品91| 五月天丁香电影| 大香蕉久久网| 精品少妇久久久久久888优播| 人妻 亚洲 视频| 久久久久视频综合| 无遮挡黄片免费观看| 免费在线观看黄色视频的| 精品少妇一区二区三区视频日本电影| 母亲3免费完整高清在线观看| 王馨瑶露胸无遮挡在线观看| 国产亚洲欧美在线一区二区| 窝窝影院91人妻| 久久中文看片网| 国产精品久久电影中文字幕 | 黄片小视频在线播放| 亚洲国产欧美在线一区| tube8黄色片| 后天国语完整版免费观看| 下体分泌物呈黄色| 纵有疾风起免费观看全集完整版| 老司机亚洲免费影院| 免费高清在线观看日韩| 精品少妇内射三级| 午夜福利免费观看在线| 一区二区三区国产精品乱码| 19禁男女啪啪无遮挡网站| 91大片在线观看| 啦啦啦免费观看视频1| 无人区码免费观看不卡 | 亚洲精品乱久久久久久| 免费黄频网站在线观看国产| 女人精品久久久久毛片| 亚洲少妇的诱惑av| 老司机午夜十八禁免费视频| 天天躁夜夜躁狠狠躁躁| 别揉我奶头~嗯~啊~动态视频| 久久国产精品影院| 男女边摸边吃奶| 国产在视频线精品| 蜜桃在线观看..| 日本vs欧美在线观看视频| 国产日韩欧美视频二区| 大陆偷拍与自拍| www.精华液| 国产精品久久久久久精品古装| 国产免费现黄频在线看| 欧美成人免费av一区二区三区 | 又黄又粗又硬又大视频| 国产免费av片在线观看野外av| 69精品国产乱码久久久| 中文字幕av电影在线播放| 丁香六月天网| 国产精品 国内视频| 亚洲精品久久午夜乱码| 国产深夜福利视频在线观看| 大陆偷拍与自拍| 人妻久久中文字幕网| 热re99久久精品国产66热6| 成人亚洲精品一区在线观看| 一区二区三区乱码不卡18| 丁香欧美五月| 久久九九热精品免费| 99九九在线精品视频| 精品一区二区三卡| 久久香蕉激情| 亚洲自偷自拍图片 自拍| 欧美变态另类bdsm刘玥| 欧美人与性动交α欧美精品济南到| 香蕉久久夜色| 高清视频免费观看一区二区| 国产片内射在线| 在线观看免费视频网站a站| 精品午夜福利视频在线观看一区 | 成人黄色视频免费在线看| 女同久久另类99精品国产91| 亚洲人成电影观看| 精品国产亚洲在线| 欧美精品一区二区免费开放| 欧美亚洲日本最大视频资源| 在线观看免费午夜福利视频| www.自偷自拍.com| 精品一区二区三区四区五区乱码| 国产在线一区二区三区精| 99热国产这里只有精品6| 国产高清视频在线播放一区| 国产熟女午夜一区二区三区| 好男人电影高清在线观看| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 1024香蕉在线观看| 欧美一级毛片孕妇| 天天躁夜夜躁狠狠躁躁| 一边摸一边抽搐一进一小说 | 成年人午夜在线观看视频| 99re6热这里在线精品视频| 老司机午夜十八禁免费视频| av免费在线观看网站| 岛国毛片在线播放| 亚洲成国产人片在线观看| 日韩一区二区三区影片| 国产成人精品在线电影| 在线播放国产精品三级| 久久人人97超碰香蕉20202| 欧美精品人与动牲交sv欧美| 精品国内亚洲2022精品成人 | 超碰97精品在线观看| 欧美激情高清一区二区三区| 黄片播放在线免费| 免费少妇av软件| 日韩有码中文字幕| 黑人巨大精品欧美一区二区mp4| 一夜夜www| 免费日韩欧美在线观看| 欧美在线一区亚洲| 色老头精品视频在线观看| 啪啪无遮挡十八禁网站| 在线观看免费视频日本深夜| 欧美成狂野欧美在线观看| 精品一区二区三区视频在线观看免费 | 18禁观看日本| 日本vs欧美在线观看视频| 亚洲成人免费av在线播放| 中文字幕色久视频| 国产高清激情床上av| 色尼玛亚洲综合影院| 1024香蕉在线观看| 日日夜夜操网爽| 色94色欧美一区二区| 亚洲国产欧美网| 丝袜美腿诱惑在线| 757午夜福利合集在线观看| 亚洲七黄色美女视频| 手机成人av网站| 极品教师在线免费播放| 亚洲色图综合在线观看| 亚洲,欧美精品.| 日韩欧美一区二区三区在线观看 | 一进一出好大好爽视频| 男女边摸边吃奶| 欧美性长视频在线观看| 久久精品人人爽人人爽视色| 伦理电影免费视频| 国产xxxxx性猛交| 日韩有码中文字幕| 在线观看免费视频日本深夜| 精品福利永久在线观看| 美国免费a级毛片| 日韩三级视频一区二区三区| 看免费av毛片| 亚洲精品美女久久久久99蜜臀| 伦理电影免费视频| 精品一区二区三卡| 色综合欧美亚洲国产小说| 午夜精品国产一区二区电影| 国产有黄有色有爽视频| 啦啦啦在线免费观看视频4| 欧美激情 高清一区二区三区| 欧美一级毛片孕妇| 无遮挡黄片免费观看| 三级毛片av免费| 三上悠亚av全集在线观看| 久久久久国内视频| 欧美人与性动交α欧美软件| 日本黄色日本黄色录像| 欧美国产精品一级二级三级| tube8黄色片| 精品国产乱码久久久久久小说| 午夜激情av网站| 免费女性裸体啪啪无遮挡网站| 91大片在线观看| 国产淫语在线视频| 成人国产av品久久久| 天天添夜夜摸| 国产99久久九九免费精品| 热99久久久久精品小说推荐| 首页视频小说图片口味搜索| 日韩免费av在线播放| 国产av又大| 伊人久久大香线蕉亚洲五| 国产区一区二久久| 亚洲精品久久成人aⅴ小说| 久久午夜综合久久蜜桃| 视频区欧美日本亚洲| 超碰成人久久| 嫁个100分男人电影在线观看| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| h视频一区二区三区| 国产在线精品亚洲第一网站| 91成年电影在线观看| 国产老妇伦熟女老妇高清| 日韩大码丰满熟妇| 久久人人97超碰香蕉20202| 少妇的丰满在线观看| 亚洲av欧美aⅴ国产| 天天躁日日躁夜夜躁夜夜| 亚洲中文av在线| 美女福利国产在线| 亚洲一码二码三码区别大吗| 热re99久久精品国产66热6| 国产麻豆69| 国产日韩欧美在线精品| 久久影院123| 亚洲av日韩精品久久久久久密| 国产成+人综合+亚洲专区| 一级片免费观看大全| 国产精品免费视频内射| 在线观看免费午夜福利视频| 日韩视频在线欧美| 淫妇啪啪啪对白视频| 久久国产精品大桥未久av| 国产精品久久久久久精品电影小说| 久久性视频一级片| 国产野战对白在线观看| 久久久久精品国产欧美久久久| 国产男靠女视频免费网站| 极品人妻少妇av视频| 欧美av亚洲av综合av国产av| 国产精品 欧美亚洲| 日韩视频在线欧美| 亚洲国产中文字幕在线视频| 黄色视频不卡| 不卡av一区二区三区| 十八禁人妻一区二区| 亚洲 国产 在线| 成年女人毛片免费观看观看9 | 性高湖久久久久久久久免费观看| 欧美国产精品va在线观看不卡| 丝袜人妻中文字幕| 嫩草影视91久久| 国产亚洲欧美精品永久| 亚洲精品美女久久久久99蜜臀| 久久精品亚洲精品国产色婷小说| 免费人妻精品一区二区三区视频| 新久久久久国产一级毛片| 在线播放国产精品三级| 一个人免费在线观看的高清视频| 亚洲午夜精品一区,二区,三区| 成人黄色视频免费在线看| tube8黄色片| 黑人欧美特级aaaaaa片| 日韩欧美免费精品| 精品国产超薄肉色丝袜足j| 国产极品粉嫩免费观看在线| 精品一品国产午夜福利视频| 国产免费视频播放在线视频| 国产精品.久久久| 一边摸一边抽搐一进一出视频| 波多野结衣一区麻豆| 在线观看人妻少妇| 久久久久久久久免费视频了| 国产精品av久久久久免费| 日韩一区二区三区影片| 亚洲免费av在线视频| 久久毛片免费看一区二区三区| 美女福利国产在线| 亚洲精品自拍成人| 女人高潮潮喷娇喘18禁视频| 久久人妻av系列| 日韩欧美免费精品| 国产一区二区三区视频了| 日韩中文字幕视频在线看片| 亚洲久久久国产精品| 免费在线观看视频国产中文字幕亚洲| 如日韩欧美国产精品一区二区三区| 男女午夜视频在线观看| 国产在线精品亚洲第一网站| 51午夜福利影视在线观看| 国产一区二区三区视频了| 精品欧美一区二区三区在线| 国产一区二区三区综合在线观看| 色老头精品视频在线观看| 精品熟女少妇八av免费久了| 成年动漫av网址| 午夜福利视频精品| 又大又爽又粗| 91成人精品电影| 制服人妻中文乱码| 男女之事视频高清在线观看| 飞空精品影院首页| 法律面前人人平等表现在哪些方面| av有码第一页| 水蜜桃什么品种好| 免费不卡黄色视频| 日本精品一区二区三区蜜桃| 热99国产精品久久久久久7| 男人操女人黄网站| 99re在线观看精品视频| 久久国产亚洲av麻豆专区| 久久久久久久久久久久大奶| 国产在线一区二区三区精| 国产精品一区二区在线不卡| 久久人妻福利社区极品人妻图片| 天堂俺去俺来也www色官网| 中文字幕色久视频| 19禁男女啪啪无遮挡网站| 国产精品99久久99久久久不卡| 一区二区av电影网| 精品乱码久久久久久99久播| 色婷婷久久久亚洲欧美| 午夜精品久久久久久毛片777| 一个人免费看片子| 免费在线观看日本一区| av超薄肉色丝袜交足视频| 自线自在国产av| 91大片在线观看| 精品免费久久久久久久清纯 | 女性生殖器流出的白浆| 在线观看免费视频日本深夜| 久久久久视频综合| 国产精品九九99| 一级a爱视频在线免费观看| 色婷婷av一区二区三区视频| 男女床上黄色一级片免费看| 色老头精品视频在线观看| 国产精品秋霞免费鲁丝片| 日韩欧美三级三区| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区三| 一区二区三区乱码不卡18| 一区在线观看完整版| 午夜福利乱码中文字幕| 咕卡用的链子| 麻豆国产av国片精品| 日韩欧美一区视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 免费观看a级毛片全部| 国产免费av片在线观看野外av| 成人永久免费在线观看视频 | 99re在线观看精品视频| 亚洲成人手机| 欧美成狂野欧美在线观看| 一区二区日韩欧美中文字幕| 一本综合久久免费| av福利片在线| 久久中文看片网| 曰老女人黄片| 最新在线观看一区二区三区| 19禁男女啪啪无遮挡网站| 麻豆乱淫一区二区| 大码成人一级视频| 国产精品一区二区在线观看99|