wafertemp.m

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

function B = wafertemp(N,Rw)

% wafertemp.m A nonlinear algebraic equation model object constructor, derived
%             from the naemodel class. Called as
%
%             A = wafertemp(N,Rw)
%
% INPUT PARAMETERS
%
%             N : number of quadrature points
%            Rw : wafer radius, (m)

% wafer parameters

delz  = 0.004;     % wafer thickness (m)
b     = 5.677e-8;  % Boltzmann constant
emiss = 0.7;       % wafer emissivity

k = 269-0.585*400+3.75e-4*400^2; % Si thermal conductivity, W/m/K

% quadgrid info

S = quadgrid('cyln',N,'r',[0 Rw]);

Dr  = linearoperator(S,'d', 'r');
DDr = linearoperator(S,'dd','r');

% Solution initial guess

T = scalarfield(S,300);

% Default radiation flux profile

Q = scalarfield(S,5000);

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

var   = assocarray({'T' T});
param = assocarray({ 'Dr' Dr 'DDr' DDr 'Q' Q ...
                     'delz' delz 'b' b 'emiss' emiss 'k' k});

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

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