Hello! I am trying to solve the 1D Fourier's equation with temperature dependent thermal properties but I do not know how to get the temperature results from pdepe at each time step and use them to evaluate the material properties for the next time step. Can you help me sort this out? Find below the example.
In this case, for example, the thermal diffisivity is temperature dependent accordint to the law:
diffr=3.34E-10*T^4 - 9.37E-07*T^3 + 9.38E-04*T^2 - 4.80E-01*T + 2.09E+02
function [c,f,s] = pdex1pde(x,t,u,dudx,diffr,Tfront)
c = 1/diffr;
f = dudx;
s = 0;
Many thanks for the help!