rooftemp.m

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

function B = rooftemp(reactortype)

% rooftemp.m A nonlinear algebraic equation model object constructor,
%          derived from the naemodel class. Called as
%
%          A = rooftemp(reactortype)
%
% INPUT PARAMETERS
%
%         reactortype % 'aixtron','thomasswann'

switch reactortype
    
    case 'aixtron'
        disp('aixtron geom')
        Tsus = 1200+273; % K
        kRoof = 0.06; % W/(m K) - value for carbon foam
                      % which is approximated by the kAir
        Lrxr = 2.54/100; % susceptor/roof spacing m
        Lins = 2.54/100; % insulation thickness m
        Esus = 0.91; % SiC emissivity
        Eroof = Esus;
    case 'thomasswann'
        disp('thomas swann geom')
        Tsus = 1040+273; % K
        kRoof = 25; % W/(m K) - value for stainless steel
        Lrxr = 1/100; % wafer/showerhead spacing m
        Lins = 5/100; % showerhead effective thickness m
        Esus = 0.91; % SiC emissivity
        Eroof = 0.2; % stainless steel emissivity
end

sig = 5.67e-8;      % J/(sec*m^2*K^4)
kH2 = 0.185;        % W/(m K)

% Solution initial guess

T = (Tsus+300)/2;

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

var   = assocarray({'T' T});
param = assocarray({ 'Tsus' Tsus 'kRoof' kRoof 'Lrxr' Lrxr ...
    'Lins' Lins 'Esus' Esus 'Eroof' Eroof 'sig' sig 'kH2' kH2 });

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

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