Rate of heat loss with two thermal conductivity values

[flying fish]

Strangely, this is for a Mat Sci Class. It is probably not as complicated as I am making it out to be, but the book has no suggestions on how to handle this problem.

Homework Statement

Calculate rate of heat loss through a fireclay refractory wall of wall with thickness .1m. Furnace operates at 1000C and outside is at 100C. Thermal conductivity of the clay at 1000C is 1.5, and 1.1 at 100C

Homework Equations

Rate of heat loss per square meter = q = -kdT/dx
dT = 900C
dx = .1m
k(1000) = 1.5
k(100) = 1.1
(q is what I am trying to find)

The Attempt at a Solution

After realizing that there are two thermal conductivity values, so I could not just plug in the numbers...I figured I would assume a linear relationship between K and Temperature, so that k = 1.1 + .4T/900 - 100(.4/900). But I'm not sure if that gets me anywhere as for solving the problem. Or...could I just take the average k value to be (1.1+1.5)/2 and use that? I'm not sure what the appropriate way to go about this problem is...

 

[Mapes]

I agree that a good first approximation would use the average thermal conductivity to give

Q=-\frac{k_1+k_2}{2}\frac{dT}{dx}

A better approximation can be found by using the exact governing equation. You can find this by doing an energy balance on a control volume, if you're familiar with that. You'd have

-k(T)\frac{dT}{dx}

going in and

-k(T)\frac{dT}{dx}+\frac{d}{dx}\left[-k(T)\frac{dT}{dx}\right]dx

coming out. (The second term is the first expression in a Taylor series expansion of the first term.) An energy balance gives

\frac{dk(T)}{dT}\left(\frac{dT}{dx}\right)^2+k\frac{d^2T}{dx^2}=0

which is the governing equation. You would likely have to solve this numerically. (If k is varies linearly with T, there might be an exact solution; I just don't know it offhand.)

 

[flying fish]

Thanks! I have a feeling that the professor is looking for the simple average solution, since this was part of a question from a Mat Sci test, on a chapter that just briefly touched on thermal properties of materials.

But it is nice to know some of the uses for these other things I am learning in some of my classes (like Taylor series and math processing software)!

 

[HugoR]

Solution

Hello,

Mapes, I would like to continue in developing a discrete equation to solve numerically. What is the next step and what would be the solution of this equation.

<br /> \frac{dk(T)}{dT}\left(\frac{dT}{dx}\right)^2+k\frac{d^2T}{dx^2}=0<br />

After reading around, one have to discretized the above equation, and identify the initial and limits conditions. Your help would be appreciated.

Thank you.
HR

 

[Mapes]

Exactly, you'll discretize as (if I did this right)

\frac{dk(T)}{dT}\left(\frac{T_\mathrm{R}-T_\mathrm{L}}{2\Delta x}\right)^2+k(T)\frac{T_\mathrm{R}+T_\mathrm{L}-2T}{(\Delta x)^2}=0

T=\frac{T_\mathrm{R}+T_\mathrm{L}}{2}+\frac{dk(T)}{dT}\frac{(T_\mathrm{R}-T_\mathrm{L})^2}{8k(T)}

Solve for T at each node and iterate until the solution converges.

I just did this for the case of k(T) varying linearly with T and interestingly enough, the exact solution (11,700 W) is the same as if you use the average value of k = 1.3. So flying fish really had the right idea earlier.

 

[HugoR]

Mapes,

Thank you for your inputs. I don't really understand what you mean by solving for T at each node and iterate until the solution converges.

Please tell me if I am wrong:

solving for T means,

<br /> T - \frac{T_\mathrm{R}+T_\mathrm{L}}{2}+\frac{dk(T)} {dT}\frac{(T_\mathrm{R}-T_\mathrm{L})^2}{8k(T)} = 0<br />

at each node means,
discretize the domain which is the wall thickness.
Let say
Tl | .T1 .T2 .T3 .T4 | Tr
Tl and Tr are known.

Here is my logic:
I can find an equation for k(T) varying linearly with T,

<br /> k(T) = 0.0004 T + 1.0556<br />

then,

<br /> \frac{dk(T)}{dT} = 00004<br />

I don't know what to do next to solve (using MATLAB),

<br /> Q=-k(T)\frac{dT}{dx}<br />

it look straight forward but I am stock.

Thank you if you can help me more.
HR

 

[Mapes]

You're totally on the right track. You have several simultaneous equations (one for each node):

T_1 - \frac{T_\mathrm{2}+1000^\circ\mathrm{C}}{2}+\frac{dk} {dT}\frac{(T_\mathrm{2}-1000^\circ\mathrm{C})^2}{8k(T_1)} = 0

T_2 - \frac{T_\mathrm{3}+T_\mathrm{1}}{2}+\frac{dk} {dT}\frac{(T_\mathrm{3}-T_\mathrm{1})^2}{8k(T_2)} = 0

T_3 - \frac{T_\mathrm{4}+T_\mathrm{2}}{2}+\frac{dk} {dT}\frac{(T_\mathrm{4}-T_\mathrm{2})^2}{8k(T_3)} = 0

T_4 - \frac{100^\circ\mathrm{C}+T_\mathrm{3}}{2}+\frac{dk} {dT}\frac{(100^\circ\mathrm{C}-T_\mathrm{3})^2}{8k(T_4)} = 0

So you see, the left and right temperatures are relative to the node, not absolute. You can have MATLAB solve these equations by making a matrix of the coefficients and inverting it, but you'll need to iterate because k is a function of T.

Your equation for Q is exactly right, and you need to discretize it too. You should get the same Q between any two nodes (this is a good check on your results).

 

[HugoR]

Mapes,

It is clearer for me. Thank you for your precise and valuable inputs.

I solved the equations using MATLAB (for fun) and I found:

1 node, q = 11,700 W (2 iterations)
11 nodes, q = 10,200 W (115 iterations)
101 nodes, q = 9,935 W (6124 iterations)
1001 nodes, q = 9,901 W (125713 iterations)
10001 nodes, q = 9,899 W (too much iterations)
100001 nodes, q = 9,898 W (way too much interations even with less precision, residue = 0.01)