residual.m

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

function A = residual(A)

% residual.m Nonlinear equation model object residual update.
%
% INPUT/OUTPUT PARAMETERS
%        A : an object derived from the neqmodel class
%
% Copyright (c) by Raymond A. Adomaitis, 1998-2003

unpack(A,'var')
unpack(A,'param')

% Interior residual

RT = kwaf*DDr*T + kwaf*DDz*T;

% Boundary condition residuals

% symmetry at center

RT  = setbval(RT,Dr*T,'r','min');

% radiation from edge

RT  = setbval(RT, kwaf*Dr*T ...
    - Esat*Ewaf*alph*sig*(Tsus^4-T^4),'r','max');

% radiation from bottom

RT  = setbval(RT, -kwaf*Dz*T ...
    - Esat*Ewaf*alph*sig*(Tsus^4-T^4),'z','min');

% radiation/conduction to reactor top

RT = setbval(RT, -kwaf*Dz*T ...
    - kH2*(T-Troof)/Lrxr ...
    - Eroof*Ewaf*alph*sig*(T^4-Troof^4),'z','max');

% Wafer bow model

delT = getbval(T,'z','min') - getbval(T,'z','max');
Ry = DDr*y - betaSapp*delT/Lw;

Ry = setbval(Ry,Dr*y,'r','min');
Ry = setbval(Ry,y,'r','max');

A = set(A,'resid',{RT, Ry});