数据处理与数值分析论.docx

数据处理与数值计算结课论文学 院 班 级 年 级 姓 名 指导教师 2014年12月29日1数据处理与数值计算结课论文摘 要化学反应工程主要研究反应动力学和反应器，反应动力学是反应器设计与分析的基础。工业反应器中除了化学反应外，还伴随有反应组分的传递、反应热的吸收或排放、反应物料的流动与混合等物理过程，极其错综复杂，其研究方法以数学模拟为主。现在在解决化学反应工程实际问题的过程中常常遇到对数学模型求解和学学模型参数估值问题。常规方法难以或不能对其进行处理，通过数据处理与数值计算这门课程的学习，我基本掌握了插值，数值微分，数值积分，最小二乘法，（非）线性方程（组）解法，常微分方程（组）的初值和边值问题。并且对计算机软件在化工生产问题中的应用有所了解，现结合具体实例介绍如下。利用最小二乘法确定动力学方程环氧乙烷可杀灭细菌（及其内孢子）、霉菌及真菌，因此可用于消毒一些不能耐受高温消毒的物品。主要用于制造其他各种溶剂(如溶纤剂等),稀释剂,非离子型表面活性剂,合成洗涤剂、抗冻剂、消毒剂、增韧剂和增塑剂等。与纤维素发生羟乙基化可合成得水溶性树脂(其环氧乙烷含量约75%)。还可用作熏蒸剂、涂料增稠剂、乳化剂、胶黏剂和纸张上浆剂等。乙烯的深度氧化是乙烯环氧化生产环氧乙烷的工业生产中重要的副反应。对反应选择性有较大影响。在铂催化剂上，乙烯深度氧化的动力学方程可表示为r=kpApB/(1+KBpB)2式中pA、pB分别表示乙烯及氧的分压。在473K等温下的实验数据如表1:表1 实验数据 PA*103MPaPB*103MPar*104mol/g*min18.990 3.230 0.672 214.220 3.000 1.072 38.860 4.080 0.598 48.320 2.030 0.713 54.370 0.890 0.610 67.750 1.740 0.834 77.750 1.820 0.828 86.170 1.730 0.656 96.130 1.730 0.694 106.980 1.560 0.791 112.870 1.060 0.418 试求该温度下的反应速率常数k和吸附平衡常数KB。【1】本题计划用最小二乘原理解决，最小二乘法的原理暨通过使数据的残差平方和最小寻找数据的最佳函数匹配。但是因为此模型公式非线性，不能直接套用线性最小二乘法，需要对公式进行线性化处理：Kr=kPAPB(1+KBPB)2 PAPBr=1/k+KBPB/k直线方程形式为y=b+a*x。即有：y=PAPBrx=PBa=KB/kb=1k利用excel计算结果如下表2：表2 最小二乘法数据表xi*103yixi*yi*103xi2*1051 3.230 0.657 2.123 1.0433 2 3.000 0.631 1.892 0.9000 3 4.080 0.777 3.172 1.6646 4 2.030 0.487 0.988 0.4121 5 0.890 0.253 0.225 0.0792 6 1.740 0.402 0.700 0.3028 7 1.820 0.413 0.751 0.3312 8 1.730 0.403 0.698 0.2993 9 1.730 0.391 0.676 0.2993 10 1.560 0.371 0.579 0.2434 11 1.060 0.270 0.286 0.1124 SUM22.870 5.055 12.090 5.6875 b=i=111xi*yi- 111*（I=111xi）*（i=111yi）/i=111xi2-111*(i=111xI)2a=y-bx由表中数据有：b=169.511482a=0.107098726k=1/b2KB=a*k有：k=87.18291638KB=1582.759093乙烯深度氧化动力学方程可以用公式:r=87.183PAPB（1+1582.759PB)2表示对回归方程进行分析：如下表3表3 原始数据与回归数据误差PA*103MPaPB*103MPar*104mol/g*minr*104mol/g*min*104mol/g*min18.993.230.6720.67761506900.005615068969 214.2231.0721.12558237610.053582376119 38.864.080.5980.5666578595-0.031342140461 48.322.030.7130.82959876410.116598764077 54.370.890.610.5844587582-0.025541241806 67.751.740.8340.83424588030.000245880301 77.751.820.8280.8165868004-0.011413199621 86.171.730.6560.66595405450.009954054523 96.131.730.6940.6616366863-0.032363313740 106.981.560.7910.7888186693-0.002181330664 112.871.060.4180.3699029320-0.048097067953 由表中数据可得：残差平方和Q=i=111(ri-ri)2=2.172*10-10剩余标准差 S=111-2*(ri-ri)2=4.913*10-6相关指数 R2=1-(ri-ri)2(ri-ri)2=0.9221利用数值积分进行反应器设计自从1928年Diels和Alder报道环戊二烯与顺丁烯二酸酐的环加成反应后，Diels-Alder反应几十年来一直吸引着有机化学家们的广泛兴趣。这个反应为合成六元环化合物提供了一条简单的途径。不仅产率高，而且反应的立体专一性和定位选择性强，成为有机合成中一个十分重要的反应。在药物、天然化合物、萜类化合物的合成中得到了广泛应用。在活塞流反应器中绝热进行丁二烯和乙烯合成环已烯反应C4H6 + C2H4C6H10( A) ( B) ( R)该反应为气相反应，反应速率方程为：r = kCACB，k= 3.16107 exp(13840 /T)，L/( mols)进料为丁二烯与乙烯的等摩尔混合物，温度为440。操作压力1.013105 Pa。该反应的热效应Hr =1.256 105kJ/mol。假定各气体的热容为常数，CpA=154，CpB=85.6，CpR =249，单位为J /( molK)。要求丁二烯的转化率达12%，试计算空时。【1】此过程为绝热变温变容反应空时：=VrQ0=CA0*0XAdXAk*CA*CB=1CA0*0XA1+A*yA0*XA2*dXAk1-XA2*(T0T)2出口温度：T=T0+*XA已知：yA0=0.5，CA0=CB0，A=1-21=-1选取入口温度440为基准温度，DHr=-1.256*105kJ/mol*K是在440下的热效应，取1mol物质A为基准，有：A + B R 1 1 0 21-XA 1-XA XA 2-XA所以，当XA=0.12时，有yA=yB=1-0.122-0.12=0.468 yR=0.12/(2-0.12)=0.064 Cpt=154*0.468+85.6*0.468+229*0.064=128.1J/mol*K=yA0*-HrCpt1+A*yA0*XA=0.5*1.256*105128.1*1-0.5*0.12=521.5CA0=PA0RT=1.013*105*0.58.314*103*713=8.5544*10-3kmol/m3将数据带入得： =18.544*10-3*00.121-0.5XA2*dXA3.16*107*e13840713+521.5XA*1-XA2*713713+521.5*XA2上式可用数值积分求得数值解。数值分析中，数值积分是计算定积分数值的方法和理论。在数学分析中，给定函数的定积分的计算不总是可行的。许多定积分不能用已知的积分公式得到精确值。数值积分是利用黎曼积分等数学定义，用数值逼近的方法近似计算给定的定积分值。利用复化辛普森公式：Sn=h3*fa+fb+2i=1n-1fx2i+4i=1nfx2i-1f(x)= 1-0.5x23.16*107*e13840713+521.5x*1-x2*713713+521.5x2 0,0.12设步长为h=0.01，则区间0,0.12等分为12份，n=b-a2*h=6。excel计算结果如表4：表4 复化辛普森公式数值积分xif（xi）Sn0.00 8.4840618.4840610.5267360.01 7.54306815.086140.02 6.72067226.882690.03 6.00040412.000810.04 5.36827621.47310.05 4.8123779.6247530.06 4.32254317.290170.07 3.8900847.7801680.08 3.50755214.030210.09 3.1685536.3371060.10 2.86758411.470330.11 2.5998995.1997980.12 2.3614022.361402=1cA0*Sn=0.5267360.008544=61.6471s较文献参考值69.5735s偏低，可能是数值积分步长过大，也可能是积分方法存在误差，可以采用精度更高的公式方法：例如以分段四次插值为基础的复化柯特斯公式：Cn=h907fx0+32k=1nfxk-34+