wafer2d.m

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

function B = wafer2d

% wafer2d.m A nonlinear algebraic equation model object constructor,
%          derived from the naemodel class. Called as
%
%          A = wafer2d
%

%         thickness = 350; % microns
%         diameter = 2;    % in

RTobj = rooftemp('thomasswann');
RTobj = newton(RTobj);
Troof = gettemp(RTobj)

unpack(RTobj,'param') % to get kH2, Lrxr, Tsus

Lrxr = Lrxr*1000; % convert to mm
kH2 = kH2/1000;   % J/(sec*mm*K)

ESiC = 0.91;         % unitless
sig = 5.67e-8;       % J/(sec*m^2*K^4)
kSiC = 300;          % J/(sec*m*K) (approx)
kSapp = 10;          % J/(sec*m*K) (8 at 600C; 9 at 700C; about 11 at 800C)
betaSapp = 5.6e-6;   % expansion coeff  1/K
thickness = 350;     % wafer thickness, microns
Lw = thickness/10^6; % convert to m
alph = 0.75;         % 1-transmittivity
Rad = 2*2.54/100;    % wafer radius

kwaf = kSapp;

Esat  = ESiC;
Eroof = 0.4;

% Number of quadrature points

nr = 20;
nz = 10;

R = quadgrid('cyln',nr,'r',[0 Rad]);
S = R*quadgrid('slab',nz,'z',[0 Lw]);

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

% Solution initial guess

T = scalarfield(S,rand(nr,nz)+Tsus);
y = scalarfield(R,0);

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

var   = assocarray({'T' T 'y' y});
param = assocarray({ 'Dr' Dr 'DDr' DDr 'Dz' Dz 'DDz' DDz ...
                     'Troof' Troof 'Tsus' Tsus 'sig' sig ...
                     'alph' alph 'Ewaf' ESiC ...
                     'Esat' Esat 'Eroof' Eroof ...
                     'kwaf' kwaf 'kH2' kH2 'Lrxr' Lrxr ...
                     'Lw' Lw 'betaSapp' betaSapp ...
                 });

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

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