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

    刮膜式分子蒸餾傳質(zhì)模型及其仿真

    2015-10-13 12:45:24李慧王珂鑫
    化工學(xué)報(bào) 2015年3期
    關(guān)鍵詞:輥筒膜式液膜

    李慧,王珂鑫

    ?

    刮膜式分子蒸餾傳質(zhì)模型及其仿真

    李慧,王珂鑫

    (長春工業(yè)大學(xué)電氣與電子工程學(xué)院,吉林長春 130012)

    分子蒸餾器中加入刮膜器可以增強(qiáng)分子蒸餾過程中的傳質(zhì)效率。以輥筒式刮膜器作為研究對(duì)象,利用進(jìn)料流量、液體濃度、刮膜電機(jī)轉(zhuǎn)速以及蒸發(fā)器尺寸等主要因素之間的結(jié)構(gòu)關(guān)系,建立刮膜式分子蒸餾的數(shù)學(xué)模型。通過將實(shí)驗(yàn)參數(shù)代入到該模型可以快速計(jì)算出分子蒸餾后流出濃度,并將計(jì)算出流出濃度與實(shí)際濃度相比較,誤差范圍縮小到1.5%~2.5%之間,為尋找最佳控制條件提供了參考。將該模型通過MATLAB依次調(diào)整每一個(gè)參數(shù),研究每一個(gè)結(jié)構(gòu)參數(shù)對(duì)分子蒸餾過程的影響,通過研究仿真圖像的變化趨勢(shì),使分子蒸餾器的設(shè)計(jì)可以精確到具體數(shù)值,為刮膜式分子蒸餾器的設(shè)計(jì)提供了依據(jù)。

    分子蒸餾;刮膜式分子蒸餾裝置;輥筒式刮膜器;數(shù)學(xué)模型

    引 言

    分子蒸餾是在高真空或惰性氣體的條件下,利用分子的平均自由程的差別,在低溫(相對(duì)于沸點(diǎn))的狀態(tài)下進(jìn)行的一種高效的液-液分離技術(shù)[1-2],如圖1所示。由于該技術(shù)具有蒸餾溫度低、受熱時(shí)間短、分離程度高、系統(tǒng)能耗低等特點(diǎn),更有利于分離高沸點(diǎn)、高黏度物質(zhì)和熱穩(wěn)定差的有機(jī)混合物。

    刮膜式分子蒸餾器是對(duì)降膜式分子蒸餾器的有效改進(jìn),在蒸餾器內(nèi)部設(shè)置了一個(gè)轉(zhuǎn)動(dòng)的刮膜器,使物料在蒸發(fā)器表面上刮成厚度均勻的液膜,其結(jié)構(gòu)簡(jiǎn)單,分離效果好,得到了廣泛應(yīng)用。雖然刮膜器的加入提高了傳質(zhì)傳熱效率,但也增加了物料流動(dòng)的復(fù)雜性。如果采用安裝傳感器檢測(cè)物料流動(dòng),卻受到傳感器尺寸和刮膜器旋轉(zhuǎn)的影響,很難實(shí)施,而玻璃制成的刮膜式分子蒸餾器為觀察物料流動(dòng)狀況提供了依據(jù)。但是這種實(shí)驗(yàn)器材尺寸有限,蒸餾量較小,在實(shí)際生產(chǎn)應(yīng)用中存在尺度變換的問題,而目前的刮膜式分子蒸餾器的設(shè)計(jì)還停留在依靠工程師的工作經(jīng)驗(yàn)的基礎(chǔ)上。因此研究刮膜器對(duì)液體流動(dòng)的影響狀況,通過各個(gè)參數(shù)間的函數(shù)關(guān)系建立合理的數(shù)學(xué)模型,可以為刮膜式分子蒸餾裝置的設(shè)計(jì)提供參考。

    國外刮膜式分子蒸餾技術(shù)研究起步較早,在研究刮膜器對(duì)液體流動(dòng)影響上取得了很多成果。McKelvey等[3-5]針對(duì)頭波的流體結(jié)構(gòu)進(jìn)行了研究,確定頭波對(duì)刮膜式分子蒸餾過程中分離效率快慢的影響。Micov等[6-7]對(duì)刮膜式分子蒸餾過程的分離因素和液膜組成變化進(jìn)行研究,并建立一維的數(shù)學(xué)模型。Mckenna[8]為了降低刮膜式分子蒸餾器的功耗,討論了各個(gè)參數(shù)對(duì)功耗的影響,為建立刮膜式分子蒸餾器的數(shù)學(xué)模型提供了重要依據(jù)。李國兵等[9]和金晶等[10]先后運(yùn)用Runge-Kutta法和平均隱式差分法對(duì)頭波和液膜變化方程進(jìn)行數(shù)值求解,并且得到很好的穩(wěn)定性,但是在計(jì)算過程中差分網(wǎng)格數(shù)每前進(jìn)一步需要求解4次頭波離散化方程,計(jì)算量較大。許松林等[11-13]利用FLUENT軟件模擬刮膜器對(duì)物料流動(dòng)的影響,來研究蒸發(fā)器內(nèi)部運(yùn)行狀況。FLUENT軟件雖然能形象地描繪出液體流動(dòng)速度矢量,液膜厚度變化,但是分子蒸餾過程中由于分子蒸餾導(dǎo)致物料的濃度變化,必然導(dǎo)致物料流動(dòng)速度和流動(dòng)方向發(fā)生變化,這是FLUENT軟件無法仿真的,所以還需要進(jìn)行數(shù)值分析才能更合理地確定分子蒸餾過程的變化。

    本文根據(jù)前人的研究結(jié)果,通過蒸發(fā)器內(nèi)部尺寸、液體濃度、進(jìn)料量以及其他相關(guān)參數(shù)結(jié)合刮膜器結(jié)構(gòu)特點(diǎn)建立刮膜式分子蒸餾的傳質(zhì)模型。并將實(shí)際實(shí)驗(yàn)參數(shù)代入該傳質(zhì)模型進(jìn)行數(shù)值計(jì)算,與實(shí)際餾出濃度相比較,從而確定模型的可靠性和建模過程中所產(chǎn)生的模型誤差范圍。然后逐一調(diào)整該模型參數(shù),通過餾出濃度的變化曲線分析每個(gè)參數(shù)對(duì)蒸餾過程的影響效果。最后將優(yōu)化后的實(shí)驗(yàn)參數(shù)代入模型仿真與未優(yōu)化的實(shí)驗(yàn)參數(shù)相比較,從而找到更為合適的實(shí)驗(yàn)參數(shù),為分子蒸餾裝置的設(shè)計(jì)和理論研究提供參考。

    1 蒸發(fā)器內(nèi)部結(jié)構(gòu)

    刮膜式分子蒸餾器是在降膜式分子蒸餾器的基礎(chǔ)上改進(jìn)而來,液體除了靠自身重力外還在刮板的作用下不斷混合更新,目前比較常見的4種刮膜器[14]如圖2所示。

    根據(jù)文獻(xiàn)[13]通過FLUENT軟件進(jìn)行的二維仿真結(jié)果分析,圓形轉(zhuǎn)子(即輥筒式)和傾斜刮板(即防滑刮刀和鉸鏈刮刀)比垂直刮板(即鐘擺式)更適用于分子蒸餾過程。而輥筒式所形成的液膜厚度相同且分布均勻不易形成斷流,傳質(zhì)效率更高,所以應(yīng)用范圍更為廣泛。因此選擇輥筒式刮膜器進(jìn)行建模更具實(shí)用性和代表性。

    圖3為刮膜式分子蒸餾器的內(nèi)部結(jié)構(gòu)。旋轉(zhuǎn)軸和安裝在旋轉(zhuǎn)軸上的輥筒共同組成刮膜器,當(dāng)旋轉(zhuǎn)軸轉(zhuǎn)動(dòng)時(shí)輥筒在離心力的作用下緊貼在蒸發(fā)器表面,物料在輥筒與蒸發(fā)面共同的作用下在蒸發(fā)器表面形成一層薄薄的液膜。當(dāng)液體進(jìn)入足夠多時(shí),輥筒前端會(huì)形成頭波。頭波與之前形成的液膜不斷地進(jìn)行混合更新,這種不斷混合更新有利于液膜內(nèi)濃度分布均勻,并且隨著轉(zhuǎn)動(dòng)輥筒自身旋轉(zhuǎn)進(jìn)一步促進(jìn)頭波內(nèi)液體的更新,增加傳熱和傳質(zhì)效果。此外在蒸發(fā)面和冷凝面之間維持在真空或者充入稀有氣體可以使分子更容易從液膜逸出[15-16]。因此刮膜器的引入加快了分子的相對(duì)運(yùn)動(dòng)速度,提高了分子蒸餾的效率。

    由于分子蒸餾是分子運(yùn)動(dòng)的理論應(yīng)用,溫度只是影響分子運(yùn)動(dòng)的快慢,溫度高,液體分子運(yùn)動(dòng)劇烈,黏度降低,流動(dòng)速度加快,分子蒸餾快,反之亦然。液體在向下流動(dòng)過程中除受到重力作用外不受其他外力作用,流動(dòng)速度主要取決于液體自身溫度、黏度和壓強(qiáng)。由于液膜是在輥筒和蒸發(fā)器共同作用下形成的一層極薄的液體,流動(dòng)過程中性質(zhì)不發(fā)生變化,可以根據(jù)Bernoulli方程計(jì)算出液體向下流動(dòng)速度。

    式中,是液體的密度,是液體的流動(dòng)速度,是重力加速度,是液體所處高度,是液體所受壓強(qiáng),是恒值常數(shù)。

    由于流入蒸發(fā)器內(nèi)液體很緩慢,1的速度可以忽略,即1=0,蒸發(fā)器內(nèi)壓強(qiáng)處處相等,即1=2,式(2)可以化簡(jiǎn)為

    式中,D是蒸發(fā)器入料口到輥筒前端的距離,2是液體最初流入輥筒流速。

    刮膜器繞軸進(jìn)行圓周運(yùn)動(dòng),其中刮膜器外緣的水平切向速度t推動(dòng)液體水平流動(dòng),在刮板外緣前端形成頭波

    t=pw(4)

    式中,是刮膜電機(jī)轉(zhuǎn)速,是蒸發(fā)器的直徑(實(shí)際為刮膜器外緣的運(yùn)動(dòng)直徑)。

    在蒸發(fā)面上選擇一指定點(diǎn),以第一次輥筒經(jīng)過該點(diǎn)開始計(jì)時(shí),液膜與頭波進(jìn)行混合更新,之后液膜繼續(xù)進(jìn)行分子蒸餾并且液膜狀態(tài)保持相對(duì)穩(wěn)定,直到下一個(gè)輥筒經(jīng)過該點(diǎn)計(jì)時(shí)結(jié)束,在液膜這段相對(duì)穩(wěn)定的時(shí)間被稱為間隔時(shí)間。

    式中,b是刮膜器上輥筒分支數(shù)。

    同樣根據(jù)Bernoulli方程,液體軸向速度vd,即

    式中,Hs是一個(gè)刮膜器旋轉(zhuǎn)一個(gè)間隔時(shí)間,液體下降高度。

    Bernoulli方程在理想流體中適用,但是實(shí)際流體中是有黏度的,因此,用式(6)計(jì)算vd時(shí)尚需要進(jìn)行修正,即

    式中,是恩氏黏度,恩氏黏度是流出200 ml液體所需時(shí)間與20℃流出同體積蒸餾水所需時(shí)間之比,20℃的蒸餾水是最接近理想流體的。

    刮膜器旋轉(zhuǎn)一個(gè)間隔的時(shí)間所形成的液膜高度等于液體流速與間隔時(shí)間相乘,由于高度變化,液體流速必然發(fā)生變化,但是時(shí)間間隔極短,液體流動(dòng)速度很緩慢,高度變化很小,所以假設(shè)在這段液膜內(nèi)速度變化忽略不計(jì),即

    刮膜器旋轉(zhuǎn)一周形成段的數(shù)量正好等于刮膜器分支數(shù)(即刮板的數(shù)量),并且把同一水平高度的b個(gè)液膜作為一級(jí)。這樣就可以利用刮膜器分支數(shù)計(jì)算出每段液膜在蒸發(fā)器表面上的距離s。

    蒸發(fā)器的總高度,即

    2 模型的建立

    在分子蒸餾過程中,由于刮膜電機(jī)的高速旋轉(zhuǎn),很多傳感器無法安裝到蒸發(fā)器內(nèi),所以如液膜厚度、物料流動(dòng)速度、物料溫度、黏度等參數(shù)無法直接檢測(cè),這都造成通過蒸發(fā)器內(nèi)部的參數(shù)直接建立傳質(zhì)模型的困難性,而蒸發(fā)器外部的參數(shù)由于獲取方便且相對(duì)準(zhǔn)確可靠,如物料進(jìn)入蒸餾器和流出蒸餾器時(shí)的濃度、流量以及刮膜電機(jī)轉(zhuǎn)速等參數(shù),所以為間接推導(dǎo)建立數(shù)學(xué)模型提供了重要的依據(jù)[17]。

    為了進(jìn)一步推導(dǎo)傳質(zhì)模型,假設(shè)每一個(gè)輥筒所過之處頭波和液膜都得到充分混合。以第級(jí)為例,列出該級(jí)蒸發(fā)器表面的質(zhì)量平衡公式為[18]

    式中,Qi為第級(jí)的流入流量,Ci為第級(jí)的流入濃度,Qo為第級(jí)的流出流量,Co為第級(jí)流出濃度,qv為第級(jí)的蒸餾速率。

    在這個(gè)模型中,假設(shè)進(jìn)入第級(jí)的溶液濃度等于從第-1級(jí)頭波流出的濃度。測(cè)量每一級(jí)的流出濃度,目的是為了估算出蒸餾速率。值得注意的是,每一級(jí)的流出濃度是每一段的平均濃度,因?yàn)檎麴s溶液一直處于攪拌當(dāng)中。

    在推導(dǎo)傳質(zhì)方程時(shí),液膜厚度非常薄,與蒸餾器的直徑相比可以近似為一個(gè)平面。這樣可以列出液膜表面蒸餾物在第段質(zhì)量平衡方程[9]

    式中,是質(zhì)量擴(kuò)散系數(shù),大多數(shù)情況下質(zhì)量擴(kuò)散系數(shù)大于靜止?fàn)顟B(tài)下溶液的真正的分子擴(kuò)散系數(shù);是從液膜表面到蒸發(fā)面的距離。

    當(dāng)開始蒸餾的液膜被刮板完全攪拌時(shí),它的濃度正好等于這個(gè)段的流出的濃度。起始條件和邊界條件[19]為

    C=Co,=0 (13)

    由于分子蒸餾的緣故,在液膜表面上液體的濃度為零,而隨著液膜深入液體濃度逐漸增加,在蒸發(fā)面和液體交界處的液體濃度幾乎不發(fā)生變化,為該段的流出濃度。

    式(13)和式(14)中,C雖然只是一個(gè)值,但是卻可以用、和三者組合變量表示為

    將式(15)代回到式(12)整理為

    求解式(16)可得到

    在液膜表面選取一點(diǎn),根據(jù)式(17)可以算出該點(diǎn)的蒸餾流量為

    根據(jù)每個(gè)段的長度定義,當(dāng)刮板以恒定速度繞軸旋轉(zhuǎn)時(shí),式中可以替換為

    式中,是液膜中一個(gè)點(diǎn)沿蒸發(fā)器表面到最近的刮膜器分支水平距離,通過該段的表面積可以算出該段的蒸餾流量

    通過該段的蒸餾流量可以算出整個(gè)第級(jí)的總蒸餾速率

    為了便于計(jì)算,引入一個(gè)參數(shù)用來代替qv,令

    進(jìn)一步假設(shè)在該段的流入流出的體積流速相對(duì)不變,質(zhì)量平衡方程可改寫為

    由此,已知初始濃度0i,就可以迭代推導(dǎo)出相關(guān)的級(jí)的濃度和最終濃度Co

    Komori等[5]在實(shí)驗(yàn)中發(fā)現(xiàn),液體在頭波的交換速率隨著垂直兩個(gè)刮板之間距離不斷增加,當(dāng)達(dá)到1.5~2 cm時(shí)交換速率不發(fā)生改變。同樣在輥筒中也采用豎排小輥筒設(shè)計(jì),所以每個(gè)輥筒之間有一個(gè)間隔,如圖3所示。在刮膜器旋轉(zhuǎn)時(shí),會(huì)在每一個(gè)間隔之間形成一條液帶。液帶的長度與一段液膜的長度相當(dāng),即

    b=s(25)

    液帶的形狀近似于半圓柱,體積大小可用液帶的間隔時(shí)間乘以液帶的流量算出,即

    液帶的高度可以根據(jù)體積計(jì)算出來,即

    此時(shí),參數(shù)可以寫成

    將參數(shù)代入式(24)中,就組成了刮膜式分子蒸餾器的數(shù)學(xué)傳質(zhì)模型。

    3 模型仿真

    式(29)是建立的刮膜式分子蒸餾的數(shù)學(xué)模型,作為數(shù)學(xué)模型在建模的過程中必然忽略一些次要因素,所以需要通過與實(shí)際實(shí)驗(yàn)相結(jié)合的方式來驗(yàn)證模型的可靠性,并求出大概的誤差范圍。故需將多組相同的實(shí)驗(yàn)參數(shù)代入模型計(jì)算仿真,將模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果相比較,求出相對(duì)誤差。并通過模型參數(shù)優(yōu)化后求出的計(jì)算結(jié)果利用相對(duì)誤差來估計(jì)實(shí)際結(jié)果的大概范圍。

    3.1 數(shù)學(xué)模型與實(shí)際結(jié)果比較

    通過HA221-40-11型超臨界CO2萃取裝置分別萃取4份五味子原油。所得量分別為:五味子原油(1) 175.0 g、五味子原油(2) 208.0 g、五味子原油(3) 205.0 g、五味子原油(4) 252.4 g。

    將4份原油分別通過MDS-80型二級(jí)刮膜式分子蒸餾設(shè)備進(jìn)行分子蒸餾,記錄實(shí)驗(yàn)室數(shù)據(jù)如表1。

    表1 五味子原油分子蒸餾實(shí)驗(yàn)數(shù)據(jù)

    ①A drop of oil volume is about 0.05 ml.

    ② Model calculation is needed to convert unit speed.

    其中第一級(jí)為薄膜蒸發(fā)器主要用于進(jìn)行脫水除雜,從第一級(jí)輕餾物流出的是水等雜質(zhì),重餾物作為第二級(jí)的進(jìn)料,直接流入第二級(jí)蒸發(fā)器。第二級(jí)為分子蒸餾器用于提純除酸,其中輕餾物是提純的精油,重餾物是剩余腳子。

    MDS-80型二級(jí)刮膜式分子蒸餾設(shè)備的設(shè)備參數(shù)如表2所示。分別檢驗(yàn)4份原油的溶質(zhì)濃度和第二級(jí)剩余腳子的溶質(zhì)濃度。將實(shí)驗(yàn)數(shù)據(jù)代入數(shù)學(xué)模型中,計(jì)算結(jié)果如表3所示。

    表2 MDS-80型二級(jí)刮膜式分子蒸餾設(shè)備參數(shù)

    表3 計(jì)算結(jié)果與實(shí)際結(jié)果相比較

    在表3中,第1份和第2份原油中進(jìn)料流量相同但刮膜電機(jī)轉(zhuǎn)速不同,轉(zhuǎn)速越慢,重餾物中輕質(zhì)成分的濃度越高。第1份和第3份原油中進(jìn)料流量相同,刮膜電機(jī)也基本相同,重餾物中輕質(zhì)成分的濃度基本相同。第1份和第4份原油中刮膜電機(jī)轉(zhuǎn)速相同但進(jìn)料流量不同,進(jìn)料流量越大,重餾物中輕質(zhì)成分的濃度越高。

    通過實(shí)際重質(zhì)餾出物的濃度與計(jì)算出的餾出濃度相比較,分析模型建立過程忽略的一些因素和假設(shè),可能產(chǎn)生誤差的主要原因如下。

    (1)實(shí)驗(yàn)過程中液體含有黏度,分子蒸餾設(shè)備內(nèi)壁可能殘留液體,由于殘留液體中重餾物含量相對(duì)于輕餾物含量較多,餾出物輕質(zhì)成分檢測(cè)時(shí)濃度高于計(jì)算結(jié)果。

    (2)建模過程中,液體下滑速度是根據(jù)每級(jí)流入點(diǎn)到入料口的垂直距離計(jì)算出來,并假設(shè)該級(jí)內(nèi)速度不發(fā)生變化,因此造成計(jì)算出的流出濃度在局部范圍內(nèi)突然略高于實(shí)際值。

    3.2 蒸餾器內(nèi)部結(jié)構(gòu)對(duì)分子蒸餾過程的影響

    以五味子原油(1)第二級(jí)分子蒸餾器的實(shí)驗(yàn)數(shù)據(jù)為依據(jù),通過依次調(diào)整每個(gè)設(shè)備參數(shù),從而確定每個(gè)參數(shù)對(duì)分子蒸餾過程的影響效果,進(jìn)而為刮膜式分子蒸餾器的設(shè)計(jì)提供參考依據(jù)。

    3.2.1 入料口到輥筒頂端距離 通常設(shè)計(jì)過程中入料口到輥筒頂端會(huì)有一段距離,物料進(jìn)入蒸發(fā)器內(nèi)沿表面下流會(huì)有一定流動(dòng)速度,這個(gè)速度就是液體參與分子蒸餾的初始流速。初始流速對(duì)分子蒸餾的影響如圖4所示。

    由圖4餾出濃度曲線可以看出不同轉(zhuǎn)速的分離程度。當(dāng)曲線趨近于直線時(shí),可以確定出最佳轉(zhuǎn)速區(qū)間。

    仿真結(jié)果表明:入料口到滾筒頂端的距離越大,液體流速隨著不斷下流越來越快,分子蒸餾沒有充分進(jìn)行,液體就已經(jīng)流出蒸餾器,造成分離效率低,輕餾物殘留較多。因此在設(shè)計(jì)分子蒸餾器時(shí)入料口到輥筒頂端的距離盡量縮短,入料時(shí)盡可能減慢流速。

    3.2.2 進(jìn)料流量 進(jìn)料流量的多少?zèng)Q定分子蒸餾器處理速度的快慢,但是進(jìn)料量的多少也受制分子蒸餾分離效果的高低,如圖5所示。

    仿真結(jié)果表明:進(jìn)料量過小,較低的刮膜電機(jī)轉(zhuǎn)速就可以得到很好的分離效果,但是刮膜電機(jī)旋轉(zhuǎn)時(shí),輥筒靠離心力緊貼蒸發(fā)器內(nèi)表面,流量過低就會(huì)造成輥筒與蒸發(fā)器產(chǎn)生摩擦,造成設(shè)備損耗。而進(jìn)料量過大,需要很高的刮膜電機(jī)轉(zhuǎn)速才能達(dá)到理想的分離效果,對(duì)刮膜電機(jī)要求較高。

    3.2.3 刮膜器分支數(shù) 刮膜器分支數(shù)影響間隔時(shí)間的大小,從而影響分子蒸餾的效果,如圖6所示。

    仿真結(jié)果表明:刮膜器分支數(shù)越少,意味著液膜更新時(shí)間越長,在相同轉(zhuǎn)速下,傳質(zhì)效率也隨著刮膜器分支數(shù)的減少而降低。

    3.2.4 刮膜器直徑 刮膜器直徑即蒸發(fā)器直徑,直徑大小影響液膜長度和面積,進(jìn)而影響分子蒸餾的分離效果,如圖7所示。

    仿真結(jié)果表明:刮膜器直徑越大液膜表面越大,相同間隔時(shí)間處理的液體量越大,傳質(zhì)效率高,同時(shí)刮膜電機(jī)的轉(zhuǎn)速調(diào)節(jié)范圍大,對(duì)分離效果具有很好的控制性。

    3.2.5 蒸發(fā)器高度 蒸發(fā)器高度直接決定液膜的級(jí)數(shù),對(duì)分子蒸餾產(chǎn)生影響,如圖8所示。

    仿真結(jié)果表明:蒸發(fā)器的高度越高分子蒸餾時(shí)間越長,傳質(zhì)效果越明顯,分離效果越明顯。但通常體積越大受空間限制,此外反應(yīng)時(shí)間過長也是高度過高所致。

    3.3 分子蒸餾器的參數(shù)優(yōu)化

    以五味子原油(1)為參照,將表1和表2里的五味子原油(1)實(shí)驗(yàn)參數(shù)代入模型中計(jì)算,并根據(jù)以上各個(gè)參數(shù)的變化對(duì)分子蒸餾過程的影響,調(diào)整設(shè)備實(shí)驗(yàn)參數(shù),使分子蒸餾過程達(dá)到最佳效果。

    由于選用的設(shè)備是小型刮膜式分子蒸餾器,本次參數(shù)優(yōu)化根據(jù)原有兩個(gè)蒸發(fā)器中實(shí)驗(yàn)參數(shù)選擇更合適的參數(shù)為依據(jù)。在進(jìn)料流量不變的條件下,蒸發(fā)器直徑較小,所以只將原分支數(shù)由3支調(diào)整為4支。蒸發(fā)器的直徑取最大值0.132 m,蒸發(fā)器高度也取最大值0.444 m,而入料口到刮膜器頂端的距離取最小值0.027 m,進(jìn)行數(shù)值仿真,仿真曲線如圖9所示。

    通過仿真結(jié)果可以看出,參數(shù)優(yōu)化后的餾出濃度明顯下降,并且當(dāng)轉(zhuǎn)速達(dá)到3 r·s-1以上時(shí)濃度已經(jīng)基本不發(fā)生變化。所以最佳的轉(zhuǎn)速為180 r·s-1左右,此時(shí)餾出濃度為249.1 mg·L-1。根據(jù)相對(duì)誤差范圍1.5%~2.5%,如果根據(jù)參數(shù)優(yōu)化后的設(shè)備進(jìn)行分子蒸餾,實(shí)際餾出的濃度為249.1~255.3 mg·L-1。

    4 結(jié) 論

    (1)根據(jù)刮膜式分子蒸餾器的刮膜器結(jié)構(gòu)特點(diǎn),通過進(jìn)料流量、液體濃度、刮膜器轉(zhuǎn)速、蒸發(fā)器尺寸等重要因素建立數(shù)學(xué)模型,并利用MATLAB進(jìn)行仿真與實(shí)際測(cè)量的結(jié)果相比較,發(fā)現(xiàn)模型計(jì)算的輸出濃度低于實(shí)際測(cè)量結(jié)果,誤差范圍是1.5%~2.5%,造成誤差的原因主要是蒸餾器內(nèi)殘留液體沒有計(jì)算在實(shí)際測(cè)量中和模型建立過程中部分參數(shù)被簡(jiǎn)化。

    (2)通過依次調(diào)整傳質(zhì)模型中的參數(shù),確定分子蒸餾器中各個(gè)設(shè)備參數(shù)對(duì)分子蒸餾過程的影響效果。縮短入料口到輥筒頂端的距離,可以減慢進(jìn)料流速,能夠增加液體在蒸餾器內(nèi)時(shí)間,增加傳質(zhì)效率。刮膜器分支數(shù)越多,直徑越大,分子蒸餾的分離效果越好。此外進(jìn)料流量和蒸發(fā)器高度要適中,過大或者過高都會(huì)影響刮膜電機(jī)對(duì)分子蒸餾過程的控制。

    (3)根據(jù)這些參數(shù)的變化規(guī)律,對(duì)原有小型分子蒸餾器的設(shè)備參數(shù)進(jìn)行調(diào)整優(yōu)化,并根據(jù)傳質(zhì)模型計(jì)算出餾出濃度與轉(zhuǎn)速之間的關(guān)系曲線,確定最佳轉(zhuǎn)速,再根據(jù)模型的相對(duì)誤差范圍求出實(shí)際的餾出濃度的范圍。

    符 號(hào) 說 明

    C——常數(shù) Cni——第n級(jí)的液體流入濃度 Cno——第n級(jí)的液體流出濃度 D——蒸發(fā)器直徑,m Fv——液膜表面一點(diǎn)的蒸餾流量,g·m-3·s-1 g——重力加速度,9.8 m·s-2 H——蒸發(fā)器的高度,m Hb——液帶高度,m Hns——第n級(jí)液膜的高度,m h——液體所處高度,m Dh——入料口到輥筒頂端的距離,m Lb——液帶的長度,m Ls——第n段液膜的長度,m Nb——刮膜器分支數(shù) p——液體所受壓強(qiáng),Pa Qni——第n級(jí)的液體流入流量,cm3·s-1 Qno——第n級(jí)的液體流出流量,cm3·s-1 qnv——第n級(jí)的蒸餾速率,g·s-1 ——第n段的蒸餾速率,g·s-1 s——第n段液膜面積,m2 Vb——液帶體積,m3 vnd——第n級(jí)液體垂直速度,m·s-1 vt——刮膜器旋轉(zhuǎn)切速度,m·s-1 dn——計(jì)算參數(shù) h——、和三者組合函數(shù) h-——h的積分 l——擴(kuò)散系數(shù) x——恩氏黏度,°E r——液體密度,g·cm-3 t——間隔時(shí)間,s w——刮膜電機(jī)轉(zhuǎn)速,r·s-1 上角標(biāo) F——最終級(jí)數(shù) n——液膜的級(jí)數(shù) 0——0級(jí) 下角標(biāo) b——液帶 i——流入 o——流出 s——液膜

    References

    [1] Hang Zhijun(杭智軍), Ying Anguo(應(yīng)安國),Wu Hao(武浩).The progress of molecular distillation technology application and research in the pharmaceutical industry [J].(化工生產(chǎn)與技術(shù)), 2012, 19(2): 37-40

    [2] Park S C, Burden D K, Nathanson G M. Surfactant control of gas transport and reactions at the surface of sulfuric acid [J]., 2009, 42(2): 379-387

    [3] James M, McKelvey G, Vincent Sharps Jr. Fluid transport in thin film polymer processors [J]., 1979, 19(9): 652-659

    [4] Komori S, Takata K, Murakami Y. Flow structure and mixing mechanism in an agitated thin-film evaporator [J]., 1988, 21(6): 639-644

    [5] Komori S,Takata K, Murakami Y, Satoru K. Flow and mixing characteristics in an agitated thin film evaporator with vertically aligned blades [J]., 1989, 22(4): 346-350

    [6] Micov M, Lutiˇsan J. Balance equations for molecular distillation [J]., 1997, 32(18): 3051-3066

    [7] Nguyen A D, Goffic F L. Limits of wiped film short-path distillation [J]., 1997, 52(16): 2661-2666

    [8] Mckenna Timothy F. Design model of a wiped film evaporator. Applications to the devolatilisation of polymer melts [J]., 1995(3): 453-467

    [9] Li Guobing(李國兵), Xu Chunjian(許春建), Zhou Ming(周明). Mathematic model of mass transfer and heat transfer of wiped film molecular distillation [J].(天津大學(xué)學(xué)報(bào)), 2003, 36(3): 261-266

    [10] Jin Jing(金晶), Liu Keping(劉克平), Tang Zhonghe(唐重和), Sun Tianfang(孫天放). Mathematic model of the wiped film molecular distillation [J].(長春工業(yè)大學(xué)學(xué)報(bào)), 2012, 33(1): 32-37

    [11] Xu Songlin(許松林), Wang Ling(王玲), Guo Kai(郭凱). Residence time distribution of liquid in wiped film molecular distillation [J].(天津大學(xué)學(xué)報(bào)), 2012, 29(1): 70-74

    [12] Wang Yanfei(王燕飛), Xu Songlin(許松林). Simulation of computational fluid dynamics for molecular distillation process [J].()(化學(xué)工程), 2010, 38(1): 30-33

    [13] Wang Xiaoming(王小明), Xu Songlin(許松林), Xu Chang(徐暢). Study on the effect of wiper structure on the liquid film flow condition in wiped film molecular distillation apparatus [J].(高校化學(xué)工程學(xué)報(bào)), 2012, 26(1): 13-18

    [14] Li Hongzhi(李鴻枝), Jiang Lin(江林), Cao Yaping(曹亞平), Zhu Ruohua(朱若華). Development of molecular distillation equipment and application in sample preparation [J].現(xiàn)代科學(xué)儀器), 2011(2): 61-67

    [15] Qi Chunhua(齊春華), Xing Yulei(邢玉雷), Kang Quan(康權(quán)), Xu Ke(徐克), Feng Houjun(馮厚軍). Pilot test of a 30 t·d-1desalination device with multi-effect and operated at low temperature by evaporating seawater [J].(化工學(xué)報(bào)), 2013, 64(8): 3023-3030

    [16] Xie Guo(謝果), Zheng Hongfei(鄭宏飛), Xiong Jianyin(熊建銀). Natural convection heat transfer and fresh water yield in vertical plate cavity with film evaporation and raw water reservoir [J].(化工學(xué)報(bào)), 2012, 63(8): 2405-2410

    [17] Zhuo Zhen(卓震), Chao Jianwei(巢建偉), Xu Jinwen(許進(jìn)文), Zhang Jing(張晶). Rotating scraped film type molecular distiller [J].(輕工機(jī)械), 2005(4): 134-138

    [18] He Xiaohua(賀小華), Jin Rucong(金如聰), Yang Yuqiang(楊玉強(qiáng)), Tang Ping(唐平). Numerical simulation based on mathematical model of distillation process in thin film evaporator and parameters analysis [J].(食品與機(jī)械), 2009, 25(2): 93-98

    [19] Bird R B, Stewart W E, Lightfoot E N. Transport Phenomena [M]. 2nd ed .Toronto: Wiley, 2002: 114-117

    Mass transfer model and simulation of wiped film molecular distillation

    LI Hui,WANG Kexin

    School of Electrical & Electronic EngineeringChangchun University of TechnologyChangchunJilinChina

    The wiper can enhance mass transfer efficiency of molecular distillation device. Based on the structural relationship between feed flow rate, liquid concentration, motor speed, evaporator size and other factors, a mathematical model of wiped film molecular distillation was established with the roller wiper as research object. By substituting the experimental parameters into the model, the effluent concentrations after molecular distillation could be calculated quickly, and the results were compared with actual concentrations with error rangereduced to 1.5%—2.5%. The model provided a reference for optimal control conditions. The parameters could be adjusted by MATLAB. The influence of each structure parameter on the molecular distillation process was studied, and knowledge of changing tendency of simulation image helped to accurately determine specific detail in the design of molecular distillation device, providing a basis for wiped film molecular distillation device design.

    molecular distillation;wiped film molecular distillation;roller wiper;mathematical model

    2014-06-23.

    WANG Kexin, 370767835@qq.com

    10.11949/j.issn.0438-1157.20140941

    TQ 028.31

    A

    0438—1157(2015)03—1026—09

    國家自然科學(xué)基金項(xiàng)目(61374138);吉林省科技支撐計(jì)劃項(xiàng)目(20130206030GX)。

    2014-06-23收到初稿,2014-12-08收到修改稿。

    聯(lián)系人:王珂鑫。第一作者:李慧(1973—),女,教授。

    supported by the National Natural Science Foundation of China (61374138) and the Jilin Key Technology R&D Program (20130206030GX)。

    猜你喜歡
    輥筒膜式液膜
    影響虹吸式烘干輥筒表面溫度的因素及其解決措施
    考慮軸彎曲的水潤滑軸承液膜建模方法
    煤氣與熱力(2022年2期)2022-03-09 06:29:30
    高空高速氣流下平板液膜流動(dòng)與破裂規(guī)律
    膜式燃?xì)獗碛?jì)量信息采集技術(shù)
    液膜破裂對(duì)PCCS降膜的影響*
    熱塑性塑料在線壓花裝置
    一種改進(jìn)的橡膠機(jī)械設(shè)備中空式冷卻輥筒回水裝置
    流量均勻性對(duì)膜式燃?xì)獗頇z定的影響研究
    膜式燃?xì)獗硎局嫡`差快速檢定裝置設(shè)計(jì)
    两个人免费观看高清视频| 精品卡一卡二卡四卡免费| 久久久久网色| 韩国av在线不卡| a 毛片基地| 成年人午夜在线观看视频| 精品卡一卡二卡四卡免费| 日日爽夜夜爽网站| 婷婷色综合大香蕉| 人成视频在线观看免费观看| 国产精品秋霞免费鲁丝片| 欧美少妇被猛烈插入视频| 国产精品久久久久成人av| 久久久久久久大尺度免费视频| 两个人免费观看高清视频| 在线看a的网站| 婷婷色综合大香蕉| 国产精品久久久久久av不卡| 18禁动态无遮挡网站| 深夜精品福利| 2022亚洲国产成人精品| 精品一区在线观看国产| 一级a爱视频在线免费观看| 高清在线视频一区二区三区| 熟女电影av网| 精品午夜福利在线看| 丰满少妇做爰视频| 久久久a久久爽久久v久久| 成人亚洲精品一区在线观看| 亚洲伊人久久精品综合| 亚洲综合精品二区| 一级爰片在线观看| 中文字幕人妻丝袜制服| 中文精品一卡2卡3卡4更新| 久久精品国产综合久久久| 国产欧美日韩一区二区三区在线| 精品人妻熟女毛片av久久网站| 亚洲精品自拍成人| 国产综合精华液| 亚洲欧美日韩另类电影网站| 亚洲第一青青草原| 欧美人与性动交α欧美精品济南到 | 欧美最新免费一区二区三区| 99热网站在线观看| 高清欧美精品videossex| 亚洲国产精品国产精品| 老司机影院毛片| 亚洲成国产人片在线观看| 啦啦啦在线免费观看视频4| 一区二区三区激情视频| 不卡视频在线观看欧美| 人人妻人人澡人人看| 免费观看性生交大片5| 另类精品久久| 伦理电影大哥的女人| 日韩中字成人| 亚洲,一卡二卡三卡| 久久久久久免费高清国产稀缺| a级片在线免费高清观看视频| 男女边摸边吃奶| 日本猛色少妇xxxxx猛交久久| 免费日韩欧美在线观看| 欧美在线黄色| 满18在线观看网站| www.熟女人妻精品国产| 亚洲一区中文字幕在线| 建设人人有责人人尽责人人享有的| 亚洲国产最新在线播放| 国产免费一区二区三区四区乱码| 在线观看免费高清a一片| 亚洲欧洲国产日韩| 免费观看a级毛片全部| 国产 一区精品| 精品国产一区二区久久| 2018国产大陆天天弄谢| 国产综合精华液| 夜夜骑夜夜射夜夜干| 少妇人妻精品综合一区二区| 亚洲成人一二三区av| 涩涩av久久男人的天堂| 秋霞伦理黄片| 亚洲精品美女久久av网站| 免费高清在线观看日韩| 亚洲人成电影观看| 可以免费在线观看a视频的电影网站 | 99热国产这里只有精品6| 久久99热这里只频精品6学生| 不卡av一区二区三区| 免费看av在线观看网站| 免费大片黄手机在线观看| 少妇被粗大猛烈的视频| 美女福利国产在线| 亚洲欧美日韩另类电影网站| 国产成人午夜福利电影在线观看| 天天操日日干夜夜撸| av福利片在线| 99久久综合免费| 国产在线一区二区三区精| 国产成人精品在线电影| av又黄又爽大尺度在线免费看| 国产精品久久久久成人av| 巨乳人妻的诱惑在线观看| 国产成人精品无人区| 久久久久视频综合| av在线播放精品| 水蜜桃什么品种好| 男女午夜视频在线观看| 午夜福利在线免费观看网站| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久成人av| 欧美97在线视频| 欧美精品人与动牲交sv欧美| 男女高潮啪啪啪动态图| 亚洲国产成人一精品久久久| 国产精品.久久久| 日韩制服骚丝袜av| 久久久久精品性色| 国产精品不卡视频一区二区| 精品国产一区二区久久| 丝袜脚勾引网站| 免费看不卡的av| 母亲3免费完整高清在线观看 | 男人操女人黄网站| 日本色播在线视频| 亚洲国产av新网站| 婷婷色综合大香蕉| 亚洲av男天堂| 一区二区av电影网| 亚洲精品国产色婷婷电影| 黄色毛片三级朝国网站| 国产 精品1| 老熟女久久久| 国产精品二区激情视频| 亚洲国产精品一区二区三区在线| av在线观看视频网站免费| 亚洲一级一片aⅴ在线观看| av国产精品久久久久影院| 国产片特级美女逼逼视频| 午夜精品国产一区二区电影| 国产av一区二区精品久久| 老司机影院毛片| 日韩中字成人| 国产片特级美女逼逼视频| 久久久欧美国产精品| 亚洲色图 男人天堂 中文字幕| 中文字幕av电影在线播放| 亚洲伊人色综图| 男女免费视频国产| 青青草视频在线视频观看| 亚洲美女搞黄在线观看| 夫妻性生交免费视频一级片| 色哟哟·www| 黄网站色视频无遮挡免费观看| av有码第一页| 亚洲欧洲日产国产| 在线亚洲精品国产二区图片欧美| 天天操日日干夜夜撸| 丰满迷人的少妇在线观看| 欧美日韩亚洲高清精品| 中文字幕另类日韩欧美亚洲嫩草| 午夜免费鲁丝| 91成人精品电影| 老鸭窝网址在线观看| 两个人看的免费小视频| 久久人人97超碰香蕉20202| 日本-黄色视频高清免费观看| 欧美国产精品一级二级三级| 国产深夜福利视频在线观看| 午夜日本视频在线| 午夜老司机福利剧场| 精品一品国产午夜福利视频| 国产极品粉嫩免费观看在线| 国产精品久久久久成人av| 欧美老熟妇乱子伦牲交| 一级毛片 在线播放| 欧美日韩精品网址| 亚洲经典国产精华液单| 欧美日韩亚洲国产一区二区在线观看 | 制服丝袜香蕉在线| 国产日韩一区二区三区精品不卡| 久久久久久久大尺度免费视频| 叶爱在线成人免费视频播放| 久久综合国产亚洲精品| 国产成人精品一,二区| av在线观看视频网站免费| 91精品伊人久久大香线蕉| 午夜影院在线不卡| 国产极品粉嫩免费观看在线| 少妇熟女欧美另类| 色94色欧美一区二区| 有码 亚洲区| 亚洲五月色婷婷综合| 国产欧美日韩一区二区三区在线| 老汉色av国产亚洲站长工具| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 日韩一区二区视频免费看| 日本av免费视频播放| 精品国产国语对白av| 秋霞伦理黄片| 色吧在线观看| 另类精品久久| 国产女主播在线喷水免费视频网站| 在线观看国产h片| 亚洲美女搞黄在线观看| 麻豆av在线久日| 免费观看无遮挡的男女| 蜜桃在线观看..| 亚洲国产av新网站| 最近手机中文字幕大全| www.自偷自拍.com| 国产不卡av网站在线观看| 国产男女超爽视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 久久99精品国语久久久| 国产一区二区在线观看av| 国产精品免费大片| 成人18禁高潮啪啪吃奶动态图| 丝袜脚勾引网站| 黄色视频在线播放观看不卡| 黑人欧美特级aaaaaa片| av女优亚洲男人天堂| 亚洲视频免费观看视频| 欧美日韩一级在线毛片| 卡戴珊不雅视频在线播放| 人妻少妇偷人精品九色| 国产精品99久久99久久久不卡 | 国产精品一二三区在线看| 宅男免费午夜| 伊人久久国产一区二区| 日本黄色日本黄色录像| 婷婷色麻豆天堂久久| 欧美另类一区| 97精品久久久久久久久久精品| av在线老鸭窝| 人人妻人人澡人人爽人人夜夜| 69精品国产乱码久久久| 日韩成人av中文字幕在线观看| 久久久久久免费高清国产稀缺| 国产又色又爽无遮挡免| 亚洲精品日韩在线中文字幕| 色哟哟·www| 久久久久久久国产电影| 汤姆久久久久久久影院中文字幕| 欧美人与性动交α欧美软件| 大香蕉久久网| 我的亚洲天堂| 少妇精品久久久久久久| 亚洲国产av新网站| 久久久久久人人人人人| 一区二区三区精品91| 亚洲成人手机| 不卡视频在线观看欧美| 中国国产av一级| 精品亚洲成a人片在线观看| 亚洲精品久久久久久婷婷小说| 人妻人人澡人人爽人人| 色94色欧美一区二区| 中文字幕最新亚洲高清| 久热这里只有精品99| 亚洲av.av天堂| 久久久久久人妻| 夫妻性生交免费视频一级片| 欧美精品国产亚洲| 伊人久久大香线蕉亚洲五| 成人免费观看视频高清| 精品国产乱码久久久久久男人| 亚洲精品国产色婷婷电影| 91成人精品电影| av免费观看日本| 最近2019中文字幕mv第一页| 成年美女黄网站色视频大全免费| 午夜日韩欧美国产| 亚洲精品美女久久av网站| 免费av中文字幕在线| 日本91视频免费播放| 日本午夜av视频| 丝袜脚勾引网站| 人体艺术视频欧美日本| 观看美女的网站| 青春草视频在线免费观看| 精品久久久久久电影网| 如何舔出高潮| 国产 一区精品| 亚洲欧美一区二区三区国产| 免费观看在线日韩| 精品少妇内射三级| 黄片小视频在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 欧美成人精品欧美一级黄| 精品人妻偷拍中文字幕| 免费日韩欧美在线观看| 黄频高清免费视频| 性色avwww在线观看| www.av在线官网国产| 制服人妻中文乱码| 亚洲视频免费观看视频| 欧美另类一区| 久久久精品免费免费高清| 亚洲成国产人片在线观看| 天天躁夜夜躁狠狠躁躁| 免费观看a级毛片全部| 99热全是精品| 丝袜人妻中文字幕| 少妇 在线观看| 亚洲av日韩在线播放| 国产免费一区二区三区四区乱码| 亚洲经典国产精华液单| 在线观看免费高清a一片| 亚洲,欧美,日韩| 欧美中文综合在线视频| 男人添女人高潮全过程视频| 中文精品一卡2卡3卡4更新| 精品人妻一区二区三区麻豆| 日本vs欧美在线观看视频| 欧美日韩亚洲高清精品| 欧美日韩精品成人综合77777| 三级国产精品片| 一级片免费观看大全| 春色校园在线视频观看| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 亚洲男人天堂网一区| 精品一区在线观看国产| 国产又色又爽无遮挡免| 成人影院久久| 国产免费又黄又爽又色| 日韩中字成人| 免费日韩欧美在线观看| 十八禁网站网址无遮挡| 亚洲国产av影院在线观看| 久久国内精品自在自线图片| 天天影视国产精品| 免费观看在线日韩| 国产97色在线日韩免费| 我要看黄色一级片免费的| 久久精品熟女亚洲av麻豆精品| 在线看a的网站| 黑人欧美特级aaaaaa片| 国产精品国产三级国产专区5o| 午夜影院在线不卡| 欧美人与善性xxx| 久久久a久久爽久久v久久| 国产精品av久久久久免费| 亚洲综合色惰| 亚洲色图综合在线观看| 美国免费a级毛片| 宅男免费午夜| 日本vs欧美在线观看视频| 国产人伦9x9x在线观看 | 男女无遮挡免费网站观看| 丝袜美腿诱惑在线| 美女xxoo啪啪120秒动态图| 国产亚洲av片在线观看秒播厂| 美女高潮到喷水免费观看| 久久这里只有精品19| 人妻少妇偷人精品九色| 日韩一区二区视频免费看| 欧美精品一区二区免费开放| 亚洲色图综合在线观看| 精品少妇一区二区三区视频日本电影 | 国产综合精华液| 天天影视国产精品| 成人手机av| 亚洲婷婷狠狠爱综合网| 成人午夜精彩视频在线观看| 欧美精品人与动牲交sv欧美| 久久韩国三级中文字幕| 一级爰片在线观看| 街头女战士在线观看网站| 夜夜骑夜夜射夜夜干| 欧美日韩综合久久久久久| 国产1区2区3区精品| 老汉色av国产亚洲站长工具| 日日摸夜夜添夜夜爱| 国产欧美亚洲国产| 国产1区2区3区精品| 五月开心婷婷网| 久久国产亚洲av麻豆专区| 精品酒店卫生间| 国精品久久久久久国模美| 欧美激情极品国产一区二区三区| 啦啦啦在线免费观看视频4| 99re6热这里在线精品视频| 色视频在线一区二区三区| 99国产精品免费福利视频| 亚洲情色 制服丝袜| kizo精华| 亚洲情色 制服丝袜| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美成人精品欧美一级黄| 观看av在线不卡| 国产欧美日韩一区二区三区在线| 国产极品粉嫩免费观看在线| 97在线人人人人妻| 亚洲一级一片aⅴ在线观看| 免费黄频网站在线观看国产| av片东京热男人的天堂| tube8黄色片| 亚洲av福利一区| 天天影视国产精品| 欧美另类一区| 伦精品一区二区三区| 久久综合国产亚洲精品| 大香蕉久久网| 一区二区三区乱码不卡18| 国产av码专区亚洲av| 欧美 亚洲 国产 日韩一| 久久人妻熟女aⅴ| 中文乱码字字幕精品一区二区三区| 成人毛片60女人毛片免费| 亚洲一区中文字幕在线| 最新的欧美精品一区二区| 青春草国产在线视频| 久久99精品国语久久久| 在线亚洲精品国产二区图片欧美| 日韩制服骚丝袜av| 巨乳人妻的诱惑在线观看| 有码 亚洲区| 午夜福利影视在线免费观看| 18禁国产床啪视频网站| 中文字幕色久视频| 丰满饥渴人妻一区二区三| 日韩熟女老妇一区二区性免费视频| 亚洲欧洲日产国产| 精品亚洲乱码少妇综合久久| 2022亚洲国产成人精品| 久久国内精品自在自线图片| 高清视频免费观看一区二区| 七月丁香在线播放| 中文字幕制服av| 久久久久久久久久久久大奶| 久久97久久精品| 午夜免费观看性视频| 日韩中字成人| 久久久久精品性色| 久久久久国产网址| 欧美另类一区| 久久久a久久爽久久v久久| 日韩中字成人| 免费日韩欧美在线观看| 亚洲国产日韩一区二区| 国产精品 欧美亚洲| 久久精品国产鲁丝片午夜精品| 大香蕉久久成人网| 国产精品熟女久久久久浪| 九草在线视频观看| 亚洲精品国产色婷婷电影| 国产野战对白在线观看| 观看av在线不卡| 曰老女人黄片| 性高湖久久久久久久久免费观看| 最近手机中文字幕大全| 91午夜精品亚洲一区二区三区| 丰满饥渴人妻一区二区三| 久久精品国产自在天天线| 十八禁高潮呻吟视频| 少妇的丰满在线观看| 久久久精品免费免费高清| 人人澡人人妻人| 黄色一级大片看看| 超碰成人久久| 最新的欧美精品一区二区| 如何舔出高潮| 一本—道久久a久久精品蜜桃钙片| 美女中出高潮动态图| 中文字幕最新亚洲高清| 国产一区二区激情短视频 | 久久人妻熟女aⅴ| 狂野欧美激情性bbbbbb| 日本av免费视频播放| 天天躁夜夜躁狠狠躁躁| 777米奇影视久久| 精品视频人人做人人爽| 黑人巨大精品欧美一区二区蜜桃| av免费在线看不卡| 国产av一区二区精品久久| 久久综合国产亚洲精品| 色网站视频免费| 欧美日韩精品网址| 婷婷色综合大香蕉| 一级爰片在线观看| 9色porny在线观看| 人人妻人人澡人人爽人人夜夜| 日本欧美国产在线视频| 一区二区三区精品91| 韩国精品一区二区三区| 国产精品免费大片| 日韩,欧美,国产一区二区三区| 99九九在线精品视频| 男的添女的下面高潮视频| 欧美亚洲日本最大视频资源| 色婷婷av一区二区三区视频| 亚洲人成网站在线观看播放| 少妇猛男粗大的猛烈进出视频| 大话2 男鬼变身卡| 人妻少妇偷人精品九色| av福利片在线| 欧美97在线视频| 久久综合国产亚洲精品| 国产精品av久久久久免费| 一级a爱视频在线免费观看| 中文欧美无线码| av在线app专区| 日本猛色少妇xxxxx猛交久久| 国产探花极品一区二区| 2021少妇久久久久久久久久久| 欧美日韩成人在线一区二区| 久久久久久人妻| 亚洲精品国产一区二区精华液| 亚洲男人天堂网一区| 欧美日韩国产mv在线观看视频| 母亲3免费完整高清在线观看 | 黑人欧美特级aaaaaa片| 亚洲成人av在线免费| 一级a爱视频在线免费观看| 免费在线观看黄色视频的| 免费少妇av软件| 国产综合精华液| 国产有黄有色有爽视频| 黄网站色视频无遮挡免费观看| 国产一区二区在线观看av| 久久久精品免费免费高清| 免费观看无遮挡的男女| 人成视频在线观看免费观看| 欧美日韩亚洲高清精品| 欧美精品一区二区大全| 免费看av在线观看网站| 少妇的丰满在线观看| 亚洲男人天堂网一区| 日韩免费高清中文字幕av| 只有这里有精品99| 97在线人人人人妻| 国产日韩一区二区三区精品不卡| 久久久精品国产亚洲av高清涩受| 国产淫语在线视频| 啦啦啦视频在线资源免费观看| 亚洲精品aⅴ在线观看| 成人二区视频| 国产精品人妻久久久影院| 国产不卡av网站在线观看| 9191精品国产免费久久| 久久99一区二区三区| 国产精品蜜桃在线观看| 精品久久蜜臀av无| 亚洲国产日韩一区二区| 曰老女人黄片| 国产国语露脸激情在线看| 欧美97在线视频| 国精品久久久久久国模美| 亚洲精品aⅴ在线观看| 狠狠精品人妻久久久久久综合| 大片电影免费在线观看免费| 欧美日本中文国产一区发布| 老司机亚洲免费影院| 美女主播在线视频| 日韩精品免费视频一区二区三区| 国产成人免费无遮挡视频| 日韩人妻精品一区2区三区| 国产深夜福利视频在线观看| 国产人伦9x9x在线观看 | 欧美精品亚洲一区二区| 国产麻豆69| 亚洲精品一二三| 婷婷色综合www| 91国产中文字幕| 午夜影院在线不卡| 国精品久久久久久国模美| 伊人亚洲综合成人网| 老熟女久久久| 丰满迷人的少妇在线观看| 天堂中文最新版在线下载| 久热久热在线精品观看| 精品久久蜜臀av无| 国产欧美亚洲国产| 一级a爱视频在线免费观看| 日本vs欧美在线观看视频| 欧美黄色片欧美黄色片| 五月天丁香电影| 99热国产这里只有精品6| 亚洲人成网站在线观看播放| 美国免费a级毛片| 亚洲色图 男人天堂 中文字幕| 日韩电影二区| 99久久综合免费| 在线观看免费日韩欧美大片| 国产xxxxx性猛交| 美女福利国产在线| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久av美女十八| 日本爱情动作片www.在线观看| 欧美精品人与动牲交sv欧美| 亚洲国产精品国产精品| 欧美日韩视频高清一区二区三区二| 国产成人精品久久二区二区91 | 天天躁狠狠躁夜夜躁狠狠躁| 啦啦啦中文免费视频观看日本| 国产精品 国内视频| 涩涩av久久男人的天堂| 久久青草综合色| 一级片免费观看大全| 国产日韩欧美亚洲二区| 久久99热这里只频精品6学生| 男的添女的下面高潮视频| 97在线人人人人妻| 校园人妻丝袜中文字幕| 美女大奶头黄色视频| 亚洲伊人色综图| 欧美日韩视频高清一区二区三区二| 天天操日日干夜夜撸| 女人高潮潮喷娇喘18禁视频| 黄色怎么调成土黄色| 久久久亚洲精品成人影院| 丝瓜视频免费看黄片| 久久ye,这里只有精品|