• <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
    一级,二级,三级黄色视频| 美女福利国产在线| 9191精品国产免费久久| 韩国精品一区二区三区| av在线老鸭窝| 黄频高清免费视频| 久久久久视频综合| 丰满人妻熟妇乱又伦精品不卡| 免费不卡黄色视频| 青春草亚洲视频在线观看| 热re99久久国产66热| a 毛片基地| 亚洲欧美激情在线| 欧美激情 高清一区二区三区| cao死你这个sao货| 欧美人与性动交α欧美精品济南到| 伦理电影免费视频| 欧美国产精品va在线观看不卡| 老汉色∧v一级毛片| 男人舔女人的私密视频| 成年人免费黄色播放视频| 欧美精品人与动牲交sv欧美| 国产一区二区三区av在线| 最近最新中文字幕大全免费视频| 一区二区三区精品91| 久久精品久久久久久噜噜老黄| 黄色怎么调成土黄色| 黑人巨大精品欧美一区二区mp4| 日韩电影二区| 母亲3免费完整高清在线观看| 91麻豆av在线| 在线 av 中文字幕| 久热这里只有精品99| 久久99一区二区三区| 成人国产av品久久久| 一区二区三区精品91| 18禁黄网站禁片午夜丰满| 91老司机精品| 女人被躁到高潮嗷嗷叫费观| 视频区欧美日本亚洲| 免费在线观看黄色视频的| 色综合欧美亚洲国产小说| 久久午夜综合久久蜜桃| 亚洲成人免费电影在线观看| 午夜福利,免费看| 超色免费av| 国产精品欧美亚洲77777| 搡老乐熟女国产| 9色porny在线观看| 一本—道久久a久久精品蜜桃钙片| 一级片'在线观看视频| 正在播放国产对白刺激| 窝窝影院91人妻| 我要看黄色一级片免费的| 日本五十路高清| 老司机影院毛片| 久久这里只有精品19| 丝袜美足系列| 欧美精品高潮呻吟av久久| 黄频高清免费视频| 欧美中文综合在线视频| 亚洲欧美精品综合一区二区三区| 精品一品国产午夜福利视频| 又大又爽又粗| 91老司机精品| 国产深夜福利视频在线观看| 亚洲av日韩精品久久久久久密| 免费在线观看黄色视频的| 真人做人爱边吃奶动态| 丝袜喷水一区| 欧美激情高清一区二区三区| 国产极品粉嫩免费观看在线| 欧美另类一区| 王馨瑶露胸无遮挡在线观看| 国产精品国产三级国产专区5o| 欧美精品高潮呻吟av久久| 男女之事视频高清在线观看| 人人妻人人澡人人看| 国产片内射在线| 午夜91福利影院| 亚洲国产精品成人久久小说| 午夜两性在线视频| 日韩制服丝袜自拍偷拍| 老汉色av国产亚洲站长工具| 免费在线观看黄色视频的| 国产精品欧美亚洲77777| 精品卡一卡二卡四卡免费| 下体分泌物呈黄色| 老熟妇仑乱视频hdxx| 在线观看免费视频网站a站| 国产一区二区三区综合在线观看| 性高湖久久久久久久久免费观看| 精品福利永久在线观看| 亚洲人成电影免费在线| 法律面前人人平等表现在哪些方面 | 亚洲五月色婷婷综合| 亚洲全国av大片| 国产真人三级小视频在线观看| 日韩,欧美,国产一区二区三区| 亚洲av欧美aⅴ国产| 日本一区二区免费在线视频| 男女高潮啪啪啪动态图| a级毛片黄视频| 另类亚洲欧美激情| 欧美日韩一级在线毛片| 无限看片的www在线观看| 国产在线观看jvid| 男女边摸边吃奶| 天堂8中文在线网| 在线亚洲精品国产二区图片欧美| 欧美xxⅹ黑人| 精品一区二区三区四区五区乱码| 国产高清国产精品国产三级| 中文精品一卡2卡3卡4更新| 黄片播放在线免费| 男人爽女人下面视频在线观看| 久久久久久人人人人人| 黄色a级毛片大全视频| 高清欧美精品videossex| 国产真人三级小视频在线观看| 日本a在线网址| 久久久国产一区二区| a级毛片黄视频| 国产一卡二卡三卡精品| 最黄视频免费看| 午夜日韩欧美国产| 天天影视国产精品| 亚洲精品日韩在线中文字幕| 免费一级毛片在线播放高清视频 | 亚洲精品粉嫩美女一区| 首页视频小说图片口味搜索| 亚洲欧美日韩另类电影网站| 国产亚洲欧美在线一区二区| 亚洲综合色网址| 下体分泌物呈黄色| 在线精品无人区一区二区三| 狂野欧美激情性xxxx| 国产精品一区二区精品视频观看| 国产在线免费精品| 91麻豆av在线| 丝袜在线中文字幕| 国产精品av久久久久免费| 大片电影免费在线观看免费| 欧美日韩黄片免| 国产精品.久久久| 老司机午夜十八禁免费视频| 亚洲精品第二区| 久久久久精品人妻al黑| 精品久久久久久电影网| 国产精品免费视频内射| 免费高清在线观看日韩| 国产高清国产精品国产三级| 国产野战对白在线观看| a级毛片在线看网站| 亚洲精品久久午夜乱码| 青春草视频在线免费观看| 久久性视频一级片| 亚洲欧美清纯卡通| 91麻豆精品激情在线观看国产 | 国产又爽黄色视频| 男女边摸边吃奶| 国精品久久久久久国模美| 在线观看人妻少妇| 亚洲美女黄色视频免费看| 男男h啪啪无遮挡| 久久国产精品男人的天堂亚洲| 亚洲精品日韩在线中文字幕| 精品少妇内射三级| 国产欧美日韩一区二区三 | 丰满饥渴人妻一区二区三| 岛国毛片在线播放| 欧美中文综合在线视频| 成人av一区二区三区在线看 | 建设人人有责人人尽责人人享有的| 成年人免费黄色播放视频| 18禁裸乳无遮挡动漫免费视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品国产色婷婷电影| 狠狠精品人妻久久久久久综合| 亚洲精品中文字幕一二三四区 | 动漫黄色视频在线观看| 国产免费福利视频在线观看| 亚洲国产精品一区三区| 91字幕亚洲| 永久免费av网站大全| 人人澡人人妻人| 天天躁夜夜躁狠狠躁躁| 久久精品成人免费网站| 高潮久久久久久久久久久不卡| 亚洲国产成人一精品久久久| 久久精品人人爽人人爽视色| 香蕉国产在线看| 啦啦啦视频在线资源免费观看| 午夜福利视频在线观看免费| 最黄视频免费看| 免费日韩欧美在线观看| 欧美日本中文国产一区发布| 成年女人毛片免费观看观看9 | 精品一区二区三区av网在线观看 | 午夜福利一区二区在线看| 久久香蕉激情| 秋霞在线观看毛片| 免费观看av网站的网址| 真人做人爱边吃奶动态| 日韩欧美一区视频在线观看| 女警被强在线播放| 亚洲三区欧美一区| 大型av网站在线播放| 国产成人av激情在线播放| 国产野战对白在线观看| av福利片在线| 亚洲欧美激情在线| 国产亚洲欧美在线一区二区| 后天国语完整版免费观看| 高潮久久久久久久久久久不卡| 97人妻天天添夜夜摸| 国产日韩欧美在线精品| 国产精品av久久久久免费| 18禁黄网站禁片午夜丰满| 亚洲中文日韩欧美视频| 丝瓜视频免费看黄片| 十八禁网站免费在线| 国产av又大| 高清在线国产一区| 亚洲熟女毛片儿| 蜜桃在线观看..| av超薄肉色丝袜交足视频| 成年女人毛片免费观看观看9 | 国产成人精品在线电影| 日本av手机在线免费观看| 男女午夜视频在线观看| 欧美人与性动交α欧美精品济南到| a级毛片黄视频| 女性生殖器流出的白浆| 欧美精品av麻豆av| 激情视频va一区二区三区| 久久亚洲精品不卡| 男女边摸边吃奶| 日日爽夜夜爽网站| 精品乱码久久久久久99久播| 精品欧美一区二区三区在线| 每晚都被弄得嗷嗷叫到高潮| 又大又爽又粗| 国产成+人综合+亚洲专区| 亚洲av美国av| 一级毛片电影观看| 日本五十路高清| 成年人黄色毛片网站| 国产精品熟女久久久久浪| 亚洲中文av在线| √禁漫天堂资源中文www| 日本五十路高清| 制服诱惑二区| 国产精品国产三级国产专区5o| 男女之事视频高清在线观看| 成人影院久久| 青春草视频在线免费观看| 国产成人av教育| 亚洲三区欧美一区| 亚洲精品日韩在线中文字幕| 手机成人av网站| 国产av精品麻豆| 高潮久久久久久久久久久不卡| avwww免费| 日本撒尿小便嘘嘘汇集6| 性色av一级| 亚洲av日韩在线播放| 中文字幕高清在线视频| 亚洲精品粉嫩美女一区| 精品少妇久久久久久888优播| 国产淫语在线视频| 亚洲欧美成人综合另类久久久| 久久国产精品影院| 男女下面插进去视频免费观看| 精品卡一卡二卡四卡免费| 99re6热这里在线精品视频| 午夜成年电影在线免费观看| 久久狼人影院| 精品久久蜜臀av无| 丝袜喷水一区| 国产区一区二久久| 精品人妻1区二区| 久9热在线精品视频| 侵犯人妻中文字幕一二三四区| 精品国产乱码久久久久久男人| 啦啦啦视频在线资源免费观看| 宅男免费午夜| 欧美另类亚洲清纯唯美| 啦啦啦免费观看视频1| 亚洲欧美成人综合另类久久久| 精品国产超薄肉色丝袜足j| 国产成人精品无人区| 午夜91福利影院| 亚洲精品国产精品久久久不卡| 91大片在线观看| 国产淫语在线视频| 欧美另类亚洲清纯唯美| 国产淫语在线视频| 熟女少妇亚洲综合色aaa.| 国产日韩欧美亚洲二区| 欧美日韩视频精品一区| 不卡av一区二区三区| 国产欧美日韩精品亚洲av| 亚洲免费av在线视频| 欧美日韩av久久| 亚洲男人天堂网一区| 国产国语露脸激情在线看| 99国产精品一区二区三区| 国产一卡二卡三卡精品| 免费在线观看黄色视频的| 热99国产精品久久久久久7| 宅男免费午夜| 操美女的视频在线观看| 午夜两性在线视频| 亚洲视频免费观看视频| 97在线人人人人妻| 欧美精品高潮呻吟av久久| 亚洲人成电影免费在线| 国产精品熟女久久久久浪| 午夜福利,免费看| 操美女的视频在线观看| 妹子高潮喷水视频| 狠狠婷婷综合久久久久久88av| 亚洲 欧美一区二区三区| 国产真人三级小视频在线观看| 日韩视频一区二区在线观看| 一本色道久久久久久精品综合| 夜夜骑夜夜射夜夜干| 免费在线观看完整版高清| 欧美日韩福利视频一区二区| 黑人操中国人逼视频| 日韩电影二区| 又大又爽又粗| 少妇猛男粗大的猛烈进出视频| 午夜影院在线不卡| 成人手机av| 欧美日韩国产mv在线观看视频| 国产一区有黄有色的免费视频| 午夜福利在线观看吧| 50天的宝宝边吃奶边哭怎么回事| 久久久国产欧美日韩av| 在线观看免费日韩欧美大片| 午夜日韩欧美国产| 日韩精品免费视频一区二区三区| 亚洲天堂av无毛| 午夜福利影视在线免费观看| 久久久精品免费免费高清| 久9热在线精品视频| 国产精品1区2区在线观看. | 天堂8中文在线网| 免费黄频网站在线观看国产| 久久性视频一级片| 欧美乱码精品一区二区三区| 亚洲视频免费观看视频| 亚洲伊人久久精品综合| 国产欧美日韩一区二区三区在线| 日韩有码中文字幕| 久久久久久久国产电影| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 亚洲色图 男人天堂 中文字幕| 啦啦啦啦在线视频资源| 国产av又大| 蜜桃国产av成人99| 亚洲人成电影免费在线| av线在线观看网站| 两人在一起打扑克的视频| av线在线观看网站| 一本—道久久a久久精品蜜桃钙片| 国产亚洲av片在线观看秒播厂| 亚洲av男天堂| 在线观看一区二区三区激情| 午夜激情久久久久久久| 少妇被粗大的猛进出69影院| 久久国产精品人妻蜜桃| 亚洲成av片中文字幕在线观看| 黄色a级毛片大全视频| 99re6热这里在线精品视频| 午夜免费鲁丝| 久久久国产精品麻豆| 99国产精品99久久久久| 亚洲九九香蕉| 色老头精品视频在线观看| 交换朋友夫妻互换小说| 又紧又爽又黄一区二区| 成年人黄色毛片网站| 黄色 视频免费看| 成人黄色视频免费在线看| 另类精品久久| bbb黄色大片| 老熟女久久久| 日韩大片免费观看网站| 啦啦啦免费观看视频1| 2018国产大陆天天弄谢| 91麻豆精品激情在线观看国产 | 午夜福利在线免费观看网站| 老司机深夜福利视频在线观看 | 丰满饥渴人妻一区二区三| 两性夫妻黄色片| 12—13女人毛片做爰片一| 亚洲一卡2卡3卡4卡5卡精品中文| 99久久人妻综合| 极品人妻少妇av视频| 久久久国产欧美日韩av| 美女扒开内裤让男人捅视频| 国产又色又爽无遮挡免| 一级毛片女人18水好多| 日韩视频一区二区在线观看| 五月天丁香电影| av网站免费在线观看视频| 欧美日韩av久久| 国产黄色免费在线视频| 悠悠久久av| 精品熟女少妇八av免费久了| videos熟女内射| www.精华液| 两个人免费观看高清视频| 久久午夜综合久久蜜桃| 老汉色av国产亚洲站长工具| 99久久99久久久精品蜜桃| 人人妻人人澡人人爽人人夜夜| 天堂8中文在线网| 日韩精品免费视频一区二区三区| 亚洲国产av影院在线观看| 秋霞在线观看毛片| 国产亚洲精品久久久久5区| 亚洲精品av麻豆狂野| 日韩欧美一区视频在线观看| 日韩 欧美 亚洲 中文字幕| 日本vs欧美在线观看视频| 国产成人系列免费观看| 久久久精品94久久精品| a级片在线免费高清观看视频| 国产激情久久老熟女| 国产一区有黄有色的免费视频| 日本wwww免费看| 国产精品久久久久久精品古装| 69精品国产乱码久久久| 精品亚洲乱码少妇综合久久| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三 | 久久精品熟女亚洲av麻豆精品| 男女国产视频网站| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲 | 9色porny在线观看| 成人亚洲精品一区在线观看| 免费日韩欧美在线观看| 超色免费av| 午夜精品国产一区二区电影| 可以免费在线观看a视频的电影网站| 午夜免费观看性视频| 国产人伦9x9x在线观看| 日韩三级视频一区二区三区| 日韩免费高清中文字幕av| 欧美精品人与动牲交sv欧美| 丝袜美足系列| 黄色视频,在线免费观看| 91九色精品人成在线观看| 欧美xxⅹ黑人| 亚洲欧美日韩另类电影网站| 一进一出抽搐动态| 国产又爽黄色视频| 黄片播放在线免费| 美女脱内裤让男人舔精品视频| 老司机亚洲免费影院| a在线观看视频网站| 亚洲国产成人一精品久久久| 老司机午夜福利在线观看视频 | 国产精品.久久久| 亚洲自偷自拍图片 自拍| 麻豆av在线久日| 午夜视频精品福利| 久久久久国产一级毛片高清牌| 在线观看一区二区三区激情| 亚洲国产欧美在线一区| 亚洲精品国产区一区二| av在线app专区| tocl精华| 热re99久久精品国产66热6| 亚洲精品自拍成人| 又紧又爽又黄一区二区| 亚洲欧美精品自产自拍| 天天躁日日躁夜夜躁夜夜| 黄色视频不卡| 久久人人爽av亚洲精品天堂| 国产xxxxx性猛交| 国产熟女午夜一区二区三区| 国产亚洲一区二区精品| 免费看十八禁软件| 另类精品久久| 国产精品亚洲av一区麻豆| 亚洲欧美精品自产自拍| 另类精品久久| 日韩大片免费观看网站| 亚洲,欧美精品.| 99久久99久久久精品蜜桃| 美女大奶头黄色视频| 婷婷成人精品国产| 精品久久蜜臀av无| 99国产精品免费福利视频| 国产在视频线精品| 免费人妻精品一区二区三区视频| 99久久人妻综合| 黑人猛操日本美女一级片| av有码第一页| 男人爽女人下面视频在线观看| 亚洲熟女精品中文字幕| 波多野结衣av一区二区av| 91大片在线观看| 天天添夜夜摸| 又紧又爽又黄一区二区| 少妇人妻久久综合中文| 国产欧美日韩一区二区三 | 涩涩av久久男人的天堂| 久久午夜综合久久蜜桃| 俄罗斯特黄特色一大片| 一区二区三区激情视频| 桃花免费在线播放| 免费在线观看完整版高清| 午夜91福利影院| 最新的欧美精品一区二区| 成人国产一区最新在线观看| 黑人猛操日本美女一级片| 五月开心婷婷网| 老司机在亚洲福利影院| 亚洲美女黄色视频免费看| tube8黄色片| 两个人看的免费小视频| 国产亚洲精品久久久久5区| 丰满饥渴人妻一区二区三| 国产97色在线日韩免费| 美女高潮喷水抽搐中文字幕| 中文字幕精品免费在线观看视频| 嫁个100分男人电影在线观看| 精品亚洲乱码少妇综合久久| 制服人妻中文乱码| 久久人妻福利社区极品人妻图片| 精品国产超薄肉色丝袜足j| 18禁裸乳无遮挡动漫免费视频| 美女高潮喷水抽搐中文字幕| 老司机深夜福利视频在线观看 | videos熟女内射| 精品第一国产精品| 欧美日韩亚洲综合一区二区三区_| 黄网站色视频无遮挡免费观看| 少妇人妻久久综合中文| 在线观看一区二区三区激情| 午夜久久久在线观看| 久久久精品免费免费高清| 国产日韩一区二区三区精品不卡| 成年人免费黄色播放视频| 亚洲精品国产av成人精品| 国产精品av久久久久免费| 久久久久网色| 免费人妻精品一区二区三区视频| 动漫黄色视频在线观看| 成人黄色视频免费在线看| 国产在线一区二区三区精| 99久久精品国产亚洲精品| 精品国产一区二区三区四区第35| √禁漫天堂资源中文www| 99国产精品免费福利视频| 精品少妇久久久久久888优播| 久久国产精品大桥未久av| 无限看片的www在线观看| 狠狠婷婷综合久久久久久88av| 久久精品亚洲熟妇少妇任你| 大型av网站在线播放| 国产日韩欧美视频二区| 在线观看免费午夜福利视频| 国产精品影院久久| 久久久久视频综合| 精品视频人人做人人爽| 午夜久久久在线观看| 深夜精品福利| 99国产精品99久久久久| 少妇裸体淫交视频免费看高清 | 人人妻人人澡人人看| 国产精品久久久av美女十八| 亚洲精品国产一区二区精华液| 黄色片一级片一级黄色片| 亚洲五月色婷婷综合| 午夜精品久久久久久毛片777| 伊人久久大香线蕉亚洲五| 亚洲伊人久久精品综合| 国产成人系列免费观看| 欧美精品av麻豆av| 高清视频免费观看一区二区| 日日爽夜夜爽网站| 美女脱内裤让男人舔精品视频| 女性被躁到高潮视频| 日韩大码丰满熟妇| 欧美 日韩 精品 国产| 一区二区av电影网| 亚洲第一av免费看| 丰满人妻熟妇乱又伦精品不卡| 国产精品二区激情视频| netflix在线观看网站| 啪啪无遮挡十八禁网站| 99热全是精品| 女警被强在线播放| 在线观看人妻少妇| av天堂久久9| 成人免费观看视频高清| 高清欧美精品videossex| 又大又爽又粗| 中文字幕制服av| 十分钟在线观看高清视频www| 午夜福利乱码中文字幕| 91老司机精品| 大陆偷拍与自拍| 精品福利观看| 69精品国产乱码久久久| 丝袜喷水一区|