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

    計算流體力學(xué)在水處理膜過程中的應(yīng)用

    2014-09-25 03:44:46張雅琴侯立安
    中國工程科學(xué) 2014年7期
    關(guān)鍵詞:膜分離傳質(zhì)極化

    張雅琴 ,張 林 ,侯立安 ,2

    (1.浙江大學(xué)化學(xué)工程與生物工程學(xué)系,杭州 310027;2.第二炮兵后勤科學(xué)技術(shù)研究所,北京 100011)

    計算流體力學(xué)在水處理膜過程中的應(yīng)用

    張雅琴1,張 林1,侯立安1,2

    (1.浙江大學(xué)化學(xué)工程與生物工程學(xué)系,杭州 310027;2.第二炮兵后勤科學(xué)技術(shù)研究所,北京 100011)

    本文介紹了計算流體力學(xué)(CFD)在膜分離過程模擬中的基本原理,對CFD在壓力膜驅(qū)動過程組件設(shè)計、傳質(zhì)和膜污染過程研究方面的應(yīng)用和CFD在具有相變膜過程(滲透汽化、膜蒸餾)以及其他膜過程中的應(yīng)用分別進行了綜述,最后對CFD在膜分離技術(shù)研究中的應(yīng)用前景進行了展望。

    計算流體力學(xué);膜分離;壓力驅(qū)動;滲透汽化;膜蒸餾

    1 前言

    計算流體力學(xué)(CFD)是以Navier-Stokes方程為基礎(chǔ),把原來在時間和空間上連續(xù)的物理量用一系列有限個離散點上的變量集合來替代,得到流場內(nèi)各個位置上的基本物理量(如速度、壓力、溫度、濃度等)的分布,以及這些物理量隨時間的變化情況。CFD技術(shù)涉及計算機、物理、流體力學(xué)、數(shù)值計算、可視化技術(shù)等多學(xué)科,這些學(xué)科的知識綜合起來,可以提供建立流體流動模型的方式和方法,其計算方法主要包括有限元法、有限體積法、有限差分法[1]。與實驗手段相比,CFD技術(shù)具有精確、效率高、成本低,可以模擬真實及理想條,后處理技術(shù)較完善,便于分析計算結(jié)果等優(yōu)點。

    膜技術(shù)被稱為“21世紀的水處理技術(shù)”,在水處理領(lǐng)域得到了廣泛的應(yīng)用。然而隨著對膜技術(shù)研究的逐步深入,人們發(fā)現(xiàn)由于受物理模型實驗條件的限制,一些膜組件結(jié)構(gòu)設(shè)計問題往往只能夠定性或者半定量進行,深入至機理層面的系統(tǒng)而量化的研究受到阻礙;此外,對于影響膜分離效率和穩(wěn)定性的傳質(zhì)、傳熱機理以及膜污染機理的研究也受限于實驗條件。CFD技術(shù)可以很好地克服這些問題,得到的結(jié)果定量化準確性高,可為實驗設(shè)計和組件優(yōu)化提供科學(xué)依據(jù),對于組件內(nèi)部的流動狀態(tài)可以進行可視化模擬。因此,越來越多的膜技術(shù)領(lǐng)域研究者嘗試將CFD技術(shù)引入到膜技術(shù)的研究中[2,3]。

    2 模擬原理

    常規(guī)流體過程的CFD是在質(zhì)量、動量和能量三大守恒定律的基礎(chǔ)上,通過解流體力學(xué)三大方程得到對變量的表述。而膜分離過程除了有流體流動還涉及到溶質(zhì)和溶劑分離的問題,因此CFD應(yīng)用于膜分離技術(shù)需要在三大守恒方程的基礎(chǔ)上,同時考慮溶質(zhì)傳遞方程。若要準確模擬膜分離過程中流體的水力學(xué)狀況,就要針對每個具體模型的特殊性,建立起相符的邊界條件和初始條件以及可以準確描述膜上下游界面的傳質(zhì)、傳熱模型。在模擬膜分離過程中,常用的邊界條件如表1所示[4]。對于存在隔網(wǎng)的膜組件,由于隔網(wǎng)單元具有周期性,會導(dǎo)致膜組件內(nèi)部流場的周期性分布,故一些研究者為了減小計算量采用了周期性邊界條件,大大減小了計算強度。此外,關(guān)于膜邊界條件的描述也不斷完善,從早期將膜視為不可滲透壁面,到后來將膜作為多孔壁面進行處理;也有一些學(xué)者將膜過程的傳質(zhì)、傳熱模型與CFD進行耦合作為膜的邊界條件。

    表1 膜分離過程中常用的邊界條件Table 1 Typical boundary condition in membrane process

    3 CFD在膜技術(shù)中的應(yīng)用現(xiàn)狀

    在過去的數(shù)年里,CFD在水處理膜過程領(lǐng)域應(yīng)用的研究主要集中在壓力驅(qū)動膜過程(超濾、納濾、反滲透)的組件優(yōu)化設(shè)計以及流場運動情況(濃差極化和膜污染)方面。最近,將CFD技術(shù)用于具有相變過程的膜分離技術(shù)傳質(zhì)、傳熱過程以及其他新型膜過程(生物膜反應(yīng)器、正滲透等)的研究報道也逐漸多起來。

    3.1 壓力驅(qū)動膜過程

    3.1.1 CFD在膜污染機理研究中的應(yīng)用

    膜污染是造成壓力驅(qū)動膜組件運行失常和限制其進一步廣泛應(yīng)用的主要因素。廣義膜污染不僅包括膜孔阻塞引起的不可逆污染,還包括由濃差極化和凝膠層形成的可逆污染。當(dāng)截留物(主要包括無機鹽、有機物、膠體以及細菌等)在膜表面沉淀與積累,使水透過膜的阻力增加,妨礙了膜面上的溶質(zhì)擴散,從而導(dǎo)致膜產(chǎn)生量和水質(zhì)下降。在膜污染研究中,CFD研究結(jié)果具有精確、量化和更具說服力的特點。

    CFD技術(shù)應(yīng)用于膜污染機理研究時,需要將CFD模型與膜污染模型耦合求解。早期的研究將膜視為簡單的平面,研究流場以及傳質(zhì)對膜污染形成的影響。Kapellos[5]、Schulenburg[6]等在研究中將膜視為多孔介質(zhì),模擬了膜污染的形成。

    然而,將滲透膜假設(shè)為不可滲透的膜平面或者處理為多空介質(zhì)與真實膜性能存在很大的差異性,因此,后來有研究者將膜的傳質(zhì)模型與CFD耦合,用來模擬真實的膜污染和傳質(zhì)過程。其中,Pinho[7]等將納濾傳質(zhì)模型與CFD模型進行耦合,并且在納濾傳質(zhì)模型中考慮溶質(zhì)—溶劑—膜之間的相互作用力,通過CFD模擬的邊界層定義,同時將滲透通量實驗數(shù)據(jù)進行二維模擬獲得截留率和流體濃度分布狀況,模擬結(jié)果與實驗結(jié)果吻合良好。

    Ahmad[8,9]等利用薄膜理論定義了滲透膜邊界條件,利用非平衡熱力學(xué)傳質(zhì)模型定義了滲透膜通量,對卷式納濾膜滲透過程進行了二維模擬。模擬中他們認為上壁面是無滑移的不可滲透壁面,只將下壁面作為可滲透的壁面(見圖1)。采用用戶自定義函數(shù)UDF將納濾傳質(zhì)模型和CFD邊界條件進行耦合,以此模擬膜表面濃度變化曲線,克服了CFD商業(yè)軟件的限制,能夠較準確地預(yù)測膜表面的溶質(zhì)濃度變化。然而,該研究把上膜面作為非滲透性膜以減小CFD計算的復(fù)雜程度,這樣的假設(shè)與實際卷式膜組件中的上、下均為滲透性膜的情況不一致,并且該模擬在卷式膜組件的幾何模型上也作了較大簡化,所以模擬結(jié)果具有一定的局限性。因此,建立更加完整和真實的物理模型和幾何模型,結(jié)合CFD計算模擬是深入研究卷式膜過程的重要途徑。

    圖1 薄膜理論中的邊界條件和濃差極化厚度[8]Fig.1 Boundary condition and concentation in modified film theroy[8]

    上膜邊界條件(非滲透性):

    式(1)中,u為x方向的速度,m/s;v為y方向的速度,m/s;mA為溶質(zhì)質(zhì)量分數(shù)。

    下膜邊界條件(納濾膜):

    式(2)中,Jv為滲透通量,m/s;DAB為雙向擴散系數(shù),m2/s;mAc為模擬區(qū)域緊挨納濾膜表面的單元格質(zhì)心處的溶質(zhì)質(zhì)量分數(shù);R為截留率;R′為透鹽率,1-R;δc為下膜表面單元與相鄰單元質(zhì)心之間的距離,m。

    其中,Jv由Spiegler-Kedem模型確定:

    式(3)中,Lp為水力學(xué)系數(shù),m/(Pa?s);P為跨膜壓力,Pa;mAw為壁面附近的溶質(zhì)質(zhì)量分數(shù);γ為滲透壓常數(shù),Pa;σ為折射系數(shù)。

    我們發(fā)現(xiàn),雖然近年來CFD模擬在膜污染研究中取得了一些進展,但是由于膜污染過程的復(fù)雜性,目前的模擬結(jié)果還難以較為真實地模擬膜污染過程。同樣,對于CFD在納濾和反滲透壓力驅(qū)動的膜過程傳質(zhì)研究中,膜傳質(zhì)模型的選擇主要集中在非平衡熱力學(xué)方面,不能全面、準確地描述不同體系的真實傳質(zhì)過程。因此需要針對各類膜污染以及傳質(zhì)機理等建立有效的CFD耦合模型,對機理進行深入研究。

    3.1.2 CFD在膜組件隔網(wǎng)優(yōu)化設(shè)計研究中的應(yīng)用

    在實際工程操作中,人們經(jīng)常在膜組件中設(shè)置隔網(wǎng)作為湍流強化器,增加流體湍流強度,產(chǎn)生二級流,增大膜表面的剪切力,提高傳質(zhì)效率,從而減小濃差極化和膜污染現(xiàn)象。同時由于隔網(wǎng)的存在,增加了流體阻力,引起壓降增加,加大了能耗,并且可能出現(xiàn)局部溝流、死區(qū)等,如何優(yōu)化膜組件中的隔網(wǎng)成為研究者們感興趣的問題。

    一些研究者認為隔網(wǎng)的距離對流道中流體流動以及漩渦具有較大影響,并對此進行了相關(guān)研究。例如Cao[10]等用CFD模擬了湍流狀態(tài)下橫截面方向設(shè)置原料間隔器對場分布的影響,并給出速度場,剪切力,剪切力場和湍動能分布圖。發(fā)現(xiàn)合理設(shè)置原料間隔器及其之間的距離對提高膜性能有效:適當(dāng)?shù)販p小原料間隔器之間的距離,可以減小剪切力峰值間的距離,增加漩渦,從而提高膜表面的傳質(zhì);但另一方面,會增大壓力降,從而增加成本。因此他們建議從考慮合理的原料間隔器的間距和適當(dāng)?shù)貕航翟O(shè)置的思路來優(yōu)化原料間隔器位置的設(shè)置。

    此外,人們發(fā)現(xiàn)除了間隔器距離,隔網(wǎng)的形狀同樣會對膜組件內(nèi)流體流動產(chǎn)生一定影響。Ahmad[11]等采用CFD技術(shù)對矩形、圓形和三角形隔網(wǎng)細絲(見圖2)的流道中均產(chǎn)生不穩(wěn)定漩渦的情況進行研究,模擬結(jié)果發(fā)現(xiàn)在相同的雷諾數(shù)下三種不同形狀的隔網(wǎng)產(chǎn)生漩渦的情況有所不同。

    圖2 隔網(wǎng)細絲形狀和尺寸[11]Fig.2 Shape and dimension of spacer filaments[11]

    目前,在膜組件的隔網(wǎng)優(yōu)化設(shè)計方面,研究者們大部分的研究集中于形狀簡單的單層隔網(wǎng)的CFD研究,對于形狀復(fù)雜的雙層隔網(wǎng)或者多層隔網(wǎng)的研究則較少。

    3.2 CFD在滲透汽化和膜蒸餾技術(shù)研究中的應(yīng)用

    與壓力驅(qū)動膜過程不同,具有相變的膜過程存在著明顯的溫差現(xiàn)象,過程中的傳質(zhì)與傳熱共同作用導(dǎo)致了溫差極化。因此,具有相變的膜過程性能和效率受濃差極化和溫差極化的共同影響,在相應(yīng)CFD研究中,研究傳質(zhì)過程的同時也需要對過程的傳熱進行研究。

    3.2.1 CFD在滲透汽化技術(shù)研究中的應(yīng)用

    滲透汽化分離技術(shù)是利用液體混合物中組分在致密膜中溶解擴散性能不同實現(xiàn)分離的膜技術(shù)。具有分離效率高、設(shè)備簡單、能耗低等優(yōu)點,近年來在水處理應(yīng)用中也備受關(guān)注。

    在滲透汽化理論傳質(zhì)研究中普遍認為,與原料側(cè)和膜本身的傳質(zhì)阻力相比,氣相側(cè)的傳質(zhì)阻力可以忽略,并且由于滲透汽化具有相變的特點,因此在建立CFD計算模型時只考慮膜上游側(cè)液體的傳質(zhì)情況,并主要依賴于溶解-擴散機理。

    Sean[12]等模擬滲透汽化平板膜器的狹窄流道內(nèi)雙組分的傳質(zhì)。在建立物理模型時假設(shè):溶質(zhì)溶解于膜的過程速度很快;與液體邊界層的傳質(zhì)阻力相比,來自膜的傳質(zhì)阻力很小;膜滲透側(cè)氣相的傳質(zhì)忽略不計。從而只對膜上游流道建模計算,采用膜上游面的擬一階化學(xué)反應(yīng)機理來描述溶質(zhì)濃差極化現(xiàn)象,并將此作為計算的一個邊界條件。采用Fluent6.0軟件模擬膜器內(nèi)的流動,得到流道內(nèi)的速度分布和狹窄膜通道濃度邊界層的溶質(zhì)濃度分布。他們還研究了膜器中擋板高度對增加傳質(zhì)的影響,模擬結(jié)果發(fā)現(xiàn),滲透汽化過程的平均傳質(zhì)系數(shù)隨著擋板高度的增加而增加。

    在滲透汽化傳熱研究方面,Gert[13]等利用CFD模擬陶瓷滲透汽化膜過程中的溫差極化現(xiàn)象,結(jié)果表明溫差極化會導(dǎo)致明顯的黏度梯度,因此嚴重影響傳質(zhì)效率。此外,Gert[14]還利用CFX4.2軟件包計算三維管式滲透汽化膜器的二級流分布和溫度分布,并與超音速計算機X線斷層攝影術(shù)測量的溫度分布相比較,計算結(jié)果與測量數(shù)據(jù)吻合良好。

    從以上CFD在滲透汽化膜過程中的研究中不難發(fā)現(xiàn),大部分研究者只考慮原料側(cè)的操作參數(shù)以及原料側(cè)傳質(zhì)過程的研究,忽略了滲透側(cè)的傳質(zhì)阻力,然而滲透側(cè)的傳質(zhì)阻力卻真實存在。因此,我們需要建立完整的滲透汽化過程物理模型進行CFD模擬,才能對過程進行更加全面和準確地認識和分析。

    3.2.2 CFD在膜蒸餾技術(shù)研究中的應(yīng)用

    膜蒸餾作為一種新型的膜分離技術(shù),在海水淡化、廢水處理等方面具有廣泛研究。根據(jù)對膜冷凝側(cè)蒸汽冷凝方式的不同,膜蒸餾可分為:直接接觸式膜蒸餾(DCMD)、氣隙式膜蒸餾(AGMD)、氣掃式膜蒸餾(SGMD)和真空膜蒸餾(VMD)。

    膜蒸餾是一種熱驅(qū)動過程,其推動力是疏水膜材料兩側(cè)的蒸汽壓差。膜蒸餾過程除了受濃差極化的影響,溫差極化更大程度上限制了該過程的效率。因此,在CFD研究中常常需要將傳質(zhì)與傳熱過程進行耦合。

    CFD在膜蒸餾研究中的應(yīng)用主要集中在直接DCMD,這是因為相比于其他三種膜蒸餾過程,DCMD組件無需額外的冷凝設(shè)備,過程傳質(zhì)阻力小,因此被廣泛研究。Shirazian[15]等利用CFD對DCMD中水蒸汽在微孔膜中的傳遞過程進行模擬,主要考察了不同操作參數(shù)如氣體流速、液體流速以及氣相濕度等對膜過程性能的影響。從模擬結(jié)果可以看出,氣體流速越大,膜蒸餾的效率越高,這是由于大的氣體流速能夠避免過程達到熱力學(xué)平衡。此外,他們還研究了膜結(jié)構(gòu)的影響,發(fā)現(xiàn)膜材料的孔隙率和曲率之比對膜性能沒有顯著影響。

    Yu[16]利用CFD對層流狀態(tài)下的DCMD過程中傳熱和傳質(zhì)過程做了模擬研究。將潛熱耦合進熱量平衡方程建立了二維傳熱模型,結(jié)合Navies-Stokes方程,對逆流流型下的熱邊界層厚度、膜壁溫度、溫度極化系數(shù)(TPC)、局部換熱系數(shù)、局部質(zhì)量通量以及熱效率進行預(yù)測,膜通量和膜面溫度相關(guān)的中空纖維膜實驗結(jié)果吻合。Shirazian[17]等也用CFD對錯流平板膜蒸餾裝置的傳質(zhì)與傳熱過程進行了模擬。他們將膜孔內(nèi)的努森擴散和黏性流模型與CFD模型耦合,能較為準確地模擬了微孔膜內(nèi)的流體流動狀況。模擬結(jié)果表明,在靠近膜表面區(qū)域,溫度差異性非常顯著,這是由于在膜表面存在著溫度邊界層(見圖3)。

    圖3 直接接觸式膜蒸餾過程中的溫度分布[17]Fig.3 Temperature distribution in direct-contact membrane distillation[17]

    此外,CFD在真空膜蒸餾過程研究中的應(yīng)用也逐漸變多。Tang[18]等采用CFD的離散相模型模擬了VMD濃縮NaCl水溶液的過程。模擬中將膜材料作為多孔介質(zhì)處理,用UDF定義了水蒸汽蒸發(fā)過程,研究了多孔材料內(nèi)部的質(zhì)量傳遞和熱量傳遞過程。此外,他們還通過模擬計算得到了不同溫度下的氣相體積分數(shù)。在該模擬中,作者合理地將微孔膜材料定義為多孔介質(zhì),一定程度上體現(xiàn)了膜的滲透性能。另一方面,他們采用CFD中的氣液離散相模型有效地模擬了料液溫度對VMD過程傳質(zhì)和滲透效率的影響。

    目前,CFD模擬在膜蒸餾過程傳質(zhì)、傳熱以及組件優(yōu)化設(shè)計方面取得了一些進展,但不難發(fā)現(xiàn)這些研究大部分集中在DCMD、VMD次之。而在另外兩種膜蒸餾過程,AGMD和SGMD中的研究則十分少見。因此,結(jié)合CFD模擬進行相應(yīng)過程的研究對深入和全面了解膜蒸餾過程是十分有必要的。

    在膜蒸餾組件優(yōu)化設(shè)計方面,CFD的應(yīng)用主要集中在組件中的隔網(wǎng)以及組件形狀對流體流動影響的研究中。Shakaib[19]等用CFD考察了有隔網(wǎng)的膜蒸餾流道中瞬時流動和溫度分布。在雷諾數(shù)較大的情況下,瞬變流動出現(xiàn)在隔網(wǎng)后面,并且沿著流道逐漸消失。這種非穩(wěn)態(tài)行為導(dǎo)致了局部溫度和換熱系數(shù)隨著時間變化。Yang[20]等利用CFD考察了一系列不同形狀的中空纖維膜蒸餾組件對DCMD性能的影響,模擬結(jié)果表明:與傳統(tǒng)的圓柱型中空纖維膜組件相比,齒輪型纖維膜組件具有最大的溫差極化系數(shù)和膜通量,交替波浪型中空纖維膜組件次之。這種改造之后的膜組件之所以能提高性能,是由于其結(jié)構(gòu)能夠產(chǎn)生強烈的二級流,從而提高膜表面的湍動和迅速更新。

    從上述研究我們可以發(fā)現(xiàn),無論是添加隔網(wǎng)還是改變組件的彎曲度,都能夠在一定程度上增加組件內(nèi)部流體流動的湍動情況,這是提高膜蒸餾過程傳質(zhì)和傳熱的重要途徑。在膜蒸餾強化傳質(zhì)的研究中,除了膜組件優(yōu)化設(shè)計外,超聲、鼓泡以及非穩(wěn)態(tài)流動等方式均能增加一定的流體湍動情況。然而,單純地實驗現(xiàn)象不能深入了解這些方式對流體流動的具體作用,因此CFD技術(shù)在這方面的研究仍然存在很大的應(yīng)用潛力。

    3.3 CFD在其他膜過程中的應(yīng)用

    除了上述的膜分離技術(shù)外,近年來隨著膜技術(shù)的發(fā)展,其他膜分離技術(shù),如膜生物反應(yīng)器、正滲透等膜分離技術(shù)也得到廣泛的研究和進一步的實際應(yīng)用。在這些過程機理和組件設(shè)計的研究中,CFD模擬同樣具有很大作用并具有相關(guān)的應(yīng)用。

    于艷[21]等用歐拉模型研究了玻璃纖維管式膜生物反應(yīng)器內(nèi)的水力學(xué)特征,定量得到了膜面水流流速、膜面剪切力、質(zhì)量流率等參數(shù),實測流速與模擬的流速在膜單元流場混合充分處較為吻合。研究了不同布置高度膜組件膜面流速和剪切力分布情況,得出膜組件距離底部曝氣器有一個最佳區(qū)域,在此位置時水力循環(huán)良好。

    近年來,正滲透作為一種新型的膜分離技術(shù)得到越來越多的關(guān)注。與反滲透過程類似,正滲透過程性能同樣受到濃差極化現(xiàn)象的影響;但與反滲透不同的是,反滲透過程只在進料側(cè)受到濃差極化的影響,而正滲透過程中膜兩側(cè)均存在不可忽略的濃差極化,包括膜材料有效層表面的濃差極化以及多空支撐層內(nèi)部的濃差極化。因此,CFD在正滲透研究中需要建立完整的非對稱膜兩側(cè)的區(qū)域模型進行模擬。Minkyu[22]、Gruber[23]等均采用CFD對正滲透過程進行了完整的模擬。

    4 CFD在膜技術(shù)領(lǐng)域研究應(yīng)用的總結(jié)與展望

    CFD技術(shù)因其自身優(yōu)點在膜技術(shù)方面的研究涉及的范圍逐漸擴大,前期主要集中在比較成熟的反滲透、超濾、微濾、納濾等壓力驅(qū)動膜分離技術(shù)的模擬;近年來,CFD在新型的膜技術(shù)如滲透汽化、膜蒸餾以及膜生物反應(yīng)器中的應(yīng)用也逐漸增多。CFD技術(shù)在膜技術(shù)研究中雖然能較好地解決一些問題,但是同樣存在一些不足和發(fā)展的空間。

    1)CFD是一門理論與實際相結(jié)合的交叉科學(xué),對研究者的知識水平和計算機能力要求高。對于實際問題的模擬,依賴于初始條件和操作條件的確定,如設(shè)置不當(dāng)將導(dǎo)致模擬結(jié)果不準確。因此,不能單獨依靠CFD模擬完成研究工作,需要與實驗相結(jié)合,模擬結(jié)果由實驗數(shù)據(jù)驗證和校正。因此,粒子圖像測速技術(shù)(PIV)、PDA等實驗流體技術(shù)需要與CFD模擬技術(shù)相輔相成。

    2)在膜傳質(zhì)機理和膜污染機理方面,應(yīng)該更加深入研究相關(guān)機理與CFD模型耦合,進一步探究濃差極化和污染層的形成,提出相應(yīng)的緩解措施,提高膜通量,提升工藝穩(wěn)定性。

    3)目前大部分研究集中在膜組件內(nèi)部的流體力學(xué)以及膜表面流動狀況,在超濾、微濾等多孔材料內(nèi)部的流動模擬較少。因此,可以結(jié)合相應(yīng)CFD模型和物理模型開展這方面的研究,在一定程度上指導(dǎo)膜材料的設(shè)計研究。

    4)一些湍動強化過程在增強膜過程的傳遞效率研究中應(yīng)用廣泛,因此可以利用CFD技術(shù)對這些耦合過程進行研究,明確強化過程對流體流動的影響。

    5)由于3D物理模型和計算模型的建立比較復(fù)雜,計算量也非常大,因此目前3D模擬主要集中在膜組件優(yōu)化設(shè)計方面;而在膜傳質(zhì)傳熱以及膜污染現(xiàn)象方面,3D模擬計算仍然比較缺乏。

    [1]王福軍.計算流體動力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社有限公司,2004.

    [2]員文權(quán),楊慶峰.計算流體力學(xué)在反滲透膜分離中的應(yīng)用[J].化工進展,2008,27(9):1357-1362.

    [3]何娟娟,黃衛(wèi)星,肖澤儀,等.計算流體力學(xué)在膜分離技術(shù)中的應(yīng)用[J].化工裝備技術(shù),2006,27(2):14-16.

    [4]Fimbres-Weihs G A,Wiley D E.Review of 3D CFD modeling of flow and mass transfer in narrow spacer-filled channels in membrane modules[J].Chemical Engineering and Processing:Process Intensification,2010,49(7):759-781.

    [5]Kapellos George E,Alexiou Terpsichori S,Payatakes Alkiviades C.Hierarchical simulato,r of biofilm growth and dynamics in granular porous materials[J].Advances in water resources,2007,30(6):1648-1667.

    [6]von der Schulenburg D A,Pintelon T R R,Picioreanu C,et al.Three dimensional simulations of biofilm growth in porous media[J].Aiche Journal,2009,55(2):494-504.

    [7]Geraldes Vítor,Semi?o Viriato,de Pinho Maria Norberta.Flow and mass transfer modelling of nanofiltration[J].Journal of Membrane Science,2001,191(1):109-128.

    [8]Ahmad AL,Lau K K,Bakar M Z,et al.Integrated CFD simulation of concentration polarization in narrow membrane channel[J].ComputersandChemicalEngineering,2005,29(10):2087-2095.

    [9]Ahmad A L,Lau K K,Abu Bakar M Z.Impact of different spacer filament geometries on concentration polarization control in narrow membrane channel[J].Journal of Membrane Science,2005,262(1):138-152.

    [10]Cao Z,Wiley D E,F(xiàn)ane A G.CFD simulations of net-type turbulence promoters in a narrow channel[J].Journal of Membrane Science,2001,185(2):157-176.

    [11]Ahmad A L,Lau K K.Impact of different spacer filaments geometries on 2D unsteady hydrodynamics and concentration polarization in spiral wound membrane channel[J].Journal of Membrane Science,2006,286(1):77-92.

    [12]Peng Ming,Vane Leland M,Liu Sean X.Numerical simulation of concentration polarization in a pervaporation module[J].Separation Science Technology,2005,39(6):1239-1257.

    [13]Van der Gulik Gert-Jan S,Janssen,R E G,Wijers J G,et al.Hydrodynamics in a ceramic pervaporation membrane reactor for resin production[J].Chemical Engineering Science,2001,56(2):371-379.

    [14]van der Gulik Gert-Jan S,Wijers Johan G,Keurentjes Jos T F.Measurement of 2D-temperature distributions in a pervaporation membrane module using ultrasonic computer tomography and comparison with computational fluid dynamics calculations[J].Journal of Membrane Science,2002,204(1):111-124.

    [15]Shirazian Saeed,Ashrafizadeh Seyed N.3D modeling and simulation of mass transfer in vapor transport through porous membranes[J].Chemical Engineering and Technology,2013,36(1):177-185.

    [16]Yu Hui,Yang Xing,Wang Rong,et al.Numerical simulation of heat and mass transfer in direct membrane distillation in a hollow fiber module with laminar flow[J].Journal of Membrane Science,2011,384(1):107-116.

    [17]Ghadiri Mehdi,F(xiàn)akhri Safoora,Shirazian Saeed.Modeling and CFD simulation of water desalination using nanoporous membrane contactors[J].Industrial&Engineering Chemistry Research,2013,52(9):3490-3498.

    [18]Tang Na,Zhang Huanju,Wang Wei.Computational fluid dynamics numerical simulation of vacuum membrane distillation for aqueous NaCl solution[J].Desalination,2011,274(1):120-129.

    [19]Shakaib M,Hasani S M F,Haque M Ehtesham-ul,et al.A CFD study of heat transfer through spacer channels of membrane distillation modules[J].Desalination and Water Treatment,2013,51(16-18):3662-3674.

    [20]Yang Xing,Yu Hui,Wang Rong,et al.Optimization of microstructured hollow fiber design for membrane distillation applications using CFD modeling[J].Journal of Membrane Science,2012,421:258-270.

    [21]于 艷,樊耀波,徐國良,等.計算流體力學(xué)對膜生物反應(yīng)器水力學(xué)特征的模擬研究[J].膜科學(xué)與技術(shù),2011,31(4):9-15.

    [22]Park Minkyu,Lee Ji Jung,Lee Sangho,et al.Determination of a constant membrane structure parameter in forward osmosis processes[J].JournalofMembraneScience,2011,375(1):241-248.

    [23]Gruber M F,Johnson C J,Tang C Y,et al.Computational fluid dynamics simulations of flow and concentration polarization in forward osmosis membrane systems[J].Journal of Membrane Science,2011,379(1):488-495.

    Computational fluid dynamics applied to membrane processes for water
    treatment

    Zhang Yaqin1,Zhang Lin1,Hou Li’an1,2
    (1.Department of Chemical and Biological Engineering,Zhejiang University,Hangzhou 310027,China;2.Institute for Logistic Science and Technology of the Second Artillery,Beijing 100011,China)

    In this paper,the basic principle of computational fluid dynamics(CFD)in membrane separation process were introduced,and the application of CFD mechanism of membrane pollution,membrane module design in pressure-driven membrane processes,as well as other membrane process,such as pervaporation,membrane distillation were presented.The prospect of application of CFD in membrane process was discussed.

    CFD;membrane separation;pressure-driven;pervaporation;membrane distillation

    TQ028.8

    A

    1009-1742(2014)07-0047-06

    2014-05-05

    國家自然科學(xué)基金重點資助項目(51238006)

    張 林,1972年出生,男,安徽當(dāng)涂縣人,教授,主要研究方向為膜科學(xué)與技術(shù);E-mail:linzhang@ziu.edu.cn

    猜你喜歡
    膜分離傳質(zhì)極化
    認知能力、技術(shù)進步與就業(yè)極化
    探析膜分離技術(shù)在環(huán)境工程中的應(yīng)用及發(fā)展
    化香樹果序多酚膜分離的動力學(xué)研究
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    基于PWM控制的新型極化電源設(shè)計與實現(xiàn)
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實驗研究
    煤顆粒熱解的傳熱傳質(zhì)分析
    計算物理(2014年1期)2014-03-11 17:00:36
    膜分離技術(shù)應(yīng)用于喉咽清口服液純化工藝的研究
    中成藥(2014年9期)2014-02-28 22:29:09
    中文字幕高清在线视频| 亚洲aⅴ乱码一区二区在线播放 | 精品欧美一区二区三区在线| 天天一区二区日本电影三级| 国产爱豆传媒在线观看 | 真人一进一出gif抽搐免费| 国产激情久久老熟女| 中文字幕精品免费在线观看视频| 97人妻精品一区二区三区麻豆 | 高清在线国产一区| 757午夜福利合集在线观看| 村上凉子中文字幕在线| 欧美性猛交╳xxx乱大交人| 在线观看免费午夜福利视频| 真人做人爱边吃奶动态| 国产精品综合久久久久久久免费| 亚洲精品美女久久久久99蜜臀| 视频在线观看一区二区三区| 中文字幕最新亚洲高清| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久久久免费视频了| 老司机靠b影院| 免费电影在线观看免费观看| 真人一进一出gif抽搐免费| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| 国产在线观看jvid| 最近最新中文字幕大全免费视频| 成人亚洲精品av一区二区| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久av美女十八| 母亲3免费完整高清在线观看| www国产在线视频色| 久久久久免费精品人妻一区二区 | 在线观看午夜福利视频| 亚洲人成网站高清观看| 亚洲全国av大片| 变态另类丝袜制服| 中文字幕精品亚洲无线码一区 | 欧美黄色片欧美黄色片| 别揉我奶头~嗯~啊~动态视频| 久久久久久久久免费视频了| 国产一区二区三区视频了| 88av欧美| 后天国语完整版免费观看| 精品第一国产精品| 老司机靠b影院| АⅤ资源中文在线天堂| 两人在一起打扑克的视频| 久久久久久久久中文| 亚洲全国av大片| 琪琪午夜伦伦电影理论片6080| 久久香蕉激情| 香蕉久久夜色| 一区二区日韩欧美中文字幕| 最近最新免费中文字幕在线| 淫妇啪啪啪对白视频| 精品国产亚洲在线| 久久狼人影院| 午夜a级毛片| 一卡2卡三卡四卡精品乱码亚洲| 特大巨黑吊av在线直播 | 精品福利观看| 精品国产超薄肉色丝袜足j| 人人妻人人澡欧美一区二区| 免费电影在线观看免费观看| 亚洲中文av在线| 欧美日韩亚洲综合一区二区三区_| av视频在线观看入口| 一本一本综合久久| 中出人妻视频一区二区| 国产精品野战在线观看| 国产精品98久久久久久宅男小说| 亚洲国产精品sss在线观看| 久久婷婷人人爽人人干人人爱| 亚洲成av片中文字幕在线观看| 国产精品日韩av在线免费观看| 免费看美女性在线毛片视频| 亚洲精华国产精华精| 老司机午夜十八禁免费视频| 国产人伦9x9x在线观看| 成人18禁高潮啪啪吃奶动态图| 91老司机精品| 免费一级毛片在线播放高清视频| 无限看片的www在线观看| 黄网站色视频无遮挡免费观看| 日韩欧美国产一区二区入口| 精品国产一区二区三区四区第35| 久久99热这里只有精品18| 国产一区二区在线av高清观看| 国产亚洲欧美98| 丰满人妻熟妇乱又伦精品不卡| 夜夜躁狠狠躁天天躁| 曰老女人黄片| 日本成人三级电影网站| 一本一本综合久久| 高清在线国产一区| 不卡av一区二区三区| 欧美不卡视频在线免费观看 | 看片在线看免费视频| 老司机深夜福利视频在线观看| 亚洲专区国产一区二区| 麻豆久久精品国产亚洲av| 国产亚洲精品久久久久久毛片| 神马国产精品三级电影在线观看 | 国产一卡二卡三卡精品| 淫秽高清视频在线观看| 国产精品 国内视频| 午夜老司机福利片| 日本在线视频免费播放| 日日摸夜夜添夜夜添小说| 白带黄色成豆腐渣| 国产精品影院久久| a级毛片在线看网站| 久久国产精品人妻蜜桃| 国产高清视频在线播放一区| 男女之事视频高清在线观看| 久久热在线av| 2021天堂中文幕一二区在线观 | 亚洲精品中文字幕在线视频| 国产精品爽爽va在线观看网站 | 男人舔女人的私密视频| 欧美zozozo另类| 精品久久久久久久久久久久久 | 免费搜索国产男女视频| 国产aⅴ精品一区二区三区波| 每晚都被弄得嗷嗷叫到高潮| 日韩视频一区二区在线观看| 老司机福利观看| 国产精品1区2区在线观看.| 88av欧美| 久久草成人影院| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 日韩有码中文字幕| 欧美不卡视频在线免费观看 | 桃红色精品国产亚洲av| 在线观看舔阴道视频| 久久国产精品男人的天堂亚洲| 午夜福利在线观看吧| 人妻丰满熟妇av一区二区三区| 日韩免费av在线播放| 国产亚洲精品久久久久5区| 法律面前人人平等表现在哪些方面| 久久精品91蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 午夜成年电影在线免费观看| www.自偷自拍.com| 欧美成人一区二区免费高清观看 | 日本在线视频免费播放| 久久精品人妻少妇| 在线观看午夜福利视频| 国产伦一二天堂av在线观看| 黄片小视频在线播放| 亚洲人成电影免费在线| 国产精品亚洲美女久久久| 久久精品国产亚洲av香蕉五月| 在线观看免费日韩欧美大片| 欧美 亚洲 国产 日韩一| 香蕉av资源在线| 可以在线观看的亚洲视频| 中文字幕精品免费在线观看视频| 一二三四社区在线视频社区8| 国产日本99.免费观看| 亚洲一区中文字幕在线| 美女高潮到喷水免费观看| 亚洲成国产人片在线观看| 日本成人三级电影网站| 制服丝袜大香蕉在线| 国产高清videossex| 亚洲全国av大片| 亚洲av熟女| 黄片播放在线免费| www日本黄色视频网| 久久九九热精品免费| 国产成人精品久久二区二区免费| 在线永久观看黄色视频| 国产视频内射| 亚洲人成网站在线播放欧美日韩| 成人亚洲精品一区在线观看| 天堂影院成人在线观看| 亚洲国产毛片av蜜桃av| 在线观看免费日韩欧美大片| 久久久久国产精品人妻aⅴ院| 亚洲精品色激情综合| av片东京热男人的天堂| 999精品在线视频| 国产伦人伦偷精品视频| 精品国产乱子伦一区二区三区| 国产野战对白在线观看| 一边摸一边抽搐一进一小说| 韩国av一区二区三区四区| 久久婷婷人人爽人人干人人爱| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产精品影院| а√天堂www在线а√下载| 黄色a级毛片大全视频| 国产野战对白在线观看| 在线观看免费午夜福利视频| 亚洲人成网站在线播放欧美日韩| 最近最新中文字幕大全免费视频| 亚洲午夜理论影院| 国产精品二区激情视频| 国产久久久一区二区三区| 亚洲全国av大片| 欧美在线一区亚洲| 久久 成人 亚洲| 草草在线视频免费看| 热99re8久久精品国产| 成年版毛片免费区| 嫩草影院精品99| 欧美大码av| 午夜影院日韩av| 一本久久中文字幕| 精品福利观看| 国产精品爽爽va在线观看网站 | 国产在线观看jvid| 亚洲五月婷婷丁香| 日日夜夜操网爽| 天天躁狠狠躁夜夜躁狠狠躁| 一夜夜www| 91麻豆av在线| 国产成人精品久久二区二区免费| 亚洲精品中文字幕在线视频| 久久精品影院6| 99久久国产精品久久久| 一级毛片女人18水好多| 亚洲成人免费电影在线观看| 欧美黄色片欧美黄色片| 真人一进一出gif抽搐免费| 啦啦啦韩国在线观看视频| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 精品一区二区三区四区五区乱码| 一本综合久久免费| 国产不卡一卡二| 亚洲五月天丁香| 制服人妻中文乱码| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 国产亚洲精品av在线| 免费电影在线观看免费观看| 日本三级黄在线观看| 男女午夜视频在线观看| 三级毛片av免费| 一边摸一边做爽爽视频免费| 久久亚洲真实| 欧美精品啪啪一区二区三区| 久久久精品欧美日韩精品| xxx96com| 丰满的人妻完整版| 老司机福利观看| 97碰自拍视频| 亚洲人成电影免费在线| 欧美亚洲日本最大视频资源| 大香蕉久久成人网| 久热爱精品视频在线9| 免费看美女性在线毛片视频| 亚洲av片天天在线观看| 久久香蕉国产精品| 我的亚洲天堂| 欧美黄色片欧美黄色片| 午夜久久久久精精品| 天天躁狠狠躁夜夜躁狠狠躁| 女警被强在线播放| 一进一出抽搐gif免费好疼| 99国产精品一区二区蜜桃av| 搞女人的毛片| 90打野战视频偷拍视频| 黄色毛片三级朝国网站| 久久久水蜜桃国产精品网| 国产成人啪精品午夜网站| 人人妻人人澡欧美一区二区| 中文字幕人妻熟女乱码| 亚洲成人久久性| 亚洲精华国产精华精| 欧美乱色亚洲激情| 99国产综合亚洲精品| 免费观看精品视频网站| 男女床上黄色一级片免费看| 法律面前人人平等表现在哪些方面| 国产亚洲精品第一综合不卡| x7x7x7水蜜桃| 韩国精品一区二区三区| 人人妻人人澡人人看| 精品国产亚洲在线| 久久精品成人免费网站| 亚洲国产精品合色在线| 日本a在线网址| 亚洲国产欧美网| 亚洲熟妇中文字幕五十中出| 50天的宝宝边吃奶边哭怎么回事| 黄片小视频在线播放| 99热这里只有精品一区 | 欧美乱色亚洲激情| 丝袜美腿诱惑在线| 黑人操中国人逼视频| 亚洲精品美女久久av网站| 伦理电影免费视频| 两个人免费观看高清视频| 精品久久蜜臀av无| 丁香欧美五月| 婷婷丁香在线五月| 美女免费视频网站| 可以免费在线观看a视频的电影网站| 韩国精品一区二区三区| 久久这里只有精品19| 国产高清视频在线播放一区| 91国产中文字幕| 变态另类成人亚洲欧美熟女| 日日干狠狠操夜夜爽| 韩国av一区二区三区四区| 美女高潮到喷水免费观看| 我的亚洲天堂| 天堂√8在线中文| 日本一本二区三区精品| 中文字幕人妻丝袜一区二区| 国内精品久久久久精免费| 宅男免费午夜| 欧美大码av| 少妇被粗大的猛进出69影院| √禁漫天堂资源中文www| 欧美在线黄色| 亚洲男人的天堂狠狠| 亚洲国产精品sss在线观看| 99久久精品国产亚洲精品| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 老汉色∧v一级毛片| 亚洲五月婷婷丁香| 亚洲一区二区三区色噜噜| 久久久久九九精品影院| 在线观看www视频免费| 又大又爽又粗| 亚洲人成网站在线播放欧美日韩| 大香蕉久久成人网| 久久精品91蜜桃| 国产成人影院久久av| 国产精品免费视频内射| 久久久久国内视频| 一边摸一边做爽爽视频免费| 91国产中文字幕| 亚洲中文字幕日韩| 两性夫妻黄色片| 日本成人三级电影网站| 在线观看免费日韩欧美大片| 亚洲一区二区三区色噜噜| 在线观看免费日韩欧美大片| 色哟哟哟哟哟哟| 91国产中文字幕| 夜夜看夜夜爽夜夜摸| 丝袜人妻中文字幕| 在线观看66精品国产| 久久天躁狠狠躁夜夜2o2o| 亚洲av五月六月丁香网| 久久精品成人免费网站| 日韩中文字幕欧美一区二区| 一本一本综合久久| 日韩免费av在线播放| 日韩有码中文字幕| 婷婷精品国产亚洲av在线| 亚洲人成伊人成综合网2020| 欧美日本亚洲视频在线播放| 免费看日本二区| 成人三级黄色视频| 午夜久久久在线观看| avwww免费| 欧美中文综合在线视频| 在线永久观看黄色视频| 国产精品一区二区免费欧美| 19禁男女啪啪无遮挡网站| 日日夜夜操网爽| 女生性感内裤真人,穿戴方法视频| 啦啦啦韩国在线观看视频| 在线天堂中文资源库| 人人澡人人妻人| 中文字幕人妻丝袜一区二区| 免费一级毛片在线播放高清视频| 亚洲狠狠婷婷综合久久图片| 久久精品亚洲精品国产色婷小说| 欧美成人性av电影在线观看| 婷婷精品国产亚洲av| 亚洲五月天丁香| 日本成人三级电影网站| tocl精华| 日本免费一区二区三区高清不卡| 午夜老司机福利片| 亚洲欧美日韩无卡精品| 欧洲精品卡2卡3卡4卡5卡区| 宅男免费午夜| 成人一区二区视频在线观看| 精品久久久久久久毛片微露脸| 久久精品国产综合久久久| 俄罗斯特黄特色一大片| 日韩精品中文字幕看吧| 99riav亚洲国产免费| 成人特级黄色片久久久久久久| 成人手机av| 成人av一区二区三区在线看| aaaaa片日本免费| 日韩成人在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 久久精品aⅴ一区二区三区四区| www.熟女人妻精品国产| 99精品欧美一区二区三区四区| 日韩大尺度精品在线看网址| 国产精品日韩av在线免费观看| 一区二区三区国产精品乱码| 欧美激情高清一区二区三区| 国产乱人伦免费视频| 日韩大码丰满熟妇| 久久久久久人人人人人| 精品一区二区三区四区五区乱码| 亚洲第一av免费看| 国产成人啪精品午夜网站| 少妇被粗大的猛进出69影院| 高潮久久久久久久久久久不卡| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久人妻蜜臀av| 国产精品精品国产色婷婷| 国产免费男女视频| 1024视频免费在线观看| 免费在线观看影片大全网站| 国内精品久久久久久久电影| 中文字幕人成人乱码亚洲影| 国产片内射在线| 亚洲精品色激情综合| 亚洲精品av麻豆狂野| 亚洲狠狠婷婷综合久久图片| 亚洲一区中文字幕在线| 久久性视频一级片| 欧美+亚洲+日韩+国产| 性欧美人与动物交配| 国内揄拍国产精品人妻在线 | 久久狼人影院| 两个人看的免费小视频| videosex国产| 岛国在线观看网站| 一区二区三区精品91| 日韩免费av在线播放| 亚洲免费av在线视频| 亚洲国产看品久久| 国产片内射在线| 日韩欧美一区二区三区在线观看| 久久青草综合色| 亚洲熟妇熟女久久| 无限看片的www在线观看| 欧美在线一区亚洲| 妹子高潮喷水视频| xxxwww97欧美| 久久精品aⅴ一区二区三区四区| 精品第一国产精品| or卡值多少钱| 老司机在亚洲福利影院| 好男人在线观看高清免费视频 | 国产久久久一区二区三区| 在线免费观看的www视频| 看片在线看免费视频| 制服丝袜大香蕉在线| av免费在线观看网站| 后天国语完整版免费观看| 男人舔奶头视频| 色综合欧美亚洲国产小说| 国产91精品成人一区二区三区| 亚洲精品av麻豆狂野| 中文字幕另类日韩欧美亚洲嫩草| 黄色成人免费大全| 在线观看日韩欧美| 18禁黄网站禁片免费观看直播| 级片在线观看| 成人免费观看视频高清| 成人国产综合亚洲| 1024手机看黄色片| 久久国产乱子伦精品免费另类| 国产精品 欧美亚洲| 国产欧美日韩精品亚洲av| 国产熟女xx| 欧美黑人巨大hd| 亚洲成人久久性| 亚洲成国产人片在线观看| 国产精品久久久人人做人人爽| 啦啦啦观看免费观看视频高清| 俄罗斯特黄特色一大片| 精品不卡国产一区二区三区| 午夜老司机福利片| 一区二区日韩欧美中文字幕| 窝窝影院91人妻| 午夜亚洲福利在线播放| 夜夜看夜夜爽夜夜摸| 很黄的视频免费| 亚洲片人在线观看| 欧美最黄视频在线播放免费| 男女之事视频高清在线观看| 欧美成人午夜精品| 午夜激情福利司机影院| 欧美黑人精品巨大| 黄色毛片三级朝国网站| 久久久久久久精品吃奶| 成人av一区二区三区在线看| 婷婷丁香在线五月| 长腿黑丝高跟| 99国产精品一区二区三区| 在线观看免费日韩欧美大片| 国产男靠女视频免费网站| svipshipincom国产片| 亚洲精品国产精品久久久不卡| 中文字幕精品免费在线观看视频| 男人的好看免费观看在线视频 | 俄罗斯特黄特色一大片| 女同久久另类99精品国产91| 日韩一卡2卡3卡4卡2021年| 国产av一区二区精品久久| 黄色毛片三级朝国网站| 国产精品久久久久久人妻精品电影| 亚洲精品国产区一区二| 国产v大片淫在线免费观看| 国产精品免费一区二区三区在线| 中文资源天堂在线| 这个男人来自地球电影免费观看| 99国产极品粉嫩在线观看| 日韩欧美在线二视频| 国产一区二区在线av高清观看| 女性被躁到高潮视频| 欧美色视频一区免费| 97人妻精品一区二区三区麻豆 | 午夜福利一区二区在线看| 久久国产亚洲av麻豆专区| 亚洲国产精品成人综合色| 久久久久久人人人人人| 亚洲男人的天堂狠狠| 久久香蕉国产精品| 国产精品99久久99久久久不卡| 18禁黄网站禁片午夜丰满| 国产主播在线观看一区二区| 亚洲精品av麻豆狂野| 免费人成视频x8x8入口观看| 日韩免费av在线播放| 99热这里只有精品一区 | 自线自在国产av| 亚洲国产中文字幕在线视频| av有码第一页| 国产精品98久久久久久宅男小说| 国产成人一区二区三区免费视频网站| 在线永久观看黄色视频| 欧美乱码精品一区二区三区| 午夜免费激情av| 久久国产乱子伦精品免费另类| 亚洲人成网站高清观看| 日本a在线网址| 亚洲成人久久性| 国产亚洲欧美精品永久| 变态另类丝袜制服| 国语自产精品视频在线第100页| 人人妻人人澡人人看| 国产精品久久久人人做人人爽| 亚洲午夜理论影院| 在线观看免费日韩欧美大片| 一级a爱片免费观看的视频| 啦啦啦韩国在线观看视频| 免费搜索国产男女视频| 神马国产精品三级电影在线观看 | 在线观看www视频免费| 久久草成人影院| 亚洲精品粉嫩美女一区| av免费在线观看网站| 免费在线观看成人毛片| 欧美黄色片欧美黄色片| 757午夜福利合集在线观看| 日韩有码中文字幕| 一进一出好大好爽视频| 欧美乱色亚洲激情| 国产真实乱freesex| 国产主播在线观看一区二区| 国产亚洲精品av在线| 给我免费播放毛片高清在线观看| 欧美在线黄色| 亚洲中文日韩欧美视频| 黄片大片在线免费观看| 亚洲第一av免费看| 一区二区三区国产精品乱码| 亚洲av熟女| 亚洲电影在线观看av| 一个人免费在线观看的高清视频| 欧美日韩福利视频一区二区| 久久久久久人人人人人| 99久久精品国产亚洲精品| 97人妻精品一区二区三区麻豆 | 91大片在线观看| 天堂√8在线中文| 国产亚洲精品一区二区www| 欧美成人性av电影在线观看| 精品少妇一区二区三区视频日本电影| 国产精品一区二区三区四区久久 | 97超级碰碰碰精品色视频在线观看| 99国产极品粉嫩在线观看| 免费看日本二区| 国产真人三级小视频在线观看| 午夜久久久久精精品| 美女高潮到喷水免费观看| 一本大道久久a久久精品| 麻豆av在线久日| 亚洲一区二区三区色噜噜| 日本撒尿小便嘘嘘汇集6| 日韩欧美一区视频在线观看| 高潮久久久久久久久久久不卡| 亚洲性夜色夜夜综合| 欧美日韩精品网址| 亚洲色图 男人天堂 中文字幕| 久久欧美精品欧美久久欧美| 成人亚洲精品av一区二区| 欧美午夜高清在线| 麻豆一二三区av精品| 最近最新中文字幕大全免费视频| 成人特级黄色片久久久久久久| 三级毛片av免费|