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

    撞擊流反應(yīng)器內(nèi)渦特性研究進(jìn)展

    2021-11-30 07:40:32張建偉高偉峰馮穎張一凡董鑫
    化工進(jìn)展 2021年11期
    關(guān)鍵詞:渦的射流流場

    張建偉,高偉峰,馮穎,張一凡,董鑫

    (沈陽化工大學(xué)機(jī)械與動力工程學(xué)院,遼寧 沈陽 110142)

    在石油化工、煤化工等過程工業(yè)中,時(shí)常伴有混合過程,提高混合效率一直是國內(nèi)外學(xué)者的研究熱點(diǎn)[1]。撞擊流作為一種新型高效的過程強(qiáng)化技術(shù)[2-3],被廣泛用于超細(xì)粉體的制備[4]、煤氣化[5]等工業(yè)生產(chǎn)中,有著優(yōu)于傳統(tǒng)同類設(shè)備的性能。撞擊流反應(yīng)器良好的混合特性與速度場、壓力場、渦量場等特性密切相關(guān)。近年來,關(guān)于流場渦結(jié)構(gòu)特性成為研究熱點(diǎn)?;诖耍疚膶ψ矒袅骰旌闲阅?、渦結(jié)構(gòu)特點(diǎn)以及湍流流場的本征正交分解(POD)分析進(jìn)行綜述,對優(yōu)化撞擊流反應(yīng)器結(jié)構(gòu)、提高反應(yīng)器混合效果奠定了理論基礎(chǔ)。

    1 撞擊流反應(yīng)器內(nèi)結(jié)構(gòu)與混合特性

    撞擊流(impinging stream)是由蘇聯(lián)學(xué)者Elperin[6]最初提出構(gòu)想,其原理如圖1 所示[7]。兩股均相或非均相流體通過噴嘴在反應(yīng)器中心形成高湍動、高物相濃度的撞擊區(qū),由此強(qiáng)化傳質(zhì)傳熱。撞擊流反應(yīng)器內(nèi)流場分為三個(gè)區(qū)域:一是噴嘴到撞擊面之間的射流區(qū);二是兩股相向流動的流體發(fā)生撞擊稱為撞擊區(qū);三是撞擊后流體由軸向轉(zhuǎn)為徑向流動,稱為徑向射流區(qū)[7-9]。其中,撞擊區(qū)內(nèi)流體的速度梯度大,撞擊駐點(diǎn)附近和徑向射流兩側(cè)形成了較大尺度的渦,渦核位置隨時(shí)間而發(fā)生移動,同時(shí)噴嘴附近產(chǎn)生渦量較大的渦。

    圖1 撞擊流反應(yīng)器流場區(qū)域[7]

    在撞擊流反應(yīng)器湍流流場內(nèi)有大量不同尺度的渦旋,大渦旋具有脈動周期長、有效能量大的特點(diǎn)[10]。湍流的混合可分為宏觀混合、介觀混合和微觀混合3個(gè)尺度。在撞擊區(qū),渦旋之間發(fā)生相互作用,較大尺度的渦分裂成小的渦旋,實(shí)現(xiàn)能量傳遞,同時(shí)達(dá)到宏觀混合。在宏觀混合中,以整個(gè)反應(yīng)器為研究對象,為直觀反映混合效果,常用反應(yīng)物轉(zhuǎn)化率[11]、離集指數(shù)Xs[12]、離析度(IOS)[13]等參數(shù)來表征反應(yīng)器的混合效果。如圖2所示[14]的混合過程,流體初始進(jìn)入反應(yīng)器內(nèi)混合尚不完全,如圖2(a)所示,此時(shí)IOS=0.3;隨著混合的持續(xù)進(jìn)行IOS 進(jìn)一步降低,混合效果進(jìn)一步改善,如圖2(b)所示,此刻IOS=0.1;當(dāng)IOS<0.05時(shí)認(rèn)為完全混合,如圖2(c)所示,整個(gè)過程流場逐漸趨于無序。宏觀混合初期,渦“生長”占主導(dǎo)地位,該過程渦的空間尺寸隨著流向距離的增加而增加,邊界處的流體“深入”到軸向射流內(nèi),由于軸向射流的能量輸入使渦的旋轉(zhuǎn)速度增加,促使渦核及附近流體運(yùn)輸?shù)綔u邊界,從而實(shí)現(xiàn)混合。渦把較大的微團(tuán)剪切成更小的微團(tuán),通過小微團(tuán)界面之間的渦流擴(kuò)散,直到達(dá)到Kolmogorov 尺度時(shí)稱為介觀混合[15]。在該尺度下,由于渦旋的運(yùn)輸作用使微團(tuán)改變其固有位置,不同組分的流體進(jìn)行相互纏繞,從而達(dá)到混合目的。同時(shí)介觀混合尺度下渦的融合作用也不可忽略,渦的融合是渦內(nèi)部流體對等組合的過程(渦核之間的融合、邊界與邊界融合),其內(nèi)部流體密度發(fā)生變化,同時(shí)新形成的渦不斷卷吸周圍流體,整個(gè)流場的密度也因此趨于一致。微觀混合是發(fā)生短距離的分子擴(kuò)散,使混合物料分子達(dá)到高度均一性的隨機(jī)分布。在分析微觀混合時(shí),僅考慮渦旋內(nèi)的分子擴(kuò)散,實(shí)際上由于撞擊流內(nèi)小渦旋內(nèi)黏性切應(yīng)力占主導(dǎo)地位,速度梯度大,離散作用不可忽略。圖2(a)為撞擊初期,此時(shí)兩股流體尚未充分混合,流場內(nèi)的渦尺度較大;圖2(b)則是兩股流體撞擊一定時(shí)間后大尺度渦相互剪切、卷吸產(chǎn)生大量的小尺度渦,使混合效果得到初步改善;撞擊流反應(yīng)器內(nèi)的渦由于黏性作用,不斷發(fā)生剪切,使反應(yīng)器內(nèi)的流體更趨近于無序,如圖2(c)所示,此刻認(rèn)為達(dá)到微觀混合[11,15]。

    圖2 撞擊流反應(yīng)器內(nèi)混合過程[14]

    眾多學(xué)者利用不同手段研究實(shí)驗(yàn)參數(shù)對撞擊流反應(yīng)器混合效果的影響,如Lyu 等[16]通過Xs反映T形通道的混合效果,發(fā)現(xiàn)大雷諾數(shù)對混合有較大影響。Zhang 等[17]發(fā)現(xiàn)T 形反應(yīng)器中的吞噬流階段,隨Re的增大,混合效果也更為理想。表1為實(shí)驗(yàn)參數(shù)對混合的影響結(jié)果匯總。

    從表1可知,流速和流速比是影響混合的重要條件,增大流速可以增加流場的湍動程度,減小撞擊面兩側(cè)渦的脫落周期和撞擊駐點(diǎn)振蕩周期,進(jìn)而加大對整個(gè)流場的擾動,從而加速流場的混合。減小流速比可有效改善流場內(nèi)部的混合效果,尤其是等流速時(shí),使撞擊面保持在反應(yīng)器中部,通過撞擊面周期的偏轉(zhuǎn)和渦周期性的生成與脫落,使撞擊面兩側(cè)的流體充分接觸進(jìn)而改善混合。除此之外,反應(yīng)器的結(jié)構(gòu)尺寸對流場的混合存在影響,噴嘴直徑越小,混合效果越好,L/D=3 時(shí)混合效果存在峰值。

    表1 實(shí)驗(yàn)參數(shù)對混合的影響

    總之,在不同的混合尺度下,渦有助于整個(gè)流場的混合,可通過渦空間尺度增大的渦運(yùn)輸作用、渦內(nèi)部流體對等融合的渦融合作用、撞擊面兩側(cè)周期性的渦脫落作用和撞擊區(qū)域大尺度渦受剪切力破裂成小尺度渦的渦破裂作用4個(gè)方面解釋。

    2 撞擊流流場渦特性

    在湍流場中存在大量不同尺度的渦,且不同尺度的渦之間存在非線性作用力,這使得湍流流動狀態(tài)更復(fù)雜[24]。研究人員主要通過模擬和實(shí)驗(yàn)的方法來探究渦結(jié)構(gòu)的復(fù)雜性。渦結(jié)構(gòu)特性常常受反應(yīng)器的結(jié)構(gòu)類型、物料相態(tài)和雷諾數(shù)的影響,因此,渦結(jié)構(gòu)也呈現(xiàn)出不同的特點(diǎn)。

    2.1 撞擊流渦特性實(shí)驗(yàn)

    為客觀真實(shí)地反映出流場內(nèi)的渦特性,研究人員通過激光多普勒技術(shù)(LDV)、平面激光誘導(dǎo)熒光技術(shù)(PLIF)、粒子圖像技術(shù)(PIV)等實(shí)驗(yàn)手段對渦形成的位置、渦系的轉(zhuǎn)變、脫落機(jī)理進(jìn)行研究分析。Zhang 等[25]通過PLIF 研究發(fā)現(xiàn)非定常吞噬流中渦合并的振蕩現(xiàn)象顯著改善混合效果。大量研究人員對渦特性進(jìn)行實(shí)驗(yàn)研究,其研究成果見表2。

    表2 渦特性實(shí)驗(yàn)研究發(fā)展

    三維LDV、PIV技術(shù)通過測量流場的速度信號分析渦的運(yùn)動特性,其優(yōu)勢為:對一點(diǎn)(或區(qū)域)的速度信號持續(xù)測量,方便研究渦體速度變化規(guī)律。PLIF 技術(shù)則是通過濃度信號角度測量流場內(nèi)的渦特性,尤其適合描述渦對混合的影響以及渦體形態(tài)演變過程。通過實(shí)驗(yàn)研究發(fā)現(xiàn)撞擊流反應(yīng)器內(nèi)最典型的渦是噴嘴兩側(cè)產(chǎn)生一對Dean 渦,該渦具有尺度大、對混合效果影響明顯等特點(diǎn)。其原理可以理解為:流體經(jīng)過噴嘴后因發(fā)生撞擊使得流體改變原有的運(yùn)動方向,形成逆壓梯度,撞擊區(qū)域壓力增加、流速降低,流體在離心力作用下轉(zhuǎn)為徑向擴(kuò)散,再與反應(yīng)器壁面撞擊后形成產(chǎn)生二次回流,形成渦對。由于Dean 渦對具有較大的尺度和較高的角速度,其渦核中心的壓力更小,通過卷吸渦邊界流體使其與渦核內(nèi)部流體充分混合,該區(qū)域流場也更趨于無序。同時(shí),流場中渦體之間會發(fā)生碰撞、相互卷吸,甚至渦系還會發(fā)生變化,渦結(jié)構(gòu)特性更加復(fù)雜,也使得流場分析起來較為困難。為此,研究人員利用數(shù)值模擬手段進(jìn)一步對渦特性進(jìn)行分析。

    2.2 撞擊流渦特性數(shù)值模擬

    歐拉法和拉格朗日法是現(xiàn)階段流場研究的主流方法。拉格朗日方法適用于描述渦結(jié)構(gòu)的產(chǎn)生、演化、消亡。在對渦結(jié)構(gòu)模擬分析時(shí),Haller[34]應(yīng)用Lyapunov指數(shù)法,觀測到三維流動含拉格朗日性質(zhì)的渦結(jié)構(gòu)。楊越[35]根據(jù)渦面場方法,得出渦面場的等值面是由渦線所構(gòu)成的渦面。同時(shí),渦面場也適用研究高對稱流與壁流動[36]中的渦結(jié)構(gòu)生成與演化過程。拉格朗日方法可以準(zhǔn)確描述渦面在轉(zhuǎn)捩過程中的幾何變形與渦重聯(lián),及其對阻力生成等重要物理過程[37-38]。由于歐拉法獲得的是瞬時(shí)場,不涉及時(shí)間序列,更適合描述渦的尺度與形態(tài)。因此,在分析渦特性時(shí)歐拉法被科研人員廣泛接受。典型的有基于渦矢量的渦線識別法;基于特征值的Q判據(jù)、Δ 判據(jù)、λ2判據(jù)、λci判據(jù)等識別方法。目前,有多種模擬研究手段對撞擊流反應(yīng)器流場特性進(jìn)行研究,如大渦模擬LES、改進(jìn)直接數(shù)值模擬DSMC、格子Boltzmann 方法,研究人員根據(jù)研究對象的不同來選擇合適的研究方法,具體見表3。

    表3 渦特性模擬研究進(jìn)展

    研究人員通過數(shù)值模擬的方式研究流場內(nèi)的渦特性,得出渦的脫落機(jī)理。如圖4所示[44]為渦脫落的3 個(gè)部分,其中a 部分流體被吸入并不斷增長形成前脫落旋渦,一部分被輸送的流體b切斷了對流的剪切層環(huán)流的供應(yīng),而流體c以逆時(shí)針循環(huán)的方式供給新形成的渦流,且該渦流強(qiáng)度在一段時(shí)間內(nèi)持續(xù)加強(qiáng)。并且用來表征渦脫落強(qiáng)度的量綱為1環(huán)流強(qiáng)度(Γω,max)隨Re的增大而增大。對流場進(jìn)行模擬分析利于研究人員直觀了解渦的結(jié)構(gòu)形式,進(jìn)而方便分析流場與渦的特征。對湍流流場模擬分析時(shí),得到了流場中典型的馬蹄形渦(圖4[40]),其特點(diǎn)是具有較大的尺度、在反應(yīng)器的邊界產(chǎn)生、運(yùn)動方向搖擺不定。

    圖3 馬蹄形渦結(jié)構(gòu)[40]

    圖4 渦脫落機(jī)理[44]

    2.3 不同種類撞擊流反應(yīng)器的渦特性

    近些年,研究人員發(fā)現(xiàn)渦的形成不僅與Re有關(guān),還與反應(yīng)器的種類、尺寸有很大的關(guān)聯(lián)度[45]。隨著研究的深入,發(fā)現(xiàn)Re對渦度的大小和渦的數(shù)量也有影響[46]。通過大量文獻(xiàn)總結(jié)出不同類型撞擊流反應(yīng)器渦結(jié)構(gòu)特性,并對其匯總,見表4。

    從表4中可以看出,無論是T形反應(yīng)器、受限撞擊流反應(yīng)器還是微小型撞擊流反應(yīng)器,渦結(jié)構(gòu)特點(diǎn)有幾處共性:①對于湍流來說,隨著雷諾數(shù)的增加,撞擊區(qū)域變薄,流場的湍動能更大;②針對層流而言,不同結(jié)構(gòu)的撞擊流反應(yīng)器都有渦發(fā)生的臨界值,且臨界值和反應(yīng)器幾何尺寸有關(guān),隨著Re的增加流場上方出現(xiàn)渦,Re繼續(xù)增大,渦的尺度加大,整個(gè)過程混合效果逐漸加強(qiáng)。

    表4 不同種類撞擊流反應(yīng)器渦特性的研究進(jìn)展

    2.4 撞擊流反應(yīng)器多相流渦特性

    在撞擊流研究領(lǐng)域中,除液-液相之外,還有氣-液、氣-固等多相撞擊流反應(yīng)器,研究者首先提出的是以氣體為連續(xù)相的撞擊流反應(yīng)器[56],液-液撞擊流反應(yīng)器在20 世紀(jì)90 年代后才逐步成為研究熱點(diǎn)。郭棟等[57]對氣-液式撞擊流反應(yīng)器的流場特性進(jìn)行數(shù)值分析,研究表明由于進(jìn)口速度的增大使駐點(diǎn)位置產(chǎn)生的渦更多,因此該位置的速度與壓強(qiáng)梯度更大,混合效果也得到顯著改善。Liu 等[58]對氣固兩相流中顆粒旋轉(zhuǎn)進(jìn)行數(shù)值分析。氣體流動如圖5所示,研究顯示,射流的兩側(cè)形成方向相反的氣體旋渦并攜有固體顆粒,這使得顆粒沉積大大減少。撞擊區(qū)內(nèi)的顆粒受相互碰撞和渦的共同作用,顆粒濃度隨氣體入口速度的增加而降低。蔣貴豐等[59]利用本征正交分解(POD)對氣固兩相撞擊流結(jié)構(gòu)進(jìn)行分析。結(jié)果顯示,撞擊流徑向流動形成一對反向旋轉(zhuǎn)的渦結(jié)構(gòu)。加入固體顆粒后,導(dǎo)致軸向產(chǎn)生脈動結(jié)構(gòu),平動能向湍動能轉(zhuǎn)化,大尺度的渦也被破壞成小尺度的湍流結(jié)構(gòu),同時(shí)耗散作用增強(qiáng)。

    在多相撞流反應(yīng)器內(nèi),由于多相物性存在差異,相對密度大的粒子會發(fā)生沉積,撞擊區(qū)域兩側(cè)形成具有較大速度梯度反向旋轉(zhuǎn)的渦,反應(yīng)器內(nèi)的大尺度的、高角速度的渦通過卷積作用使發(fā)生沉積的粒子再次達(dá)到湍動區(qū),強(qiáng)化多相間接觸,即提高粒子的接觸時(shí)間和有效碰撞概率;多相撞擊流反應(yīng)器中存在粒子使原有流場分布規(guī)律發(fā)生改變,軸向射流部分呈脈動狀態(tài),撞擊面兩側(cè)的渦變得更加無序,大尺度的渦破碎也導(dǎo)致規(guī)律性渦系結(jié)構(gòu)消失,同時(shí)在撞擊駐點(diǎn)附近的渦數(shù)量增多,速度、壓力梯度變大,可以強(qiáng)化反應(yīng)器內(nèi)混合效果。

    3 湍流場的POD模態(tài)分析

    在湍流流場中,存在不同尺度相干結(jié)構(gòu)、非相干結(jié)構(gòu)和平均流。通過POD 分析捕捉湍流流場的大尺度相干結(jié)構(gòu)[60],再利用低維動力系統(tǒng),得到其動力學(xué)特性[61]。其優(yōu)點(diǎn)為通過時(shí)間模式替換空間模式,使計(jì)算量大幅減少。在分析結(jié)果中,占主導(dǎo)地位的模態(tài)通常代表流場中較高能量的大尺度結(jié)構(gòu),其特征值可以表示出該階模態(tài)所含能量所占湍動能的比例,也可以是速度場在對應(yīng)特征空間投影的平均能量[62]。根據(jù)Lamotte[63]、Mallor[64]等的研究,在撞擊流反應(yīng)器內(nèi),前4階模態(tài)能量比較突出,并且隨階數(shù)的增加特征值迅速衰減,湍動能降低。實(shí)際流動可視為各界模態(tài)以一定權(quán)重比在時(shí)空域中的組合[65]。Saha 等[66]研究表明隨著模態(tài)階數(shù)的增加,含能比與變化率不斷降低,其能量比如圖6 所示[67],高階模態(tài)對總動能的貢獻(xiàn)很低,對混合的影響效果很小。因此,在對流場進(jìn)行POD 分析時(shí),一般重點(diǎn)研究前4階模態(tài)。當(dāng)前,針對射流領(lǐng)域?qū)u特性的研究主要集中在圓柱射流、平板射流、撞擊射流等形式,POD分析方法在該領(lǐng)域取得了豐碩成果。

    圖6 前若干階本征模態(tài)函數(shù)的能量累積貢獻(xiàn)[67]

    3.1 圓柱射流的POD分析

    圓柱射流與撞擊流有很大的相似之處,射流在經(jīng)過圓柱壁面后發(fā)生的尾跡擾動,與撞擊流的徑向射流區(qū)的流動方式接近。此外,流體在撞擊圓柱后形成的交變尾跡渦使圓柱發(fā)生振動。研究人員通過POD 方法對圓柱射流渦特性進(jìn)行分析,其研究成果詳情見表5。

    表5 圓柱射流的POD研究進(jìn)展

    圓柱繞流流場中的渦與撞擊流場的渦相似,同樣具有周期性,在初始射流失穩(wěn)觸發(fā)的情況下,自由射流渦形成了對稱的二次渦結(jié)構(gòu);不同的是產(chǎn)生的渦周期和渦本身會發(fā)生自身撞擊。對撞擊流反應(yīng)器而言,要優(yōu)化反應(yīng)器的結(jié)構(gòu)充分利用二次渦,提高反應(yīng)器的混合效率。

    3.2 平板射流渦特性的POD分析

    撞擊流最初被Powell[72]認(rèn)為是一種鏡像模型,雖沒有關(guān)注到兩股相向撞擊流的相互作用,但平板射流部分類似撞擊流的鏡像模型,可通過鏡像模型對平板射流渦特性進(jìn)行初步研究。使用POD 方法通過捕捉大尺度渦能量分析渦的運(yùn)動特性、區(qū)分渦系結(jié)構(gòu)、探究脫落頻率等。近年來平板射流的POD分析研究進(jìn)展匯總見表6。

    表6 平板射流的POD研究進(jìn)展

    如圖7 所示,1 階模態(tài)的渦能量明顯高于2 階模態(tài),渦的尺度也更大。對于不同參數(shù)的平板射流而言,渦特征存在一定的共性:一是渦主要發(fā)生在噴嘴出口處;二是射流在撞擊平板時(shí)會形成一對同向渦。差異性主要體現(xiàn)在不同的Re、Ls/Ds、Hs/Ds的渦系結(jié)構(gòu)和空間形態(tài)存在區(qū)別。

    圖7 弦厚比Ls/Ds=1時(shí),不同雷諾數(shù)所對應(yīng)第1、2階模態(tài)[73]

    3.3 撞擊流渦特性的POD分析

    撞擊流的渦特性復(fù)雜多樣,主要受反應(yīng)器的結(jié)構(gòu)型式(如水平對置、雙層撞擊流反應(yīng)器等)、研究工況(如Re、L/D)和物料相態(tài)(氣液、液液、氣固等)的影響。其中,張建偉等[79]首次采用POD方法對撞擊流混合器內(nèi)渦量場進(jìn)行實(shí)驗(yàn)探究,結(jié)果如圖8所示,發(fā)現(xiàn)1階模態(tài)結(jié)構(gòu)宏觀表現(xiàn)為撞擊面兩側(cè)大尺度地反向旋轉(zhuǎn)的一對渦結(jié)構(gòu),并確定了同種動量比下,L/D=3 時(shí)1 階模態(tài)能量占比最高,其結(jié)果如圖9所示。栗晶等[80]對氣固兩相撞擊流反應(yīng)器進(jìn)行POD 分析,發(fā)現(xiàn)顆粒的加入使原有的軸向射流中渦的流向尺度增大。

    圖8 流場速度矢量圖[79]

    圖9 不同L/D下流場能量的變化[79]

    為進(jìn)一步探究反應(yīng)器結(jié)構(gòu)形式對撞擊流反應(yīng)器內(nèi)流場渦結(jié)構(gòu)的影響,張建偉等[10]對雙層撞擊流反應(yīng)器進(jìn)行研究。以上下噴嘴雷諾數(shù)之比K作為變量,在K=1 時(shí)混合效果最好,通過對稱流場發(fā)現(xiàn),L/D=3 時(shí)1 階模態(tài)占比最高,占總能量的23.9%,其混合效果要優(yōu)于不對稱場。對稱撞擊優(yōu)于不對撞擊,原因是對稱射流撞擊后形成的徑向射流速度較大,對下層噴嘴形成的徑向射流速度發(fā)展產(chǎn)生促進(jìn),產(chǎn)生的渦旋較多,渦旋強(qiáng)度增大,提高混合效率。同時(shí)通過POD 方法對速度場進(jìn)行分解,從能量的角度來表征渦的消散,大尺度渦有更高的強(qiáng)度和能量,在撞擊破裂成小尺度渦后強(qiáng)度和湍動能降低,但此時(shí)流場具有較高的脈動特性。隨后,小尺度渦由于流體的黏性作用而發(fā)生耗散。

    通過對不同結(jié)構(gòu)撞擊流裝置進(jìn)行模態(tài)分析可知,改變反應(yīng)器結(jié)構(gòu)和噴嘴間距、進(jìn)口流量都能使1 階模態(tài)能量增大,能量占比增加[81],提高混合效果。

    4 目前研究工作進(jìn)展

    對水平對置撞擊流反應(yīng)器的渦結(jié)構(gòu)進(jìn)行模擬分析,分別得到徑向Q值(圖10)、渦量ω(圖11)。根據(jù)張建偉課題組[79]之前工作可推知,反應(yīng)器的徑向速度對分布徑向Q值、渦量分布有著較大影響,并且趨勢相同呈雙峰分度。Q值在噴嘴邊緣達(dá)到最大,因出口邊緣流體的速度與噴嘴壁面垂直,極易產(chǎn)生渦,因此呈雙峰形式。由于流體噴嘴出口壁面處的剪切力大,進(jìn)而產(chǎn)生渦強(qiáng)度較大,因此呈雙峰分布,受出口影響導(dǎo)致左峰高于右峰。Y在±30mm附近時(shí),兩股流體在撞擊面撞擊后發(fā)生合流,再次出現(xiàn)較大強(qiáng)度的渦。圖12 為不同高度界面過中線沿X方向的分布曲線。從圖中可以看出在撞擊面兩側(cè)出現(xiàn)渦量峰值,且在Y=10mm 時(shí)渦量最大,沿Y軸向上渦量峰值逐漸減小,說明渦體的能量在持續(xù)衰落。但在Y=0mm 處出現(xiàn)大幅度衰落,因?yàn)樵撐恢脼閮晒闪黧w的射流區(qū),兩股流體速度大,削弱了切向力,即產(chǎn)生的渦量也較小。

    圖10 徑向Q值

    圖11 徑向渦量

    圖12 軸向渦量分布曲線

    5 結(jié)語與展望

    從撞擊流反應(yīng)器的混合機(jī)理入手,通過對流場渦特性進(jìn)行分析,總結(jié)出流場內(nèi)渦的生成條件、影響渦量因素、渦的脫落周期等。再結(jié)合三種不同射流形式的POD 模態(tài)分析得到了流場內(nèi)渦的演變規(guī)律以及對混合的影響。未來可以在以下幾個(gè)方向進(jìn)一步對撞擊流反應(yīng)器內(nèi)的渦特性進(jìn)行研究。

    (1)通過優(yōu)化反應(yīng)器結(jié)構(gòu)和改變外部條件來改變反應(yīng)器內(nèi)的微觀混合效果。如采用多層、多向流撞擊和加入超聲激勵(lì)等,都更有利于流場小尺度渦的形成,從而引發(fā)撞擊區(qū)域劇烈湍動,強(qiáng)化流場的傳質(zhì)傳熱的效率。

    (2)在對流場進(jìn)行模擬分析時(shí),可以對網(wǎng)格結(jié)構(gòu)、數(shù)學(xué)模型等進(jìn)行優(yōu)化。反應(yīng)器的網(wǎng)格質(zhì)量關(guān)系到計(jì)算結(jié)果的精確度和準(zhǔn)確性,良好的網(wǎng)格質(zhì)量可以減少計(jì)算機(jī)的計(jì)算量同時(shí)具有較高的精度。在劃分網(wǎng)格時(shí),結(jié)合大數(shù)據(jù)、智能算法等手段優(yōu)化網(wǎng)格,可以提高網(wǎng)格的精度和質(zhì)量,進(jìn)而提高計(jì)算結(jié)果的精確度,保證模擬結(jié)果的可靠性。再結(jié)合PIV、PLIF、VIV等實(shí)驗(yàn)方法來揭示流動混合機(jī)理。

    (3)對于反應(yīng)器內(nèi)的渦結(jié)構(gòu)的形成、長大、消失機(jī)理的理解還不夠深刻,限制了撞擊流技術(shù)的發(fā)展、應(yīng)用和新型裝置的研發(fā)。對此可以嘗試通過非線性分析手段,完善撞擊流的湍流渦特性理論。

    符號說明

    D—— 噴嘴直徑,mm

    D0—— 圓柱直徑,mm

    Ds—— 弓形平板直徑,mm

    Hs—— 弓形平板高度,mm

    h—— 圓柱高度,mm

    IOS —— 離析度

    K—— 上下噴嘴雷諾數(shù)之比

    L—— 噴嘴間距,mm

    L0—— 行程長度,mm

    Ls—— 弓形平板弦長,mm

    Q—— 渦強(qiáng)度,s-2

    Re—— 雷諾數(shù)

    t—— 時(shí)間,s

    Xs—— 離集指數(shù)

    ω—— 渦量,s-1

    Γω,max—— 量綱為1環(huán)流強(qiáng)度

    猜你喜歡
    渦的射流流場
    深海逃逸艙射流注水均壓過程仿真分析
    低壓天然氣泄漏射流擴(kuò)散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計(jì)算
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動研究
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    溝槽對湍流邊界層中展向渦影響的實(shí)驗(yàn)研究
    開縫圓柱縫隙傾斜角對脫落渦的影響
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計(jì)分析
    三角柱體對機(jī)翼根部馬蹄渦的影響
    基于瞬態(tài)流場計(jì)算的滑動軸承靜平衡位置求解
    免费少妇av软件| 国产免费av片在线观看野外av| 久久国产精品人妻蜜桃| 色婷婷久久久亚洲欧美| 日本五十路高清| 免费av毛片视频| 日韩欧美免费精品| 极品人妻少妇av视频| 久久人人爽av亚洲精品天堂| 97人妻天天添夜夜摸| 韩国精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 嫩草影视91久久| 亚洲久久久国产精品| av视频在线观看入口| 91国产中文字幕| 成人三级做爰电影| 操美女的视频在线观看| 在线天堂中文资源库| 午夜老司机福利片| 无人区码免费观看不卡| 母亲3免费完整高清在线观看| 美国免费a级毛片| 亚洲九九香蕉| 日韩大码丰满熟妇| 首页视频小说图片口味搜索| 日本五十路高清| 真人做人爱边吃奶动态| 精品第一国产精品| 亚洲avbb在线观看| 国产欧美日韩一区二区精品| 欧美成狂野欧美在线观看| 亚洲人成电影观看| 黄色片一级片一级黄色片| 嫁个100分男人电影在线观看| 精品欧美国产一区二区三| 国产精品秋霞免费鲁丝片| 午夜亚洲福利在线播放| 国产精品,欧美在线| 欧美日韩黄片免| 91字幕亚洲| 热re99久久国产66热| 人人妻人人澡人人看| 国产精品一区二区在线不卡| av免费在线观看网站| 搡老岳熟女国产| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色 | 精品久久久久久久人妻蜜臀av | 久久久国产精品麻豆| 啪啪无遮挡十八禁网站| 国产av一区在线观看免费| av电影中文网址| 自线自在国产av| 19禁男女啪啪无遮挡网站| 成人手机av| 国产麻豆成人av免费视频| 精品第一国产精品| 国产视频一区二区在线看| 一区二区日韩欧美中文字幕| 亚洲精品久久国产高清桃花| 最近最新中文字幕大全电影3 | 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆 | 亚洲五月天丁香| www国产在线视频色| 99国产精品一区二区三区| 国产一区二区三区视频了| 成人三级黄色视频| 日本a在线网址| 欧美一级毛片孕妇| 欧美乱色亚洲激情| 国产av在哪里看| 精品久久久久久久毛片微露脸| 久久久精品国产亚洲av高清涩受| 人妻丰满熟妇av一区二区三区| 欧美精品亚洲一区二区| 日韩欧美一区二区三区在线观看| 中文字幕另类日韩欧美亚洲嫩草| 老司机午夜十八禁免费视频| 嫁个100分男人电影在线观看| 久久人人97超碰香蕉20202| 亚洲国产欧美日韩在线播放| 老汉色∧v一级毛片| ponron亚洲| 男人舔女人的私密视频| 欧美+亚洲+日韩+国产| 无遮挡黄片免费观看| 日韩欧美一区视频在线观看| 一区福利在线观看| 亚洲熟妇熟女久久| 亚洲狠狠婷婷综合久久图片| 91成年电影在线观看| 99国产精品一区二区蜜桃av| tocl精华| 精品久久蜜臀av无| 91大片在线观看| 国产精品99久久99久久久不卡| 午夜a级毛片| 久久久久久久久久久久大奶| 久久久久精品国产欧美久久久| 国产视频一区二区在线看| 99热只有精品国产| 波多野结衣一区麻豆| 精品一区二区三区av网在线观看| 精品人妻1区二区| 18禁国产床啪视频网站| 午夜免费激情av| 最近最新中文字幕大全电影3 | 一本大道久久a久久精品| 亚洲电影在线观看av| 深夜精品福利| 亚洲欧美一区二区三区黑人| 长腿黑丝高跟| 欧美精品亚洲一区二区| 亚洲国产看品久久| 精品久久久久久,| 国产av一区二区精品久久| 在线免费观看的www视频| АⅤ资源中文在线天堂| 亚洲av成人av| 久久精品影院6| АⅤ资源中文在线天堂| 久久人妻熟女aⅴ| 好男人在线观看高清免费视频 | 日韩一卡2卡3卡4卡2021年| 91大片在线观看| 亚洲最大成人中文| 老汉色∧v一级毛片| 两个人免费观看高清视频| 搡老妇女老女人老熟妇| 多毛熟女@视频| 亚洲国产高清在线一区二区三 | or卡值多少钱| 国产精品永久免费网站| 动漫黄色视频在线观看| 国产av精品麻豆| 精品午夜福利视频在线观看一区| 侵犯人妻中文字幕一二三四区| 1024香蕉在线观看| 国内精品久久久久久久电影| 久久久久国产精品人妻aⅴ院| 亚洲成av片中文字幕在线观看| 一级片免费观看大全| 国产成人啪精品午夜网站| 两性午夜刺激爽爽歪歪视频在线观看 | 我的亚洲天堂| 欧美av亚洲av综合av国产av| а√天堂www在线а√下载| 久久伊人香网站| 变态另类成人亚洲欧美熟女 | 亚洲国产欧美网| 亚洲国产精品sss在线观看| 国产精品爽爽va在线观看网站 | 高清在线国产一区| 嫩草影视91久久| 精品免费久久久久久久清纯| 夜夜躁狠狠躁天天躁| 亚洲 国产 在线| 国产伦人伦偷精品视频| 成人亚洲精品av一区二区| 亚洲欧美激情在线| 变态另类丝袜制服| 久久性视频一级片| 欧美绝顶高潮抽搐喷水| 亚洲片人在线观看| 国产成人av教育| 国产熟女xx| 一a级毛片在线观看| 国产精华一区二区三区| 一本综合久久免费| 国产精品久久久久久精品电影 | 国内久久婷婷六月综合欲色啪| 一边摸一边抽搐一进一出视频| 18禁国产床啪视频网站| 男人操女人黄网站| 欧美一级a爱片免费观看看 | 欧美日韩黄片免| 国产一区二区在线av高清观看| 夜夜夜夜夜久久久久| 一级a爱视频在线免费观看| 亚洲少妇的诱惑av| 男女床上黄色一级片免费看| 久久久国产成人精品二区| 国产精品98久久久久久宅男小说| 亚洲无线在线观看| 久久人妻av系列| 丝袜在线中文字幕| 黄频高清免费视频| 成人欧美大片| 国产主播在线观看一区二区| 97超级碰碰碰精品色视频在线观看| 国产精品野战在线观看| 一进一出抽搐gif免费好疼| 国产欧美日韩一区二区三| 12—13女人毛片做爰片一| 亚洲国产精品999在线| 日本精品一区二区三区蜜桃| 免费人成视频x8x8入口观看| 黄色片一级片一级黄色片| 成人18禁在线播放| 搡老熟女国产l中国老女人| 一边摸一边做爽爽视频免费| 久久欧美精品欧美久久欧美| 婷婷六月久久综合丁香| 性色av乱码一区二区三区2| 麻豆国产av国片精品| 欧美中文综合在线视频| 色综合婷婷激情| 一级黄色大片毛片| 美女大奶头视频| 亚洲av第一区精品v没综合| 中文字幕人妻熟女乱码| 免费在线观看日本一区| 超碰成人久久| 18禁黄网站禁片午夜丰满| 日本三级黄在线观看| 国产1区2区3区精品| 亚洲自偷自拍图片 自拍| 亚洲最大成人中文| 女性被躁到高潮视频| 激情视频va一区二区三区| 免费在线观看影片大全网站| 亚洲熟妇中文字幕五十中出| 久久香蕉国产精品| 波多野结衣一区麻豆| 亚洲欧洲精品一区二区精品久久久| 一本大道久久a久久精品| 亚洲少妇的诱惑av| 叶爱在线成人免费视频播放| 在线观看日韩欧美| 欧美另类亚洲清纯唯美| 91九色精品人成在线观看| 欧美日韩黄片免| 两个人免费观看高清视频| 99国产综合亚洲精品| 可以在线观看的亚洲视频| 欧美最黄视频在线播放免费| 精品国产国语对白av| 美女大奶头视频| 成年人黄色毛片网站| 日韩精品中文字幕看吧| 国产不卡一卡二| 桃红色精品国产亚洲av| 欧美 亚洲 国产 日韩一| 午夜影院日韩av| 青草久久国产| 9热在线视频观看99| 亚洲av日韩精品久久久久久密| 亚洲男人的天堂狠狠| 亚洲全国av大片| 久久精品91无色码中文字幕| 午夜日韩欧美国产| 99精品欧美一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲人成网站在线播放欧美日韩| 高清黄色对白视频在线免费看| 国产99久久九九免费精品| 国产高清激情床上av| 国产精品美女特级片免费视频播放器 | 国产亚洲精品第一综合不卡| 在线观看日韩欧美| 欧美精品亚洲一区二区| 满18在线观看网站| 中文字幕精品免费在线观看视频| 午夜激情av网站| 久久香蕉精品热| 国产蜜桃级精品一区二区三区| 欧美性长视频在线观看| 99精品在免费线老司机午夜| 久久久久久免费高清国产稀缺| 在线观看一区二区三区| 制服诱惑二区| 欧美乱色亚洲激情| 国产私拍福利视频在线观看| 69精品国产乱码久久久| 亚洲国产精品久久男人天堂| 亚洲国产高清在线一区二区三 | 久热爱精品视频在线9| 美女高潮喷水抽搐中文字幕| 99精品久久久久人妻精品| 欧美乱码精品一区二区三区| 啦啦啦免费观看视频1| 亚洲aⅴ乱码一区二区在线播放 | 久久天躁狠狠躁夜夜2o2o| 国产精品永久免费网站| 亚洲成人免费电影在线观看| 久久久久久久午夜电影| 国产亚洲欧美精品永久| 女人爽到高潮嗷嗷叫在线视频| 成人三级黄色视频| 黄片大片在线免费观看| 丝袜美足系列| 波多野结衣高清无吗| 嫩草影院精品99| 极品人妻少妇av视频| 亚洲自偷自拍图片 自拍| 国产91精品成人一区二区三区| 99香蕉大伊视频| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区蜜桃| 色综合站精品国产| 精品一区二区三区四区五区乱码| 脱女人内裤的视频| 久久婷婷成人综合色麻豆| 国产精品香港三级国产av潘金莲| 极品教师在线免费播放| 国产成人欧美在线观看| 久久国产亚洲av麻豆专区| 午夜视频精品福利| 亚洲精品美女久久久久99蜜臀| 国产精品日韩av在线免费观看 | 成年版毛片免费区| av在线播放免费不卡| 91av网站免费观看| 欧美一级毛片孕妇| 老司机深夜福利视频在线观看| 精品国产一区二区三区四区第35| 两个人视频免费观看高清| 黄色丝袜av网址大全| 动漫黄色视频在线观看| videosex国产| 久久精品国产99精品国产亚洲性色 | 欧美最黄视频在线播放免费| 亚洲天堂国产精品一区在线| 国产色视频综合| 一区二区三区国产精品乱码| 欧美在线黄色| 国产亚洲av嫩草精品影院| 成人精品一区二区免费| 久久欧美精品欧美久久欧美| 久久人妻福利社区极品人妻图片| 日韩欧美国产在线观看| 免费看a级黄色片| 18美女黄网站色大片免费观看| 黄片播放在线免费| 久久久久亚洲av毛片大全| 久久人妻福利社区极品人妻图片| 欧美另类亚洲清纯唯美| 午夜两性在线视频| 88av欧美| 视频区欧美日本亚洲| 国产不卡一卡二| 国产区一区二久久| 在线观看免费午夜福利视频| 久久久久久久久免费视频了| 久久婷婷成人综合色麻豆| 男女床上黄色一级片免费看| 好男人在线观看高清免费视频 | 黄色a级毛片大全视频| 国产欧美日韩一区二区三区在线| 国内毛片毛片毛片毛片毛片| 多毛熟女@视频| 久久久久久人人人人人| 男女之事视频高清在线观看| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 久久久久久久午夜电影| 精品久久久久久成人av| 亚洲精品国产区一区二| 国产精品精品国产色婷婷| 日韩欧美在线二视频| 变态另类成人亚洲欧美熟女 | 久久伊人香网站| 嫩草影视91久久| 久久人妻熟女aⅴ| 欧美久久黑人一区二区| 成人特级黄色片久久久久久久| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 一级a爱视频在线免费观看| 久久久久久国产a免费观看| 国产97色在线日韩免费| 亚洲一区中文字幕在线| 亚洲av成人不卡在线观看播放网| 亚洲av美国av| 精品久久蜜臀av无| 亚洲成人久久性| 久久久久久久精品吃奶| 国产在线精品亚洲第一网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成av人片免费观看| 91成年电影在线观看| 久久人妻福利社区极品人妻图片| 一边摸一边抽搐一进一小说| 制服丝袜大香蕉在线| 久久精品国产亚洲av高清一级| 少妇的丰满在线观看| 97人妻天天添夜夜摸| 欧美激情极品国产一区二区三区| 成人18禁高潮啪啪吃奶动态图| 色综合婷婷激情| 涩涩av久久男人的天堂| 亚洲片人在线观看| 天天添夜夜摸| 久久久精品国产亚洲av高清涩受| 亚洲av片天天在线观看| 午夜福利高清视频| 免费人成视频x8x8入口观看| 十八禁人妻一区二区| 午夜福利免费观看在线| 久久 成人 亚洲| 在线天堂中文资源库| 亚洲情色 制服丝袜| 在线永久观看黄色视频| 国内久久婷婷六月综合欲色啪| 亚洲精品久久成人aⅴ小说| 色综合欧美亚洲国产小说| 此物有八面人人有两片| 国产麻豆成人av免费视频| 天天躁夜夜躁狠狠躁躁| 亚洲人成伊人成综合网2020| 亚洲国产精品sss在线观看| 久久久久久亚洲精品国产蜜桃av| 99re在线观看精品视频| 午夜久久久在线观看| 国产精品久久久人人做人人爽| 午夜福利成人在线免费观看| 国产激情欧美一区二区| 午夜福利,免费看| 欧美日韩亚洲国产一区二区在线观看| 亚洲第一av免费看| 熟女少妇亚洲综合色aaa.| 欧美日本亚洲视频在线播放| 母亲3免费完整高清在线观看| 波多野结衣巨乳人妻| 亚洲情色 制服丝袜| 亚洲人成电影免费在线| bbb黄色大片| 日韩国内少妇激情av| 一级a爱片免费观看的视频| 美女高潮到喷水免费观看| 欧美色视频一区免费| 可以在线观看毛片的网站| 国产高清激情床上av| 99久久国产精品久久久| 亚洲 欧美 日韩 在线 免费| 男女之事视频高清在线观看| 亚洲欧美日韩高清在线视频| avwww免费| 后天国语完整版免费观看| 午夜激情av网站| 日韩免费av在线播放| 欧美亚洲日本最大视频资源| 禁无遮挡网站| 欧美成狂野欧美在线观看| 老司机在亚洲福利影院| 亚洲av第一区精品v没综合| 精品高清国产在线一区| 精品一品国产午夜福利视频| 久久久久国产一级毛片高清牌| 无限看片的www在线观看| 女性生殖器流出的白浆| 50天的宝宝边吃奶边哭怎么回事| 黄色女人牲交| 亚洲免费av在线视频| 国产视频一区二区在线看| cao死你这个sao货| 757午夜福利合集在线观看| 国产一区二区三区在线臀色熟女| 精品久久久久久久久久免费视频| 国产91精品成人一区二区三区| 久久久久精品国产欧美久久久| av视频免费观看在线观看| 久久久久久亚洲精品国产蜜桃av| 免费看a级黄色片| 别揉我奶头~嗯~啊~动态视频| 亚洲av五月六月丁香网| 最近最新中文字幕大全免费视频| 国产熟女xx| 成人免费观看视频高清| 精品国产乱子伦一区二区三区| 国产精品免费一区二区三区在线| 欧美日韩亚洲国产一区二区在线观看| 人人妻人人澡欧美一区二区 | 丁香六月欧美| 国产xxxxx性猛交| 亚洲avbb在线观看| 国产又爽黄色视频| 亚洲欧美激情综合另类| 亚洲精品一区av在线观看| av片东京热男人的天堂| 亚洲自拍偷在线| 三级毛片av免费| 每晚都被弄得嗷嗷叫到高潮| 色精品久久人妻99蜜桃| 天天一区二区日本电影三级 | 久久国产精品影院| 99国产精品一区二区蜜桃av| 国产精品国产高清国产av| 久久亚洲真实| 熟妇人妻久久中文字幕3abv| 成在线人永久免费视频| 中文字幕高清在线视频| 香蕉丝袜av| 亚洲av五月六月丁香网| 母亲3免费完整高清在线观看| 亚洲人成网站在线播放欧美日韩| 一本久久中文字幕| 麻豆成人av在线观看| 啪啪无遮挡十八禁网站| 激情视频va一区二区三区| 免费搜索国产男女视频| 亚洲国产精品合色在线| 中亚洲国语对白在线视频| 亚洲精品一卡2卡三卡4卡5卡| 最近最新免费中文字幕在线| 97人妻天天添夜夜摸| 高清毛片免费观看视频网站| 悠悠久久av| 99在线视频只有这里精品首页| 老汉色∧v一级毛片| 99香蕉大伊视频| 好男人电影高清在线观看| 夜夜夜夜夜久久久久| 国产真人三级小视频在线观看| 国产精品亚洲一级av第二区| 黄色女人牲交| 中文字幕最新亚洲高清| www.999成人在线观看| 色综合婷婷激情| 满18在线观看网站| 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| 女人被躁到高潮嗷嗷叫费观| 成年人黄色毛片网站| 国产精品亚洲一级av第二区| 亚洲人成电影观看| 一边摸一边抽搐一进一出视频| 国产精品久久久久久人妻精品电影| av电影中文网址| 国产精品自产拍在线观看55亚洲| 真人做人爱边吃奶动态| 国产一区二区激情短视频| 99国产精品99久久久久| 91av网站免费观看| 免费在线观看视频国产中文字幕亚洲| 国产亚洲欧美在线一区二区| 麻豆一二三区av精品| 长腿黑丝高跟| 久久 成人 亚洲| 亚洲无线在线观看| 亚洲av片天天在线观看| 亚洲无线在线观看| 国产熟女xx| 黄色视频,在线免费观看| 黑人操中国人逼视频| 亚洲九九香蕉| 久久精品国产综合久久久| 看免费av毛片| 色综合亚洲欧美另类图片| 亚洲七黄色美女视频| 99久久99久久久精品蜜桃| 天堂动漫精品| 欧美 亚洲 国产 日韩一| 日韩高清综合在线| 涩涩av久久男人的天堂| 高清在线国产一区| 一个人免费在线观看的高清视频| 中文字幕精品免费在线观看视频| 欧美成狂野欧美在线观看| 波多野结衣高清无吗| 日本在线视频免费播放| 午夜福利免费观看在线| 九色亚洲精品在线播放| 国产国语露脸激情在线看| 一a级毛片在线观看| 亚洲 欧美一区二区三区| 69av精品久久久久久| 男人舔女人的私密视频| 欧美中文综合在线视频| 久久久久久大精品| 国产精品亚洲一级av第二区| 桃色一区二区三区在线观看| 亚洲男人天堂网一区| 久久精品亚洲熟妇少妇任你| 久久久久九九精品影院| 天天躁狠狠躁夜夜躁狠狠躁| 精品国产乱子伦一区二区三区| 少妇裸体淫交视频免费看高清 | 男女下面插进去视频免费观看| 麻豆久久精品国产亚洲av| 午夜福利成人在线免费观看| 国产欧美日韩精品亚洲av| 欧美激情久久久久久爽电影 | 一级a爱片免费观看的视频| 一级毛片女人18水好多| 国产成人精品久久二区二区91| 亚洲专区中文字幕在线| 国产高清videossex| 俄罗斯特黄特色一大片| 黄色 视频免费看| 法律面前人人平等表现在哪些方面| 久久久国产精品麻豆| 亚洲成人免费电影在线观看| 精品久久蜜臀av无| 女性生殖器流出的白浆| 亚洲五月色婷婷综合| 久久久久久亚洲精品国产蜜桃av| 亚洲男人天堂网一区| 老司机在亚洲福利影院| 99久久精品国产亚洲精品| 欧美黑人欧美精品刺激| 老熟妇乱子伦视频在线观看| 国产国语露脸激情在线看| 91麻豆精品激情在线观看国产| 欧美 亚洲 国产 日韩一| 国产精品美女特级片免费视频播放器 | 久久精品国产亚洲av香蕉五月| 精品高清国产在线一区| 两性夫妻黄色片| 欧美在线黄色| 精品高清国产在线一区| 制服诱惑二区| 国产高清激情床上av|