function main
clear all
clc
global n G dl M M1 b1 b2 C4 C6 C7 C8 C9 a
T=273.2+112.13; %反应器内平均温度,K
RMI=1.1; %氯化氢乙炔分子比
u=0.062*3.5; %气体线速度,m/s
C_A=exp(1.1835+1.055e-3*T)*4.184e3; %乙炔热容,J/kg/K
C_H=0.192*4.184e3; %氯化氢热容,J/Kg/K
C_V=exp(2.238+1.256e-3*T)*4.184*1000/62.5; %氯乙烯热容,J/Kg/K
L_A=exp(9.6084+1.413*log(T))*0.1; %乙炔导热系数,W/m
L_H=exp(-7.574+0.9769*log(T))*0.1; %氯化氢导热系数,W/m/K
L_V=0.018; %氯乙烯导热系数,W/m/K
r_A=0.0115*(T/353.2)^0.97*0.001; %乙炔黏度,Pa.s
r_H=0.0165*(T/353.2)^0.95*0.001; %氯化氢黏度,pa.s
r_V=8.601e-3*(T/250.1)^0.