pfrlength.m

matlab在化学工程中的应用实例程序,对大家研究学习有指导作用。

function PFRlength
% 等温管式反应器的管长计算
%
%   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

global PA0 eps

d = 0.1;                % 反应器内经，m
F0 = 10;                % 进气速率，kmol/h
P = 1;                  % 反应压力，atm
yA0 = 2/3;              % A的初始摩尔分率
PA0 = yA0 * P;          % A的初始分压，atm
delta = -1/2;
eps = yA0 * delta;      % eps: 膨胀因子
xA = [0.1:0.1:0.5];     % A的摩尔转化率

% 计算不同转化率下的管长
for i = 1:length(xA)
    I = quadl(@func,0,xA(i));   % 用自适应Lobatto法求积,积分限为[0,xA]
    L(i) = 4*yA0*F0/(pi*d^2)*I; % L: 管长,m
end
xAe = fzero(@Rate,0.5);         % 平衡转化率

% 结果显示
fprintf('Results:\n\txA:\t')
fprintf('\t%.2f',xA)
fprintf('\n\tL (m):')
fprintf('\t%.2f',L)
fprintf('\n\txAe = %.3f',xAe)

% ------------------------------------------------------------------
function f = func(xA)  
f = 1./Rate(xA);

% ------------------------------------------------------------------
function r = Rate(xA)
global PA0 eps
PA = PA0 * (1-xA)./(1+eps*xA);
PS = PA0/2 * xA./(1+eps*xA);
r = 116 * (PA.*PA - PS/1.27);