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

    高梯度磁場(chǎng)提升單纖維捕集PM2.5性能的機(jī)理

    2019-07-31 07:41:52張儷安刁永發(fā)莊加瑋周發(fā)山沈恒根
    中國(guó)環(huán)境科學(xué) 2019年7期
    關(guān)鍵詞:磁化率磁場(chǎng)強(qiáng)度磁性

    張儷安,刁永發(fā),莊加瑋,周發(fā)山,沈恒根

    高梯度磁場(chǎng)提升單纖維捕集PM2.5性能的機(jī)理

    張儷安,刁永發(fā)*,莊加瑋,周發(fā)山,沈恒根

    (東華大學(xué)環(huán)境科學(xué)與工程學(xué)院,上海 201620)

    以鋼鐵廠和有色金屬行業(yè)排放的PM2.5為研究對(duì)象,基于離散相模型DPM(Discrete Phase Model),并加入U(xiǎn)DF自定義編程,研究高梯度磁場(chǎng)下不同入口風(fēng)速、顆粒粒徑、外磁場(chǎng)強(qiáng)度、磁性纖維磁感應(yīng)強(qiáng)度以及磁化率對(duì)捕集效率的影響,并結(jié)合顆粒運(yùn)動(dòng)軌跡和受力情況對(duì)其進(jìn)行分析.結(jié)果表明:當(dāng)0.5μm£p£2.5μm,=0.1m/s時(shí),利用高梯度磁場(chǎng)(=0.1T,=0.06T)可以使單纖維捕集PM2.5的效率提高為原來(lái)的4.23倍,得出磁性纖維周圍存在2個(gè)引力區(qū)和2個(gè)斥力區(qū).同時(shí),在高梯度磁場(chǎng)中磁性纖維對(duì)PM2.5的捕集效率隨入口風(fēng)速呈先減小后趨于平穩(wěn)的規(guī)律;而捕集效率隨粉塵粒徑呈先增大后減小的規(guī)律.當(dāng)p=1.0μm時(shí)的捕集效率提升最大,無(wú)論是外磁場(chǎng)強(qiáng)度還是磁性纖維磁感應(yīng)強(qiáng)度,磁性纖維對(duì)顆粒的捕集效率與場(chǎng)強(qiáng)都呈一次函數(shù)關(guān)系,效率增長(zhǎng)率B>H;隨著顆粒磁化率的增加,磁性纖維對(duì)顆粒的捕集呈現(xiàn)兩段線性增長(zhǎng)規(guī)律,前后兩段效率增長(zhǎng)率1>2.當(dāng)顆粒經(jīng)過(guò)高梯度磁場(chǎng)區(qū)域時(shí),入口風(fēng)速、粉塵粒徑、場(chǎng)強(qiáng)對(duì)運(yùn)動(dòng)軌跡影響較大,而磁化率對(duì)運(yùn)動(dòng)軌跡影響較小.

    高梯度磁場(chǎng);PM2.5;捕集效率;磁化率

    隨著鋼鐵和有色金屬行業(yè)的發(fā)展,其生產(chǎn)過(guò)程中會(huì)產(chǎn)生大量PM2.5,由于數(shù)量多、沉降速度慢、比表面積大、可作為其它污染物的載體,對(duì)人體健康和大氣環(huán)境質(zhì)量的影響很大[1-2].

    為了更好的研究纖維對(duì)顆粒的捕集規(guī)律,國(guó)內(nèi)外研究者從機(jī)理上系統(tǒng)的研究了單纖維對(duì)顆粒物捕集.前期階段研究多集中在經(jīng)驗(yàn)公式和半經(jīng)驗(yàn)公式[3-5]的研究.隨后,為了更好的研究單纖維對(duì)顆粒的捕集規(guī)律,其中一些學(xué)者從不同工況的角度較為系統(tǒng)的考察了影響單纖維捕集顆粒的因素[6-8],這些研究深刻揭示了單纖維在捕集顆粒過(guò)程中隨工況改變的變化規(guī)律.而另一部分學(xué)者則專注于研究顆粒在纖維周圍積聚、沉積位點(diǎn)以及三維樹(shù)枝狀顆粒沉積物的形成過(guò)程[9-11].同時(shí),在研究方法上,DEM (Discrete Element Model)方法的使用[12-13]不僅可以分析單纖維在流場(chǎng)中的放置方式、表面黏附系數(shù)、顆粒直徑以及雷諾數(shù)對(duì)捕集效率等傳統(tǒng)工況的影響,而且顆粒的可視化可以更加清晰的觀察粉塵顆粒被纖維捕集過(guò)程.

    通過(guò)總結(jié)發(fā)現(xiàn),當(dāng)前對(duì)于單纖維捕集顆粒的研究主要體現(xiàn)在經(jīng)驗(yàn)公式、工況參數(shù)以及纖維沉積過(guò)程的研究,相關(guān)研究無(wú)法解決PM2.5等微細(xì)顆粒由于存在穿透窗口而難以捕集的問(wèn)題[14].因此,需要在原有的基礎(chǔ)上開(kāi)發(fā)新的技術(shù).

    現(xiàn)階段高梯度磁場(chǎng)技術(shù)在處理細(xì)顆?;蛭⒓?xì)顆粒弱磁性物料方面因具有分離能力高,結(jié)構(gòu)簡(jiǎn)單和維護(hù)成本低等優(yōu)勢(shì)而被廣泛使用.Svoboda等[15]對(duì)高梯度磁場(chǎng)分離進(jìn)行研究,通過(guò)增加外磁場(chǎng)、磁介質(zhì)磁化強(qiáng)度可有效提高顆粒物的分離效率.Ravnik等[16]通過(guò)計(jì)算顆粒在流體和磁場(chǎng)力作用下的運(yùn)動(dòng)軌跡,驗(yàn)證了在窄通道中進(jìn)行高梯度磁選的可行性.Baik等[17]對(duì)高梯度磁場(chǎng)研究結(jié)果表明,在流體中作用于顆粒的磁場(chǎng)力與磁通密度和磁場(chǎng)梯度成正比,且捕集顆粒能力方面,高梯度磁場(chǎng)系統(tǒng)要強(qiáng)于永磁性系統(tǒng).Zheng等[18]研究表明高梯度磁場(chǎng)下橢圓截面矩陣和圓形截面矩陣對(duì)于微米級(jí)顆粒捕集效果明顯.通過(guò)總結(jié)發(fā)現(xiàn),雖然利用高梯度磁場(chǎng)對(duì)微細(xì)顆粒的捕集已有研究,但是利用該技術(shù)在纖維捕集粉塵領(lǐng)域的研究卻鮮見(jiàn)報(bào)道,尤其是在鋼鐵和有色金屬行業(yè),由于所排放的粉塵含有鐵磁性物質(zhì)[19]而更容易被捕集.因此為了解決目前單纖維模型存在預(yù)測(cè)與實(shí)驗(yàn)值差別較大的問(wèn)題,通過(guò)UDF編程的方法來(lái)提出改進(jìn)模型,以彌補(bǔ)研究的不足.

    其次在研究高梯度磁場(chǎng)下纖維捕集粉塵顆粒的過(guò)程中,磁性纖維可直接通過(guò)紡絲或基體纖維的物理、化學(xué)改性制備.磁性纖維在工業(yè)使用中可分離含塵煙氣中的鐵磁性物質(zhì),通過(guò)后期的間歇性電磁振打方式實(shí)現(xiàn)回收[20];且磁性濾料配合磁性吸附劑可實(shí)現(xiàn)對(duì)燃煤煙氣中痕量Hg0[21]的脫除,在除塵的同時(shí)凈化有毒污染氣體;此外,外磁場(chǎng)的加入可使磁性纖維形成非均勻的致密磁場(chǎng),增加磁場(chǎng)強(qiáng)度.因此,本研究通過(guò)對(duì)高梯度磁場(chǎng)下磁性纖維對(duì)PM2.5的捕集過(guò)程為研究對(duì)象,建立捕集PM2.5的物理模型.基于CFD-DPM方法,加入U(xiǎn)DF自定義編程對(duì)PM2.5在高梯度磁場(chǎng)中的被捕集規(guī)律進(jìn)行數(shù)值模擬,計(jì)算在不同工況下PM2.5的運(yùn)動(dòng)軌跡,考察了入口風(fēng)速、粉塵粒徑、外磁場(chǎng)強(qiáng)度和磁性纖維磁場(chǎng)強(qiáng)度以及顆粒磁化率對(duì)捕集效率的影響,旨在為高梯度磁場(chǎng)下磁性纖維對(duì)PM2.5的捕集優(yōu)化設(shè)計(jì)提供理論指導(dǎo).

    1 計(jì)算模型

    1.1 氣固兩相流模型

    對(duì)于內(nèi)部氣?固兩相流動(dòng)進(jìn)行數(shù)值模擬需先計(jì)算氣相場(chǎng).可采用標(biāo)準(zhǔn)k-ε模型、穩(wěn)態(tài)及不可壓縮模型進(jìn)行數(shù)值模擬.控制方程(連續(xù)性方程、動(dòng)量方程)[22]如下:

    文中用二階迎風(fēng)格式SIMPLE算法對(duì)離散化動(dòng)量方程進(jìn)行壓力速度耦合求解,并將連續(xù)性方程及動(dòng)量方程在直角坐標(biāo)系,,方向上的收斂殘差設(shè)定在10?6以內(nèi).一旦獲得穩(wěn)定的氣相流場(chǎng),就將顆粒從進(jìn)口以面射流源形式注入計(jì)算區(qū)域.假定顆粒在流場(chǎng)中做無(wú)旋運(yùn)動(dòng),顆粒的運(yùn)動(dòng)平衡方程表達(dá)式[23]如下:

    式中:vp,v分別為顆粒運(yùn)動(dòng)速度和流體速度, m/s;FD為流體的曳力,N;g為重力加速度,m/s2;ρp,ρa(bǔ)分別為顆粒和空氣的密度,kg/m3;FM為通過(guò)UDF編程在高梯度磁場(chǎng)中所受的磁場(chǎng)力,N;Fother為受到的一些可忽略的力,N.

    圖1中,通過(guò)建立極坐標(biāo)系將M沿徑向和切向進(jìn)行分解,如下[24]:

    式中:0為真空磁導(dǎo)率,2.256′10-6;p為磁化率;和分別為外磁場(chǎng)強(qiáng)度和磁性纖維的磁感應(yīng)強(qiáng)度,T;和分別為極坐標(biāo)下的極徑和極角.

    雖然作用在顆粒上的力相當(dāng)復(fù)雜,會(huì)受到壓力梯度力、Basset力、Magnus力、Saffman升力等一系列的作用力[25-26].但是在單纖維捕集顆粒模型中,由于各力在適用條件下對(duì)顆粒的影響很小,為了簡(jiǎn)化計(jì)算,均可忽略不計(jì)[27].

    1.2 單纖維捕集效率的計(jì)算方法

    式中:為單纖維捕集的效率;in為入口通入的顆粒數(shù);out為出口逃逸的顆粒數(shù).

    1.3 邊界條件設(shè)置

    圖2 計(jì)算區(qū)域及邊界條件設(shè)置示意

    具體的邊界條件如圖2所示,以單纖維捕集顆粒模型為例,具體邊界條件設(shè)置如下:計(jì)算域入口邊界設(shè)為速度進(jìn)口,出口邊界設(shè)為壓力出口;在模型計(jì)算中,將纖維橫截面面積與假想控制面的比值作為填充率對(duì)待[28],即=f2/2,同時(shí)根據(jù)Davies結(jié)合實(shí)驗(yàn)結(jié)果得出的結(jié)論,在0.6%~30%的范圍內(nèi)都是正確的[3].當(dāng)入口高度=5.3f時(shí),則填充率=0.035,因此模型的長(zhǎng)、寬、高分別定為240,120,80μm是可行的;纖維的直徑為15μm;計(jì)算區(qū)域四周的邊界則根據(jù)單纖維結(jié)構(gòu)特點(diǎn)設(shè)為對(duì)稱邊界條件;纖維表面邊界設(shè)為無(wú)滑移邊界條件;圖3為高梯度磁場(chǎng)的磁場(chǎng)形式,高梯度磁場(chǎng)的形成則是在均勻的背景磁場(chǎng)中填充飽和聚磁介質(zhì)以產(chǎn)生高磁場(chǎng)梯度磁場(chǎng);圖4為鋼鐵廠生產(chǎn)過(guò)程中產(chǎn)生粉塵的XRD圖譜,主要成分為CaCO3,當(dāng)2=35.42°和2=44.14°,分別出現(xiàn)了Fe3O4和Fe的特征峰,說(shuō)明鋼鐵廠排放的粉塵具有一定的鐵磁性,容易被磁化.

    圖3 高梯度磁場(chǎng)內(nèi)圓柱形磁纖維周圍磁場(chǎng)示意

    圖4 鋼鐵廠產(chǎn)生粉塵的XRD圖譜

    1.4 網(wǎng)格獨(dú)立性檢驗(yàn)

    為了去除網(wǎng)格數(shù)量對(duì)數(shù)值模擬計(jì)算準(zhǔn)確性的影響,對(duì)模型進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證,計(jì)算不同網(wǎng)格密度下的壓力損失,模擬結(jié)果如圖5所示.隨著網(wǎng)格密度的增加,單纖維模型結(jié)構(gòu)的壓力損失逐漸增大,且網(wǎng)格數(shù)為14萬(wàn)、55萬(wàn)、90萬(wàn)左右時(shí),單纖維模型結(jié)構(gòu)的壓力損失隨入口風(fēng)速變規(guī)律一致,與Davies壓力損失經(jīng)驗(yàn)公式(6)的對(duì)比誤差都在5%范圍內(nèi). 選取其中一種工況進(jìn)行效率計(jì)算,當(dāng)速度=0.2m/s,p=2.5μm時(shí),單纖維捕集效率與Davies效率經(jīng)驗(yàn)公式(7)的誤差分別為8.51%、3.50%、1.50%,后兩者的誤差都在5%范圍內(nèi),根據(jù)網(wǎng)格數(shù)量和誤差綜合考慮選取55萬(wàn)的網(wǎng)格用于數(shù)值模型的計(jì)算,且本文的計(jì)算模型網(wǎng)格采用的是六面體結(jié)構(gòu)化網(wǎng)格.

    1.5 與經(jīng)驗(yàn)公式對(duì)比

    為了驗(yàn)證單纖維結(jié)構(gòu)捕集粉塵顆粒數(shù)值模擬的準(zhǔn)確性,計(jì)算了單纖維結(jié)構(gòu)模型的捕集效率,并與實(shí)驗(yàn)計(jì)算公式進(jìn)行對(duì)比,由圖6,7可知,過(guò)濾效率的誤差在10%范圍內(nèi).粒徑較小時(shí)模擬值與經(jīng)驗(yàn)公式相差較大,這是因?yàn)镈avies效率經(jīng)驗(yàn)公式(7)只考慮了顆粒的攔截作用和慣性碰撞,而當(dāng)粉塵顆粒0.5μm

    式中:P為直接碰撞系數(shù);St為斯托克斯數(shù);0為計(jì)算修正后的捕集效率;為計(jì)算捕集效率;為填充率.

    圖6 不同粒徑下單纖維捕集效率計(jì)算

    圖7 不同St下單纖維捕集效率計(jì)算

    計(jì)算結(jié)果與經(jīng)驗(yàn)公式間的誤差基本可接受,表明采用該方法進(jìn)行單纖維捕集粉塵顆粒的研究可行.

    2 結(jié)果與討論

    2.1 高梯度磁場(chǎng)下PM2.5顆粒的運(yùn)動(dòng)規(guī)律

    當(dāng)PM2.5被氣體攜帶進(jìn)入高梯度磁場(chǎng)時(shí),由于纖維的長(zhǎng)度遠(yuǎn)大于纖維的直徑,因此不考慮各流場(chǎng)和磁場(chǎng)中各物理量沿磁性纖維軸向的變化.PM2.5剛進(jìn)入時(shí),分布均勻.隨著顆粒的繼續(xù)移動(dòng),當(dāng)粉塵顆粒進(jìn)入高梯度磁場(chǎng)區(qū)域時(shí),由于主要受到磁場(chǎng)力、曳力、重力以及布朗力的共同作用,運(yùn)動(dòng)狀態(tài)發(fā)生改變.通過(guò)顆粒的運(yùn)動(dòng)軌跡8(a)和8(b)對(duì)比可知,在高梯度磁場(chǎng)中,磁性纖維周圍存在2個(gè)引力區(qū)和2個(gè)斥力區(qū).這與孫仲元等[28]的磁選理論研究一致.通過(guò)與傳統(tǒng)流場(chǎng)形式下單纖維對(duì)顆粒的捕集進(jìn)行對(duì)比,可知當(dāng)加入高梯度磁場(chǎng)時(shí),單纖維捕集顆粒能力明顯提高.

    圖8 高梯度磁場(chǎng)下磁性纖維捕集顆粒的運(yùn)動(dòng)軌跡

    2.2 高梯度磁場(chǎng)下磁場(chǎng)對(duì)粉塵捕集效率的影響

    2.2.1 入口風(fēng)速對(duì)捕集效率的影響 如圖9可知,磁性纖維對(duì)PM2.5顆粒的捕集效率隨入口風(fēng)速呈先減小后趨于平穩(wěn)的規(guī)律.隨著入口風(fēng)速的增加,磁性纖維對(duì)粉塵顆粒捕集效率逐漸減小.

    圖9 高梯度磁場(chǎng)下捕集效率與入口風(fēng)速的關(guān)系

    圖10中,其余軌跡圖中速度與磁場(chǎng)方向都與本圖相同.當(dāng)p=1.0μm時(shí),隨著入口風(fēng)速的增加,在引力區(qū)落在磁性纖維上的顆粒數(shù)逐漸減小,且斥力區(qū)的“空腔”縮小.這是因?yàn)?當(dāng)顆粒接近磁性纖維時(shí),雖然在斥力區(qū)所受的斥力相同,但是顆粒速度越大,運(yùn)動(dòng)狀態(tài)越不易改變,在斥力區(qū)運(yùn)動(dòng)的距離越長(zhǎng),導(dǎo)致斥力區(qū)的“空腔”縮小.當(dāng)顆粒進(jìn)入引力區(qū)時(shí),由于氣流越大對(duì)顆粒的攜帶能力越強(qiáng),在引力區(qū)被流體帶走的顆粒越多,被引力區(qū)吸引捕集的顆粒就越少.同時(shí),氣流速度越大,相應(yīng)的顆粒在磁場(chǎng)中的作用時(shí)間越短,顆粒所受的磁場(chǎng)作用效果就會(huì)在一定程度上減弱,因此捕集效率減少.

    圖10 不同入口風(fēng)速下磁性纖維捕集顆粒軌跡圖(H=0.1T,B=0.06T,cp=0.025)

    2.2.2 粉塵粒徑對(duì)捕集效率的影響 由圖11可知,當(dāng)0.5μm£p£2.5μm,=0.1m/s時(shí),隨著粉塵粒徑的增加,磁性纖維對(duì)PM2.5的捕集效率呈現(xiàn)先增大后減小的趨勢(shì),利用高梯度磁場(chǎng)(=0.1T,=0.06T)可以使單纖維捕集PM2.5的效率提高為原來(lái)的4.23倍.當(dāng)p=1.0μm,此時(shí)的捕集效率提高最大.通過(guò)粉塵顆粒的運(yùn)動(dòng)軌跡圖12可知,當(dāng)粉塵粒徑較小時(shí),氣流攜帶顆粒能力強(qiáng),被氣流攜帶的顆粒在經(jīng)過(guò)引力區(qū)時(shí)很少被磁性纖維捕獲,因此捕集效率較低.隨著顆粒的粒徑增大,在引力區(qū)磁場(chǎng)對(duì)顆粒的磁力增大,磁性纖維對(duì)顆粒捕集效率提高,隨著粉塵顆粒的繼續(xù)增大,顆粒在斥力區(qū)受到斥力作用同步增強(qiáng),在接近纖維的過(guò)程中由于受斥力作用而遠(yuǎn)離磁性纖維.顆粒粒徑越大,遠(yuǎn)離纖維距離越大,即“空腔”增大.當(dāng)顆粒再經(jīng)過(guò)引力區(qū)時(shí),由于磁場(chǎng)強(qiáng)度會(huì)隨著與纖維間距的增大而減弱,顆粒與纖維間距的增加使得顆粒所受的磁場(chǎng)力減小,再加上流場(chǎng)的作用,在引力區(qū)顆粒受引力運(yùn)動(dòng)方向幾乎與流場(chǎng)垂直,受流場(chǎng)影響很大,此時(shí)很難再被纖維捕集,捕集效率降低.

    2.2.3 磁場(chǎng)強(qiáng)度對(duì)捕集效率的影響 由圖13(a)可知,磁性纖維對(duì)粉塵顆粒的捕集效率與磁性纖維磁感應(yīng)強(qiáng)度呈一次函數(shù)關(guān)系.隨著磁性纖維磁感應(yīng)強(qiáng)度增加,磁性纖維對(duì)PM2.5的捕集作用逐漸增強(qiáng),與此同時(shí),當(dāng)施加的外磁場(chǎng)強(qiáng)度不同時(shí),如圖13(b)所示,磁性纖維對(duì)粉塵顆粒的捕集效率與外磁場(chǎng)強(qiáng)度同樣呈一次函數(shù)關(guān)系.隨著外磁場(chǎng)強(qiáng)度的增加,磁性纖維對(duì)粉塵的捕集作用逐漸增強(qiáng),通過(guò)擬合得出圖13(a)和13(b)中兩個(gè)一次函數(shù)的平均增長(zhǎng)斜率分別為137.54和12.79,即捕集效率的增長(zhǎng)率B>H,磁性纖維磁感應(yīng)強(qiáng)度對(duì)PM2.5顆粒捕集效率的影響程度要大于外磁場(chǎng)強(qiáng)度的影響程度.

    圖11 高梯度磁場(chǎng)下捕集效率與粉塵粒徑的關(guān)系

    圖12 不同粒徑下磁性纖維捕集顆粒軌跡(H=0.1T,B=0.06T,=0.025)

    圖13 高梯度磁場(chǎng)下捕集效率與磁場(chǎng)強(qiáng)度的關(guān)系

    (a)磁性纖維磁感應(yīng)強(qiáng)度; (b)外加磁場(chǎng)強(qiáng)度

    通過(guò)粉塵顆粒的運(yùn)動(dòng)軌跡(圖 14)可知,當(dāng)=0.2m/s,p=1.0μm時(shí),隨著外磁場(chǎng)強(qiáng)度的增加,在引力區(qū)中落入磁性纖維表面的顆粒增多.主要原因是當(dāng)外磁場(chǎng)增加時(shí),磁場(chǎng)范圍和強(qiáng)度增大,且磁性纖維附近的磁場(chǎng)梯度增加,粉塵顆粒所受的磁場(chǎng)力相應(yīng)增強(qiáng),因此捕集效率增大.

    圖14 不同外磁場(chǎng)強(qiáng)度下磁性纖維捕集顆粒軌跡圖(B=0.06T,=0.025)

    圖15 高梯度磁場(chǎng)下捕集效率與顆粒磁化率的關(guān)系

    2.2.4 顆粒的磁化率對(duì)捕集效率的影響 當(dāng)PM2.5進(jìn)入高梯度磁場(chǎng)時(shí),粉塵顆粒被磁化,繼而受到磁場(chǎng)力的作用.由圖15可知,當(dāng)=0.2m/s,p=1.0μm時(shí),隨著粉塵顆粒磁化率的增加,磁性纖維對(duì)粉塵顆粒的捕集規(guī)律呈現(xiàn)出兩段線性增加規(guī)律,且起始段線性增加斜率要大于后段線性增加的斜率.這是因?yàn)楫?dāng)PM2.5磁化率較小時(shí),此時(shí)顆粒磁性雖然很弱,但是相對(duì)于無(wú)磁性的PM2.5,磁場(chǎng)力的出現(xiàn)會(huì)使得顆粒在磁場(chǎng)中所受的磁場(chǎng)力顯著變化,因此增長(zhǎng)的斜率較大.隨著磁化率的繼續(xù)增加,此時(shí)顆粒的磁性同樣增強(qiáng),但磁場(chǎng)力的變化相比之前減弱,因此增長(zhǎng)的斜率減小.通過(guò)PM2.5的運(yùn)動(dòng)軌跡(圖16)可知,當(dāng)=0.2m/s,p=1.0μm時(shí),隨著顆粒磁化率的增加,在引力區(qū)中被磁性纖維表面捕集顆粒數(shù)逐漸增加.根據(jù)“空腔”的大小可知磁化率的大小對(duì)運(yùn)動(dòng)軌跡影響較小.

    圖16 不同磁化率下磁性纖維捕集顆粒軌跡圖(H=0.1T,B=0.06T)

    3 結(jié)論

    3.1 PM2.5在高梯度磁場(chǎng)中運(yùn)動(dòng)時(shí),磁場(chǎng)力對(duì)于PM2.5的運(yùn)動(dòng)軌跡影響較大,當(dāng)0.5μm£p£2.5μm,=0.1m/s時(shí),利用高梯度磁場(chǎng)(=0.1T,=0.06T)可以使單纖維捕集PM2.5的效率提高為原來(lái)的4.23倍,且在磁性纖維周圍存在2個(gè)引力區(qū)和2個(gè)斥力區(qū).

    3.2 當(dāng)粉塵粒徑、磁場(chǎng)強(qiáng)度、磁化率一定時(shí),在高梯度磁場(chǎng)中,磁性纖維對(duì)PM2.5的捕集效率隨入口風(fēng)速呈先減小后趨于平穩(wěn)的規(guī)律,當(dāng)入口風(fēng)速、磁場(chǎng)強(qiáng)度、磁化率一定時(shí),磁性纖維對(duì)PM2.5的捕集效率隨粉塵粒徑的增加呈先增大后減小的趨勢(shì),當(dāng)p= 1.0μm時(shí)的捕集效率提高最大.

    3.3 無(wú)論是外磁場(chǎng)強(qiáng)度還是磁性纖維磁感應(yīng)強(qiáng)度,當(dāng)入口風(fēng)速、顆粒粒徑一定時(shí),磁性纖維對(duì)粉塵顆粒的捕集效率與磁場(chǎng)強(qiáng)度都呈一次函數(shù)的關(guān)系,隨著磁場(chǎng)強(qiáng)度的增加,捕集效率增加,且增長(zhǎng)斜率B>H.

    3.4 當(dāng)入口風(fēng)速、粉塵粒徑、磁場(chǎng)強(qiáng)度一定時(shí),隨著PM2.5磁化系數(shù)的增加,磁性纖維對(duì)PM2.5的捕集規(guī)律呈現(xiàn)兩段線性增加規(guī)律,且開(kāi)始段線性增加的斜率要大于后段,磁化率對(duì)顆粒的運(yùn)動(dòng)軌跡影響較小.

    [1] Wang W, Tang D, Liu H. Research on current pollution status and pollution characteristics of PM2.5in china [J]. Research of environmental sciences, 2000,13(1):1-4.

    [2] 邢瓊予,戴啟立,畢曉輝,等.我國(guó)中西部典型城市PM2.5中痕量金屬的時(shí)空分布特征和健康影響[J]. 中國(guó)環(huán)境科學(xué), 2019,39(2):574- 582.Xing Q Y, Dai Q L, Bi X H, et al. Temporal-spatial variation and health effects of trace metals in PM2.5in four central-cities of China [J]. China Environmental Science, 2019,39(2):574-582.

    [3] Davies C N. Air filtration [M]. London: Academic press, 1973:237.

    [4] Endo Y, Chen D R, Pui D Y H. Effects of particle polydispersity and shape factor during dust cake loading on air filters [J]. Powder Technology, 1998,98(3):241-249.

    [5] Badel E, Delisee C, Lux J. 3D structural characterisation, deformation measurements and assessment of low-density wood fibreboard under compression: The use of X-ray microtomography [J]. Composites Science and Technology, 2008,68(7/8):1654-1663.

    [6] 顧叢匯,呂士武,李 瑞,等.纖維對(duì)PM2.5過(guò)濾性能的影響 [J]. 化工學(xué)報(bào), 2014,65(6):2137-2147. Gu C H, Lv S W, Li R, et al. Influence of fiber on filtration performance for PM2.5[J]. CIES Journal, 2014,65(6):2137-2148.

    [7] Shu Z Y, Qian F P, Huang X P. Numerical study on the effect of wall parameters on single fiber trapping efficiency [J]. Journal of Filtration and Separation, 2017,27(1):1-5.

    [8] 楊 林,付海明,李 杰,等.單纖維對(duì)慣性顆粒穩(wěn)態(tài)過(guò)濾捕集效率的數(shù)值模擬分析 [J]. 東華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014,40(3):345- 349. Yang L, Fu H M, Li J, et al. Numerical simulation analysis of stationary collection efficiency of inertial particles on single fiber surface [J]. Journal of DongHua University (Nature Science), 2014, 40(3):345-349.

    [9] Kasper G, Schollmeier S, Meyer J. Structure and density of deposits formed on filter fibers by inertial particle deposition and bounce [J]. Journal of Aerosol Science, 2010,41(12):1167-1182.

    [10] Hosseini S A, Tafreshi H V. Modeling particle-loaded single fiber efficiency and fiber drag using ANSYS–Fluent CFD code [J]. Computers and Fluids, 2012,66:157-166.

    [11] Bang J, Yoon W. Stochastic analysis of a collection process of submicron particles on a single fiber accounting for the changes in flow field due to particle collection [J]. Journal of Mechanical Science and Technology, 2014,28(9):3719-3732.

    [12] Li S Q, Marshall J S. Discrete element simulation of micro-particle deposition on a cylindrical fiber in an array [J]. Journal of Aerosol Science, 2007,38(10):1031-1046.

    [13] Yue C, Zhang Q, Zhai Z Q. Numerical simulation of the filtration process in fibrous filters using CFD-DEM method [J]. Journal of Aerosol Science, 2016,101:174-187.

    [14] 熊桂龍,李水清,陳 晟,等.增強(qiáng)PM2.5脫除的新型電除塵技術(shù)的發(fā)展 [J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2015,35(9):2217-2213. Xiong G L, Li S Q, Chen M, et al. Development of advanced electrostatic precipitation technologies for reducing PM2.5emissions from coal-fired power plants [J]. Proceedings ofthe CSEE, 2015, 35(9):2217-2223.

    [15] Svoboda J, Fujita T. Recent developments in magnetic methods of material separation [J]. Minerals Engineering, 2003,16(9):785-792.

    [16] Ravnik J, Hriber M. High gradient magnetic particle separation in viscous flows by 3D BEM [J]. Computational Mechanics, 2013,51(4): 465-474.

    [17] Baik S K , Ha D W , Kwon J M , et al. Magnetic force on a magnetic particle within a high gradient magnetic separator [J]. Physica C Superconductivity, 2013,484(484):333-337.

    [18] Zheng X, Wang Y, Lu D. Investigation of the particle capture of elliptic cross-sectional matrix for high gradient magnetic separation [J]. Powder Technology, 2016,297:303-310.

    [19] 吳勝利,張風(fēng)杰,張建良,等.鋼廠含鋅粉塵基本物性及其成球性能研究[J]. 環(huán)境工程, 2015,33(7):90-95. Wu S L, Zhang F J, Zhang J L, et al. Basic characteristics and balling performance of Zinc-Bearing dust in steel plant [J]. Environmental Engineering, 2015,33(7):90-95.

    [20] 趙長(zhǎng)遂.燃煤可吸入顆粒物的脫除方法:中國(guó), CN201172006 [P]. 2007-11-14. Zhao C S. Method for the removal of inhalable particulate matter from coal burning: Chana, CN201172006 [P]. 2007-11-14.

    [21] 趙 艷,周發(fā)山,刁永發(fā),等.磁性載銀濾料對(duì)燃煤煙氣中Hg0的吸附性能[J]. 環(huán)境工程學(xué)報(bào), 2017,11(12):6339-6344. Zhao Y, Zhou F S, Diao Y F, et al. Adsorption performance of magnetic silver-loaded filter material on Hg0from coal combustion flue gas [J]. Chinese Journal of Environmental Engineering, 2017, 11(12):6339-6344.

    [22] Qian F, Huang N, Zhu X, et al. Numerical study of the gas–solid flow characteristic of fibrous media based on SEM using CFD–DEM [J]. Powder Technology, 2013,249:63-70.

    [23] 朱紅軍,林元華,謝龍漢. Fluent12流體分析及工程仿真 [M]. 北京:清華大學(xué)出版社, 2011:186-188. Zhu H J, Lin Y H, Xie L H. Fluid analysis and engineering simulation of Fluent 12 [M]. Beijing: Tsinghua University Press, 2011:186-188.

    [24] 王發(fā)揮,王秀芳,張 丹.高梯度磁場(chǎng)中多根磁介質(zhì)捕集磁性微粒的數(shù)值模擬 [J]. 金屬礦山, 2010,39(2):103-106. Wang F H, Wang X F, Zhang D. Numerical simulation for multiple magnetic medium in entrapping magnetism particles in high gradient magnetic field [J]. Metal Mine, 2010,39(2):103-106.

    [25] 林金賢,林 棋,婁 晨,等.氣固兩相流管道流動(dòng)阻力特性數(shù)值模擬 [J]. 油氣儲(chǔ)運(yùn), 2014,33(1):32-41. Lin J X, Lin Q, Lou C, et al. Numerical simulation of bow drag characteristics of gas-solid two-phase bow in pipeline [J]. Oil and Gas Storage and Transportation, 2014,33(1):32-41.

    [26] Tripathy S K, Bhoja S K, Raghu Kumar C, et al. A short review on hydraulic classification and its development in mineral industry [J]. Powder Technology, 2015,270:205-220.

    [27] 朱 輝,付海明,亢燕銘.單纖維過(guò)濾阻力與慣性捕集效率數(shù)值分析[J]. 中國(guó)環(huán)境科學(xué), 2017,37(4):1298-1306. Zhu H, Fu H M, Kang Y M. Numerical analysis of pressure drop and inertial collection efficiency of a single fiber [J]. China Environmental Science, 2017,37(4):1298-1306.

    [28] 雷陳磊.單纖維非穩(wěn)態(tài)過(guò)濾的捕集效率及其壓降特性研究 [D]. 上海:東華大學(xué), 2015. Lei C L. Research on collection efficiency and pressure drop of a particle-loaded single fiber [D]. Shanghai: DongHua University, 2015.

    [29] 孫中元.磁選理論 [M]. 長(zhǎng)沙:中南大學(xué)出版社, 2007. Sun Z Y. Magnetic separation theory [M]. Changsha: Central South University Press, 2007:220-225.

    The mechanism of high gradient magnetic field improving the performance of single fiber capture PM2.5.

    ZHANG Li-an, DIAO Yong-fa*, ZHUANG Jia-wei, ZHOU Fa-shan, SHEN Heng-gen

    (College of Environmental Science and Engineering, Dong Hua University, Shanghai 201620, China)., 2019,39(7):2765~2773

    PM2.5emissions from iron and steel and nonferrous industries as object, based on the discrete particle model DPM and UDF custom programming, The effects of inlet velocity, particle diameter, external magnetic field strength, magnetic fiber magnetic flux density and magnetic susceptibility on the capture efficiency were studied and analyzed by combination between the particle motion trajectory and force. The results showed that when 0.5μm£p£2.5μm,=0.1m/s, the efficiency of single fiber to capture PM2.5could be improved by using high gradient magnetic field(=0.1T,=0.06T)4.23times. It was also found that there were two gravitational zones and two repulsive zones around the magnetic fibers. At the same time, in the high gradient magnetic field, the capture efficiency of the magnetic fiber to PM2.5decreased first and then stabilized with the inlet velocity; while the capture efficiency increased first and then decreased with the particle diameter of the dust. Whenp=1.0μm, the increasing of capture efficiency was maximized at this time; whether it was the external magnetic field strength or the magnetic fiber magnetic flux density, the capturing efficiency of magnetic fiber to particles was a linear function of the field intensity, with the efficiency growth rate ofB>H; with the increase of the magnetic susceptibility of the particles, the capturing of particles by magnetic fibers presented a two-stage linear growth law, and the growth rate1>2. When the particles passed through the high gradient magnetic field, the inlet velocity, dust particle diameter and field strength had a great influence on the motion trajectory, while the magnetic susceptibility had little effect on the motion trajectory.

    high gradient magnetic field;PM2.5;capturing efficiency;magnetic susceptibility

    X513

    A

    1000-6923(2019)07-2765-09

    張儷安(1990-),男,安徽淮北市人,東華大學(xué)博士研究生,主要從事PM2.5顆粒的磁團(tuán)聚研究.發(fā)表論文1篇.

    2018-12-20

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2018YFC0705300);中央高?;究蒲袠I(yè)務(wù)費(fèi)重點(diǎn)項(xiàng)目(2232017A-09)

    * 責(zé)任作者, 教授, diaoyongfa@dhu.edu.cn

    猜你喜歡
    磁化率磁場(chǎng)強(qiáng)度磁性
    關(guān)于醫(yī)用磁共振成像系統(tǒng)(MRI)磁場(chǎng)強(qiáng)度建標(biāo)
    一種永磁種子磁化機(jī)的設(shè)計(jì)
    自制磁性螺絲刀
    超高商業(yè)大廈內(nèi)部磁場(chǎng)強(qiáng)度的研究與分析
    科技資訊(2017年3期)2017-03-25 08:54:13
    磁性離子交換樹(shù)脂的制備及其對(duì)Cr3+的吸附
    基于超拉普拉斯分布的磁化率重建算法
    巖(礦)石標(biāo)本磁化率測(cè)定方法試驗(yàn)及認(rèn)識(shí)
    一種新型磁性指紋刷的構(gòu)思
    淺談對(duì)磁場(chǎng)強(qiáng)度H和磁感應(yīng)強(qiáng)度B的認(rèn)識(shí)
    溫度對(duì)不同初始狀態(tài)ising模型磁化強(qiáng)度和磁化率的影響
    河南科技(2014年18期)2014-02-27 14:15:23
    免费观看a级毛片全部| 天美传媒精品一区二区| 女性被躁到高潮视频| 自线自在国产av| 欧美人与善性xxx| 欧美精品一区二区大全| 亚州av有码| 在线观看免费高清a一片| 国产69精品久久久久777片| 日韩精品免费视频一区二区三区 | 夜夜骑夜夜射夜夜干| 自线自在国产av| 亚洲国产日韩一区二区| 国产精品成人在线| 美女xxoo啪啪120秒动态图| 婷婷色综合www| 久久久久久久久久久丰满| 欧美人与善性xxx| 免费高清在线观看视频在线观看| 国产精品免费大片| 水蜜桃什么品种好| 免费看不卡的av| 天天操日日干夜夜撸| 精品一区二区免费观看| 午夜福利视频精品| 乱人伦中国视频| 亚洲欧美日韩东京热| 久久99精品国语久久久| 国产亚洲最大av| 简卡轻食公司| 一级毛片aaaaaa免费看小| 夜夜看夜夜爽夜夜摸| www.色视频.com| 黄色毛片三级朝国网站 | 一区二区av电影网| 丰满少妇做爰视频| 高清在线视频一区二区三区| 波野结衣二区三区在线| 色吧在线观看| 成年人午夜在线观看视频| 欧美精品高潮呻吟av久久| 亚洲中文av在线| 精品久久久久久久久亚洲| 久久综合国产亚洲精品| 特大巨黑吊av在线直播| 三级经典国产精品| 日韩大片免费观看网站| 亚洲精品日本国产第一区| tube8黄色片| 黑丝袜美女国产一区| 久久精品久久久久久噜噜老黄| 嫩草影院新地址| 婷婷色av中文字幕| 国产精品.久久久| 69精品国产乱码久久久| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 中文字幕人妻丝袜制服| 精品久久久久久久久av| 一级黄片播放器| 久久久欧美国产精品| a级毛片在线看网站| 爱豆传媒免费全集在线观看| 色婷婷久久久亚洲欧美| 亚洲图色成人| 亚洲电影在线观看av| 日本vs欧美在线观看视频 | 免费看光身美女| 自拍偷自拍亚洲精品老妇| 91久久精品国产一区二区三区| 国产精品久久久久久精品电影小说| 精品亚洲成国产av| 亚洲电影在线观看av| 国产精品国产三级国产av玫瑰| 午夜91福利影院| 国产爽快片一区二区三区| 日韩,欧美,国产一区二区三区| 免费高清在线观看视频在线观看| 青青草视频在线视频观看| 少妇 在线观看| a级一级毛片免费在线观看| 欧美三级亚洲精品| √禁漫天堂资源中文www| 女性被躁到高潮视频| 欧美97在线视频| 亚洲电影在线观看av| 亚洲一级一片aⅴ在线观看| 午夜老司机福利剧场| 乱人伦中国视频| 欧美最新免费一区二区三区| 日韩三级伦理在线观看| 欧美日韩av久久| 欧美97在线视频| 国产女主播在线喷水免费视频网站| 各种免费的搞黄视频| 我要看黄色一级片免费的| 免费黄频网站在线观看国产| 国产成人精品无人区| 777米奇影视久久| 亚洲四区av| 免费少妇av软件| 在线看a的网站| 精品午夜福利在线看| 国内精品宾馆在线| 国产午夜精品一二区理论片| 亚洲精品久久久久久婷婷小说| 免费久久久久久久精品成人欧美视频 | 18禁裸乳无遮挡动漫免费视频| 午夜福利影视在线免费观看| 一级毛片aaaaaa免费看小| 午夜福利网站1000一区二区三区| 亚洲欧美日韩卡通动漫| 十八禁高潮呻吟视频 | av国产精品久久久久影院| 丝袜在线中文字幕| 精品熟女少妇av免费看| 亚洲内射少妇av| 国产黄色视频一区二区在线观看| 国产欧美亚洲国产| 一级毛片 在线播放| 亚洲欧美清纯卡通| 大又大粗又爽又黄少妇毛片口| a级毛片在线看网站| 99九九线精品视频在线观看视频| 亚洲国产精品国产精品| 久久免费观看电影| 免费观看无遮挡的男女| 久久精品熟女亚洲av麻豆精品| 亚洲精品乱码久久久久久按摩| 美女xxoo啪啪120秒动态图| 久久久久久久大尺度免费视频| h日本视频在线播放| 一级毛片 在线播放| 亚洲国产色片| 久久鲁丝午夜福利片| 亚洲av成人精品一二三区| 国产在视频线精品| 亚洲人成网站在线观看播放| 久久精品久久久久久噜噜老黄| 插阴视频在线观看视频| 精品视频人人做人人爽| 26uuu在线亚洲综合色| 能在线免费看毛片的网站| 黄色欧美视频在线观看| 国产真实伦视频高清在线观看| 啦啦啦视频在线资源免费观看| 精品人妻一区二区三区麻豆| 在线观看av片永久免费下载| 日本av手机在线免费观看| 午夜精品国产一区二区电影| 午夜福利网站1000一区二区三区| √禁漫天堂资源中文www| 日韩亚洲欧美综合| 精品少妇内射三级| 在线观看av片永久免费下载| 这个男人来自地球电影免费观看 | 我要看日韩黄色一级片| 大香蕉97超碰在线| 亚洲国产精品国产精品| 亚洲,欧美,日韩| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| 99久国产av精品国产电影| 久久精品国产a三级三级三级| 久久狼人影院| 久久久久久久久久人人人人人人| 欧美精品人与动牲交sv欧美| 国产69精品久久久久777片| 亚洲av不卡在线观看| 十八禁高潮呻吟视频 | 国产精品无大码| 各种免费的搞黄视频| 久久午夜综合久久蜜桃| a 毛片基地| 精品久久久噜噜| 成人亚洲欧美一区二区av| 精华霜和精华液先用哪个| 成人影院久久| 熟女人妻精品中文字幕| 国产色爽女视频免费观看| 久久女婷五月综合色啪小说| 国产一区二区三区av在线| 成人影院久久| 久久99热6这里只有精品| 三级经典国产精品| 日韩av免费高清视频| 热re99久久精品国产66热6| 波野结衣二区三区在线| 久久影院123| 黑人巨大精品欧美一区二区蜜桃 | 波野结衣二区三区在线| 亚洲国产精品一区三区| 丁香六月天网| 日本午夜av视频| 欧美 亚洲 国产 日韩一| 成年女人在线观看亚洲视频| 哪个播放器可以免费观看大片| 日韩免费高清中文字幕av| 午夜免费男女啪啪视频观看| 婷婷色综合大香蕉| a 毛片基地| 国产成人freesex在线| 国产成人91sexporn| 国产精品久久久久久精品电影小说| 人人妻人人爽人人添夜夜欢视频 | 久久99热6这里只有精品| 成人亚洲欧美一区二区av| 国产精品三级大全| 午夜激情福利司机影院| 妹子高潮喷水视频| 日韩免费高清中文字幕av| 国产黄色视频一区二区在线观看| 亚洲不卡免费看| 久久久a久久爽久久v久久| 伦理电影大哥的女人| 精品国产一区二区久久| 亚洲国产最新在线播放| 麻豆乱淫一区二区| 少妇被粗大的猛进出69影院 | 欧美bdsm另类| 噜噜噜噜噜久久久久久91| 日本黄色片子视频| 欧美日本中文国产一区发布| 亚洲国产精品一区二区三区在线| 日韩精品有码人妻一区| av国产精品久久久久影院| 妹子高潮喷水视频| 十八禁网站网址无遮挡 | 女的被弄到高潮叫床怎么办| 内射极品少妇av片p| 成人黄色视频免费在线看| 能在线免费看毛片的网站| 美女中出高潮动态图| √禁漫天堂资源中文www| 久久久久久久大尺度免费视频| 我的老师免费观看完整版| 日日爽夜夜爽网站| 午夜福利,免费看| 国产精品人妻久久久影院| 另类精品久久| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 九草在线视频观看| 久久影院123| 国产日韩欧美在线精品| www.色视频.com| 天堂俺去俺来也www色官网| 你懂的网址亚洲精品在线观看| 9色porny在线观看| 久久久久久久久久久丰满| 免费人妻精品一区二区三区视频| 人妻少妇偷人精品九色| 肉色欧美久久久久久久蜜桃| 热re99久久国产66热| 两个人免费观看高清视频 | 秋霞伦理黄片| 日韩中文字幕视频在线看片| 亚洲国产精品国产精品| 国产欧美日韩综合在线一区二区 | 天堂中文最新版在线下载| 国产日韩一区二区三区精品不卡 | 免费看不卡的av| 女人精品久久久久毛片| 亚洲人成网站在线播| 一级,二级,三级黄色视频| 成人无遮挡网站| 欧美+日韩+精品| 亚洲内射少妇av| 日日啪夜夜爽| 久久精品久久久久久噜噜老黄| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| 国模一区二区三区四区视频| 黄色一级大片看看| 综合色丁香网| 久久av网站| 国产男女内射视频| 黄色欧美视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久久久免| 日本与韩国留学比较| 久久久久网色| 午夜福利在线观看免费完整高清在| 亚洲,一卡二卡三卡| 午夜视频国产福利| 九色成人免费人妻av| 精品亚洲乱码少妇综合久久| 免费播放大片免费观看视频在线观看| 久久久久精品性色| 99re6热这里在线精品视频| 欧美日韩在线观看h| 老熟女久久久| 九草在线视频观看| 99热网站在线观看| 18禁裸乳无遮挡动漫免费视频| 一级av片app| 色94色欧美一区二区| 精品一区在线观看国产| 日韩一区二区三区影片| 色5月婷婷丁香| 日日啪夜夜撸| 亚洲高清免费不卡视频| 日韩欧美一区视频在线观看 | 最近中文字幕高清免费大全6| 在线观看一区二区三区激情| 久久99精品国语久久久| 亚洲,一卡二卡三卡| 中文字幕亚洲精品专区| 成人亚洲欧美一区二区av| 夫妻午夜视频| 国产精品伦人一区二区| 欧美日韩综合久久久久久| 国产真实伦视频高清在线观看| 国产av国产精品国产| 丝袜喷水一区| 久久免费观看电影| 亚洲情色 制服丝袜| 亚洲精品456在线播放app| 欧美最新免费一区二区三区| 五月玫瑰六月丁香| 夜夜骑夜夜射夜夜干| 亚洲av不卡在线观看| 国产探花极品一区二区| 少妇人妻久久综合中文| av不卡在线播放| 赤兔流量卡办理| 国产一区二区三区av在线| 99九九在线精品视频 | 97超视频在线观看视频| 国产美女午夜福利| 久久久久国产网址| 99久久中文字幕三级久久日本| 亚洲欧美一区二区三区国产| 国产成人freesex在线| 中国三级夫妇交换| 狂野欧美白嫩少妇大欣赏| 黄色一级大片看看| 又粗又硬又长又爽又黄的视频| 天堂俺去俺来也www色官网| 久久久久久久久久人人人人人人| 午夜久久久在线观看| 久久精品国产a三级三级三级| 亚洲熟女精品中文字幕| 晚上一个人看的免费电影| 欧美亚洲 丝袜 人妻 在线| 黄色视频在线播放观看不卡| 国产精品一区二区性色av| 亚洲欧美清纯卡通| 欧美三级亚洲精品| 国产片特级美女逼逼视频| 全区人妻精品视频| 性色avwww在线观看| 国产高清不卡午夜福利| 亚洲国产最新在线播放| 免费av不卡在线播放| 91精品国产国语对白视频| 内射极品少妇av片p| 久久影院123| 亚洲精华国产精华液的使用体验| 亚洲内射少妇av| 又黄又爽又刺激的免费视频.| 精品国产一区二区久久| 久久精品国产鲁丝片午夜精品| 国产欧美日韩综合在线一区二区 | 2022亚洲国产成人精品| 三级国产精品片| 国内少妇人妻偷人精品xxx网站| 在现免费观看毛片| 久久热精品热| 亚洲精品第二区| 日韩人妻高清精品专区| 国产精品99久久99久久久不卡 | 免费观看无遮挡的男女| 老司机影院成人| 国产精品99久久久久久久久| 久久久久久久久大av| 下体分泌物呈黄色| 高清视频免费观看一区二区| 亚洲av男天堂| 韩国高清视频一区二区三区| 九九在线视频观看精品| 亚洲情色 制服丝袜| 成人美女网站在线观看视频| 91aial.com中文字幕在线观看| 中文字幕精品免费在线观看视频 | 精品国产国语对白av| 亚洲精华国产精华液的使用体验| 成人18禁高潮啪啪吃奶动态图 | 免费av不卡在线播放| 国产成人a∨麻豆精品| 国产精品.久久久| 啦啦啦视频在线资源免费观看| 黄色一级大片看看| 美女内射精品一级片tv| 黄色日韩在线| 我要看黄色一级片免费的| 日本免费在线观看一区| 国产老妇伦熟女老妇高清| 三级国产精品欧美在线观看| 亚洲国产色片| 日韩人妻高清精品专区| 热re99久久精品国产66热6| 街头女战士在线观看网站| 亚洲精品中文字幕在线视频 | 亚洲国产精品成人久久小说| 在线观看av片永久免费下载| 免费在线观看成人毛片| 国产永久视频网站| 亚洲国产色片| 免费人妻精品一区二区三区视频| 99热这里只有是精品在线观看| 国产极品天堂在线| 免费看日本二区| 深夜a级毛片| 极品教师在线视频| 国产欧美日韩综合在线一区二区 | 亚洲精品亚洲一区二区| 99热全是精品| 国产黄色免费在线视频| 丝袜在线中文字幕| 久久久久精品久久久久真实原创| 日本91视频免费播放| 国产精品人妻久久久影院| 亚洲丝袜综合中文字幕| 在线观看国产h片| 熟女av电影| 精品少妇久久久久久888优播| 人妻制服诱惑在线中文字幕| 一级毛片我不卡| 肉色欧美久久久久久久蜜桃| 最近2019中文字幕mv第一页| 国产精品福利在线免费观看| 亚洲精品国产av蜜桃| 欧美bdsm另类| 中文在线观看免费www的网站| 亚洲av福利一区| 国产欧美日韩精品一区二区| 天堂8中文在线网| 久久久久国产网址| 日韩免费高清中文字幕av| 边亲边吃奶的免费视频| 啦啦啦啦在线视频资源| 少妇人妻 视频| 久久久久久久久久久丰满| 亚洲自偷自拍三级| kizo精华| 91精品伊人久久大香线蕉| 成人亚洲欧美一区二区av| 老司机亚洲免费影院| 两个人的视频大全免费| 天堂俺去俺来也www色官网| 91aial.com中文字幕在线观看| 久久久久久人妻| 桃花免费在线播放| 午夜免费观看性视频| 国产高清国产精品国产三级| 黄色配什么色好看| 卡戴珊不雅视频在线播放| 99视频精品全部免费 在线| 一本一本综合久久| 一区二区三区精品91| 91精品国产九色| 久久精品国产a三级三级三级| 日韩中字成人| 高清在线视频一区二区三区| 精品久久国产蜜桃| 麻豆成人av视频| 久久韩国三级中文字幕| 亚洲欧美精品自产自拍| 久久人人爽人人片av| 免费观看a级毛片全部| 国产乱人偷精品视频| 亚洲国产欧美在线一区| 久久久欧美国产精品| 青青草视频在线视频观看| 欧美日本中文国产一区发布| 亚洲va在线va天堂va国产| 交换朋友夫妻互换小说| 日本黄色片子视频| 日本wwww免费看| 国产成人精品一,二区| 丰满乱子伦码专区| 黄色配什么色好看| 亚洲综合精品二区| 91午夜精品亚洲一区二区三区| 成人二区视频| 七月丁香在线播放| 免费看不卡的av| 日本av手机在线免费观看| 色94色欧美一区二区| 在线 av 中文字幕| 蜜臀久久99精品久久宅男| 亚洲av不卡在线观看| 99久久中文字幕三级久久日本| 国产在线男女| 啦啦啦中文免费视频观看日本| 免费av不卡在线播放| a 毛片基地| 国产一区二区三区av在线| 涩涩av久久男人的天堂| 激情五月婷婷亚洲| 成人亚洲精品一区在线观看| 国产精品嫩草影院av在线观看| 亚洲av中文av极速乱| 亚洲av不卡在线观看| 亚洲伊人久久精品综合| 少妇的逼水好多| 久久久久久久精品精品| 永久网站在线| 国产精品国产三级国产专区5o| 久久久久精品久久久久真实原创| 亚洲第一av免费看| 一级二级三级毛片免费看| 欧美最新免费一区二区三区| 大香蕉97超碰在线| 深夜a级毛片| 亚洲精品一区蜜桃| av一本久久久久| 在线精品无人区一区二区三| 国模一区二区三区四区视频| 日日啪夜夜撸| 自线自在国产av| 有码 亚洲区| 精品少妇黑人巨大在线播放| 春色校园在线视频观看| 这个男人来自地球电影免费观看 | 岛国毛片在线播放| 在线 av 中文字幕| 国产精品成人在线| 少妇精品久久久久久久| 国产日韩一区二区三区精品不卡 | 国产女主播在线喷水免费视频网站| 丰满饥渴人妻一区二区三| 草草在线视频免费看| av视频免费观看在线观看| 亚洲欧美日韩卡通动漫| 国产免费一级a男人的天堂| 在线观看www视频免费| 亚洲综合色惰| 国产伦理片在线播放av一区| 欧美日韩国产mv在线观看视频| 日韩av不卡免费在线播放| 久热久热在线精品观看| 久久久国产一区二区| 99国产精品免费福利视频| 我要看黄色一级片免费的| 最新的欧美精品一区二区| 涩涩av久久男人的天堂| 我要看日韩黄色一级片| 国产av精品麻豆| 99热6这里只有精品| 亚洲欧美成人综合另类久久久| 午夜视频国产福利| 国产免费视频播放在线视频| 久久人妻熟女aⅴ| 国产片特级美女逼逼视频| 日韩欧美 国产精品| av国产久精品久网站免费入址| 人妻制服诱惑在线中文字幕| 国产精品女同一区二区软件| 亚洲欧美成人精品一区二区| 久久久久人妻精品一区果冻| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 午夜视频国产福利| 成人18禁高潮啪啪吃奶动态图 | 久久久久久人妻| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一区蜜桃| xxx大片免费视频| 青春草视频在线免费观看| 少妇的逼好多水| 午夜免费观看性视频| 9色porny在线观看| 简卡轻食公司| 免费高清在线观看视频在线观看| 妹子高潮喷水视频| 欧美日韩视频精品一区| 高清毛片免费看| 亚洲国产欧美在线一区| 99国产精品免费福利视频| 国产片特级美女逼逼视频| 精品国产乱码久久久久久小说| 80岁老熟妇乱子伦牲交| 99久久综合免费| 男女边摸边吃奶| 国产亚洲av片在线观看秒播厂| 久久影院123| 日本黄大片高清| 国产精品久久久久久久久免| 国产精品欧美亚洲77777| 少妇高潮的动态图| 韩国av在线不卡| www.色视频.com| 欧美激情国产日韩精品一区| 夜夜骑夜夜射夜夜干| 成人黄色视频免费在线看| 久久这里有精品视频免费| 18禁在线无遮挡免费观看视频| 久久99蜜桃精品久久| 精品久久久噜噜| 男男h啪啪无遮挡| 日日啪夜夜撸| 免费少妇av软件| 久久 成人 亚洲| 2022亚洲国产成人精品| 又爽又黄a免费视频| 久久久久久久精品精品| 91久久精品电影网| 国产日韩欧美亚洲二区| 亚洲av电影在线观看一区二区三区| 在线看a的网站| 内地一区二区视频在线| 美女xxoo啪啪120秒动态图| 国产免费福利视频在线观看|