?? wafer2d.m
字號:
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; % inRTobj = rooftemp('thomasswann');RTobj = newton(RTobj);Troof = gettemp(RTobj)unpack(RTobj,'param') % to get kH2, Lrxr, TsusLrxr = Lrxr*1000; % convert to mmkH2 = kH2/1000; % J/(sec*mm*K)ESiC = 0.91; % unitlesssig = 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/Kthickness = 350; % wafer thickness, micronsLw = thickness/10^6; % convert to malph = 0.75; % 1-transmittivityRad = 2*2.54/100; % wafer radiuskwaf = kSapp;Esat = ESiC;Eroof = 0.4;% Number of quadrature pointsnr = 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 guessT = scalarfield(S,rand(nr,nz)+Tsus);y = scalarfield(R,0);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% set up an instance of the twodww objectvar = 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);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -