Hello !
I'm trying to implement a Phase Change Material in PDE Toolbox in a 2D model, so that I build a model and I need to use Temp-dependent heat capacity and conductivity. I define the capacity as follow and it works well
d = @(region,state) 900.*(4250 +(15000-4250).*exp(-((21.7-state.u)./3).^2)); %Cp*rho
...specifyCoefficients(modele,'m',0,'d',d,...,'face',1);
But to define conductivity I need a stattement condition on temperature. So I wrote a function in c.m
function y = c(region,state)if state.u < 21.7 y = 0.25;else y = 0.20;end
And I call the function in the main program as follow ..
cond = @cspecifyCoefficients(modele,'m',0,'d',d,'c',cond,...,'face',1);
But when I run the program , I get the following error:
Error using formGlobalKF2DCoefficient evaluation function, "c", was requested to calculate coefficients at 276 locations so should have returned amatrix with 276 columns. Instead it returned a matrix with 1 columns.
Does someone know how to implement a statement condition on my coefficients ?
Best Answer