pfrconv.m

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

function PFRconv
% 已知等温管式反应器的管长,计算出口转化率
%
%   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
L = 5;                  % 反应管长，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 = fzero(@Func4Fzero,0.3,[],d,L,F0,yA0);         % 转化率
fprintf('所求的转化率为：\txA = %.4f',xA)

% ------------------------------------------------------------------
function f = Func4Fzero(xA,d,L,F0,yA0)  
I = quadl(@Func,0,xA);      % 用自适应Lobatto法求积,积分限为[0,xA]
f = L - 4*yA0*F0/(pi*d^2)*I;

% ------------------------------------------------------------------
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);