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

    柱形充液室內(nèi)多股燃?xì)馍淞髁鲌?chǎng)的氣體與液體兩相流場(chǎng)演化特性

    2016-11-09 10:08:44趙嘉俊余永剛
    兵工學(xué)報(bào) 2016年10期
    關(guān)鍵詞:柱形噴孔空腔

    趙嘉俊, 余永剛

    (南京理工大學(xué) 能源與動(dòng)力工程學(xué)院, 江蘇 南京 210094)

    ?

    柱形充液室內(nèi)多股燃?xì)馍淞髁鲌?chǎng)的氣體與液體兩相流場(chǎng)演化特性

    趙嘉俊, 余永剛

    (南京理工大學(xué) 能源與動(dòng)力工程學(xué)院, 江蘇 南京 210094)

    基于水下槍炮發(fā)射的工程背景,對(duì)火藥燃燒形成的多股高壓燃?xì)馍淞鲝哪M彈頭頭部噴入柱形充液室的過(guò)程進(jìn)行了實(shí)驗(yàn)觀測(cè)和數(shù)值模擬研究。實(shí)驗(yàn)觀測(cè)采用高速錄像記錄,數(shù)值模型采用流體體積函數(shù)模型來(lái)描述液體環(huán)境中多個(gè)燃?xì)馀莸臄U(kuò)展,數(shù)值模擬結(jié)果與實(shí)驗(yàn)觀測(cè)結(jié)果相符。數(shù)值仿真結(jié)果表明:噴孔軸線上的燃?xì)鈮毫?jīng)歷下降和上升的過(guò)程后達(dá)到穩(wěn)定值,空腔附近的水流場(chǎng)由于燃?xì)庾饔贸霈F(xiàn)壓力峰,隨后逐漸衰減;在燃?xì)饪涨粩U(kuò)展的過(guò)程中,離噴孔距離越遠(yuǎn)的截面上燃?xì)? 水流場(chǎng)的壓力分布越均勻,截面平均壓力隨著空腔擴(kuò)展逐漸下降;隨著時(shí)間的增長(zhǎng),截面上氣體組分分?jǐn)?shù)可達(dá)70%.

    兵器科學(xué)與技術(shù); 氣體與液體相互作用; 多股燃?xì)馍淞鳎?湍流摻混; 數(shù)值計(jì)算

    0 引言

    一種水下槍炮的發(fā)射方式是將其直接浸沒(méi)在水中發(fā)射,此時(shí)身管內(nèi)充滿水,彈丸出膛之前需要推動(dòng)其前方的液體運(yùn)動(dòng)。由于水的密度約為空氣的800倍,膛內(nèi)壓力隨著彈丸速度的增大急劇上升而易引發(fā)膛炸。借助導(dǎo)彈氣幕式發(fā)射的模式,通過(guò)在彈丸前端導(dǎo)入火藥燃燒形成的多股燃?xì)馍淞鱽?lái)排開(kāi)水柱,避免彈丸直接在水中運(yùn)動(dòng)。采用這種新方法可望降低槍炮在水下發(fā)射過(guò)程中燃燒室的壓力,提高發(fā)射安全性。為分析這種發(fā)射方式的作用機(jī)理,需要研究柱形充液室內(nèi)多股燃?xì)獾臄U(kuò)展及其形成Taylor空腔內(nèi)部的壓力分布特性。

    液體環(huán)境中多股燃?xì)獾臄U(kuò)展與火箭水下點(diǎn)火燃?xì)馀莸臄U(kuò)展相似,均為高溫高壓燃?xì)馍淞髋c液體相互作用問(wèn)題。趙世平等[1]實(shí)驗(yàn)研究了發(fā)射管內(nèi)充滿水時(shí)導(dǎo)彈的發(fā)射過(guò)程,表明燃?xì)鈬娙胨袝r(shí)由于水的不可壓縮性會(huì)形成一個(gè)較高的初始膛壓峰值。湯龍生等[2]對(duì)燃?xì)鈹U(kuò)展初期水中形成的壓力波進(jìn)行實(shí)驗(yàn)研究,表明該壓力波的壓力峰值衰減較快,對(duì)水介質(zhì)場(chǎng)的影響較小。為了深入了解燃?xì)獾牧鲃?dòng)狀態(tài),國(guó)內(nèi)外研究者采用數(shù)值方法來(lái)模擬液體環(huán)境中的燃?xì)鈹U(kuò)展。早期的研究采用零維等壓泡模型[3-4]來(lái)描述燃?xì)獾牧鲃?dòng)狀態(tài),隨后又出現(xiàn)一維[5]和二維模型[6]來(lái)描述燃?xì)狻T谇叭说幕A(chǔ)上,文獻(xiàn)[7-8]采用三維歐拉方程描述燃?xì)鉅顟B(tài),基于歐拉- 拉格朗日方法耦合燃?xì)? 水流場(chǎng)導(dǎo)彈水下發(fā)射的燃?xì)鈹U(kuò)展。Lindau等[9]基于歐拉- 歐拉方法提出了包含燃?xì)?、水和水汽的三相流?jì)算模型并對(duì)水下燃?xì)馔七M(jìn)過(guò)程進(jìn)行了模擬。同樣基于歐拉- 歐拉方法,曹嘉怡等[10-11]、王建儒等[12]和陳煥龍等[13]采用Mixture混合模型描述燃?xì)? 水的流場(chǎng),而向敏等[14]、唐嘉寧等[15-16]則采用流體體積函數(shù)(VOF)模型描述燃?xì)馀莸陌l(fā)展。相比于Mixture混合模型,VOF模型能獲得明顯的氣體與液體(簡(jiǎn)稱氣液)界面,適合描述燃?xì)鈹U(kuò)展初期氣液未完全混合的狀態(tài)。劉傳龍等[17]分析了發(fā)射過(guò)程中氣泡彈性對(duì)彈體的影響。

    對(duì)于以上研究,燃?xì)馍淞骶峭ㄈ腴_(kāi)放的液體環(huán)境,對(duì)于受限空間內(nèi)燃?xì)馍淞髋c液體相互作用問(wèn)題,主要是以整裝式液體炮為工程背景,莽珊珊等[18-19]和薛曉春等[20-21]分別開(kāi)展了單股和雙股燃?xì)馍淞髋c充液室液體相互作用特性研究。而對(duì)于水下槍炮發(fā)射這一工程背景,采用多股燃?xì)鈱?shí)時(shí)排水這種方法需要先研究錐狀分布的多股燃?xì)馍淞髟谥纬湟菏覂?nèi)擴(kuò)展以及氣液摻混問(wèn)題,該問(wèn)題未見(jiàn)報(bào)道。本文基于靜態(tài)模擬裝置的實(shí)驗(yàn),采用VOF模型對(duì)柱形充液室內(nèi)的多股燃?xì)鈹U(kuò)展進(jìn)行模擬。分析了多股燃?xì)馍淞鞯臄U(kuò)展特性,主要分析多個(gè)Taylor空腔內(nèi)部的壓力分布特征。

    1 數(shù)理模型

    1.1物理模型

    為研究柱形充液室內(nèi)多股燃?xì)馀潘O(shè)計(jì)了靜態(tài)模擬裝置來(lái)進(jìn)行實(shí)驗(yàn)?;鹚幵诿荛]高壓燃燒室內(nèi)燃燒形成燃?xì)?,燃?xì)鈮毫_(dá)到閾值時(shí)沖破密封膜片并通過(guò)模擬彈丸頭部的多個(gè)噴孔噴入充滿水的透明觀察室,由高速攝像機(jī)同步記錄燃?xì)馍淞魉纬蒚aylor空腔的演化過(guò)程。為研究多股燃?xì)獾臄U(kuò)展及氣液相互摻混的過(guò)程,對(duì)觀察室內(nèi)的區(qū)域進(jìn)行建模計(jì)算,計(jì)算模型和實(shí)驗(yàn)裝置的對(duì)比如圖1所示。模擬彈丸頭部開(kāi)有5個(gè)直徑分別為3 mm的噴孔,中心1個(gè),其余4個(gè)平均分布在錐形側(cè)面,側(cè)面噴孔與水平方向有45°夾角。觀察室內(nèi)徑55 mm,噴嘴頂部到觀察室頂部距離為127 mm,Photron Ultima APX i2高速攝像機(jī)拍攝頻率為4 000幀/s. 在模擬多股燃?xì)鈹U(kuò)展時(shí),計(jì)算模型的尺寸與模擬實(shí)驗(yàn)的參數(shù)保持一致,并在建模過(guò)程中做出如下假設(shè):

    1)多股燃?xì)馍淞鲊娙胫纬湟菏沂且粋€(gè)三維非穩(wěn)態(tài)過(guò)程,噴孔附近的雷諾數(shù)遠(yuǎn)高于臨界值,屬于完全湍流狀態(tài)。

    2)假設(shè)燃?xì)鉃榭蓧嚎s氣體,其狀態(tài)方程由理想狀態(tài)方程描述。

    3)實(shí)驗(yàn)持續(xù)時(shí)間僅幾毫秒,將不考慮氣液界面上的化學(xué)反應(yīng)及相變過(guò)程。

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

    基于上述假設(shè),可得控制方程如下:

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

    (1)

    式中:ρ代表混合密度,ρ=αgρg+αlρl,下標(biāo)g和l分別代表燃?xì)夂退?,αg、αl為相體積分?jǐn)?shù);u代表混合速度,u=(αgρgug+αlρlul)/ρ.

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

    (2)

    式中:p為靜壓;μ=αgμg+αlμl為混合動(dòng)力黏性系數(shù);μt為湍流黏性系數(shù);g為重力加速度(m/s2)。

    3)能量方程:

    (3)

    式中:ke為有效熱傳導(dǎo)系數(shù);

    (4)

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

    pg=ρgRTg,

    (5)

    式中:R為通用氣體常數(shù),值為8.31 J/(mol·K)。

    5)湍流模型。對(duì)于完全湍流狀態(tài),文獻(xiàn)[15-16]和薛曉春等[21]采用標(biāo)準(zhǔn)湍流模型來(lái)模擬燃?xì)獾呐蛎洝6鴮?duì)于柱形射流,Evgenevna等[22]認(rèn)為用可實(shí)現(xiàn)的k-ε湍流模型處理湍流會(huì)得到比標(biāo)準(zhǔn)湍流模型更好的結(jié)果。湍動(dòng)能k和湍流脈動(dòng)耗散率ε以及μt為

    Gk+Gb-ρε-YM,

    (6)

    (7)

    (8)

    式中:Gk為速度梯度產(chǎn)生的湍流動(dòng)能;Gb為浮力產(chǎn)生的湍流動(dòng)能;YM為可壓縮性的影響;常數(shù)依照文獻(xiàn)[23]設(shè)置:C1ε=1.44,C2=1.9,σk=1.0,σε=1.2.

    1.3初始條件與邊界條件

    圖2為計(jì)算模型邊界示意圖。初始狀態(tài)下圓柱形觀察室內(nèi)充滿液體:ρ0=ρl,T0=Tl,p0=pl. 依照實(shí)驗(yàn)研究條件,入口為高溫高壓燃?xì)猓鋮?shù)為燃?xì)鈪?shù):ρi=ρg,pi=pg,Ti=Tg,出口為大氣環(huán)境,其參數(shù)為大氣環(huán)境參數(shù):po=p∞,To=T∞. 壁面定義為絕熱且無(wú)滑移。

    圖2 計(jì)算模型邊界示意圖Fig.2 Schematic diagram of boundary in numerical model

    2 結(jié)果及分析

    2.1兩相流動(dòng)

    數(shù)值模型的計(jì)算區(qū)域?yàn)閷?shí)驗(yàn)觀測(cè)的觀察室區(qū)域,基于多孔噴嘴的對(duì)稱性,將計(jì)算區(qū)域簡(jiǎn)化為四分之一圓柱,網(wǎng)格平行于流動(dòng)方向并在壁面有加密處理,網(wǎng)格總數(shù)約為48萬(wàn),并通過(guò)無(wú)關(guān)性驗(yàn)證,采用FLUENT軟件來(lái)進(jìn)行數(shù)值模擬。針對(duì)模擬實(shí)驗(yàn)工況,入口處燃?xì)饪倝簽? MPa,溫度為2 000 K進(jìn)行數(shù)值模擬。圖3為0.5 ms時(shí)刻數(shù)值計(jì)算與實(shí)驗(yàn)觀測(cè)的比較圖,實(shí)驗(yàn)中高速攝像機(jī)從兩個(gè)相互垂直的方向來(lái)記錄燃?xì)獾臄U(kuò)展過(guò)程,圖中給出這兩個(gè)方向數(shù)值計(jì)算和實(shí)驗(yàn)觀測(cè)的結(jié)果??梢?jiàn)在t=0.5 ms時(shí)刻,燃?xì)馔ㄟ^(guò)模擬彈丸頂部的多個(gè)噴孔進(jìn)入液體形成多個(gè)單獨(dú)的Taylor空腔。圖4為1.0 ms時(shí)刻數(shù)值計(jì)算與實(shí)驗(yàn)觀測(cè)的對(duì)比圖,到t=1.0 ms時(shí)刻,側(cè)面噴孔形成的燃?xì)馍淞饕炎矒粲^察室壁面,側(cè)面射流的形狀及其氣體流動(dòng)方向均發(fā)生改變。在高速攝像機(jī)的記錄結(jié)果中,可以看到不同燃?xì)馍淞髦g依然存在間隙。圖5為3.0 ms時(shí)刻數(shù)值計(jì)算與實(shí)驗(yàn)觀測(cè)的對(duì)比圖,燃?xì)馍淞餍纬傻目涨痪蛏蠑U(kuò)展,側(cè)面射流和中心射流出現(xiàn)了聚并的過(guò)程,但頭部依然保持較大的空隙。圖6為射流頭部位移值比較圖,s為射流頭部在豎直方向上的位移,可見(jiàn)數(shù)值計(jì)算與實(shí)驗(yàn)測(cè)量的燃?xì)馍淞黝^部位移值相差較小。結(jié)合圖3~圖6的對(duì)比可以看出,數(shù)值計(jì)算結(jié)果反映了柱形充液室內(nèi)多股燃?xì)馍淞鲾U(kuò)展的主要特征。

    圖3 t=0.5 ms時(shí)數(shù)值計(jì)算與實(shí)驗(yàn)觀測(cè)對(duì)比圖Fig.3    Comparison between simulated and observed results for t=0.5 ms

    圖4 t=1.0 ms時(shí)數(shù)值計(jì)算與實(shí)驗(yàn)觀測(cè)對(duì)比圖Fig.4    Comparison between simulated and observed results for t=1.0 ms

    圖5 t=3.0 ms時(shí)數(shù)值計(jì)算與實(shí)驗(yàn)觀測(cè)對(duì)比圖Fig.5    Comparison between simulated and observed results for t=3.0 ms

    圖6 射流頭部位移值的比較Fig.6 Comparison of head displacements

    2.2流場(chǎng)壓力

    在水下武器的發(fā)射過(guò)程中,膛內(nèi)壓力是影響彈丸運(yùn)動(dòng)的主要因素,因此主要分析柱形充液室中多股燃?xì)馍淞鞯膲毫Ψ植继卣?。圖7為沿噴孔軸線的靜壓變化曲線圖,h為對(duì)應(yīng)噴孔軸線上的距離。由圖7(a)可見(jiàn)中心燃?xì)馍淞鬏S線上的壓力變化主要集中噴孔附近燃?xì)庥蓢娍讎姵鲋蟆H細(xì)鈴亩嗫讎娮靽姵鰰r(shí),靜壓約為2 MPa,由于噴孔處燃?xì)忪o壓遠(yuǎn)遠(yuǎn)高于環(huán)境靜壓,噴孔附近形成膨脹區(qū)域和壓縮區(qū)域,從而靜壓先下降、再上升。經(jīng)歷膨脹壓縮區(qū)域后,燃?xì)膺_(dá)到相對(duì)穩(wěn)定的壓力值,該壓力值隨著中心射流的發(fā)展而逐漸下降,由t=1.0 ms時(shí)刻的0.7 MPa下降到t=3.0 ms時(shí)刻的0.3 MPa. 在高速燃?xì)鈬娙氲退僖后w中時(shí),會(huì)產(chǎn)生“水錘效應(yīng)”,氣液交界面附近將出現(xiàn)一個(gè)壓力峰,即圖7(a)中t=0.5 ms靜壓曲線出現(xiàn)的壓力高峰,該峰值迅速下降直至消失。圖7(b)為側(cè)面燃?xì)馍淞鞯撵o壓變化曲線,由于側(cè)面燃?xì)饬鲃?dòng)方向在撞擊壁面后發(fā)生改變,因此側(cè)面射流在h=25 mm后的靜壓分布為沿壁面方向。與中心燃?xì)馍淞飨嗨?,?cè)面射流靜壓經(jīng)歷下降和上升的過(guò)程后達(dá)到了穩(wěn)定值,同時(shí)在t=0.5 ms時(shí)刻也出現(xiàn)壓力波動(dòng)。不同的是,側(cè)面射流在t=1.0 ms和t=1.5 ms時(shí)刻也出現(xiàn)了壓力峰。

    圖7 沿噴孔軸線的靜壓變化曲線Fig.7 Distribution of static pressure along nozzle axis

    圖8為不同時(shí)刻的壓力分布和相分布云圖,下半圖為相分布圖(淺色代表液體,深色代表燃?xì)?,其中x代表所選截面(直徑為55 mm)到噴嘴的軸向距離,r代表徑向距離。由圖8(a)可見(jiàn),在t=1.0 ms時(shí)刻,不同截面上的壓力分布有所不同。在x=10 mm截面上,中心燃?xì)馍淞鞲浇鼔毫ο鄬?duì)較高,壓力值較低的區(qū)域在側(cè)面射流內(nèi)部。在x=20 mm和x=30 mm截面上,中心燃?xì)馍淞鞲浇筛邏簠^(qū)域轉(zhuǎn)變?yōu)榈蛪簠^(qū)域,高壓區(qū)域?yàn)閭?cè)面射流附近的液體區(qū)域。結(jié)合圖4和圖8(a)可見(jiàn),側(cè)面射流撞擊壁面后頭部形狀發(fā)生變化,高壓區(qū)域處在燃?xì)馀c壁面之間,可見(jiàn)此時(shí)的高壓可能是由燃?xì)鈱?duì)壁面附近的液體作用而形成,這也解釋了圖7(b)中t=1.0 ms的壓力曲線上觀測(cè)到的壓力峰。在t=1.0 ms時(shí)刻,不同截面的平均壓力相差不大,約為0.7 MPa. 圖8(b)為t=3.0 ms時(shí)刻的壓力和相分布云圖,此時(shí)燃?xì)庹紦?jù)了柱形空間內(nèi)的大部分區(qū)域,在x=20 mm處還出現(xiàn)了相鄰側(cè)面射流聚并的現(xiàn)象。此時(shí)截面上氣相和液相區(qū)域之間的壓力差相對(duì)于t=1.0 ms時(shí)刻要小,截面的平均壓力由0.7 MPa下降到0.4 MPa. 可見(jiàn)隨著燃?xì)馍淞鞯臄U(kuò)展,氣液流場(chǎng)壓力趨向均勻分布的同時(shí)在逐漸降低。

    圖8 不同時(shí)刻的壓力分布和相分布云圖Fig.8 The static pressure and phase contours at different moment

    為反映氣體通道尺寸隨時(shí)間的變化特性,針對(duì)不同截面上的空氣所占面積百分?jǐn)?shù)進(jìn)行統(tǒng)計(jì)。圖9為x=10 mm,x=20 mm,x=30 mm和x=60 mm 4個(gè)截面上氣體所占百分?jǐn)?shù)隨時(shí)間的變化圖。由圖9可見(jiàn),不同截面上氣體的面積百分?jǐn)?shù)均隨著時(shí)間的增長(zhǎng)逐步增大。到t=3.0 ms時(shí)刻,x=10 mm,x=20 mm和x=30 mm 3個(gè)截面上的氣體所占面積百分?jǐn)?shù)可達(dá)70%.x=60 mm截面上的氣體百分?jǐn)?shù)快速增大。同時(shí)可見(jiàn),x=10 mm,x=20 mm和x=30 mm 3個(gè)截面上的氣體百分?jǐn)?shù)在t=2.0 ms時(shí)刻后增長(zhǎng)趨于平緩,這表明燃?xì)庠谶@些截面上的徑向擴(kuò)展已趨于穩(wěn)定。

    圖9 氣體組分分?jǐn)?shù)隨時(shí)間的變化圖Fig.9 Change of gas fraction over time

    3 結(jié)論

    本文利用數(shù)值計(jì)算方法模擬了柱形充液室內(nèi)多股燃?xì)馍淞鞯臄U(kuò)展,比較了模擬計(jì)算與實(shí)驗(yàn)觀測(cè)結(jié)果,分析了燃?xì)馀c水兩相流場(chǎng)的壓力分布規(guī)律,得到以下結(jié)論:

    1)在多股燃?xì)馍淞鲾U(kuò)展的模擬中,燃?xì)馍淞餍纬傻腡aylor空腔在初期相互獨(dú)立,相鄰側(cè)面射流隨著射流的擴(kuò)展發(fā)生匯聚。射流頭部位移的模擬結(jié)果與實(shí)驗(yàn)觀測(cè)結(jié)果吻合,數(shù)值結(jié)果能較好地描述兩相流動(dòng)的主要特征。

    2)燃?xì)鈬娙胨笮纬啥鄠€(gè)Taylor空腔,空腔內(nèi)氣流壓力經(jīng)歷膨脹區(qū)的下降和壓縮區(qū)的上升后達(dá)到穩(wěn)定值。在燃?xì)鈹U(kuò)展初期和側(cè)面燃?xì)馍淞髯矒舯诿孢^(guò)程中,Taylor空腔附近的水流場(chǎng)出現(xiàn)壓力峰,峰值逐漸減小最后消失。

    3)在Taylor空腔擴(kuò)展的過(guò)程中,燃?xì)馀c水流場(chǎng)不同截面上的壓力分布隨著截面離噴孔距離的增加而趨向均勻。不同截面上壓力的平均值在同一時(shí)刻相差不大,但截面平均壓力隨著Taylor空腔的擴(kuò)展而逐漸下降。

    4)不同截面上的氣體百分?jǐn)?shù)隨著時(shí)間的增長(zhǎng)逐步增大。到t=3.0 ms時(shí)刻,不同截面上的氣體所占面積百分?jǐn)?shù)可達(dá)70%.

    結(jié)論表明,多股燃?xì)馍淞骺稍谂懦耐瑫r(shí)形成低壓空氣區(qū)域,從而改善火炮的水下發(fā)射環(huán)境,為水下火炮的高速發(fā)射提供理論依據(jù)。

    References)

    [1]趙世平, 李江, 何國(guó)強(qiáng), 等. 固體燃?xì)獍l(fā)生器動(dòng)力模擬水下發(fā)射試驗(yàn)研究[J]. 固體火箭技術(shù), 2006, 29(1):5-8.

    ZHAO Shi-ping, LI Jiang, HE Guo-qiang, et al. Investigation on simulated underwater projectile test by solid propellant gas generator[J]. Journal of Solid Rocket Technology, 2006, 29(1):5-8.(in Chinese)

    [2]湯龍生, 劉宇, 吳智鋒, 等. 水下超聲速燃?xì)馍淞鳉馀莸纳L(zhǎng)及壓力波傳播特性實(shí)驗(yàn)研究[J]. 推進(jìn)技術(shù), 2011, 32(3):417-420.

    TANG Long-sheng, LIU Yu, WU Zhi-feng, et al. Experimental study on characteristics of bubble growth and pressure wave propagation by supersonic gas jets under water[J]. Journal of Propulsion Technology, 2011, 32(3):417-420.(in Chinese)

    [3]魯傳敬, 陳方, 樊泓, 等. 導(dǎo)彈水下點(diǎn)火的流體動(dòng)力研究[J]. 航空學(xué)報(bào), 1992,13(4):124-130.

    LU Chuan-jing, CHEN Fang, FAN Hong, et al. The fluid dynamic research on the under-water ignition of missile[J]. Acta Aeronautica et Astronautica Sinica, 1992, 13(4):124-130.(in Chinese)

    [4]黃建春, 葉取源, 朱世權(quán). 不同發(fā)射深度下導(dǎo)彈水下點(diǎn)火氣水流體動(dòng)力計(jì)算[J]. 應(yīng)用力學(xué)學(xué)報(bào), 1994, 11(3):19-24, 138.

    HUANG Jian-chun, YE Qu-yuan, ZHU Shi-quan. Gas-water dynamic calculation for the underwater ignition of a missile at different depths[J]. Chinese Journal of Applied Mechanics, 1994, 11(3):19-24, 138.(in Chinese)

    [5]王誠(chéng), 葉取源, 何友聲. 導(dǎo)彈水下發(fā)射燃?xì)馀萦?jì)算[J]. 應(yīng)用力學(xué)學(xué)報(bào), 1997, 14(3):3-9, 147.

    WANG Cheng, YE Qu-yuan, HE You-sheng. Calculation of an exhausted gas cavity behind an underwater launched missile[J]. Chinese Journal of Applied Mechanics, 1997, 14(3):3-9, 147.(in Chinese)

    [6]仲峰泉, 陸夕云, 莊禮賢. 火箭水下發(fā)射復(fù)雜流場(chǎng)的近似數(shù)值模擬[J]. 宇航學(xué)報(bào), 2000, 21(2): 1-7.

    ZHONG Feng-quan, LU Xi-yun, ZHUANG Li-xian. Numerical simulation of the complex flow field for rocket launch under water[J]. Journal of Astron Autics, 2000, 21(2):1-7.(in Chinese)

    [7]程用勝, 劉樺. 噴管不同擺角對(duì)水下發(fā)射導(dǎo)彈受力的影響[J]. 水動(dòng)力學(xué)研究與進(jìn)展A輯, 2007, 22(1):83-92.

    CHENG Yong-sheng, LIU Hua. The effects of nozzle deflection angles on a missile launched underwater[J]. Journal of Hydrodynamics, Ser. A, 2007, 22(1):83-92.(in Chinese)

    [8]Cheng Y S, Liu H. Mathematical modeling of fluid flows for underwater missile launch[J]. Journal of Hydrodynamics, Ser. B, 2006, 18(3):492-497.

    [9]Lindau J W, Venkateswaran S, Kunz R F, et al. Multiphase computations for underwater propulsive flows[C]∥16th AIAA Computational Fluid Dynamics Conference. Orlando, FL, US:AIAA, 2003.

    [10]曹嘉怡, 魯傳敬, 李杰, 等. 潛射導(dǎo)彈水下垂直自拋發(fā)射過(guò)程研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展A輯, 2006, 21(6):752-759.

    CAO Jia-yi, LU Chuan-jing, LI Jie, et al. Research on the vertical launching of direct ignition underwater missile[J]. Journal of Hydrodynamics, Ser. A, 2006, 21(6):752-759.(in Chinese)

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

    CAO Jia-yi, LU Chuan-jing, LI Jie, et al. Research on dynamics characteristics of underwater supersonic gas jets[J]. Journal of Hydrodynamics, Ser. A, 2009, 24(5):575-582.(in Chinese)

    [12]王建儒, 趙仕廠. 水下固體火箭發(fā)動(dòng)機(jī)尾流場(chǎng)計(jì)算[J]. 固體火箭技術(shù), 2007, 30(5): 388-391.

    WANG Jian-ru, ZHAO Shi-chang. Computations for solid rocket moter tail flow under water[J]. Journal of Solid Rocket Technology, 2007, 30(5):388-391.(in Chinese)

    [13]陳煥龍, 王檸, 劉華坪, 等. 不同發(fā)射深度下噴管燃?xì)馍淞魈匦匝芯縖J]. 水動(dòng)力學(xué)研究與進(jìn)展A輯, 2012, 27(6):659-666.

    CHEN Huan-long, WANG Ning, LIU Hua-ping, et al. Investigation of nozzle gas jet characteristics with different launch depth underwater[J]. Chinese Journal of Hydrodynamics, Ser. A, 2012, 27(6):659-666.(in Chinese)

    [14]向敏, 吳雄, 張為華, 等. 水下固體發(fā)動(dòng)機(jī)尾流場(chǎng)數(shù)值仿真[J]. 推進(jìn)技術(shù), 2009, 30(4): 479-483.

    XIANG Min, WU Xiong, ZHANG Wei-hua, et al. Numerical simulation for underwater solid motor tail flow[J]. Journal of Propulsion Technology, 2009, 30(4):479-483.(in Chinese)

    [15]唐嘉寧, 李世鵬, 王寧飛. 水下固體火箭發(fā)動(dòng)機(jī)的負(fù)推力現(xiàn)象研究[J]. 固體火箭技術(shù), 2012, 35(3):325-329, 343.

    TANG Jia-ning, LI Shi-peng, WANG Ning-fei. Study of the negative thrust of the underwater solid rocket engines[J]. Journal of Solid Rocket Technology, 2012, 35(3):325-329, 343.(in Chinese)

    [16]Tang J, Wang N, Shyy W. Flow structures of gaseous jets injected into water for underwater propulsion[J]. Acta Mechanica Sinica, 2011, 27(4):461-472.

    [17]劉傳龍, 張宇文, 王亞?wèn)|, 等. 考慮適配器彈性的潛射導(dǎo)彈出筒載荷特性研究[J]. 兵工學(xué)報(bào), 2015, 36(2):379-384.

    LIU Chuan-long, ZHANG Yu-wen, WANG Ya-dong, et al. Investigation into load characteristics of submarine-launched missile being ejected from launch tube considering the adapter elasticity[J]. Acta Armamentarii, 2015, 36(2):379-384.(in Chinese)

    [18]莽珊珊, 余永剛. 邊界形狀對(duì)高壓受限射流擴(kuò)展穩(wěn)定性影響[J]. 推進(jìn)技術(shù), 2011, 32(3): 334-338.

    MANG Shan-shan, YU Yong-gang. Effects of chamber wall shape on high-pressure gas jet expansion stability in bulk-loaded liquid[J]. Journal of Propulsion Technology, 2011, 32(3):334-338.(in Chinese)

    [19]莽珊珊, 余永剛. 高壓燃?xì)馍淞髟谡b液體中擴(kuò)展過(guò)程的實(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]. Explosion and Shock Waves, 2011, 31(3):300-305.(in Chinese)

    [20]薛曉春, 余永剛, 張琦, 等. 漸擴(kuò)邊界形狀對(duì)雙股燃?xì)馍淞鲾U(kuò)展特性影響的實(shí)驗(yàn)研究[J]. 彈道學(xué)報(bào), 2013, 25(2): 44-47.

    XUE Xiao-chun, YU Yong-gang, ZHANG Qi, et al. Experimental study on effects of stepped-wall boundary on expansion characteristic of twin combustion-gas Jets[J]. Journal of Ballistics, 2013, 25(2):44-47.(in Chinese)

    [21]薛曉春, 余永剛, 張琦. 雙束燃?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]. Acta Armamentarii, 2013, 34(6):669-677.(in Chinese)

    [22]Evgenevna I E, Evgenevna I T, Viktorovich B P. Analysis of the application of turbulence models in the calculation of supersonic gas jet[J]. American Journal of Applied Sciences, 2014, 11(11):1914-1920.

    [23]Shih T H, Liou W W, Shabbir A, et al. A newk-εeddy viscosity model for high Reynolds number turbulent flows[J]. Computers & Fluids, 1995, 24(3):227-238.

    The Evolution of Gas-liquid Flow of Multiple Combustion Gas Jets in a Cylindrical Water-filled Chamber

    ZHAO Jia-jun, YU Yong-gang

    (School of Energy and Power Engineering, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China)

    The injection process of multiple high-pressure combustion gas jets is observed experimentally and simulated numerically during gun firing underwater. The gas is generated by gunpowder burning and injected into a cylindrical water-filled chamber through simulated bullet. In experiment, a high-speed camera system is used to record the injecting process of air. The volume-of-fluid method is used to describe the expansion of the gas cavities in water. The numerical results show a good agreement with experimental observations. Numerical simulations show that the gas pressure along the axis of orifice reaches a steady value after declining and rising. A pressure peak appears in the water near the gas cavity under the effect of combustion gas, and then is gradually attenuated. During the expansion of the gas cavities, the pressure distribution of gas-water flow field on a cross section away from nozzle is relatively uniform, and the average pressure on the section decreases with the expansion in the cavity. The area fraction of gas on the section could be up to 70% over time.

    ordnance science and technology; gas-water interaction; multiple combustion gas jet; turbulent mixing; numerical simulation

    2016-01-06

    國(guó)家自然科學(xué)基金項(xiàng)目(11372139)

    趙嘉俊(1990—),男,博士研究生。E-mail: njustzhaojiajun@126.com;

    余永剛(1963—),男,教授,博士生導(dǎo)師。E-mail:yygnjust801@163.com

    O358

    A

    1000-1093(2016)10-1852-08

    10.3969/j.issn.1000-1093.2016.10.012

    猜你喜歡
    柱形噴孔空腔
    基于邊光滑有限元法的二維復(fù)合彈性空腔聲振特性分析
    柴油機(jī)噴油嘴變截面噴孔內(nèi)壁粗糙度影響研究
    基于CFD的噴嘴結(jié)構(gòu)參數(shù)對(duì)各孔內(nèi)部流動(dòng)特性影響研究
    非柱形容器的壓力和壓強(qiáng)
    從“柱形凸透鏡成像”實(shí)驗(yàn)感受體驗(yàn)教學(xué)的魅力
    空腔參數(shù)對(duì)重力壩穩(wěn)定的影響分析
    前置污水去油池
    前置污水去油池
    基于Fluent的空氣射流切削式反循環(huán)鉆頭參數(shù)優(yōu)化
    鉆探工程(2015年11期)2015-01-01 02:53:50
    溪洛渡水電站GIL柱形絕緣子局部放電原因分析
    欧美日韩av久久| 成年人免费黄色播放视频| 51午夜福利影视在线观看| bbb黄色大片| 国产精品一区二区免费欧美 | 亚洲激情五月婷婷啪啪| 大码成人一级视频| 国产免费一区二区三区四区乱码| 每晚都被弄得嗷嗷叫到高潮| 欧美黄色片欧美黄色片| av网站免费在线观看视频| 十分钟在线观看高清视频www| 日日夜夜操网爽| 国产精品99久久99久久久不卡| 在线观看一区二区三区激情| 亚洲av电影在线进入| 别揉我奶头~嗯~啊~动态视频 | 国产男女超爽视频在线观看| 18禁黄网站禁片午夜丰满| 国产99久久九九免费精品| 亚洲三区欧美一区| 1024香蕉在线观看| a在线观看视频网站| www.精华液| 在线十欧美十亚洲十日本专区| 新久久久久国产一级毛片| 色94色欧美一区二区| 最黄视频免费看| 久久人人爽av亚洲精品天堂| 国产精品1区2区在线观看. | 一级黄色大片毛片| 美女午夜性视频免费| 国产男女超爽视频在线观看| 亚洲国产av影院在线观看| 日日摸夜夜添夜夜添小说| 国产成人欧美| 老汉色av国产亚洲站长工具| 桃红色精品国产亚洲av| 国产精品一区二区在线不卡| 国产亚洲av高清不卡| 国产免费av片在线观看野外av| 天天操日日干夜夜撸| 日韩中文字幕视频在线看片| 一区二区三区乱码不卡18| 嫩草影视91久久| 深夜精品福利| 亚洲七黄色美女视频| 亚洲黑人精品在线| 搡老岳熟女国产| 在线天堂中文资源库| 国产精品.久久久| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美一区二区三区久久| 侵犯人妻中文字幕一二三四区| 乱人伦中国视频| 欧美变态另类bdsm刘玥| 国产三级黄色录像| 久久精品国产亚洲av高清一级| 国产一区二区三区av在线| 欧美黄色淫秽网站| 久久热在线av| 欧美人与性动交α欧美软件| 12—13女人毛片做爰片一| 啦啦啦免费观看视频1| 日本黄色日本黄色录像| 日本a在线网址| 国产精品成人在线| 精品欧美一区二区三区在线| 亚洲情色 制服丝袜| 免费在线观看视频国产中文字幕亚洲 | 老司机影院成人| 国产真人三级小视频在线观看| 亚洲欧美精品自产自拍| 成人三级做爰电影| a 毛片基地| 热99久久久久精品小说推荐| 大香蕉久久成人网| 国产99久久九九免费精品| 狠狠精品人妻久久久久久综合| 国产在线免费精品| 国产麻豆69| 国产免费一区二区三区四区乱码| 国产日韩欧美亚洲二区| 亚洲七黄色美女视频| 亚洲国产中文字幕在线视频| 国产精品1区2区在线观看. | 国产一级毛片在线| 亚洲av日韩精品久久久久久密| bbb黄色大片| 欧美日韩国产mv在线观看视频| 精品人妻一区二区三区麻豆| 满18在线观看网站| 欧美少妇被猛烈插入视频| 国产日韩欧美视频二区| 国产不卡av网站在线观看| 亚洲七黄色美女视频| 人妻 亚洲 视频| 三上悠亚av全集在线观看| 欧美在线黄色| 天天躁狠狠躁夜夜躁狠狠躁| 天天添夜夜摸| 一区二区三区精品91| 久久ye,这里只有精品| 精品少妇一区二区三区视频日本电影| 国产主播在线观看一区二区| 亚洲精品粉嫩美女一区| 国产精品久久久人人做人人爽| av天堂在线播放| 人妻一区二区av| 免费高清在线观看日韩| 桃红色精品国产亚洲av| 精品欧美一区二区三区在线| 国产高清国产精品国产三级| 91九色精品人成在线观看| 97人妻天天添夜夜摸| a 毛片基地| 亚洲国产成人一精品久久久| 精品人妻1区二区| 热99久久久久精品小说推荐| 肉色欧美久久久久久久蜜桃| 热99久久久久精品小说推荐| 宅男免费午夜| 黄色视频,在线免费观看| 国产黄色免费在线视频| 一本—道久久a久久精品蜜桃钙片| av一本久久久久| 啦啦啦在线免费观看视频4| 大片免费播放器 马上看| 亚洲五月色婷婷综合| videos熟女内射| 久久精品人人爽人人爽视色| 日韩人妻精品一区2区三区| 国产成人免费观看mmmm| 精品第一国产精品| 亚洲成人手机| 18禁观看日本| 国产深夜福利视频在线观看| 91麻豆精品激情在线观看国产 | 久久性视频一级片| 日韩中文字幕视频在线看片| 真人做人爱边吃奶动态| 国产亚洲一区二区精品| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 99精国产麻豆久久婷婷| 国产福利在线免费观看视频| 一进一出抽搐动态| 日日夜夜操网爽| 12—13女人毛片做爰片一| 性高湖久久久久久久久免费观看| 欧美精品一区二区免费开放| 国产精品一区二区免费欧美 | 欧美日韩福利视频一区二区| 真人做人爱边吃奶动态| 天堂俺去俺来也www色官网| 亚洲精品成人av观看孕妇| 巨乳人妻的诱惑在线观看| 日本a在线网址| svipshipincom国产片| 精品一区二区三区四区五区乱码| 亚洲国产中文字幕在线视频| 欧美国产精品va在线观看不卡| 老司机靠b影院| 在线 av 中文字幕| 国产亚洲av片在线观看秒播厂| 99精国产麻豆久久婷婷| 日韩大片免费观看网站| 91精品国产国语对白视频| 水蜜桃什么品种好| 在线观看一区二区三区激情| 一区二区三区激情视频| 免费在线观看影片大全网站| 亚洲精品久久成人aⅴ小说| 女人久久www免费人成看片| 亚洲精品日韩在线中文字幕| 制服人妻中文乱码| 亚洲欧美精品综合一区二区三区| 欧美精品高潮呻吟av久久| 亚洲性夜色夜夜综合| 一进一出抽搐动态| 日日爽夜夜爽网站| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| a在线观看视频网站| 精品一区在线观看国产| av免费在线观看网站| 亚洲国产看品久久| 国产成人精品久久二区二区91| 亚洲第一欧美日韩一区二区三区 | 久久影院123| 秋霞在线观看毛片| 99re6热这里在线精品视频| 18在线观看网站| xxxhd国产人妻xxx| 搡老岳熟女国产| 国产av又大| 日韩免费高清中文字幕av| 欧美午夜高清在线| 女人精品久久久久毛片| 国产精品国产三级国产专区5o| 免费高清在线观看日韩| 久久国产精品人妻蜜桃| 日本91视频免费播放| 久久精品亚洲av国产电影网| 国产精品久久久久成人av| 高清视频免费观看一区二区| 中文精品一卡2卡3卡4更新| 国产无遮挡羞羞视频在线观看| 国产亚洲精品久久久久5区| 亚洲情色 制服丝袜| 亚洲国产毛片av蜜桃av| 制服诱惑二区| a级毛片黄视频| 日韩,欧美,国产一区二区三区| 亚洲第一欧美日韩一区二区三区 | 一个人免费看片子| 欧美人与性动交α欧美精品济南到| 亚洲欧美清纯卡通| 一二三四社区在线视频社区8| 国产成人欧美在线观看 | 水蜜桃什么品种好| 国产精品自产拍在线观看55亚洲 | 2018国产大陆天天弄谢| 三上悠亚av全集在线观看| 老鸭窝网址在线观看| 久久国产亚洲av麻豆专区| 国产av一区二区精品久久| 国产精品秋霞免费鲁丝片| 又紧又爽又黄一区二区| 老司机影院毛片| 91大片在线观看| 91成人精品电影| 欧美日韩视频精品一区| 亚洲一区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 久9热在线精品视频| 老司机亚洲免费影院| 欧美日本中文国产一区发布| 法律面前人人平等表现在哪些方面 | 亚洲av欧美aⅴ国产| 亚洲成人手机| 在线精品无人区一区二区三| 法律面前人人平等表现在哪些方面 | 免费人妻精品一区二区三区视频| 免费日韩欧美在线观看| 午夜福利免费观看在线| 久久久久久久久免费视频了| avwww免费| 极品人妻少妇av视频| 国产精品国产三级国产专区5o| 国产伦人伦偷精品视频| 久久国产亚洲av麻豆专区| kizo精华| 2018国产大陆天天弄谢| 极品少妇高潮喷水抽搐| 美女主播在线视频| 一二三四在线观看免费中文在| 99精国产麻豆久久婷婷| 一级毛片精品| 欧美激情 高清一区二区三区| 首页视频小说图片口味搜索| 好男人电影高清在线观看| 欧美亚洲 丝袜 人妻 在线| 欧美人与性动交α欧美软件| 国产精品九九99| 成人影院久久| 一区福利在线观看| 欧美黄色淫秽网站| 91字幕亚洲| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| 五月开心婷婷网| 国产日韩一区二区三区精品不卡| 成年动漫av网址| e午夜精品久久久久久久| 亚洲国产欧美日韩在线播放| 免费高清在线观看视频在线观看| 亚洲国产精品成人久久小说| 最近中文字幕2019免费版| 脱女人内裤的视频| 91麻豆精品激情在线观看国产 | 一二三四社区在线视频社区8| 国产成人啪精品午夜网站| 精品人妻1区二区| 国产熟女午夜一区二区三区| 韩国高清视频一区二区三区| 久久精品成人免费网站| 美女国产高潮福利片在线看| 欧美激情高清一区二区三区| 国产免费视频播放在线视频| 高清视频免费观看一区二区| av天堂久久9| 精品视频人人做人人爽| 一本—道久久a久久精品蜜桃钙片| 18禁观看日本| 麻豆国产av国片精品| 亚洲欧美精品综合一区二区三区| 国产成人精品久久二区二区91| 亚洲欧美一区二区三区久久| 久久久国产精品麻豆| 久久久久久久久免费视频了| 1024香蕉在线观看| 久久久精品94久久精品| 啦啦啦免费观看视频1| 操美女的视频在线观看| 国产深夜福利视频在线观看| 老司机深夜福利视频在线观看 | 欧美激情久久久久久爽电影 | 国产精品亚洲av一区麻豆| 久久久久久久大尺度免费视频| 免费人妻精品一区二区三区视频| 亚洲欧美日韩另类电影网站| 精品福利永久在线观看| 国产精品熟女久久久久浪| 国产在线视频一区二区| 精品欧美一区二区三区在线| 男人爽女人下面视频在线观看| 欧美成狂野欧美在线观看| 无限看片的www在线观看| 成年人免费黄色播放视频| 久久精品人人爽人人爽视色| 午夜激情久久久久久久| 国产亚洲一区二区精品| 51午夜福利影视在线观看| 久久久久久久久免费视频了| 亚洲av成人不卡在线观看播放网 | 国产区一区二久久| 午夜福利影视在线免费观看| 91麻豆精品激情在线观看国产 | 亚洲av成人一区二区三| 人妻 亚洲 视频| 人人妻人人澡人人看| 国产亚洲精品一区二区www | 免费在线观看视频国产中文字幕亚洲 | 国产精品 国内视频| 国精品久久久久久国模美| 欧美97在线视频| 亚洲色图 男人天堂 中文字幕| 我要看黄色一级片免费的| 亚洲成人国产一区在线观看| 老司机午夜十八禁免费视频| 狂野欧美激情性bbbbbb| 视频区欧美日本亚洲| 欧美另类亚洲清纯唯美| 美女视频免费永久观看网站| 一个人免费在线观看的高清视频 | 亚洲中文字幕日韩| 国产一区二区激情短视频 | 美女中出高潮动态图| 青草久久国产| 午夜91福利影院| 亚洲成人国产一区在线观看| 亚洲av男天堂| 亚洲专区国产一区二区| 成人亚洲精品一区在线观看| 午夜久久久在线观看| 亚洲精品国产区一区二| 中文字幕人妻丝袜制服| 亚洲av电影在线观看一区二区三区| 少妇粗大呻吟视频| 大型av网站在线播放| 国精品久久久久久国模美| 精品视频人人做人人爽| av福利片在线| 男人爽女人下面视频在线观看| 国产欧美亚洲国产| 国产精品免费大片| 黄色 视频免费看| 久久精品国产亚洲av高清一级| 精品一区二区三区av网在线观看 | 亚洲专区国产一区二区| 国产成人精品在线电影| 男女国产视频网站| 动漫黄色视频在线观看| 成年美女黄网站色视频大全免费| 欧美大码av| 99国产极品粉嫩在线观看| 成年美女黄网站色视频大全免费| 日本av手机在线免费观看| tube8黄色片| 一区二区三区精品91| 亚洲专区字幕在线| 一本一本久久a久久精品综合妖精| 日本一区二区免费在线视频| 老司机在亚洲福利影院| www.av在线官网国产| 国产男女内射视频| 午夜视频精品福利| 久久久久久久久久久久大奶| 大片免费播放器 马上看| 久久精品国产亚洲av香蕉五月 | 国产欧美日韩一区二区三 | 超碰成人久久| 欧美日本中文国产一区发布| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久成人aⅴ小说| bbb黄色大片| 欧美成人午夜精品| 性色av乱码一区二区三区2| av天堂久久9| 巨乳人妻的诱惑在线观看| 一进一出抽搐动态| 黄网站色视频无遮挡免费观看| 亚洲欧美精品自产自拍| 精品卡一卡二卡四卡免费| 在线永久观看黄色视频| 国产国语露脸激情在线看| 欧美精品亚洲一区二区| 精品人妻熟女毛片av久久网站| 一本色道久久久久久精品综合| 性色av一级| 51午夜福利影视在线观看| 国精品久久久久久国模美| 人妻 亚洲 视频| 一区二区av电影网| 侵犯人妻中文字幕一二三四区| 99国产综合亚洲精品| 我的亚洲天堂| 免费少妇av软件| 欧美日韩视频精品一区| 久久精品亚洲熟妇少妇任你| 女性生殖器流出的白浆| 国产人伦9x9x在线观看| 欧美在线黄色| 久久精品久久久久久噜噜老黄| 久热爱精品视频在线9| 日韩视频一区二区在线观看| 50天的宝宝边吃奶边哭怎么回事| 久久国产精品男人的天堂亚洲| 免费在线观看黄色视频的| 秋霞在线观看毛片| 伊人亚洲综合成人网| 成年美女黄网站色视频大全免费| 脱女人内裤的视频| av视频免费观看在线观看| 国产成人免费观看mmmm| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品日韩在线中文字幕| 久久精品亚洲熟妇少妇任你| 性高湖久久久久久久久免费观看| 欧美精品av麻豆av| 国产在线观看jvid| 国产欧美日韩综合在线一区二区| 国产亚洲av高清不卡| 在线看a的网站| 国精品久久久久久国模美| 久久毛片免费看一区二区三区| 日韩一卡2卡3卡4卡2021年| 麻豆乱淫一区二区| 99国产精品一区二区蜜桃av | 美女高潮喷水抽搐中文字幕| 日本猛色少妇xxxxx猛交久久| 成年av动漫网址| 香蕉国产在线看| av免费在线观看网站| 免费观看av网站的网址| 99久久综合免费| av视频免费观看在线观看| 亚洲精品乱久久久久久| 老熟妇仑乱视频hdxx| a在线观看视频网站| 久久国产精品影院| 69av精品久久久久久 | 黄色视频在线播放观看不卡| 最近中文字幕2019免费版| 亚洲精品av麻豆狂野| 成年动漫av网址| 2018国产大陆天天弄谢| 一本综合久久免费| 老司机影院毛片| 热re99久久国产66热| 亚洲精品美女久久av网站| 久久午夜综合久久蜜桃| 欧美97在线视频| 狂野欧美激情性xxxx| 精品亚洲成国产av| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 欧美日韩视频精品一区| 热re99久久精品国产66热6| 亚洲精品久久久久久婷婷小说| 午夜激情av网站| 国产福利在线免费观看视频| 亚洲专区字幕在线| 亚洲性夜色夜夜综合| 99国产精品99久久久久| 老汉色av国产亚洲站长工具| 久久影院123| 他把我摸到了高潮在线观看 | 久久女婷五月综合色啪小说| 大香蕉久久成人网| 成人18禁高潮啪啪吃奶动态图| 青春草视频在线免费观看| 各种免费的搞黄视频| 国产日韩欧美亚洲二区| 精品国产一区二区久久| 国产成人欧美| 如日韩欧美国产精品一区二区三区| 一区二区三区激情视频| 亚洲av国产av综合av卡| 国产精品免费视频内射| 亚洲精品日韩在线中文字幕| 久热爱精品视频在线9| 欧美在线一区亚洲| 色老头精品视频在线观看| 高清黄色对白视频在线免费看| 一本大道久久a久久精品| 国产精品av久久久久免费| 9191精品国产免费久久| 美女高潮喷水抽搐中文字幕| 9色porny在线观看| 亚洲男人天堂网一区| 伦理电影免费视频| 免费少妇av软件| 亚洲综合色网址| 午夜福利乱码中文字幕| www.精华液| 人妻 亚洲 视频| 亚洲欧洲精品一区二区精品久久久| 在线永久观看黄色视频| 亚洲欧洲精品一区二区精品久久久| 日韩欧美一区视频在线观看| 夫妻午夜视频| 久久狼人影院| 亚洲国产中文字幕在线视频| 午夜福利乱码中文字幕| 国产高清videossex| 久久毛片免费看一区二区三区| av一本久久久久| 亚洲,欧美精品.| 国产成人啪精品午夜网站| av视频免费观看在线观看| 久久久久精品国产欧美久久久 | 久久精品亚洲av国产电影网| 国产精品久久久久久精品古装| 国产av又大| 老熟妇仑乱视频hdxx| 国产精品秋霞免费鲁丝片| 亚洲成人免费av在线播放| 国产成人欧美在线观看 | 精品高清国产在线一区| 一个人免费在线观看的高清视频 | 国产真人三级小视频在线观看| 午夜激情久久久久久久| 一区二区日韩欧美中文字幕| av天堂久久9| 午夜福利,免费看| 久久青草综合色| 91成年电影在线观看| 国产成+人综合+亚洲专区| 亚洲精品中文字幕一二三四区 | 亚洲精品国产区一区二| 最黄视频免费看| 精品亚洲成a人片在线观看| 各种免费的搞黄视频| 精品人妻在线不人妻| 国产亚洲午夜精品一区二区久久| 80岁老熟妇乱子伦牲交| www.精华液| 午夜福利,免费看| 国精品久久久久久国模美| 国产男女超爽视频在线观看| 亚洲自偷自拍图片 自拍| 欧美精品一区二区大全| 亚洲欧美激情在线| 青草久久国产| 男人舔女人的私密视频| 十八禁人妻一区二区| 12—13女人毛片做爰片一| 女警被强在线播放| 97在线人人人人妻| 成人国语在线视频| 欧美av亚洲av综合av国产av| 熟女少妇亚洲综合色aaa.| 国产又色又爽无遮挡免| av天堂久久9| 老汉色av国产亚洲站长工具| 亚洲国产欧美网| 国产亚洲精品久久久久5区| 曰老女人黄片| 成在线人永久免费视频| 精品久久蜜臀av无| av在线老鸭窝| 国产亚洲午夜精品一区二区久久| 亚洲成人国产一区在线观看| av在线播放精品| 成人免费观看视频高清| 9191精品国产免费久久| 亚洲av日韩在线播放| 国产精品99久久99久久久不卡| 免费人妻精品一区二区三区视频| 亚洲精品国产av蜜桃| 十分钟在线观看高清视频www| 精品一区二区三卡| 夫妻午夜视频| 一二三四社区在线视频社区8| 王馨瑶露胸无遮挡在线观看| 国产伦人伦偷精品视频| 狂野欧美激情性bbbbbb| 久久香蕉激情| 99久久人妻综合| 中文字幕色久视频| 国产日韩欧美在线精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲男人天堂网一区| 色播在线永久视频| 永久免费av网站大全| 欧美成人午夜精品| 免费在线观看黄色视频的| 国产精品麻豆人妻色哟哟久久| 国产精品免费大片| 在线精品无人区一区二区三| 欧美乱码精品一区二区三区| a级毛片在线看网站| 欧美日韩亚洲高清精品|