1、課 程 設 計 說 明 書 武 漢 工 程 大 學 化工(hugng)與制藥學院課程設計說明書課題(kt)名稱 專業(zhuny)班級 學生學號 學生姓名 學生成績 指導教師 課題工作時間 武漢工程大學化工與制藥學院課程設計任務書專業(zhuny) 班級(bnj) 學生(xu sheng)姓名 發題時間： 年 月 日課題名稱1萬噸/年產鄰二甲苯催化氧化制鄰苯二甲酸酐反應器的設計課題條件（文獻資料、儀器設備、指導力量）在過量空氣存在情況下，以五氧化二釩為催化劑對鄰二甲苯氧化制備鄰苯二甲酸酐產品，反應氣中的苯酐產品通過兩級冷卻得以回收，第一級得到液體產品，第二級得到固體產品，最后，對粗品采用精餾的

2、方法進行精制。反應器為列管式固定床反應器，管內裝球型催化劑，管外采用熔鹽循環以移走放應熱，并副產高壓水蒸氣。反應產物以氣體的形式從反應器移出，進入氣體冷卻器，用高壓水冷卻，并副產中壓蒸汽。粗品采用精餾的方法進行精制得到產品鄰苯二甲酸酐。鄰二甲苯的轉化率為100，轉化為苯酐的選擇性為79.生產時間按每年8000h計算。設計任務（含實驗、分析、計算、繪圖、論述等內容）使用迭代方法來確定滿足所有關于鄰二甲苯轉化率和溫度分布要求的最優設計。由計算機模擬得到的設計參數包括：完全反應所需的管長，為有效控制溫度所需的管徑，能將溫差減小到可接受程度的冷卻劑流量以及為有效散熱和防止反應器內出現熱點所需的總傳熱系

3、數。（設計參數包括：反應管直徑長度、催化劑裝填高度、殼程換熱面積、壓力降、冷卻劑流量、沿管程溫度分布、沿管程組成分布、傳熱速率、反應物及產物和冷卻劑的流速。要求：編寫鄰苯二甲酸酐催化反應器模擬計算的FORTRAN程序，并在計算機上計算出上述各項設計參數；用CAD畫出鄰二甲苯氧化制鄰苯二甲酸酐的工藝流程圖及催化反應器的外型圖；設計說明書邏輯清楚，層次分明，書寫工整，獨立完成；附簡短的中、英文摘要。設計所需技術參數序號 名稱 反應器進料組成（） 反應器出口氣體組成（）1 鄰二甲苯 9.5 2 氧氣 19.0 5.33 氮氣 71.5 71.54 苯酐 10.55 馬來酸酐 0.66 CO2 6.6

4、7 水 5.58 雜質 0.19 反應原料進口溫度205反應溫度307反應器出口溫度307熔鹽組成40NaNO2,7%NaNO3 ,53%KNO3(Wt%)熔鹽入口溫度371熔鹽出口溫度385反應壓力175kPa設計(shj)說明書內容設計(shj)任務書題目(tm)名稱，目錄及頁碼前言包括設計依據、主要內容、特點等按設計任務順序說明結語包括設計體會、收獲、評述、建議、致謝等參考文獻主要技術符號說明進度計劃設計動員，下達任務 0.5天收集資料，閱讀教材，理順設計思路 0.51.5天編程，調試，設計計算 35天整理設計資料，撰寫設計說明書 2天指導教師審查，答辯 2天指導教師簽名： 教研室主任簽

5、名： 年 月 日 年 月 日化工與制藥學院課程設計綜合成績評定表學生姓名學生班級設計題目 鄰二甲苯催化氧化制鄰苯二甲酸酐反應器的設計指導教師評語指導教師簽字：年 月 日答辯記錄答辯組成員簽字： 記錄人：年 月 日成績綜合評定欄設計情況答辯情況項 目權重分值項 目權重分值1、計算和繪圖能力351、回答問題能力202、綜合運用專業知識能力102、表述能力（邏輯性、條理性）103、運用計算機能力和外語能力104、查閱資料、運用工具書的能力55、獨立完成設計能力56、書寫情況（文字能力、整潔度）5綜合成績指導(zhdo)教師簽名： 學科(xuk)部主任簽名： 年 月 日 年 月 日摘 要化學(huxu

6、)反應器是整個化工(hugng)工藝流程的核心(hxn)設備，是化工生產中實現化學物質轉化的必要工序，被譽為化工廠的“心臟”。本次課程設計是綜合運用已學課程知識，完成以一個單元設備反應器設計為主的一次實踐性教學，是連接理論和實際的重要橋梁，在整個教學過程中起到培養學生實際能力的重要作用。使用迭代方法來確定滿足所有關于鄰二甲苯轉化率和溫度分布要求的最優設計。運用計算機模擬并計算得到管長、管徑、壓力降、冷卻劑流量等設計參數。反應器設計要求學生更加熟悉化工工程設計的基本內容。宗旨是使學生能夠熟練掌握化工生產工藝設計理念、設計程序和設計方法，鍛煉學生分析問題、解決問題的能力。培養學生綜合理論知識解決實

7、際問題的能力，訓練學生的計算能力和思考問題的能力，訓練學生查閱技術資料、選用公式、搜集數據和制圖的能力。關鍵詞：催化反應器；設計；鄰二甲苯AbstractChemical reactor is the core of the entire chemical process equipment, chemical production processes is necessary to achieve the transformation of chemical substances, chemical plant known as the heart. The course has been

8、designed to be integrated use of knowledge courses, complete with a unit device once the reactor design based practical teaching, is to connect theory and practice important bridge, play an important practical ability of students in the teaching process effect. Using an iterative method to determine

9、 satisfy all about the o-xylene conversion rate and temperature distribution optimum design requirements. The use of computer simulation and calculate the tube length, diameter, pressure drop, coolant flow and other design parameters. Reactor design requires students to become more familiar with the

10、 basic content of chemical engineering. Aim is to enable students to master the chemical production process design, design process and design methods, training students to analyze problems, problem-solving abilities. Students comprehensive theoretical knowledge and ability to solve practical problem

11、s, training of students of computing power and the ability to think, and training students ability to access technical data, the choice of formulas, data collection and mapping.Keywords： HYPERLINK /dict_result.aspx?searchword=%e5%82%ac%e5%8c%96%e5%8f%8d%e5%ba%94%e5%99%a8&tjType=sentence&style=&t=cat

12、alytic+reactors catalytic reactor； HYPERLINK /dict_result.aspx?searchword=%e8%ae%be%e8%ae%a1&tjType=sentence&style=&t=design design； HYPERLINK javascript:showjdsw(jd_t,j_) o-xylene 目 錄 TOC o 1-3 h z u HYPERLINK l _Toc409361909 摘 要 PAGEREF _Toc409361909 h I HYPERLINK l _Toc409361910 Abstract PAGEREF

13、_Toc409361910 h II HYPERLINK l _Toc409361911 前 言 PAGEREF _Toc409361911 h 1 HYPERLINK l _Toc409361912 第1章 文獻綜述 PAGEREF _Toc409361912 h 2 HYPERLINK l _Toc409361913 1.1設計內容 PAGEREF _Toc409361913 h 2 HYPERLINK l _Toc409361914 1.2催化反應器的分類 PAGEREF _Toc409361914 h 2 HYPERLINK l _Toc409361915 1.3設計原理 PAGERE

14、F _Toc409361915 h 2 HYPERLINK l _Toc409361916 第2章反應器設計參數的計算 PAGEREF _Toc409361916 h 4 HYPERLINK l _Toc409361918 2.1物質的理基礎數據 PAGEREF _Toc409361918 h 4 HYPERLINK l _Toc409361919 2.2物料衡算 PAGEREF _Toc409361919 h 4 HYPERLINK l _Toc409361920 2.3計算機模擬 PAGEREF _Toc409361920 h 5 HYPERLINK l _Toc409361921 2.4

15、管數 PAGEREF _Toc409361921 h 6 HYPERLINK l _Toc409361922 2.5殼程換熱面積 PAGEREF _Toc409361922 h 6 HYPERLINK l _Toc409361923 2.6壓力降 PAGEREF _Toc409361923 h 7 HYPERLINK l _Toc409361933 2.7.1管程壓力降 PAGEREF _Toc409361933 h 7 HYPERLINK l _Toc409361934 2.7.2殼程壓力降 PAGEREF _Toc409361934 h 7 HYPERLINK l _Toc40936193

16、5 2.7冷卻劑流量 PAGEREF _Toc409361935 h 8 HYPERLINK l _Toc409361936 2.8沿管程溫度分布 PAGEREF _Toc409361936 h 9 HYPERLINK l _Toc409361937 2.9沿管程組成分布 PAGEREF _Toc409361937 h 10 HYPERLINK l _Toc409361938 2.10傳熱速率 PAGEREF _Toc409361938 h 11 HYPERLINK l _Toc409361939 第3章 設計結果匯總 PAGEREF _Toc409361939 h 12 HYPERLINK

17、l _Toc409361940 第4章 結論與建議 PAGEREF _Toc409361940 h 14 HYPERLINK l _Toc409361941 致 謝 PAGEREF _Toc409361941 h 15 HYPERLINK l _Toc409361942 參考文獻 PAGEREF _Toc409361942 h 16 HYPERLINK l _Toc409361943 附 錄 PAGEREF _Toc409361943 h 17武漢工程大學本科課程設計設計前 言化學(huxu)反應器是整個化工(hugng)工藝流程的核心(hxn)設備，是化工生產中實現化學物質轉化的必要工序，被

18、譽為化工廠的“心臟”。本次課程設計是綜合運用化學反應工程 化工過程模擬 化工工藝學 化工設計及開發和化工制圖等課程知識，完成以一個單元設備反應器設計為主的一次性實踐教學，是連接理論和實際的重要橋梁，在整個教學過程中起到培養學生實際能力的重要作用。本次設計使用迭代方法來確定滿足所有關于鄰二甲苯轉化率和溫度分布要求的最優設計。由計算機模擬得到的設計參數包括：完全反應所需的管長，為有效控制溫度所需的管徑，能將溫差減小到可接受程度的冷卻劑流量以及為有效散熱和防止反應器內出現熱點所需的總傳熱系數。（設計參數包括：反應管直徑長度、催化劑裝填高度、殼程換熱面積、壓力降、冷卻劑流量、沿管程溫度分布、沿管程組成

19、分布、傳熱速率、反應物及產物和冷卻劑的流速。）反應器課程設計要求學生更加熟悉化工工程設計的基本內容。宗旨是使學生能夠熟練掌握化工生產工藝設計理念、設計程序和設計方法，鍛煉學生分析問題、解決問題的能力。培養學生綜合理論知識解決實際問題的能力，訓練學生的計算能力和思考問題的能力，訓練學生查閱技術資料、選用公式、搜集數據和制圖的能力。因自己學識有限，設計中一定有很多疏漏和錯誤之處，懇請老師同學批評指正。在此，感謝老師的耐心指導使得我的設計工作得以圓滿完成！第1章 文獻(wnxin)綜述設計(shj)內容化學(huxu)反應器是整個化工工藝流程的核心，是實現化學物質轉化的必要工序。反應器的設計一般包括

20、下列內容：根據反應過程的化學基礎和生產工藝的基本要求，進行反應器的選型設計；根據化學反應與有關流體力學、熱量、質量傳遞過程綜合宏觀反應動力學，計算反應器的結構尺寸；反應器的機械設計、穩定性分析。催化反應器的分類催化反應器是一個在其中發生催化反應的容器或罐，可由各種材料制成，材料的選擇取決于反應發生時的溫度和壓力以及所要求的抗腐蝕性，催化反應器也可能有多種類型，這取決于所用催化劑的種類以及反應物和產物的物理特性。催化反應器可按各種方法分類， 根據催化作用的性質分類，可分為均相催化反應器和非均相催化反應器。如果催化劑并不形成與反應物或產物分離的單獨相，則催化作用稱為均相的。如果催化劑形成與反應物分

21、離的單獨相，則催化作用稱為非均相的。根據在反應器中的催化劑顆粒是否移動，非均相反應器可分為：固定床反應器，顆粒固定在位置上稱為密相固定床；移動床反應器，顆粒堆積在緩慢移動的密相床中，排出某些結垢的顆粒并加入某些新鮮的顆粒；流化床反應器，催化劑顆粒被向上流動的氣體托住如流化床一樣；懸浮床反應器，催化劑顆粒懸浮在液體中。 設計原理在過量空氣存在情況下，以五氧化二釩為催化劑對鄰二甲苯氧化制備鄰苯二甲酸酐產品，反應氣中的苯酐產品通過兩級冷卻得以回收，第一級得到液體產品，第二級得到固體產品，最后，對粗品采用精餾的方法進行精制。反應器為列管式固定床反應器，管內裝球型催化劑，管外采用熔鹽循環以移走放應熱，并

22、副產高壓水蒸氣。反應產物以氣體的形式從反應器移出，進入氣體冷卻器，用高壓水冷卻，并副產中壓蒸汽。粗品采用精餾的方法進行精制得到產品鄰苯二甲酸酐。鄰二甲苯的轉化率為100，轉化為苯酐的選擇性為79.生產時間按每年8000h計算。反應的方程式為：主反應：副反應： 使用迭代方法來確定滿足所有關于鄰二甲苯轉化率和溫度分布要求的最優設計。由計算機模擬得到的設計參數包括：完全反應所需的管長，為有效控制溫度所需的管徑，能將溫差(wnch)減小到可接受程度的冷卻劑流量以及為有效散熱和防止反應器內出現熱點所需的總傳熱系數。第2章反應器設計(shj)參數的計算(j sun)物質(wzh)的理基礎數據由chemBl

23、ink化學品數據庫查得反應物與生成物的基本信息與理化性質如表2.1所示。表2.1物質的理化性質名稱鄰二甲苯苯酐馬來酸酐英文名o-XylenePhthalic anhydrideMaleic Anhydride分子結構分子式C8H10C8H4O3C4H2O3分子量106.16148.1298.06密度0.8791.50.1g/cm3熔點-26-23 C129-132 C54C沸點143-145 C284 C202.2222C閃點31 C152 C103.3333C溶解性6 g/L (20 C)物料衡算根據設計任務所給定生產要求，已知鄰二甲苯的轉化率為100，轉化為苯酐的選擇性為79，生產時間按每

24、年8000h計算。進行物料衡算，計算出反應器出口氣體的總流量，以及進出口各組分的質量流量，為計算反應器的設計參數提供依據。如表2.2所示。出口氣體的總流量為：110780000.105=11904.76kg/h表2.2反應器進出口物料組成及流量名稱反應器進口組成/%進口質量流率kg/h反應器出口組成/%出口質量流率kg/h鄰二甲苯9.501130.9524O219.002261.90485.30630.9524N271.508511.904871.508511.9048苯酐10.501250.0000馬來酸酐0.6071.4286CO26.60785.7143水5.50654.7619雜質0.

25、1011.9048計算機模擬(mn)計算機模擬(mn)，使用迭代方法來確定(qudng)滿足所有關于鄰二甲苯轉化率和溫度分布要求的最優設計。圖2.1溫度、轉化率與l的關系圖關系由計算機模擬結果，為有效控制溫度所需的管徑取為214mm.管數反應器所需要的列管數（正方形排列，管心距PT=32mm）為：N=110780000.34=3677根所以最終排管密度為一根管0.032m2=977根管/m2最小管束面積為367797799%=3.81m2由公式(gngsh)D24=(NPT)2得最小管束(gunsh)直徑為D=436770.03223.1499%=2.75m3m殼程換熱面積(min j)使用K

26、ern(1950)的方法，用雷諾數和普朗特數來確定殼程換熱面積：A=DiCBPT其中A：殼程錯流面積，m2；Di：殼內徑，m，Di =3m；C：列管間距，m， C=30-21=9mm=0.009mB：折流板間距，m，一般取殼徑的（0.21）倍，取B=30.2=0.6mPT：管心距，m，取32mm，即PT=0.032m因此，殼程換熱面積為A=DiCBPT=3(30-21)10-30.60.032=0.51m2壓力降管程壓力降通過催化劑的壓降，由填充床的方程計算：PL=1501-DP+1.75G1- GDPg其中為填充床層空隙率，=0.405；為氣相粘度，=3.3410-5Pa.S；DP為顆粒直徑

27、，DP=6.0mm；G為質量流量，G=2.60kg/(m2.s)； 為氣相平均密度，=0.77 kg/m3 L為床層高度，L=4m所以P4=1501-0.4053.3410-5610-3+1.752.601-0.405 0.4052.60610-30.779.81P=1.505KPa管程壓降有三部分：催化劑床層、絲網和支撐以及(yj)出入口損失。催化劑層壓降：P0=1.0757KPa絲網(s wn)和支撐：P1=0.3 KPa出入口損失(snsh)：P2=0.1 KPa管程總壓降P=1.505+0.3+0.1=1.905KPa殼程壓力降Kern給出了計算管式換熱器殼程壓降的關聯式：P=fGs2

28、Ds(N+1)2gDe其中f 為摩擦因數；Gs為熔鹽質量流量；Ds為殼層內徑，Ds=3m； N折流板數，N=6； 為熔鹽密度，=1851kg/m3； Ds為有效管徑，Ds=0.02079 ；為管程數，=1當Re0500時適用f0=5Re0-0.228則摩擦因數f=512687.67-0.228=0.5799所以P=0.57991020.1423(5+1)218519.80.020791=16802.54Pa=16.802kPa冷卻劑流量由公式UA0T=msCpsTs1-Ts2設U初值為70W/(m2.), T=T1-T2=307-205=102tm=t1-t2lnt1t2=371-205-38

29、5-307ln371-205385-307=116.51A0=HrUtm=1399510370116.51=1715.51m2m=UA0TcpsTs1-Ts2=701715.511021.61103(385-371)=520.27kg/sG=mA=520.270.51=1020.14kg/(m2.s)u=G=1020.141851=0.55m/sPr(377)=6.29Re=du=Gd=1020.140.024751.9910-3=12687.67查Re與jH的關系(gun x)線圖，得到(d do) jH=42因此(ync)Nu=jHPr0.33=426.290.33=77.06 De=0.

30、99d0=0.990.021=0.02079mmh0=NukDe=77.060.5150.03168=1252.64w/(m2.)hk=2kdln(d0di)=2250.032+0.0252ln（3225）=7106.80w/(m2.)1U=1hi+1h0+1hk=174.85+11252.64+17106.80則總傳熱系數為U =69.93 w/(m2)本次迭代以初值U的值相差為70-69.9370100%=0.1%所以相差小，迭代結束。因此熔鹽流量為520.27kg/s沿管程溫度分布Calderbank推導了一段反應管區間內的能量平衡： T-Tw=-e-l-e-l其中=BSM-G； =Hr

31、y0cp；=2URGcp；T：反應氣體溫度，；Tw：管壁溫度，；y0：二甲苯摩爾分率；Hr：鄰二甲苯的摩爾放熱量，Hr=13395KJ/mol；Cp：反應物的熱容，cp=1.19KJ/(kg.)R：列管半徑U-總傳熱系數，U70W/(m2.)。考慮了循環導熱鹽后，在一段反應管區間內的熱平衡如下T1TrmcpdT+my0Hrx2-x1=T1T2mcpdT+UATmy0Hrx2-x1=mcpT2-Tr+UAT其中(qzhng)m為質量(zhling)流量；T為列管平均溫度 T=0.8Tgas+0.2Twall； A為換熱面積(min j)T=T-Tsalt需要迭代確定。用同一區間的能量平衡可以計算

32、出鹽的溫度變化UAT=mscpsTs1-Ts2其中ms為熔鹽質量流量，kg/s; cps為導熱鹽熱容cps=1.61KJ/(kg.) =BSM-G=15000.836.329359.24.65 =Hry0cp=133950.0284.651.19=1465.57 =2URGcp=2700.01059359.21.191.20 R=0.0212=0.0105沿管程組成分布假定反應具有一級反應動力學特征，速率常數與溫度無關dxdl=-k1-x x=1-exp(-L)=BSM-G其中B：催化劑床層的裝填密度， B=1500kg/m3S：指前因子，由Calderbank(1974)實驗測定，值為0.8

33、M-：進料的平均分子質量；G：表觀氣體流量kg(m2.h)M-=0.905MWair+0.095MWox=0.90529+ 0.095106=36.32G=3.24kg/h40.0252m=9359.2kg(m2.h)=2.6kg(m2.s)=BSM-G=15000.836.329359.24.65傳熱(chun r)速率管程換熱系數(xsh)Kreith和Bohn(1986)給出了基于(jy)努塞爾特數（Nusselt）的管程換熱系數關聯式：Nu=0.203(RePr)0.33+(0.220Re)0.8Pr0.4其中 Re=GD(1-) ； Pr=cpk ； hi=NukD ；則Re=GD(

34、1-) =2.600.0213.3410-7（1-0.405）=2.75105Pr=cpk=1.193.3410-70.0499=7.9710-6Nu=0.203(RePr)0.33+(0.220Re)0.8Pr0.4=0.203（2.751057.9710-6）0.33+0.220（2.75105）0.8（7.9710-6）0.4=0.2031.2956+0.220205.139645.39hi=NukD=45.390.04990.021=107.86W/(m2.)第3章 設計結果(ji gu)匯總設計(shj)數據(shj)匯總表如下所示。表3.1設計結果匯總表項目內容反應器反應器類型管式

35、固定床催化反應器方向垂直操作方式連續進料反應物100%鄰二甲苯，空氣空氣與烴的比率9.5:1烴負荷340g/(h.管)總流量39523.80（kg/h）溫度205壓力175kpa催化劑類型五氧化二釩V2O5鄰二甲苯轉化率100%苯酐的選擇性79%苯酐產率1.10kg/kg（鄰二甲苯）載體6mm瓷球支撐方式金屬絲網和夾環列管數目3677長度4.0m催化劑填充高度3.9m內徑21mm外徑25mm規格BWG12(高規格)傳熱面積0.51m2程數單程，向上列管排列方式正方形直排列管間距7mm壓降2.221kpa入口溫度205出口溫度370殼程流體導熱鹽流體組成（質量百分數）40%NaNO2,7%NaN

36、O3,53%KNO2流量215057.23(kg/h)入口溫度371出口溫度385每管移出的熱量2783.3w傳熱系數82.8697w/(內徑5000mm程數單程折流板3壓降2.221kpa第4章 結論(jiln)與建議這次(zh c)課程設計讓我再一次感受(gnshu)到再完美的想象也比不上一次嘗試，理論知識再扎實，也只有通過實際的操作和驗證才能正真領悟理論知識與工程實際的聯系與區別。剛拿到設計任務書的時候，連續看了幾天的書，卻無從下手。沒有完全理清楚設計思路，就開始嘗試用Aspen Plus模擬，每次模擬到中間的時候就感覺缺少條件進行不下去，換了各種方法模擬，還是不行。經過反思之后，我又重

37、新回到設計任務書。我就把設計任務書重新看了幾遍，然后就去查閱文獻，找相關的設計方法，流程模擬的教程。所以我就按照設計任務書上的進度計劃，先進行設計計算。在設計方案的過程中，通過查閱資料，最大的感觸就是書本上學習到的知識只是冰山一角，在實際生產過程中有許多問題不是僅靠書本知識就能解決的。總體的設計雖然比較繁瑣，但還是學到了很多東西。做的比較匆忙，可能存在一些不足之處，之后會繼續改進的。致 謝最后(zuhu)，感謝(gnxi)學校(xuxio)開設了課程設計這門課程，讓我有了一次能夠學以致用的機會，此次設計不僅鞏固了所學的化學反應工程知識更極大拓寬了我的知識面，讓我認識了實際化工生產過程和理論的聯

38、系和差別，這對將來的畢業設計無疑將起到重要的作用。但因自己學識有限，設計中一定有很多疏漏和錯誤之處，懇請老師同學批評指正。最難忘的是得到了老師的熱心指導使得我的設計工作得以圓滿完成，*老師耐心的輔導和幫助使我這次的課程設計事半功倍！在此，向他們表示衷心的感謝！參考文獻吳元欣，丁一剛，劉生鵬. 化學(huxu)反應(fnyng)工程M. 北京(bi jn)：化學工業出版社，2010.6吳元欣，張珩. 反應工程簡明教程M. 北京：高等教育出版社，2013.田文德，汪梅，王英龍. 化工過程計算機輔助設計基礎M. 北京：化學工業出版社，2012.8孫蘭義. 化工流程模擬實訓：Aspen Plus教程M

39、. 北京：化學工業出版社，2012.熊杰明, 楊索和. HYPERLINK 10/opac/openlink.php?title=Aspen+Plus%E5%AE%9E%E4%BE%8B%E6%95%99%E7%A8%8B Aspen Plus實例教程M. 北京：化學工業出版社，2013.尹先清. 化工設計M. 北京：石油工業出版社，2006.11鄒蘭，閻傳智. 化工工藝工程設計M. 成都：成都科技大學出版社，1998.8米鎮濤. 化學工藝學M. 北京：化學工業出版社，2006.3譚榮偉. 化工設計CAD繪圖快速入門M. 北京：化學工業出版社:，2014.7附 錄附錄(fl)1計算機程序模擬(

40、mn)： FORTRAN清單(qngdn)Program to model temperature and composition profilesIn catalytic tubular fixed bed reactor for the partialOxidation of ortho-xylene to phthalic anbydrideVariable names refer to equations derived in Appendix B.DOUBLE PRECISION X(0:90),T(0:90),TS(0:90),TW(0:90),ENGTH(0:90),FLOW,S

41、_FLOW,AREA,CP,S_CP,DEL_T(2),TR1,R2,DENS,MW,G,Y0,RADIUS,TOTLEN,EFFLEN,ALPHA,TA,GAMMA,FACTOR,S,HC,U,HEAT_R.WTEGER I,KDefine output file for results.PEN(FILE=“REACTOR.DAT”,UNIT=3.5,STATUS=“OLD”)Physical properties commom to both catalysts.DNS=1500DIUS=0.0125EA=0.003927AT_R=12868HC=210/1000CP=1.19S_CP=1

42、.61Physical properties specific to the LAR catalysts.FLOW=0.00100G=9813.862Y0=0.0952MW=30.9Physical properties specific to trational catalystsFLOW=0.00139G=10170Y0=0.0356MW=29.8Define constantsALPHA=DENS*S*MW/GBETA=HEAT_R*Y0*ALPHA/CPGAMMA=2*3600*HC/(RADTUS*G*CP)PRINT*,LAR CATALYSTPRINT*,COUNTER-CURRENT SALT FLOWPRINT*,GRADUATED CATALYST ACTIVITYPRINT*Parameters to be inputted by the user.PRINT*,LENGTH OF TUBULAR R