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

    圓盤反應(yīng)器液膜表面更新數(shù)值模擬

    2015-08-20 07:31:42鄧斌戴干策
    化工學(xué)報 2015年4期
    關(guān)鍵詞:附壁液膜圓盤

    鄧斌,戴干策

    (華東理工大學(xué)化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗(yàn)室,上海 200237)

    引 言

    圓盤反應(yīng)器中水平軸旋轉(zhuǎn)帶動垂直圓盤,從釜體內(nèi)拖曳液體,在圓盤表面形成不斷更新的液膜。圓盤反應(yīng)器早期用于廢水處理及血液中氧傳遞。20世紀(jì)60年代應(yīng)用于合成纖維聚酯,并獲得成功[1]。90年代以來,應(yīng)用范圍進(jìn)一步擴(kuò)大,包括檸檬酸生產(chǎn)[2]、光催化[3]、精細(xì)化工[4]、石油回收[5]、微生物燃料電池[6]等。工業(yè)應(yīng)用推動了圓盤反應(yīng)器的研究。Bintanja 等[7]、Zeevalkink 等[8]、Afanasiev 等[5]研究了圓盤表面膜厚的影響因素,獲得了膜厚分布規(guī)律。Zeevalkink 等[9]、Kubsad 等[10]考察了圓盤反應(yīng)器的傳質(zhì)特性。王良生[11]、周賢爵[12]對聚酯工業(yè)中圓盤反應(yīng)器的持液量和軸功率等進(jìn)行了系統(tǒng)研究。盡管對圓盤反應(yīng)器的流動與傳質(zhì)研究有所進(jìn)展,但液膜表面更新這一重要特性至今研究甚少。

    表面更新概念始于Danckwerts[13],考察氣液界面?zhèn)髻|(zhì),認(rèn)為液相主體內(nèi)的湍流旋渦或流體微團(tuán)遷移到相界面,進(jìn)行物質(zhì)和能量交換,一段時間后界面微團(tuán)流體被新的微團(tuán)取代,進(jìn)入主體內(nèi),完成更新。這一概念從20世紀(jì)50年代初期提出就曾經(jīng)引起廣泛關(guān)注。20世紀(jì)60~70年代,對低黏物系,Perlmutter[14]、Harriott[15]、Koppel[16]、Chung 等[17]諸多學(xué)者對表面更新理論提出了不同程度的修改,概念上都更為合理,但是對更新頻率、旋渦表面停留時間這些關(guān)鍵參數(shù)仍未給出簡要的計(jì)算方法。

    對中高黏度物系,表面更新特性更為重要。Murakami 等[18]在1972年根據(jù)拉絲模型估算聚酯圓盤反應(yīng)器內(nèi)液膜的表面更新頻率為0.8~80 s-1。Biesenberger 等[19]在1983年提出隨機(jī)表面更新模型和選擇性表面更新模型應(yīng)用于聚酯圓盤反應(yīng)器。潘勤敏[20]在1987年采用示蹤方法分析圓盤反應(yīng)器內(nèi)界面更新,提出以液相微元表面停留時間的倒數(shù)作為更新頻率,同時得到更新頻率S=1/tr=2N,認(rèn)為表面更新正比于轉(zhuǎn)速。

    周賢爵[12]在2001年研究籠式反應(yīng)器流體動力學(xué)時,將表面更新頻率與液膜變形聯(lián)系,利用激光多普勒實(shí)驗(yàn)測定圓盤液膜速度分布,計(jì)算更新頻率,并給出量綱1 關(guān)聯(lián)式。這是一種合理而且有意義的嘗試,但該研究未能探討圓盤結(jié)構(gòu)與表面更新的關(guān)系。實(shí)際上,聚酯圓盤表面開窗,因而存在兩種類型的液膜[21]:一種是附壁膜(wall-bounded film),僅一個自由面;另一種是自由膜(free film),無壁面支撐,存在兩個自由面。兩種液膜行為上的差異,特別是表面更新,有待研究。

    本研究在分析低黏/高黏物系表面更新的基礎(chǔ)上對圓盤層流液膜建立表面更新頻率與變形率的關(guān)系,通過數(shù)值模擬圓盤反應(yīng)器流場特性計(jì)算更新頻率并探討相關(guān)影響因素,從表面更新的角度將自由膜與附壁膜對比研究,進(jìn)一步揭示圓盤表面開窗對過程強(qiáng)化的意義,并為圓盤結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供量化分析的基礎(chǔ)。

    1 表面更新理論分析

    1.1 液膜表面更新與變形率

    Danckwerts 表面更新理論中,液體表面視為不同停留時間、不同濃度液體微團(tuán)“鑲嵌”而成,液體表面周期性地被充分混合的主體相流體微團(tuán)(與氣液界面接觸的湍流渦)取代。

    單位時間內(nèi)表面被新鮮液體取代的次數(shù)即為表面更新頻率,它可用黏性耗散范圍內(nèi)湍流渦的特征參數(shù)旋渦特征速度ue、旋渦長度尺寸l定義[22]

    在表面更新拉伸模型中[23]以新鮮表面產(chǎn)生的速率作為表面更新的頻率,因此

    式中,r為膜表面平均交換速率,m2·s-1;A為液膜與氣相接觸的總面積,m2。

    對于層流液膜(含波動層流),不存在湍流渦,混合與表面更新依賴剪切變形與拉伸變形,表面更新頻率由變形率決定。下面試就一維、二維液膜的表面更新做簡單推演。

    1.2 一維液膜流動表面更新

    如圖1所示,一寬為W的無限長平板,當(dāng)t=0時,液膜以初速度u0在重力作用下加速向下流動,經(jīng)過L達(dá)到充分發(fā)展,此時液膜表面流速為u。

    圖1 液膜一維流動Fig.1 One-dimensional flow of liquid film

    參照式(1)、式(2),表面更新頻率定義為

    式中,u、L分別為流場的特征速度和特征長度。

    將式(3)分子、分母同乘寬度W,則有

    式中,ΔA為Δt時間內(nèi)由拉伸流動速度變化所致液膜增加的面積,即 (u-u0)WΔt。

    于是所得的液膜表面更新頻率為

    當(dāng)L→ 0時,局部表面更新頻率為SP

    1.3 二維液膜流動與表面更新

    對二維液膜拉伸流動,如圖2所示。液膜在A1B1處沿x方向的速度是u,在C1D1處沿x方向的速度是液膜在A1D1處沿y方向的速度是v,在B1C1處沿y方向的速度是不考慮u隨y的變化和v隨x的變化,經(jīng)過δt時間后,由流體質(zhì)點(diǎn)組成的流體線C1D1較A1B1多移動距離,B1C1較A1D1向下多移動距離(這里為負(fù)),因此A1B1C1D1的新位置A1B2C2D2仍為矩形,但邊長不同,微團(tuán)發(fā)生了形變,其特征量是單位時間內(nèi)長度的相對變化,為線變形速率,以ε表示。

    圖2 液膜二維流動Fig.2 Two-dimensional flow of liquid film

    在x方向有

    在y方向有

    對于不可壓縮流體,連續(xù)性方程為

    從以上分析可知,液膜的更新頻率與液膜表面的線變形速率直接相關(guān)。由此二維局部更新頻率定義為

    圓盤表面液膜的平均更新頻率則為

    2 計(jì)算方法

    根據(jù)圓盤反應(yīng)器內(nèi)氣液兩相的流動特點(diǎn)和接觸方式,采用Hirt 等[24]提出的VOF 模型進(jìn)行計(jì)算。該方法的基本原理:通過分析網(wǎng)格單元中流體和網(wǎng)格體積比函數(shù)確定自由面,追蹤流體界面的變化,而非追蹤自由液面上質(zhì)點(diǎn)的運(yùn)動。VOF 方法可以處理自由面重入等強(qiáng)非線性現(xiàn)象,所需計(jì)算時間短、存儲量少。若液相體積分?jǐn)?shù)α=1,說明該單元全部為指定相流體所占據(jù);若α=0,該單元內(nèi)無指定相流體;當(dāng)0<α<1時,該單元稱為交界面單元。

    對于三維流動,VOF 模型的連續(xù)方程、動量方程、體積分?jǐn)?shù)連續(xù)方程分別為

    并有

    因表面張力作用引起的動量方程的源項(xiàng)為

    式中,ρ為密度,u為速度矢量,μ為流體動力黏度,αi為第i相的體積分?jǐn)?shù),κ為表面曲率,下角標(biāo)1、2 分別代表空氣和液體。

    計(jì)算中將圓盤和軸表面的邊界設(shè)置為旋轉(zhuǎn)壁面,壁面邊界條件為無滑移,自由面通過求解VOF模型方程得到。本研究不考慮液位變化的影響,圓盤浸沒深度始終保持為圓盤的半徑,即保持液位與轉(zhuǎn)軸在同一平面。

    采用Gambit 進(jìn)行網(wǎng)格劃分,ANSYS FLUENT商業(yè)軟件求解。圓盤反應(yīng)器網(wǎng)格劃分和計(jì)算區(qū)域如 圖3所示,均為六面體網(wǎng)格。為了減少網(wǎng)格數(shù),節(jié)約計(jì)算時間,根據(jù)圓盤兩側(cè)流場關(guān)于圓盤對稱的特點(diǎn),對于普通圓盤只將圓盤一側(cè)的區(qū)間選取為計(jì)算域,如圖4所示,并進(jìn)行驗(yàn)證,見表1。對于開窗圓盤,兩側(cè)均為計(jì)算區(qū)域??紤]到研究重點(diǎn)是圓盤表面液膜更新,對圓盤附近區(qū)域網(wǎng)格按邊界層加密方法進(jìn)行加密,保證平均膜厚內(nèi)有10 個網(wǎng)格,網(wǎng)格數(shù)控制在10 萬~20 萬個之間,兼顧計(jì)算時間和計(jì)算精度。計(jì)算過程中改變圓盤半徑、轉(zhuǎn)速、結(jié)構(gòu)和黏度等參數(shù),本研究模擬范圍為9<Re<3 00,取層流模型。借有限體積法離散控制方程;壁面服從無滑移條件,壓力-速度耦合方程的求解采用PISO 方法;對氣液界面的追蹤采用精度較高的 Geo-Reconstruct 界面重構(gòu)技術(shù)??紤]重力和表面張力影響,參考壓力為標(biāo)準(zhǔn)大氣壓。根據(jù)操作條件的不同,計(jì)算的時間步長在10-5~10-4s 之間,待計(jì)算穩(wěn)定(膜厚和功率波動小于5%)再提取數(shù)據(jù)。

    圖3 對稱與非對稱圓盤的網(wǎng)格劃分和計(jì)算域Fig.3 Mesh generation and computational domain

    圖4 CFD 計(jì)算模擬單盤基本外形參數(shù)Fig.4 Computational simulation parameters of single disc

    3 模型驗(yàn)證

    從網(wǎng)格無關(guān)性、液相流型、膜厚和液膜速度4個方面進(jìn)行模型驗(yàn)證。網(wǎng)格無關(guān)性檢驗(yàn)見表1,不同網(wǎng)格數(shù)的非對稱網(wǎng)格、膜厚和功率相差均小于3%。垂直旋轉(zhuǎn)圓盤表面部分浸沒于液體中,液相流型為經(jīng)典的Stewartson[25]流型,與相同狀態(tài)下的CFD 模擬結(jié)果比較如圖5、圖6所示,兩者流型 一致。

    圖5 Stewartson 流型俯視圖Fig.5 Top view of Stewartson flow pattern

    圖6 CFD 模擬跡線俯視圖Fig.6 Top view of trace of CFD simulation (y=-10 mm)

    Afanasiev 等[5]在2008年通過有限元方法求解四階非線性偏微分方程研究圓盤表面液膜的分布規(guī)律,發(fā)現(xiàn)半徑一定時膜厚隨角度增加而減小,最后趨于恒定膜厚h∞,并得到了恒定膜厚的表達(dá)式。本研究模擬半徑R=150 mm、μ=1000 mPa·s、N=90 r·min-1的旋轉(zhuǎn)圓盤,當(dāng)r=100 mm 時膜厚隨角度的變化如圖7所示,膜厚隨角度增加而減小,h∞(膜厚隨角度不再減小時的厚度)接近Afanasiev 等的研究結(jié)果。

    膜厚實(shí)驗(yàn)驗(yàn)證采用黏度為253 mPa·s 的甘油 水溶液,轉(zhuǎn)軸水平,圓盤豎直,圓盤半徑75 mm,盤面1/2 浸沒在水槽中,轉(zhuǎn)速140 r·min-1。水槽足夠大,壁效應(yīng)可以忽略。實(shí)驗(yàn)過程中,圓盤轉(zhuǎn)速穩(wěn)定后,用德國Precitec Optronik GmbH 公司的CHR E測量系統(tǒng)進(jìn)行非接觸式在線動態(tài)測量圓盤表面局部膜厚,所用探頭型號RB200050。探頭的光線在液膜表面形成的光點(diǎn)直徑12 μm,可認(rèn)為是單點(diǎn)測試。液膜厚度方向的分別率為0.1 μm,采集頻率為1000 Hz,膜厚穩(wěn)定后取平均值。

    表1 對稱與非對稱圓盤的模擬參數(shù)和結(jié)果Table 1 Simulation parameters and results of symmetric and asymmetric disc

    圖7 r=100 mm 時膜厚切向分布Fig.7 Tangential distribution of film thickness at r=100 mm

    圖8 θ=π/2 時膜厚徑向分布Fig.8 Radial distribution of film thickness at θ=π/2

    實(shí)驗(yàn)主要測量了圓盤表面液膜厚度徑向分布。以液面作為x軸。固定角度為π/2,分別取半徑r為25、35、45、55、65、68 mm 的6 個點(diǎn)進(jìn)行膜厚測量。CFD 模擬參數(shù)完全按照實(shí)際測量工況設(shè)計(jì),盡量減少不必要的誤差,數(shù)據(jù)采集和處理也同實(shí)驗(yàn)一致,取相同位置的膜厚數(shù)據(jù)進(jìn)行對比,兩者的對比結(jié)果如圖8所示。

    對比實(shí)驗(yàn)和CFD 計(jì)算結(jié)果可知,在圓盤表面膜厚較薄或較厚的地方實(shí)驗(yàn)測量結(jié)果和CFD 計(jì)算結(jié)果吻合較好,誤差小于0.2 mm;在膜厚由薄向厚過渡的地方兩者誤差小于0.5 mm。這是因?yàn)樵谶^渡區(qū)域膜厚梯度較大,測量探頭位置的移動和網(wǎng)格劃分都會造成兩者誤差增加。文中最大誤差為0.5 mm,能夠滿足計(jì)算分析要求。

    James[26]對平板抽出過程中液膜內(nèi)速度分布進(jìn)行了理論推導(dǎo),認(rèn)為在液膜內(nèi)沿平板方向存在向上和向下的速度分布,如圖9所示。圓盤旋轉(zhuǎn)時,液膜上升段與平板抽出模型類似,因此將上升段液 膜內(nèi)速度分布與理論解進(jìn)行對比,如圖9和圖10所示。

    圖9 理論推導(dǎo)液膜內(nèi)速度分布Fig.9 Distribution of velocity theoretically

    圖10 CFD 模擬上升段液膜內(nèi)速度分布Fig.10 Distribution of velocity by CFD

    總體而言,CFD 模擬結(jié)果與文獻(xiàn)和實(shí)驗(yàn)對比,表明VOF 模型既能從宏觀角度模擬液相流型特征也能從微觀角度模擬圓盤表面膜厚變化,因此本研究數(shù)值模擬的結(jié)果是可信的。

    4 結(jié)果與討論

    4.1 圓盤液膜更新頻率空間分布

    表面更新與速度分布密切相關(guān),因此先對圓盤表面速度分布進(jìn)行分析。對圓盤半徑R=150mm、液體黏度μ=1000 mPa·s、轉(zhuǎn)速N=90 r·min-1進(jìn)行數(shù)值模擬,研究圓盤表面液膜的速度分布。以網(wǎng)格內(nèi)液相分率VOF=0.99 等值面作為液膜自由面,即氣液接觸面。速度矢量分布如圖11所示,根據(jù)速 度方向和大小的差異將圓盤表面分為起始區(qū)、加速區(qū)、勻速區(qū)。

    圖11 表面液膜速度矢量圖Fig.11 Liquid film velocity vector of rotating disc

    圖12 液膜更新頻率分布Fig.12 Distribution of surface renewal frequency

    圖13 液膜更新頻率隨角度和半徑的變化Fig.13 Relationship between liquid surface renewal frequency and angle and radius

    4.2 圓盤液膜平均更新頻率及其影響因素

    對于工業(yè)過程,更加關(guān)注液膜更新頻率的宏觀特征,因此先對圓盤平均更新頻率進(jìn)行研究,考察轉(zhuǎn)速、半徑等因素的影響。

    將平均更新頻率S0量綱1 化可得

    式中,k為比例系數(shù)。

    對R和ω進(jìn)行量綱分析,得到:a=-0 .10,b=0.39。

    圖14 液膜平均更新頻率擬合Fig.14 Linear fitting of average liquid film surface renewal frequency

    擬合結(jié)果表明S0∝ω1.188,增加轉(zhuǎn)速顯著強(qiáng)化更新。液膜平均更新頻率S0∝R0.678,增加半徑更有利于加快液膜更新。在研究的黏度范圍內(nèi),液膜平均更新頻率S0∝μ-0.39,液膜更新頻率對黏度敏感,隨黏度增加急劇減小。

    4.3 圓盤結(jié)構(gòu)與表面更新

    圖15 圓盤結(jié)構(gòu)特征Fig.15 Diagram of different discs structure

    對圓盤半徑均為75 mm、液體黏度μ=1000 mPa·s、轉(zhuǎn)速N=90 r·min-1進(jìn)行模擬,物系和操作條件不變,分別從開窗與不開窗、開窗形狀和開窗大小3 個方面研究圓盤結(jié)構(gòu)與表面更新的關(guān)系。為了方便表述,扇形窗圓盤記為圓盤A,是聚酯工業(yè)中最常見的開窗圓盤;普通圓盤記為圓盤B;開6 個大圓形窗的圓盤記為圓盤C;開12 個小圓形窗的圓盤記為圓盤D。結(jié)構(gòu)分別如圖15所示,3 個開窗圓盤具有相同的開窗面積。

    4.3.1 開窗圓盤與普通圓盤液膜更新對比 普通圓盤B 與開窗圓盤A 比較,表面更新頻率分布如圖16和圖17所示。由于開窗區(qū)域的存在,圓盤液膜表面更新頻率明顯加快,特別是在加速區(qū)和加速區(qū)與勻速區(qū)的過渡區(qū)域,開窗對于表面更新強(qiáng)化作用明顯。計(jì)算得到圓盤B 平均液膜更新頻率SB=1.68 s-1;圓盤A 平均液膜更新頻率SA=2.35 s-1,比圓盤B 高40%。

    開窗區(qū)域自由膜更新頻率如圖18所示,3 個開窗區(qū)域自由膜平均更新頻率為SAF=2.63 s-1。自由膜更新頻率沿角度依次減小,這與普通圓盤附壁膜更新頻率分布規(guī)律類似。自由膜上部更新頻率明顯 大于下部,因?yàn)樵陂_窗區(qū)域沒有壁面作用,重力作用明顯,液膜向下加速,更新頻率加快,液膜在窗口下降一段距離后達(dá)到勻速,更新頻率減慢。

    圖16 圓盤B 表面更新頻率分布Fig.16 Distribution of disc B surface renewal frequency

    圖17 圓盤A 表面更新頻率分布Fig.17 Distribution of disc A surface renewal frequency

    定量比較,圓盤表面半徑r=45 mm 更新頻率隨角度θ的變化規(guī)律如圖19所示。半徑r=45 mm的圓剛好經(jīng)過開窗區(qū)域,開窗區(qū)域自由膜更新頻率明顯高于對應(yīng)的普通圓盤附壁膜。普通圓盤附壁膜更新頻率在加速區(qū)出現(xiàn)最大值之后一直減小;開窗圓盤在每個窗口位置自由膜更新頻率都會出現(xiàn)一個峰值,峰值隨角度增加而減小。并且自由膜周圍的附壁膜更新頻率也得到了提高,圓盤A 附壁膜平均更新頻率為SAW=2.28 s-1,比圓盤B 高35%。

    圖18 圓盤A 自由膜表面更新Fig.18 Free liquid film renewal frequency of disk A

    圖19 半徑r=45 mm 時更新頻率隨角度的變化規(guī)律Fig.19 Relationship between liquid surface renewal frequency and angle (r=45 mm)

    圖20 θ=π 時液膜速度變化對比Fig.20 Comparison of liquid film velocity (θ=π)

    圖21 圓盤A 膜厚分布(z=0.061 m 為圓盤壁面)Fig.21 Distribution of disc A film thickness

    分析認(rèn)為,開窗區(qū)域沒有壁面的影響,使得重力作用效果明顯,液膜向下加速,重力方向(-y方 向)速度加快,如圖20所示(θ=π,開窗區(qū)域r=30~50 mm),液膜從窗口上緣r=50 mm 處開始加速,到達(dá)窗口下緣r=30 mm 處液膜積累,出現(xiàn)減速,到達(dá)壁面后又開始加速。根據(jù)質(zhì)量守恒定律,液膜加速使膜厚變薄,如圖21所示,開窗區(qū)域液膜厚度顯著減小,內(nèi)層液體暴露出來,更新加快。

    4.3.2 開窗形狀與液膜更新的關(guān)系 圓盤C 表面更新分布如圖22所示。計(jì)算發(fā)現(xiàn),開窗區(qū)域更新頻率較圓盤A 更快,圓盤C 開窗區(qū)域自由膜平均更新頻率SCF=3.26 s-1,高出平均值36%,高出圓盤A 自由膜平均更新頻率24%,SAF=2.63 s-1。但是圓盤C附壁膜平均更新頻率SCW=1.99 s-1,比圓盤A 低13%,SAW=2.28 s-1。最終圓盤C 液膜表面平均更新頻率為SC=2.37 s-1,接近圓盤A,SA=2.35 s-1。

    圖22 圓盤C 表面更新頻率分布Fig.22 Distribution of disk C surface renewal frequency

    圖23 半徑r=45mm 時更新頻率隨角度的變化規(guī)律Fig.23 Relationship between surface renewal and angle

    圓盤A、B、C 半徑r=45 mm 時液膜更新頻率隨角度的變化規(guī)律如圖23所示。圓盤C 開窗區(qū)域自由膜的更新頻率比圓盤A 快,但是開窗區(qū)域壁面附近附壁膜的更新頻率沒有因?yàn)殚_窗區(qū)域的存在而明顯增加,而是接近實(shí)心圓盤B,這點(diǎn)不同于圓盤A。因此圓盤C 附壁膜更新頻率小于圓盤A,最終的結(jié)果是兩者平均更新頻率相差不大??梢酝茰y,當(dāng)增加圓盤開孔率后,開圓形窗的優(yōu)勢會顯現(xiàn)出來。

    4.3.3 開窗大小與液膜更新的關(guān)系 對比開窗圓盤C 和D,圓盤D 中每個小圓窗的面積為圓盤C 中大圓窗的1/2,兩者總的開窗面積相同。計(jì)算發(fā)現(xiàn),自由膜平均更新頻率SDF=3.28 s-1,與圓盤C 接近;附壁膜更新頻率SDW=2.82 s-1,比圓盤C 高42%??傮w而言,圓盤D 液膜平均更新頻率SD=2.91 s-1,比圓盤B 高73%,比圓盤A 和C 高22%。

    在計(jì)算范圍內(nèi)結(jié)果顯示,開窗大小對于自由膜更新頻率影響不顯著(單窗面積不宜過?。?,增加整體開孔率作用更突出。適當(dāng)減小開窗區(qū)域間距,可強(qiáng)化附壁膜的更新頻率。

    4.3.4 4 種圓盤液膜表面更新對比 將4 種圓盤自由膜與附壁膜和液膜平均更新頻率分別作圖比較,如圖24所示。

    通過比較,可以得到以下規(guī)律。

    (1)圓盤A 與B 比較,圓盤A 液膜平均更新頻率、自由膜和附壁膜更新頻率比圓盤B 附壁膜分別高40%、57%和36%。開窗區(qū)域自由膜更新頻率顯著提高。

    (2)圓盤A 與C 比較,圓盤C 的自由膜更新頻率比A 高24%,附壁膜更新頻率比A 低13%,平均更新頻率相當(dāng)。圓形窗口形成的自由膜比扇形窗自由膜更新頻率快,但是對附壁膜影響較小。

    圖24 圓盤結(jié)構(gòu)對液膜更新的影響Fig.24 Relationship between disc structures and liquid film surface renewal frequency

    (3) 圓盤C 與D 比較,兩者自由膜更新頻率 相當(dāng),圓盤D 附壁膜更新頻率比C 高42%,因此平均更新頻率比C 高23%。適當(dāng)減小開窗區(qū)域的間距能夠有效強(qiáng)化自由膜對附壁膜的影響。

    5 結(jié) 論

    表面更新是探討相際傳質(zhì)的基本屬性,本研究建議更新頻率與變形率相關(guān)聯(lián),通過數(shù)值模擬探討圓盤液膜更新,得到以下結(jié)論。

    (1)普通圓盤液膜局部更新頻率沿盤面呈一定的分布;隨半徑增加到一定值后不再變化。隨角度先增加后減小,最大值出現(xiàn)在加速區(qū)。更新頻率較其他區(qū)域高150%,較平均值高75%。

    (2)普通圓盤液膜平均更新頻率,隨轉(zhuǎn)速S∝ω1.2變化,隨半徑S∝R0.7變化。

    (3)開窗圓盤表面存在自由膜和附壁膜,這是兩種不同性質(zhì)的液膜。開窗區(qū)域自由膜厚度較薄,更新較快,同時能激化周圍壁面液膜更新,兩種液膜更新頻率的量化比較揭示了開窗圓盤強(qiáng)化機(jī)理,也證明開窗是一種強(qiáng)化措施。數(shù)值模擬可以對結(jié)構(gòu)化圓盤做優(yōu)化設(shè)計(jì)。

    符 號 說 明

    A——?dú)庖航佑|面積,m2

    Ca——Capillary 數(shù),Ca=μω R/σ

    Fr——Froude 數(shù),F(xiàn)r=ω2R/g

    g——重力加速度,m·s-2

    h——液膜厚度,mm

    l——旋渦長度尺寸,m

    N——圓盤轉(zhuǎn)速,r·min-1

    R——圓盤半徑,mm

    Re——Reynolds 數(shù),Re=

    r——某點(diǎn)位置半徑,mm

    SAW——圓盤A 表面附壁膜平均更新頻率,s-1

    SA0——圓盤A 表面自由膜平均更新頻率,s-1

    SP——液膜局部更新頻率,s-1

    S0——液膜平均更新頻率,s-1

    t——時間,s

    ue——旋渦特征速度,m·s-1

    α——相體積分率

    θ——角度,rad

    ρ——密度,kg·m-3

    σ——表面張力,N·m

    下角標(biāo)

    A,B,C,D——圓盤

    F——自由膜

    W——附壁膜

    0——平均值

    [1]Woo B,Choi K Y,Song K H.Melt polycondensation of bisphenol a polycarbonate by forced gas sweeping process (Ⅱ):Continuous rotating-disk reactor [J].Industrial & Engineering Chemistry Research,2001,40 (16):3459-3466

    [2]Sakurai A,Imai H,Takenaka Y,et al.Simulation of citric acid production by rotating disk contactor [J].Biotechnol.Bioeng.,1997,56 (6):689-696

    [3]Li K,Yang C,Xu Y,et al.Effect of inorganic anions on Rhodamine B removal under visible light irradiation using Bi2O3/Ti rotating disk reactor [J].Chemical Engineering Journal,2012,211/212:208-215

    [4]Vasiliadou I A,Molina R,Martínez F,et al.Experimental and modeling study on removal of pharmaceutically active compounds in rotating biological contactors [J].Journal of Hazardous Materials,2014,274:473-482

    [5]Afanasiev K,Münch A,Wagner B.Thin film dynamics on a vertically rotating disk partially immersed in a liquid bath [J].Applied Mathematical Modelling,2008,32 (9):1894-1911

    [6]Suor D,Ma J,Wang Z,et al.Enhanced power production from waste activated sludge in rotating-cathode microbial fuel cells:the effects of aquatic worm predation [J].Chemical Engineering Journal,2014,248:415-421

    [7]Bintanja H H J,van der Erve J J V M,Boelhouwer C.Oxygen transfer in a rotating disc treatment plant [J].Water Research,1975,9 (12):1147-1153

    [8]Zeevalkink J A,Kelderman P,Boelhouwer C.Liquid film thickness in a rotating disc gas-liquid contactor [J].Water Research,1978,12 (8):577-581

    [9]Zeevalkink J A,Kelderman P,Visser D C,et al.Physical mass transfer in a rotating disc gas-liquid contactor [J].Water Research,1979,13 (9):913-919

    [10]Kubsad V,Chaudhari S,Gupta S K.Model for oxygen transfer in rotating biological contactor [J].Water Research,2004,38 (20):4297-4304

    [11]Wang Liangsheng (王良生).The liquid moment and mass transfer on the rotating disc-ring [D].Shanghai:East China University of Science and Technology,1999

    [12]Zhou Xianjue (周賢爵).Flow behavior of high viscosity fluid in a cage-like reactor [D].Shanghai:East China University of Science and Technology,2001

    [13]Danckwerts P V.Significance of liquid-film coefficients in gas absorption [J].Industrial & Engineering Chemistry,1951,43 (6):1460-1467

    [14]Perlmutter D D.Surface-renewal models in mass transfer [J].Chemical Engineering Science,1961,16 (3/4):287-296

    [15]Harriott P.A random eddy modification of the penetration theory [J].Chemical Engineering Science,1962,17 (3):149-154

    [16]Koppel L B P R.Statistical models for surface renewal in heat and mass transfer (Ⅳ):Wall to fluidized bed heat transfer coefficients [J].AIChE J.,1970,16 (3):464-471

    [17]Benjamin T F,Chung L T F A.General mathematical models of transport processes with and without chemical reactions [J].Canadian Journal of Chemical Engineering,1971,3 (49):340-345

    [18]Murakami Y,Fujimoto K,Kakimoto S,et al.On a high viscosity polymer finisher apparatus with two agitator axes having multidisks [J].Journal of Chemical Engineering of Japan,1972,5 (3):257-263

    [19]Biesenberger J A,Sebastian D H.Principles of Polymerization Engineering [M].New York:Wiley,1983

    [20]Pan Qinmin ( 潘 勤 敏).The gas-liquid mass transfer and devolatilization during the process of polymerization [D].Hangzhou:Zhejiang University,1987

    [21]Dai Gance (戴干策),Zhou Xianjue (周賢爵).The liquid film moment on rotating disc-ring PET finishers [J].Chemical Engineering(China) (化學(xué)工程),2002,30 (S):211-216

    [22]Banerjee S,Scott D S,Rhodes E.Mass transfer to falling wavy liquid films in turbulent flow [J].Industrial & Engineering Chemistry Fundamentals,1968,7 (1):22-27

    [23]Jajuee B,Margaritis A,Karamanev D,et al.Application of surface-renewal-stretch model for interface mass transfer [J].ChemicalEngineering Science,2006,61 (12):3917-3929

    [24]Hirt C W,Nichols B D.Volume of fluid (VOF) method for the dynamics of free boundaries [J].Journal of Computational Physics,1981,39 (1):201-225

    [25]Stewartson K.On the flow between two rotating coaxial disks [J].Mathematical Proceedings of the Cambridge Philosophical Society,1953,49 (2):333-341

    [26]James O Wilkes.Fluid Mechanics for Chemical Engineers [M].USA:Prentice Hall PTR,1999:280-285

    猜你喜歡
    附壁液膜圓盤
    純磨玻璃結(jié)節(jié)的HRCT征象對非附壁型浸潤性肺腺癌的診斷價值
    考慮軸彎曲的水潤滑軸承液膜建模方法
    高空高速氣流下平板液膜流動與破裂規(guī)律
    液膜破裂對PCCS降膜的影響*
    圓盤鋸刀頭的一種改進(jìn)工藝
    石材(2020年6期)2020-08-24 08:27:00
    掘進(jìn)工作面附壁風(fēng)筒控塵特性研究
    單位圓盤上全純映照模的精細(xì)Schwarz引理
    奇怪的大圓盤
    基于Profibus-DP的圓盤澆鑄控制系統(tǒng)的應(yīng)用
    卵巢交界性黏液性腫瘤伴惡變及附壁肉瘤結(jié)節(jié)1 例及文獻(xiàn)復(fù)習(xí)
    日本-黄色视频高清免费观看| 26uuu在线亚洲综合色| 国产精品人妻久久久影院| 亚洲国产精品一区三区| 麻豆精品久久久久久蜜桃| 丝袜喷水一区| 大香蕉久久网| 欧美人与性动交α欧美精品济南到 | 久久国产精品大桥未久av| 一级,二级,三级黄色视频| av又黄又爽大尺度在线免费看| 毛片一级片免费看久久久久| 亚洲精品在线美女| 国产伦理片在线播放av一区| av在线观看视频网站免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费视频播放在线视频| 久久精品国产亚洲av天美| 1024香蕉在线观看| 黄片无遮挡物在线观看| 欧美日韩视频高清一区二区三区二| 水蜜桃什么品种好| 亚洲国产看品久久| 制服丝袜香蕉在线| 各种免费的搞黄视频| 9191精品国产免费久久| 亚洲伊人色综图| 在线天堂中文资源库| av在线观看视频网站免费| 男女国产视频网站| 亚洲精品在线美女| 伦理电影大哥的女人| 两性夫妻黄色片| xxx大片免费视频| 国产一区二区 视频在线| 我要看黄色一级片免费的| 99热全是精品| 伊人久久国产一区二区| 精品亚洲乱码少妇综合久久| 亚洲国产毛片av蜜桃av| 秋霞伦理黄片| 欧美人与善性xxx| 久久久国产精品麻豆| 老鸭窝网址在线观看| 成年女人毛片免费观看观看9 | 极品人妻少妇av视频| 搡老乐熟女国产| 97人妻天天添夜夜摸| 日韩电影二区| freevideosex欧美| av天堂久久9| 一级a爱视频在线免费观看| 久久狼人影院| 精品国产露脸久久av麻豆| 欧美人与善性xxx| 久久这里有精品视频免费| 两个人免费观看高清视频| 日韩电影二区| 九草在线视频观看| 成年美女黄网站色视频大全免费| 满18在线观看网站| 狠狠婷婷综合久久久久久88av| 国产精品一区二区在线观看99| 肉色欧美久久久久久久蜜桃| 久久久久久久久久久免费av| 女人久久www免费人成看片| 欧美激情 高清一区二区三区| 精品亚洲乱码少妇综合久久| 欧美av亚洲av综合av国产av | 九色亚洲精品在线播放| 人妻人人澡人人爽人人| 一本大道久久a久久精品| 日韩欧美精品免费久久| 狂野欧美激情性bbbbbb| a级毛片黄视频| 狠狠精品人妻久久久久久综合| 精品国产一区二区三区四区第35| 在线天堂中文资源库| 男女边摸边吃奶| 黑人猛操日本美女一级片| 亚洲欧美精品综合一区二区三区 | 色吧在线观看| 丝袜美足系列| 欧美在线黄色| 国产女主播在线喷水免费视频网站| av又黄又爽大尺度在线免费看| 校园人妻丝袜中文字幕| 咕卡用的链子| 最近2019中文字幕mv第一页| 老汉色∧v一级毛片| 老女人水多毛片| 成人午夜精彩视频在线观看| 日本黄色日本黄色录像| 色网站视频免费| av.在线天堂| 狠狠婷婷综合久久久久久88av| 男女边摸边吃奶| 国产探花极品一区二区| 欧美黄色片欧美黄色片| 赤兔流量卡办理| 蜜桃在线观看..| 国产又色又爽无遮挡免| 大香蕉久久网| 久久久久久久久久久久大奶| 两性夫妻黄色片| 日韩三级伦理在线观看| av在线播放精品| 男人舔女人的私密视频| 亚洲综合色惰| 国产免费一区二区三区四区乱码| 人妻 亚洲 视频| 高清黄色对白视频在线免费看| av福利片在线| 男女高潮啪啪啪动态图| 成人18禁高潮啪啪吃奶动态图| 综合色丁香网| av卡一久久| 日韩一区二区三区影片| 女人被躁到高潮嗷嗷叫费观| 男女午夜视频在线观看| 精品视频人人做人人爽| 美女中出高潮动态图| 欧美激情 高清一区二区三区| 欧美激情极品国产一区二区三区| 国产熟女午夜一区二区三区| 亚洲欧美清纯卡通| 午夜免费男女啪啪视频观看| 国产极品粉嫩免费观看在线| 亚洲国产色片| kizo精华| 最新中文字幕久久久久| 春色校园在线视频观看| 欧美黄色片欧美黄色片| 18禁裸乳无遮挡动漫免费视频| 曰老女人黄片| 日韩,欧美,国产一区二区三区| 视频区图区小说| 亚洲成人手机| 久久狼人影院| 免费在线观看黄色视频的| 国产在视频线精品| 尾随美女入室| 中国三级夫妇交换| 久久精品夜色国产| 国产白丝娇喘喷水9色精品| 桃花免费在线播放| 九色亚洲精品在线播放| 色94色欧美一区二区| 爱豆传媒免费全集在线观看| 午夜激情av网站| 亚洲久久久国产精品| 视频区图区小说| 久久国内精品自在自线图片| 亚洲男人天堂网一区| 最新中文字幕久久久久| 一二三四在线观看免费中文在| 亚洲男人天堂网一区| 亚洲美女视频黄频| 97精品久久久久久久久久精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 啦啦啦视频在线资源免费观看| 国产乱人偷精品视频| 成人毛片60女人毛片免费| 亚洲av电影在线观看一区二区三区| 精品少妇一区二区三区视频日本电影 | 亚洲成av片中文字幕在线观看 | 一级毛片黄色毛片免费观看视频| 九草在线视频观看| 欧美日韩精品成人综合77777| 伦理电影免费视频| 一级a爱视频在线免费观看| 日韩制服丝袜自拍偷拍| 日韩在线高清观看一区二区三区| 波野结衣二区三区在线| 两性夫妻黄色片| 久久ye,这里只有精品| 热re99久久精品国产66热6| 女性生殖器流出的白浆| 伦理电影免费视频| 91成人精品电影| 交换朋友夫妻互换小说| 免费观看无遮挡的男女| 日韩欧美一区视频在线观看| freevideosex欧美| 日韩欧美精品免费久久| 久久精品国产亚洲av涩爱| 国产一区二区在线观看av| 美女国产高潮福利片在线看| 亚洲成色77777| 国产精品三级大全| 久热久热在线精品观看| 亚洲成人av在线免费| 熟女少妇亚洲综合色aaa.| 国产精品免费视频内射| www.av在线官网国产| 一级毛片 在线播放| 国产精品99久久99久久久不卡 | 中文欧美无线码| 欧美日韩国产mv在线观看视频| 午夜福利乱码中文字幕| av片东京热男人的天堂| 久久久精品区二区三区| 大话2 男鬼变身卡| 亚洲婷婷狠狠爱综合网| 18+在线观看网站| 精品一区二区三区四区五区乱码 | 国产片内射在线| 国产爽快片一区二区三区| 亚洲av日韩在线播放| 亚洲精品国产色婷婷电影| 自线自在国产av| 伊人亚洲综合成人网| 国精品久久久久久国模美| 伊人久久国产一区二区| 亚洲精品,欧美精品| 色播在线永久视频| 热99久久久久精品小说推荐| 国产一区二区在线观看av| 最近2019中文字幕mv第一页| 色94色欧美一区二区| 又粗又硬又长又爽又黄的视频| 国产精品一二三区在线看| 国产成人精品在线电影| 欧美精品av麻豆av| 国产成人精品一,二区| 久久99精品国语久久久| 99久国产av精品国产电影| 国产精品二区激情视频| 欧美少妇被猛烈插入视频| 久久久久久人妻| 日韩一本色道免费dvd| av卡一久久| 少妇被粗大的猛进出69影院| 成人国产av品久久久| 午夜福利视频在线观看免费| 亚洲av日韩在线播放| 亚洲欧洲日产国产| 在线 av 中文字幕| 美国免费a级毛片| 国产成人av激情在线播放| 男的添女的下面高潮视频| 一本色道久久久久久精品综合| 久久久久精品性色| 青春草国产在线视频| 黄色配什么色好看| 777米奇影视久久| www日本在线高清视频| 欧美97在线视频| 精品国产露脸久久av麻豆| 精品久久蜜臀av无| 亚洲五月色婷婷综合| 久久久久国产网址| 性少妇av在线| 又黄又粗又硬又大视频| 两个人免费观看高清视频| 国产精品一国产av| av卡一久久| 国产免费又黄又爽又色| 2022亚洲国产成人精品| 国产精品欧美亚洲77777| 欧美日韩av久久| 精品人妻在线不人妻| 满18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 日本wwww免费看| 欧美 亚洲 国产 日韩一| 国产精品二区激情视频| 黑人猛操日本美女一级片| 亚洲伊人色综图| 街头女战士在线观看网站| 免费在线观看视频国产中文字幕亚洲 | 亚洲综合色网址| 女的被弄到高潮叫床怎么办| 亚洲精品久久久久久婷婷小说| 亚洲欧美清纯卡通| 国产高清不卡午夜福利| 91在线精品国自产拍蜜月| 少妇被粗大的猛进出69影院| 成年女人毛片免费观看观看9 | 精品少妇内射三级| 久久精品aⅴ一区二区三区四区 | 日韩欧美一区视频在线观看| 宅男免费午夜| 欧美日韩成人在线一区二区| 国产高清国产精品国产三级| 最近的中文字幕免费完整| 精品人妻一区二区三区麻豆| 精品少妇黑人巨大在线播放| 国产精品国产三级专区第一集| 国产精品 欧美亚洲| 久久久久精品性色| 成年女人在线观看亚洲视频| 中文字幕人妻丝袜制服| 多毛熟女@视频| 激情视频va一区二区三区| 在线精品无人区一区二区三| 亚洲熟女精品中文字幕| 亚洲色图综合在线观看| 国产亚洲精品第一综合不卡| 日本黄色日本黄色录像| 国产乱人偷精品视频| 久久女婷五月综合色啪小说| 日韩精品有码人妻一区| 高清欧美精品videossex| 国产成人午夜福利电影在线观看| 男人操女人黄网站| 亚洲婷婷狠狠爱综合网| 国产精品香港三级国产av潘金莲 | 国产老妇伦熟女老妇高清| 尾随美女入室| 国产麻豆69| 亚洲欧美中文字幕日韩二区| 伦精品一区二区三区| 久久久精品国产亚洲av高清涩受| 老司机亚洲免费影院| 欧美老熟妇乱子伦牲交| 中文欧美无线码| 午夜91福利影院| 久久亚洲国产成人精品v| 少妇人妻久久综合中文| 欧美 亚洲 国产 日韩一| 少妇被粗大的猛进出69影院| 日韩制服骚丝袜av| 丝瓜视频免费看黄片| 成人黄色视频免费在线看| 黄网站色视频无遮挡免费观看| 日韩人妻精品一区2区三区| 亚洲av免费高清在线观看| 欧美变态另类bdsm刘玥| 少妇被粗大猛烈的视频| 一区二区三区乱码不卡18| 秋霞伦理黄片| 最新中文字幕久久久久| 伦理电影免费视频| 波多野结衣一区麻豆| 亚洲欧洲精品一区二区精品久久久 | 午夜福利乱码中文字幕| 免费女性裸体啪啪无遮挡网站| 日韩精品有码人妻一区| 人人妻人人添人人爽欧美一区卜| 日本爱情动作片www.在线观看| 人妻人人澡人人爽人人| 99久久人妻综合| 亚洲av在线观看美女高潮| 国产爽快片一区二区三区| 99久久人妻综合| 国产成人aa在线观看| 午夜影院在线不卡| 国产免费福利视频在线观看| 欧美日韩视频精品一区| 超色免费av| 日本-黄色视频高清免费观看| 午夜av观看不卡| 天堂8中文在线网| 丝瓜视频免费看黄片| 久久国产精品大桥未久av| 99香蕉大伊视频| 99久久人妻综合| 国产不卡av网站在线观看| 伊人亚洲综合成人网| 熟女少妇亚洲综合色aaa.| 好男人视频免费观看在线| 涩涩av久久男人的天堂| 久久国产精品大桥未久av| 亚洲成国产人片在线观看| www日本在线高清视频| 亚洲久久久国产精品| 成人黄色视频免费在线看| 国精品久久久久久国模美| 国产精品国产三级国产专区5o| 考比视频在线观看| 国产精品无大码| 午夜免费鲁丝| 国产无遮挡羞羞视频在线观看| 老司机影院毛片| 亚洲av成人精品一二三区| 久久国产亚洲av麻豆专区| 中文字幕制服av| 亚洲中文av在线| 七月丁香在线播放| 国语对白做爰xxxⅹ性视频网站| 男女午夜视频在线观看| 国产有黄有色有爽视频| 免费在线观看黄色视频的| 热99国产精品久久久久久7| 午夜福利视频在线观看免费| 国产有黄有色有爽视频| 五月伊人婷婷丁香| 欧美 亚洲 国产 日韩一| 国产日韩欧美在线精品| 女人久久www免费人成看片| 视频区图区小说| 中文天堂在线官网| 女人精品久久久久毛片| 国产在视频线精品| 欧美av亚洲av综合av国产av | 丝袜脚勾引网站| 夫妻性生交免费视频一级片| 亚洲成国产人片在线观看| 日韩免费高清中文字幕av| 免费在线观看完整版高清| 亚洲伊人久久精品综合| 免费看av在线观看网站| 亚洲精华国产精华液的使用体验| 午夜av观看不卡| 一级黄片播放器| 国产探花极品一区二区| 伊人亚洲综合成人网| 极品人妻少妇av视频| 精品少妇久久久久久888优播| 中文乱码字字幕精品一区二区三区| 18禁动态无遮挡网站| 97在线人人人人妻| 人体艺术视频欧美日本| 国产黄频视频在线观看| 少妇人妻久久综合中文| 9色porny在线观看| 少妇熟女欧美另类| 考比视频在线观看| 亚洲一码二码三码区别大吗| av福利片在线| 少妇猛男粗大的猛烈进出视频| 一边摸一边做爽爽视频免费| 亚洲精品一区蜜桃| 熟女电影av网| 中文字幕人妻熟女乱码| 我的亚洲天堂| 久久精品夜色国产| 天天躁狠狠躁夜夜躁狠狠躁| 母亲3免费完整高清在线观看 | 永久网站在线| 国产极品粉嫩免费观看在线| 日本爱情动作片www.在线观看| 人妻少妇偷人精品九色| 久久久久网色| av在线老鸭窝| 国产一区二区激情短视频 | 夫妻性生交免费视频一级片| 9色porny在线观看| 婷婷成人精品国产| 免费在线观看完整版高清| 一本—道久久a久久精品蜜桃钙片| 午夜激情av网站| 在线精品无人区一区二区三| 国产精品免费大片| 999精品在线视频| 91午夜精品亚洲一区二区三区| 日日摸夜夜添夜夜爱| 久久久精品免费免费高清| 亚洲激情五月婷婷啪啪| 看免费av毛片| 国产精品国产av在线观看| 日韩不卡一区二区三区视频在线| 91成人精品电影| 中文字幕人妻熟女乱码| 女人高潮潮喷娇喘18禁视频| 日韩制服骚丝袜av| 男女无遮挡免费网站观看| 大码成人一级视频| 欧美老熟妇乱子伦牲交| 丰满乱子伦码专区| 国产又色又爽无遮挡免| 精品亚洲成a人片在线观看| 香蕉国产在线看| 在线观看www视频免费| 精品国产乱码久久久久久小说| 亚洲国产色片| 亚洲色图 男人天堂 中文字幕| 亚洲精品第二区| a级毛片在线看网站| 国产无遮挡羞羞视频在线观看| 免费高清在线观看视频在线观看| 欧美成人精品欧美一级黄| 91aial.com中文字幕在线观看| 伊人久久大香线蕉亚洲五| 五月伊人婷婷丁香| 欧美 日韩 精品 国产| 99久久中文字幕三级久久日本| 99久久人妻综合| 最新中文字幕久久久久| av在线老鸭窝| 大码成人一级视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄片播放在线免费| 日韩在线高清观看一区二区三区| 美女国产视频在线观看| 伊人久久大香线蕉亚洲五| 国产熟女午夜一区二区三区| 色婷婷久久久亚洲欧美| 电影成人av| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 大码成人一级视频| 亚洲精品久久成人aⅴ小说| 9色porny在线观看| 国产毛片在线视频| 欧美日韩国产mv在线观看视频| av网站免费在线观看视频| 日韩 亚洲 欧美在线| 日本av免费视频播放| 美女视频免费永久观看网站| 各种免费的搞黄视频| 欧美老熟妇乱子伦牲交| 校园人妻丝袜中文字幕| 日本-黄色视频高清免费观看| 女性被躁到高潮视频| 欧美精品一区二区免费开放| 免费看av在线观看网站| 国产97色在线日韩免费| 午夜福利视频在线观看免费| 国产精品一区二区在线不卡| 欧美老熟妇乱子伦牲交| 欧美另类一区| 人人妻人人澡人人爽人人夜夜| 天堂8中文在线网| 人体艺术视频欧美日本| 最近手机中文字幕大全| 精品福利永久在线观看| 国产成人a∨麻豆精品| 成人午夜精彩视频在线观看| 日韩三级伦理在线观看| 亚洲精品国产av蜜桃| 男女边摸边吃奶| 久久国产亚洲av麻豆专区| 蜜桃在线观看..| av免费在线看不卡| 亚洲精品在线美女| 国产av精品麻豆| 91久久精品国产一区二区三区| 亚洲三区欧美一区| 熟女av电影| 婷婷色综合大香蕉| 在线 av 中文字幕| 国产精品二区激情视频| 国产成人91sexporn| 精品久久久精品久久久| 高清在线视频一区二区三区| 国产av码专区亚洲av| 制服人妻中文乱码| 国产一级毛片在线| 亚洲中文av在线| 日本vs欧美在线观看视频| 一级黄片播放器| 国产成人午夜福利电影在线观看| 成人手机av| 一区二区三区四区激情视频| 久久精品国产鲁丝片午夜精品| 黑丝袜美女国产一区| 激情五月婷婷亚洲| 黑人巨大精品欧美一区二区蜜桃| 99国产综合亚洲精品| 伦理电影免费视频| 久久国产精品男人的天堂亚洲| 精品久久蜜臀av无| 国产av一区二区精品久久| 高清av免费在线| 国产成人午夜福利电影在线观看| av福利片在线| 99热全是精品| 亚洲精品日韩在线中文字幕| 91精品伊人久久大香线蕉| 不卡视频在线观看欧美| 成人18禁高潮啪啪吃奶动态图| 亚洲精华国产精华液的使用体验| 99热网站在线观看| 久久国产精品男人的天堂亚洲| 美女国产视频在线观看| 欧美97在线视频| av网站在线播放免费| 精品人妻在线不人妻| 色哟哟·www| 熟女av电影| 中文精品一卡2卡3卡4更新| 精品亚洲成a人片在线观看| 香蕉国产在线看| 99热全是精品| 最新中文字幕久久久久| 国产爽快片一区二区三区| 亚洲综合色惰| 99久久人妻综合| 啦啦啦在线观看免费高清www| av在线老鸭窝| 久久久精品94久久精品| 一区二区日韩欧美中文字幕| 国产高清国产精品国产三级| 日韩在线高清观看一区二区三区| 天美传媒精品一区二区| 在线观看人妻少妇| 欧美日韩视频精品一区| 亚洲一区中文字幕在线| 亚洲一码二码三码区别大吗| 一区二区三区四区激情视频| 两个人免费观看高清视频| av片东京热男人的天堂| 精品亚洲成国产av| 日韩制服骚丝袜av| 精品人妻一区二区三区麻豆| 精品久久久久久电影网| 婷婷成人精品国产| 国产成人精品久久二区二区91 | 五月天丁香电影| 国产又色又爽无遮挡免| 日韩精品有码人妻一区| 亚洲四区av| 人妻 亚洲 视频| 五月伊人婷婷丁香| 曰老女人黄片| 一级毛片黄色毛片免费观看视频| www.精华液| 伊人久久大香线蕉亚洲五| 成人亚洲欧美一区二区av| av一本久久久久| 国产在线一区二区三区精| 亚洲一码二码三码区别大吗| 黑人欧美特级aaaaaa片|