北京化工大学2006一一2007学年第一学期 《计算化学》期末考试试卷标准答案 第一部分填空题（共30分，每空一分） 1.数值计算、非数值计算。 2.建立数学模型。 3.二分法：①函数x)在区间[a,b]内连续，②函数在a、b两点之值a)与b) 异号：牛顿迭代法：①∫”x)有明确表示式且容易计算②容易确定根的初值xo: 二分法：压力对化学反应中反应物平衡转化率的影响。 4.线性方程组的系数矩阵为对角占优势。 5.插值法。 6.一阶导数，由c~1数据求化学反应速率。 7.y=+h+e;y=6+hx+b3++bxn+s。最小二乘法。 8.x)函数形式己知，但其积分不能表示成初等函数的闭合形式，x)函数形 式未知，但其离散数据表已给出。(1)用插值程序求任意点的函数值。(2)用 Simpson求积程序计算[a,b]区间中离散点下的面积。 9.Euler法、Runge-.Kutta法。化学动力学。 10.量子化学/结构化学。 11,不用求导，甚至没有目标函数表达式时也可使用。 I2.Monte Carlo法、分子动力学方法。 13.Chemwindow/Chem Office Origin,Gaussian,Materials studio. 14.稳态模拟、动态模拟
1 北京化工大学 2006——2007 学年第一学期 《计算化学》期末考试试卷标准答案 第一部分 填空题（共 30 分，每空一分） 1．数值计算、非数值计算。 2．建立数学模型。 3．二分法：①函数 f(x)在区间[a,b]内连续，②函数在 a、b 两点之值 f(a)与 f(b) 异号；牛顿迭代法：①f ’(x)有明确表示式且容易计算②容易确定根的初值 x0； 二分法：压力对化学反应中反应物平衡转化率的影响。 4．线性方程组的系数矩阵为对角占优势。 5．插值法。 6．一阶导数，由 c~t 数据求化学反应速率。 7．y=ax+b+ε； = + + + ⋅⋅⋅+ + ε m m y b b x b x b x 0 1 1 2 2 。最小二乘法。 8．f(x)函数形式已知，但其积分不能表示成初等函数的闭合形式，f(x)函数形 式未知，但其离散数据表已给出。（1）用插值程序求任意点的函数值。（2）用 Simpson 求积程序计算[a,b]区间中离散点下的面积。 9．Euler 法、Runge-Kutta 法。化学动力学。 10．量子化学／结构化学。 11．不用求导，甚至没有目标函数表达式时也可使用。 12．Monte Carlo 法、分子动力学方法。 13．Chemwindow/Chem Office ，Origin，Gaussian，Materials studio。 14．稳态模拟、动态模拟
第二部分上机题（共70分） A 1.计算原理（化学原理和计算方法）(15分) 带式板应图装型东：业=5式大业 式中，Fo为反应物A的进气流速，”A为以浓度表示的反应速率。 反应管长为 L=4F 1 dx 其中反应物A的进气流速FA0=xAoF%=1/1+0.5)F=2F/3(5分) kAL-A2 2AKA-1 XA=0 no 0 n。/2 XA no(1-xA)noxA/2 no/2 可得P4=20-2p 3-XA pa 故，-=p-P (1) TAe =TAp /RT
2 第二部分 上机题（共 70 分） A 1． 计算原理（化学原理和计算方法）（15 分） 等温管式反应器模型为： A 0 A, R A0 A d 1 V x r F x c ∫ = 式中， FA0 为反应物 A 的进气流速， c rA, 为以浓度表示的反应速率。 反应管长为 A 0 A, 2 A0 A d 4 1 x D r F L x c ∫ = π 其中反应物 A 的进气流速 FA0 = xA0F0 = 1/(1+ 0.5)F0 = 2F0 / 3 （5 分） 2A A2 kA,1 kA,-1 0 xA = n0 0 n0 / 2 Ax (1 ) 0 A n − x n0 xA / 2 n0 / 2 可得 p x x p A A A 3 - 2(1− ) = p x x p A A A 3 - 2 = 故 p x x p k x x k k p k p t p r p A A A,-1 2 2 A 2 A A,1 A,-1 A 2 A,1 A A A, (3 - ) 3 - 4(1 ) d d 2 − − = = − = − （1） rA,c = rA, p / RT I
所以反应管长L πD20r (2) πD2 (3)(5分 由于被积函数 L有确定的形式（式(1）的倒数)，所以可直接利用Si即som求积 分的方法计算可得到(3)式S的值，其中调用自编子程序计算函数， 1的值，然后代入(2) 式，可得不同的转化率时反应管长。因积分下限0，为了防止出错，计算时取0.0001。(5 分) 程序框图(20分) 开始 输入：数据点数F0,EPS,P,T,XA0,XA①，(I=1,N)(5分) A=0.000L,B=XA①，(=1，N) (5分) 调用Simpson法子程序计算积分S (用自编子程序计算函数值)，L=4F,RTS/πD2 (5分) 输出：不同的转化率时反应管长 (5分) 结束 3. 源程序(20分) PROGRAM MAIN (10分) DOUBLE PRECISION EPS.P.XA.A.B.S2.L.D.FO.FAO.XAO DIMENSION XA(5) EXTERNALF OPEN(6.FILE-2007TESTAOUT TXU STATUS-UNKNOWN)
3 所以反应管长 A 0 A, 2 A0 A 0 A, 2 A0 A A d 4 1 d 4 1 x D r F RT x D r F L x p x c ∫ ∫ = = π π （2） 令 A 0 A, A d 1 x r S x p ∫ = （3） （5 分） 由于被积函数 p rA, 1 有确定的形式（式（1）的倒数），所以可直接利用 Simpson 求积 分的方法计算可得到（3）式 S 的值，其中调用自编子程序计算函数 p rA, 1 的值，然后代入（2） 式，可得不同的转化率时反应管长。因积分下限 0，为了防止出错，计算时取 0.0001。（5 分） 2． 程序框图（20 分） 3． 源程序（20 分） PROGRAM MAIN （10 分） DOUBLE PRECISION EPS,P,XA,A,B,S2,L,D,F0,FA0,XA0 DIMENSION XA(5) EXTERNAL F OPEN(6,FILE='2007TESTAOUT.Txt',STATUS='UNKNOWN') 开始 输入：数据点数 F0,EPS，P，T，XA0，XA (I),（I=1，N） （5 分） A=0.0001, B= XA (I),（I=1，N） （5 分） 调用 Simpson 法子程序计算积分 S （用自编子程序计算函数值）， 2 A0 L = 4F RTS /πD （5 分） 输出：不同的转化率时反应管长 L （5 分） 结束
EPS-10E-6 ATA XA0.10.020,0.30,0.40.0.50/ 0) 00 FAO XA0 R=8.314 D010=l.5 A=0.0001 R-YACD A B EPS S2.F) S2FA043.14159265D/DRT WRITE(6 I XAOL 11 FORMAT(X.XA-.F12.4.4X.L-F12.4) 10 CONTINUE END FUNCTION FCXA) (5分) DOUBLE PRECISION XA.K1.K2.P D-1 0E5 DATA K1/8565E-5/K2/6.833/ F=1K1*4*(1-XA)*(1-XA)产p*p-K2*XA*3-XA)P(3-XA)*2) RETURN END SUBROUTINE SIMP(A.B.EPS.S2.F) (5分) DOUBLE PRECISION A.B.H.EPS.S2.A1 H=B-A T1=0.5H*(F(A+F(B)) N=l 91=0 10 8=0 D020K=1,N A1=A+K-0.5*H S=S+E(A1 T2=0.5*(T1+H*S) S2=T2+T2.T1V30 IF(ABS(S2).LE.1.0)D=ABS(S2-S1) If(ABS(S2)GT LONTHEN D=ABS((S2-S1VS2) END IE IF(D.LT.EPS)GOTO50 NEN+N H=0.5H T1=T2 S1=S2 GOTO 10 ￥
4 EPS=1.0E-6 DATA XA/0.10,0.20,0.30,0.40,0.50/ D=0.10 P=1.0E5 F0=10.0 XA0=1/(1+0.5) FA0=F0*XA0 R=8.314 T=911.0 DO 10 I=1,5 A=0.0001 B=XA(I) IF(B.EQ.0)STOP CALL SIMP(A,B,EPS,S2,F) L=S2*FA0*4/3.14159265/D/D*R*T WRITE(6,11) XA(I),L 11 FORMAT(3X,'XA=',F12.4,4X,'L=',F12.4) 10 CONTINUE END FUNCTION F(XA) （5 分） DOUBLE PRECISION XA,K1,K2,P P=1.0E5 DATA K1/8.565E-5/,K2/6.833/ F=1/((K1*4*(1-XA)*(1-XA)*P*P-K2*XA*(3-XA)*P)/((3-XA)**2)) RETURN END SUBROUTINE SIMP(A,B,EPS,S2,F) （5 分） DOUBLE PRECISION A,B,H,EPS,S2,A1 H=B-A T1=0.5*H*(F(A)+F(B)) N=1 S1=0 10 S=0 DO 20 K=1,N A1=A+(K-0.5)*H 20 S=S+F(A1) T2=0.5*(T1+H*S) S2=T2+(T2-T1)/3.0 IF(ABS(S2).LE.1.0)D=ABS(S2-S1) IF(ABS(S2).GT.1.0)THEN D=ABS((S2-S1)/S2) END IF IF(D.LT.EPS)GOTO 50 N=N+N H=0.5*H T1=T2 S1=S2 GOTO 10
50 RETURN END 运行结果(15分) 1.8765(每个3分) m 42 4000 13.5401 5000 58596 x≤0.4 当x>0.4以后，所需反应器管长剧增。因此设计反应器时应选取 ◇ 1.计算原理（化学原理和计算方法）(15分) 设反应达平衡时生成C0的物质的量为xmol,有 CH4+H2O CO+3H2 起始物质的量/mol 12 00 平衡物质的量mol 1-x2-x x 3x 反应达平衡时系统中总物质的量为：1-x+2-x+x+3x=3+2x 故平衡时系统中各物质的组成为： 1-x 2-x yal=3+2x' yH0=3+2x x 3x yc0=3+2x m=3+2x 式中：y为物质的量分数。 (5分) 代入平衡关系式得： K”=KP x(3x)3「 =1.247 (1) (1-x)2-x)(3+2x)p° (1)式经整理后得 x.(3x)3 -1.247=0 (2) 1-x)2-x)(3+2x)p° 令：f(x)= x-(3x)3 p -1.247 (5分) (1-x)2-x)(3+2x)p° 调用二分法求解方程f(x)=0
5 50 RETURN END 4． 运行结果（15 分） xA= .1000 L/m= 1.8765（每个 3 分） xA= .2000 L/m= 4.2871 xA= .3000 L/m= 7.6705 xA= .4000 L/m= 13.5401 xA= .5000 L/m= 58.5966 由计算结果可知：当 0.4 xA > 以后，所需反应器管长剧增。因此设计反应器时应选取 0.4 xA ≤ 。 B 1．计算原理（化学原理和计算方法）（15 分） 设反应达平衡时生成 CO 的物质的量为 x mol，有 CH4+H2O CO+3H2 起始物质的量/mol 1 2 0 0 平衡物质的量/mol 1-x 2-x x 3x 反应达平衡时系统中总物质的量为： 1− x + 2 − x + x + 3x = 3 + 2x 故平衡时系统中各物质的组成为： x x y x x y x x y x x y 3 2 3 , 3 2 3 2 2 , 3 2 1 2 4 2 CO H CH H O + = + = + − = + − = 式中：y 为物质的量分数。 （5 分） 代入平衡关系式得： 1.247 (1 )(2 ) (3 2 ) (3 ) ( ) 2 3 = ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − − + ⋅ = = ∑ θ θ θ x p p x x x x n p p K K B n （1） （1） 式经整理后得 1.247 0 (1 )(2 ) (3 2 ) (3 ) 2 3 − = ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − − + ⋅ θ x p p x x x x （2） 令： 1.247 (1 )(2 ) (3 2 ) (3 ) ( ) 2 3 − ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − − + ⋅ = θ x p p x x x x f x （5 分） 调用二分法求解方程 f (x) = 0