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

    搗固焦?fàn)t高壓氨水消煙的數(shù)值模擬

    2015-12-07 02:35:28
    關(guān)鍵詞:數(shù)值模擬

    ?

    搗固焦?fàn)t高壓氨水消煙的數(shù)值模擬

    陳光1a,陳賽1a,鄭光明2,董璐1b,閆丹1a,王雙1a

    (1.安徽工業(yè)大學(xué)a.能源與環(huán)境學(xué)院,b.建筑工程學(xué)院,安徽馬鞍山243002; 2.安徽同興環(huán)保工程技術(shù)有限責(zé)任公司,安徽馬鞍山238100)

    摘要:搗固焦?fàn)t采用高壓氨水噴射消煙技術(shù)控制裝煤產(chǎn)生的煙塵。為研究裝煤時(shí)氨水噴射速度與炭化室產(chǎn)氣量之間的關(guān)系,基于工藝過程和Fluent軟件將某鋼廠4.3 m搗固焦?fàn)t引射排氣系統(tǒng)簡(jiǎn)化成二維模型,采用Mixture模型對(duì)氨水噴射過程進(jìn)行模擬。結(jié)果表明:氨水噴射速度保持不變時(shí),隨著炭化室產(chǎn)氣量增加,炭化室機(jī)側(cè)爐門的氣體由向室內(nèi)流動(dòng)逐漸轉(zhuǎn)變?yōu)橄蚴彝饬鲃?dòng);裝煤時(shí)集氣管煙氣量為0.6 m3/s時(shí),為保證機(jī)側(cè)爐門氣體靜止,所需氨水噴射速度為28.4 m/s;炭化室機(jī)側(cè)爐門氣體保持靜止時(shí),隨著炭化室產(chǎn)氣量增加,氨水噴射速度隨之幾乎呈線性增加。

    關(guān)鍵詞:搗固焦?fàn)t;高壓氨水;Mixture模型;數(shù)值模擬

    搗固焦?fàn)t因其能夠增加高揮發(fā)分和弱黏結(jié)煤用量、提高焦炭質(zhì)量等優(yōu)勢(shì)已成為煉焦技術(shù)的重點(diǎn)發(fā)展方向之一。裝煤煙塵治理成為制約搗固焦?fàn)t工藝進(jìn)一步發(fā)展的重要因素[1-2]。搗固焦?fàn)t采用側(cè)裝煤,裝煤時(shí)由于機(jī)側(cè)爐口長(zhǎng)時(shí)間處于敞開狀態(tài),炭化室內(nèi)產(chǎn)生的煙塵量大,壓力較高,易造成煙塵和荒煤氣從炭化室機(jī)側(cè)爐門處逸出,造成環(huán)境污染。若消煙處理效果不佳會(huì)對(duì)大氣環(huán)境及生產(chǎn)工作區(qū)域造成重度空氣污染[3-6]。

    近年來,高壓氨水噴射消煙裝煤技術(shù)成為焦?fàn)t煙塵治理的發(fā)展方向。高壓氨水噴射消煙原理是在橋管部位噴射高壓氨水,氨水速度較大,以使上升管和炭化室形成較大吸力,可使上升管根部產(chǎn)生約200 Pa的負(fù)壓,在機(jī)側(cè)爐門處的壓力處于零壓或者微負(fù)壓,將裝煤時(shí)產(chǎn)生的大量煙氣導(dǎo)入炭化室,可有效地減少裝煤時(shí)機(jī)側(cè)爐口逸出的煙塵量。但是若高壓氨水噴射卷吸煙氣所形成的負(fù)壓不足或者負(fù)壓過大,則會(huì)造成機(jī)側(cè)爐門處室內(nèi)煙塵外冒或者室外空氣灌入,而這兩種情況在實(shí)際生產(chǎn)中均要避免。因此,如何匹配上升管處高壓氨水噴射的速度與炭化室裝煤時(shí)產(chǎn)氣量之間的關(guān)系,是解決該問題的關(guān)鍵。目前有學(xué)者對(duì)搗固焦?fàn)t煙塵處理進(jìn)行了相關(guān)的研究,羅時(shí)政等[7]開發(fā)了炭化室壓力自動(dòng)調(diào)節(jié)系統(tǒng),并與高壓氨水噴射結(jié)合,解決了結(jié)焦末期出現(xiàn)負(fù)壓的問題,具有一定效果。趙洪波等[8]在集氣管上增加高壓氨水噴灑管道,利用高壓氨水與除塵車聯(lián)合操作,達(dá)到除塵效果。邵豐等[9]探討了高壓氨水噴嘴結(jié)構(gòu)、安裝位置、直徑及壓力對(duì)煙塵疏導(dǎo)效果的作用和規(guī)律,為煙塵疏導(dǎo)工藝的設(shè)計(jì)提供依據(jù)。賈仕文[10]采取增加爐頭煙塵回收系統(tǒng)、設(shè)置集氣管壓力自動(dòng)調(diào)節(jié)系統(tǒng),解決搗固焦?fàn)t煙塵污染。獨(dú)立碳化室直接引射到集氣管容易產(chǎn)生安全危險(xiǎn),采用聯(lián)通管將裝煤時(shí)候的煙氣聯(lián)通到隔壁的高溫碳化室,然后再?gòu)倪@個(gè)碳化室的上升管通過氨水引射到集氣管中。

    本文將某鋼廠搗固焦?fàn)t的碳化室-聯(lián)通管-隔壁碳化室-橋管-集氣管和氨水引射裝置所構(gòu)成的系統(tǒng)簡(jiǎn)化成二維模型,基于FLUENT軟件,運(yùn)用數(shù)值模擬的方法,對(duì)高壓氨水噴射速度與炭化室產(chǎn)氣量之間的關(guān)系進(jìn)行研究。在保證機(jī)側(cè)爐門處氣體靜止的前提下,探討高壓氨水噴射速度與炭化室產(chǎn)氣量之間的關(guān)系,其結(jié)果對(duì)于防止室內(nèi)煙氣的逸出或者室外冷空氣的灌入具有一定的實(shí)際意義。

    1 幾何模型的創(chuàng)建及網(wǎng)格劃分

    1.1幾何模型創(chuàng)建

    模型基于某焦化廠的高度為4.3 m搗固焦?fàn)t。該廠搗固焦?fàn)t采用高壓氨水噴射消煙裝煤,在裝煤時(shí)利用U型管導(dǎo)煙車將該炭化室煙氣導(dǎo)入到相鄰處于結(jié)焦末期的炭化室,進(jìn)而在上升管高壓氨水高速噴射作用下,煙氣被卷吸入集氣管,如圖1所示。

    由于炭化室尺寸較大,因此將其簡(jiǎn)化為二維模型進(jìn)行研究,如圖2所示。模型的簡(jiǎn)化思路是盡量將簡(jiǎn)化后管路的阻力損失及流量與現(xiàn)場(chǎng)測(cè)試一致。對(duì)于4.3 m搗固焦?fàn)t,炭化室長(zhǎng)×寬×高為14.0 m× 0.5 m×4.3 m,其中煤餅距離爐頂?shù)母叨?.3 m,本文取炭化室沿長(zhǎng)度方向的垂直截面及上升管集氣管的垂直截面進(jìn)行研究,其中U型連通管在現(xiàn)場(chǎng)是水平連接兩炭化室,本模型在由三維簡(jiǎn)化成二維的過程中,將U型連通管簡(jiǎn)化為圖2中所示的方式,在簡(jiǎn)化過程中保證局部阻力損失不變。具體的簡(jiǎn)化思路如下:搗固焦?fàn)t在裝煤時(shí),炭化室的焦側(cè)為關(guān)閉狀態(tài),可認(rèn)為煙氣靜止,因此只取炭化室長(zhǎng)度的一半進(jìn)行研究,對(duì)整個(gè)管路的阻力影響不大,即圖中炭化室長(zhǎng)度為7.0 m;煤餅高度為4.0 m,炭化室的凈高則為0.3 m,煙氣只在煤餅上方流動(dòng),因此圖中簡(jiǎn)化的炭化室寬度即為炭化室的凈高0.3 m;為了模擬炭化室煙氣的產(chǎn)生,簡(jiǎn)化了一個(gè)2.0 m長(zhǎng)的產(chǎn)氣口,模擬時(shí)保證該產(chǎn)氣口的流量與實(shí)際煙氣發(fā)生量相同,由于炭化室的尺寸較大,產(chǎn)生的煙氣流速很低,因此在炭化室內(nèi)煙氣阻力損失較少,經(jīng)過測(cè)算此處阻力只占整個(gè)管路阻力的5%左右,進(jìn)行如此簡(jiǎn)化,對(duì)整個(gè)管路的阻力影響并不大;由于U型連通管中有密封裝置、調(diào)節(jié)閥、彎頭等,此處管路的阻力很大,如果按照實(shí)際的U型連通管的截面當(dāng)量直徑,則在簡(jiǎn)化成二維模型后阻力損失較小,經(jīng)過測(cè)算,截面為實(shí)際截面的1/4時(shí),阻力損失與實(shí)際接近,因此取U型連通管的寬度為0.050 m;圖中U型連通管的長(zhǎng)度為現(xiàn)場(chǎng)U型連通管的實(shí)際長(zhǎng)度,即為2.6 m;對(duì)于相鄰處于結(jié)焦末期的炭化室的簡(jiǎn)化類似于裝煤炭化室;上升管處的氨水噴射裝置簡(jiǎn)化為一個(gè)矩形噴射口,噴口直徑為0.015 m,氨水噴射入口位于上升管與集氣管之間,向右下方傾斜45°,且圖中將集氣管進(jìn)行延長(zhǎng),防止在計(jì)算時(shí)產(chǎn)生回流,改善模擬的穩(wěn)定性。

    1.2網(wǎng)格劃分

    文中使用ICEM CFD軟件進(jìn)行網(wǎng)格劃分,采用三角形網(wǎng)格。由于模型中尺寸相差較大,因此采用非均勻網(wǎng)格劃分,對(duì)U型導(dǎo)管、氨水噴射出口等區(qū)域進(jìn)行網(wǎng)格局部加密,以在網(wǎng)格盡量少的情況下保證網(wǎng)格精度。由于模型尺寸較大,因此選取局部網(wǎng)格如圖3所示。

    本模型同時(shí)進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證,在網(wǎng)格數(shù)量分別為3×104,6×104,12×104,20×104在同一工況下考察爐門出口處速度的大小,對(duì)比發(fā)現(xiàn)網(wǎng)格數(shù)為6×104,12×104,20×104之間的差值僅為1%,而3×104與6×104之間的差值為8%,因此綜合考慮可選取網(wǎng)格數(shù)為6×104進(jìn)行最終的模擬計(jì)算。

    2 數(shù)學(xué)模型的建立

    2.1數(shù)學(xué)方程描述

    高壓氨水噴射入,由于氨水速度過高,產(chǎn)生對(duì)周圍煙氣的卷吸,此原理類似于射流泵,因此本文采用文獻(xiàn)[11]中射流泵的模擬方法,基于歐拉方法的兩相流Mixture模型。雖然炭化室中煙氣的產(chǎn)氣量少,流速較低,但是高壓氨水噴射入管道中后流速很高,同時(shí)氨水對(duì)煙氣的卷吸過程均為湍流,因此流動(dòng)過程為湍流,本文采用RNG k-ε湍流模型。

    2.1.1Mixture模型的建立

    Mixture模型可用于模擬各相有不同速度的多相流,但是在小空間尺寸上需假定局部保持平衡,而且相與相之間具有相當(dāng)強(qiáng)的耦合[12],因此該模型可用于射流過程中氣液兩相流動(dòng)的數(shù)值計(jì)算。

    2.1.2連續(xù)性方程

    式中:ρm為混合密度,kg/m3;υm為質(zhì)量平均速度,m/s;下標(biāo)m代表混合相。

    2.1.3動(dòng)量方程

    混合物模型的動(dòng)量方程可以通過對(duì)氣液兩相各自的動(dòng)量方程求和得到

    式中:F為體積力,N;μm為混合相動(dòng)力黏度,N·s/㎡;g為重力加速度,m/s2;αp為氣泡相所占的體積分?jǐn)?shù);αq為流體相所占的體積分?jǐn)?shù);ρq為流體相密度,kg/m3;ρp為氣泡相密度,kg/m3;υdr,q為二次流體相的漂移速度,m/s;υdr,p為二次氣泡相的漂移速度,m/s;下標(biāo)q代表流體相,p代表氣泡相。

    2.1.4氣相體積分?jǐn)?shù)方程

    氣泡相對(duì)于混合速度有一個(gè)漂移速度。根據(jù)氣泡相的連續(xù)性方程,可以得到氣泡相體積分?jǐn)?shù)方程

    2.1.5湍流模型

    采用RNG k-ε模型。與標(biāo)準(zhǔn)k-ε模型相比,RNG k-ε模型通過修正湍流粘度,考慮了平均流動(dòng)中的旋轉(zhuǎn)及旋流流動(dòng)情況[13]。

    k方程

    ε方程

    2.2邊界條件及收斂標(biāo)準(zhǔn)

    2.2.1裝煤炭化室煙氣量計(jì)算

    研究對(duì)象是某鋼廠的2×49孔的4.3 m搗固焦?fàn)t,需要對(duì)裝煤時(shí)產(chǎn)生的煙氣量進(jìn)行計(jì)算,進(jìn)而確定模擬的邊界條件。

    根據(jù)文獻(xiàn)[14]中炭化室裝煤煙氣量的計(jì)算方法,確定本炭化室裝煤煙氣量。4.3 m搗固焦?fàn)t裝煤量為21 t/孔,煉焦過程中的煤氣發(fā)生量為320 m3/t,結(jié)焦時(shí)間炭化室寬為500 mm的取24 h計(jì)算,由此得結(jié)焦時(shí)煙氣產(chǎn)量V0=0.078 m3?s,對(duì)于裝煤過程中產(chǎn)生的煤氣量約為結(jié)焦過程中煤氣發(fā)生量的7倍,則裝煤時(shí)本孔炭化室的產(chǎn)氣量V1=7×V0=0.55 m3?s。因裝煤時(shí)相鄰炭化室處于結(jié)焦末期,產(chǎn)氣量相對(duì)較少,約為初裝煤時(shí)的50%,則相鄰炭化室的產(chǎn)氣量V2=0.039 m3/s,裝煤時(shí)集氣管的流量V=V1+V2=0.6 m3/s。

    2.2.2邊界條件及收斂標(biāo)準(zhǔn)設(shè)置

    將產(chǎn)氣口、氨水噴口邊界設(shè)置為速度入口邊界,產(chǎn)氣量取值為0.2,0.4,0.6,0.8,1.0 m3/s等一系列值,氨水噴射速度取值為20,25,30,35,40 m/s等一系列值,入口湍流取值按水力直徑及湍流強(qiáng)度給定;機(jī)側(cè)爐門邊界設(shè)置為壓力出口邊界,且壓力為大氣壓,即為零壓;集氣管出口設(shè)置為自由出流邊界。

    收斂標(biāo)準(zhǔn):連續(xù)性方程、湍流方程、混合相方程各殘差均10-3。

    3 計(jì)算結(jié)果分析

    3.1流場(chǎng)模擬分析

    取氨水噴射速度為20 m/s,炭化室產(chǎn)氣量為0.2 m3/s,對(duì)該工況模擬得到速度矢量圖分布(圖4)和機(jī)側(cè)爐門及氨水噴射出口局部速度矢量圖(圖5)。

    從圖4可以看出,由于氨水噴出速度較高,射入上升管中時(shí),由于氣、液之間的較大的速度差,以及粘性力作用,導(dǎo)致液體周圍的煙氣被卷吸,產(chǎn)生負(fù)壓,從而使炭化室內(nèi)的煙氣向上升管流動(dòng),由于此時(shí)該炭化室煙氣量較小,導(dǎo)致產(chǎn)生的負(fù)壓過大,導(dǎo)致機(jī)側(cè)爐門的氣流向室內(nèi)流動(dòng)。從圖5可以看到在機(jī)側(cè)爐門處氣體的流動(dòng)方向向內(nèi),表明室外冷風(fēng)在向室內(nèi)灌入。

    取氨水噴射速度為20 m/s,炭化室產(chǎn)氣量為0.4 m3/s,對(duì)該工況模擬得到壓力云圖分布(圖6)和側(cè)爐門局部速度矢量圖(圖7)。

    從圖6可以看出,若產(chǎn)氣量增大到0.4 m3/s時(shí),對(duì)機(jī)側(cè)爐門處的局部速度矢量場(chǎng)進(jìn)行分析。由于此時(shí)該炭化室的產(chǎn)氣量較大,導(dǎo)致高速氨水卷吸氣流所形成的負(fù)壓無法帶走大量的煙氣量,因此導(dǎo)致炭化室機(jī)側(cè)爐門的煙氣是向外流動(dòng)。從圖7可以看出機(jī)側(cè)爐門處的氣流方向向外,表明炭化室內(nèi)的煙氣在向室外逸出。

    3.2模擬結(jié)果分析

    通過數(shù)值模擬得出的計(jì)算結(jié)果,將在不同氨水噴射速度情況下改變炭化室不同產(chǎn)氣量所得到機(jī)側(cè)爐門速度變化置于一張圖中進(jìn)行比較,如圖8所示(圖中約定負(fù)號(hào)表示機(jī)側(cè)爐門氣流方向向內(nèi))。

    以氨水噴射速度為20 m/s且保持不變時(shí)為例,從圖8可以看出:產(chǎn)氣量在0.2 m3/s時(shí),機(jī)側(cè)爐門速度為-2.14 m/s,為負(fù)值,表明氣體流向炭化室內(nèi),這是由于炭化室的產(chǎn)氣量較小,高壓氨水噴射產(chǎn)生的負(fù)壓較大,導(dǎo)致室外氣體向內(nèi)流入;當(dāng)炭化室產(chǎn)氣速度增加為0.4 m3/s時(shí),機(jī)側(cè)爐門的速度為-0.005 5 m/s,接近于零且為負(fù)值,此時(shí)表明機(jī)側(cè)爐門口為微負(fù)壓,符合實(shí)際生產(chǎn)的工作要求;當(dāng)炭化室產(chǎn)氣量繼續(xù)增大為0.6,0.8,1.0 m3/s時(shí),炭化室的產(chǎn)氣量增加,但是此時(shí)高壓氨水噴射形成的負(fù)壓已經(jīng)不足以將炭化室的氣體吸走,因此室內(nèi)氣體向外流動(dòng),機(jī)側(cè)爐門處為正壓,且隨著產(chǎn)氣量的增加機(jī)側(cè)爐門處氣體向外的流動(dòng)速度也增加,分別為1.24,3.22,4.23 m3/s。在不同氨水噴射速度下,炭化室的產(chǎn)氣速度發(fā)生變化時(shí),機(jī)側(cè)爐門的氣體流動(dòng)均表現(xiàn)為類似的情況:在一定氨水噴射速度下,隨著炭化室內(nèi)產(chǎn)氣量的增加,機(jī)側(cè)爐門的速度逐漸增大,且由負(fù)值逐漸變?yōu)檎担砻鳈C(jī)側(cè)爐門的氣體由向內(nèi)逐漸變化為向外流動(dòng)。

    由圖8可看出:對(duì)于計(jì)算的某工況下炭化室裝煤時(shí)的煙氣量為0.6 m3/s,當(dāng)氨水噴射速度小于25 m/s時(shí),機(jī)側(cè)爐門處的氣體速度為正值,表明此時(shí)氨水噴射速度較低,卷吸煙氣不能產(chǎn)生足夠負(fù)壓,因此煙氣流出炭化室外;當(dāng)氨水噴射速度增大至30 m/s以上時(shí),機(jī)側(cè)爐門處的氣體速度為負(fù)值,表明此時(shí)高速氨水噴射卷吸煙氣所形成的負(fù)壓足夠大,導(dǎo)致機(jī)側(cè)爐門處的氣體吸入炭化室內(nèi)。由此可斷定當(dāng)產(chǎn)氣量為0.6 m3/s時(shí),高壓氨水噴射應(yīng)控制速度在25~30 m/s范圍內(nèi),根據(jù)計(jì)算得氨水噴射速度為28.4 m/s。此外,從圖8可以得出當(dāng)炭化室機(jī)側(cè)爐門處氣流速度為0 m/s時(shí),在不同的氨水噴射速度下對(duì)應(yīng)的炭化室產(chǎn)氣量,數(shù)據(jù)見圖9。

    從圖9可以看出,當(dāng)炭化室產(chǎn)氣速度增大時(shí),氨水噴射速度也同步增加,且二者之間的變化幾乎呈現(xiàn)為線性變化。當(dāng)炭化室產(chǎn)氣量為0.4 m3/s,對(duì)應(yīng)的氨水噴射速度為20 m/s時(shí),此時(shí)機(jī)側(cè)爐門氣體速度為零,即煙氣靜止;當(dāng)炭化室產(chǎn)氣量為0.53 m3/s時(shí),氨水噴射速度為25 m/s,隨著炭化室產(chǎn)氣量增加,氨水噴射速度也相應(yīng)增加,這是因?yàn)楫?dāng)煙氣量增加時(shí),此時(shí)炭化室內(nèi)的壓力增大,為了保證炭化室機(jī)側(cè)爐門處氣流靜止,需要增大氨水噴射速度,來增強(qiáng)在氨水的卷吸能力,以在炭化室形成更大的負(fù)壓,進(jìn)而將產(chǎn)生的煙氣帶走。因此在實(shí)際的生產(chǎn)過程中,要根據(jù)產(chǎn)氣量的變化來調(diào)節(jié)氨水噴射的速度,以保證在機(jī)側(cè)爐門處既沒有煙氣外冒,也沒有空氣灌入。

    4 結(jié) 論

    通過對(duì)高壓氨水噴射消煙裝置的二維模型簡(jiǎn)化,采用基于歐拉方法的Mixture模型對(duì)氨水在煙氣中噴射的兩相流過程模擬,得出高壓氨水噴射速度與炭化室不同產(chǎn)氣量之間的關(guān)系,以及相對(duì)應(yīng)的機(jī)側(cè)爐門處氣體流動(dòng)的變化規(guī)律。

    (1)在一定氨水噴射速度下,隨著炭化室內(nèi)產(chǎn)氣量的增加,機(jī)側(cè)爐門處氣體由向內(nèi)逐漸轉(zhuǎn)變?yōu)橄蛲饬鲃?dòng)。

    (2)對(duì)某工況下的裝煤時(shí)產(chǎn)氣量為0.6 m3/s,為保證機(jī)側(cè)爐門氣體靜止,根據(jù)計(jì)算得到氨水噴射速度為28.4 m/s。

    (3)當(dāng)炭化室產(chǎn)氣速度增大時(shí),氨水噴射速度也隨之增加,且二者之間的變化幾乎呈現(xiàn)為線性變化。

    參考文獻(xiàn):

    [1]馬建安.搗固焦?fàn)t裝煤煙塵治理技術(shù)的開發(fā)應(yīng)用[J].煤化工, 2009, 37(6):21-23.

    [2]李秀珍,孟令春,寧芳青,等.焦?fàn)t集氣管壓力解耦控制系統(tǒng)[J].安徽工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2008, 25(2):155-158.

    [3]張赟.搗固型焦?fàn)t裝煤煙氣治理技術(shù)的改進(jìn)[J].機(jī)械管理開發(fā),2006(5):53-55.

    [4]陳怡,陳新.搗固焦?fàn)t推焦裝煤消煙除塵淺談[J].江蘇冶金,2008(6):72-76.

    [5]余子權(quán).中型焦?fàn)t高壓氨水消煙裝煤[J].燃料與化工,1994, 25(4):174-177.

    [6]許伙暖.高壓氨水消煙裝煤技術(shù)的應(yīng)用[J].冶金叢刊,1994(6):33-34.

    [7]羅時(shí)政,劉成雷,許明.焦?fàn)t炭化室壓力自動(dòng)調(diào)節(jié)技術(shù)的應(yīng)用[J].燃料與化工,2014, 2:33-35.

    [8]趙洪波,馮昭.高壓氨水在搗固焦?fàn)t裝煤消煙中的應(yīng)用[J].中國(guó)高新技術(shù)企業(yè),2009,16:26-27.

    [9]邵豐,肖正浩,劉洪春,等.搗固焦?fàn)t高壓氨水煙塵輸導(dǎo)效果的影響因素[J].燃料與化工,2011(5):25-28.

    [10]賈仕文. 5_5 m搗固焦?fàn)t煙塵治理系統(tǒng)改善技術(shù)[J].山東化工,2015, 44(7):115-116.

    [11]王濤,李紅,馮世峰,等.基于Mixture模型對(duì)自吸噴灌泵自吸過程的氣液兩相流模擬[C]//現(xiàn)代節(jié)水高效農(nóng)業(yè)與生態(tài)灌區(qū)建設(shè)(上). 2010.

    [12] Kumar S, Bergada J M.Axial piston pump gooved slipper analysis by CFD simulation of three-dimensional N-S equation in cylindrical coordinates[J]. Computers and Fluids, 2007, 28(1):648-663.

    [13]王福軍.計(jì)算流體動(dòng)力學(xué)分析[M].北京:清華大學(xué)出版社,2004:201.

    [14]王恭鐸.焦?fàn)t裝煤煤煙塵治理的煙氣量計(jì)算[J].燃料與化工,2002, 33(1):14-17.

    責(zé)任編輯:丁吉海

    Numerical Simulation of High Pressure Ammonia Liquor for Smokeless of Stamping Coke Oven

    CHEN Guang1a, CHEN Sai1a, ZHENG Guangming2, DONG Lu1b, YAN Dan1a, WANG Shuang1a

    (1.a Schoolof EnergyandEnvironmentEngineering,bSchoolofCivilEngineering&Architecture,Anhui University of Technology, Ma'anshan 243032, China; 2.Anhui Tongxing Environmental Protecting Engineering Co. Ltd., Ma'anshan 238100, China)

    Abstract:The technology of high pressure ammonia liquor for smokeless was used to control the dust emission in stamping coke oven. In order to study the relationship between the ammonia injection rate and the carbonization chamber gas production, a 4.3 m stamping coke oven equipment of a steel plant was simplified as a two-dimensional model based on the proceeds and Fluent software, the two-phase process of high pressure ammonia liquor injection in rising pipe transmitting smoke was simulated by using Mixture model. The results show that: while the ammonia injection speed is a constant, as chamber gas production increases, the flow direction of the gas flowing in chamber by the side door gradually transforms to outdoor; When the production output of chamber gas is 0.6 m3/s, in order to ensure that the side door gas velocity is zero, ammonia injection rate required is 28.4 m/s; In order to ensure that the velocity of the gas in chamber side door is zero, ammonia injection speed increases almost linearly with the increase of chamber gas production.

    Key words:stamping coke oven; high pressure ammonia liquor; mixture model; numerical simulation

    作者簡(jiǎn)介:陳光(1963-),男,河北玉田人,博士,教授,主要研究方向?yàn)闊峁だ碚撆c系統(tǒng)節(jié)能。

    基金項(xiàng)目:安徽省自然科學(xué)研究項(xiàng)目(KJ2011A056)

    收稿日期:2015-05-21

    文章編號(hào):1671-7872(2015)-03-0245-06

    doi:10.3969/i.issn.1671-7872.2015.03.009

    文獻(xiàn)標(biāo)志碼:A

    中圖分類號(hào):TP229

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    欧美日韩中文字幕国产精品一区二区三区 | 啪啪无遮挡十八禁网站| 少妇 在线观看| 啦啦啦在线免费观看视频4| 国产欧美日韩综合在线一区二区| 中亚洲国语对白在线视频| 极品少妇高潮喷水抽搐| 日韩欧美国产一区二区入口| 757午夜福利合集在线观看| 色综合婷婷激情| 国产精品 欧美亚洲| 在线观看免费午夜福利视频| 国产精华一区二区三区| 久久精品91无色码中文字幕| 国产av精品麻豆| 亚洲,欧美精品.| 一a级毛片在线观看| 亚洲美女黄片视频| 91成人精品电影| 日韩欧美国产一区二区入口| 热re99久久精品国产66热6| 欧美日韩av久久| av天堂久久9| 亚洲avbb在线观看| 国产成人免费观看mmmm| 别揉我奶头~嗯~啊~动态视频| 少妇裸体淫交视频免费看高清 | 女人精品久久久久毛片| 欧美成狂野欧美在线观看| 一本综合久久免费| 亚洲一区二区三区欧美精品| 免费看a级黄色片| 国产高清视频在线播放一区| 天堂俺去俺来也www色官网| 91精品国产国语对白视频| 国产成人av激情在线播放| 精品久久久久久久久久免费视频 | 欧美 日韩 精品 国产| 久久久国产一区二区| 亚洲五月天丁香| 另类亚洲欧美激情| 免费少妇av软件| 欧美 日韩 精品 国产| 久热爱精品视频在线9| 50天的宝宝边吃奶边哭怎么回事| av片东京热男人的天堂| 日本wwww免费看| 色94色欧美一区二区| 国产精品香港三级国产av潘金莲| 国产一区二区三区视频了| 国产视频一区二区在线看| 亚洲avbb在线观看| 婷婷丁香在线五月| 高清av免费在线| 最新在线观看一区二区三区| 可以免费在线观看a视频的电影网站| 免费久久久久久久精品成人欧美视频| www.熟女人妻精品国产| 免费在线观看黄色视频的| 亚洲av欧美aⅴ国产| 亚洲,欧美精品.| 亚洲第一av免费看| 身体一侧抽搐| 又黄又粗又硬又大视频| 亚洲色图综合在线观看| www.精华液| 大码成人一级视频| 精品无人区乱码1区二区| 亚洲av熟女| 美女高潮到喷水免费观看| 精品一品国产午夜福利视频| 亚洲精品中文字幕在线视频| 看黄色毛片网站| 夜夜躁狠狠躁天天躁| 18禁裸乳无遮挡动漫免费视频| 久久久国产成人精品二区 | 久久精品成人免费网站| 69av精品久久久久久| 国产精品自产拍在线观看55亚洲 | 久久久久视频综合| 国产片内射在线| 精品熟女少妇八av免费久了| 久久精品成人免费网站| 国产成人精品久久二区二区91| 成在线人永久免费视频| 久久久久精品人妻al黑| 免费在线观看完整版高清| 免费在线观看完整版高清| 99国产综合亚洲精品| 久久人妻熟女aⅴ| 视频在线观看一区二区三区| 一级作爱视频免费观看| 美女扒开内裤让男人捅视频| 日日摸夜夜添夜夜添小说| 国产免费男女视频| 曰老女人黄片| 高清欧美精品videossex| 两性午夜刺激爽爽歪歪视频在线观看 | 老汉色∧v一级毛片| 久久人人97超碰香蕉20202| 在线观看www视频免费| 黑人操中国人逼视频| 丝袜美足系列| 亚洲精品乱久久久久久| av超薄肉色丝袜交足视频| 亚洲视频免费观看视频| 老汉色av国产亚洲站长工具| 人人妻人人爽人人添夜夜欢视频| 国产午夜精品久久久久久| avwww免费| 久久精品aⅴ一区二区三区四区| 在线视频色国产色| 欧美成人午夜精品| 亚洲一区二区三区不卡视频| 久久国产亚洲av麻豆专区| 久久人妻av系列| 国产片内射在线| 国产欧美日韩一区二区精品| 在线观看免费日韩欧美大片| 中文字幕人妻丝袜一区二区| 热99久久久久精品小说推荐| 日韩欧美一区二区三区在线观看 | 天堂√8在线中文| 床上黄色一级片| 国产高清视频在线观看网站| 99热只有精品国产| 国产亚洲欧美在线一区二区| 淫秽高清视频在线观看| 国产精品三级大全| 久久久久久久久大av| 亚洲av日韩精品久久久久久密| 99精品欧美一区二区三区四区| 中文字幕人妻丝袜一区二区| 久久精品影院6| 国产亚洲欧美在线一区二区| 51午夜福利影视在线观看| 亚洲七黄色美女视频| 男人的好看免费观看在线视频| 亚洲五月天丁香| 亚洲精品久久国产高清桃花| 黄片大片在线免费观看| 亚洲欧美日韩卡通动漫| 久久久久久久精品吃奶| 桃红色精品国产亚洲av| xxxwww97欧美| 亚洲欧美激情综合另类| 夜夜夜夜夜久久久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩高清在线视频| 黄色丝袜av网址大全| 国产精品电影一区二区三区| 欧美高清成人免费视频www| 亚洲中文日韩欧美视频| 欧美不卡视频在线免费观看| 天堂网av新在线| 老司机深夜福利视频在线观看| 国产精品98久久久久久宅男小说| 国产一级毛片七仙女欲春2| 1000部很黄的大片| 午夜福利欧美成人| 国产精品免费一区二区三区在线| 乱人视频在线观看| 亚洲最大成人手机在线| 久久久久精品国产欧美久久久| 少妇的逼水好多| 亚洲精品一卡2卡三卡4卡5卡| 日本免费一区二区三区高清不卡| 欧美bdsm另类| 哪里可以看免费的av片| 欧美成人a在线观看| 97超级碰碰碰精品色视频在线观看| 少妇的丰满在线观看| 尤物成人国产欧美一区二区三区| 人妻久久中文字幕网| 一区二区三区国产精品乱码| 婷婷精品国产亚洲av| 国内少妇人妻偷人精品xxx网站| 在线视频色国产色| 亚洲av成人av| 欧美在线一区亚洲| 美女 人体艺术 gogo| 亚洲不卡免费看| 亚洲午夜理论影院| 精品一区二区三区人妻视频| 国产成人系列免费观看| 欧美一区二区国产精品久久精品| 国产成人影院久久av| 少妇的丰满在线观看| 一区二区三区国产精品乱码| 午夜福利在线观看免费完整高清在 | 在线观看日韩欧美| 在线播放国产精品三级| 国产亚洲精品一区二区www| 日本黄色视频三级网站网址| 在线视频色国产色| 国产一区二区激情短视频| 久久久久久久精品吃奶| 性色avwww在线观看| bbb黄色大片| 老司机午夜福利在线观看视频| 日韩大尺度精品在线看网址| 国产毛片a区久久久久| 欧美激情久久久久久爽电影| av黄色大香蕉| 国产精品99久久久久久久久| 国内精品美女久久久久久| 日韩人妻高清精品专区| 国产成人福利小说| 中文字幕熟女人妻在线| 成年女人永久免费观看视频| 欧美bdsm另类| 亚洲18禁久久av| 国产精品日韩av在线免费观看| or卡值多少钱| 欧美日韩福利视频一区二区| 九色国产91popny在线| 757午夜福利合集在线观看| 老熟妇仑乱视频hdxx| 亚洲激情在线av| 女人十人毛片免费观看3o分钟| 久久草成人影院| 欧美中文综合在线视频| 1024手机看黄色片| 日韩av在线大香蕉| 久久久久久久久大av| 欧美又色又爽又黄视频| 国内少妇人妻偷人精品xxx网站| av福利片在线观看| 成人国产一区最新在线观看| 村上凉子中文字幕在线| 精华霜和精华液先用哪个| 久久亚洲精品不卡| 国产精品一区二区三区四区久久| 岛国视频午夜一区免费看| 88av欧美| 色综合站精品国产| 国产精品一及| 久久精品国产自在天天线| 国产老妇女一区| 欧美+亚洲+日韩+国产| 成人三级黄色视频| 欧美黑人欧美精品刺激| 在线观看一区二区三区| 欧美不卡视频在线免费观看| 精品不卡国产一区二区三区| 成年版毛片免费区| 精品不卡国产一区二区三区| 国产视频一区二区在线看| 91av网一区二区| 一级a爱片免费观看的视频| 美女免费视频网站| 亚洲精品粉嫩美女一区| 观看美女的网站| 精品久久久久久久末码| 1024手机看黄色片| 国产亚洲欧美98| 男女那种视频在线观看| 国产精品永久免费网站| 日本在线视频免费播放| 黑人欧美特级aaaaaa片| 黑人欧美特级aaaaaa片| 欧美精品啪啪一区二区三区| 亚洲 欧美 日韩 在线 免费| 大型黄色视频在线免费观看| 久久婷婷人人爽人人干人人爱| 少妇裸体淫交视频免费看高清| 岛国在线免费视频观看| 在线免费观看不下载黄p国产 | 乱人视频在线观看| 国产高清有码在线观看视频| 国产亚洲精品久久久com| 97超级碰碰碰精品色视频在线观看| 亚洲在线自拍视频| 每晚都被弄得嗷嗷叫到高潮| 久久中文看片网| 亚洲av日韩精品久久久久久密| 日韩av在线大香蕉| 免费在线观看亚洲国产| 亚洲人成网站高清观看| 亚洲国产精品久久男人天堂| 成人精品一区二区免费| 在线观看免费视频日本深夜| 国产av在哪里看| 国产精品日韩av在线免费观看| 99精品欧美一区二区三区四区| 舔av片在线| 久久精品91无色码中文字幕| 欧美日韩一级在线毛片| 精品午夜福利视频在线观看一区| av专区在线播放| 国产探花极品一区二区| www日本在线高清视频| 欧美日韩黄片免| 人人妻人人澡欧美一区二区| 国产真人三级小视频在线观看| 午夜福利欧美成人| 国产精品一区二区三区四区久久| 九九热线精品视视频播放| 嫩草影院精品99| 一边摸一边抽搐一进一小说| 精品欧美国产一区二区三| 日本撒尿小便嘘嘘汇集6| 国产97色在线日韩免费| 久久久国产成人精品二区| 丰满人妻熟妇乱又伦精品不卡| 亚洲在线自拍视频| 国产探花极品一区二区| 手机成人av网站| 国产精品久久久久久久电影 | av国产免费在线观看| 免费av毛片视频| 国产一区二区激情短视频| 在线观看美女被高潮喷水网站 | 国产欧美日韩一区二区精品| 欧美黄色淫秽网站| 国产成人a区在线观看| 国产探花极品一区二区| 老熟妇仑乱视频hdxx| 久久这里只有精品中国| 两个人视频免费观看高清| 精品欧美国产一区二区三| 桃色一区二区三区在线观看| 亚洲av熟女| 国产精品亚洲一级av第二区| 我的老师免费观看完整版| 国产成人系列免费观看| 美女大奶头视频| 丰满人妻一区二区三区视频av | 男人舔女人下体高潮全视频| АⅤ资源中文在线天堂| 亚洲aⅴ乱码一区二区在线播放| 免费搜索国产男女视频| 国产精品久久视频播放| 观看美女的网站| 欧美一区二区亚洲| 久久久国产成人免费| 亚洲人成伊人成综合网2020| 国内久久婷婷六月综合欲色啪| 一级毛片女人18水好多| 亚洲无线观看免费| 成人特级黄色片久久久久久久| 麻豆成人av在线观看| 日韩欧美 国产精品| 欧美+亚洲+日韩+国产| 色综合亚洲欧美另类图片| 丰满乱子伦码专区| 欧美日韩乱码在线| 一本久久中文字幕| 日本一二三区视频观看| 老司机午夜十八禁免费视频| 亚洲在线自拍视频| 欧美高清成人免费视频www| 欧美日韩综合久久久久久 | 亚洲欧美精品综合久久99| 国产精品久久久久久久久免 | 久久欧美精品欧美久久欧美| 国内揄拍国产精品人妻在线| 亚洲精品粉嫩美女一区| 久久久久久国产a免费观看| 欧美成狂野欧美在线观看| 国产精品久久电影中文字幕| 国产 一区 欧美 日韩| 免费看光身美女| 亚洲专区国产一区二区| 怎么达到女性高潮| 天堂网av新在线| 免费在线观看日本一区| 三级男女做爰猛烈吃奶摸视频| 欧美三级亚洲精品| 国产一级毛片七仙女欲春2| 日韩精品青青久久久久久| 18禁黄网站禁片免费观看直播| 久久久久久国产a免费观看| 亚洲av免费在线观看| 90打野战视频偷拍视频| 亚洲五月天丁香| 特大巨黑吊av在线直播| 亚洲第一电影网av| 国产亚洲精品一区二区www| 一区二区三区高清视频在线| 国产免费av片在线观看野外av| 变态另类丝袜制服| 欧美成人一区二区免费高清观看| 国产亚洲欧美98| 亚洲成av人片在线播放无| 九九在线视频观看精品| 女人高潮潮喷娇喘18禁视频| 亚洲片人在线观看| 精品久久久久久久毛片微露脸| 久久性视频一级片| 嫁个100分男人电影在线观看| 在线观看日韩欧美| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品乱码久久久v下载方式 | 国产精品电影一区二区三区| 国产精品乱码一区二三区的特点| 久久国产乱子伦精品免费另类| 亚洲内射少妇av| 午夜精品久久久久久毛片777| 国产又黄又爽又无遮挡在线| 老司机深夜福利视频在线观看| 国产熟女xx| 波多野结衣高清无吗| 一二三四社区在线视频社区8| 欧美激情在线99| 亚洲国产色片| 丁香六月欧美| 男人舔奶头视频| 色综合站精品国产| 久久久色成人| 可以在线观看毛片的网站| 欧美日韩精品网址| 18禁黄网站禁片免费观看直播| 国产一区二区激情短视频| 女人十人毛片免费观看3o分钟| 少妇人妻精品综合一区二区 | 亚洲精品在线观看二区| 国产精品一及| 观看美女的网站| 俺也久久电影网| 级片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 在线观看日韩欧美| 特大巨黑吊av在线直播| 我的老师免费观看完整版| 亚洲av第一区精品v没综合| 日本一本二区三区精品| 丁香欧美五月| 欧美性感艳星| 中文字幕人妻熟人妻熟丝袜美 | 此物有八面人人有两片| 亚洲av一区综合| 欧美乱妇无乱码| 国产精品永久免费网站| 香蕉久久夜色| 国内少妇人妻偷人精品xxx网站| 亚洲无线在线观看| 国产视频内射| 久久久久免费精品人妻一区二区| 午夜日韩欧美国产| 日本在线视频免费播放| 欧美最新免费一区二区三区 | 亚洲av熟女| 国产v大片淫在线免费观看| 夜夜躁狠狠躁天天躁| 黄色日韩在线| 日韩欧美一区二区三区在线观看| 欧美黑人巨大hd| 黄片大片在线免费观看| 精品日产1卡2卡| 熟妇人妻久久中文字幕3abv| 亚洲乱码一区二区免费版| 美女免费视频网站| 欧美在线一区亚洲| 亚洲久久久久久中文字幕| 色老头精品视频在线观看| 丁香欧美五月| 欧美黄色片欧美黄色片| 国产精品久久久久久人妻精品电影| 欧美日韩亚洲国产一区二区在线观看| 国产伦一二天堂av在线观看| 亚洲专区中文字幕在线| 神马国产精品三级电影在线观看| 亚洲av电影在线进入| 国产欧美日韩一区二区精品| 国产精品电影一区二区三区| 老汉色∧v一级毛片| 十八禁网站免费在线| 中亚洲国语对白在线视频| 国产亚洲欧美98| 国产欧美日韩一区二区精品| 真实男女啪啪啪动态图| 国产午夜福利久久久久久| 色综合亚洲欧美另类图片| 久久九九热精品免费| 日韩欧美在线二视频| 亚洲av成人精品一区久久| 亚洲 欧美 日韩 在线 免费| 他把我摸到了高潮在线观看| 亚洲精华国产精华精| 国产久久久一区二区三区| 激情在线观看视频在线高清| 日韩成人在线观看一区二区三区| 69人妻影院| 99热这里只有精品一区| 成年版毛片免费区| 女同久久另类99精品国产91| 3wmmmm亚洲av在线观看| 熟女电影av网| 国产精品久久久久久精品电影| 男人舔女人下体高潮全视频| 亚洲av成人精品一区久久| 国产乱人伦免费视频| 黄片小视频在线播放| 国产亚洲精品久久久com| 真人做人爱边吃奶动态| www日本黄色视频网| 亚洲黑人精品在线| 精品一区二区三区视频在线 | 观看美女的网站| 国产精品永久免费网站| 身体一侧抽搐| 嫩草影院精品99| 桃色一区二区三区在线观看| 香蕉久久夜色| 又粗又爽又猛毛片免费看| 欧美成人性av电影在线观看| 日本一本二区三区精品| 俄罗斯特黄特色一大片| www国产在线视频色| av天堂中文字幕网| 免费看a级黄色片| 国产精品女同一区二区软件 | 精品久久久久久久久久免费视频| 在线观看舔阴道视频| 啦啦啦免费观看视频1| 99在线人妻在线中文字幕| 亚洲精品456在线播放app | 毛片女人毛片| 亚洲欧美精品综合久久99| av天堂中文字幕网| 99国产精品一区二区蜜桃av| 麻豆一二三区av精品| 国产欧美日韩精品亚洲av| 一本精品99久久精品77| avwww免费| 最新在线观看一区二区三区| 999久久久精品免费观看国产| 69av精品久久久久久| 少妇熟女aⅴ在线视频| 日本精品一区二区三区蜜桃| 级片在线观看| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 欧美日韩亚洲国产一区二区在线观看| 女同久久另类99精品国产91| 午夜福利18| av在线天堂中文字幕| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产高清视频在线播放一区| 欧美一区二区国产精品久久精品| tocl精华| 国产毛片a区久久久久| 88av欧美| 性欧美人与动物交配| 人人妻人人看人人澡| 中文字幕av成人在线电影| 亚洲欧美一区二区三区黑人| 97超级碰碰碰精品色视频在线观看| 日韩国内少妇激情av| 欧美+亚洲+日韩+国产| 91在线精品国自产拍蜜月 | 毛片女人毛片| 丁香六月欧美| 亚洲精品日韩av片在线观看 | 久久久久久久午夜电影| 美女 人体艺术 gogo| 精品久久久久久久久久免费视频| 大型黄色视频在线免费观看| 日韩av在线大香蕉| 国产色婷婷99| 97碰自拍视频| 女生性感内裤真人,穿戴方法视频| 亚洲成a人片在线一区二区| 国产精品一区二区三区四区免费观看 | 欧美国产日韩亚洲一区| 日韩欧美 国产精品| 亚洲aⅴ乱码一区二区在线播放| 欧美绝顶高潮抽搐喷水| 两个人视频免费观看高清| 一本综合久久免费| 中文字幕熟女人妻在线| 少妇的逼水好多| 99精品欧美一区二区三区四区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 级片在线观看| 他把我摸到了高潮在线观看| 国产高清videossex| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类 | 黄色丝袜av网址大全| 国产主播在线观看一区二区| 欧美极品一区二区三区四区| 欧美av亚洲av综合av国产av| 久久久国产成人精品二区| 国产男靠女视频免费网站| 女人被狂操c到高潮| 成人三级黄色视频| 蜜桃亚洲精品一区二区三区| 少妇人妻一区二区三区视频| 日本a在线网址| 午夜精品一区二区三区免费看| 婷婷精品国产亚洲av在线| 国产伦人伦偷精品视频| 精品人妻一区二区三区麻豆 | 国产又黄又爽又无遮挡在线| 亚洲一区二区三区色噜噜| 香蕉丝袜av| 国产精品野战在线观看| 久久久久久久久大av| 色综合欧美亚洲国产小说| 国产高清视频在线播放一区| 成人午夜高清在线视频| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久一区二区三区 | 非洲黑人性xxxx精品又粗又长| 午夜两性在线视频| 99国产精品一区二区蜜桃av| 中文资源天堂在线| 亚洲内射少妇av| 国产又黄又爽又无遮挡在线| 久久欧美精品欧美久久欧美| a级毛片a级免费在线| 亚洲国产精品久久男人天堂| 免费看a级黄色片|