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

    超聲波霧化強(qiáng)化氨水降膜吸收特性的數(shù)值研究

    2022-10-18 02:45:32何韞玉李舒宏
    關(guān)鍵詞:吸收器降膜霧化器

    何韞玉 李舒宏

    (東南大學(xué)能源與環(huán)境學(xué)院, 南京 210096)

    在國家建立綠色低碳循環(huán)發(fā)展經(jīng)濟(jì)體系,促進(jìn)經(jīng)濟(jì)社會(huì)發(fā)展全面綠色轉(zhuǎn)型的總目標(biāo)下,提升制冷系統(tǒng)性能、降低能耗可以為解決我國資源、環(huán)境、生態(tài)問題作出巨大貢獻(xiàn).目前應(yīng)用最廣泛的制冷系統(tǒng)是壓縮式制冷系統(tǒng),但傳統(tǒng)壓縮式制冷系統(tǒng)中使用的氫氟氯碳化物(HCFCs)和氯氟烴(CFCs)制冷劑對環(huán)境的破壞作用很大,且消耗大量電能.氨吸收式制冷以其直接利用低品位熱源(如太陽能、工業(yè)廢熱余熱)來獲得冷量等優(yōu)勢得到越來越多的關(guān)注[1].

    氨水吸收式制冷技術(shù)已有多年的研究,但目前使用的氨吸收式制冷系統(tǒng)的性能系數(shù)較低.為了充分發(fā)揮氨吸收式制冷系統(tǒng)的優(yōu)勢,擴(kuò)大其應(yīng)用領(lǐng)域,世界各國都在研究和開發(fā)提高其性能的方法.吸收器是影響氨水吸收式制冷機(jī)性能系數(shù)和設(shè)備小型化的重要部件,吸收過程的好壞直接影響整個(gè)制冷循環(huán)的性能.目前大多數(shù)吸收式制冷系統(tǒng)中所采用的吸收器結(jié)構(gòu)為降膜吸收,然而,傳統(tǒng)降膜吸收器的吸收效率不高,出口氨水濃度只能達(dá)到出口飽和狀態(tài)濃度的20%~50%[2].其降膜過程汽液接觸的時(shí)間有限,吸收不夠充分.

    為了提高吸收式制冷系統(tǒng)中吸收器的吸收效率,許多學(xué)者進(jìn)行了相關(guān)研究[3-5].目前提高吸收器性能的方法有3種:對降膜管進(jìn)行表面處理、添加表面活性劑及外加力場[6].表面處理是通過改變降膜管的形狀和結(jié)構(gòu)來增加表面粗糙度,從而加劇流體的擾動(dòng),提高傳熱傳質(zhì)系數(shù)[7-10];許多研究者[11-15]發(fā)現(xiàn)添加少量添加劑(表面活性劑、納米顆粒等)有助于界面湍流,提高傳熱傳質(zhì)系數(shù).此外,相關(guān)領(lǐng)域研究者[16-18]發(fā)現(xiàn)利用電磁場等外加力場可以改善氨-水吸收式的制冷性能.

    超聲波是一種典型的利用外力提高傳熱傳質(zhì)速率的方法,已在溶液除濕、萃取及多相化學(xué)反應(yīng)等領(lǐng)域得到應(yīng)用.陸志等[19]提出在溶液除濕系統(tǒng)中引入超聲波霧化技術(shù)的概念.邊爭等[20]對超聲霧化液體除濕空調(diào)系統(tǒng)的性能進(jìn)行了實(shí)驗(yàn)研究;李曦等[21]研究了超聲霧化液體除濕系統(tǒng)對室內(nèi)空氣品質(zhì)的影響;Yang等[22-25]將超聲波霧化器應(yīng)用于溶液除濕,除濕速率和除濕效率提升可達(dá)31.0%和24.6%.王俐等[26]驗(yàn)證了超聲波霧化技術(shù)的可行性,其理論研究指出超聲波技術(shù)應(yīng)用于溶液除濕系統(tǒng)可以大大降低系統(tǒng)所需的液氣比.Laugier等[27]研究了壓力和超聲波對氣體溶解度和傳質(zhì)的影響,發(fā)現(xiàn)在氣體誘導(dǎo)下,隨著壓力的增加,傳質(zhì)速率增加了11倍.霧化除濕過程本質(zhì)上是一個(gè)絕熱吸收過程.液滴從空氣中捕集水蒸氣的機(jī)理與吸收器中氨水溶液捕集氨蒸汽的機(jī)理相似.因此,可以認(rèn)為超聲霧化與降膜吸收器相結(jié)合可以作為一種提高氨水降膜吸收傳質(zhì)能力的新方法.

    高頻超聲(1.5~2 MHz)對溶液有良好的霧化效果[28],超聲波的存在可以促進(jìn)表面更新效應(yīng),引起界面湍流.超聲霧化產(chǎn)生的微米級(jí)液滴懸浮在吸收降膜管外空間,可以大大增加汽液接觸面積,在有限體積內(nèi)提高吸收器的吸收能力.將降膜溶液部分轉(zhuǎn)化為霧滴,可以在不增大吸收器尺寸或改變降膜管形狀的情況下增加足夠的傳熱傳質(zhì)面積.通過增加液滴的數(shù)量密度和減小液滴的尺寸,可以優(yōu)化吸收器的有效傳熱傳質(zhì)面積.然而,霧化吸收的放熱會(huì)使蒸汽的溫度升高,降低降膜過程的傳質(zhì)勢,接觸面積的增大和傳質(zhì)勢的減小是對立的.因此在降膜吸收器內(nèi)耦合霧化吸收過程,通過降膜管對蒸汽進(jìn)行冷卻,可以改善超聲霧化的不利影響.

    本文通過數(shù)值模擬研究超聲波霧化對降膜吸收性能的影響,并為后續(xù)實(shí)驗(yàn)和垂直管外降膜吸收器的設(shè)計(jì)和優(yōu)化提供一定的理論指導(dǎo).以垂直管外降膜吸收器為研究對象,提出在降膜管外壁上安裝超聲波霧化器的新方法.超聲波霧化器通過與降膜管連接的霧化擋板固定,酒杯形霧化擋板可使剩余未被霧化的氨水溶液沿?fù)醢逋獗诶^續(xù)流向降膜管.本文基于降膜吸收模型和離散相模型建立了新型吸收器中蒸汽側(cè)含有離散相的降膜吸收傳熱傳質(zhì)的數(shù)值模型,從理論上分析了超聲波霧化器對降膜吸收過程的影響.通過定量計(jì)算論證了超聲波霧化對吸收傳質(zhì)效果的增強(qiáng)作用,得到了超聲波霧化器的最佳安裝位置.在此基礎(chǔ)上,進(jìn)一步對比分析了裝有超聲波霧化器的降膜吸收器和傳統(tǒng)降膜吸收器內(nèi)部的濃度場和溫度場分布情況.

    1 物理模型

    1.1 新型降膜吸收器的結(jié)構(gòu)

    為了充分發(fā)揮超聲波霧化器的霧化效果,提出在降膜管外壁上安裝超聲波霧化器的新方法.新型降膜吸收器的結(jié)構(gòu)如圖1所示.吸收器由37根1 m長的不銹鋼降膜管組成,降膜管管長1 m,壁厚1.2 mm,外徑25.4 mm,內(nèi)徑23 mm,縱向管間距和橫向管間距分別為63.5和55 mm,降膜管錯(cuò)排放置.如圖1所示,冷卻水由冷卻水入口流入殼體,經(jīng)過冷卻水分布槽分配均勻,流入各降膜管內(nèi)自下而上流動(dòng),氨蒸汽從吸收器側(cè)壁氨氣入口側(cè)送入吸收器腔體,氨水稀溶液由吸收器頂部入口流入,在布液器作用下形成液膜并沿管壁流動(dòng)形成穩(wěn)定降膜,液膜與吸收器腔體中的氨氣進(jìn)行傳熱和傳質(zhì)過程,在汽液接觸過程中液膜濃度自上而下連續(xù)地增大.

    在每根降膜管的外壁上沿氨水溶液的流向安裝高頻環(huán)形超聲波霧化器.超聲波霧化器由與降膜管相連的霧化擋板固定,擋板外壁呈酒杯弧形,能使得未被霧化的氨水吸收溶液繼續(xù)沿?fù)醢逋獗诹飨蚪的す?并沿降膜管繼續(xù)流動(dòng).超聲波霧化器產(chǎn)生高頻振蕩,將沿降膜管外壁流動(dòng)的溶液霧化成細(xì)小霧滴,霧滴擴(kuò)散到吸收器腔內(nèi)吸收氨蒸汽.吸收器殼體底部設(shè)有濃溶液集液槽,集液槽與濃溶液出口連通.空間中的霧滴和液膜在吸收氨氣后匯集進(jìn)入吸收器底部的溶液集液槽,二者混合形成濃氨水溶液,通過濃溶液出口流出吸收器殼體.

    1—冷卻水出口;2—布液盤;3—霧化擋板;4—濃溶液集液槽;5—冷卻水入口;6—冷卻水分布槽;7—濃溶液出口;8—氨氣入口;9—降膜管;10—超聲波霧化器;11—稀溶液入口;12—集水槽;13—降膜;14—霧滴

    為了計(jì)算超聲波霧化器對吸收性能的影響,在理論制冷量為2.5 kW、發(fā)生溫度為105 ℃、蒸發(fā)溫度為-15 ℃、冷凝溫度為30 ℃的條件下,計(jì)算出吸收器的入口設(shè)計(jì)條件如表1所示.設(shè)計(jì)霧滴直徑ddo為超聲波霧化器頻率f的函數(shù),其計(jì)算公式如下:

    (1)

    式中,ρd為霧滴的密度,kg/m3;σd為液相的表面張力系數(shù),N/m.

    表1 吸收器入口參數(shù)

    1.2 吸收強(qiáng)化效果評價(jià)指標(biāo)

    為了定量表征超聲波霧化的優(yōu)化效果,將是否安裝有超聲波霧化器的吸收器分別命名為新型吸收器和傳統(tǒng)吸收器.為了描述超聲波霧化器對吸收過程的增強(qiáng)效果,定義了以下幾個(gè)評價(jià)指標(biāo).

    對于相同幾何尺寸的吸收器,吸收效果可以通過氨的吸收速率直接表示.單管氨吸收速率N定義如下:

    N=Ms,o+Md,o-Ms,i

    (2)

    式中,Ms,o為吸收器底部液膜的質(zhì)量流量;Md,o為吸收器底部霧滴的質(zhì)量流量;Ms,i為吸收器入口溶液的質(zhì)量流量.

    局部氨吸收速率n定義如下:

    ndL=js,vdAs,v+jd,vdAd,v

    (3)

    式中,js,v為液膜與氨蒸汽之間的局部傳質(zhì)速率;jd,v為霧滴與氨蒸汽之間的局部傳質(zhì)速率;L為長度;As,v表示液膜與氨蒸汽之間的傳質(zhì)面積;Ad,v表示霧滴與氨蒸汽之間的傳質(zhì)面積.

    局部傳質(zhì)速率比β定義為新型吸收器與傳統(tǒng)吸收器在降膜管同一位置的局部氨吸收速率之比, 即

    (4)

    式中,nnew、ntra分別表示新型吸收器和傳統(tǒng)吸收器的局部氨吸收速率.

    2 數(shù)學(xué)模型

    假設(shè)氨水溶液、氨蒸汽和冷卻劑水均勻分布在降膜管中,建立了傳熱傳質(zhì)計(jì)算的數(shù)值模型.通過網(wǎng)格無關(guān)性檢驗(yàn)將降膜管沿長度離散化為1 000個(gè)單元,并采用有限差分法求解數(shù)學(xué)模型.圖2為傳熱傳質(zhì)計(jì)算單元模型.圖中,ms、md、mv、mw分別表示液膜、霧滴、氨蒸汽和冷卻水的質(zhì)量流量;Ts、Td、Tv、Tw分別表示液膜、霧滴、氨蒸汽和冷卻水的

    圖2 降膜管單元控制體示意圖

    單元體溫度;qs,v表示液膜與氨蒸汽之間的局部傳熱速率;qd,v表示霧滴與氨蒸汽之間的局部傳熱速率.降膜溶液在降膜管外壁面呈降膜向下流動(dòng),與氨蒸汽在蒸汽壓差的作用下于汽液交界面處發(fā)生吸收反應(yīng).吸收過程中產(chǎn)生的吸收熱使溶液溫度升高,為不影響吸收效果,該熱量被降膜管內(nèi)的冷卻水帶走,被霧化的霧滴在吸收器腔體吸收氨蒸汽,降膜吸收與霧化吸收同時(shí)進(jìn)行.

    2.1 控制方程

    在建立數(shù)學(xué)模型的過程中,為簡化模型方便求解對模型進(jìn)行如下假設(shè):① 溶液為牛頓型流體,不可壓縮,考慮吸收過程中溶液物性的變化;② 考慮液膜厚度變化;③ 霧滴為標(biāo)準(zhǔn)等徑球體,霧滴內(nèi)部溫度分布均勻;④ 在平行于降膜方向上,液滴的速度不為零,其他方向上的速度為零;⑤ 采用雙膜理論計(jì)算傳質(zhì)過程;⑥ 忽略通過吸收器壁的熱量損失,忽略降膜管和冷卻劑水的熱性能變化;⑦ 氨氣由吸收器側(cè)壁送入吸收器腔體,進(jìn)入的氨氣被完全吸收;⑧ 流動(dòng)是一維的并處于穩(wěn)態(tài).

    2.1.1 質(zhì)量守恒方程

    安裝有超聲波霧化器的單元中,降膜側(cè)的質(zhì)量守恒方程可由下式表示:

    ms,out=(1-a)ms,in

    (5)

    式中,ms,out表示流出該單元體的液膜質(zhì)量流量;ms,in表示流入該單元體的液膜質(zhì)量流量;a為霧化率.

    未安裝有超聲波霧化器的單元降膜側(cè)的質(zhì)量守恒方程由下式表示:

    (6)

    式中,dL為單元體長度.

    式(6)中js,v可以表示為

    (7)

    蒸汽側(cè)的傳質(zhì)系數(shù)Kv可由下式求得:

    (8)

    式中,Hv為蒸汽側(cè)對流換熱系數(shù);Dv表示氨氣擴(kuò)散系數(shù);αv表示氨氣的導(dǎo)熱系數(shù);cp,v表示氨氣的比定壓熱容;ρv表示氨氣的密度.

    降膜溶液側(cè)傳質(zhì)系數(shù)Ks采用Yih[29]提出的一個(gè)適用于較大范圍的降膜吸收器液側(cè)傳質(zhì)系數(shù)方程計(jì)算:

    (9)

    式中,Refilm表示液膜雷諾數(shù);Ds表示液膜擴(kuò)散系數(shù);ρs表示液膜的密度;ms表示液膜的質(zhì)量流量;μs表示液膜的黏度;g表示重力加速度.

    蒸汽側(cè)對流換熱系數(shù)Hv可表示為

    (10)

    式中,Ka表示氨蒸汽的導(dǎo)熱系數(shù).

    安裝有超聲波霧化器的單元中,霧滴側(cè)的質(zhì)量守恒方程可由下式表示:

    md,out=ams,in

    (11)

    式中,md,out表示流出該單元體的霧滴質(zhì)量流量.

    未安裝有超聲波霧化器的單元霧滴的質(zhì)量守恒方程由下式表示:

    (12)

    根據(jù)經(jīng)典的含液滴的多相流理論,將氨蒸汽作為連續(xù)相介質(zhì),液滴作為稀相顆粒處理.離散相模型將液滴視為標(biāo)準(zhǔn)球體,采用集總參數(shù)法處理單個(gè)液滴,忽略了液滴內(nèi)部的濃度和溫度等物理量梯度.

    單元內(nèi)霧滴和蒸汽之間的傳質(zhì)速率jd,v由下式表示:

    (13)

    每個(gè)單元霧滴與氨蒸汽的傳質(zhì)面積dAd,v可由下式計(jì)算得到:

    (14)

    式中,ud為液滴速度.

    每個(gè)單元的霧滴直徑dd可由下式計(jì)算得到:

    (15)

    式中,y表示霧滴的數(shù)量.

    式(13)中霧滴側(cè)傳質(zhì)系數(shù)Kd可表示為

    (16)

    式中,Kd,i為霧滴的內(nèi)部傳質(zhì)系數(shù),可以用Incropera等[30]提出的公式計(jì)算;Kd,o為霧滴的外部傳質(zhì)系數(shù),可用Amokrane等[31]提出的公式計(jì)算;Hms表示亨利系數(shù).在式(14)中,液滴速度ud可以用牛頓第二定律計(jì)算,其計(jì)算如下式所示:

    (17)

    阻力fd由下式計(jì)算:

    (18)

    由于液滴是球形的,故阻力系數(shù)Cd可采用球形阻力系數(shù)來計(jì)算.

    蒸汽側(cè)質(zhì)量守恒方程可表示為

    mvdL=-js,vdAs,v-jd,vdAd,v

    (19)

    吸收器出口溶液由霧滴和降膜側(cè)溶液混合而成,其質(zhì)量守恒方程如下式所示:

    Ms,out+Md,out=Mout

    (20)

    式中,Ms,out表示吸收器底部液膜的質(zhì)量流量;Md,out表示吸收器底部霧滴的質(zhì)量流量;Mout表示吸收器出口溶液質(zhì)量流量.

    2.1.2 組分守恒方程

    在安裝超聲波霧化器的單元和未安裝超聲波霧化器的單元中,降膜側(cè)溶液的組分守恒方程分別如下式所示:

    (21)

    (22)

    式中,ws表示液膜中氨氣的質(zhì)量分?jǐn)?shù).

    霧滴側(cè)的組分守恒方程如下式所示:

    (23)

    式中,wd表示霧滴中氨氣的質(zhì)量分?jǐn)?shù).

    吸收器出口溶液的組分守恒方程由下式表示:

    Ms,outws,out+Md,outwd,out=Moutwout

    (24)

    式中,ws,out表示吸收器底部液膜的氨質(zhì)量分?jǐn)?shù);wd,out表示吸收器底部霧滴的氨質(zhì)量分?jǐn)?shù);wout表示吸收器出口溶液氨質(zhì)量分?jǐn)?shù).

    2.1.3 能量守恒方程

    在安裝超聲波霧化器的單元和未安裝超聲波霧化器的單元中,降膜側(cè)的能量守恒方程分別由下式給出:

    (25)

    (26)

    式中,hs、hv分別表示液膜、氨蒸汽的焓;qs,w表示液膜與冷卻水之間的傳熱速率;As,w表示液膜與冷卻水之間的傳熱面積.

    液膜與氨蒸汽之間的傳熱速率qs,v可以由下式求得:

    qs,v=Hs,v(Tv-Ts)

    (27)

    蒸汽與霧滴之間的換熱系數(shù)Hs,v可以由下式表示:

    (28)

    式(28)中,蒸汽側(cè)對流傳熱系數(shù)Hv由式(10)計(jì)算;液膜側(cè)對流換熱系數(shù)Hs可由下式求得:

    (29)

    式中,ρs、cp,s、αs、Ds分別表示溶液的密度、比定壓熱容、擴(kuò)散系數(shù)和導(dǎo)熱系數(shù).

    qs,w可表示為

    qs,w=Hs,w(Tw-Ts)

    (30)

    冷卻水與降膜側(cè)溶液之間的傳熱系數(shù)Hs,w可根據(jù)下式計(jì)算:

    (31)

    式中,At表示降膜管的平均面積;At,i表示按降膜管內(nèi)徑計(jì)算得到的面積;;dn表示降膜管的外徑;di表示降膜管的內(nèi)徑;kt表示降膜管的導(dǎo)熱系數(shù).降膜側(cè)溶液與降膜管之間的對流換熱系數(shù)Hs,t可根據(jù)Gesellschaft[32]提出的相關(guān)方程計(jì)算,冷卻水與降膜管之間的對流換熱系數(shù)Hw,t可以用Colburn公式計(jì)算.

    霧滴側(cè)和蒸汽側(cè)的能量守恒方程分別由下式表示:

    (32)

    mvhvdL=-(js,vhv+qs,v)dAs,v-(jd,vhv+

    qd,v)dAd,v+qcddAcd

    (33)

    式中,hd表示霧滴的焓;qcd表示相鄰單元間蒸汽的導(dǎo)熱速率;Acd表示相鄰單元間蒸汽的傳質(zhì)面積.

    其中,霧滴與氨蒸汽之間的傳熱速率qd,v可由下式計(jì)算:

    qd,v=Hd,v(Tv-Td)

    (34)

    式中,Hd,v表示霧滴與氨蒸汽之間的傳熱系數(shù).

    努塞爾數(shù)Nud,v由下式表示:

    (35)

    式中,μg、μtp分別表示霧滴和氨蒸汽的黏度系數(shù).

    吸收器出口處的能量守恒方程可表示為

    Ms,ouths,out+Md,outhd,out=Mouthout

    (36)

    式中,hs,out表示吸收器底部液膜的焓;hd,out表示吸收器底部霧滴的焓;hout表示吸收器出口溶液的焓.

    3 模型求解與驗(yàn)證

    3.1 模型求解

    由于模型涉及多元多相流動(dòng),且蒸汽側(cè)存在離散的液滴,其直接求解過程存在較大的計(jì)算量.本文將蒸汽側(cè)和溶液側(cè)通過相界面劃分,降膜溶液側(cè)為二元溶液傳熱傳質(zhì)問題,蒸汽側(cè)為含離散相的傳熱傳質(zhì)問題,二者求解均有成熟的算法.本文建立的數(shù)學(xué)模型考慮了吸收過程中溶液物性的變化,對氨水溶液的物性計(jì)算采用了Schulz氨水狀態(tài)方程,其以吉布斯函數(shù)形式給出,與實(shí)驗(yàn)結(jié)果偏差較小.選取吸收器中的一根降膜管作為計(jì)算對象,對整個(gè)求解區(qū)域進(jìn)行均勻的網(wǎng)格劃分,計(jì)算網(wǎng)格的尺寸通過網(wǎng)格獨(dú)立性校驗(yàn)確定,從而保證計(jì)算精度的同時(shí)節(jié)約計(jì)算資源,沿降膜方向?qū)⒂?jì)算域劃分為1 000個(gè)單元,降膜管內(nèi)部沿法向劃分為100層,對偏微分方程組進(jìn)行離散化,差分方程式中對流項(xiàng)的計(jì)算使用QUICK法.對擴(kuò)散通量,使用中心差分法.使用matlab編程計(jì)算求解非線性方程,入口條件如表2所示,氨吸收速率浮動(dòng)小于10-8kg/s時(shí)停止迭代計(jì)算.

    表2 降膜吸收模型驗(yàn)證結(jié)果

    3.2 模型驗(yàn)證

    安裝有超聲波霧化器的新型吸收器內(nèi)部進(jìn)行的是蒸汽側(cè)含有不飽和霧滴的降膜吸收過程,本文建立的數(shù)學(xué)模型是將傳統(tǒng)降膜吸收模型和霧滴離散相模型通過蒸汽側(cè)進(jìn)行耦合,在超聲波霧化器安裝高度以上進(jìn)行的是傳統(tǒng)降膜吸收過程,在超聲波霧化器安裝高度之下同時(shí)進(jìn)行降膜側(cè)溶液傳熱傳質(zhì)和霧滴傳熱傳質(zhì)過程.降膜吸收已有成熟的模型可以描述吸收過程的傳熱傳質(zhì),在離散相模型的基礎(chǔ)上建立了霧化吸收的數(shù)值模型,液滴速度、溫度的概率密度函數(shù)是由動(dòng)力學(xué)方程確定的,離散相系統(tǒng)的平均方程也是基于動(dòng)力學(xué)方程.將汽液之間的傳熱傳質(zhì)作為汽液控制方程的源相,這是計(jì)算液滴-流體相互作用的一種常用方法.當(dāng)霧滴進(jìn)入蒸汽時(shí),將蒸汽視為連續(xù)相介質(zhì),霧滴視為稀相顆粒.根據(jù)含液滴的經(jīng)典多相流理論,將液滴視為標(biāo)準(zhǔn)球體,采用集總參數(shù)法處理單個(gè)液滴,忽略了液滴內(nèi)部的濃度和溫度等物理量梯度.利用文獻(xiàn)[33]的實(shí)驗(yàn)數(shù)據(jù)對本文的降膜吸收模型進(jìn)行了分析驗(yàn)證,驗(yàn)證結(jié)果如表2所示.在相同的入口條件下本文模型的模擬結(jié)果與文獻(xiàn)[33]的實(shí)驗(yàn)結(jié)果吻合良好.本文在建立數(shù)學(xué)模型和求解的過程中盡量減少模擬和實(shí)際的偏差,考慮了氨水溶液物性的變化和液膜厚度變化等因素.

    4 結(jié)果與討論

    圖3給出了超聲波霧化器安裝位置對單管氨吸收率N和吸收器出口溶液氨質(zhì)量分?jǐn)?shù)wout的影響.當(dāng)超聲波霧化器安裝在距降膜管頂部1 m處時(shí),吸收過程相當(dāng)于傳統(tǒng)降膜吸收過程.從圖3可以看出,隨著超聲波霧化器安裝高度的降低,單管氨吸收率N和吸收器出口溶液氨質(zhì)量分?jǐn)?shù)先升高而后降低,在降膜管從上至下0.45 m處達(dá)到最大值,此時(shí)單管總氨吸收率N為0.028 1 g/s,出口溶液氨質(zhì)量分?jǐn)?shù)為37.5%.傳統(tǒng)吸收器中總氨吸收率N和出口溶液氨質(zhì)量分?jǐn)?shù)分別為0.024 4 g/s和36.3%.與傳統(tǒng)吸收器相比,在最佳位置安裝超聲波霧化器的新型吸收器氨吸收率提高了15.1%.

    從圖3中可看出,超聲波霧化器對吸收效果的提升程度隨著其安裝位置自上向下移動(dòng)呈現(xiàn)先增大后減小的趨勢.若超聲波霧化器安裝在降膜管靠上部分,則霧滴的運(yùn)動(dòng)行程和停留時(shí)間變長,霧滴還未落于底部就已達(dá)到飽和狀態(tài),當(dāng)霧滴飽和后其溫度高于氨蒸汽的溫度,故在霧滴吸收達(dá)到飽和后直至落于底部前的這段行程里,霧滴將不斷地向氨蒸汽釋放熱量,此部分熱量來源于霧滴與氨蒸汽的溫差,并非霧滴絕熱吸收氨蒸汽釋放的熱量.故若超聲波霧化器安裝位置靠上,霧滴行程變長,其向氨蒸汽釋放的熱量會(huì)增大而吸收氨蒸汽的質(zhì)量不會(huì)增大.同時(shí),其對氨蒸汽的升溫影響也會(huì)降低降膜過程的傳質(zhì)推動(dòng)力,故降膜管中總氨吸收速率會(huì)減小.若超聲波霧化器安裝在降膜管靠下位置,則霧滴在吸收器腔體內(nèi)的運(yùn)動(dòng)行程和停留時(shí)間將會(huì)縮短,霧滴還未充分吸收氨氣就落于吸收器底部,其傳質(zhì)過程受到限制,故而導(dǎo)致霧滴吸收氨氣量減少.

    圖3 N和wout隨不同超聲波霧化器位置的變化關(guān)系

    圖4(a)表示安裝有超聲波霧化器的新型吸收器與傳統(tǒng)吸收器的局部傳質(zhì)率比β沿吸收器高度的變化趨勢.冷卻水在降膜管內(nèi)自下而上流動(dòng),氨水溶液沿降膜管外壁自上而下流動(dòng).從圖中可以看出在超聲波霧化器安裝位置,局部傳質(zhì)速率比驟增,新型吸收器內(nèi)局部氨吸收速率顯著提高,這是由于在吸收過程中霧滴產(chǎn)生的瞬間,汽液接觸面積顯著增大,從而局部吸收傳質(zhì)速率迅速提高.圖4(b)為局部氨吸收速率n沿吸收器高度的變化趨勢,在超聲波霧化器安裝位置,由于降膜側(cè)溶液部分被霧化,液膜質(zhì)量流量急劇下降,導(dǎo)致液膜局部氨吸收速率驟降.同時(shí),液膜厚度的減薄會(huì)使冷卻水與液膜之間的傳熱熱阻減小,使得液膜吸收氨氣的熱量能迅速傳遞給冷卻水,液膜溫度迅速下降,雖然液膜溫度隨吸收進(jìn)行下降而液膜濃度隨吸收進(jìn)行升高,但此時(shí)溫度下降對傳質(zhì)推動(dòng)力的影響大于濃度升高對其的影響,故液膜局部吸收速率總體增加,其在持續(xù)下降之前有一個(gè)短暫的回升.隨著液膜沿管壁向下流動(dòng)并逐漸趨于飽和,濃度升高對其的影響逐漸大于溫度下降對傳質(zhì)推動(dòng)力的影響,其傳質(zhì)推動(dòng)力不斷下降,局部氨吸收速率逐漸降低.由于霧化吸收類似于絕熱過程,隨著吸收進(jìn)行其溫度和濃度不斷升高,吸收推動(dòng)勢逐漸減小,故其局部氨吸收速率沿吸收器高度下降逐漸減小.

    (a) 局部傳質(zhì)速率比β

    (b)局部傳質(zhì)速率n

    圖5給出了安裝有超聲波霧化器的新型吸收器和傳統(tǒng)吸收器中沿降膜管高度的氨質(zhì)量分?jǐn)?shù)分布.wnew和wtra分別表示新型吸收器和傳統(tǒng)吸收器中局部氨質(zhì)量分?jǐn)?shù).新型吸收器中局部氨質(zhì)量分?jǐn)?shù)wnew可以表示為

    (37)

    式中,wd,new為局部霧滴氨質(zhì)量分?jǐn)?shù);ws,new為局部降膜溶液氨質(zhì)量分?jǐn)?shù).

    圖5 氨質(zhì)量分?jǐn)?shù)隨降膜管長度的變化曲線

    從圖5中可以觀察到,傳統(tǒng)吸收器出口溶液氨質(zhì)量分?jǐn)?shù)為36.3%;新型吸收器中底部降膜側(cè)溶液氨質(zhì)量分?jǐn)?shù)為38.6%,底部霧滴氨質(zhì)量分?jǐn)?shù)為35.3%,二者混合后出口溶液氨質(zhì)量分?jǐn)?shù)為37.5%.霧滴產(chǎn)生后,wnew相較于wtra增加得更快,這充分證明了霧化吸收過程顯著增加了汽液傳質(zhì)面積,提高了吸收傳質(zhì)效果.由于霧滴產(chǎn)生時(shí)降膜側(cè)液膜的質(zhì)量流量和厚度減小,因此液膜吸收氨氣速率顯著增加,圖5中降膜管長度大于0.45 m后液膜氨質(zhì)量分?jǐn)?shù)曲線斜率增大.此外,由于降膜側(cè)溶液吸收氨氣放出的熱量由冷卻水帶走,而霧滴吸收為絕熱吸收,故霧滴吸收氨氣速率沿降膜管下降得更快,霧滴氨質(zhì)量分?jǐn)?shù)曲線的增長率相較于液膜氨質(zhì)量分?jǐn)?shù)曲線的增長率下降得更快.

    超聲波作用于溶液會(huì)產(chǎn)生熱效應(yīng).圖6表示新型吸收器和傳統(tǒng)吸收器中液膜、霧滴、氨蒸汽的溫度沿吸收器高度的變化曲線.新型吸收器出口溶液溫度為27.04 ℃,而傳統(tǒng)吸收器出口溶液溫度為24.41 ℃,新型吸收器中的出口溶液溫度相較于傳統(tǒng)吸收器提高了2.6 ℃.由此可見霧化過程對出口溶液的溫升影響較小;此外,新型吸收器中底部氨蒸汽溫度比傳統(tǒng)吸收器低0.8 ℃.新型吸收器冷卻水出口溫度比傳統(tǒng)吸收器高0.2 ℃.由圖6可知,與傳統(tǒng)吸收器相比,新型吸收器中霧滴產(chǎn)生后液膜的溫度下降速度更快,這是由于液膜厚度的減薄,冷卻水與液膜之間的傳熱增強(qiáng).由于霧化吸收將熱量釋放到氨蒸汽中,導(dǎo)致氨蒸汽溫度在霧滴產(chǎn)生的瞬間驟升.

    但是,由于20 ℃的氨蒸汽是由吸收器側(cè)壁送入吸收器腔體的,并非由頂部送入,因此新型吸收器中蒸汽的溫度變化趨勢與霧滴溫度變化趨勢不一致,側(cè)送入的蒸汽同時(shí)受到液膜和霧滴的加熱作用,由于液膜的溫度沿降膜管逐漸下降,新型吸收器中降膜管下部液膜溫度低于氨蒸汽溫度,液膜可以冷卻蒸汽,霧滴加熱蒸汽,且液膜溫度下降的速率大于霧滴溫度上升速率,故液膜和霧滴釋放給蒸汽的總熱量隨高度下降而減小,靠近降膜管底部的液膜和霧滴釋放給氨蒸汽的熱量相比于上部減小了,氨蒸汽溫度沿管高度呈現(xiàn)下降的趨勢;由于冷卻水的質(zhì)量流量大,故超聲波霧化對冷卻水溫度的溫升影響較小,新型吸收器和傳統(tǒng)吸收器中冷卻水出口溫度分別為22.4和22.2 ℃.綜合來看,超聲波霧化過程對吸收器內(nèi)部的溫度升高影響較小,霧化過程對降膜吸收過程的不利溫升影響與其強(qiáng)化吸收的效果相比可以忽略不計(jì).

    5 結(jié)論

    1) 本文通過數(shù)值模擬證明了超聲波霧化器對降膜吸收性能的提升效果.在本文設(shè)計(jì)參數(shù)下,超聲波霧化器的最佳安裝高度為降膜管從上至下0.45 m處.隨著超聲波霧化器安裝高度的下降,超聲波霧化對降膜吸收的提升效果先增大后減小.

    2) 當(dāng)超聲波霧化器安裝在最佳位置,其他工況為設(shè)計(jì)工況時(shí),單管氨吸收速率N提升了15.1%,吸收器出口吸收溶液氨質(zhì)量分?jǐn)?shù)由36.3%提高到37.5%.

    3) 超聲波霧化過程產(chǎn)生的熱效應(yīng)對降膜吸收的溫升影響較小,與其強(qiáng)化吸收的效果相比可以忽略不計(jì).與傳統(tǒng)吸收器相比,安裝有超聲波霧化器的新型吸收器出口溶液溫度提高了2.6 ℃;新型吸收器中底部氨蒸汽溫度比傳統(tǒng)吸收器中低0.8 ℃,其冷卻水出口溫度比傳統(tǒng)吸收器高0.2 ℃.

    猜你喜歡
    吸收器降膜霧化器
    腔式寬光譜高吸收比標(biāo)準(zhǔn)器研制及測試分析
    霧化器在聚氯乙烯糊樹脂生產(chǎn)中的應(yīng)用
    中國氯堿(2021年10期)2021-12-21 06:18:08
    燒堿裝置降膜蒸發(fā)器的腐蝕檢查
    霧化器晾干裝置的制作與使用
    波浪能點(diǎn)吸收器結(jié)構(gòu)設(shè)計(jì)與數(shù)值優(yōu)化
    能源吸收器
    基于CST的紅外吸收器特性分析
    板式降膜蒸發(fā)器激光焊蜂窩狀波紋傳熱板
    中國造紙(2015年7期)2015-12-16 12:40:48
    水平管外R404A降膜蒸發(fā)傳熱的實(shí)驗(yàn)研究
    基于ARM9的低頻超聲霧化器驅(qū)動(dòng)電源優(yōu)化設(shè)計(jì)及仿真
    国产一区亚洲一区在线观看| 亚洲久久久久久中文字幕| 青青草视频在线视频观看| 爱豆传媒免费全集在线观看| 亚洲综合色惰| 国产精品福利在线免费观看| 看黄色毛片网站| 国产亚洲一区二区精品| 午夜福利在线观看吧| 美女内射精品一级片tv| av一本久久久久| 2018国产大陆天天弄谢| 国产精品一区二区三区四区久久| 99热网站在线观看| 国产女主播在线喷水免费视频网站 | 99久国产av精品国产电影| 日韩成人av中文字幕在线观看| 午夜激情久久久久久久| 哪个播放器可以免费观看大片| 国产精品国产三级专区第一集| 99re6热这里在线精品视频| 最近最新中文字幕大全电影3| 99久国产av精品| 亚洲欧美清纯卡通| 国产精品三级大全| 看非洲黑人一级黄片| 久久久久久久大尺度免费视频| 欧美激情在线99| 亚洲av电影在线观看一区二区三区 | 国产探花在线观看一区二区| 久热久热在线精品观看| 国内精品宾馆在线| 91在线精品国自产拍蜜月| 亚洲最大成人av| 国产精品三级大全| 爱豆传媒免费全集在线观看| 美女主播在线视频| 日韩人妻高清精品专区| 高清午夜精品一区二区三区| 日本熟妇午夜| 亚洲欧美成人精品一区二区| 日本一本二区三区精品| 亚洲精品国产av蜜桃| 免费在线观看成人毛片| 精品一区二区三区视频在线| 欧美xxxx黑人xx丫x性爽| 欧美性感艳星| 九色成人免费人妻av| 国产 亚洲一区二区三区 | 自拍偷自拍亚洲精品老妇| 麻豆久久精品国产亚洲av| 三级经典国产精品| 高清视频免费观看一区二区 | 国产伦精品一区二区三区视频9| 亚洲熟女精品中文字幕| 午夜日本视频在线| av在线蜜桃| 欧美最新免费一区二区三区| 国产午夜精品论理片| 日韩视频在线欧美| 伊人久久国产一区二区| 国产一级毛片在线| 三级毛片av免费| 久久国内精品自在自线图片| 卡戴珊不雅视频在线播放| 免费大片18禁| 舔av片在线| 国产黄a三级三级三级人| 国产成人免费观看mmmm| 亚洲欧洲国产日韩| 亚洲熟女精品中文字幕| 男女边摸边吃奶| 大话2 男鬼变身卡| 国产老妇女一区| 色综合色国产| 亚洲精华国产精华液的使用体验| 成人综合一区亚洲| 亚洲人成网站高清观看| 蜜桃久久精品国产亚洲av| 欧美丝袜亚洲另类| 欧美xxxx性猛交bbbb| 中文乱码字字幕精品一区二区三区 | 麻豆精品久久久久久蜜桃| 久久精品夜色国产| 亚洲精品日本国产第一区| 精品久久久精品久久久| 免费观看性生交大片5| 偷拍熟女少妇极品色| 蜜臀久久99精品久久宅男| 欧美不卡视频在线免费观看| 老师上课跳d突然被开到最大视频| 免费黄网站久久成人精品| 久久草成人影院| 美女高潮的动态| 99热这里只有是精品在线观看| 欧美xxxx性猛交bbbb| 51国产日韩欧美| 成年人午夜在线观看视频 | 丝瓜视频免费看黄片| 精品久久久久久成人av| 超碰97精品在线观看| 久久精品国产亚洲网站| 国产黄色视频一区二区在线观看| 老司机影院成人| 国产精品嫩草影院av在线观看| 最近手机中文字幕大全| 午夜老司机福利剧场| 性插视频无遮挡在线免费观看| 卡戴珊不雅视频在线播放| 亚洲婷婷狠狠爱综合网| 91精品伊人久久大香线蕉| 日韩 亚洲 欧美在线| 51国产日韩欧美| 国产在线一区二区三区精| 国产探花极品一区二区| 人妻制服诱惑在线中文字幕| 狠狠精品人妻久久久久久综合| 美女被艹到高潮喷水动态| 国产伦理片在线播放av一区| 美女脱内裤让男人舔精品视频| 国产 一区精品| h日本视频在线播放| 最近手机中文字幕大全| 一区二区三区乱码不卡18| 亚洲精品影视一区二区三区av| 老司机影院成人| 91精品伊人久久大香线蕉| 大香蕉97超碰在线| 97精品久久久久久久久久精品| 亚洲精品第二区| 国产美女午夜福利| 青春草国产在线视频| 又爽又黄无遮挡网站| 如何舔出高潮| 99热这里只有精品一区| 亚洲精品乱码久久久v下载方式| 欧美日韩亚洲高清精品| 熟女人妻精品中文字幕| 国产精品嫩草影院av在线观看| 亚洲国产色片| 最近2019中文字幕mv第一页| 日韩一区二区三区影片| 91精品国产九色| 国产91av在线免费观看| 乱码一卡2卡4卡精品| 国产精品伦人一区二区| 99热这里只有精品一区| 天堂影院成人在线观看| 精品久久久噜噜| 精品久久久噜噜| 亚洲人成网站在线观看播放| av在线亚洲专区| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区免费观看| 黄色日韩在线| 国产男女超爽视频在线观看| 国产毛片a区久久久久| 白带黄色成豆腐渣| 人体艺术视频欧美日本| 天堂中文最新版在线下载 | 亚洲国产精品成人综合色| 亚洲四区av| 又爽又黄无遮挡网站| 午夜福利成人在线免费观看| 亚洲av在线观看美女高潮| 伊人久久国产一区二区| 成人综合一区亚洲| 免费看不卡的av| 成人高潮视频无遮挡免费网站| av天堂中文字幕网| 大陆偷拍与自拍| 免费观看av网站的网址| av网站免费在线观看视频 | 七月丁香在线播放| 九九久久精品国产亚洲av麻豆| 国产精品国产三级国产av玫瑰| 午夜激情久久久久久久| 午夜激情久久久久久久| 免费观看无遮挡的男女| 人妻少妇偷人精品九色| 婷婷色av中文字幕| 久久久久久久亚洲中文字幕| 午夜免费观看性视频| 精品酒店卫生间| .国产精品久久| 美女脱内裤让男人舔精品视频| 国产有黄有色有爽视频| 亚洲av中文av极速乱| 精品国产三级普通话版| 国产精品国产三级国产av玫瑰| 久久人人爽人人爽人人片va| 最近视频中文字幕2019在线8| 久久久久久久午夜电影| 嫩草影院新地址| 男人和女人高潮做爰伦理| 特大巨黑吊av在线直播| 天堂俺去俺来也www色官网 | 美女被艹到高潮喷水动态| 亚洲欧美成人精品一区二区| 十八禁网站网址无遮挡 | 日产精品乱码卡一卡2卡三| 亚洲av成人精品一二三区| 亚洲熟妇中文字幕五十中出| 一级二级三级毛片免费看| 少妇裸体淫交视频免费看高清| 中国国产av一级| 嫩草影院入口| 人妻少妇偷人精品九色| 一级av片app| 日韩电影二区| 三级国产精品欧美在线观看| 七月丁香在线播放| 毛片一级片免费看久久久久| 亚洲国产成人一精品久久久| 亚洲欧美日韩东京热| 国产精品人妻久久久久久| 少妇的逼好多水| 亚洲精品色激情综合| 国国产精品蜜臀av免费| 观看美女的网站| 亚洲精品影视一区二区三区av| 国产老妇伦熟女老妇高清| 成人鲁丝片一二三区免费| 欧美高清性xxxxhd video| 大又大粗又爽又黄少妇毛片口| 午夜免费男女啪啪视频观看| 久久久久国产网址| 国产免费福利视频在线观看| 日韩人妻高清精品专区| 免费观看a级毛片全部| 精品一区二区免费观看| 欧美精品国产亚洲| 乱人视频在线观看| 六月丁香七月| 水蜜桃什么品种好| 久久草成人影院| 黄色欧美视频在线观看| 国产成人一区二区在线| 亚洲在久久综合| 三级经典国产精品| 国内少妇人妻偷人精品xxx网站| 国产黄色免费在线视频| 自拍偷自拍亚洲精品老妇| 精品一区二区免费观看| 国产亚洲最大av| 欧美日韩国产mv在线观看视频 | 国产精品人妻久久久久久| 欧美高清成人免费视频www| 国产毛片a区久久久久| 天堂中文最新版在线下载 | 丝袜喷水一区| 97超碰精品成人国产| 波野结衣二区三区在线| 亚洲精品aⅴ在线观看| 又大又黄又爽视频免费| a级毛片免费高清观看在线播放| 美女内射精品一级片tv| 少妇人妻精品综合一区二区| av线在线观看网站| 欧美日本视频| 国产午夜福利久久久久久| 亚洲国产精品sss在线观看| 777米奇影视久久| 男女视频在线观看网站免费| 久久精品夜色国产| 国产乱来视频区| 日韩成人伦理影院| av天堂中文字幕网| 国产一区二区在线观看日韩| 亚洲精品久久午夜乱码| 青春草国产在线视频| 成人二区视频| 久久久欧美国产精品| 国产伦在线观看视频一区| 欧美最新免费一区二区三区| 免费高清在线观看视频在线观看| 午夜激情久久久久久久| 一级毛片我不卡| 熟妇人妻不卡中文字幕| 国产在线男女| 中文天堂在线官网| 国产成人免费观看mmmm| 午夜激情久久久久久久| 一本久久精品| 欧美3d第一页| 亚洲av中文字字幕乱码综合| 国产女主播在线喷水免费视频网站 | 国产成年人精品一区二区| 大片免费播放器 马上看| 免费看日本二区| 国产精品99久久久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲av福利一区| 欧美成人a在线观看| 成年av动漫网址| 日韩视频在线欧美| 国产精品一区二区在线观看99 | 我的女老师完整版在线观看| 69人妻影院| 午夜精品国产一区二区电影 | 欧美3d第一页| 国产毛片a区久久久久| 精品国产三级普通话版| 毛片一级片免费看久久久久| 九九爱精品视频在线观看| 国产白丝娇喘喷水9色精品| 日本午夜av视频| av卡一久久| 午夜激情欧美在线| 中文在线观看免费www的网站| 国产毛片a区久久久久| 国产午夜精品久久久久久一区二区三区| 久久久久久国产a免费观看| 最近手机中文字幕大全| 亚洲人成网站在线播| 欧美区成人在线视频| 亚洲av免费在线观看| av卡一久久| 成年av动漫网址| 国产伦一二天堂av在线观看| 男人舔奶头视频| 精品久久国产蜜桃| 午夜激情福利司机影院| 日韩视频在线欧美| 亚洲成人精品中文字幕电影| 亚洲av成人av| 97热精品久久久久久| 久久久久久久国产电影| 99视频精品全部免费 在线| 看免费成人av毛片| 亚洲欧美精品专区久久| 搡老妇女老女人老熟妇| 欧美极品一区二区三区四区| 精品不卡国产一区二区三区| 国产 一区精品| 深爱激情五月婷婷| 国产精品蜜桃在线观看| 欧美激情久久久久久爽电影| 色综合亚洲欧美另类图片| 老司机影院成人| 欧美一级a爱片免费观看看| 国产精品福利在线免费观看| 亚洲综合精品二区| 欧美潮喷喷水| 国产精品福利在线免费观看| 六月丁香七月| 日韩电影二区| 国产精品无大码| 大片免费播放器 马上看| 亚洲av成人精品一区久久| 一区二区三区高清视频在线| 青春草亚洲视频在线观看| 亚洲av不卡在线观看| 精品熟女少妇av免费看| 国产亚洲午夜精品一区二区久久 | 国产精品三级大全| 人妻少妇偷人精品九色| av专区在线播放| 亚洲激情五月婷婷啪啪| 久久精品久久精品一区二区三区| 国产精品伦人一区二区| 国产探花极品一区二区| 91aial.com中文字幕在线观看| 亚洲国产精品成人综合色| 亚洲av成人精品一区久久| 亚洲精品aⅴ在线观看| 乱人视频在线观看| 亚洲av不卡在线观看| 国产精品.久久久| 日韩伦理黄色片| 丝袜美腿在线中文| 日本三级黄在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲av.av天堂| 午夜福利视频1000在线观看| 真实男女啪啪啪动态图| 亚洲电影在线观看av| 久久人人爽人人爽人人片va| 高清欧美精品videossex| 丝袜喷水一区| av在线蜜桃| 日本免费a在线| 日本三级黄在线观看| 永久免费av网站大全| 亚洲欧洲日产国产| 在线免费观看不下载黄p国产| 日韩在线高清观看一区二区三区| 国产午夜精品久久久久久一区二区三区| 色哟哟·www| 亚洲精品成人av观看孕妇| 80岁老熟妇乱子伦牲交| 久久久午夜欧美精品| 亚洲av一区综合| 日本三级黄在线观看| 肉色欧美久久久久久久蜜桃 | 搡女人真爽免费视频火全软件| 欧美日韩国产mv在线观看视频 | 国产不卡一卡二| 午夜福利在线观看吧| 久久久久久久亚洲中文字幕| 男的添女的下面高潮视频| 国产亚洲最大av| 国产精品伦人一区二区| 看十八女毛片水多多多| 色吧在线观看| 午夜爱爱视频在线播放| 日韩在线高清观看一区二区三区| 国产一区二区三区av在线| 午夜老司机福利剧场| 成年女人在线观看亚洲视频 | 91久久精品电影网| 丰满乱子伦码专区| 免费看av在线观看网站| 亚洲最大成人中文| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院入口| av线在线观看网站| 日本黄大片高清| 青青草视频在线视频观看| 国产久久久一区二区三区| 欧美激情国产日韩精品一区| 日日啪夜夜撸| 2018国产大陆天天弄谢| 国产淫片久久久久久久久| 麻豆成人av视频| 久久这里只有精品中国| 国语对白做爰xxxⅹ性视频网站| 日韩电影二区| 中文欧美无线码| 毛片女人毛片| 国产精品人妻久久久久久| 男人舔女人下体高潮全视频| 国模一区二区三区四区视频| 99久久人妻综合| av在线蜜桃| 欧美性感艳星| 欧美极品一区二区三区四区| 白带黄色成豆腐渣| xxx大片免费视频| 99热这里只有是精品在线观看| 国产成人freesex在线| 国产精品国产三级国产av玫瑰| 欧美丝袜亚洲另类| 熟女电影av网| 久久精品久久精品一区二区三区| 久久久久久久亚洲中文字幕| 蜜臀久久99精品久久宅男| 91精品伊人久久大香线蕉| a级一级毛片免费在线观看| 日韩中字成人| 男女视频在线观看网站免费| 亚洲三级黄色毛片| 少妇被粗大猛烈的视频| 国产精品美女特级片免费视频播放器| 国内精品一区二区在线观看| 精品久久国产蜜桃| 亚洲天堂国产精品一区在线| 亚洲精品日韩av片在线观看| 麻豆精品久久久久久蜜桃| 国产午夜福利久久久久久| 日韩成人av中文字幕在线观看| 80岁老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 又爽又黄无遮挡网站| ponron亚洲| 一本久久精品| 国产精品.久久久| 一个人免费在线观看电影| 国产久久久一区二区三区| 久久热精品热| 六月丁香七月| 免费av观看视频| 色综合亚洲欧美另类图片| 午夜激情欧美在线| 久久99热6这里只有精品| 中国美白少妇内射xxxbb| 亚洲精品色激情综合| 中文字幕制服av| 日韩一本色道免费dvd| 国产精品久久久久久av不卡| 午夜视频国产福利| 亚洲av成人精品一二三区| 日韩 亚洲 欧美在线| 精品酒店卫生间| 久久久久久国产a免费观看| 亚洲不卡免费看| 亚洲精品日韩在线中文字幕| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 丰满乱子伦码专区| 午夜精品一区二区三区免费看| 久久久久久久久大av| 欧美高清性xxxxhd video| 精品人妻视频免费看| 秋霞在线观看毛片| 国产精品综合久久久久久久免费| 男插女下体视频免费在线播放| 免费无遮挡裸体视频| 搡老妇女老女人老熟妇| 亚洲无线观看免费| 综合色av麻豆| 日韩欧美一区视频在线观看 | 亚洲欧美精品专区久久| 国产片特级美女逼逼视频| 亚洲成色77777| 午夜亚洲福利在线播放| 午夜免费观看性视频| 亚洲av.av天堂| 午夜福利成人在线免费观看| 男人舔奶头视频| 精品欧美国产一区二区三| 国产男女超爽视频在线观看| 天美传媒精品一区二区| 国产又色又爽无遮挡免| 亚洲精品色激情综合| 99久久人妻综合| 黄片无遮挡物在线观看| 在线免费十八禁| 熟女电影av网| 国产成人a∨麻豆精品| kizo精华| 亚洲av成人av| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 国产麻豆成人av免费视频| 身体一侧抽搐| 18禁在线无遮挡免费观看视频| 国产精品不卡视频一区二区| 九色成人免费人妻av| 国产色爽女视频免费观看| 久久这里只有精品中国| 国产精品久久久久久久久免| 免费高清在线观看视频在线观看| 麻豆成人午夜福利视频| 免费观看av网站的网址| 嫩草影院入口| 色综合色国产| 啦啦啦中文免费视频观看日本| 国产亚洲精品av在线| 麻豆国产97在线/欧美| 国产色婷婷99| 美女国产视频在线观看| 久久久久久久大尺度免费视频| 国产亚洲5aaaaa淫片| 亚洲成人久久爱视频| 亚洲欧美精品自产自拍| 一级毛片我不卡| 白带黄色成豆腐渣| 69人妻影院| 久久久久久久久久成人| 青春草亚洲视频在线观看| 精品久久久久久久人妻蜜臀av| 国产男人的电影天堂91| 一级爰片在线观看| 国产精品久久久久久精品电影| 熟妇人妻不卡中文字幕| 欧美日韩精品成人综合77777| 久久精品久久久久久久性| 婷婷色麻豆天堂久久| 国产成人aa在线观看| 国产精品综合久久久久久久免费| 好男人在线观看高清免费视频| 免费观看在线日韩| 国产在视频线精品| 在线a可以看的网站| 久久午夜福利片| 黄色日韩在线| 18禁在线无遮挡免费观看视频| 成年版毛片免费区| 国产黄片视频在线免费观看| 日韩国内少妇激情av| 视频中文字幕在线观看| 亚洲性久久影院| 偷拍熟女少妇极品色| 色5月婷婷丁香| 国产高清国产精品国产三级 | 九色成人免费人妻av| 中文字幕久久专区| 国产亚洲91精品色在线| 日韩视频在线欧美| 国产男人的电影天堂91| 精品久久久噜噜| 国产亚洲91精品色在线| 精品一区二区三卡| 国产成人精品婷婷| 成人一区二区视频在线观看| 亚洲国产欧美人成| 欧美高清成人免费视频www| 亚洲最大成人av| 亚洲无线观看免费| 搡老妇女老女人老熟妇| 欧美不卡视频在线免费观看| 性插视频无遮挡在线免费观看| 欧美激情久久久久久爽电影| 卡戴珊不雅视频在线播放| 国产男女超爽视频在线观看| 成人无遮挡网站| 你懂的网址亚洲精品在线观看| 男插女下体视频免费在线播放| 精品国产一区二区三区久久久樱花 | 激情 狠狠 欧美| 色5月婷婷丁香| 一夜夜www| 看十八女毛片水多多多| 欧美成人一区二区免费高清观看| 夫妻性生交免费视频一级片| 99久久精品国产国产毛片| 成人毛片a级毛片在线播放| 欧美xxxx性猛交bbbb| 国产av在哪里看| 麻豆精品久久久久久蜜桃| 日韩欧美精品v在线| 18禁动态无遮挡网站| 日本wwww免费看| 国产高清三级在线|