radcoolblock.m

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

function B = radcoolblock

% radcoolblock.m A nonlinear algebraic equation model object constructor, derived
%                from the naemodel class. Called as
%
%                A = radcoolblock
%
% INPUT PARAMETERS
%                none

% Number of quadrature points

nx = 20;
ny = 10;

S = quadgrid('slab',nx,'x',[0 1])*quadgrid('slab',ny,'y',[0 0.3]);

Dx  = linearoperator(S,'d', 'x');
DDx = linearoperator(S,'dd','x');
Dy  = linearoperator(S,'d', 'y');
DDy = linearoperator(S,'dd','y');

% Solution initial guess

T = scalarfield(S,rand(nx,ny));

n = valsize(T);

% BC value at y = 1

Tbc = scalarfield(S,1);

plot(T), title('Current solution')
drawnow

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

var   = assocarray({'T' T});
param = assocarray({ 'Dx' Dx 'DDx' DDx 'Dy' Dy 'DDy' DDy ...
                     'Tbc' Tbc });

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

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