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

    高壓甲烷射流沖擊預(yù)混火焰過(guò)程中渦量的變化

    2023-09-26 07:34:08梁曉杰
    關(guān)鍵詞:渦量量值前鋒

    雷 艷,王 瑩,仇 滔,梁曉杰

    (北京工業(yè)大學(xué) 環(huán)境與生命學(xué)部,北京 100124)

    隨著日益嚴(yán)格的低碳化要求,清潔代用燃料的應(yīng)用成為發(fā)動(dòng)機(jī)技術(shù)研究熱點(diǎn),天然氣是十分有前景的內(nèi)燃機(jī)替代燃料[1].天然氣高壓直噴技術(shù)可以有效提高發(fā)動(dòng)機(jī)熱效率并降低排放,天然氣高壓噴射進(jìn)入氣缸,并被預(yù)噴燃料形成的火焰引燃.天然氣射流與引燃火焰間的相互作用決定了燃燒過(guò)程.因此,研究氣體射流沖擊火焰有非常重要的實(shí)際意義[2].

    在火災(zāi)消防領(lǐng)域橫掠風(fēng)對(duì)擴(kuò)散火焰是有影響的,針對(duì)環(huán)境變量風(fēng)速,蔣依[3]以甲醇汽油為試驗(yàn)燃料,了解有風(fēng)條件下甲醇汽油燃燒發(fā)展的全過(guò)程,得出風(fēng)速對(duì)火焰傾角、火焰高度和火焰尺寸的影響機(jī)制及規(guī)律.莊磊[4]以航空煤油為燃料,進(jìn)行了4種不同直徑油池的燃燒試驗(yàn),試驗(yàn)表明:在有風(fēng)條件下,小直徑油池的燃燒速率隨風(fēng)速的增加呈非單調(diào)變化.李權(quán)威等[5]在低湍流風(fēng)洞試驗(yàn)條件下(0~3.95m/s的縱向通風(fēng)環(huán)境)對(duì)邊長(zhǎng)為4~9cm的正方形乙醇池火進(jìn)行燃燒試驗(yàn),試驗(yàn)表明:油池接受熱反饋的量和通風(fēng)引起的對(duì)流蒸發(fā)是影響池火燃燒速率的原因.Kuang等[6]對(duì)交叉流中的池火燃燒速率進(jìn)行試驗(yàn),當(dāng)橫流風(fēng)速較小時(shí),燃燒速率隨橫流風(fēng)速單調(diào)增大.在航空發(fā)動(dòng)機(jī)領(lǐng)域,超燃沖壓發(fā)動(dòng)機(jī)燃燒室內(nèi)的燃燒就涉及到射流與已燃火焰的相互作用,其中燃料混合、火焰的穩(wěn)定等問(wèn)題都是研究人員感興趣的領(lǐng)域[7].在鍋爐燃燒方面,馬啟磊等[8]在不同二次風(fēng)風(fēng)門(mén)開(kāi)度工況下對(duì)一臺(tái)600MW超臨界前后墻旋流對(duì)沖鍋爐進(jìn)行數(shù)值模擬,模擬表明:風(fēng)門(mén)開(kāi)度越小,爐膛出口NOx質(zhì)量分?jǐn)?shù)越大,鍋爐燃燒效率下降.宋景慧等[9]發(fā)現(xiàn)在總二次風(fēng)風(fēng)量不變情況下,隨著燃盡風(fēng)風(fēng)量占二次風(fēng)總風(fēng)量比例的增大,爐膛出口截面上NOx排放質(zhì)量分?jǐn)?shù)逐漸降低.

    而在內(nèi)燃機(jī)高壓直噴領(lǐng)域的不同之處在于氣流對(duì)火焰是邊射流、邊混合和邊燃燒的模式,學(xué)者們對(duì)高壓射流與火焰之間的關(guān)系做出了一系列研究.Li等[10]研究了不同噴射參數(shù)對(duì)先導(dǎo)點(diǎn)火直噴天然氣發(fā)動(dòng)機(jī)性能和燃燒特性的影響,發(fā)現(xiàn)可通過(guò)延緩和縮短噴射時(shí)間、縮短柴油噴射脈沖寬度和增加噴射壓力來(lái)提高最大放熱率.Lei等[11-12]開(kāi)展高壓甲烷射流對(duì)層流火焰作用的試驗(yàn),得出了高壓甲烷射流一旦沖擊預(yù)混火焰,火焰環(huán)境的壓力和溫度會(huì)改變,并且甲烷噴射時(shí)刻與點(diǎn)火時(shí)刻的時(shí)間間隔不同則火焰的形成演變過(guò)程就不同,隨著時(shí)間間隔的增加,預(yù)混層流火焰等效半徑增大,這樣更能成功地引燃甲烷射流,后以湍流燃燒模式繼續(xù)發(fā)展.王濤[13]對(duì)直噴天然氣發(fā)動(dòng)機(jī)缸內(nèi)火焰面發(fā)展進(jìn)行數(shù)值模擬研究,得出燃料噴射定時(shí)通過(guò)影響點(diǎn)火時(shí)刻缸內(nèi)混合氣當(dāng)量比分布情況,來(lái)決定缸內(nèi)火焰?zhèn)鞑コ跗诘闹鲗?dǎo)燃燒形式和火焰?zhèn)鞑ニ俣龋M(jìn)而影響缸內(nèi)火焰面面積的增大速度,最終影響發(fā)動(dòng)機(jī)整體的燃燒性能.以上研究均表明,氣流對(duì)火焰的燃燒速率和效率、傳播速度及火焰形態(tài)均有明顯影響.內(nèi)燃機(jī)內(nèi)噴射條件不同,射流對(duì)預(yù)混火焰的影響不同,并且均從宏觀的角度來(lái)分析這一變化,微觀上渦量也可用來(lái)分析流場(chǎng)中的氣體流動(dòng).

    以渦量作為觀測(cè)量來(lái)描述整個(gè)流場(chǎng)中的旋渦運(yùn)動(dòng)也在多個(gè)領(lǐng)域應(yīng)用[14].渦是湍流的一種基本結(jié)構(gòu),理清渦的形成和研究如何采用可行的控制手段對(duì)渦旋起主導(dǎo)作用的流動(dòng)實(shí)施有效控制是重要的[15].蔣時(shí)澤等[16]發(fā)現(xiàn)半球陣列式渦流發(fā)生器與水平線的夾角越大,尾跡區(qū)域渦的大小和分布受到前緣區(qū)域渦的影響越大.雷越[17]發(fā)現(xiàn)飛機(jī)起飛時(shí)側(cè)風(fēng)條件下,來(lái)流速度低、離地間隙小的情況下易生成地面渦,其吸入進(jìn)氣道后會(huì)造成總壓損失.潘衛(wèi)軍等[18]研究表明,側(cè)風(fēng)速度對(duì)尾渦的耗散作用并不呈正比.王維軍等[15]通過(guò)渦動(dòng)力學(xué)診斷可以準(zhǔn)確捕捉到離心泵葉輪、蝸殼中渦結(jié)構(gòu)的變化,給離心泵渦控制、流動(dòng)減阻提供理論依據(jù).郭運(yùn)珍[19]研究得出,在螺旋進(jìn)氣道上方加裝1/4尺寸的擋板時(shí),表現(xiàn)出了較好的綜合效果.因此,對(duì)于進(jìn)一步分析直噴天然氣發(fā)動(dòng)機(jī)缸內(nèi)氣體流動(dòng),提高其燃燒效率來(lái)說(shuō),研究渦量是有必要的,且在此方面的研究也較少.

    直噴天然氣發(fā)動(dòng)機(jī)缸內(nèi)燃燒過(guò)程存在一個(gè)基本問(wèn)題就是射流與預(yù)混火焰的相互作用.為了探索此問(wèn)題,前期試驗(yàn)發(fā)現(xiàn)了當(dāng)高壓甲烷射流作用于甲烷預(yù)混火焰時(shí),存在一種現(xiàn)象,即甲烷預(yù)混火焰在射流作用下繼續(xù)膨脹發(fā)展,并引燃甲烷射流產(chǎn)生噴射火焰.整個(gè)燃燒過(guò)程中,彈內(nèi)層流火焰和湍流火焰并存[11].并且前期也討論了彈內(nèi)湍流燃燒形成的原因,即高壓射流擾動(dòng)產(chǎn)生的渦量為湍流燃燒的形成提供適宜的環(huán)境,當(dāng)量比適宜燃燒(先傳質(zhì));當(dāng)射流驅(qū)動(dòng)火焰表面到達(dá)這一區(qū)域時(shí),再加上合適的溫度(再傳熱)形成湍流燃燒[20].渦量是度量旋渦運(yùn)動(dòng)的重要參數(shù),前期均未更深一步討論彈內(nèi)其場(chǎng)的分布及其對(duì)火焰有何影響.由于試驗(yàn)無(wú)法直接獲得彈內(nèi)渦量場(chǎng)的變化,筆者從仿真的角度進(jìn)一步分析高壓甲烷射流沖擊預(yù)混火焰過(guò)程中動(dòng)量參數(shù)渦量的分布規(guī)律以及其對(duì)火焰的影響;實(shí)現(xiàn)對(duì)自由射流模式下彈內(nèi)氣體流動(dòng)以及預(yù)混火焰模式和射流燃燒模式下彈內(nèi)氣體燃燒過(guò)程的有效預(yù)測(cè),并對(duì)控制內(nèi)燃機(jī)總體工作性能的優(yōu)化具有重要意義.

    1 理論分析及數(shù)值方法

    1.1 渦量的定義及計(jì)算

    渦和渦量是流體力學(xué)中兩個(gè)最基本的概念,渦指的是具有渦旋運(yùn)動(dòng)的流體區(qū)域,而渦量是描述漩渦運(yùn)動(dòng)的重要物理量[14],其量化了流體在特定區(qū)域的旋轉(zhuǎn)運(yùn)動(dòng),數(shù)學(xué)上,速度矢量的旋度產(chǎn)生渦量,渦量的計(jì)算見(jiàn)公式(1).

    式中:?代表哈密頓算子;i、j和k表示x、y和z三個(gè)基向量;xi、xj和xk表示x、y和z軸分量;v為流體速度矢量;υi、υj和υk分別表示速度的x、y和z軸分量;收斂計(jì)算渦量的大小為,見(jiàn)公式(2).

    1.2 模型建立及驗(yàn)證

    筆者基于前期完成的高壓甲烷射流沖擊預(yù)燃火焰試驗(yàn)的結(jié)果[11-12]建立數(shù)值模型,試驗(yàn)系統(tǒng)主要由定容燃燒彈裝置、進(jìn)/排氣系統(tǒng)、高速紋影系統(tǒng)和同步控制系統(tǒng)組成,如圖1所示.在射流沖擊火焰試驗(yàn)過(guò)程中,點(diǎn)火針先點(diǎn)燃預(yù)噴甲烷,經(jīng)過(guò)噴射延時(shí)1ms后,甲烷開(kāi)始噴射.

    圖1 試驗(yàn)系統(tǒng)原理示意Fig.1 Schematic diagram of the test system

    定容燃燒彈內(nèi)部容積為圓柱體,點(diǎn)火針、噴嘴均位于其軸線上,且噴嘴出口距點(diǎn)火針上方40mm.為觀察氣體流動(dòng)及整個(gè)燃燒過(guò)程(3ms內(nèi)),且保證彈壁對(duì)內(nèi)部場(chǎng)無(wú)影響,仿真計(jì)算區(qū)域是一個(gè)直徑d1=100mm、高度l1=100mm的圓柱體(與定容燃燒彈同軸),噴嘴出口位于其上表面,如圖1中紅色虛框線所示;建立的三維模型見(jiàn)圖2a,噴孔直徑d0=1mm,高度l0=10mm.

    圖2 模型示意Fig.2 Schematic of the model

    對(duì)比分析自由射流模式、預(yù)混燃燒模式和射流燃燒模式3種模式,研究射流沖擊火焰過(guò)程中渦量的變化特性.射流方向?yàn)閆軸正方向,坐標(biāo)原點(diǎn)(0,0)設(shè)置在噴嘴出口處,如圖2b所示.表1為3種模式的參數(shù)設(shè)置.設(shè)置點(diǎn)火時(shí)刻(預(yù)混燃燒模式和射流燃燒模式)為時(shí)間原點(diǎn),甲烷噴射延時(shí)為1ms(射流燃燒模式、自由射流模式中甲烷噴射時(shí)刻均為1ms).根據(jù)射流燃燒過(guò)程中甲烷射流發(fā)展過(guò)程定義不同時(shí)刻為t1、tin、tc、tout和t2;而在自由射流模式和預(yù)混燃燒模式中其僅代指對(duì)應(yīng)時(shí)刻,無(wú)其他含義.

    表1 3種模式參數(shù)設(shè)置Tab.1 Parameter setting of three modes

    仿真模型網(wǎng)格進(jìn)行了局部加密.一是根據(jù)實(shí)際情況進(jìn)行局部固定加密,由于噴嘴內(nèi)部和噴嘴出口10mm處氣體流速大,與周?chē)鷼怏w流速存在較大速度梯度,因而將其進(jìn)行6級(jí)圓臺(tái)固定加密;采用source/sink模型以0.02J能量點(diǎn)火的方式進(jìn)行點(diǎn)火,火焰前鋒處存在較大溫度梯度,將其進(jìn)行5級(jí)球形固定加密;二是自適應(yīng)加密,高壓射流與預(yù)混火焰的相互作用過(guò)程中彈內(nèi)的溫度場(chǎng)、速度場(chǎng)和壓力場(chǎng)有較大變化,因此,以溫度、速度和壓力為指標(biāo)對(duì)彈內(nèi)流場(chǎng)進(jìn)行5級(jí)自適應(yīng)加密.兩種加密方法的使用在節(jié)約計(jì)算資源的同時(shí)提高了計(jì)算精度.

    求解計(jì)算時(shí),對(duì)物理模型的燃燒模型進(jìn)行簡(jiǎn)化假設(shè),采用CONVERGE中SAGE燃燒模型,湍流模型采用雷諾平均納維斯托克斯(RANS)與重整化群(RNG k-ξ)相結(jié)合的方式;定容燃燒彈內(nèi)的氣體均假設(shè)為理想氣體;定容燃燒彈彈壁假設(shè)為絕熱;整個(gè)過(guò)程不考慮熱輻射.

    反應(yīng)機(jī)理采用GRI3.0版本機(jī)理[21],其包含53個(gè)組分、325個(gè)基元反應(yīng),能夠較吻合地模擬出不同燃燒模式下的甲烷燃燒.

    圖3顯示網(wǎng)格無(wú)關(guān)性結(jié)果,分別選取2、4和8mm的基本網(wǎng)格尺寸,并將3種基本網(wǎng)格尺寸下層流火焰半徑計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)做比較;考慮到計(jì)算準(zhǔn)確性和時(shí)間的低成本,其基本尺寸選為4mm.

    圖3 網(wǎng)格無(wú)關(guān)性驗(yàn)證Fig.3 Verification of grid independence

    前期對(duì)富燃工況下高壓射流沖擊預(yù)燃火焰進(jìn)行了數(shù)值模擬,其中建立的仿真模型得到較好的驗(yàn)證,將仿真得到的層流火焰半徑和貫穿距分別與試驗(yàn)數(shù)據(jù)相比,相對(duì)誤差保持在5%左右,量級(jí)基本一致,并且以甲烷密度為標(biāo)準(zhǔn)對(duì)仿真流場(chǎng)形態(tài)變化與試驗(yàn)做了對(duì)比,兩者得到的流場(chǎng)形態(tài)相似,具有較高計(jì)算精度,本次仿真繼續(xù)使用之前研究模型[20].

    2 結(jié)果與討論

    射流燃燒模式下,高壓射流對(duì)預(yù)混火焰具有較強(qiáng)影響,分別分析自由射流模式和預(yù)混燃燒模式的動(dòng)量參數(shù)變化,并詳細(xì)分析射流燃燒模式下,高壓射流沖擊預(yù)混火焰過(guò)程中渦量變化特征.

    2.1 自由射流模式與預(yù)混燃燒模式

    自由射流模式是只向充有甲烷空氣混合氣(當(dāng)量比為1)的定容燃燒彈內(nèi)噴射高壓甲烷且無(wú)點(diǎn)火;預(yù)混燃燒模式是充有甲烷空氣混合氣(當(dāng)量比為1)的定容燃燒彈內(nèi)只點(diǎn)火點(diǎn)燃預(yù)混氣體,并且無(wú)氣體噴射.圖4為自由射流模式、預(yù)混燃燒模式下的速度.自由射流模式下高壓甲烷射流速度沿Z軸方向不斷下降,這是因?yàn)樯淞魇艿奖尘皻怏w的阻力,使得甲烷自噴嘴噴出后貫穿速度快速下降后逐漸趨于平穩(wěn);預(yù)混火焰速度基本保持不變,約為5.4m/s.

    圖4 自由射流模式和預(yù)混燃燒模式下的速度Fig.4Velocity under conditions of free jet and premixed combustion mode

    圖5為自由射流模式下t1、tin、tc、tout和t2時(shí)刻定容燃燒彈內(nèi)的渦量及速度矢量.氣體射流形成的主渦區(qū)具有對(duì)稱(chēng)性,射流前鋒外圍存在尺度不一的小渦區(qū),不同時(shí)刻其分布的位置也不盡相同.高壓射流噴入后其前端在擠壓背景氣體的同時(shí),噴嘴入口處由于初始流束速度較大,壓力較小,此處與噴嘴遠(yuǎn)端區(qū)域存在壓差,進(jìn)而產(chǎn)生卷吸作用,擠壓和卷吸的共同作用使氣體以Z軸為軸心向四周擴(kuò)散,形成渦量為正的旋渦(橢圓標(biāo)記處),并且兩者的共同作用也會(huì)在射流前鋒與Z軸的交點(diǎn)處形成分叉,分叉會(huì)隨著時(shí)間而增大.由此可得出彈內(nèi)渦量場(chǎng)的變化與速度場(chǎng)具有一定的關(guān)系.

    圖5 自由射流模式渦量及速度矢量Fig.5Vorticity cloud diagram and velocity vector diagram of free jet model

    圖6為自由射流模式軸線距離上截面的平均渦量.t1=1.2ms時(shí),射入時(shí)間較短,受射流影響渦量在噴嘴出口10mm處具有較高水平,且此時(shí)的射流對(duì)噴嘴遠(yuǎn)端氣體沒(méi)有影響.因此,30~100mm渦量幾乎為0;隨著時(shí)間的增加,高壓射流對(duì)彈內(nèi)氣體的影響范圍變大,表現(xiàn)在tin=1.4ms時(shí)30mm之前的區(qū)域渦量較大;tc=1.6ms、tout=1.8ms和t2=3.0ms時(shí)分別在40、50和70mm之前的區(qū)域渦量逐漸變大,并且存在峰值,此峰值會(huì)隨著時(shí)間的增加而變得更大.發(fā)展一段時(shí)間后,渦量出現(xiàn)峰值的原因是氣體逆流(Z軸負(fù)方向)產(chǎn)生旋渦.

    圖6 自由射流模式軸線距離上截面的平均渦量Fig.6 Mean vorticity of the free jet mode at the distance from the axis to the cross section

    圖7為自由射流模式軸線距離上截面的平均壓力.射流噴入時(shí),高速氣體會(huì)擠壓彈內(nèi)背景氣體,以貫穿距處為界限,其后的壓力會(huì)有所增加;由于壓差的存在使部分射流氣體在此處的流動(dòng)方向發(fā)生偏轉(zhuǎn),即向Z軸負(fù)方向流動(dòng),再加上卷吸作用.因此,氣體會(huì)在貫穿距這一分界線的前端形成旋渦,進(jìn)而渦量在貫穿距前端會(huì)增加,出現(xiàn)峰值.

    圖7 自由射流模式軸線距離上截面的平均壓力Fig.7 Average pressure of free jet mode at the distance from the axis to the section

    圖8為預(yù)混燃燒火焰不同時(shí)刻渦量及速度矢量.初始時(shí)刻,預(yù)混火焰呈球形向四周發(fā)展,表面較光滑,火焰前鋒面外部的未燃區(qū)和內(nèi)部的已燃區(qū)渦量均較小,約為50s-1;而只有火焰前鋒面這一化學(xué)反應(yīng)區(qū)渦量較大,約為300s-1,其值是彈內(nèi)其他區(qū)域的6倍;隨著預(yù)混燃燒時(shí)間的增加,光滑球形火焰面上會(huì)產(chǎn)生裂紋和胞狀結(jié)構(gòu),表面上褶皺的形成使此處渦區(qū)擴(kuò)大,渦量有所增加,這一規(guī)律可從圖8上直觀地看出.

    圖9為預(yù)混燃燒模式軸線距離上截面平均渦量,可定量表達(dá)出渦量隨時(shí)間及軸線位置的變化規(guī)律,圖中d1、din、dc、dout和d2分別表示對(duì)應(yīng)時(shí)刻預(yù)混火焰直徑.反應(yīng)時(shí)間越短,球形火焰半徑越小,火焰體積和表面積也均處于較小狀態(tài),截面上火焰前鋒處反應(yīng)區(qū)較小,因而平均渦量值越小,且火焰前鋒處渦量值與其他截面處相差不大.火焰前鋒隨著時(shí)間的增加而增大,使截面上的平均渦量值增大,t2時(shí),火焰前鋒處與其他截面相差較明顯,約為15s-1,且渦量值以點(diǎn)火處(Z=40mm)為中心,呈對(duì)稱(chēng)分布.

    圖9 預(yù)混燃燒模式軸線距離上截面的平均渦量Fig.9Average vorticity of premixed combustion mode from axis to cross section

    2.2 射流燃燒模式

    圖10顯示射流燃燒的發(fā)展過(guò)程中OH基質(zhì)量分?jǐn)?shù)的變化.預(yù)混火焰先形成,1.0ms后甲烷開(kāi)始噴射,射流前鋒面向前發(fā)展,沖擊預(yù)混火焰迎風(fēng)面,從火焰背風(fēng)面沖出離開(kāi)火焰,期間火焰不斷發(fā)展擴(kuò)大.整個(gè)射流作用預(yù)混火焰的過(guò)程分為射流前鋒面未接觸火焰階段(0~1.3ms)、射流前鋒面進(jìn)入火焰階段(1.4~1.8ms)以及射流前峰面離開(kāi)火焰階段(1.9~3.0ms).

    圖10 甲烷質(zhì)量分?jǐn)?shù)為0~0.25的三維火焰面X=0截面上OH基質(zhì)量分?jǐn)?shù)Fig.10 Mass fraction of the OH cloud map on the X=0 section of the three-dimensional flame surface in the methane mass fraction of 0—0.25

    圖11為射流燃燒模式下定容燃燒彈內(nèi)的渦量等值線及速度矢量.射流噴入后,定容燃燒彈內(nèi)的氣體流動(dòng)與自由射流模式一致,均在射流貫穿距的前端形成渦.在射流前鋒面進(jìn)入火焰階段,射流作用于火焰,渦量增加.

    圖11 射流燃燒模式X=0截面上選取時(shí)刻渦量等值線及速度矢量Fig.11 Vorticity contour line and velocity vector diagram at selected time on the section of jet combustion mode X=0

    不同階段的射流貫穿距l(xiāng)jf定義如表2所示,可以根據(jù)ljf求出不同階段的速度.圖12為射流燃燒、自由射流和預(yù)混燃燒3個(gè)模式的速度對(duì)比,在射流前鋒面未接觸火焰階段和射流前鋒面離開(kāi)火焰階段,射流燃燒和自由射流模式的平均貫穿速度幾乎相等;但射流前鋒面進(jìn)入火焰階段射流燃燒的平均貫穿速度約是預(yù)混燃燒火焰面發(fā)展速度的10倍,而自由射流平均貫穿速度約是預(yù)混燃燒火焰面發(fā)展速度的8倍,射流燃燒的平均貫穿速度比自由射流的大11.4m/s,且大于自由射流平均貫穿速度與預(yù)混燃燒火焰面速度之和.這說(shuō)明射流前鋒面進(jìn)入火焰階段速度的增加受動(dòng)量參數(shù)渦量變化的影響.

    表2 貫穿距數(shù)據(jù)定義Tab.2 Definition of penetration distance

    圖12 3個(gè)模式速度對(duì)比Fig.12 Comparison of velocity in three modes

    圖13a為不同Z軸位置(20、40、60和80mm)處各工況截面平均渦量值隨時(shí)間的變化,圖13b為40mm處平均渦量值隨時(shí)間的變化,圖13c為t=3.0ms時(shí)渦量在軸線位置上的變化.在射流未接觸預(yù)混火焰時(shí)(t<1.3ms),自由射流貫穿距與射流燃燒貫穿距重合.因此,兩種工況的平均渦量相等,圖13a中表現(xiàn)為射流前鋒面未接觸火焰階段上兩種模式的渦量曲線重合.

    圖13 不同截面位置的平均渦量值Fig.13 Average vorticity values at different cross-sectional locations

    由圖13a可以看出,射流燃燒模式截面平均渦量大于自由射流模式下的值.在射流前鋒面進(jìn)入火焰階段,射流燃燒模式4個(gè)截面處的平均渦量值分別在1.3、1.3、1.6和1.7ms開(kāi)始大于自由射流模式,軸線距離越遠(yuǎn)的區(qū)域,甲烷射流對(duì)其擾動(dòng)開(kāi)始的時(shí)間越晚.并且在20、60和80mm處兩模式平均渦量值相差較40mm處的小.在40mm處、1.6ms時(shí),射流燃燒模式平均渦量值與自由射流模式差值最大,約為260s-1,是此時(shí)此截面預(yù)混燃燒平均渦量值6s-1的43倍,如圖13b所示.射流燃燒模式,射流前鋒面進(jìn)入火焰階段渦量的大幅增加導(dǎo)致其平均貫穿速度的增加,射流燃燒不是自由射流與預(yù)混燃燒物理動(dòng)量上的簡(jiǎn)單耦合.

    在射流前鋒面離開(kāi)火焰階段,隨著時(shí)間的增加,射流流束貫穿速度減小,射流燃燒模式火焰前峰面向四周擴(kuò)大,加大其周?chē)鷼怏w的擾動(dòng),因而各截面上平均渦量值均呈增加趨勢(shì);但自由射流模式下、同一時(shí)刻,渦量隨著軸線距離先增加達(dá)到峰值,隨后驟然降低.當(dāng)t=3ms時(shí),射流燃燒和自由射流模式的平均渦量最大值均出現(xiàn)在60mm處;從圖13c可以看出,60mm截面處射流燃燒平均渦量值約為自由射流模式平均渦量值的2倍.同樣佐證了射流燃燒不是自由射流與預(yù)混燃燒物理動(dòng)量上的簡(jiǎn)單耦合這一規(guī)律.

    3 結(jié) 論

    射流沖擊火焰的全過(guò)程被分為3個(gè)階段,通過(guò)詳細(xì)分析自由射流模式、預(yù)混燃燒模式和射流燃燒模式3個(gè)階段的渦量,得出以下結(jié)論:

    (1) 自由射流模式下,同一時(shí)刻,渦量隨著軸線距離先增加達(dá)到峰值,隨后驟然降低;其出現(xiàn)峰值的原因是射流主體擠壓遠(yuǎn)端背景氣體使兩者之間產(chǎn)生壓差,射流氣體回流產(chǎn)生旋渦.

    (2) 預(yù)混燃燒模式下,只有球形火焰前鋒處的渦量處于較高水平,其值約為300s-1,是彈內(nèi)其他區(qū)域的6倍.

    (3) 射流燃燒模式射流前鋒面進(jìn)入火焰階段的貫穿速度約是預(yù)混燃燒火焰面發(fā)展速度的10倍,而自由射流平均貫穿速度約是預(yù)混燃燒火焰面發(fā)展速度的8倍;在40mm處、1.6ms時(shí),射流燃燒模式平均渦量值與自由射流模式差值最大,約為260s-1,是此時(shí)此截面預(yù)混燃燒模式平均渦量值6s-1的43倍,此階段速度的增加受動(dòng)量參數(shù)渦量的影響;在射流前鋒面進(jìn)入火焰階段和射流前鋒面離開(kāi)火焰兩個(gè)階段,射流燃燒不是自由射流與預(yù)混燃燒物理動(dòng)量上的簡(jiǎn)單耦合.

    猜你喜歡
    渦量量值前鋒
    多元向量值區(qū)域和加權(quán)風(fēng)險(xiǎn)值
    含沙空化對(duì)軸流泵內(nèi)渦量分布的影響
    廣州市番禺區(qū)石碁鎮(zhèn)前鋒小學(xué)作品集
    ViVi美眉(2021年10期)2021-01-25 23:51:49
    基于QAR數(shù)據(jù)的碳當(dāng)量值適航符合性驗(yàn)證方法
    帶有中心值的量值的公差表示
    山東冶金(2018年5期)2018-11-22 05:12:28
    自由表面渦流動(dòng)現(xiàn)象的數(shù)值模擬
    籃球的由來(lái)
    旋量值函數(shù)的Plemelj公式
    航態(tài)對(duì)大型船舶甲板氣流場(chǎng)的影響
    The application of numerical simulation of delta wing with blunt leading edge using RANS/LES hybrid method
    日韩av在线大香蕉| 亚洲国产毛片av蜜桃av| 日本欧美视频一区| 国产精品av久久久久免费| 亚洲国产精品成人综合色| 一进一出抽搐动态| 亚洲在线自拍视频| 亚洲av五月六月丁香网| 亚洲国产毛片av蜜桃av| 狠狠狠狠99中文字幕| 久久这里只有精品19| svipshipincom国产片| 久久精品亚洲熟妇少妇任你| 亚洲少妇的诱惑av| 精品久久久久久久毛片微露脸| 波多野结衣一区麻豆| 天堂影院成人在线观看| 成人特级黄色片久久久久久久| 日韩精品免费视频一区二区三区| 亚洲国产精品成人综合色| 老司机午夜福利在线观看视频| 99久久99久久久精品蜜桃| 悠悠久久av| 国产精品香港三级国产av潘金莲| 久久人妻熟女aⅴ| 一级作爱视频免费观看| 久久青草综合色| 国产精品亚洲一级av第二区| 少妇被粗大的猛进出69影院| 在线观看午夜福利视频| 久久久久精品国产欧美久久久| 久久久久国产精品人妻aⅴ院| 男女午夜视频在线观看| 国产伦人伦偷精品视频| 午夜久久久在线观看| 日韩国内少妇激情av| 日本五十路高清| 亚洲中文av在线| 一本久久中文字幕| 91老司机精品| 国产亚洲av高清不卡| 一区福利在线观看| 亚洲国产精品成人综合色| 欧美成人免费av一区二区三区| 搡老熟女国产l中国老女人| 亚洲国产欧美日韩在线播放| 成人免费观看视频高清| 亚洲成a人片在线一区二区| 免费一级毛片在线播放高清视频 | 亚洲人成77777在线视频| 亚洲人成电影观看| 99精品欧美一区二区三区四区| 久久久久国内视频| 久久国产精品人妻蜜桃| 动漫黄色视频在线观看| 免费少妇av软件| 亚洲九九香蕉| 久久久国产成人精品二区| 精品久久久精品久久久| 国产精品久久久av美女十八| av视频在线观看入口| 精品无人区乱码1区二区| 一级毛片女人18水好多| 丝袜在线中文字幕| 国产黄a三级三级三级人| 亚洲五月婷婷丁香| 脱女人内裤的视频| 精品人妻1区二区| or卡值多少钱| 久久久国产成人精品二区| 熟妇人妻久久中文字幕3abv| 在线视频色国产色| 日韩欧美免费精品| 久久精品国产99精品国产亚洲性色 | 亚洲欧洲精品一区二区精品久久久| 亚洲色图av天堂| 国产主播在线观看一区二区| 久久久久久久精品吃奶| 国产片内射在线| 久久国产精品影院| 99热只有精品国产| 91国产中文字幕| 精品一区二区三区av网在线观看| 悠悠久久av| 国产一区二区三区综合在线观看| 免费在线观看亚洲国产| 熟妇人妻久久中文字幕3abv| 国产主播在线观看一区二区| 禁无遮挡网站| 国内精品久久久久精免费| 妹子高潮喷水视频| 亚洲av电影在线进入| 精品无人区乱码1区二区| 精品卡一卡二卡四卡免费| 日韩欧美国产一区二区入口| 一夜夜www| 久久草成人影院| 亚洲一区中文字幕在线| 一进一出好大好爽视频| 一本久久中文字幕| 精品高清国产在线一区| 日韩成人在线观看一区二区三区| 色播亚洲综合网| 午夜久久久在线观看| 变态另类成人亚洲欧美熟女 | 国产激情欧美一区二区| 久久人人精品亚洲av| 日韩视频一区二区在线观看| 人妻久久中文字幕网| 日韩精品中文字幕看吧| 宅男免费午夜| 热99re8久久精品国产| 丝袜美足系列| 久久狼人影院| 母亲3免费完整高清在线观看| 国产激情久久老熟女| 国语自产精品视频在线第100页| 1024香蕉在线观看| 人人妻人人澡人人看| 丝袜美足系列| 少妇被粗大的猛进出69影院| 精品熟女少妇八av免费久了| 一区二区三区国产精品乱码| 亚洲成人精品中文字幕电影| 色精品久久人妻99蜜桃| 久久青草综合色| 99国产极品粉嫩在线观看| av在线天堂中文字幕| 国产区一区二久久| 中文亚洲av片在线观看爽| 日韩欧美三级三区| 性少妇av在线| av免费在线观看网站| 男女做爰动态图高潮gif福利片 | 国产精品99久久99久久久不卡| 两个人免费观看高清视频| 中亚洲国语对白在线视频| 久久热在线av| 日韩免费av在线播放| 久久久国产精品麻豆| e午夜精品久久久久久久| 欧美绝顶高潮抽搐喷水| 午夜免费观看网址| 午夜精品久久久久久毛片777| 91在线观看av| 婷婷六月久久综合丁香| 亚洲色图 男人天堂 中文字幕| 欧美最黄视频在线播放免费| 亚洲第一av免费看| 欧美乱妇无乱码| 国产精品电影一区二区三区| 嫩草影视91久久| 精品国产乱码久久久久久男人| 狂野欧美激情性xxxx| 日韩精品中文字幕看吧| 亚洲精品美女久久av网站| x7x7x7水蜜桃| 国产精品99久久99久久久不卡| 日韩欧美三级三区| aaaaa片日本免费| 不卡av一区二区三区| 精品一品国产午夜福利视频| 黄色片一级片一级黄色片| 韩国av一区二区三区四区| 久久国产精品男人的天堂亚洲| 这个男人来自地球电影免费观看| 亚洲国产欧美网| 久久国产精品男人的天堂亚洲| 丝袜人妻中文字幕| 欧美最黄视频在线播放免费| 国产又色又爽无遮挡免费看| 女警被强在线播放| 9热在线视频观看99| 91麻豆精品激情在线观看国产| 久久香蕉激情| 一边摸一边抽搐一进一出视频| 宅男免费午夜| 久久久久国产一级毛片高清牌| 黄色成人免费大全| 在线天堂中文资源库| 十分钟在线观看高清视频www| 老熟妇乱子伦视频在线观看| 亚洲av美国av| 校园春色视频在线观看| 亚洲精品中文字幕在线视频| 99国产精品一区二区蜜桃av| 国产成人av激情在线播放| 男人舔女人的私密视频| 亚洲国产精品成人综合色| 国产在线精品亚洲第一网站| 午夜影院日韩av| 1024视频免费在线观看| 亚洲欧美激情在线| 精品久久久久久久毛片微露脸| www.精华液| 色哟哟哟哟哟哟| 亚洲狠狠婷婷综合久久图片| 久久影院123| 自线自在国产av| 国产伦人伦偷精品视频| 欧美最黄视频在线播放免费| 色综合亚洲欧美另类图片| 欧美日韩福利视频一区二区| 一级黄色大片毛片| 亚洲精品美女久久久久99蜜臀| 国产成人一区二区三区免费视频网站| 一区二区三区激情视频| 欧美成人性av电影在线观看| 免费在线观看日本一区| 男女做爰动态图高潮gif福利片 | 中文亚洲av片在线观看爽| 91成年电影在线观看| 久久草成人影院| 国产一区二区三区视频了| 大陆偷拍与自拍| 亚洲国产中文字幕在线视频| 亚洲精品在线观看二区| 桃红色精品国产亚洲av| 色婷婷久久久亚洲欧美| 亚洲精品粉嫩美女一区| 一卡2卡三卡四卡精品乱码亚洲| 99久久久亚洲精品蜜臀av| 国产精品精品国产色婷婷| www.精华液| 自线自在国产av| 电影成人av| 看免费av毛片| 精品国产一区二区三区四区第35| 国产麻豆成人av免费视频| av免费在线观看网站| 女性生殖器流出的白浆| √禁漫天堂资源中文www| 大型黄色视频在线免费观看| 一级片免费观看大全| 免费在线观看日本一区| 久久精品国产清高在天天线| 国产精品亚洲av一区麻豆| 夜夜看夜夜爽夜夜摸| 久久久国产精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 亚洲伊人色综图| 久久精品成人免费网站| 亚洲va日本ⅴa欧美va伊人久久| 国产成年人精品一区二区| 亚洲人成伊人成综合网2020| 亚洲精品国产区一区二| 日本 欧美在线| 精品高清国产在线一区| 1024香蕉在线观看| 成人三级做爰电影| av中文乱码字幕在线| 一本久久中文字幕| 国产精品98久久久久久宅男小说| 欧美日韩乱码在线| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 美女 人体艺术 gogo| 国产精品1区2区在线观看.| 视频区欧美日本亚洲| 久久久久精品国产欧美久久久| 欧美日本亚洲视频在线播放| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 国产三级黄色录像| xxx96com| 欧美在线黄色| www国产在线视频色| 国产成人欧美| 一区二区三区精品91| 无限看片的www在线观看| 天堂动漫精品| 两个人看的免费小视频| 免费少妇av软件| 精品国产一区二区久久| 国产一级毛片七仙女欲春2 | 欧美大码av| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 国产精品1区2区在线观看.| 成在线人永久免费视频| 看黄色毛片网站| 色综合亚洲欧美另类图片| 久久精品国产99精品国产亚洲性色 | 色综合亚洲欧美另类图片| 桃色一区二区三区在线观看| 亚洲伊人色综图| 青草久久国产| 久久国产精品影院| 人成视频在线观看免费观看| 香蕉久久夜色| 嫩草影院精品99| 人人澡人人妻人| 黄色女人牲交| 亚洲男人的天堂狠狠| 嫩草影视91久久| 精品欧美国产一区二区三| 欧美不卡视频在线免费观看 | 91精品三级在线观看| 熟女少妇亚洲综合色aaa.| 天天躁夜夜躁狠狠躁躁| 国产熟女午夜一区二区三区| 999精品在线视频| 亚洲精华国产精华精| 国产精品乱码一区二三区的特点 | 岛国在线观看网站| 国产亚洲精品久久久久5区| 国内精品久久久久精免费| 久久天堂一区二区三区四区| 熟女少妇亚洲综合色aaa.| 午夜福利18| 一区在线观看完整版| 亚洲av五月六月丁香网| 啪啪无遮挡十八禁网站| 韩国av一区二区三区四区| 国产亚洲欧美在线一区二区| 制服人妻中文乱码| 国产主播在线观看一区二区| 国产成人啪精品午夜网站| 国产精品久久久久久亚洲av鲁大| 久久久精品国产亚洲av高清涩受| 免费在线观看影片大全网站| 久久九九热精品免费| 国产精品99久久99久久久不卡| 极品教师在线免费播放| 亚洲一区二区三区不卡视频| 淫秽高清视频在线观看| 一级毛片高清免费大全| 变态另类丝袜制服| 怎么达到女性高潮| 国产一区二区三区综合在线观看| 国产成人av激情在线播放| 91麻豆精品激情在线观看国产| 欧美另类亚洲清纯唯美| 精品国产一区二区三区四区第35| 欧美日韩亚洲国产一区二区在线观看| 国产又爽黄色视频| 精品日产1卡2卡| 亚洲精品在线观看二区| 精品无人区乱码1区二区| 少妇粗大呻吟视频| 别揉我奶头~嗯~啊~动态视频| 999久久久精品免费观看国产| 亚洲五月天丁香| 色婷婷久久久亚洲欧美| а√天堂www在线а√下载| 久久精品国产综合久久久| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲va日本ⅴa欧美va伊人久久| 国产精品爽爽va在线观看网站 | 亚洲自偷自拍图片 自拍| 亚洲av日韩精品久久久久久密| 国产成人av激情在线播放| 亚洲,欧美精品.| 日本在线视频免费播放| 午夜福利一区二区在线看| 欧美中文综合在线视频| 精品国产一区二区三区四区第35| 亚洲国产中文字幕在线视频| 久久国产精品男人的天堂亚洲| 国产单亲对白刺激| 国产成人欧美在线观看| 美女午夜性视频免费| 精品午夜福利视频在线观看一区| 亚洲少妇的诱惑av| www国产在线视频色| 人人妻人人澡欧美一区二区 | 久久精品国产亚洲av香蕉五月| www国产在线视频色| 村上凉子中文字幕在线| 夜夜夜夜夜久久久久| 免费看a级黄色片| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 日韩有码中文字幕| 精品电影一区二区在线| 黄色片一级片一级黄色片| 久久久久久人人人人人| 很黄的视频免费| 日韩精品青青久久久久久| 人人妻人人澡人人看| 午夜福利高清视频| 无人区码免费观看不卡| 日韩欧美国产一区二区入口| 日韩一卡2卡3卡4卡2021年| 12—13女人毛片做爰片一| 亚洲色图 男人天堂 中文字幕| 不卡一级毛片| 99精品久久久久人妻精品| 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| 日本五十路高清| 国产亚洲欧美98| 韩国精品一区二区三区| 国产97色在线日韩免费| 美女大奶头视频| 国产精品亚洲av一区麻豆| 亚洲成人精品中文字幕电影| 大码成人一级视频| 国产主播在线观看一区二区| 亚洲欧洲精品一区二区精品久久久| 亚洲激情在线av| 国产精品一区二区在线不卡| 91精品三级在线观看| 高潮久久久久久久久久久不卡| 天堂动漫精品| 男男h啪啪无遮挡| 90打野战视频偷拍视频| 欧美 亚洲 国产 日韩一| 日本免费a在线| 欧美绝顶高潮抽搐喷水| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲精品一区二区www| 亚洲av电影在线进入| 美国免费a级毛片| 免费观看人在逋| 岛国在线观看网站| 一级,二级,三级黄色视频| 亚洲熟女毛片儿| a在线观看视频网站| 国产一区二区三区视频了| 好男人在线观看高清免费视频 | 女性生殖器流出的白浆| 在线观看舔阴道视频| 男女之事视频高清在线观看| 一边摸一边抽搐一进一出视频| 色精品久久人妻99蜜桃| 色哟哟哟哟哟哟| 欧美日韩福利视频一区二区| 村上凉子中文字幕在线| 日韩欧美在线二视频| 欧美老熟妇乱子伦牲交| 日韩大尺度精品在线看网址 | 色老头精品视频在线观看| 免费搜索国产男女视频| 变态另类成人亚洲欧美熟女 | 久久久久亚洲av毛片大全| 极品人妻少妇av视频| 国产亚洲精品第一综合不卡| 国产亚洲欧美精品永久| 午夜福利,免费看| 精品国内亚洲2022精品成人| 男女下面插进去视频免费观看| 亚洲最大成人中文| 一个人观看的视频www高清免费观看 | 亚洲一区二区三区色噜噜| 69精品国产乱码久久久| 免费看a级黄色片| 中文字幕色久视频| 久久婷婷人人爽人人干人人爱 | or卡值多少钱| 久久精品亚洲熟妇少妇任你| 黄频高清免费视频| 99久久99久久久精品蜜桃| 国产精品日韩av在线免费观看 | 欧美另类亚洲清纯唯美| 人人妻人人澡欧美一区二区 | 久久久国产欧美日韩av| 女人精品久久久久毛片| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频| 久久久久久免费高清国产稀缺| 欧美日韩乱码在线| 午夜日韩欧美国产| 美女高潮到喷水免费观看| 天堂影院成人在线观看| 男女做爰动态图高潮gif福利片 | 老司机深夜福利视频在线观看| 亚洲av美国av| 午夜久久久在线观看| 久久人人97超碰香蕉20202| 色精品久久人妻99蜜桃| 国产伦一二天堂av在线观看| 变态另类丝袜制服| 久久热在线av| 精品国产美女av久久久久小说| 狂野欧美激情性xxxx| 国产精品久久久久久亚洲av鲁大| 中文字幕人妻熟女乱码| 亚洲中文字幕日韩| 母亲3免费完整高清在线观看| 搞女人的毛片| 欧美激情高清一区二区三区| 男男h啪啪无遮挡| 国产一区二区三区综合在线观看| 亚洲avbb在线观看| 在线视频色国产色| 很黄的视频免费| 好男人电影高清在线观看| 成人国语在线视频| 国产91精品成人一区二区三区| 中文字幕久久专区| 叶爱在线成人免费视频播放| 老熟妇乱子伦视频在线观看| 欧美成狂野欧美在线观看| av中文乱码字幕在线| 熟女少妇亚洲综合色aaa.| 亚洲成av片中文字幕在线观看| 久久久久久久久中文| 亚洲精品中文字幕一二三四区| 精品不卡国产一区二区三区| 一进一出抽搐动态| 日韩精品免费视频一区二区三区| 亚洲国产中文字幕在线视频| 美女免费视频网站| 亚洲成国产人片在线观看| 乱人伦中国视频| 日本vs欧美在线观看视频| 久久久久亚洲av毛片大全| 日本vs欧美在线观看视频| 大型黄色视频在线免费观看| 无人区码免费观看不卡| 免费一级毛片在线播放高清视频 | 黄片大片在线免费观看| 久久精品91蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美+亚洲+日韩+国产| av福利片在线| 男人舔女人的私密视频| 日本一区二区免费在线视频| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区| 97人妻精品一区二区三区麻豆 | 青草久久国产| 日韩免费av在线播放| 日本一区二区免费在线视频| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 久久婷婷人人爽人人干人人爱 | 色精品久久人妻99蜜桃| 国产精品亚洲av一区麻豆| 国产又色又爽无遮挡免费看| 女人精品久久久久毛片| 日韩成人在线观看一区二区三区| 精品欧美一区二区三区在线| 午夜a级毛片| 国内久久婷婷六月综合欲色啪| 午夜福利一区二区在线看| 午夜久久久在线观看| av欧美777| 中文字幕最新亚洲高清| 免费av毛片视频| 久久久久久久久免费视频了| 人人妻人人爽人人添夜夜欢视频| 黄色丝袜av网址大全| 欧美午夜高清在线| 最近最新中文字幕大全电影3 | 亚洲自拍偷在线| 久久精品国产99精品国产亚洲性色 | 国产精品久久久久久亚洲av鲁大| 两人在一起打扑克的视频| 精品久久久久久,| 亚洲成a人片在线一区二区| 亚洲少妇的诱惑av| 国产成人av教育| 亚洲精品在线观看二区| 亚洲人成伊人成综合网2020| 亚洲专区中文字幕在线| 在线观看免费午夜福利视频| 欧美日本亚洲视频在线播放| 中出人妻视频一区二区| 成人三级黄色视频| 免费久久久久久久精品成人欧美视频| 国产av一区二区精品久久| 色尼玛亚洲综合影院| 男女床上黄色一级片免费看| 夜夜躁狠狠躁天天躁| 99在线人妻在线中文字幕| 国产亚洲精品av在线| 久久人妻福利社区极品人妻图片| 亚洲最大成人中文| 搡老妇女老女人老熟妇| 欧美日韩瑟瑟在线播放| 狠狠狠狠99中文字幕| 侵犯人妻中文字幕一二三四区| 国产精华一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 日韩有码中文字幕| 久久久久亚洲av毛片大全| 中出人妻视频一区二区| 两个人免费观看高清视频| xxx96com| 人人妻人人爽人人添夜夜欢视频| 午夜免费鲁丝| 久久亚洲精品不卡| 天堂影院成人在线观看| 十八禁网站免费在线| 一级毛片高清免费大全| 国产亚洲精品久久久久5区| 又黄又爽又免费观看的视频| 亚洲五月天丁香| 丝袜在线中文字幕| 女警被强在线播放| 国产精品影院久久| 欧美一级毛片孕妇| 久久亚洲精品不卡| 亚洲情色 制服丝袜| 19禁男女啪啪无遮挡网站| 色播在线永久视频| 高清黄色对白视频在线免费看| 国产精品国产高清国产av| av在线天堂中文字幕| 亚洲色图av天堂| 老鸭窝网址在线观看| 一级黄色大片毛片| 后天国语完整版免费观看| 一区二区三区高清视频在线| 人人妻人人澡人人看| 欧美日本亚洲视频在线播放| 三级毛片av免费| 久久人人精品亚洲av| 熟女少妇亚洲综合色aaa.| 亚洲第一青青草原| 日日爽夜夜爽网站| 亚洲七黄色美女视频| 一边摸一边做爽爽视频免费|