pfrspacetime.m

来自「《Matlab在化学工程中的应用》该书是我见过的一本相当有用的一本书 具有很高」· M 代码 · 共 38 行

function PFRSpaceTime
% PFRSpaceTime.m
% 等温管式反应器的管长计算
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/21 $

clear all
clc

k = 0.01;
T = 488.15;     % K
yA0 = 0.5;      % 进料气中A的摩尔分率
P = 5;          % atm
R = 0.082;      % 理想气体常数
CA0 = yA0*P/(R*T);  % mol/liter
epsilon = 1;
xAf = 0.8;

% the interval: [a, b], the step: d
a = 0;
b = xAf;
d = (b-1)/100; 

% 自适应Lobatto求积法(Adaptive Lobatto quadrature, high order)
I = quadl(@func,a,b,[],[],CA0,k,epsilon);
tau = CA0*I;
fprintf('\tSpace time: %.1f (s)',tau)

% ------------------------------------------------------------------
function f = func(xA,CA0,k,epsilon)  
CA = CA0*(1-xA)./(1+epsilon*xA);
rate = k*sqrt(CA);
f = 1./rate;