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

    渦輪槳攪拌槽內(nèi)湍流特性的V3V實驗及大渦模擬

    2017-11-22 10:20:05施乃進(jìn)周勇軍鮑蘇洋辛偉陶蘭蘭
    化工學(xué)報 2017年11期
    關(guān)鍵詞:尾渦大渦渦量

    施乃進(jìn),周勇軍,鮑蘇洋,辛偉,陶蘭蘭

    (1南京工業(yè)大學(xué)機械與動力工程學(xué)院,江蘇 南京 211816;2江蘇省特種設(shè)備安全監(jiān)督檢驗研究院國家化工設(shè)備質(zhì)量監(jiān)督檢驗中心,江蘇 蘇州 215600)

    渦輪槳攪拌槽內(nèi)湍流特性的V3V實驗及大渦模擬

    施乃進(jìn)1,周勇軍1,鮑蘇洋1,辛偉2,陶蘭蘭2

    (1南京工業(yè)大學(xué)機械與動力工程學(xué)院,江蘇 南京 211816;2江蘇省特種設(shè)備安全監(jiān)督檢驗研究院國家化工設(shè)備質(zhì)量監(jiān)督檢驗中心,江蘇 蘇州 215600)

    分別用體三維速度測量技術(shù)(volumetric three-component velocimetry measurements,V3V)和大渦模擬(large eddy simulation,LES)方法對渦輪槳攪拌槽內(nèi)流場進(jìn)行研究,發(fā)現(xiàn)在完全湍流狀態(tài)下,渦輪槳攪拌槽內(nèi)流場的量綱1相平均速度及湍動能分布同Reynolds數(shù)無關(guān)。用V3V方法實現(xiàn)了Rushton槳葉附近三維流場的重構(gòu);探討尾渦的三維結(jié)構(gòu)及運動規(guī)律;分析了葉片后方 30°截面軸向、徑向和環(huán)向速度沿徑向分布規(guī)律。用V3V實驗結(jié)果對比了2D-PIV(particle image velocimetry)數(shù)據(jù)中的尾渦渦對位置和渦量,渦對位置吻合度較好,但2D-PIV中渦量較V3V小37.5%;通過大渦模擬得到完整的尾渦結(jié)構(gòu),發(fā)現(xiàn)在葉片上邊緣后側(cè)存在一個和尾渦形成方式相同但不成對出現(xiàn)的渦結(jié)構(gòu);將大渦模擬結(jié)果和2D-PIV及V3V實驗結(jié)果對比發(fā)現(xiàn),大渦模擬在速度分布及尾渦運動軌跡方面均同實驗結(jié)果吻合較好,表明大渦模擬能較好地預(yù)測渦輪槳攪拌槽內(nèi)流場。

    V3V;PIV;大渦模擬;尾渦;渦輪槳

    引 言

    渦輪槳攪拌槽廣泛應(yīng)用于化工、生物、石化、食品和制藥等過程工業(yè)[1-2]。攪拌槽內(nèi)流場結(jié)構(gòu)十分復(fù)雜。攪拌槳通過槳葉附近流體將機械能逐漸傳遞至攪拌槽內(nèi),葉片產(chǎn)生的尾渦是這一區(qū)域流場的顯著特征,攪拌槳葉周圍的尾渦有利于流體進(jìn)行充分的混合[3-4]。因此,有效地認(rèn)識渦輪槳攪拌槽的流場特性及尾渦性質(zhì),對于攪拌槽的設(shè)計及優(yōu)化有著重大意義。

    目前,國內(nèi)外對渦輪槳攪拌槽內(nèi)流場的實驗研究較多,主要集中在2D-PIV和三維PIV研究。程先明等[5]用 2D-PIV方法實驗研究了相同功率下不同葉片長度Rushton槳的流場,研究了槳葉離底距離對尾渦特性的影響;Chara等[6]使用2D-PIV技術(shù),研究了Reynolds數(shù)對Rushton槳攪拌槽內(nèi)尾渦分布的影響;Sharp等[7]用PIV方法研究發(fā)現(xiàn)Rushton攪拌槳在尾渦結(jié)構(gòu)處嚴(yán)重趨于各向異性;Escudie等[8-9]采用三維PIV方法研究了渦輪攪拌槳的尾渦特性,發(fā)現(xiàn)三維 PIV的渦尺寸較大,但尾渦發(fā)展軌跡與2D-PIV結(jié)果相似;宋戈[10]用2D-PIV和三維PIV結(jié)合的相位解析測量方法研究了渦輪槳攪拌槽內(nèi)流場的湍流特性,發(fā)現(xiàn)標(biāo)準(zhǔn)Rushton槳射流區(qū)和尾渦結(jié)構(gòu)區(qū)流場偏離各向同性,且Reynolds數(shù)對各向異性的影響不大。

    由于V3V方法較二維PIV方法先進(jìn),能測量到攪拌槽內(nèi)流場的環(huán)向速度,彌補了長期以來對環(huán)向速度實驗測量的空白。Troolin等[11-12]將V3V實驗方法首先應(yīng)用于渦輪槳攪拌槽內(nèi)流場研究,并取得一定成果。他們分析了渦旋軌跡并且計算了渦流強度,對渦旋的3種識別方法進(jìn)行了比較,驗證了V3V方法的可行性;分析了葉端旋渦的結(jié)構(gòu)和徑向射流的高斯特性,揭示了葉輪附近高剪切區(qū)的湍流產(chǎn)生原理。鮑蘇洋等[13]采用V3V方法對槳葉附近湍流各向同性假設(shè)進(jìn)行了分析,發(fā)現(xiàn)槳葉區(qū)和尾渦所在位置湍動能被各向同性假設(shè)近似法高估了 25%~33%,槳葉區(qū)和尾渦所在位置趨向于各向異性。

    近年來有研究者開始使用大渦模擬方法研究攪拌槽內(nèi)三維流場的湍流性能。張艷紅等[14]、苗一等[15]、李志鵬等[16]采用LES方法研究了渦輪槳攪拌槽內(nèi)單循環(huán)流動特性,得到了攪拌槽內(nèi)流場的瞬時流場,并初步探討了尾渦的發(fā)展規(guī)律,與實驗結(jié)果對比,發(fā)現(xiàn)LES方法可以準(zhǔn)確預(yù)測攪拌槽內(nèi)湍流流動的非穩(wěn)態(tài)及周期性脈動;Bao等[17]用2D-PIV和LES方法對比研究了渦輪槳攪拌槽內(nèi)流場,結(jié)果表明湍動能受尾渦位置影響;Li等[18]分別用三維PIV方法和LES方法,對Rushton攪拌槳等3種不同槳型下攪拌槽內(nèi)的流場特性進(jìn)行了對比研究,發(fā)現(xiàn)尾渦在槳葉上方形成橢圓形尾渦,尾渦所在位置為高湍動能區(qū)域。

    2D-PIV方法測量區(qū)域較大,但無法實現(xiàn)對環(huán)向速度的測量,只能根據(jù)各向同性假設(shè)預(yù)測流場的環(huán)向速度。V3V方法可以對140 mm×140 mm×100 mm的空間區(qū)域內(nèi)三維流場進(jìn)行測量,實現(xiàn)對三維流場的重構(gòu),有利于理解尾渦促進(jìn)流體混合過程。雖然有涉及渦輪槳攪拌槽內(nèi)尾渦鑒別方法的 V3V及 2D-PIV實驗研究,但未實現(xiàn)渦輪槳攪拌槽內(nèi)尾渦的三維結(jié)構(gòu)圖的實驗重構(gòu)及大渦模擬,未涉及流場的V3V實驗研究及LES方法對比研究。

    本文采用V3V技術(shù),分析了渦輪槳攪拌槽內(nèi)槳葉附近流場的流動特性,深入研究尾渦的結(jié)構(gòu)和發(fā)展規(guī)律,并驗證 LES方法在流場預(yù)測方面與 V3V及 2D-PIV實驗方法的差距,為模擬方法的改進(jìn)提供依據(jù),以期為工業(yè)放大提供參考。

    1 實驗裝置和方法

    本實驗在內(nèi)徑T=420 mm的圓柱形透明攪拌槽內(nèi)進(jìn)行,攪拌槽主體由有機玻璃制成,底部為不銹鋼標(biāo)準(zhǔn)橢圓封頭,配置單層標(biāo)準(zhǔn)六直葉渦輪槳,槳葉直徑 D=0.48T,離底距離 C=T/2,無擋板條件,實驗體系為單相水體系,液面高度H=550 mm,如圖1所示。槳葉轉(zhuǎn)速N=1.67 s-1,Re=66007。為了減弱攪拌槽曲面光束折射對測量產(chǎn)生的影響,圓形攪拌槽外設(shè)置有同樣工作介質(zhì)的方形玻璃槽,液面與攪拌槽內(nèi)水面持平。

    采用美國TSI公司V3V測試系統(tǒng)測量攪拌槽內(nèi)流場,脈沖激光器的能量為 380 mJ。V3V能夠?qū)w區(qū)域完成“瞬態(tài)”測量和“長時間”的測量,本實驗用V3V對長140 mm、寬140 mm、深100 mm的區(qū)域進(jìn)行測量,如圖2所示。3D相機探頭是由3個400像素的CCD相機組成,3個相機以“共面-三角”形式集成在同一個面板上,相機的重合拍攝區(qū)域為V3V的測量區(qū)域。測量區(qū)域中每個粒子的圖像同時被3個相機捕捉,每個粒子在3個相機成像并構(gòu)成三角形。對三角形的分析得到粒子的空間位置。利用針孔成像的原理分析三角形,從而得到了粒子的空間位置的速度矢量。示蹤粒子為空心玻璃珠,平均直徑約為50 μm,密度與水相同。示蹤粒子在激光器拍攝時,選擇要求避光性較好的實驗室,避免外界光線對粒子的影響。

    圖1 攪拌槽結(jié)構(gòu)Fig.1 Scheme of stirred tank

    2 大渦模擬

    2.1 大渦模型

    LES模擬是介于Reynolds平均法和直接模擬法之間的一種湍流數(shù)值模擬方法,對非穩(wěn)態(tài)流場具有較好的預(yù)測效果。LES模型利用濾波函數(shù)將流場中的渦劃分為大、小尺度的渦,對大尺度的渦直接計算N-S方程,用亞格子尺度應(yīng)力處理小尺度的渦[19]。大渦模擬方法由于采用非穩(wěn)態(tài)的計算方法,對計算機運算能力的要求相對較高,計算時間較長。在大渦模擬中,對不可壓縮流體,通過連續(xù)性方程和動量方程進(jìn)行空間濾波后,可得控制方程如下[20]

    其中,τij由亞格子應(yīng)力模型求取,一般采用標(biāo)準(zhǔn)Smagorinsky-Lily模型,但引入的模型參數(shù)CS會對湍流參數(shù)(湍動能耗散率)產(chǎn)生影響。本工作采用改進(jìn)的動態(tài)Smagorinsky-Lily模型[21],該模型認(rèn)為子應(yīng)力與局部流體應(yīng)變率呈正比,表達(dá)式如下

    其中,κ為卡門常數(shù),一般取 0.42;在動態(tài)Smagorinsky-Lily模型中CS是變化的,通過局部流場信息計算,從而更真實地反映流場情況。

    2.2 數(shù)值求解方法

    為了便于實驗測試結(jié)果與模擬結(jié)果相比較,本文的模擬模型與實驗?zāi)P拖嘁恢?。渦輪槳攪拌槽結(jié)構(gòu)較復(fù)雜,采用非結(jié)構(gòu)化四面體網(wǎng)格對模型進(jìn)行離散化,并對槳葉和軸附近網(wǎng)格進(jìn)行加密處理,其中槳葉區(qū)網(wǎng)格尺寸為1~2 mm,其余區(qū)域網(wǎng)格尺寸為4~6 mm,網(wǎng)格結(jié)構(gòu)如圖3所示。槳葉的旋轉(zhuǎn)運動通過多重參考系(MRF)方法實現(xiàn)。大渦模擬中動量方程的離散采用中心差分格式,時間推進(jìn)采用二階精度的隱式格式,壓力與速度的耦合采用 PISO算法。標(biāo)準(zhǔn) k-?模型計算時分別采用二階迎風(fēng)格式和SIMPLE算法。計算時先利用標(biāo)準(zhǔn)k-?模型計算初始流場,然后再進(jìn)行非穩(wěn)態(tài)的大渦模擬[22]。每個槳葉旋轉(zhuǎn)周期內(nèi)取36個時間步,每個時間步內(nèi)速度殘差均收斂至10-4。約20個槳葉旋轉(zhuǎn)周期后,流場達(dá)到近似穩(wěn)態(tài),然后進(jìn)行樣本數(shù)據(jù)的收集、統(tǒng)計和后處理等工作。

    圖2 V3V的測量區(qū)域Fig.2 Measurement area of V3V

    圖3 計算域網(wǎng)格示意圖Fig.3 Grid of computational regional

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

    本文通過攪拌功率準(zhǔn)數(shù)和網(wǎng)格數(shù)相關(guān)關(guān)系判定網(wǎng)格無關(guān)性[23]。表1是同一個模型不同網(wǎng)格數(shù)情況下功率準(zhǔn)數(shù)及其相對偏差。從表中可發(fā)現(xiàn),當(dāng)網(wǎng)格數(shù)達(dá)到95萬時,功率準(zhǔn)數(shù)的偏差已經(jīng)小于3%,再增加網(wǎng)格數(shù)到105萬對模擬結(jié)果的影響不大,綜合考慮本研究選擇網(wǎng)格數(shù)為 95萬的網(wǎng)格尺寸模型對流場進(jìn)行大渦模擬。

    表1 網(wǎng)格數(shù)與功率準(zhǔn)數(shù)及其相對偏差Table 1 Effect of grid size on power number and deviation

    3 實驗結(jié)果與討論

    3.1 Reynolds數(shù)對流場特性的影響

    本實驗通過調(diào)整轉(zhuǎn)速改變 Reynolds數(shù)的方法來實現(xiàn)對不同Reynolds數(shù)情況下流場的PIV結(jié)果分析,攪拌槽內(nèi)流體不同轉(zhuǎn)速下的Reynolds數(shù)如表2所示,選取槳葉后θ=30°,2z/W=0位置上流體速度作為研究對象。圖4為不同Reynolds數(shù)V3V量綱1相位平均速度和湍動能分布,發(fā)現(xiàn)標(biāo)準(zhǔn)化徑向速度、軸向速度、環(huán)向速度和湍動能分布基本不隨Reynolds數(shù)的增大而改變,這同Zhou等[24]、Yoon等[25]和Li等[16,18]的完全湍流狀態(tài)下(Re>40000),相同攪拌槽內(nèi)流場的量綱1相位平均速度及湍動能分布同 Reynolds數(shù)無關(guān)結(jié)論十分吻合,而同 Yoon等[26]的結(jié)論不一致,這主要是由于他們的實驗未達(dá)到完全湍流狀態(tài)。湍動能只有在尾渦位置(r/R=1.2和1.6)處變化較大,而其余位置吻合非常好,這主要是由于尾渦中心處傾向于各向異性。因此,在完全湍流狀態(tài)下,攪拌槽內(nèi)量綱1相位平均速度及湍動能分布同Reynolds數(shù)無關(guān),這將為攪拌器在工業(yè)放大中的設(shè)計及改進(jìn)提供重要參考。

    表2 不同轉(zhuǎn)速下的Reynolds數(shù)Table 2 Reynolds number under different rotation speeds

    3.2 整體流型分析

    3.2.1 三維流場V3V重構(gòu) 本實驗在Reynolds數(shù)為66007的情況下,對長140 mm、寬140 mm、深100 mm的三維流場區(qū)域進(jìn)行V3V測量,配對得約18000個矢量。圖5為速度等值面和含速度矢量的水平切面圖。由圖可知,在x*=0.4~0.8的渦輪槳葉片區(qū)域,速度分布十分均勻。只有x*=1.2~1.4的射流區(qū)速度變化顯著,其中0.8Vtip等值面非常小,說明高速射流區(qū)非常小,衰減過程十分明顯。圖6為垂直切片及3個速度等值面主視圖,從圖中可以清楚地看出最大速度分布在x*=1.2~1.3,y*=0.5~1區(qū)域內(nèi),而不是槳葉中心高度 2z/W=0處呈對稱分布,這表明射流區(qū)流體有明顯的向上軸向速度,這是由于本實驗得到的為標(biāo)準(zhǔn)的雙循環(huán)流型,流體從葉片端部射流排出后稍微偏上運動。其中0.8Vtip等值面(紅色)被0.67Vtip等值面(黃色)包絡(luò)其中,說明此區(qū)域即為前槳葉的高速射流區(qū)。槳葉寬度位置附近0.42Vtip(藍(lán)色)等值面沿y*=0.5基本對稱,并且較其余位置遠(yuǎn)離槳葉區(qū),這是由于槳葉寬度流體主要受槳葉射流影響,徑向速度較大;而其余位置流體是受射流區(qū)流體的裹挾及槳葉的卷吸綜合作用的結(jié)果,徑向速度u相對較小,運動較慢。

    圖4 不同Reynolds數(shù)V3V量綱1相位平均速度和湍動能分布(θ=30°,2z/W=0)Fig.4 V3V normalized phase-resolved velocity and TKE profiles for different Re(θ=30°,2z/W=0)

    圖5 速度等值面和水平切面圖(Re=66007)Fig.5 Isosurface of normalized velocity and three horizontal slices(Re=66007)

    3.2.2 流場速度分布 圖7為槳葉附近區(qū)域不同軸向位置 V3V 速度分布(Re=66007,θ=30°)。從圖可看出,在r/R=0.6~0.9區(qū)域,徑向速度、軸向速度和環(huán)向速度大小與軸向位置無關(guān),因為此區(qū)域?qū)儆跇~葉片覆蓋區(qū)域,流場主要受葉片影響。而在r/R>0.9位置,流體受射流區(qū)及邊界條件影響增大,尤其是軸向速度,在r/R=0.9位置處即受到雙循環(huán)流型影響,在2z/W=-1高度出現(xiàn)軸向向上速度,而2z/W=1高度卻出現(xiàn)軸向向下速度。通過對比3幅圖可以得出,渦輪槳攪拌槽內(nèi)流場徑向射流較強,射流區(qū)域中心位置和槳葉卷吸位置處速度突變明顯,越遠(yuǎn)離該區(qū)域的流體,速度變化越平緩。

    圖6 垂直切片及3個速度等值面主視圖Fig.6 Front view of a vertical slice and three isosurfaces of normalized velocity

    圖8為θ=30°時V3V及2D-PIV速度場云圖對比。由圖可看出V3V及2D-PIV速度場云圖吻合非常好,圖左下角均出現(xiàn)速度較大區(qū)域;槳葉前端射流有明顯向上速度,且最大速度區(qū)域出現(xiàn)在r/R=1.2~1.3,2z/W=0.5~1區(qū)域,此位置亦為尾渦所在位置。槳葉右上方流體速度較大,且該位置速度矢量圖可以判斷出,在右上方為一個大型循環(huán)流,這是由槳葉形狀和離底距離決定的。由于V3V測量的是流體的實際速度(u、v、w),而 2D-PIV只測量得到u、v,缺少環(huán)向速度w,且由于葉片尖端射流區(qū)環(huán)向速度較大,因此在r/R=1.3~1.4,2z/W=1位置附近云圖同 2D-PIV出現(xiàn)差別,這主要由2D-PIV測量方法本身的缺陷引起的。

    圖9為LES與實驗結(jié)果的標(biāo)準(zhǔn)化徑向、軸向速度對比(2z/W=0,θ=30°),從圖可以看出,大渦模擬方法與實驗結(jié)果吻合較好,即大渦模擬分別在r/R=1.2和r/R=1.7位置出現(xiàn)徑向速度極大值,在r/R=1.4附近位置出現(xiàn)軸向速度極大值,且在極大值附近速度吻合非常好,只是在遠(yuǎn)離槳葉位置處LES結(jié)果較實驗結(jié)果更平滑,而實驗結(jié)果中可以發(fā)現(xiàn)明顯的隨機脈動引起的流體速度波動,且遠(yuǎn)離槳葉位置處波動對結(jié)果影響較大,這和尾渦結(jié)構(gòu)在遠(yuǎn)離槳葉位置出現(xiàn)的結(jié)構(gòu)不連續(xù)性一致。

    圖7 不同軸向位置V3V速度分布圖(Re=66007,θ=30°)Fig.7 V3V velocity profiles for different axial location(Re=66007,θ=30°)

    3.3 尾渦分析

    圖8 V3V及2D-PIV速度場云圖對比Fig.8 Comparison of velocity profile from V3V and 2D-PIV

    圖9 標(biāo)準(zhǔn)化徑向、軸向速度對比(2z/W=0, θ=30°)Fig.9 Comparision of normalized radial and axial velocity(2z/W=0, θ=30°)

    圖10 渦量等值面和速度等值面圖(Re=66007)Fig.10 Isosurfaces of normalized vorticity and velocity(Re=66007)

    3.3.1 尾渦結(jié)構(gòu)分析 圖10為轉(zhuǎn)速N=1.67 s-1時,V3V測量的渦量等值面和速度等值面圖。由圖可以看出,尾渦渦對中間為高速射流區(qū),即尾渦渦對分布于高速射流區(qū)兩側(cè),且尾渦運動稍滯后于高速射流區(qū);相同渦量等值面,上尾渦明顯較下尾渦大,下尾渦徑向運動較上尾渦快,尾渦耗散亦較快;隨著θ的增大,尾渦不斷耗散,因此尾渦渦對結(jié)構(gòu)連續(xù)性變差,由于下尾渦運動及耗散較快,因此在尾渦即將耗散完畢處結(jié)構(gòu)連續(xù)性更差。徑向射流遇筒體壁面分流,分別沿上、下兩方向運動,在槳葉上部和下部分別形成了兩個大型的橢圓形渦旋,尾渦右上方和右下方兩渦量較大區(qū)域正是兩渦旋中心位置。

    圖11為LES尾渦結(jié)構(gòu)圖,圖中可明顯看出尾渦渦對成型、發(fā)展及耗散消失過程,上尾渦較下尾渦耗散快,下尾渦在θ=60°附近消失,上尾渦在90°才耗散完畢,且θ=0~40°范圍內(nèi),尾渦結(jié)構(gòu)、位移較一致,在θ>40°時,受湍流隨機脈動影響很大,甚至部分尾渦結(jié)構(gòu)出現(xiàn)不連續(xù)現(xiàn)象,這和圖10的V3V測量結(jié)果中尾渦結(jié)構(gòu)不連續(xù)性非常吻合。可以明顯看出下尾渦較上尾渦徑向運動快,且耗散亦較快。該結(jié)論同Li等[18]的SPIV方法得出的公認(rèn)尾渦結(jié)構(gòu)形成及發(fā)展規(guī)律基本一致,只是在尾渦耗散完畢的位置處尚存在爭議。

    圖11 大渦模擬尾渦結(jié)構(gòu)Fig.11 Visualization of trailing vorticies obtained by LES

    除了在每個葉片后面存在一個渦對外,在葉片上方或者后方還存在一個較大的渦結(jié)構(gòu),但結(jié)構(gòu)并不完整,大小不同,規(guī)律性較尾渦渦對弱,且現(xiàn)象并非像尾渦一樣成對出現(xiàn),而只是在圓盤上面區(qū)域出現(xiàn)。結(jié)合帶速度矢量的渦量云圖可以發(fā)現(xiàn),該渦成型過程和尾渦相似,只是此渦為槳葉內(nèi)側(cè)上邊緣附近流體受槳葉葉片厚度方向壓力梯度影響,流體繞著葉片上邊緣旋轉(zhuǎn),形成了一個渦結(jié)構(gòu),按照同尾渦相同的鑒別方法,該渦在θ=40°附近即耗散完畢,這表明此尾渦的結(jié)構(gòu)保持性很差;且隨著角度增大,渦沿軸向向下、徑向向外運動,渦結(jié)構(gòu)范圍r/R=0.5~0.8,軸向位置不低于圓盤。

    3.3.2 尾渦特性的對比分析 圖12為θ=30°時截面渦量分布。可以看出,V3V和2D-PIV實驗結(jié)果均出現(xiàn)兩個渦量較大區(qū)域,此位置即為尾渦渦對位置,隨著距離尾渦距離增大,尾渦渦量呈逐漸減小趨勢。且V3V實驗結(jié)果中尾渦渦對位置同2D-PIV結(jié)果吻合非常好。但是,2D-PIV結(jié)果中渦量較V3V中小37.5%,這是由于V3V可以實現(xiàn)三維速度的測量,而 2D-PIV僅能測量徑向速度及軸向速度,因此渦量計算缺少環(huán)向速度分量。

    圖12 V3V和2D-PIV渦量對比(θ=30°,Re=66007)Fig.12 Comparison of vorticity magnitude from V3V and 2D-PIV (θ=30°,Re=66007)

    目前,確定尾渦中心主要有3種方法,分別為速度法、最大渦量法及特征值法。Escudie等[8]采用這3種方法分別對尾渦中心進(jìn)行研究,發(fā)現(xiàn)3種方法得出的結(jié)果相差不大。本文采用最大渦量法來確定尾渦中心,渦量最大處即為尾渦的中心位置,研究了尾渦的運動軌跡。圖13為V3V和LES尾渦徑向運動軌跡對比。從圖中可以發(fā)現(xiàn)尾渦運動軌跡基本一致,但LES研究結(jié)果表明,由于雙循環(huán)結(jié)構(gòu)中下尾渦運動較上尾渦快,耗散亦較快,因此下尾渦在θ=60°基本耗散完畢,而上尾渦結(jié)構(gòu)一直持續(xù)到槳葉后方90°位置,但Escudie等[8]用DES方法發(fā)現(xiàn)尾渦最大持續(xù)到槳葉后方100°,Sharp等[11]用V3V方法在槳葉后方140°位置發(fā)現(xiàn)尾渦結(jié)構(gòu)。因此,渦輪槳尾渦耗散完畢位置有待繼續(xù)研究。

    圖13 2D-PIV和LES尾渦徑向運動軌跡對比Fig.13 Radial displacement of trailing vortices obtained by 2D-PIV and LES

    4 結(jié) 論

    (1)采用 V3V方法對渦輪槳攪拌槽內(nèi)槳葉附近流場進(jìn)行了實驗研究,在完全湍流狀態(tài)下(Re>40000),渦輪槳攪拌槽內(nèi)流場的量綱1相位平均速度及湍動能分布同 Reynolds數(shù)無關(guān)。對比V3V和2D-PIV的速度場云圖分布,發(fā)現(xiàn)兩種方法研究結(jié)果吻合很好,尤其是在槳葉射流區(qū)。

    (2)用V3V實驗方法實現(xiàn)了尾渦的三維重構(gòu),下尾渦運動比上尾渦稍快,且下渦量耗散亦較快,且渦對沿軸向向上運動;尾渦渦對分布于高速射流區(qū)兩側(cè),且尾渦運動稍滯后于高速射流區(qū);隨著θ的變化,尾渦不斷耗散,因此尾渦渦對結(jié)構(gòu),尤其是下尾渦連續(xù)性變差,對比 2D-PIV實驗結(jié)果中尾渦渦對位置,吻合度非常高。

    (3)對渦輪槳攪拌槽內(nèi)流場進(jìn)行了大渦模擬研究,模擬得到了非常完整的尾渦結(jié)構(gòu),發(fā)現(xiàn)下尾渦在θ=60°附近消失,上尾渦在 90°才耗散完畢,且θ=0~40°范圍內(nèi),尾渦結(jié)構(gòu)、位移較一致,在θ>40°時,部分尾渦結(jié)構(gòu)出現(xiàn)不連續(xù)現(xiàn)象。將大渦模擬結(jié)果和2D-PIV及V3V實驗結(jié)果對比發(fā)現(xiàn),大渦模擬在速度分布及尾渦運動軌跡方面均同實驗結(jié)果吻合較好,表明大渦模擬能較好預(yù)測渦輪槳攪拌槽內(nèi)流場。

    符 號 說 明

    c——離底距離,m

    D——槳葉直徑,m

    H——液面高度,m

    N——槳葉轉(zhuǎn)速,s-1

    R——渦輪槳槳葉半徑,m

    Re——Reynolds數(shù)

    r——坐標(biāo)系徑向坐標(biāo),m

    T——攪拌槽直徑,m

    u——徑向速度,m·s-1

    V——流體速度,m·s-1

    Vtip——葉端線速度,m·s-1

    v——軸向速度,m·s-1

    W——槳葉葉片寬度,m

    w——環(huán)向速度,m·s-1

    x*——x/R

    y*——2y/W

    z——坐標(biāo)系軸向坐標(biāo),m

    θ——槳葉后角度,(°)

    [1] 劉敏珊, 張麗娜, 董其伍. 渦輪槳攪拌槽內(nèi)混合特性模擬研究[J].工程熱物理學(xué)報, 2009, 30(10): 1700-1702.LIU M S, ZHANG L N, DONG Q W. Numerical research of mixing characteristics in a stirred tank with turbine impellers[J]. Journal of Engineering Thermophysics, 2009, 30(10): 1700-1702.

    [2] 周國忠, 王英琛, 施力田. 用 CFD研究攪拌槽內(nèi)的混合過程[J].化工學(xué)報, 2003, 54(7): 886-890.ZHOU G Z, WANG Y C, SHI L T. CFD study of mixing process in stirred tank[J]. Journal of Chemical Industry and Engineering (China),2003, 54(7): 886-890.

    [3] 蔡清白, 戴干策. 翼型槳葉片尾渦結(jié)構(gòu)的PIV實測與LES模擬研究[J]. 過程工程學(xué)報, 2011, (4): 541-548.CAI Q B, DAI G C. Measurement of trailing vortices around a hydrofoil impeller with PIV and large eddy simulation[J]. The Chinese Journal of Process Engineering, 2011, (4): 541-548.

    [4] SHARP K V, ADRIAN R J. PIV study of small-scale flow structure around a Rushton turbine[J]. AIChE Journal, 2001,47(4): 766-778.

    [5] 程先明, 李志鵬, 高正明, 等. 渦輪槳攪拌槽內(nèi)流動及尾渦特性研究[J]. 北京化工大學(xué)學(xué)報(自然科學(xué)版), 2009, 36(6): 16-21.CHENG X M, LI Z P, GAO Z M,et al. Characteristics of flow fields and trailing vortices in a stirred tank with a Rushton turbine[J].Journal of Beijing University of Chemical Technology(Natural Science), 2009, 36(6): 16-21.

    [6] CHARA Z, KYSELA B, KONFRST J,et al. Study of fluid flow in baffled vessels stirred by a Rushton standard impeller[J]. Applied Mathematics and Computation, 2016, (272): 614-628.

    [7] SHARP K V, ADRIAN R J. PIV study of small-scale flow structure around a Rushton turbine[J]. AIChE Journal, 2001, 47(1): 766-778.

    [8] ESCUDIE R, BOUYER D, LINE A. Characterization of trailing vortices generated by a Rushton turbine[J]. AIChE Journal, 2004,50(1): 75-86.

    [9] ESCUDIE R, LINE A. A simplified procedure to identify trailing vortices generated by a Rushton turbine[J]. AIChE Journal, 2007,53(2): 523-526.

    [10] 宋戈. 渦輪槳攪拌槽內(nèi)湍流特性的三維 PIV實驗研究[D]. 北京:北京化工大學(xué), 2012.SONG G. The experiment investigation of fluid characteristics in the stirrer tank using stereo PIV[D]. Beijing: Beijing University of Chemical Technology, 2012.

    [11] SHARP K V, HILL D, TROOLIN D,et al. Volumetric three-component velocimetry measurements of the turbulent flow around a Rushton turbine[J]. Exp. Fluids, 2010, 48(1): 167-183.

    [12] HILL D, TROOLIN D, WALTERS G,et al. Volumetric 3-component velocimetry (V3V) measurements of the turbulent flow in stirred tank reactors[C]//Proc. 14th Intern. Symp. Applications of Laser Techniques to Fluid Mechanics. Lisbon, Portugal, 2008: 1-12.

    [13] 鮑蘇洋, 周勇軍, 王璐璐, 等. 渦輪槳攪拌槽內(nèi)流場特性的 V3V實驗[J]. 化工學(xué)報, 2016, 67(11): 4580-4586.BAO S Y, ZHOU Y J, WANG L L,et al. V3V study on flow field characteristics in a stirred vessel with a Rushton turbine impeller[J].CIESC Journal, 2016, 67(11): 4580-4586.

    [14] 張艷紅, 楊超, 毛在砂. 大渦模擬攪拌槽中的液相流動[J]. 化工學(xué)報, 2007, 58(10): 2474-2479.ZHANG Y H, YANG C, MAO Z S. Large eddy simulation of turbulent flow in a stirred tank[J]. Journal of Chemical Industry and Engineering(China), 2007, 58(10): 2474-2479.

    [15] 苗一, 潘家禎, 牛國瑞, 等. 多層槳攪拌槽內(nèi)的宏觀混合特性[J].華東理工大學(xué)學(xué)報(自然科學(xué)版), 2006, (3): 357-360.MIAO Y, PAN J Z, NIU G R,et al. Mixing in stirred tanks with multiple impellers[J]. Journal of East China University of Science and Technology (Natural Science Edition), 2006, (3): 357-360.

    [16] 李志鵬, 高正明. 渦輪槳攪拌槽內(nèi)單循環(huán)流動特性的大渦模擬[J].過程工程學(xué)報, 2007, 7(5): 900-904.LI Z P, GAO Z M. Large eddy simulation of single loop flow field in a Rushton impeller stirred tank[J]. The Chinese Journal of Process Engineering, 2007, 7(5): 900-904.

    [17] BAO Y, LI Z, GAO Z. PIV experiments and large eddy simulations of single loop flow fields in Rushton turbine stirred tank[J]. Chemical Engineering Science, 2011, 66(6): 1219-1231.

    [18] LI Z P, BAO S Y. Stereo-PIV experiments and large eddy simulations of flow fields in stirred tanks with Rushton and curved-blade turbines[J]. AIChE Journal, 2013, 59(10): 3986-4003.

    [19] BAKKER A, OSHINOWO L M. Modelling of turbulence in stirred vessels using large eddy simulation[J]. Trans. IChemE:Chem.Eng.Res., 2004, 82: 1169-1178.

    [20] 張兆順, 崔桂香, 許春曉.湍流大渦數(shù)值模擬的理論和應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2009: 236-247.ZHANG Z S, CUI G X, XU C X. Theory and Application of Turbulent Large Eddy Simulation[M]. Beijing: Tsinghua University Press, 2009: 236-247.

    [21] MENEVEAU C, KATZ J. Scale-invariance and turbulence models for large-eddy simulation[J]. Annual. Rev. Fluid Mech., 2000, 32:1-32.

    [22] DERKSEN J, VAN D, AKKER E A. Large eddy simulations on the flow driven by Rushton turbine[J]. AIChE Journal, 1999, 45(2):209-221.

    [23] 周勇軍, 盧源, 陳明濠, 等. 改進(jìn)型 INTER-MIG攪拌槽內(nèi)固液懸浮特性的數(shù)值模擬[J]. 過程工程學(xué)報, 2014, (5): 744-750.ZHOU Y J, LU Y, CHEN M H, et al. Numerical simulation of solid-liquid suspension characteristics in a stirred tank with improved INTER-MIG impeller[J]. The Chinese Journal of Process Engineering,2014, (5): 744-750.

    [24] ZHOU G, KRESTA S M. Impact of tank geometry on the maximum turbulence energy dissipation rate for impellers[J]. AIChE Journal,1996, 42(9): 2476-2490.

    [25] YOON H S, HILL D F, BALACHANDAR S, et al. Reynolds number scaling of flow in a Rushton turbine stirred tank(Ⅰ): Mean flow,circular jet and tip vortex scaling[J]. Chemical Engineering Science,2005, 60(12): 3169-3183.

    [26] YOON H S, BALACHANDAR S, MAN Y H. Large eddy simulation of flow in an unbaffled stirred tank for different Reynolds numbers[J].Physics of Fluids, 2009, 21(8): 137-1.

    date:2017-05-10.

    ZHOU Yongjun, associate professor, zhouyj@njtech.edu.cn.

    V3V study and large eddy simulation of turbulence characteristics in a stirred vessel with Rushton turbine impeller

    SHI Naijin1, ZHOU Yongjun1, BAO Suyang1, XIN Wei2, TAO Lanlan2
    (1College of Mechanical and Power Engineering,Nanjing Tech University,Nanjing211816,Jiangsu,China;2National Quality Supervision & Inspection Center of Chemical Equipment,Special Equipment Safety Supervision Inspection Institute of Jiangsu Province,Suzhou215600,Jiangsu,China)

    The volumetric three-component velocimetry (V3V) measurements and the large eddy simulation (LES)were used to study the flow field near the standard Rushton impeller in a stirred vessel. The results of V3V indicate that the Reynolds number has almost no influence on the normalized phase-resolved velocity and TKE when the flow field is absolutely turbulent. The structure and movement law of trailing vortices are discussed by reconstruction of the three-dimensional flow field with V3V. The distributions of radial, axial and tangential velocities in 30°cross-section behind the impeller were analyzed. The distributions of vortex pair by V3V were in good agreement with these by 2D-PIV (particle image velocimetry), however the vorticity by V3V is about 37.5%larger. LES method was used to simulate flow field in stirred tank with Rushton turbine impeller. The velocity distribution and trailing trajectory of LES results are in good agreement with the 2D-PIV and V3V data. Complete structure of the trailing vortices had been obtained and a vortex, which is similar with trailing vortices, was discovered near the up surface of the impeller.

    V3V; PIV; LES; trailing vortex; Rushton impeller

    TQ 022

    A

    0438—1157(2017)11—4069—10

    10.11949/j.issn.0438-1157.20170594

    2017-05-10收到初稿,2017-06-27收到修改稿。

    聯(lián)系人:周勇軍。

    施乃進(jìn)(1992—),男,碩士研究生。

    猜你喜歡
    尾渦大渦渦量
    不同B-V頻率下的飛機尾渦數(shù)值模擬研究
    高空巡航階段的飛機尾渦流場演化特性研究
    含沙空化對軸流泵內(nèi)渦量分布的影響
    基于壁面射流的下?lián)舯┝鞣欠€(wěn)態(tài)風(fēng)場大渦模擬
    基于激光雷達(dá)回波的動態(tài)尾渦特征參數(shù)計算
    自由表面渦流動現(xiàn)象的數(shù)值模擬
    干擾板作用下飛機尾渦流場近地演變機理研究
    軸流風(fēng)機葉尖泄漏流動的大渦模擬
    基于大渦模擬的旋風(fēng)分離器錐體結(jié)構(gòu)影響研究
    航態(tài)對大型船舶甲板氣流場的影響
    亚洲中文av在线| 国产精品,欧美在线| 日韩国内少妇激情av| 亚洲欧美精品综合久久99| 欧美成狂野欧美在线观看| 国产一区二区三区在线臀色熟女| 日韩成人在线观看一区二区三区| 又紧又爽又黄一区二区| 男男h啪啪无遮挡| 一卡2卡三卡四卡精品乱码亚洲| 午夜精品久久久久久毛片777| 日本a在线网址| 国产精品二区激情视频| 欧美又色又爽又黄视频| 亚洲欧美精品综合一区二区三区| 中文字幕高清在线视频| 久久午夜综合久久蜜桃| 欧美一级毛片孕妇| 免费在线观看亚洲国产| 国产真实乱freesex| 亚洲五月天丁香| 久久精品国产清高在天天线| 日本黄色视频三级网站网址| 日韩中文字幕欧美一区二区| 欧美日韩瑟瑟在线播放| 99国产精品99久久久久| 国产高清videossex| 欧美一级毛片孕妇| 欧洲精品卡2卡3卡4卡5卡区| 色尼玛亚洲综合影院| 中文字幕精品亚洲无线码一区 | 亚洲九九香蕉| 女生性感内裤真人,穿戴方法视频| 国产精品免费一区二区三区在线| 国产又色又爽无遮挡免费看| 女人高潮潮喷娇喘18禁视频| 999精品在线视频| 亚洲av五月六月丁香网| 少妇熟女aⅴ在线视频| 色综合婷婷激情| 国产aⅴ精品一区二区三区波| 国产亚洲欧美精品永久| 最新在线观看一区二区三区| 精品熟女少妇八av免费久了| 嫁个100分男人电影在线观看| 亚洲国产欧美一区二区综合| 国产成人精品久久二区二区91| 久久香蕉国产精品| 老汉色av国产亚洲站长工具| 亚洲熟妇中文字幕五十中出| 1024香蕉在线观看| 国产私拍福利视频在线观看| 两个人免费观看高清视频| 中文字幕人成人乱码亚洲影| 精品久久久久久,| 天天躁狠狠躁夜夜躁狠狠躁| av中文乱码字幕在线| 此物有八面人人有两片| 欧美+亚洲+日韩+国产| 国产爱豆传媒在线观看 | 久久国产精品影院| 亚洲专区字幕在线| 一区二区三区国产精品乱码| 欧美最黄视频在线播放免费| 黄色毛片三级朝国网站| 久久久久久久久免费视频了| 欧美大码av| 国产aⅴ精品一区二区三区波| 精品午夜福利视频在线观看一区| 国产区一区二久久| 欧美中文日本在线观看视频| 日韩一卡2卡3卡4卡2021年| 日韩一卡2卡3卡4卡2021年| 老司机午夜福利在线观看视频| 久久久国产成人精品二区| 在线av久久热| 在线av久久热| 久久久久国产一级毛片高清牌| 久久精品国产99精品国产亚洲性色| 国产亚洲欧美98| 午夜久久久在线观看| 欧美国产日韩亚洲一区| 亚洲男人的天堂狠狠| 色综合欧美亚洲国产小说| 国产精品亚洲av一区麻豆| 美女扒开内裤让男人捅视频| 日本三级黄在线观看| 99热这里只有精品一区 | 淫妇啪啪啪对白视频| 男人操女人黄网站| 国产精品久久久久久人妻精品电影| www.999成人在线观看| 黄色视频不卡| 黄色视频,在线免费观看| 久久久久久久午夜电影| 亚洲av美国av| 亚洲精品粉嫩美女一区| 女性生殖器流出的白浆| 免费看美女性在线毛片视频| 搡老妇女老女人老熟妇| 18禁美女被吸乳视频| 亚洲黑人精品在线| tocl精华| 在线观看日韩欧美| 亚洲成人久久爱视频| xxxwww97欧美| 国产精品亚洲av一区麻豆| 精品久久蜜臀av无| 91麻豆av在线| 精品午夜福利视频在线观看一区| 午夜福利成人在线免费观看| 欧美中文综合在线视频| 日韩欧美在线二视频| 国产人伦9x9x在线观看| 精品一区二区三区av网在线观看| 看片在线看免费视频| 91麻豆av在线| 亚洲人成伊人成综合网2020| 在线观看免费午夜福利视频| 少妇的丰满在线观看| 色尼玛亚洲综合影院| 男女那种视频在线观看| 我的亚洲天堂| 白带黄色成豆腐渣| 999精品在线视频| 91字幕亚洲| 在线观看66精品国产| www日本在线高清视频| av电影中文网址| 欧美激情久久久久久爽电影| 一级黄色大片毛片| 一二三四在线观看免费中文在| 日本黄色视频三级网站网址| 不卡一级毛片| 亚洲精品一区av在线观看| 国产精品香港三级国产av潘金莲| www日本黄色视频网| 国产精品,欧美在线| 19禁男女啪啪无遮挡网站| 国产精品av久久久久免费| 亚洲性夜色夜夜综合| 欧美性猛交╳xxx乱大交人| 中文字幕另类日韩欧美亚洲嫩草| 久久天躁狠狠躁夜夜2o2o| 一级片免费观看大全| 亚洲国产欧美一区二区综合| 一级片免费观看大全| 一边摸一边抽搐一进一小说| 免费人成视频x8x8入口观看| 1024视频免费在线观看| 国产亚洲av嫩草精品影院| 久久久国产成人免费| 麻豆国产av国片精品| 国产视频内射| 国产精品亚洲av一区麻豆| 最近最新中文字幕大全电影3 | 久久久久久久久中文| 搞女人的毛片| 国产精品亚洲一级av第二区| 久久精品人妻少妇| 成人三级黄色视频| 久久精品国产清高在天天线| 国产精品自产拍在线观看55亚洲| 精品免费久久久久久久清纯| 久久天躁狠狠躁夜夜2o2o| 国产真实乱freesex| 免费在线观看亚洲国产| 亚洲第一欧美日韩一区二区三区| 可以免费在线观看a视频的电影网站| 熟女少妇亚洲综合色aaa.| 国产激情欧美一区二区| 搡老熟女国产l中国老女人| 免费在线观看视频国产中文字幕亚洲| 亚洲国产精品合色在线| 亚洲精品国产区一区二| 一进一出抽搐gif免费好疼| 2021天堂中文幕一二区在线观 | 操出白浆在线播放| 人妻久久中文字幕网| 久久人妻av系列| 好看av亚洲va欧美ⅴa在| 久热爱精品视频在线9| 国产一级毛片七仙女欲春2 | 麻豆成人午夜福利视频| 在线看三级毛片| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品av麻豆狂野| 人人妻人人看人人澡| 日韩国内少妇激情av| 夜夜看夜夜爽夜夜摸| 99国产极品粉嫩在线观看| 啦啦啦免费观看视频1| 色综合亚洲欧美另类图片| 听说在线观看完整版免费高清| 99久久无色码亚洲精品果冻| 亚洲av熟女| 午夜免费激情av| 一级作爱视频免费观看| 免费在线观看成人毛片| 欧美黄色淫秽网站| 午夜久久久在线观看| 成人一区二区视频在线观看| 亚洲自偷自拍图片 自拍| 黄网站色视频无遮挡免费观看| 国产av在哪里看| 国产精品久久久人人做人人爽| 欧美日韩亚洲综合一区二区三区_| av视频在线观看入口| 免费搜索国产男女视频| 别揉我奶头~嗯~啊~动态视频| 久久久久久亚洲精品国产蜜桃av| 亚洲av电影不卡..在线观看| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区mp4| 大型黄色视频在线免费观看| 国产在线观看jvid| 在线免费观看的www视频| 丝袜美腿诱惑在线| 国产成人精品久久二区二区91| 一级作爱视频免费观看| 91麻豆av在线| 欧美三级亚洲精品| 不卡一级毛片| 丁香欧美五月| 国产黄a三级三级三级人| 中文字幕另类日韩欧美亚洲嫩草| 成人午夜高清在线视频 | 一卡2卡三卡四卡精品乱码亚洲| 看黄色毛片网站| 白带黄色成豆腐渣| 国产高清视频在线播放一区| 美女午夜性视频免费| 久久亚洲精品不卡| 亚洲中文av在线| а√天堂www在线а√下载| 大香蕉久久成人网| 精品久久蜜臀av无| 脱女人内裤的视频| 国产伦一二天堂av在线观看| 欧美日韩精品网址| 亚洲真实伦在线观看| 免费在线观看完整版高清| 91av网站免费观看| 欧美色视频一区免费| 久久久久国产一级毛片高清牌| 窝窝影院91人妻| 欧美日韩瑟瑟在线播放| 18禁观看日本| 国产精品久久久久久精品电影 | 熟女少妇亚洲综合色aaa.| 中文字幕高清在线视频| 欧美成人午夜精品| 日本免费一区二区三区高清不卡| 国产又爽黄色视频| 久久久国产精品麻豆| 男人的好看免费观看在线视频 | 欧美日韩乱码在线| 久久久久九九精品影院| 91麻豆av在线| 国产精品一区二区免费欧美| 人妻丰满熟妇av一区二区三区| 美女国产高潮福利片在线看| 激情在线观看视频在线高清| 亚洲人成伊人成综合网2020| 亚洲精品美女久久久久99蜜臀| av电影中文网址| 一个人观看的视频www高清免费观看 | 精品一区二区三区av网在线观看| 手机成人av网站| 国产又爽黄色视频| 日本免费a在线| 美女 人体艺术 gogo| 黑人操中国人逼视频| 女人高潮潮喷娇喘18禁视频| 久久久久久大精品| 亚洲人成伊人成综合网2020| 国产亚洲精品av在线| 日韩免费av在线播放| 亚洲精品粉嫩美女一区| 国产真实乱freesex| 波多野结衣巨乳人妻| 午夜福利视频1000在线观看| 很黄的视频免费| 亚洲午夜理论影院| 亚洲精品一区av在线观看| 午夜福利免费观看在线| 欧美日韩一级在线毛片| 欧美最黄视频在线播放免费| 欧美亚洲日本最大视频资源| 美女高潮喷水抽搐中文字幕| 久久九九热精品免费| 久久久久久久久免费视频了| 亚洲熟女毛片儿| 国产精品二区激情视频| 国产真实乱freesex| 黄色a级毛片大全视频| 午夜两性在线视频| 亚洲久久久国产精品| 别揉我奶头~嗯~啊~动态视频| 看免费av毛片| 日韩免费av在线播放| 久久中文字幕人妻熟女| 在线av久久热| 国产成+人综合+亚洲专区| 无遮挡黄片免费观看| 国产高清有码在线观看视频 | 亚洲国产精品999在线| 久久精品国产亚洲av香蕉五月| 少妇裸体淫交视频免费看高清 | 日韩av在线大香蕉| 很黄的视频免费| 国产熟女xx| 精华霜和精华液先用哪个| 免费搜索国产男女视频| 欧美大码av| 国产爱豆传媒在线观看 | 欧美日韩亚洲综合一区二区三区_| 亚洲第一欧美日韩一区二区三区| 欧美日韩瑟瑟在线播放| 亚洲美女黄片视频| 亚洲精品中文字幕一二三四区| 国产精品久久电影中文字幕| 少妇 在线观看| 国产私拍福利视频在线观看| 后天国语完整版免费观看| 国产精品电影一区二区三区| 12—13女人毛片做爰片一| 久久久久久大精品| 十分钟在线观看高清视频www| 国产免费男女视频| 午夜福利18| 久久精品国产亚洲av高清一级| 久久久国产成人精品二区| 激情在线观看视频在线高清| 亚洲激情在线av| 好男人在线观看高清免费视频 | 久久中文看片网| 香蕉国产在线看| 窝窝影院91人妻| 亚洲精品在线观看二区| 国产伦人伦偷精品视频| www国产在线视频色| 女性被躁到高潮视频| 日本 av在线| 俺也久久电影网| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区三区视频了| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 欧美绝顶高潮抽搐喷水| 国产精品久久久人人做人人爽| 制服诱惑二区| 人人妻人人看人人澡| 美女大奶头视频| 视频在线观看一区二区三区| 精品第一国产精品| 精品电影一区二区在线| av在线天堂中文字幕| 欧美日韩亚洲国产一区二区在线观看| 757午夜福利合集在线观看| 一区二区日韩欧美中文字幕| 久久久久九九精品影院| 日日摸夜夜添夜夜添小说| 神马国产精品三级电影在线观看 | 19禁男女啪啪无遮挡网站| 免费在线观看黄色视频的| 亚洲av成人一区二区三| 成人三级黄色视频| 国产精品久久久久久人妻精品电影| 天堂影院成人在线观看| 日韩中文字幕欧美一区二区| 哪里可以看免费的av片| www国产在线视频色| 久久国产精品影院| av有码第一页| 成人国语在线视频| 一二三四在线观看免费中文在| 在线天堂中文资源库| 免费在线观看成人毛片| 亚洲av五月六月丁香网| 亚洲精品色激情综合| 中文字幕精品亚洲无线码一区 | 日韩成人在线观看一区二区三区| 欧美在线一区亚洲| 久久午夜综合久久蜜桃| 成人国产一区最新在线观看| 搡老妇女老女人老熟妇| 成人18禁在线播放| 国产99白浆流出| 俄罗斯特黄特色一大片| 人妻久久中文字幕网| 动漫黄色视频在线观看| 露出奶头的视频| 亚洲一区高清亚洲精品| 又紧又爽又黄一区二区| 19禁男女啪啪无遮挡网站| 狂野欧美激情性xxxx| 日本成人三级电影网站| 国产黄a三级三级三级人| 女人被狂操c到高潮| 亚洲国产高清在线一区二区三 | 人成视频在线观看免费观看| 亚洲在线自拍视频| 91字幕亚洲| 人成视频在线观看免费观看| 夜夜躁狠狠躁天天躁| 久久久国产欧美日韩av| 国产成人精品无人区| 国产激情欧美一区二区| 国产又爽黄色视频| 一本久久中文字幕| 又黄又粗又硬又大视频| 在线永久观看黄色视频| 操出白浆在线播放| 精品福利观看| 国产一区二区在线av高清观看| 亚洲欧美精品综合久久99| 一区二区三区激情视频| 嫩草影院精品99| 在线播放国产精品三级| 亚洲全国av大片| 啪啪无遮挡十八禁网站| 中文字幕久久专区| 欧美av亚洲av综合av国产av| a在线观看视频网站| 亚洲国产欧美网| 中文字幕人成人乱码亚洲影| 欧美日韩一级在线毛片| 日韩精品中文字幕看吧| 啦啦啦 在线观看视频| 精品久久久久久成人av| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| 亚洲精品国产一区二区精华液| 999久久久国产精品视频| 亚洲七黄色美女视频| 男女那种视频在线观看| 好男人电影高清在线观看| 夜夜爽天天搞| 99国产精品99久久久久| 国产欧美日韩一区二区三| bbb黄色大片| 婷婷精品国产亚洲av| 在线免费观看的www视频| 麻豆成人午夜福利视频| 久久久精品国产亚洲av高清涩受| 国产精品98久久久久久宅男小说| 女性生殖器流出的白浆| 久久久久精品国产欧美久久久| 国产亚洲精品第一综合不卡| 91国产中文字幕| 又黄又爽又免费观看的视频| 久久精品影院6| 丰满人妻熟妇乱又伦精品不卡| 久久欧美精品欧美久久欧美| 亚洲七黄色美女视频| 国产精品野战在线观看| 国产精品二区激情视频| 国产1区2区3区精品| 嫩草影院精品99| 岛国在线观看网站| a在线观看视频网站| 午夜福利高清视频| e午夜精品久久久久久久| 一区福利在线观看| 亚洲人成网站高清观看| 亚洲av电影在线进入| 欧美日韩瑟瑟在线播放| 黄片大片在线免费观看| 此物有八面人人有两片| 国产亚洲欧美98| 亚洲成人精品中文字幕电影| 精品久久久久久久久久久久久 | 黄色丝袜av网址大全| 亚洲精品美女久久久久99蜜臀| 9191精品国产免费久久| 色播亚洲综合网| 国产一区二区三区视频了| 国产日本99.免费观看| 亚洲精品久久成人aⅴ小说| 黄色a级毛片大全视频| 日韩视频一区二区在线观看| 黄色成人免费大全| 欧美日韩一级在线毛片| 美女免费视频网站| 国产在线精品亚洲第一网站| 亚洲九九香蕉| 一边摸一边做爽爽视频免费| x7x7x7水蜜桃| 久久精品国产亚洲av香蕉五月| 啦啦啦韩国在线观看视频| 制服丝袜大香蕉在线| 国产亚洲精品久久久久久毛片| 无人区码免费观看不卡| 男女下面进入的视频免费午夜 | 91麻豆精品激情在线观看国产| 窝窝影院91人妻| 亚洲精品中文字幕一二三四区| 亚洲国产高清在线一区二区三 | 欧美一区二区精品小视频在线| 可以免费在线观看a视频的电影网站| 久久精品国产清高在天天线| 色综合婷婷激情| 一级片免费观看大全| 老熟妇仑乱视频hdxx| 免费看日本二区| 免费电影在线观看免费观看| 在线观看舔阴道视频| 久久精品成人免费网站| 天天躁夜夜躁狠狠躁躁| 亚洲精品国产一区二区精华液| 老汉色av国产亚洲站长工具| 天天一区二区日本电影三级| 男人操女人黄网站| 97人妻精品一区二区三区麻豆 | 国产蜜桃级精品一区二区三区| 欧美最黄视频在线播放免费| 成人午夜高清在线视频 | 国产单亲对白刺激| 88av欧美| 不卡av一区二区三区| 国产亚洲av嫩草精品影院| 国产三级在线视频| 久久国产精品人妻蜜桃| 国内毛片毛片毛片毛片毛片| 亚洲电影在线观看av| 国产区一区二久久| 日韩欧美三级三区| 无限看片的www在线观看| 午夜福利视频1000在线观看| 欧美av亚洲av综合av国产av| 夜夜躁狠狠躁天天躁| 午夜亚洲福利在线播放| 在线观看66精品国产| 免费在线观看亚洲国产| 1024香蕉在线观看| 一本综合久久免费| 欧美成狂野欧美在线观看| 亚洲国产精品999在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美一区二区三区黑人| 悠悠久久av| 欧美另类亚洲清纯唯美| 久久久久久久午夜电影| 久热爱精品视频在线9| 亚洲欧美激情综合另类| 久久中文看片网| 亚洲专区字幕在线| 真人一进一出gif抽搐免费| 美女免费视频网站| 可以在线观看毛片的网站| 人成视频在线观看免费观看| 成人国产一区最新在线观看| 国产人伦9x9x在线观看| 精品国产乱码久久久久久男人| 亚洲中文字幕日韩| 久久草成人影院| 亚洲成av人片免费观看| 日韩成人在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 成人午夜高清在线视频 | 欧美日韩亚洲国产一区二区在线观看| 亚洲精品国产精品久久久不卡| 国产黄片美女视频| 精品一区二区三区四区五区乱码| 1024香蕉在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 俺也久久电影网| 午夜福利18| 免费在线观看日本一区| 一级片免费观看大全| 人人妻,人人澡人人爽秒播| 九色国产91popny在线| 一夜夜www| 日韩免费av在线播放| 国产片内射在线| 中文字幕另类日韩欧美亚洲嫩草| 丝袜在线中文字幕| 午夜成年电影在线免费观看| a级毛片a级免费在线| 亚洲中文字幕日韩| 婷婷丁香在线五月| 亚洲无线在线观看| 亚洲 欧美 日韩 在线 免费| 婷婷丁香在线五月| 午夜福利免费观看在线| 精品国产亚洲在线| 校园春色视频在线观看| 免费高清在线观看日韩| 亚洲一区二区三区色噜噜| 国产激情久久老熟女| 黄色视频,在线免费观看| 身体一侧抽搐| 国产精品 欧美亚洲| 免费看美女性在线毛片视频| 成人av一区二区三区在线看| 黄色丝袜av网址大全| 一夜夜www| 嫩草影视91久久| 国产在线精品亚洲第一网站| 国产成人精品久久二区二区91| 日韩欧美国产一区二区入口| 人人澡人人妻人| 母亲3免费完整高清在线观看| 最好的美女福利视频网| 欧美最黄视频在线播放免费| 久久精品成人免费网站| 一级a爱片免费观看的视频| 一区二区三区国产精品乱码| 色精品久久人妻99蜜桃| 亚洲国产精品sss在线观看| 久久 成人 亚洲| 在线观看免费日韩欧美大片| 午夜a级毛片| 久久久久免费精品人妻一区二区 |