radialrxrmodel.m

MDPSAS工具箱是马里兰大学开发的

function B = radialrxrmodel(A,vo,Ts)

% radialrxrmodel.m A nonlinear equation model object constructor, derived
%           from the neqmodel class. Called as
%
%          B = radialrxrmodel(A)
%
% INPUT PARAMETERS
%
%         A : a radialrxrgeom object

unpack(A)

if nargin < 2
    vo = 20;   % inlet velocity (m/s)
end
if nargin < 3
    Ts = 1000; % susceptor temperature (K)
end

xAo = 0.1;
k1  = 1e19;
E1  = 3e4;

xBo = 0;
k2  = 1e15;
E2  = 2e4;

xCo = 0;
kd  = 100;

Tamb = 300; % ambient temperature, K

h = 1; % heat transfer coefficient

T  = scalarfield(R,(Tamb+Ts)/2); % Temperature solution initial guess
xA = scalarfield(R,xAo);         % Species A profile solution initial guess
xB = scalarfield(R,xAo/2);       % Species B profile solution initial guess
xC = scalarfield(R,xAo/2);       % Species C profile solution initial guess

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% set up an instance of the rxnmodel object

var   = assocarray({ 'T' T 'xA' xA 'xB' xB 'xC' xC });
param = assocarray({ 'R' R 'RSF' RSF 'nrp' nrp ...
    'Rs' Rs 'Rp' Rp 'Rw' Rw 'Ro' Ro ...
    'vo' vo 'Tamb' Tamb 'Ts' Ts 'h' h ...
    'xAo' xAo 'k1' k1 'E1' E1 ...
    'xBo' xBo 'k2' k2 'E2' E2 ...
    'xCo' xCo 'kd' kd  ...
    });

A = naemodel(var,param);
B = struct([]);

B = class(B,'radialrxrmodel',A);