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

    新型一體化脫硫工藝反應(yīng)器內(nèi)流場(chǎng)的數(shù)值模擬研究*

    2016-03-13 05:52:12楊竟憲張海濤
    環(huán)境污染與防治 2016年12期
    關(guān)鍵詞:氣力輸送塔內(nèi)流場(chǎng)

    楊竟憲 張海濤

    (華東理工大學(xué)化工學(xué)院,上海 200237)

    新型一體化脫硫工藝反應(yīng)器內(nèi)流場(chǎng)的數(shù)值模擬研究*

    楊竟憲 張海濤#

    (華東理工大學(xué)化工學(xué)院,上海 200237)

    基于顆粒動(dòng)力學(xué)理論,利用Fluent 軟件中的氣體-顆粒兩相流體模型,對(duì)新型一體化脫硫工藝(NID)反應(yīng)器中不同工況下的三維流場(chǎng)進(jìn)行數(shù)值模擬,并結(jié)合工廠實(shí)測(cè)數(shù)據(jù)及文獻(xiàn)數(shù)據(jù)對(duì)模擬結(jié)果進(jìn)行驗(yàn)證。結(jié)果表明,F(xiàn)luent軟件的模擬數(shù)據(jù)與工廠實(shí)測(cè)數(shù)據(jù)及文獻(xiàn)數(shù)據(jù)相吻合,利用Fluent軟件模擬了NID反應(yīng)器內(nèi)顆粒速度和濃度的分布,分析了反應(yīng)器內(nèi)流場(chǎng)的形成機(jī)制以及煙氣進(jìn)速和顆粒粒徑對(duì)塔內(nèi)壓降的影響。由模擬結(jié)果可知,在煙氣進(jìn)速為18 m/s、顆粒粒徑為0.5~50.0 μm時(shí),顆粒濃度分布均勻,塔內(nèi)回流區(qū)最小,壓強(qiáng)變化穩(wěn)定,最有利于提高實(shí)際工業(yè)運(yùn)行的穩(wěn)定性與安全性。

    NID 氣力輸送 數(shù)值模擬 壓降

    新型一體化脫硫工藝(NID)是在噴霧干燥法(半干法)脫硫工藝的基礎(chǔ)上發(fā)展而來[1],其脫硫塔裝置中的反應(yīng)器屬于矩形輸送床,煙氣進(jìn)入反應(yīng)器后帶動(dòng)循環(huán)顆粒和脫硫劑豎直向上運(yùn)動(dòng),顆粒稀相輸送與脫硫反應(yīng)同時(shí)進(jìn)行。

    在公開的資料中,有關(guān)NID運(yùn)行的報(bào)道均來自于工廠脫硫裝置,主要包括了脫硫塔整體壓降和進(jìn)出口溫度等運(yùn)行數(shù)據(jù)[2-5],尚未有文獻(xiàn)對(duì)脫硫塔內(nèi)流場(chǎng)進(jìn)行模擬分析。為提高實(shí)際工業(yè)運(yùn)行過程中設(shè)備的經(jīng)濟(jì)性和安全性,需要對(duì)NID脫硫塔內(nèi)反應(yīng)器的流場(chǎng)分布進(jìn)行模擬計(jì)算,根據(jù)反應(yīng)器內(nèi)顆粒速度場(chǎng)和濃度場(chǎng)的分布,分析壓強(qiáng)場(chǎng)的分布狀況與形成機(jī)制,從而得出優(yōu)化設(shè)備運(yùn)行安全性的條件。壓降分布越均勻,設(shè)備運(yùn)行越安全。脫硫塔內(nèi)反應(yīng)器的固體進(jìn)料中,99%(質(zhì)量分?jǐn)?shù))以上為循環(huán)顆粒,脫硫反應(yīng)對(duì)流場(chǎng)分布沒有任何影響,因此本研究主要對(duì)反應(yīng)器中循環(huán)顆粒氣力輸送的物理過程進(jìn)行計(jì)算分析。

    氣力輸送是能源、冶金和化工等行業(yè)普遍使用的輸送物料的方法,而目前氣力輸送的研究多為實(shí)驗(yàn)研究[6-8]。氣體速度和顆粒粒徑均會(huì)對(duì)顆粒輸送設(shè)備的流場(chǎng)產(chǎn)生影響,工廠運(yùn)行和實(shí)驗(yàn)室小試的操作過程中既不能精確控制顆粒粒徑等運(yùn)行參數(shù),也不能對(duì)塔內(nèi)氣固相流場(chǎng)細(xì)節(jié)進(jìn)行分析[9-12]。因此,計(jì)算流體力學(xué)(CFD)成為了工業(yè)反應(yīng)器系統(tǒng)設(shè)計(jì)和參數(shù)優(yōu)化的有力手段,主要計(jì)算工具為Fluent軟件。

    利用Fluent軟件對(duì)氣體-顆粒兩相流進(jìn)行數(shù)值模擬時(shí),主要使用Euler法和Lagrange法兩種模型。國(guó)內(nèi)外很多學(xué)者都利用Fluent軟件選擇合適的模型,并對(duì)類似的兩相流場(chǎng)進(jìn)行基礎(chǔ)研究。MAROCCO等[13]采用Euler-Lagrange法模擬濕法脫硫反應(yīng)器中的壓降和出口溫度,模擬結(jié)果與實(shí)際值相比壓降相差35.3 Pa,溫度僅相差0.74 K。GOMEZ等[14]模擬了脫硫反應(yīng)器內(nèi)不同工況下氣液固各相的速度以及反應(yīng)器內(nèi)等壓線,模擬結(jié)果和工廠數(shù)據(jù)相比誤差小于5%。以上文獻(xiàn)數(shù)據(jù)均顯示了數(shù)值模擬在流場(chǎng)計(jì)算中的準(zhǔn)確性。

    本研究針對(duì)實(shí)際運(yùn)行的工業(yè)NID反應(yīng)器內(nèi)流場(chǎng)狀況進(jìn)行研究,以顆粒動(dòng)力學(xué)理論為基礎(chǔ)建立兩流體模型,即Euler法,該模型能夠全面地考慮氣固相間作用、顆粒湍動(dòng)黏度和顆粒間的碰撞作用,對(duì)流場(chǎng)細(xì)節(jié)進(jìn)行準(zhǔn)確描述。在氣力輸送過程中,影響壓降分布的主要因素是氣固速度和顆粒濃度分布,其中反應(yīng)器內(nèi)速度的大小和方向均會(huì)對(duì)顆粒速度和濃度的分布產(chǎn)生影響,因此,本研究從煙氣進(jìn)速和循環(huán)顆粒粒徑兩個(gè)參數(shù)入手,考察不同工況下顆粒在反應(yīng)器內(nèi)的速度分布、濃度分布以及反應(yīng)器內(nèi)壓強(qiáng)分布。根據(jù)模擬得出的顆粒速度矢量圖、濃度分布圖和反應(yīng)器內(nèi)壓降分布圖進(jìn)行分析,了解壓降分布的形成機(jī)制,從而為設(shè)備運(yùn)行參數(shù)的優(yōu)化提供參考。

    1 研究?jī)?nèi)容

    1.1 研究對(duì)象

    研究對(duì)象為NID脫硫塔(半干式),其物理模型如圖1所示。由圖1可見,NID脫硫塔主要包括底部U形彎頭和矩形脫硫反應(yīng)器。底部的U形彎頭呈漸開設(shè)計(jì),內(nèi)彎與外彎的距離增加20%。反應(yīng)器是垂直煙道結(jié)構(gòu),從高爐出來的待處理煙氣通過U形彎頭,在反應(yīng)器入口處與循環(huán)物料混合并向上流動(dòng),循環(huán)物料表面帶有的少量水分與高溫?zé)煔饨佑|后蒸發(fā),煙氣溫度降低。煙氣出入口尺寸均為2.000 m×4.000 m,循環(huán)物料入口為0.500 m×4.000 m,反應(yīng)器高度為17.312 m。

    圖1 NID脫硫塔物理模型Fig.1 Physical model of NID desulfurization tower

    1.2 研究方法

    本研究采用Gambit作為前處理軟件進(jìn)行建模并劃分網(wǎng)格(見圖2),將NID脫硫塔分成3個(gè)相連的部分,并對(duì)其中的U形彎頭和反應(yīng)段連接部分進(jìn)行加密處理,總網(wǎng)格數(shù)為224 240個(gè)。然后將建立好的模型在Fluent 14.5中進(jìn)行計(jì)算和分析處理。

    圖2 NID脫硫塔網(wǎng)格劃分Fig.2 Grid of NID desulfurization tower

    1.3 基本假設(shè)

    由于反應(yīng)器內(nèi)的液相僅為顆粒表面的少量水分,氣液相間的傳熱傳質(zhì)僅對(duì)脫硫塔內(nèi)溫度產(chǎn)生影響,因此塔內(nèi)流場(chǎng)可認(rèn)為是氣固兩相流,并且在計(jì)算過程中作如下假設(shè)[15]:

    (1) 循環(huán)灰顆粒是直徑相同、密度相同的圓球顆粒;

    (2) 煙氣和蒸氣流動(dòng)是三維、不可壓縮的非定常流動(dòng);

    (3) 不考慮溫度、傳質(zhì)、傳熱、相變對(duì)湍流的影響;

    (4)煙氣流動(dòng)處于旺盛湍流區(qū);

    (5) 脫硫反應(yīng)器壁為絕熱,對(duì)流及換熱只在3相之間。

    1.4 參數(shù)設(shè)置

    1.4.1 求解參數(shù)

    采用非定常的壓力基隱式求解器,湍流模型為可實(shí)現(xiàn)k-ε方程,壓力速度耦合方式選用SIMPLE算法,松弛因子采用默認(rèn)值,對(duì)控制方程的離散則采用二階迎風(fēng)格式,其中能量方程的收斂標(biāo)準(zhǔn)為10-6,其他方程的收斂標(biāo)準(zhǔn)均為10-3。

    1.4.2 基本計(jì)算參數(shù)

    本模擬涉及的物料主要包括煙氣、循環(huán)顆粒以及水分,基本參數(shù)見表1。煙氣和循環(huán)顆粒簡(jiǎn)化如下:

    (1) 入口煙氣為7種氣體和煙塵混合的氣體,由于模擬中不涉及化學(xué)反應(yīng)且忽略煙塵的體積率,因此采用空氣模擬煙氣,并改變相關(guān)參數(shù)。

    表1 模擬基本參數(shù)

    (2) 循環(huán)物料實(shí)際為粉煤灰,由于煙氣降溫的主要因素是水分,因此本模擬中用固體硅顆粒模擬粉煤灰,并改變密度值。

    (3) 為研究煙氣進(jìn)速和顆粒粒徑對(duì)流場(chǎng)的影響,選取最大和最小煙氣進(jìn)速以及不同數(shù)量級(jí)的顆粒粒徑作為設(shè)定模擬條件。

    1.4.3 邊界條件

    煙氣入口和循環(huán)物料入口均設(shè)為速度入口,出口的邊界條件設(shè)為充分發(fā)展,壁面則設(shè)置為無滑移、絕熱的墻。

    2 計(jì)算結(jié)果與討論

    2.1 模擬驗(yàn)證

    在上述設(shè)定條件下模擬得到不同工況下出口溫度、顆粒停留時(shí)間、速度場(chǎng)、濃度場(chǎng)和壓強(qiáng)場(chǎng)等數(shù)據(jù),然而最終模擬結(jié)果與實(shí)測(cè)值有一定誤差,原因如下:(1)模擬時(shí)不考慮煙氣的可壓縮性和反應(yīng)器壁的傳熱性,計(jì)算溫度偏高;(2)工廠運(yùn)行中環(huán)境溫度和壓強(qiáng)等均會(huì)影響實(shí)測(cè)數(shù)據(jù);(3)工廠運(yùn)行中的數(shù)據(jù)測(cè)量點(diǎn)和模擬計(jì)算中的觀測(cè)點(diǎn)不完全一致。

    由李超等[16]和汪建慧[17]的研究可知,煙氣進(jìn)速越大,氣固輸送量增加,顆粒加速度增大,因此脫硫塔整體壓降增大且顆粒停留時(shí)間縮短;顆粒粒徑越大時(shí),顆粒跟隨性越差,氣固間摩擦引起的壓力損失越大,脫硫塔整體壓降增大,顆粒停留時(shí)間延長(zhǎng);顆粒粒徑越大,顆粒比表面積越小,越不利于氣液間傳質(zhì)傳熱,出口溫度升高。圖3至圖5中列出了模擬的顆粒停留時(shí)間、反應(yīng)器壓降和出口溫度隨不同工況的變化,可見這些參數(shù)的變化趨勢(shì)與文獻(xiàn)所述規(guī)律相一致。

    將模擬數(shù)據(jù)與工廠實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析,最高出口溫度模擬值與實(shí)測(cè)值的相對(duì)偏差僅為4.04%,最大壓降和最小壓降的模擬結(jié)果與實(shí)側(cè)數(shù)據(jù)相比存在一定偏差,但相對(duì)偏差最大為12.60%。

    圖3 不同工況下顆粒停留時(shí)間Fig.3 Residence time under different conditions

    圖4 不同工況下反應(yīng)器的壓降Fig.4 Pressure drop across reactor under different conditions

    由蔡容容等[18]的研究可知,NID中顆粒在脫硫反應(yīng)器內(nèi)的停留時(shí)間的理論值為1.0~2.0 s,模擬結(jié)果與理論值相符(見圖(3))。綜合分析,本研究建立的模型和相應(yīng)算法具有較高的準(zhǔn)確性和較好的預(yù)測(cè)能力。

    2.2 反應(yīng)器內(nèi)速度場(chǎng)研究

    實(shí)際工業(yè)運(yùn)行中,反應(yīng)器內(nèi)顆粒物的速度與走向無法測(cè)定,而利用Fluent軟件可以模擬得到不同工況下反應(yīng)器內(nèi)顆粒的速度分布。以顆粒粒徑為900.0 μm,煙氣進(jìn)速為15 m/s的工況為例,反應(yīng)器內(nèi)顆粒的速度云圖模擬結(jié)果見圖6。

    圖5 不同工況下反應(yīng)器出口溫度Fig.5 Outlet temperature under different conditions

    圖6 反應(yīng)器內(nèi)顆粒速度云圖Fig.6 Velocity contour of particles in the desulfurization reactor

    由圖6可見,由于氣固混合的不均性導(dǎo)致反應(yīng)器內(nèi)顆粒產(chǎn)生明顯的速度分層,根據(jù)顆粒物流動(dòng)特征,將反應(yīng)器循環(huán)物料入口側(cè)的低速區(qū)進(jìn)一步劃分為5個(gè)區(qū)域,不同區(qū)域的顆粒速度矢量圖見圖7。

    由圖7可見,低速區(qū)內(nèi)存在明顯的回流區(qū),這是由于進(jìn)入反應(yīng)器內(nèi)的顆粒主要受到兩個(gè)力:自身的重力和煙氣的曳力。在左側(cè)高速區(qū)內(nèi)顆粒受到的曳力大于重力,顆粒隨煙氣向上加速并流運(yùn)動(dòng)。在右側(cè)低速區(qū)內(nèi),氣固接觸不充分,顆粒周圍氣速過低,受到的曳力小于重力,顆粒減速向上運(yùn)動(dòng),最后回落形成回流區(qū)1。隨著反應(yīng)器高度增加,回流區(qū)1外側(cè)的顆粒向內(nèi)側(cè)補(bǔ)充,此時(shí)顆粒的速度大于回流區(qū)1內(nèi)顆粒的速度,但是曳力仍小于重力,形成了回流區(qū)2?;亓鲄^(qū)2外的顆粒繼續(xù)補(bǔ)充時(shí),受到的曳力大于重力,不再產(chǎn)生回流區(qū),顆粒向上形成并流區(qū)。

    表2列出了不同工況下反應(yīng)器內(nèi)回流區(qū)域的大小,可見當(dāng)煙氣進(jìn)速增大時(shí),總體來看低速區(qū)范圍變小,回流區(qū)變窄,高度增加。而顆粒粒徑減小時(shí),跟隨性增強(qiáng),回流區(qū)內(nèi)顆粒量減小,回流區(qū)面積減小。當(dāng)粒徑為0.5 μm時(shí),反應(yīng)器內(nèi)僅存在一個(gè)回流區(qū)。

    2.3 反應(yīng)器內(nèi)顆粒濃度分布研究

    顆粒濃度的分布同樣可以反映輸送過程中顆粒的運(yùn)動(dòng)狀態(tài),圖8列出了煙氣進(jìn)速分別為15、18 m/s時(shí),不同粒徑顆粒在反應(yīng)器水平截面上的顆粒濃度(以體積分?jǐn)?shù)表征)沿高度上的變化。綜合分析可知,在1.000~7.000 m高度內(nèi),由于有回流區(qū)的存在,顆粒濃度基數(shù)較高,又由于在該范圍內(nèi)部分顆粒處于加速狀態(tài),顆粒濃度以較高的速度降低;在7.000 m以上的區(qū)域內(nèi),顆?;亓鳜F(xiàn)象消失,顆粒隨煙氣并流向上運(yùn)動(dòng),顆粒平均加速度降低,顆粒濃度出現(xiàn)降速變緩的趨勢(shì),且煙氣進(jìn)速越小、粒徑越小時(shí),顆粒濃度降速變緩的現(xiàn)象越明顯。結(jié)合表2可知,煙氣進(jìn)速相同時(shí),顆粒粒徑越小,跟隨性越好,形成的總回流區(qū)范圍也就越小,從而顆粒越早進(jìn)入并流向上運(yùn)動(dòng)段,最終顆粒濃度分布基本不隨高度變化,達(dá)到穩(wěn)定狀態(tài)??傮w看來,當(dāng)顆粒粒徑相同時(shí),煙氣進(jìn)速越大,循環(huán)顆粒進(jìn)料量成比例增加,回流區(qū)內(nèi)顆粒濃度增加,顆粒濃度達(dá)到穩(wěn)定的時(shí)間延長(zhǎng)。

    圖7 不同區(qū)域內(nèi)顆粒速度矢量圖Fig.7 Velocity vector of particles in different area of desulfurization reactor

    表2 不同工況下反應(yīng)器內(nèi)回流區(qū)面積

    圖8 不同工況下顆粒濃度隨反應(yīng)器高度的變化Fig.8 Effect of reactor height on particles concentration under different conditions

    圖9 不同工況下反應(yīng)器內(nèi)壓降隨高度的變化Fig.9 Effect of reactor height on pressure drop in the reactor under different conditions

    結(jié)合圖8對(duì)垂直方向上的顆粒濃度分布進(jìn)行分析,可見在0.500~1.000 m進(jìn)料區(qū)循環(huán)物料低速溢流進(jìn)料,局部顆粒濃度相對(duì)偏高;在1.000~7.000 m回流區(qū)內(nèi),進(jìn)入反應(yīng)器的顆粒沿著回流區(qū)外側(cè)向上運(yùn)動(dòng),同時(shí)部分回流區(qū)內(nèi)的顆粒脫離回流區(qū)域,進(jìn)入左側(cè)高速區(qū)域。因此,總體看來,回流區(qū)內(nèi)顆粒濃度隨高度逐漸降低,局部最高濃度出現(xiàn)在回流區(qū)與高速區(qū)的過渡區(qū)域內(nèi)。當(dāng)顆粒進(jìn)行并流區(qū)后(高度7.000~17.312 m),其受到的曳力不小于重力,顆粒隨煙氣并流向上運(yùn)動(dòng)進(jìn)入加速段,平均濃度下降。

    2.4 反應(yīng)器內(nèi)壓強(qiáng)場(chǎng)的研究

    反應(yīng)器內(nèi)整體壓降的大小與氣固輸送量和氣固間摩擦損失有關(guān),且回流現(xiàn)象的存在對(duì)反應(yīng)器的壓強(qiáng)分布有著明顯的影響,圖9顯示了不同工況下反應(yīng)器內(nèi)壓強(qiáng)隨高度的變化。結(jié)合圖4可以看出,當(dāng)煙氣進(jìn)速增大時(shí),脫硫塔氣固輸送量增加,且氣速與顆粒粒徑的增大均會(huì)加大氣固相間的速度差,即增大了煙氣和顆粒間摩擦力,從而使反應(yīng)器整體壓降增大。

    結(jié)合圖8分析可知,顆粒進(jìn)入反應(yīng)器后顆粒濃度高速下降,相對(duì)應(yīng)的反應(yīng)器內(nèi)壓降大幅度增加。結(jié)合速度場(chǎng)分析可知,不同工況的反應(yīng)器內(nèi)均存在著回流現(xiàn)象,且回流區(qū)內(nèi)顆粒和器壁的碰撞摩擦均會(huì)增加壓力損失,因此在該范圍內(nèi)會(huì)出現(xiàn)局部低壓。當(dāng)顆粒粒徑為500.0、900.0 μm時(shí),氣固相間速度差值最大且回流區(qū)域最廣,此時(shí)局部壓強(qiáng)低于出口壓強(qiáng),不利于設(shè)備的安全運(yùn)行。而在小粒徑顆粒的工況下,隨著煙氣進(jìn)速的降低,回流區(qū)范圍變大,塔內(nèi)壓降存在明顯的波動(dòng),不利于設(shè)備的穩(wěn)定運(yùn)行。

    實(shí)際運(yùn)行過程中,反應(yīng)器內(nèi)局部壓降過大不僅對(duì)生產(chǎn)設(shè)備有更高的要求,也會(huì)產(chǎn)生安全隱患,因此,從安全性和穩(wěn)定性考慮,高煙氣進(jìn)速(18 m/s)和小粒徑顆粒(0.5~50.0 μm)最有利于實(shí)際工業(yè)運(yùn)行中的安全生產(chǎn)。

    3 結(jié) 論

    基于顆粒動(dòng)力學(xué)理論,對(duì)NID脫硫反應(yīng)器的流場(chǎng)進(jìn)行三維模擬,得出了不同工況下反應(yīng)器出口溫度、顆粒停留時(shí)間、速度分布、濃度分布和塔內(nèi)壓強(qiáng)分布等數(shù)據(jù)。考察并分析了不同煙氣進(jìn)速及顆粒粒徑工況下,反應(yīng)器內(nèi)的顆粒流動(dòng)細(xì)節(jié)。數(shù)值模擬結(jié)果表明:

    (1) 將顆粒相做擬流體,利用Fluent軟件模擬NID脫硫反應(yīng)器的氣體-顆粒兩相流體模型,發(fā)現(xiàn)Fluent軟件對(duì)稀相氣力輸送具有很好的適應(yīng)性,模擬結(jié)果與文獻(xiàn)及工廠數(shù)據(jù)基本吻合,F(xiàn)luent軟件可以成為NID脫硫反應(yīng)器系統(tǒng)設(shè)計(jì)和優(yōu)化的有力工具。

    (2) 在不同的工況下,反應(yīng)器內(nèi)速度分層現(xiàn)象一直存在,煙氣進(jìn)速越大,分層現(xiàn)象越明顯,且在低速區(qū)域內(nèi)存在不可忽略的回流現(xiàn)象。回流現(xiàn)象的存在增加了顆粒濃度分布的不均性,對(duì)塔內(nèi)壓降分布產(chǎn)生不利影響。

    (3) 顆粒粒徑越小,反應(yīng)器內(nèi)回流區(qū)域越小,顆粒分布均勻性增加,塔內(nèi)局部低壓現(xiàn)象改善;煙氣進(jìn)速越大,塔內(nèi)輸送量增加,脫硫塔整體壓降增大,反應(yīng)器內(nèi)回流區(qū)域變窄,塔內(nèi)壓強(qiáng)變化的穩(wěn)定性增加。

    (4) 從運(yùn)行成本、操作安全性以及塔內(nèi)氣固混合均勻性等角度出發(fā),實(shí)際NID反應(yīng)器運(yùn)行中應(yīng)以18 m/s煙氣進(jìn)速和0.5~50.0 μm顆粒粒徑為宜。

    [1] 楊飏.煙氣脫硫脫硝凈化工程技術(shù)與設(shè)備[M].北京:化學(xué)工業(yè)出版社,2013.

    [2] WANG Yunjun,DUAN Yufeng,YANG Liguo,et al.Experimental research and microscopic mechanism analysis on simultaneous mercury removal by novel integrated semi-dry desulfurization systems[J].Journal of Power Engineering,2009,29(3):270-276.

    [3] WU Chengjun,DUAN Yufeng,WANG Yunjun,et al.Characteristics of mercury emission and demercurization property of NID system of a 410 t/h pulverized coal fired boiler[J].Journal of Fuel Chemistry and Technology,2008,36(5):540-544.

    [4] WANG Yunjun,DUAN Yufeng,YANG Liguo,et al.Mercury speciation and emission from the coal-fired power plant filled with flue gas desulfurization equipment[J].The Canadian Journal of Chemical Engineering,2010,88(5):867-873.

    [5] 蘇清發(fā),蔣亞兵,吳慕正,等.NID半干法脫硫灰生產(chǎn)蒸壓磚中試實(shí)驗(yàn)研究[J].環(huán)境工程學(xué)報(bào),2014,8(9):3996-4000.

    [6] 張善榮.散料輸送與貯存[M].北京:化學(xué)工業(yè)出版社,1994.

    [7] RAUTIAINEN A,STEWART G,POIKOIAINEN V,et al.An experiment study of vertical pneumatic convey[J]. Powder Technology,1999,104(2):139-150.

    [8] LI H,TOMITA Y.Particle velocity and concentration characteristics in a horizontal dilute swirling flow pneumatic conveying[J].Powder Technology,2000,107(1/2):144-152

    [9] 謝灼利,張政.氣力輸送的數(shù)值模擬研究[J].北京化工大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,28(1):22-27.

    [10] 蒲文灝,熊源泉,趙長(zhǎng)遂,等.垂直管煤粉高壓密相氣力輸送特性的模擬研究[J].中國(guó)電機(jī)工程學(xué)報(bào),2008,28(17):21-25.

    [11] 景山,蔣攀亮,胡慶元,等.垂直稀相氣力輸送壓降特征實(shí)驗(yàn)[J].過程工程學(xué)報(bào),2001,1(1):25-29.

    [12] 蒲文灝,趙長(zhǎng)遂,熊源泉,等.水平管加壓密相煤粉氣力輸送數(shù)值模擬[J].化工學(xué)報(bào),2008,59(10):2601-2607.

    [13] MAROCCO L,INZOLI F.Multiphase Euler-Lagrange CFD simulation applied to wet flue gas desulphurisation technology[J].International Journal of Multiphase Flow,2009,35(2):185-194.

    [14] GOMEZ A,FUEYO N,TOMAS A.Detailed modelling of a flue-gas desulfurisation plant[J].Computers & Chemical Engineering,2007,31(11):1419-1431.

    [15] BROWN K,KALATA W,SCHICK R.Optimization of SO2scrubber using CFD modeling[J].Procedia Engineering,2014,83:170-180.

    [16] 李超,代正華,許建良,等.多噴嘴對(duì)置式氣化爐內(nèi)顆粒停留時(shí)間分布數(shù)學(xué)模擬研究[J].高?;瘜W(xué)工程學(xué)報(bào),2011,25(3):416-422.

    [17] 汪建慧.NID循環(huán)半干法煙氣脫硫裝置脫硫效率的分析[J].環(huán)境工程,2011,29(增刊1):150-152.

    [18] 蔡容容,張衍國(guó),蒙晨瑋,等.流化床內(nèi)球形大顆粒停留時(shí)間分布的正交實(shí)驗(yàn)研究[J].中國(guó)電機(jī)工程學(xué)報(bào),2014(5):713-717.

    Numericalsimulationofflowfieldinthenovelintegrateddesulphurizationreactor

    YANGJingxian,ZHANGHaitao.

    (InstituteofChemicalEngineering,EastChinaUniversityofScienceandTechnology,Shanghai200237)

    Gas-particles two fluid model of Fluent software was introduced to study three dimensional flow behaviors of pneumatic conveying in novel integrated desulphrization (NID) reactor on the basis of the kinetic theory of granular flow. The simulation was run with different inlet velocity and particle size. The result showed that the simulated data had good agreement with the literature data and plant actual monitoring data. The simulated data includes the distribution of particle velocity and particle concentration, the formation mechanism of flow field in the reactor and the effect of inlet velocity, particle size on pressure drop of reactor were investigated. The result indicates that as the inlet velocity was 18 m/s and the particle size was 0.5~50.0 μm, the area of back flow was minmized, which optimized the stability and security of NID reactor in industrial operation.

    NID; pneumatic conveying; numerical simulation; pressure drop

    楊竟憲,女,1991年生,碩士研究生,研究方向?yàn)闊煔饷摿颉?

    。

    10.15985/j.cnki.1001-3865.2016.12.016

    編輯:丁 懷 (

    2016-07-05)

    猜你喜歡
    氣力輸送塔內(nèi)流場(chǎng)
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場(chǎng)計(jì)算
    轉(zhuǎn)杯紡排雜區(qū)流場(chǎng)與排雜性能
    氣力輸送系統(tǒng)和管鏈機(jī)輸送系統(tǒng)在石子煤輸送中的比較
    鍋爐脫硫沉積對(duì)節(jié)能環(huán)保的影響
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場(chǎng)統(tǒng)計(jì)分析
    石灰石-石膏濕法煙氣脫硫塔內(nèi)流場(chǎng)模擬及優(yōu)化分析
    高濃度高氣壓在燒結(jié)用石灰氣力輸送中的應(yīng)用
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    炭黑氣力輸送裝置主要元器件的選擇
    一種粉料氣力輸送槽
    精品日产1卡2卡| 午夜福利一区二区在线看| 国产一区二区激情短视频| 亚洲欧美一区二区三区黑人| 女人被狂操c到高潮| 色综合站精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 欧美精品亚洲一区二区| 国产日韩一区二区三区精品不卡| 最新在线观看一区二区三区| 最新在线观看一区二区三区| 操出白浆在线播放| 欧美日本中文国产一区发布| 久久精品91蜜桃| 国产精品久久电影中文字幕| 51午夜福利影视在线观看| 一区福利在线观看| 亚洲成人免费电影在线观看| 久久青草综合色| 少妇的丰满在线观看| 1024视频免费在线观看| 亚洲精品美女久久久久99蜜臀| 国产精品影院久久| 欧美日韩精品网址| 精品久久久久久久久久免费视频| 欧美日韩福利视频一区二区| 禁无遮挡网站| 丁香六月欧美| 1024视频免费在线观看| 99国产综合亚洲精品| 人人妻人人爽人人添夜夜欢视频| 手机成人av网站| 亚洲精品国产色婷婷电影| 桃色一区二区三区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 十八禁人妻一区二区| 老司机深夜福利视频在线观看| 色精品久久人妻99蜜桃| 成在线人永久免费视频| 欧美黄色片欧美黄色片| 波多野结衣高清无吗| 亚洲成人久久性| 一边摸一边抽搐一进一出视频| 久久九九热精品免费| 99国产精品一区二区三区| www.精华液| 少妇裸体淫交视频免费看高清 | 日韩精品中文字幕看吧| 日韩大尺度精品在线看网址 | 看黄色毛片网站| av片东京热男人的天堂| 国产高清有码在线观看视频 | 高清毛片免费观看视频网站| 高清黄色对白视频在线免费看| 亚洲伊人色综图| 色尼玛亚洲综合影院| 亚洲精品国产精品久久久不卡| 久久 成人 亚洲| 日本一区二区免费在线视频| 丝袜人妻中文字幕| 啪啪无遮挡十八禁网站| 国产av一区二区精品久久| 亚洲欧美激情在线| 亚洲av片天天在线观看| 国产精品 欧美亚洲| 午夜老司机福利片| 淫秽高清视频在线观看| 国产私拍福利视频在线观看| 好男人电影高清在线观看| 精品福利观看| 韩国av一区二区三区四区| 黄片播放在线免费| 午夜老司机福利片| 这个男人来自地球电影免费观看| 久久久久久久久久久久大奶| 亚洲九九香蕉| 女性生殖器流出的白浆| 久久中文看片网| 国产不卡一卡二| 国产野战对白在线观看| 久热爱精品视频在线9| 波多野结衣高清无吗| 狠狠狠狠99中文字幕| 国产一区二区激情短视频| 成人亚洲精品一区在线观看| a级毛片在线看网站| 大香蕉久久成人网| 国内精品久久久久久久电影| 桃红色精品国产亚洲av| 亚洲黑人精品在线| 亚洲国产欧美一区二区综合| 人人妻人人爽人人添夜夜欢视频| 日韩欧美国产一区二区入口| 亚洲精品国产区一区二| 亚洲最大成人中文| 老司机在亚洲福利影院| 黑人操中国人逼视频| 亚洲天堂国产精品一区在线| svipshipincom国产片| 一本久久中文字幕| 久久亚洲真实| 国产精华一区二区三区| 女人被狂操c到高潮| 一级毛片精品| 亚洲五月色婷婷综合| ponron亚洲| 男女下面进入的视频免费午夜 | 久久人人爽av亚洲精品天堂| 老司机福利观看| 琪琪午夜伦伦电影理论片6080| 在线永久观看黄色视频| 女人高潮潮喷娇喘18禁视频| 无遮挡黄片免费观看| 老司机深夜福利视频在线观看| 日日夜夜操网爽| 日韩大尺度精品在线看网址 | 国产精品野战在线观看| 老司机午夜十八禁免费视频| 午夜福利在线观看吧| 久久人妻福利社区极品人妻图片| 成人国语在线视频| 十八禁人妻一区二区| 黄频高清免费视频| 国产一区二区在线av高清观看| 两个人免费观看高清视频| 人人妻人人澡欧美一区二区 | 极品教师在线免费播放| 午夜福利免费观看在线| 亚洲av美国av| 国产一区二区三区在线臀色熟女| 搡老妇女老女人老熟妇| 女生性感内裤真人,穿戴方法视频| 欧美日本中文国产一区发布| 亚洲成人精品中文字幕电影| 无遮挡黄片免费观看| 亚洲成人免费电影在线观看| 一级a爱视频在线免费观看| 亚洲三区欧美一区| 久久久久国产精品人妻aⅴ院| 国产亚洲精品一区二区www| 免费在线观看日本一区| 最近最新中文字幕大全电影3 | 日韩欧美免费精品| 叶爱在线成人免费视频播放| 国产成人精品在线电影| www.999成人在线观看| 在线免费观看的www视频| 女人高潮潮喷娇喘18禁视频| 纯流量卡能插随身wifi吗| 亚洲美女黄片视频| 婷婷精品国产亚洲av在线| 搡老岳熟女国产| 一级a爱片免费观看的视频| 一区福利在线观看| 午夜福利欧美成人| av在线播放免费不卡| 色综合亚洲欧美另类图片| 久久久久亚洲av毛片大全| 麻豆久久精品国产亚洲av| 亚洲中文字幕一区二区三区有码在线看 | 淫妇啪啪啪对白视频| 成年人黄色毛片网站| 久9热在线精品视频| 欧美激情高清一区二区三区| 岛国视频午夜一区免费看| 天堂动漫精品| www.999成人在线观看| www.999成人在线观看| 亚洲性夜色夜夜综合| 亚洲五月天丁香| a级毛片在线看网站| 午夜亚洲福利在线播放| 欧美激情极品国产一区二区三区| 精品乱码久久久久久99久播| 日本精品一区二区三区蜜桃| 欧美老熟妇乱子伦牲交| 他把我摸到了高潮在线观看| 搡老岳熟女国产| 欧美黑人精品巨大| 少妇被粗大的猛进出69影院| 免费少妇av软件| 97人妻精品一区二区三区麻豆 | 精品日产1卡2卡| 色尼玛亚洲综合影院| 久久久久精品国产欧美久久久| 多毛熟女@视频| 欧美成人免费av一区二区三区| 法律面前人人平等表现在哪些方面| 高清在线国产一区| 夜夜躁狠狠躁天天躁| 国产精品电影一区二区三区| 99久久99久久久精品蜜桃| 美女国产高潮福利片在线看| 999久久久国产精品视频| 午夜福利18| 亚洲成人国产一区在线观看| 色综合站精品国产| 精品国产乱码久久久久久男人| 国产亚洲av嫩草精品影院| 久久午夜亚洲精品久久| 久99久视频精品免费| 亚洲免费av在线视频| 亚洲精品国产一区二区精华液| 成年人黄色毛片网站| 91麻豆av在线| 人妻丰满熟妇av一区二区三区| 又黄又粗又硬又大视频| 欧美日韩黄片免| 91字幕亚洲| 国产精品久久久久久亚洲av鲁大| 色综合亚洲欧美另类图片| 国产成人精品无人区| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美国产一区二区入口| 亚洲精品久久成人aⅴ小说| 免费少妇av软件| av在线播放免费不卡| 一二三四社区在线视频社区8| 十八禁网站免费在线| 极品教师在线免费播放| 十八禁人妻一区二区| 99久久综合精品五月天人人| 亚洲精品国产色婷婷电影| 欧美一级a爱片免费观看看 | 老司机靠b影院| 波多野结衣av一区二区av| 99精品欧美一区二区三区四区| 99国产综合亚洲精品| 久久九九热精品免费| 国产成人欧美在线观看| 一级毛片女人18水好多| 国产精品国产高清国产av| 欧美日韩精品网址| 最好的美女福利视频网| 亚洲熟妇中文字幕五十中出| 18美女黄网站色大片免费观看| 午夜老司机福利片| 国产成人啪精品午夜网站| 色综合欧美亚洲国产小说| 99精品在免费线老司机午夜| 曰老女人黄片| 91字幕亚洲| АⅤ资源中文在线天堂| 国产成人精品久久二区二区91| 亚洲人成电影观看| 中文字幕人妻熟女乱码| 亚洲国产精品久久男人天堂| netflix在线观看网站| 国产欧美日韩一区二区三区在线| 在线观看免费视频网站a站| 国产av精品麻豆| 久久午夜亚洲精品久久| 亚洲精品国产一区二区精华液| 搡老熟女国产l中国老女人| 麻豆久久精品国产亚洲av| 一进一出好大好爽视频| 99riav亚洲国产免费| 操美女的视频在线观看| 色尼玛亚洲综合影院| 无限看片的www在线观看| 国产男靠女视频免费网站| 亚洲av美国av| av视频在线观看入口| 很黄的视频免费| 日日爽夜夜爽网站| 亚洲国产高清在线一区二区三 | 久久久久久免费高清国产稀缺| 久热这里只有精品99| cao死你这个sao货| 久久久久久亚洲精品国产蜜桃av| 久久国产亚洲av麻豆专区| 久久久久久人人人人人| 国产午夜精品久久久久久| 精品国产一区二区久久| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播放欧美日韩| 日韩精品免费视频一区二区三区| 91成年电影在线观看| 精品欧美国产一区二区三| 99精品在免费线老司机午夜| 精品国产超薄肉色丝袜足j| 变态另类丝袜制服| 给我免费播放毛片高清在线观看| 中出人妻视频一区二区| 亚洲精品av麻豆狂野| 国产一级毛片七仙女欲春2 | 日韩精品青青久久久久久| 中文字幕人妻丝袜一区二区| 欧美成人一区二区免费高清观看 | 美女扒开内裤让男人捅视频| 亚洲一区中文字幕在线| 一a级毛片在线观看| 国产精品 国内视频| av天堂久久9| 亚洲中文字幕日韩| 91精品三级在线观看| 一本久久中文字幕| 亚洲第一电影网av| 一边摸一边抽搐一进一出视频| 亚洲国产看品久久| 夜夜爽天天搞| 日本三级黄在线观看| 在线国产一区二区在线| 一二三四在线观看免费中文在| 亚洲成人精品中文字幕电影| 欧美日本中文国产一区发布| a级毛片在线看网站| 大码成人一级视频| 久久国产亚洲av麻豆专区| 亚洲三区欧美一区| 日韩 欧美 亚洲 中文字幕| 久久精品人人爽人人爽视色| 两性夫妻黄色片| 午夜福利成人在线免费观看| 手机成人av网站| 日本五十路高清| 免费高清视频大片| 久久精品91无色码中文字幕| 女人被狂操c到高潮| 91精品国产国语对白视频| 黄网站色视频无遮挡免费观看| 国产精品香港三级国产av潘金莲| 男女之事视频高清在线观看| 国产国语露脸激情在线看| 国产日韩一区二区三区精品不卡| 不卡一级毛片| 高清在线国产一区| 好看av亚洲va欧美ⅴa在| 麻豆久久精品国产亚洲av| 国产欧美日韩一区二区三| 亚洲九九香蕉| 国产麻豆69| 免费在线观看视频国产中文字幕亚洲| 12—13女人毛片做爰片一| 宅男免费午夜| 欧美成人性av电影在线观看| e午夜精品久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲美女黄片视频| 91成年电影在线观看| 久9热在线精品视频| 久久人妻av系列| 伊人久久大香线蕉亚洲五| 中国美女看黄片| 50天的宝宝边吃奶边哭怎么回事| av中文乱码字幕在线| 一级,二级,三级黄色视频| 亚洲男人的天堂狠狠| 欧美成人午夜精品| 18禁国产床啪视频网站| 高清毛片免费观看视频网站| 国产片内射在线| 操出白浆在线播放| 婷婷丁香在线五月| 熟妇人妻久久中文字幕3abv| 久久久久亚洲av毛片大全| 国产一区二区激情短视频| 真人一进一出gif抽搐免费| 欧美一级a爱片免费观看看 | 亚洲一区高清亚洲精品| 色婷婷久久久亚洲欧美| 国产野战对白在线观看| 久久中文字幕一级| 99久久国产精品久久久| 我的亚洲天堂| 老司机靠b影院| 香蕉丝袜av| av天堂久久9| av中文乱码字幕在线| 亚洲国产精品成人综合色| 精品第一国产精品| 国产成人系列免费观看| 一区在线观看完整版| 欧美国产精品va在线观看不卡| 国产私拍福利视频在线观看| 黄色片一级片一级黄色片| 日韩国内少妇激情av| 在线永久观看黄色视频| 母亲3免费完整高清在线观看| 人妻久久中文字幕网| 999久久久精品免费观看国产| 成人免费观看视频高清| 在线av久久热| 亚洲欧美日韩另类电影网站| 国产成年人精品一区二区| 99久久国产精品久久久| 日韩精品中文字幕看吧| 免费在线观看完整版高清| 波多野结衣一区麻豆| 国产熟女xx| 亚洲男人天堂网一区| 国产麻豆成人av免费视频| 中文字幕高清在线视频| 国产三级在线视频| 久久久久久久久久久久大奶| www.熟女人妻精品国产| 激情视频va一区二区三区| 在线天堂中文资源库| 精品国产一区二区久久| 9191精品国产免费久久| 精品久久久久久,| 亚洲一码二码三码区别大吗| 亚洲第一电影网av| 黄色 视频免费看| 午夜福利免费观看在线| 欧美中文日本在线观看视频| 欧美日韩瑟瑟在线播放| 久久精品国产亚洲av香蕉五月| 这个男人来自地球电影免费观看| 丁香六月欧美| 亚洲精品一区av在线观看| 国产精品av久久久久免费| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av五月六月丁香网| 国产精品免费一区二区三区在线| 后天国语完整版免费观看| 亚洲专区字幕在线| 巨乳人妻的诱惑在线观看| 欧美乱妇无乱码| 搞女人的毛片| 99热只有精品国产| 久99久视频精品免费| 国产精品乱码一区二三区的特点 | 久久午夜亚洲精品久久| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| 好男人在线观看高清免费视频 | 热re99久久国产66热| 亚洲视频免费观看视频| 最新在线观看一区二区三区| 天天一区二区日本电影三级 | 亚洲精品国产一区二区精华液| 老司机靠b影院| 日韩欧美免费精品| 91大片在线观看| 亚洲欧美日韩无卡精品| 视频在线观看一区二区三区| 可以在线观看毛片的网站| 国产精品久久久av美女十八| 日韩大码丰满熟妇| 丝袜美足系列| 国产精品二区激情视频| 一级毛片高清免费大全| 在线免费观看的www视频| 久久久久九九精品影院| 黑人操中国人逼视频| 一区二区日韩欧美中文字幕| 久久香蕉激情| 国产精品亚洲美女久久久| 在线观看免费视频网站a站| 久久九九热精品免费| 欧美日本视频| 18禁黄网站禁片午夜丰满| 欧美中文综合在线视频| 国产午夜福利久久久久久| 午夜亚洲福利在线播放| 午夜免费鲁丝| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡免费网站照片 | 国产高清有码在线观看视频 | 熟妇人妻久久中文字幕3abv| 久久亚洲精品不卡| 69精品国产乱码久久久| 757午夜福利合集在线观看| 十八禁网站免费在线| 色av中文字幕| 男人舔女人的私密视频| 老鸭窝网址在线观看| 国产亚洲精品av在线| 亚洲视频免费观看视频| 一进一出抽搐动态| 曰老女人黄片| 成人免费观看视频高清| 精品日产1卡2卡| 国产高清视频在线播放一区| 两性夫妻黄色片| 女性被躁到高潮视频| 久久久国产精品麻豆| 波多野结衣巨乳人妻| 国产在线观看jvid| 1024视频免费在线观看| 国产亚洲精品综合一区在线观看 | x7x7x7水蜜桃| 少妇 在线观看| 国产成人精品在线电影| 亚洲精品国产一区二区精华液| 国产精品综合久久久久久久免费 | 亚洲免费av在线视频| 99久久久亚洲精品蜜臀av| 国产成人影院久久av| 女性被躁到高潮视频| 国产亚洲欧美在线一区二区| 亚洲欧洲精品一区二区精品久久久| 搡老岳熟女国产| 美女 人体艺术 gogo| 午夜视频精品福利| 国产亚洲av嫩草精品影院| 国语自产精品视频在线第100页| 99国产综合亚洲精品| 免费一级毛片在线播放高清视频 | 欧美成人免费av一区二区三区| 人人妻,人人澡人人爽秒播| 日本一区二区免费在线视频| 久久国产精品影院| 国产精品av久久久久免费| 中文字幕av电影在线播放| 99热只有精品国产| 欧美国产日韩亚洲一区| 久久久国产欧美日韩av| 久9热在线精品视频| 一二三四在线观看免费中文在| 黄频高清免费视频| 黄色毛片三级朝国网站| 免费av毛片视频| 国产亚洲欧美98| ponron亚洲| 99在线视频只有这里精品首页| 美女国产高潮福利片在线看| 国产av精品麻豆| 波多野结衣巨乳人妻| 欧美日本中文国产一区发布| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品999在线| 69精品国产乱码久久久| 变态另类成人亚洲欧美熟女 | 夜夜夜夜夜久久久久| 午夜精品国产一区二区电影| 中文字幕人妻熟女乱码| 久久精品91蜜桃| 中国美女看黄片| 国产精品国产高清国产av| 天天添夜夜摸| 丁香六月欧美| 99国产精品免费福利视频| 久久热在线av| 亚洲第一青青草原| 黄色丝袜av网址大全| 久久国产乱子伦精品免费另类| 精品熟女少妇八av免费久了| 亚洲精品在线美女| 午夜福利在线观看吧| 日韩国内少妇激情av| 亚洲av电影在线进入| 国内精品久久久久精免费| 午夜免费观看网址| 日韩欧美在线二视频| 美女免费视频网站| 操美女的视频在线观看| 午夜精品国产一区二区电影| 亚洲人成电影观看| 国产精品99久久99久久久不卡| 黄频高清免费视频| 91在线观看av| 亚洲无线在线观看| 欧美最黄视频在线播放免费| 又黄又爽又免费观看的视频| 中文字幕人妻丝袜一区二区| 国产高清videossex| 国产成人系列免费观看| 成人国语在线视频| 91成人精品电影| 精品久久蜜臀av无| 一区福利在线观看| 亚洲男人的天堂狠狠| 国产精品香港三级国产av潘金莲| 成人三级做爰电影| 男女下面进入的视频免费午夜 | 亚洲精品美女久久av网站| 免费观看精品视频网站| 啦啦啦 在线观看视频| 看黄色毛片网站| 国产精品av久久久久免费| 欧美乱码精品一区二区三区| 国产精品自产拍在线观看55亚洲| 亚洲一区二区三区不卡视频| 满18在线观看网站| 色老头精品视频在线观看| 中亚洲国语对白在线视频| 亚洲五月色婷婷综合| 一区二区三区精品91| 999久久久精品免费观看国产| 男人操女人黄网站| 一进一出抽搐gif免费好疼| 国产成人av教育| 人人妻,人人澡人人爽秒播| 在线av久久热| 在线观看一区二区三区| 多毛熟女@视频| 在线永久观看黄色视频| 久久影院123| 19禁男女啪啪无遮挡网站| 日韩大尺度精品在线看网址 | 欧美另类亚洲清纯唯美| 久久精品国产亚洲av香蕉五月| 成人18禁高潮啪啪吃奶动态图| 欧美日本中文国产一区发布| 男女午夜视频在线观看| 亚洲欧美精品综合久久99| 在线观看免费视频网站a站| 久久热在线av| 淫秽高清视频在线观看| 亚洲av五月六月丁香网| 久久久国产成人免费| 午夜免费成人在线视频| 可以免费在线观看a视频的电影网站| 亚洲精品久久成人aⅴ小说| 日韩欧美在线二视频| 欧美老熟妇乱子伦牲交| 国产在线精品亚洲第一网站| 91麻豆精品激情在线观看国产| 欧美久久黑人一区二区| 国产午夜精品久久久久久| 香蕉久久夜色| 欧美在线一区亚洲| 少妇 在线观看|