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

    煉鎂還原罐內(nèi)傳熱及其強化換熱分析

    2023-08-09 08:55:01張小艷柯亞萍劉浪侯東壯韓子怡劉清江
    關鍵詞:數(shù)值模擬溫度場

    張小艷 柯亞萍 劉浪 侯東壯 韓子怡 劉清江

    摘 要:在碳達峰碳中和背景下,皮江法煉鎂工藝的高能耗問題不容忽視,而造成該工藝高能耗的原因之一是還原罐內(nèi)物料傳熱效率低下。為提高還原罐的傳熱效率,對還原罐內(nèi)的傳熱過程進行數(shù)值模擬,分析球團直徑、球團初始溫度、還原罐內(nèi)徑和肋片參數(shù)對還原罐徑向傳熱性能的影響規(guī)律。結果表明:球團直徑越大、初始溫度越高、還原罐內(nèi)徑越小,越有利于熱量向還原罐內(nèi)部傳遞,其中球團直徑的變化對還原罐中心溫度變化影響最大,直徑為18,25,32 mm和39 mm的球團依次加熱7,5,4 h和3 h后,達到還原反應所需溫度(大于1 300 K);肋片長度越長、寬度越寬、數(shù)量越多,對還原罐中心溫度的提升效果越明顯,其中肋片長度的變化對其影響最大,肋片長度由65 mm以25 mm增量增至115 mm,加熱1 h時,還原罐中心溫度依次增大16.8%、9.6%。

    研究結果為皮江法煉鎂還原罐的換熱及其強化提供理論參考,也為縮短煉鎂周期提供很好的解決思路。

    關鍵詞:還原罐;溫度場;還原程度;強化換熱;數(shù)值模擬

    中圖分類號:TG 156.95

    文獻標志碼:A

    文章編號:1672-9315(2023)04-0647-11

    DOI:10.13800/j.cnki.xakjdxxb.2023.0401

    Analysis of heat transfer and enhancement in

    magnesium reduction tank

    ZHANG Xiaoyan1,KE Yaping1,LIU Lang1,2,HOU Dongzhuang3,HAN Ziyi1,LIU Qingjiang1

    (1.College of Energy Science and Engineering,Xian University of Science and Technology,Xian 710054,China;

    2.Key Laboratory of Western Mine Exploitation and Hazard Prevention,Ministry of Education,

    Xian University of Science and Technology,Xian 710054,China;

    3.

    College of Geology and Environment,Xian University of Science and Technology,Xian 710054,China)

    Abstract:Under the background of peaking carbon neutrality,the high energy consumption of the Pidgeons process of magnesium smelting cannot be ignored,and one of the reasons for the high energy consumption of the process is the low heat transfer efficiency of the materials in the reduction tank.In order to improve the heat transfer efficiency of the reduction tank,the heat transfer process in the reduction tank was numerically simulated,and the influences of pellet diameter,initial temperature,inside diameter and fin parameters on the radial heat transfer performance of the reduction tank were analyzed.The results show:the larger the pellet diameter,the higher the initial temperature and the smaller the inner diameter of the reduction tank,the more favorable the heat transfer to the inside of the reduction tank,in which the change of pellet diameter has the greatest influence on the change of temperature in the center of the reduction tank.The pellets with diameters of 18,25,32 mm and 39 mm were heated successively for 7,5,4 h and 3 h to reach the required temperature of reduction reaction(greater than 1 300 K).The longer the length,wider the width and more fin,the more obvious the improvement effect on the center temperature of the reduction tank is,among which the change of fin length has the greatest influence,the fin length increases from 65 mm to 115 mm in an increment of 25 mm.After 1 h heating,the center temperature of the reduction tank increases by 16.8% and 9.6% in turn.The results? provide a theoretical reference for the heat exchange and strengthening of the reduction tank in? Pidgeons process of magnesium smelting, and also provide a better solution for shortening the cycle of magnesium smelting.

    Key words:reduction tank;temperature field;degree of reduction;enhanced heat transfer;numerical simulation

    0 引 言

    中國豐富的鎂礦儲量為鎂產(chǎn)業(yè)的發(fā)展奠定了基礎[1-2],自2000年以來,中國的產(chǎn)鎂總量占世界生產(chǎn)總量的40%,是全球第一產(chǎn)鎂大國[3-4]。皮江法煉鎂工藝(硅熱還原法)是中國最主要的金屬鎂生產(chǎn)工藝,但該工藝存在還原周期長導致的能耗大問題[5-6],是單位能耗最高的有色冶金行業(yè)之一[7-8]。在碳達峰碳中和背景下,為進一步降低皮江法煉鎂的能耗,重點應放在優(yōu)化產(chǎn)業(yè)結構及提高能源利用效上[9],而提高能源利用率,應該從皮江法煉鎂工藝的傳熱機理等方面深入研究尋找突破口。

    目前對皮江法煉鎂工藝的研究主要集中在還原反應條件的優(yōu)化方面。唐祁峰等重點介紹碳熱還原煉鎂技術的動力學研究進展,指出目前碳熱還原法煉鎂要實現(xiàn)工業(yè)化生產(chǎn)還需解決的問題[10]。穆曉輝等介紹國內(nèi)外不同還原劑的熱法煉鎂研究進展,并分析不同還原劑對煉鎂過程的影響[11]。GUO等研究用鋁熱法從低品位菱鎂礦中直接提取鎂的還原過程,結果表明該方法的鎂提取率高于常規(guī)煉鎂工藝[12]。謝興以單質(zhì)鎂生產(chǎn)過程產(chǎn)生的熵増作為主要研究手段應用,分析熱過程和熱設備中有效能損失的動力學原因[13]。鄧軍平等通過對水泥中摻加適量鎂渣不僅提高了復合水泥的強度,還一定程度上解決了鎂渣粉塵污染問題[14]。

    BUGDAYC等研究還原劑類型和還原劑量對皮江法煉鎂工藝的影響[15]。

    WADA等開發(fā)了一種新型的微波驅(qū)動皮江法煉鎂工藝[16]。WU等采用10 Pa壓力下耦合電熱場制備金屬鎂,該方法有效地降低了通過常規(guī)方法提取鎂所需的反應裝置溫度[17]。CHE等提出一種在流動氬氣中進行的非真空生產(chǎn)鎂的新工藝,結果表明新工藝可縮短煉鎂[18]。喬錦華通過對球團中硅鐵的不同配比進行試驗,分析硅含量對粗鎂產(chǎn)量的影響[19]。

    大量學者對皮江法煉鎂工藝的研究中發(fā)現(xiàn),該工藝存在還原罐內(nèi)球團溫升速度過慢或者存在冷芯使罐內(nèi)球團反應不充分等問題,從而影響單位時間產(chǎn)鎂率和產(chǎn)鎂量。為探究其中的原因,李波分析高硅白云石產(chǎn)鎂率低的原因,并通試驗驗證皮江法煉鎂工藝的原材料采用高硅白云石是可行的[20]。ZHANG等對硅熱還原過程的化學動力學機理進行試驗研究,建立包含化學反應的三維傳熱非定常數(shù)值方法,并對此進行了驗證[21]。

    FU等提出將白云石、硅鐵、螢石和粘結劑等粉末材料混合制成球團的一種生產(chǎn)鎂的新工藝[22]。除此之外也有不少學者對煉鎂還原罐結構等進行研究。傅大學等為解決皮江法還原罐內(nèi)傳熱慢的問題,采用數(shù)值方法研究預制球團在還原罐內(nèi)的傳熱規(guī)律[23]。ZHANG等開發(fā)利用高溫碳化硅顆粒增強還原罐內(nèi)傳熱的循環(huán)利用新工藝,并對此進行數(shù)值模擬,結果表明該技術可顯著提高傳熱效率,縮短生產(chǎn)周期,降低生產(chǎn)成本[24]。鄧軍平等設計一種有利于加快物料還原反應速率的新型內(nèi)熱式煉鎂反應器,

    該反應器具有傳熱速度快,內(nèi)部溫度梯度小的特點[25]。

    在對皮江法煉鎂工藝進行降耗研究過程中,主要目標是大力提高資源和能源的利用率。目前造成皮江法煉鎂工藝高能耗的原因之一是真空熱還原階段還原罐內(nèi)物料傳熱效率低下,從而導致還原周期(一般為8~12 h)較長、產(chǎn)鎂率低等問題。為縮短煉鎂周期,提高還原罐的傳熱速率,

    以還原罐為研究對象,

    基于ANSYS仿真模擬軟件與傳熱學基本理論,對還原罐內(nèi)傳熱過程進行研究,利用溫度和還原反應程度變化,獲知各影響因素下還原罐傳熱性能的變化特性。為進一步降低皮江法煉鎂工藝的能耗,縮短還原周期,實現(xiàn)煉鎂綠色化、低碳化目標提供理論參考。

    1 數(shù)值模擬方法

    1.1 物理模型

    皮江法煉鎂過程是將還原罐內(nèi)抽成真空,然后在1 200 ℃高溫下進行固相反應還原出鎂單質(zhì),得到粗鎂。還原罐內(nèi)部傳熱過程如圖1所示,主要包括:(a)罐壁與球團間的輻射;(b)球團與球團之間的導熱;(c)罐壁與球團間的導熱;(d)球團之間的輻射和球團自身的化學反應吸熱5種形式。還原罐為壁厚30 mm,長度3 m的圓筒,核桃狀球團填滿罐內(nèi)。在大量文獻研究中,學者們多將還原罐內(nèi)球團簡化為一個整體,由于還原罐內(nèi)球團與球團之間存在輻射傳熱,將球團視作一個整體去建模會引起較大誤差,所以將球團簡化為球狀球團,通過等效體積法可計算得到球團直徑。由于還原罐長度較長,直徑較小且還原爐內(nèi)溫度均勻,認為熱量僅沿還原罐徑向傳遞,故可在軸向取單層球團厚度建立三維物理模型,如圖2所示,并且為了后續(xù)表述簡潔化,模型幾何參數(shù)尺寸符號見表1。由于傳熱過程十分復雜,為了便于模擬計算結合實際做出如下假設。

    1)罐壁向球團傳熱、球團由外向內(nèi)的傳熱均可視為一維徑向傳熱。

    2)在整個煉鎂周期內(nèi),未發(fā)生化學反應之前無內(nèi)熱源,達到反應溫度后球團持續(xù)發(fā)生還原反應并忽略鎂蒸氣對傳熱過程的影響。

    3)還原罐外側為定壁溫條件。

    1.3 計算方法及網(wǎng)格劃分

    采用

    Fluent模擬軟件對還原罐傳熱過程進行數(shù)值模擬,采用

    Solidworks軟件對還原罐及球團(還原罐外徑為330 mm,內(nèi)徑為270 mm,球團直徑25 mm)進行建模,完成建模后利用Workbench的Mesh模塊進行網(wǎng)格劃分(圖3),對罐內(nèi)球團進行布爾減操作,畫出流體域部分,整個模型的網(wǎng)格劃分采用結構性四面體和六面體網(wǎng)格,為計算精度要求對罐內(nèi)球團部分進行加密處理。同時,在還原罐傳熱過程的數(shù)值模擬中,采用壓力求解器和非穩(wěn)態(tài)計算法,開啟能量方程和DO輻射模型,為提高收斂速度與計算精度,采用有限元體積法,壓力耦合器采用Simple,松弛因子保持不變。

    1.4 網(wǎng)格無關性驗證

    在數(shù)值模擬中,網(wǎng)格數(shù)量對計算精度有一定影響,為了提高計算精度,針對還原罐內(nèi)10 h的溫升過程,分別進行了網(wǎng)格無關性驗證及步長獨立性驗證,以確保模擬結果的準確性。圖4為罐中心溫度隨時間變化。

    從圖4可以看出,網(wǎng)格數(shù)為288 665時比119 669的溫度略微升高,和211 347對比幾乎沒有區(qū)別,所以選擇網(wǎng)格數(shù)為211 347進行數(shù)值模擬計算;時間步長為2,1,0.5 s的計算結果基本一致,所以選取2 s的時間步長。

    1.5 模型驗證

    為驗證模擬計算結果和源項UDF編譯的準確性,采用文獻[27]的結果進行驗證。所驗證還原罐模型的尺寸為內(nèi)徑270 mm,外徑330 mm,壁厚30 mm,罐內(nèi)球團直徑17.8 mm,罐體和球團的物性參數(shù)及邊界條件的設置與文獻一致,其中罐體加熱溫度為1 473 K,球團初始溫度為298.15 K。根據(jù)模擬結果計算還原反應程度與文獻模擬結果進行對比,圖5是數(shù)值模擬與文獻模擬結果的對比情況。其中還原反應程度定義為已反應球團體積和總球團體積之比。通過對比可知,還原罐內(nèi)溫度模擬結果與文獻模擬結果接近,誤差在10%內(nèi),最大誤差在2 h時為8.84%,驗證了還原罐模型和模擬方法的準確性和可靠性。

    2 結果與分析

    2.1 球團直徑對還原罐內(nèi)傳熱的影響

    在還原罐外徑為330 mm,內(nèi)徑為270 mm,球團初始溫度為300 K,球團直徑分別為18,25,32,39 mm的情況下,對還原罐內(nèi)的溫度場進行數(shù)值模擬。從圖6可以看出,不同球團直徑時,還原罐中心溫度在不同時間段升溫速率差別明顯,在0~1 h內(nèi)直徑為39 mm的球團溫升速度最快,在7 h以后中心溫度逐漸趨于穩(wěn)定。當球團直徑由d=18 mm以7 mm增量依次增加至39 mm時,加熱1 h的中心溫度由659 K分別增至934,1 169,1 266 K,增幅分別為275,235 K和97 K,可見隨著球團直徑的增加,還原罐中心溫度明顯升高,但增加幅度逐漸減小,其中d=18 mm增長至d=25 mm時的增幅最大。從加熱時間來看,d=18 mm的球團在加熱7 h后才達到還原反應所要求的溫度(大于1 300 K),而d=25 mm、d=32 mm和r=39 mm的球團分別需要加熱5,4 h和3 h,可以看出球團直徑越大加熱達到化學反應溫度的時間越短。在達到反應條件后,球團就會發(fā)生還原反應進行吸熱,從而影響熱量向罐內(nèi)的傳遞,還原罐內(nèi)中心溫度的升幅逐漸減小最終趨于平穩(wěn)。

    從圖7可以看出,在相同加熱時間下,隨著球團直徑的增加,溫度分布更加均勻。在加熱1 h時,球團直徑由18 mm每增加7 mm,罐內(nèi)最大溫差(罐內(nèi)R=0 mm與R=135 mm處球團的溫差)由774 K分別減小至520,278 K和195 K,可以看出球團直徑增加會使罐內(nèi)最大溫差減小,但減幅是逐漸降低的,其中d=18 mm增長至d=25 mm時的減幅最大。此外圖中可以明顯看出,直徑為18,25,32 mm的球團在加熱2h后的溫度分布曲線分別與直徑為25,32,39 mm的球團在加熱1 h后的溫度分布曲線十分接近,可見增大球團直徑可以有效提高傳熱速率。

    2.2 球團初始溫度對還原罐內(nèi)傳熱的影響

    在還原罐外徑為330 mm,內(nèi)徑為270 mm,球團直徑為18 mm時,球團初始溫度分別為300,800,1 000 K和1 200 K的情況下對還原罐內(nèi)的溫度場進行數(shù)值模擬。從圖8可以看出,在相同加熱條件下,提高球團初始溫度,可使還原罐中心溫度更快達到反應溫度。當球團初始溫度為300? K時,加熱7 h還原罐中心溫度才能達到1 373 K,從而使最內(nèi)部球團開始反應;當球團初始溫度提高到800 K時加熱7 h,與初溫300 K加熱7 h的溫度進行對比,只是略有升高;當球團初始溫度提高到1 000 K時,加熱只需6 h就可使罐中心溫度可提升至1 370 K;當球團初始溫度提高到1 200 K時,使罐中心溫度升至1 370 K僅需加熱4 h。從以上數(shù)據(jù)可以看出,提高球團初始溫度,可以縮短還原罐內(nèi)部球團達到化學反應所需溫度的時間,但對于球團初始溫度的提升越大則效果越好,并且在初始溫度300~800 K范圍內(nèi)對縮短達到化學反應條件時間的作用可以忽略不計,但升至1 000 K以上會有明顯效果。造成該現(xiàn)象的原因是球團升溫過程中在達到反應條件后,球團就會發(fā)生還原反應而吸熱,則向內(nèi)部傳遞的熱量會低于未發(fā)生化學反應時,只有使球團初始溫度越接近化學反應溫度,才能更快縮短煉鎂周期。

    圖9顯示不同球團初始溫度下還原罐內(nèi)的溫度分布云圖,從圖9可以看出,不同初始溫度下還原罐內(nèi)的溫度分布規(guī)律,在加熱1 h時,初始溫度為1 200 K的還原罐內(nèi)高溫球團數(shù)量明顯多于初始溫度為300 K;在加熱2 h時,初始溫度為300 K的還原罐內(nèi)依然存在代表1 000 K的綠色球團,而初始溫度為1 200 K的還原罐內(nèi)最低溫度已經(jīng)達到1 300 K,在加熱4 h時,初始溫度為1 200 K的還原罐內(nèi)球團均已達到反應溫度,與初始溫度為300 K的對比很明顯可以縮短還原罐的加熱周期。

    2.3 還原罐內(nèi)徑對還原罐內(nèi)傳熱的影響

    球團初始溫度為300 K,球團直徑為25 mm,在還原罐內(nèi)徑分別為220,270,320 mm和370 mm,壁厚30 mm的情況下對還原罐內(nèi)的溫度場進行數(shù)值模擬。從圖10可以看出,還原罐內(nèi)徑的增大雖然增加了還原罐受熱面積,但同時因球團裝料量的增加使得熱阻變大,不利于熱量向罐內(nèi)的傳遞。在相同加熱時長下,還原罐直徑越小罐中心溫度越高。當還原罐內(nèi)徑由220 mm每增加50 mm,加熱1 h時,還原罐中心溫度由1 092 K依次減小至934,714,501 K,加熱3 h時,還原罐中心溫度由1 366 K依次減小至1 286,1 180,1 064 K。以還原罐中心溫度加熱到1 400 K為節(jié)點,從加熱時間來看,內(nèi)徑為220,270 mm和320 mm所需要的時間分別為4,6,8 h,而當還原罐內(nèi)徑增加至370 mm時,加熱8 h仍未達到1 400 K。從以上數(shù)據(jù)可明顯看出增加還原罐內(nèi)徑嚴重影響還原罐內(nèi)傳熱。

    從圖11可以看出,在相同加熱時間下,隨著還原罐直徑的增加,還原反應程度逐漸降低,且還原罐內(nèi)徑越大罐中心溫度達到反應溫度所需的時間越長。加熱1 h時,內(nèi)徑為220,270,320 mm和370 mm的還原罐內(nèi)還原反應程度分別為22.6%、19.2%、15.4%和11.7%,在這一時間段內(nèi)還原罐內(nèi)球團處于溫升階段,所以隨著還原罐內(nèi)徑的增加還原程度變化差異不是很明顯;還原罐內(nèi)徑由220 mm每增加50 mm,加熱4 h時,還原罐內(nèi)還原反應程度由92.1%依次減小至69%、49.11%、37.8%,由此可明顯看出增大還原罐內(nèi)徑會嚴重影響煉鎂周期。

    2.4 肋片參數(shù)對加肋片的還原罐內(nèi)傳熱的影響

    擴大還原罐直徑可提高還原罐球團填充量,但會嚴重影響煉鎂周期,若能解決傳熱慢的問題則能提高單罐產(chǎn)鎂量。針對大型還原罐傳熱慢的問題借鑒文獻[28]提出的適用于豎罐皮江法煉鎂的肋片罐結構設計方案,將其中空罐及肋片分布和數(shù)量進行了修改,肋片罐結構示意圖如圖12所示。通過數(shù)值模擬的方法研究肋片參數(shù)對還原罐內(nèi)傳熱的影響規(guī)律。與內(nèi)徑270 mm的還原罐對比,內(nèi)徑370 mm罐體的球團裝料量提升1.87倍。選擇直徑370 mm的還原罐,球團直徑為25 mm,初始溫度為300 K,分別對含有不同肋長、肋寬以及數(shù)量的還原罐內(nèi)溫度場進行數(shù)值模擬,運用控制變量法分析肋長、肋寬以及數(shù)量對還原罐傳熱的影響規(guī)律。

    從圖13可以看出,在同一加熱時長下,加肋片的罐體中心溫度均大于無肋片罐,且肋長越長罐中心溫度越高,肋長等量增加對應中心溫度也近似等量升高。與無肋片罐相比,肋寬為15 mm、肋長分別為65,90,115 mm的肋片罐,加熱1 h時,還原罐中心溫度由501 K依次增加至581,679,744 K,增量分別為80,98,65 K。隨著加熱時間的延長,各肋長的還原罐中心溫度差值逐漸減小,加熱8 h時,還原罐中心溫度依次為1 360,1 378,1 389,1 398 K,增量均小于20 K,分析認為,隨著罐內(nèi)溫度的升高,與熱源的溫差逐漸減小,傳熱速率下降,罐內(nèi)溫度也逐漸趨于平緩。

    從圖14可以看出,同一加熱時長下,有肋片的還原罐罐中心溫度均高于無肋片還原罐,且肋片越寬,罐體中心溫度越高。與無肋片罐相比,肋長為115 mm,寬度分別為10,15,20 mm的肋片罐,加熱1 h時,還原罐中心溫度由501 K依次增至691,744,801 K,增幅分別為53 K和57 K,增幅略有增大。隨著加熱時間的延長,各肋寬的還原罐中心溫度差值逐漸減小,在7 h時,差值均小于15 K。

    從圖15可以看出,在同一加熱時長下,肋片數(shù)越多,罐體中心溫度越高,且隨著加熱時間的延長,不同數(shù)量的肋片罐中心溫度差值逐漸減小。肋片數(shù)量由n=4每增加4時,加熱1h時的中心溫度由679 K分別增至769,850 K,增幅分別為90 K和82 K,可見隨著肋片數(shù)量的增加,還原罐中心溫度明顯升高,但增加幅度略微減小,其中n=4增長至n=8時的增幅最大。

    雖然肋片數(shù)的增加可有效提高罐中心溫度,但與增加肋長和肋寬相比,肋片數(shù)的增加對罐內(nèi)球團的填充量有很大影響。圖16為不同肋片數(shù)時還原罐內(nèi)的溫度分布云圖,從球團填充量來看,肋片數(shù)越多,罐內(nèi)球團填充量則越少,在肋片數(shù)n=4時,罐內(nèi)球團可以最大程度上充滿罐體,在n=8時,球團數(shù)量減小,在n=12時,球團數(shù)量最少。從溫度分布來看,在同一加熱時長下,肋片數(shù)越多則罐內(nèi)溫度分布越均勻。在加熱1 h,肋片數(shù)n=4時,處在750~800 K的低溫球團數(shù)量(藍色區(qū)域)明顯多于肋片數(shù)n=8和n=12,隨著加熱時間的延長,這一規(guī)律保持不變,均為肋片數(shù)少的罐中心溫度較低的球團數(shù)量多于肋片數(shù)較多的。在加熱6 h時,肋片數(shù)n=4的還原罐內(nèi)溫度大于1 450 K的高溫球團(紅色區(qū)域)數(shù)量明顯少于肋片數(shù)n=8和n=12,尤其是n=12時,除了肋片罐中心的球團外其余球團溫度均處于1 450 K以上。

    2.5 各影響因素下還原罐傳熱性能綜合評估

    表2為各影響因素下加熱1~3 h(煉鎂周期內(nèi)溫升階段)時還原罐中心溫度的變化情況,罐中心溫度變化值及變化率均為一個煉鎂周期(8 h)計算所得。從表2可以看出,球團直徑越大、初始溫度越高、還原罐內(nèi)徑越小,越有利于熱量向罐體內(nèi)部的傳遞。在球團直徑以7 mm增量由18 mm增加至39 mm的過程中,加熱1 h時,還原罐中心溫度依次升高34.9%、17.4%、7.4%。在球團初始溫度由300 K升高至1 200 K的過程中,由于300 K至800 K范圍內(nèi)數(shù)據(jù)差異不大,所以直接對比了300 K升高至800 K的結果。當球團初始溫度由300 K依次增至8 00,1 000,1 200 K,加熱1 h時,還原罐中心溫度依次增大了36.9%、14.8%、19.1%。在還原罐直徑由220 mm增加至370 mm的過程中,加熱1 h時,直徑每增加50 mm,還原罐中心溫度依次減小14.5%、23.6%、29.8%。

    圖17更直觀地顯示出球團直徑、球團初始溫度與還原罐中心溫度呈正相關,還原罐內(nèi)徑與還原罐中心溫度呈負相關,并且隨著加熱時間的延長,還原罐中心溫度隨球團直徑、球團初始溫度、還原罐內(nèi)徑變化而變化的幅度逐漸減小。其中對還原罐中心溫度變化影響最大的是球團直徑的變化,其次是還原罐內(nèi)徑,最后是球團初始溫度。此處需對此進行說明,圖中的球團初始溫度影響因素的數(shù)據(jù)并非等量遞增的,由于球團初始溫度300~800 K的差值較大,所以使球團初始溫度影響因素下1 h的數(shù)據(jù)高于球團直徑影響因素下1 h處的數(shù)據(jù)。

    表3為肋片罐各影響因素下加熱1~3 h時的還原罐中心溫度變化情況,從表3可以看出,肋片長度越長、寬度越寬、數(shù)量越多,越有利于熱量向罐體內(nèi)部的傳遞。在肋片長度由65 mm增加至115 mm的過程中,每增加25 mm,加熱1 h時,還

    原罐中心溫度依次增大16.8%、9.6%。在肋片寬度由10 mm增加至20 mm的過程中,每增加5 mm,加熱1 h時,還原罐中心溫度依次增大7.7%、7.6%。在肋片數(shù)由4個增加至12個的過程中,每增加4個,加熱1 h時,還原罐中心溫度依次增大13.3%、10.6%。

    圖18更直觀地顯示出肋片長度、肋片寬度以及肋片數(shù)量與還原罐中心溫度呈正相關,并且同樣呈現(xiàn)出隨著加熱時間的延長,還原罐中心溫度隨肋片長度、寬度、數(shù)量變化而變化的幅度逐漸減小。其中肋片長度的變化對還原罐中心溫度的影響最大,其次是肋片數(shù)量,肋片寬度對還原罐中心溫度的影響最小。

    2.6 優(yōu)化分析

    研究的目的是提高還原罐的傳熱效率,縮短煉鎂周期從而提高煉鎂廠的經(jīng)濟效益。通過分析結果可以看出,增大球團直徑、提高球團初始溫度、減小還原罐內(nèi)徑以及還原罐內(nèi)增生肋片的方式均能有效提高還原罐傳熱效率,在實際工業(yè)生產(chǎn)的設備改造過程中,還要考慮經(jīng)濟性。在現(xiàn)有條件下,要提高單罐產(chǎn)鎂量,首先應該選擇增大罐體直徑,在大直徑的還原罐內(nèi)壁增設肋片并選擇大直徑的球團進行還原反應,原因是在提高了單罐產(chǎn)鎂量的基礎上,增大還原罐傳熱效率,從而使煉鎂周期不變甚至縮短,這大大提高了能源利用率。

    3 結 論

    1)

    當球團直徑增加至39 mm時,加熱3 h可達到還原反應所需溫度,此外將球團初始溫度提升至1 200 K時,加熱4 h可達到1 370 K,表明提高球團直徑和初始溫度,可以使還原罐中心更快達到反應溫度。

    2)還原罐內(nèi)徑越大還原反應程度越低,達到完全反應時間也越長。還原罐內(nèi)徑由220 mm增加至370 mm,在加熱4 h時,還原反應程度由92.1%減小為37.8%。通過在還原罐內(nèi)壁加肋片可以有效解決大型還原罐傳熱慢的問題,加肋片的罐體中心溫度在同一加熱時長下均大于無肋片的,肋片長度越長、寬度越寬、數(shù)量越多,還原罐內(nèi)傳熱效率越高。

    3)綜合分析表明球團直徑的變化對還原罐中心溫度變化影響最大,其次是還原罐直徑,球團初始溫度對還原罐中心溫度影響最小。還原罐內(nèi)壁加肋片有利于熱量向罐中心傳遞,肋片長度的變化對還原罐中心溫度影響最大,其次是肋片數(shù)量,肋片寬度對還原罐中心溫度影響最小。

    4)在實際工業(yè)生產(chǎn)的設備改造過程中,考慮經(jīng)濟性的前提下,要提高單罐產(chǎn)鎂量,首先應該選擇增大罐體直徑,對其內(nèi)壁增設肋片并選擇大直徑的球團進行還原反應。

    參考文獻(References):

    [1] 于旭光,邱竹賢.鎂工業(yè)生產(chǎn)及應用的現(xiàn)狀和展望[J].材料與冶金學報,2003(3):189-192.

    YU Xuguang,QIU Zhuxian.Present situation and prospect of magnesium production and applications[J].Journal of Materials and Metallurgy,2003(3):189-192.

    [2]RAMAKRISHNAN S,KOLTUN P.Global warming impact of the magnesium produced in China using the Pidgeon process[J].Resources,Conservation and Recycling,2004,42(1):49-64.

    [3]馬登民.鎂及鎂合金產(chǎn)業(yè)發(fā)展研究[J].世界有色金屬,2019,19:165-168.

    MA Dengmin.Research on the development of magnesium and magnesium alloy industry[J].World Nonferrous Metals,2019,19:165-168.

    [4]林如海,孫前.走向新時代的中國鎂工業(yè)[J].中國有色金屬,2019(3):38-41.

    LIN Hairu.Chinese magnesium industry towards a new era[J].China Nonferrous Metals,2019(3):38-41.

    [5]朱祖武,曹黎華.國內(nèi)皮江法煉鎂技術的缺陷與改進途徑[J].南昌高專學報,2007(2):92-94.

    ZHU Zuwu,CAO Lihua.The deficiency in pidgeon process of mg-smelting & methods lmprovements[J].Journal of Nanchang College,2007(2):92-94.

    [6]GAO F,NIE Z,WANG Z,et al.Assessing environmental impact of magnesium production using Pidgeon process in China[J].Transactions of Nonferrous Metals Society of China,2008,18(3):749-754.

    [7]任虎奎,任建勛.我國金屬鎂冶煉技術現(xiàn)狀與發(fā)展趨勢[J].陜西煤炭,2011,30(1):36-37,69.

    REN Hukui,REN Jianxun.Present situation and development trend of magnesium smelting technology in China[J].

    Shaanxi Coal Mine,

    2011,30(1):36-37,69.

    [8]CHERUBINI F,RAUGEI M,ULGITAIl S.LCA of magnesium production[J].Resources,Conservation and Recycling,2008,52(8-9):1093-1100.

    [9]尚梅,李強.區(qū)域經(jīng)濟發(fā)展與能源消費間關系研究[J].西安科技大學學報,2015,35(6):787-793.

    SHANG Mei,LI Qiang.Relationship between Chinas regional economic development and energy consumption[J].Journal of Xian University of Science and Technology,2015,35(6):787-793.

    [10]唐祁峰,高家誠,陳小華,等.真空碳熱還原制鎂工藝的研究進展[J].材料導報,2014,28(23):64-67,85.

    TANG Qifeng,GAO Jiacheng,CHEN Xiaohua,et al.Research progress of magnesium production by carbothermic reduction at vacuum[J].Materials Reports,2014,28(23):64-67,85.

    [11]穆曉輝,羅洪杰,王耀武,等.熱法煉鎂還原劑研究進展[J].輕金屬,2014(2):45-49.

    MU Xiaohui,LUO Hongjie,Wang Yaowu,et al.Research progress on agent in thermal magnesium reduction process[J].Light Metals,2014(2):45-49.

    [12]GUO J,LI X,ZHANG T,et al.Comparison of extraction behavior of magnesium from magnesite/magnesia by aluminothermic process in flowing argon[J].Journal of Sustainable Metallurgy,2022,8(4):1756-1768.

    [13]謝興.新法鋁熱真空還原鎂工藝及熱力學分析[D].沈陽:東北大學,2014.

    XIE Xing.Magnesium production by aluminothermic reduction in vacuum and thermodynamic analysis[D].Shenyang:Northeastern University,2014.

    [14]鄧軍平,陳新年,郭一萍.鎂渣和礦渣對復合水泥性能的影響[J].西安科技大學學報,2008,28(4):735-739.

    DENG Junping,CHEN Xinnian,GUO Yiping.Effect of magnesium slag and slag on properties of composite cement[J].Journal of Xian University of Science and Technology,2008,28(4):735-739.

    [15]BUGDAYCI M,TURAN A,ALKAN M,et al.Effect of reductant type on the metallothermic magnesium production process[J].High Temperature Materials and Processes,2018,37(1):1-8.

    [16]WADA Y,F(xiàn)UJII S,SUZUKI E,et al.Smelting magnesium metal using a microwave pidgeon method[J].Scientific Reports,2017,7(1):46512.

    [17]WU H,ZHAO P,JING M,et al.Magnesium production by a coupled electric and thermal field[J].Vacuum,2021,183:109822.

    [18]CHE Y,ZHANG C,SONG J,et al.The silicothermic reduction of magnesium in flowing argon and numerical simulation of novel technology[J].Journal of Magnesium and Alloys,2020,8(3):752-760.

    [19]喬錦華,梁國軍.球團中硅含量對粗鎂產(chǎn)量的影響[J].冶金管理,2019(3):38-40.

    QIAO Jinhua,LIANG Guojun.Effect of silicon content in pellets on crude magnesium yield[J].China Steel Focus,2019(3):38-40.

    [20]李波,褚丙武.高硅白云石皮江法煉鎂實驗[J].有色礦冶,2015,31(6):32-35,38.

    LI Bo,CHU Bingwu.Experiment of magnesium smelting with high silica dolomite by Pidgeon method[J].Non-Ferrous Mining and Metallurgy,2015,31(6):32-35,38.

    [21]ZHANG C,CHU H,GU M,et al.Experimental and numerical investigation of silicothermic reduction process with detailed chemical kinetics and thermal radiation[J].Applied Thermal Engineering,2018,135:454-462.

    [22]FU D X,ZHANG T A,GUAN L K,et al.Magnesium production by silicothermic reduction of dolime in pre-prepared dolomite pellets[J].JOM,2016,68(12):3208-3213.

    [23]傅大學,張廷安,豆志河,等.球團預制-硅熱還原煉鎂還原罐內(nèi)傳熱[J].東北大學學報(自然科學版),2018,39(5):643-648.

    FU Daxue,ZHANG Tingan,DOU Zhihe,et al.Heat transfer in a retort for mg-extraction by pre-prepared pellet silicothermic process[J].Journal of Northeastern University(Natural Science),2018,39(5):643-648.

    [24]ZHANG C,YAN W,CHE Y,et al.Experimental and numerical studies on heat transfer enhancement in the silicothermic reduction process[J].Energy Reports,2022(8):6999-7009.

    [25]鄧軍平,王曉剛.內(nèi)熱式煉鎂工藝實驗研究[J].材料導報,2014,28(S2):179-181,204.

    DENG Junping,WANG Xiaogang.Experimental research on the inner thermal magnesium reduction process[J].Materials Reports,2014,28(S2):179-181,204.

    [26]楊康定,陳群,任建勛.煉鎂還原罐內(nèi)傳熱與化學反應的耦合特性[J].清華大學學報(自然科學版),2009,49(5):755-758.

    YANG Kangding,CHEN Qun,REN Jianxun.Coupled heat transfer and chemical reactions in magnesium production retorts[J].Journal of Tsinghua University(Science and Technology),2009,49(5):755-758.

    [27]馮俊小,張志遠.鎂還原罐內(nèi)強化換熱研究[J].北京科技大學學報,2012,34(7):830-835.

    FENG Junxiao,ZHANG Zhiyuan.Enhancement of heat transfer in a magnesium reduction retort[J].Chinese Journal of Engineering,2012,34(7):830-835.

    (責任編輯:劉潔)

    猜你喜歡
    數(shù)值模擬溫度場
    鋁合金加筋板焊接溫度場和殘余應力數(shù)值模擬
    基于紋影法的溫度場分布測量方法
    測控技術(2018年4期)2018-11-25 09:47:10
    MJS工法與凍結法結合加固區(qū)溫度場研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    F40船板鋼高溫鍛壓溫度場及組織演變模擬
    上海金屬(2016年2期)2016-11-23 05:34:33
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設計與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    欧美区成人在线视频| 亚洲av熟女| 在线免费十八禁| 国产精品蜜桃在线观看| 欧美日韩国产亚洲二区| av福利片在线观看| 久久精品久久久久久噜噜老黄 | 欧美成人午夜免费资源| 国产乱来视频区| 国产精品伦人一区二区| 国产精品蜜桃在线观看| 久久久久国产网址| 波野结衣二区三区在线| 在线观看66精品国产| 嫩草影院入口| 国产一区二区在线av高清观看| 91狼人影院| 看片在线看免费视频| 精品国内亚洲2022精品成人| 国产美女午夜福利| 国产真实乱freesex| 在现免费观看毛片| 一个人观看的视频www高清免费观看| 国产精品日韩av在线免费观看| 亚洲伊人久久精品综合 | 亚洲欧美成人综合另类久久久 | 国产一级毛片在线| 免费在线观看成人毛片| 欧美丝袜亚洲另类| 男人和女人高潮做爰伦理| 春色校园在线视频观看| 国产 一区 欧美 日韩| 久久精品国产亚洲av涩爱| 免费av不卡在线播放| 日韩国内少妇激情av| 99热网站在线观看| 国产视频首页在线观看| 亚洲精品影视一区二区三区av| 一本久久精品| 欧美色视频一区免费| 欧美丝袜亚洲另类| 久热久热在线精品观看| 国产综合懂色| 久久久久免费精品人妻一区二区| 欧美区成人在线视频| 中文欧美无线码| 亚洲综合色惰| 日本午夜av视频| 亚洲精品日韩av片在线观看| 欧美日韩在线观看h| 久久精品久久久久久噜噜老黄 | 非洲黑人性xxxx精品又粗又长| 在现免费观看毛片| 国产片特级美女逼逼视频| 91久久精品国产一区二区成人| 亚洲内射少妇av| 亚洲最大成人av| 免费观看性生交大片5| 少妇的逼好多水| 亚洲av免费在线观看| 国产真实伦视频高清在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲欧美精品专区久久| 日韩一本色道免费dvd| 熟女电影av网| 中文资源天堂在线| 尾随美女入室| 国产 一区 欧美 日韩| 欧美高清性xxxxhd video| 国产综合懂色| 色5月婷婷丁香| 国产一级毛片在线| 久久亚洲精品不卡| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美精品自产自拍| 亚洲精品日韩在线中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 国产av不卡久久| 亚洲婷婷狠狠爱综合网| 最近2019中文字幕mv第一页| 又爽又黄a免费视频| 国产免费男女视频| 亚洲中文字幕日韩| 少妇丰满av| 国产亚洲午夜精品一区二区久久 | 赤兔流量卡办理| 小说图片视频综合网站| 日本免费一区二区三区高清不卡| 日韩精品有码人妻一区| 亚洲欧洲日产国产| 一边摸一边抽搐一进一小说| 国产精品久久久久久久久免| 国产伦精品一区二区三区四那| 人人妻人人看人人澡| 欧美区成人在线视频| 少妇被粗大猛烈的视频| 床上黄色一级片| 99久久精品国产国产毛片| 女的被弄到高潮叫床怎么办| 美女被艹到高潮喷水动态| 美女脱内裤让男人舔精品视频| 亚洲成人中文字幕在线播放| 亚洲精品色激情综合| 97热精品久久久久久| 国产精品久久久久久精品电影小说 | 一边摸一边抽搐一进一小说| 欧美成人a在线观看| 51国产日韩欧美| 日本三级黄在线观看| 日韩成人伦理影院| 精品一区二区三区人妻视频| 菩萨蛮人人尽说江南好唐韦庄 | 天堂av国产一区二区熟女人妻| 亚洲国产精品久久男人天堂| 啦啦啦观看免费观看视频高清| a级一级毛片免费在线观看| 国产精品99久久久久久久久| 中文亚洲av片在线观看爽| 只有这里有精品99| 国产私拍福利视频在线观看| 国产成人精品婷婷| 嫩草影院新地址| 久久久久久大精品| 欧美激情在线99| 日本黄色片子视频| 午夜精品一区二区三区免费看| 亚洲性久久影院| 男的添女的下面高潮视频| 久久久久久久久大av| 男女视频在线观看网站免费| 91精品国产九色| 欧美激情国产日韩精品一区| 欧美日韩一区二区视频在线观看视频在线 | 国产精品无大码| 丰满乱子伦码专区| 国产精品嫩草影院av在线观看| 国产亚洲一区二区精品| 中文字幕av在线有码专区| 熟女电影av网| 亚洲欧美成人综合另类久久久 | 亚洲欧美成人精品一区二区| 久久久久久九九精品二区国产| 日本wwww免费看| 六月丁香七月| 亚洲国产精品国产精品| 免费av不卡在线播放| 男人舔女人下体高潮全视频| 一区二区三区四区激情视频| 一级av片app| 久久精品国产亚洲av涩爱| 国产精品综合久久久久久久免费| av免费在线看不卡| 国产精品久久久久久av不卡| 一卡2卡三卡四卡精品乱码亚洲| 久久久久国产网址| av国产免费在线观看| 精品免费久久久久久久清纯| 久99久视频精品免费| 夜夜爽夜夜爽视频| 亚洲成人av在线免费| 亚洲乱码一区二区免费版| 亚洲av免费在线观看| 18禁在线播放成人免费| 国产美女午夜福利| 看黄色毛片网站| 男的添女的下面高潮视频| 国产精品,欧美在线| 日日摸夜夜添夜夜爱| 亚洲欧美一区二区三区国产| 欧美一级a爱片免费观看看| 亚洲电影在线观看av| 亚洲伊人久久精品综合 | 人妻系列 视频| 日日摸夜夜添夜夜添av毛片| 永久免费av网站大全| 美女国产视频在线观看| av视频在线观看入口| 久99久视频精品免费| 丝袜喷水一区| 久久久久免费精品人妻一区二区| 男人舔奶头视频| 久久久久久国产a免费观看| 美女cb高潮喷水在线观看| 亚洲精品456在线播放app| 小蜜桃在线观看免费完整版高清| 日本黄色视频三级网站网址| 日韩大片免费观看网站 | 人妻少妇偷人精品九色| 亚洲在线自拍视频| 在线观看66精品国产| www.色视频.com| videossex国产| 天堂网av新在线| 精品不卡国产一区二区三区| 国产一区二区在线观看日韩| or卡值多少钱| 少妇人妻一区二区三区视频| 亚洲伊人久久精品综合 | 少妇熟女欧美另类| 国产激情偷乱视频一区二区| 男的添女的下面高潮视频| 七月丁香在线播放| 国产免费视频播放在线视频 | 美女高潮的动态| 秋霞伦理黄片| h日本视频在线播放| 久久精品人妻少妇| 国产真实伦视频高清在线观看| 欧美bdsm另类| 嫩草影院精品99| 好男人视频免费观看在线| 非洲黑人性xxxx精品又粗又长| 嘟嘟电影网在线观看| 直男gayav资源| 精品不卡国产一区二区三区| 免费播放大片免费观看视频在线观看 | 国产一区二区在线观看日韩| 成人特级av手机在线观看| 成年女人永久免费观看视频| 一区二区三区免费毛片| 又爽又黄无遮挡网站| 亚洲欧美精品综合久久99| 青春草视频在线免费观看| 欧美区成人在线视频| 久久久久久伊人网av| 毛片女人毛片| 日韩视频在线欧美| 小蜜桃在线观看免费完整版高清| 99在线视频只有这里精品首页| 美女被艹到高潮喷水动态| 自拍偷自拍亚洲精品老妇| 日日干狠狠操夜夜爽| 成人亚洲精品av一区二区| 亚洲色图av天堂| 久久人妻av系列| 欧美丝袜亚洲另类| 国产又色又爽无遮挡免| 少妇丰满av| 99在线视频只有这里精品首页| 成人一区二区视频在线观看| 男人舔奶头视频| av在线播放精品| 欧美成人精品欧美一级黄| 日韩在线高清观看一区二区三区| 日韩欧美精品v在线| 一二三四中文在线观看免费高清| 级片在线观看| 国产乱人视频| 淫秽高清视频在线观看| 啦啦啦啦在线视频资源| 麻豆成人av视频| 精品久久久久久久久亚洲| 深爱激情五月婷婷| 久久久久性生活片| 亚洲在线观看片| 国产毛片a区久久久久| 1000部很黄的大片| 又黄又爽又刺激的免费视频.| 老女人水多毛片| 熟女电影av网| 爱豆传媒免费全集在线观看| 欧美成人一区二区免费高清观看| 久热久热在线精品观看| 欧美高清成人免费视频www| 欧美激情久久久久久爽电影| 午夜福利在线观看免费完整高清在| 中文字幕av成人在线电影| 国产精品一及| 午夜亚洲福利在线播放| 乱系列少妇在线播放| 一卡2卡三卡四卡精品乱码亚洲| 国产在线男女| 国产熟女欧美一区二区| 午夜a级毛片| 91久久精品国产一区二区三区| 亚洲av免费高清在线观看| 国产精品永久免费网站| 白带黄色成豆腐渣| 色噜噜av男人的天堂激情| a级毛片免费高清观看在线播放| 国产精华一区二区三区| 天堂影院成人在线观看| 又粗又爽又猛毛片免费看| 欧美成人a在线观看| 国产伦在线观看视频一区| 亚洲精品色激情综合| 成年av动漫网址| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 91精品一卡2卡3卡4卡| 我的女老师完整版在线观看| 国产精品久久久久久久久免| 亚洲自拍偷在线| 国产精品熟女久久久久浪| 亚洲最大成人中文| 校园人妻丝袜中文字幕| 国产精品99久久久久久久久| 深爱激情五月婷婷| 在线播放无遮挡| 高清av免费在线| 级片在线观看| 美女国产视频在线观看| 高清毛片免费看| 两个人的视频大全免费| 男人和女人高潮做爰伦理| 波野结衣二区三区在线| av卡一久久| 亚洲在久久综合| 淫秽高清视频在线观看| 一级毛片aaaaaa免费看小| 日产精品乱码卡一卡2卡三| 男女视频在线观看网站免费| 亚洲国产日韩欧美精品在线观看| 我的女老师完整版在线观看| 午夜精品国产一区二区电影 | 亚洲av成人av| 国产麻豆成人av免费视频| 日韩,欧美,国产一区二区三区 | 久久久久久久久久成人| 偷拍熟女少妇极品色| 天堂av国产一区二区熟女人妻| 国产精品精品国产色婷婷| 99久久成人亚洲精品观看| 只有这里有精品99| 建设人人有责人人尽责人人享有的 | 成人毛片60女人毛片免费| 午夜福利在线观看吧| 亚洲在线自拍视频| 亚洲欧美日韩无卡精品| 亚洲欧美日韩卡通动漫| 精品99又大又爽又粗少妇毛片| 高清午夜精品一区二区三区| 亚洲av免费在线观看| 建设人人有责人人尽责人人享有的 | 国产成人a∨麻豆精品| 国产久久久一区二区三区| 一个人看视频在线观看www免费| 婷婷色麻豆天堂久久 | 国产精品麻豆人妻色哟哟久久 | 一个人看的www免费观看视频| 日韩欧美精品v在线| 99久久无色码亚洲精品果冻| 亚洲精品乱码久久久v下载方式| 伊人久久精品亚洲午夜| 午夜福利在线在线| 久久99热这里只频精品6学生 | 国产午夜精品论理片| 草草在线视频免费看| 99热6这里只有精品| 成年女人看的毛片在线观看| 一边摸一边抽搐一进一小说| 熟女电影av网| 久久久欧美国产精品| 边亲边吃奶的免费视频| 日韩三级伦理在线观看| 99久久精品国产国产毛片| av在线老鸭窝| 午夜福利高清视频| 国产高清国产精品国产三级 | 免费看av在线观看网站| 免费av观看视频| 精品久久久久久久久亚洲| 欧美成人a在线观看| 日本黄色片子视频| 99久久九九国产精品国产免费| 国产精品国产三级国产av玫瑰| 一本一本综合久久| 成人一区二区视频在线观看| 2022亚洲国产成人精品| 日本爱情动作片www.在线观看| 亚洲精品乱码久久久v下载方式| 欧美日韩综合久久久久久| 国产精品日韩av在线免费观看| 成人二区视频| 亚洲精品日韩av片在线观看| 日本一二三区视频观看| 人妻夜夜爽99麻豆av| 日韩高清综合在线| 国产熟女欧美一区二区| 亚洲四区av| 国产高清不卡午夜福利| 日本爱情动作片www.在线观看| 亚洲,欧美,日韩| 亚洲自偷自拍三级| 欧美潮喷喷水| 激情 狠狠 欧美| 欧美潮喷喷水| 欧美性猛交黑人性爽| 少妇的逼水好多| 91久久精品国产一区二区三区| 看片在线看免费视频| 一级毛片电影观看 | 国产成人91sexporn| 亚洲国产精品国产精品| 黑人高潮一二区| 日本三级黄在线观看| 精品久久久久久电影网 | 国产高清视频在线观看网站| 国产欧美日韩精品一区二区| 少妇被粗大猛烈的视频| 国产av一区在线观看免费| 精品国产三级普通话版| 午夜久久久久精精品| 99久久中文字幕三级久久日本| 亚洲天堂国产精品一区在线| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品乱码久久久v下载方式| 国内揄拍国产精品人妻在线| 亚洲天堂国产精品一区在线| 国产成人精品婷婷| 国产日韩欧美在线精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 观看美女的网站| 亚洲最大成人中文| 精品久久久久久久末码| 亚洲欧美精品自产自拍| 亚洲人与动物交配视频| 黄片无遮挡物在线观看| 干丝袜人妻中文字幕| 日韩欧美精品v在线| 国产精品蜜桃在线观看| 国产成人91sexporn| 真实男女啪啪啪动态图| 国产伦精品一区二区三区视频9| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 国产淫语在线视频| 免费看光身美女| 只有这里有精品99| 国产亚洲精品久久久com| 波多野结衣巨乳人妻| 高清视频免费观看一区二区 | 久久韩国三级中文字幕| 亚洲av中文字字幕乱码综合| 欧美精品国产亚洲| 成人鲁丝片一二三区免费| 精品酒店卫生间| 99在线人妻在线中文字幕| 国产单亲对白刺激| 日韩成人伦理影院| 免费看a级黄色片| 国产成人a∨麻豆精品| 看黄色毛片网站| 国产不卡一卡二| 亚洲国产精品成人综合色| av专区在线播放| 草草在线视频免费看| 久久99热6这里只有精品| 黄色一级大片看看| 99在线人妻在线中文字幕| 国产91av在线免费观看| 午夜福利高清视频| 97热精品久久久久久| 久久久久久久久久成人| 午夜视频国产福利| АⅤ资源中文在线天堂| 久久久久免费精品人妻一区二区| 免费观看性生交大片5| 国产精品一区二区在线观看99 | 老司机影院成人| 男的添女的下面高潮视频| 久久99热6这里只有精品| 亚洲成人中文字幕在线播放| 乱码一卡2卡4卡精品| av又黄又爽大尺度在线免费看 | 亚洲av男天堂| 99久久精品一区二区三区| 国产亚洲最大av| 国产精品乱码一区二三区的特点| 欧美xxxx黑人xx丫x性爽| 国产一级毛片七仙女欲春2| 乱系列少妇在线播放| 亚洲乱码一区二区免费版| 晚上一个人看的免费电影| 久久久久久久久大av| 蜜桃久久精品国产亚洲av| 亚洲乱码一区二区免费版| 久久久精品欧美日韩精品| 狠狠狠狠99中文字幕| 亚洲无线观看免费| 成人漫画全彩无遮挡| 成人综合一区亚洲| 亚洲av男天堂| 亚洲美女搞黄在线观看| 蜜臀久久99精品久久宅男| 韩国av在线不卡| 日韩av在线免费看完整版不卡| 最后的刺客免费高清国语| 国内精品一区二区在线观看| 色哟哟·www| 夜夜爽夜夜爽视频| 亚洲av男天堂| 日韩欧美国产在线观看| 久久这里只有精品中国| 午夜福利成人在线免费观看| 日韩av在线免费看完整版不卡| 精品一区二区三区视频在线| 你懂的网址亚洲精品在线观看 | 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看 | 在线观看美女被高潮喷水网站| 中文字幕久久专区| 国产成人精品婷婷| 精品久久久噜噜| 七月丁香在线播放| 观看免费一级毛片| 最近最新中文字幕免费大全7| 又爽又黄a免费视频| 日日撸夜夜添| 99热这里只有精品一区| 91aial.com中文字幕在线观看| 欧美3d第一页| 26uuu在线亚洲综合色| 久久精品国产亚洲av天美| 国产不卡一卡二| 蜜桃久久精品国产亚洲av| 又粗又硬又长又爽又黄的视频| 欧美潮喷喷水| 69av精品久久久久久| 一级黄色大片毛片| 一级黄片播放器| 亚洲自拍偷在线| 中文字幕免费在线视频6| 久久久久精品久久久久真实原创| 女的被弄到高潮叫床怎么办| 国产av不卡久久| 免费看光身美女| 日本欧美国产在线视频| 亚洲自拍偷在线| 国产色爽女视频免费观看| 99九九线精品视频在线观看视频| 精品国产一区二区三区久久久樱花 | 精品久久久久久久久久久久久| 男女边吃奶边做爰视频| 国产精品乱码一区二三区的特点| 免费观看的影片在线观看| 神马国产精品三级电影在线观看| 如何舔出高潮| 亚洲自偷自拍三级| 不卡视频在线观看欧美| 1024手机看黄色片| 国产午夜精品久久久久久一区二区三区| 欧美高清性xxxxhd video| 男女边吃奶边做爰视频| 国产精品1区2区在线观看.| 狂野欧美白嫩少妇大欣赏| 国产伦精品一区二区三区视频9| 国产伦精品一区二区三区四那| 精品免费久久久久久久清纯| 成人毛片60女人毛片免费| 特大巨黑吊av在线直播| 国产一区二区在线av高清观看| 人妻制服诱惑在线中文字幕| 国产精品嫩草影院av在线观看| 久久国产乱子免费精品| 国产探花在线观看一区二区| 一二三四中文在线观看免费高清| 国产精品综合久久久久久久免费| 欧美成人精品欧美一级黄| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av免费高清在线观看| 亚洲精品日韩av片在线观看| 中文字幕久久专区| 日韩av在线大香蕉| 日韩av在线免费看完整版不卡| 直男gayav资源| 国产精品爽爽va在线观看网站| 三级男女做爰猛烈吃奶摸视频| 国产片特级美女逼逼视频| 国产精品无大码| 欧美97在线视频| 一个人免费在线观看电影| 三级国产精品欧美在线观看| 欧美另类亚洲清纯唯美| 久久久久久大精品| 精品国产一区二区三区久久久樱花 | av专区在线播放| 嫩草影院入口| 人人妻人人澡欧美一区二区| av黄色大香蕉| 免费观看人在逋| 成人毛片60女人毛片免费| 免费黄网站久久成人精品| 国产视频首页在线观看| 成年版毛片免费区| 午夜老司机福利剧场| 国产极品精品免费视频能看的| 国内少妇人妻偷人精品xxx网站| 简卡轻食公司| 日日啪夜夜撸| 99热6这里只有精品| 日本与韩国留学比较| 女人被狂操c到高潮| 久久99热这里只有精品18| 久久人妻av系列| 国产国拍精品亚洲av在线观看| 久久这里有精品视频免费| 亚洲国产欧美在线一区| 免费人成在线观看视频色| 听说在线观看完整版免费高清| 少妇熟女aⅴ在线视频| 亚洲婷婷狠狠爱综合网| 麻豆久久精品国产亚洲av| 99热这里只有是精品在线观看| 日韩中字成人| 99热这里只有是精品50| 狂野欧美激情性xxxx在线观看| 色5月婷婷丁香| 欧美成人a在线观看| 我的老师免费观看完整版| 久久精品夜色国产| 亚洲欧美精品专区久久| 亚洲美女搞黄在线观看| 高清av免费在线| 久久久久久久久久成人|