• <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)軸承靜平衡位置求解
    炭黑氣力輸送裝置主要元器件的選擇
    一種粉料氣力輸送槽
    国产淫语在线视频| 日韩精品有码人妻一区| 中文字幕亚洲精品专区| 精品久久久久久久久av| 久久久久久九九精品二区国产| 国产成人午夜福利电影在线观看| 水蜜桃什么品种好| 国产视频首页在线观看| 日韩伦理黄色片| 国产精品人妻久久久久久| 久久精品国产亚洲网站| 一区二区三区乱码不卡18| 亚洲精品色激情综合| 成年版毛片免费区| 99热全是精品| 少妇裸体淫交视频免费看高清| 国产不卡一卡二| 夜夜看夜夜爽夜夜摸| 天堂网av新在线| 一级黄片播放器| 女人十人毛片免费观看3o分钟| 日韩制服骚丝袜av| 亚洲av中文av极速乱| 午夜视频国产福利| 成人av在线播放网站| 久久综合国产亚洲精品| 久久国产乱子免费精品| 欧美最新免费一区二区三区| 免费观看精品视频网站| 亚洲一区高清亚洲精品| 亚洲aⅴ乱码一区二区在线播放| 久久这里只有精品中国| 国产有黄有色有爽视频| 国产伦精品一区二区三区四那| 麻豆乱淫一区二区| 一二三四中文在线观看免费高清| 国产亚洲91精品色在线| 亚洲精品视频女| 国产毛片a区久久久久| 小蜜桃在线观看免费完整版高清| 免费黄网站久久成人精品| 精品久久久久久久末码| 色综合站精品国产| 又粗又硬又长又爽又黄的视频| 欧美xxⅹ黑人| 女人久久www免费人成看片| 午夜精品一区二区三区免费看| 麻豆av噜噜一区二区三区| 亚洲精品日本国产第一区| 一个人免费在线观看电影| 九色成人免费人妻av| 男人舔奶头视频| 精品国产露脸久久av麻豆 | 日本三级黄在线观看| 久久午夜福利片| 免费大片18禁| 国产黄a三级三级三级人| 97人妻精品一区二区三区麻豆| 国产成人午夜福利电影在线观看| 中文字幕免费在线视频6| 久久久色成人| 一夜夜www| 亚洲精品日韩av片在线观看| 网址你懂的国产日韩在线| 国产激情偷乱视频一区二区| 特大巨黑吊av在线直播| 毛片女人毛片| 久久精品久久精品一区二区三区| 久久久久久久久久黄片| 91aial.com中文字幕在线观看| 国产精品久久久久久久电影| 国产精品不卡视频一区二区| 尤物成人国产欧美一区二区三区| 少妇人妻精品综合一区二区| 性插视频无遮挡在线免费观看| 久久久久久久午夜电影| 免费av毛片视频| 街头女战士在线观看网站| 永久网站在线| 国产探花极品一区二区| 欧美潮喷喷水| 极品少妇高潮喷水抽搐| 三级国产精品片| 国产亚洲av嫩草精品影院| 69人妻影院| 综合色av麻豆| 在线观看一区二区三区| 成人毛片a级毛片在线播放| 国内精品美女久久久久久| 国产白丝娇喘喷水9色精品| 成人一区二区视频在线观看| 国产 一区 欧美 日韩| 久久久久网色| 日日摸夜夜添夜夜爱| 91久久精品国产一区二区三区| 精华霜和精华液先用哪个| 男女边摸边吃奶| 亚洲第一区二区三区不卡| 久久鲁丝午夜福利片| 波多野结衣巨乳人妻| videos熟女内射| 国产v大片淫在线免费观看| 一级黄片播放器| 精品久久久精品久久久| 97热精品久久久久久| 欧美成人a在线观看| 观看美女的网站| 日日撸夜夜添| 七月丁香在线播放| 日韩三级伦理在线观看| 久久精品人妻少妇| 国产一级毛片七仙女欲春2| 成人高潮视频无遮挡免费网站| av在线老鸭窝| 精品久久久久久成人av| 18禁动态无遮挡网站| 免费黄色在线免费观看| 大香蕉97超碰在线| 大片免费播放器 马上看| 国产综合精华液| 91在线精品国自产拍蜜月| 精品久久久久久久久av| 日韩av免费高清视频| 国产精品久久久久久精品电影| 久久精品久久久久久噜噜老黄| 哪个播放器可以免费观看大片| 亚洲欧洲国产日韩| 国产乱人视频| 肉色欧美久久久久久久蜜桃 | 欧美日韩亚洲高清精品| 99热这里只有是精品在线观看| 亚洲四区av| 97热精品久久久久久| 色播亚洲综合网| 成人午夜高清在线视频| 亚洲av福利一区| 亚洲精品乱久久久久久| 亚洲欧美成人综合另类久久久| 中文欧美无线码| 免费高清在线观看视频在线观看| 国内少妇人妻偷人精品xxx网站| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩东京热| 国产成人午夜福利电影在线观看| 免费人成在线观看视频色| 免费少妇av软件| 亚洲人成网站高清观看| 国产日韩欧美在线精品| 国产一级毛片在线| 亚洲成人一二三区av| 成年免费大片在线观看| 99热这里只有是精品在线观看| 成年免费大片在线观看| 欧美潮喷喷水| 尤物成人国产欧美一区二区三区| 国产v大片淫在线免费观看| 国产成人aa在线观看| 人妻少妇偷人精品九色| 一本一本综合久久| 丰满少妇做爰视频| 美女cb高潮喷水在线观看| 一个人免费在线观看电影| 青青草视频在线视频观看| 一个人观看的视频www高清免费观看| 国内精品宾馆在线| 日日干狠狠操夜夜爽| 色播亚洲综合网| 久久国内精品自在自线图片| 久久久久九九精品影院| 国产淫片久久久久久久久| 简卡轻食公司| 日韩av在线大香蕉| 国产亚洲av嫩草精品影院| 国产老妇伦熟女老妇高清| 能在线免费看毛片的网站| 日本欧美国产在线视频| 天天躁日日操中文字幕| 一区二区三区乱码不卡18| 国产精品美女特级片免费视频播放器| 久久97久久精品| 久久韩国三级中文字幕| 两个人视频免费观看高清| 久久午夜福利片| 亚洲精品成人av观看孕妇| 国产精品蜜桃在线观看| 老司机影院成人| 2021少妇久久久久久久久久久| 插阴视频在线观看视频| av国产久精品久网站免费入址| 18禁在线无遮挡免费观看视频| 18禁在线无遮挡免费观看视频| videossex国产| 久久国内精品自在自线图片| 男人和女人高潮做爰伦理| 久久久久精品性色| 国产麻豆成人av免费视频| 国产麻豆成人av免费视频| 日韩一区二区三区影片| 国产片特级美女逼逼视频| 一区二区三区乱码不卡18| 最近的中文字幕免费完整| 色尼玛亚洲综合影院| 久久久a久久爽久久v久久| 青春草视频在线免费观看| 国产精品人妻久久久久久| 精品久久久久久久人妻蜜臀av| 久久久精品免费免费高清| 久久久欧美国产精品| 少妇高潮的动态图| 人体艺术视频欧美日本| 在线天堂最新版资源| 麻豆久久精品国产亚洲av| 亚洲精品aⅴ在线观看| freevideosex欧美| eeuss影院久久| 色综合亚洲欧美另类图片| 亚洲精品色激情综合| 久久久精品94久久精品| 青青草视频在线视频观看| 91精品伊人久久大香线蕉| 精品久久久精品久久久| 久久久久精品性色| 欧美丝袜亚洲另类| 韩国高清视频一区二区三区| 欧美日韩亚洲高清精品| 国产大屁股一区二区在线视频| 成年女人在线观看亚洲视频 | 亚洲av电影不卡..在线观看| 日本三级黄在线观看| 国产一区有黄有色的免费视频 | 精品国产露脸久久av麻豆 | 国产精品av视频在线免费观看| 男女边吃奶边做爰视频| 国产一区二区亚洲精品在线观看| 国产黄色小视频在线观看| 精品久久久精品久久久| 国产精品熟女久久久久浪| 欧美3d第一页| 国产精品麻豆人妻色哟哟久久 | 国产精品一区二区三区四区久久| 有码 亚洲区| 激情 狠狠 欧美| 亚洲精品影视一区二区三区av| 在线 av 中文字幕| 超碰97精品在线观看| 国产乱人偷精品视频| videos熟女内射| 国产免费又黄又爽又色| 国产毛片a区久久久久| 男人舔奶头视频| 国产精品一区二区三区四区久久| 日韩大片免费观看网站| 欧美精品国产亚洲| 18禁在线播放成人免费| 两个人视频免费观看高清| 亚洲国产精品成人久久小说| 久久久久久久久久久丰满| h日本视频在线播放| 国产有黄有色有爽视频| 亚洲真实伦在线观看| 亚洲美女视频黄频| 色哟哟·www| 毛片一级片免费看久久久久| 能在线免费观看的黄片| 熟妇人妻不卡中文字幕| 久久这里有精品视频免费| 色综合色国产| 中文天堂在线官网| 黄片wwwwww| 色哟哟·www| 卡戴珊不雅视频在线播放| 久久久久精品性色| 国产91av在线免费观看| 久久久精品免费免费高清| 三级经典国产精品| 我的老师免费观看完整版| 黄片wwwwww| 成人特级av手机在线观看| 99久久精品国产国产毛片| 午夜免费激情av| 国产v大片淫在线免费观看| 国产探花在线观看一区二区| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 久久久精品94久久精品| 一个人观看的视频www高清免费观看| 国产美女午夜福利| 国产精品国产三级国产专区5o| 精品久久久久久久久久久久久| 国产伦理片在线播放av一区| 夜夜看夜夜爽夜夜摸| 91av网一区二区| 麻豆国产97在线/欧美| 欧美xxxx黑人xx丫x性爽| 99热全是精品| 亚洲av成人精品一区久久| 国内精品一区二区在线观看| 欧美日韩综合久久久久久| 国产午夜精品久久久久久一区二区三区| 亚洲av免费在线观看| 欧美zozozo另类| 少妇熟女欧美另类| 丝袜喷水一区| 国产成人精品福利久久| 日韩av不卡免费在线播放| 日韩大片免费观看网站| 两个人视频免费观看高清| 最近2019中文字幕mv第一页| 久久久久久久久久久丰满| 国产毛片a区久久久久| 午夜福利在线观看免费完整高清在| 午夜福利在线在线| 97在线视频观看| 亚洲国产精品成人久久小说| videos熟女内射| 大香蕉久久网| 亚洲欧美日韩无卡精品| 我要看日韩黄色一级片| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 久久国产乱子免费精品| 欧美另类一区| 好男人视频免费观看在线| 青春草国产在线视频| 成人午夜精彩视频在线观看| 大香蕉久久网| 免费看a级黄色片| 久久精品综合一区二区三区| 亚洲精品影视一区二区三区av| 国产一级毛片在线| 麻豆成人av视频| 日韩欧美 国产精品| 少妇的逼好多水| 亚洲av免费高清在线观看| 亚洲人与动物交配视频| 真实男女啪啪啪动态图| 亚洲av二区三区四区| 国产精品国产三级专区第一集| 亚洲av.av天堂| av国产久精品久网站免费入址| 欧美日韩亚洲高清精品| 中文字幕av在线有码专区| 欧美最新免费一区二区三区| 国产精品国产三级专区第一集| 九色成人免费人妻av| 婷婷色综合www| 欧美成人午夜免费资源| 亚洲最大成人中文| 街头女战士在线观看网站| 日韩强制内射视频| 成人av在线播放网站| 国产精品麻豆人妻色哟哟久久 | 中文精品一卡2卡3卡4更新| 一本久久精品| 日日啪夜夜爽| 波野结衣二区三区在线| 国产一区二区在线观看日韩| 天美传媒精品一区二区| 色视频www国产| 中文字幕av在线有码专区| 国内少妇人妻偷人精品xxx网站| 亚洲国产欧美在线一区| 最后的刺客免费高清国语| 丝瓜视频免费看黄片| 免费大片黄手机在线观看| 91av网一区二区| 久久99热这里只频精品6学生| 国产乱来视频区| 亚洲国产日韩欧美精品在线观看| 亚洲精品乱码久久久久久按摩| 夫妻午夜视频| 日韩欧美精品免费久久| 欧美性猛交╳xxx乱大交人| 最近2019中文字幕mv第一页| 免费av观看视频| 国产av在哪里看| 国产亚洲最大av| 大香蕉97超碰在线| 亚洲不卡免费看| 亚洲成色77777| 精品久久久久久电影网| 一区二区三区免费毛片| 精品人妻一区二区三区麻豆| 一个人看的www免费观看视频| 国内揄拍国产精品人妻在线| 丝袜美腿在线中文| 亚洲美女搞黄在线观看| 女人久久www免费人成看片| 又大又黄又爽视频免费| 好男人视频免费观看在线| av线在线观看网站| 亚洲性久久影院| 美女大奶头视频| 卡戴珊不雅视频在线播放| 在线观看人妻少妇| 亚洲电影在线观看av| or卡值多少钱| 国产精品伦人一区二区| 亚洲成人av在线免费| 亚洲怡红院男人天堂| 好男人在线观看高清免费视频| 夜夜爽夜夜爽视频| 亚洲av国产av综合av卡| 一级毛片久久久久久久久女| 欧美日韩亚洲高清精品| 国产精品.久久久| 国产不卡一卡二| 欧美日韩国产mv在线观看视频 | 日本色播在线视频| 国产精品久久久久久av不卡| 中文字幕免费在线视频6| 夜夜看夜夜爽夜夜摸| 91久久精品国产一区二区三区| 亚洲四区av| 人妻一区二区av| 色综合亚洲欧美另类图片| 精品欧美国产一区二区三| 欧美97在线视频| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 日日啪夜夜爽| 亚洲精品自拍成人| 一个人免费在线观看电影| 国产精品蜜桃在线观看| 欧美不卡视频在线免费观看| 久久久久久久久大av| 免费观看的影片在线观看| 日韩一区二区视频免费看| 午夜日本视频在线| 少妇裸体淫交视频免费看高清| 精品99又大又爽又粗少妇毛片| 黄色日韩在线| 国产精品一区二区三区四区久久| 人妻少妇偷人精品九色| 国产伦一二天堂av在线观看| 又大又黄又爽视频免费| 成人无遮挡网站| 中文字幕久久专区| 久久99精品国语久久久| 日韩av在线大香蕉| 久久久欧美国产精品| 亚洲国产av新网站| 亚洲成人一二三区av| 看十八女毛片水多多多| 亚洲精品乱久久久久久| 三级毛片av免费| 成人高潮视频无遮挡免费网站| 久久久久久伊人网av| 国产美女午夜福利| 伦精品一区二区三区| 国产成人91sexporn| 久99久视频精品免费| 一个人看视频在线观看www免费| 日韩制服骚丝袜av| 日韩伦理黄色片| 国产成人a∨麻豆精品| 亚洲三级黄色毛片| 26uuu在线亚洲综合色| 麻豆成人午夜福利视频| 国产精品.久久久| 色哟哟·www| 久久久久久久久久久丰满| 国产乱人视频| 成人无遮挡网站| 高清视频免费观看一区二区 | 成年免费大片在线观看| 日本-黄色视频高清免费观看| 听说在线观看完整版免费高清| 日日啪夜夜爽| 汤姆久久久久久久影院中文字幕 | 99热全是精品| 欧美最新免费一区二区三区| 一级a做视频免费观看| 亚洲国产成人一精品久久久| 免费观看av网站的网址| 国产精品精品国产色婷婷| 欧美高清性xxxxhd video| 亚洲精品aⅴ在线观看| 国产高清有码在线观看视频| 国产日韩欧美在线精品| 亚洲欧美日韩无卡精品| 日本一二三区视频观看| 免费av观看视频| 亚洲av.av天堂| 成年免费大片在线观看| 一级毛片黄色毛片免费观看视频| 激情 狠狠 欧美| videos熟女内射| 汤姆久久久久久久影院中文字幕 | 嘟嘟电影网在线观看| 寂寞人妻少妇视频99o| 观看美女的网站| 天堂网av新在线| 午夜激情久久久久久久| 欧美xxxx性猛交bbbb| 波野结衣二区三区在线| 国产亚洲最大av| 国产亚洲av片在线观看秒播厂 | 日韩av免费高清视频| 亚洲国产成人一精品久久久| 男女边摸边吃奶| 高清欧美精品videossex| 内射极品少妇av片p| 亚洲内射少妇av| 观看免费一级毛片| 婷婷色综合大香蕉| 青春草国产在线视频| 我的老师免费观看完整版| 久久午夜福利片| 国产伦在线观看视频一区| 国产一区二区亚洲精品在线观看| 乱码一卡2卡4卡精品| 美女主播在线视频| 久久精品国产亚洲网站| 亚洲av日韩在线播放| 欧美xxxx黑人xx丫x性爽| 国产午夜精品论理片| 九色成人免费人妻av| 乱码一卡2卡4卡精品| 国产一区二区在线观看日韩| 建设人人有责人人尽责人人享有的 | 成人av在线播放网站| 成人高潮视频无遮挡免费网站| a级毛片免费高清观看在线播放| 国产大屁股一区二区在线视频| 久久精品国产自在天天线| 97精品久久久久久久久久精品| 亚洲欧美一区二区三区黑人 | 成年版毛片免费区| 亚洲av电影在线观看一区二区三区 | 两个人视频免费观看高清| 97超碰精品成人国产| 亚洲av福利一区| 午夜久久久久精精品| 男女下面进入的视频免费午夜| 非洲黑人性xxxx精品又粗又长| 肉色欧美久久久久久久蜜桃 | 乱人视频在线观看| 中文字幕久久专区| 亚洲av.av天堂| 国产91av在线免费观看| 久久99蜜桃精品久久| 人妻少妇偷人精品九色| 熟女电影av网| 97热精品久久久久久| 伦理电影大哥的女人| 舔av片在线| 午夜视频国产福利| 又粗又硬又长又爽又黄的视频| 老师上课跳d突然被开到最大视频| 又大又黄又爽视频免费| 久热久热在线精品观看| 亚洲欧美精品专区久久| 99热6这里只有精品| 日韩一本色道免费dvd| 老司机影院成人| av.在线天堂| 亚洲va在线va天堂va国产| 国产精品一区二区三区四区免费观看| 国产视频首页在线观看| 成年女人在线观看亚洲视频 | 深夜a级毛片| 久久久久性生活片| 国产乱人视频| 亚洲欧美中文字幕日韩二区| 亚洲电影在线观看av| 成人无遮挡网站| 亚洲经典国产精华液单| 女人十人毛片免费观看3o分钟| 99re6热这里在线精品视频| 亚洲av电影在线观看一区二区三区 | 精品久久久久久久末码| 久久久久久久久久人人人人人人| 亚洲欧美清纯卡通| 三级男女做爰猛烈吃奶摸视频| 一区二区三区免费毛片| av福利片在线观看| 一个人观看的视频www高清免费观看| 欧美不卡视频在线免费观看| 成人午夜高清在线视频| 国产成人freesex在线| 日本av手机在线免费观看| 中文资源天堂在线| 18禁在线播放成人免费| 亚洲性久久影院| 免费看美女性在线毛片视频| 久久久久久久久久成人| 18禁在线播放成人免费| 国产午夜精品论理片| 青青草视频在线视频观看| 亚洲性久久影院| 一夜夜www| 亚洲欧美日韩东京热| 国产黄片视频在线免费观看| 免费无遮挡裸体视频| 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 美女内射精品一级片tv| 国产高清不卡午夜福利| 国模一区二区三区四区视频| 九九在线视频观看精品| 日本一二三区视频观看| 午夜福利在线观看免费完整高清在| 2021少妇久久久久久久久久久| 亚洲国产高清在线一区二区三| 亚洲最大成人手机在线| 午夜福利网站1000一区二区三区| 国产淫片久久久久久久久| ponron亚洲| 亚洲人与动物交配视频| 免费高清在线观看视频在线观看| 麻豆精品久久久久久蜜桃| 免费高清在线观看视频在线观看| 国产精品国产三级国产专区5o| 亚洲真实伦在线观看|