I'm trying to obtain the temperature at a given point of my layered 2D cross-section geometry considering thermal conductivity, convection and radiation. Inside the geometry there are horizontal and vertical edges.
The code is working fine when no convection or radiation are used. However when I include them, the results don't vary much (only on the 6th decimal) although they should do it (for radiation) since I'm considering high temperatures.
I'm defining what edges should consider the radiation and convection in the following way:
eps = @(region,state) obj.epsilon*region
outerCC = @(region,state) obj.conv_coef*region
Do I need to specify the pde coefficients for each edge or face? or the solver extracts them when using thermalBC?
In the code, I'm not using commands such as ApplyBoundaryCondition or specifyCoefficients . If I try to use specifyCoefficients, I define the coefficients in the following way:
Ta = 25;
c = length*K;
a = @(~,state) 2*obj.conv_coef + 2*obj.epsilon*obj.tm_SS.StefanBoltzmannConstant*state.u.^3;
f = 2*obj.conv_coef*Ta + 2*obj.epsilon*obj.tm_SS.StefanBoltzmannConstant*Ta^4;
d = length*Density*SpecificHeat;
but always obtain the following error:
Check for missing argument or incorrect argument data type in call to function 'specifyCoefficients'.
The model is defined as follows:
obj.tm_SS = createpde('thermal','steadystate');
gdm = [r0;r9;r1;r2;r3;r4;r5;r6;r7;r8]';
obj.g = decsg(gdm,'R0+(R9-R1)+R1+(R2-R3-R6-R7)+R3+R6+R7+(R4-R6-R7)+R6+R7+(R5-R8)+R8',['R0';'R9';'R1';'R2';'R3';'R4';'R5';'R6';'R7';'R8']');
Do you know what I'm doing wrong? or, Do you have any advice on how to procede?