正交配置求解问题.docVIP

正交配置求解问题： 运用正交配置法求解有轴向扩散的固定床反应器中催化反应的温度和浓度分布。柱形固体床反应器中催化反应的温度和浓度方程为： ＝+ =+ ︱r=1 =|r=0=0 -|r=1=Biw[T(1,z)-Tw(z)], -|r=1=0 T(r,0)=T0 , c (r,0)=c0 ＝+ =+ ︱r=1 =|r=0=0 -|r=1=Biw[T(1,z)-Tw(z)], -|r=1=0 T(r,0)=T0 , c (r,0)=c0 其中R（c，T）为催化反应的速率方程，其形式为 R(c,T)= + =+ ︱r=1 =|r=0=0 -|r=1=Biw[T(1,z)-Tw(z)], -|r=1=0 T(r,0)=T0 , c (r,0)=c0 其中R（c，T）为催化反应的速率方程，其形式为 R(c,T)= 解题思路： 应用对称的正交配置法，有下面的方程和初始条件： =+ (1- ) = + (1-) Tj(0)=T0 , cj(0)=c0 边界条件为：-AN+1,iTi=Biw(TN+1-Tw) , AN+1,ici=0 将温度和浓度的边界条件代入微分方程，消去边界值，可得2N个常微分方程，而将两边界条件的代数方程同2N个常微分方程组联合，就组成2N+2个微分代数方程组。结合正交配置系数的计算程序与常微分方程组或微分方程组求解程序，可得到反应器中的温度和浓度分布。具体做法如下： 利用对称的正交配置格式： 1、对称常微分方程程序：（collab.for,Dlsode.for） 主程序： IMPLICIT REAL*8 (A-H,O-Z) EXTERNAL FEX, JEX DIMENSION AS(19,19),BS(19,19),Q(19,19),XS(19),WS(19) DIMENSION DIF1(19),DIF2(19),DIF3(19),ROOT(19),V1(19),V2(19) DIMENSION Y(99), ATOL(99), RWORK(10920), IWORK(120) DOUBLE PRECISION YN1,YN2 COMMON /AB/ N,AS,BS COMMON /BC/YN1,YN2 C N--- FOR SYMMETRIC COLLOCATION USED FOR PARTICLE AND C M--- FOR ASYMMETRIC COLLOCATION USED FOR COLUMN N=7 IW=1 IS=2 CALL COLL(AS,BS,Q,XS,WS,19,N,IW,IS) NS=N+1 WRITE(*,*) * SYMMETRIC SITUATION: * WRITE(*,*) * POLYNOMIAL ROOTS * WRITE(*,*) (XS(I),I=1,NS) WRITE(*,*) WRITE(*,*) * A-MATRIX * DO 20 I=1,NS 20 WRITE(*,*) (AS(I,J),J=1,NS) WRITE(*,*) WRITE(*,*) * B-MATRIX * DO 30 I=1,NS 30 WRITE(*,*) (BS(I,J),J=1,NS) WRITE(*,*) WRITE(*,*) * W-MATRIX * WRITE(*,*) (WS(J),J=1,NS) C CALCULATING THE PARAMETERS OF THE PROBLEM, WHICH WILL BE USED C FOR THE DIMENSIONLESS FORM OF AND DEFINING OF THE PROBLEM. NEQ=2*N LRW=22+9*NEQ+NEQ**2 LIW=20+NEQ C INITIAL CONDITIONS DO 201 I=1,N Y(I) = 1.D0 Y(N+I)=0.D0 201 CONTINUE YN1=1.0D0 YN2=0.D0