Heat flux equation coding

Hi

Can you please help me code a 1-D heat flux equation q =- k(dT/dR).

I have a tabulated data of both T and R values.

Thank you in advance

Answer 1

[Walter Roberson]

[sorted_t, idx] = sort(T);
sorted_R = R(idx);
gr = gradient(sorted_t, sorted_R);
q = -mean(gr);

[Torsten]

I think T is a function of R. Thus sorting T with respect to R seems more plausible in this case. And why do you take the mean ?

[Walter Roberson]

@Torsten

Ah, you are right, T being a function of R makes more sense for dT/dR . I was reading T as Time

[sorted_R, idx] = sort(R);
sorted_T = T(idx);
gr = gradient(sorted_R, sorted_T);
k = -mean(gr);

The user has a table of T and R values, so the implied task is to estimate k. Once the gradients have been estimated, the task turns into a best-fit of a linear multiplier with a known set of values, and the best fit in such cases is mean.

The situation would have been different if the user had a k value already, but they don't mention that.

[Nour]

I do have the value of K (thermal conductivity). You're right, T is temperature not time and it's function of R (radius).

Thank you. I appreciate it

[Torsten]

So problem solved ?