- #1
skrat
- 748
- 8
Homework Statement
A long squared rod with sides ##l##, thermal conductivity ##\lambda ##, specific heat ##c_p## and density ##\rho ## is cooled with water with constant temperature. Inside our squared rod we have a long square heater with sides ##l/2## and is completely centered inside the rod. Heat sources per unit volume in this is equal to ##q [W/m^3]##. (When calculating the temperature, you can use only the biggest three terms of the sum)
What is the maximum possible temperature of the surrounding water , before the heater starts to melt down? Melting temperature is ##T_c##.
Homework Equations
The Attempt at a Solution
I am really not sure this is correct, but at least it is something:
So we have two squares, one with sides ##l## and the other one centered in the first one with sides ##l/2##. I have put the origin of my coordinate system in the left bottom corner of the bigger square. This may not be a really good idea but there is a good reason why I chose this option.
And the reason is, that if the origin is there, than the function that solves ##\nabla ^2T(x,y)=0## is $$T(x,y)=\sum_{n,m}A_{n,m}sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)$$ Of course the equation I am solving is not ##\nabla ^2T=0## instead I have to consider the heater! Therefore ##\nabla ^2 T(x,y)=-\frac{q}{\rho c_p D}##.
Meaning I should write ##\frac{q}{\rho c_p D}## in terms of this sum ##T(x,y)=\sum_{n,m}C_{n,m}sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)##, so let's do that: $$\frac{q}{\rho c_p D}\int_{l/4}^{3l/4}\int _{l/4}^{3l/4}sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)dxdy=C_{n,m}(\frac l 2)^2$$ Which brings me to $$C=\frac{8q}{\rho c_p D\pi ^2}\frac{1}{mn}$$ but only for (both) odd ##m## and ##n##.
So finally $$\frac{q}{\rho c_p D}=\sum_{m,n}\frac{8q}{\rho c_p D\pi ^2}\frac{1}{mn}sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)$$
Now I can start solving ##\nabla ^2 T(x,y)=-\frac{q}{\rho c_p D}## $$\sum _{m,n}A_{m,n}\frac{\pi ^2}{l^2}(m^2+n^2)sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)=-\sum_{m,n}\frac{8q}{\rho c_p D\pi ^2}\frac{1}{mn}sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)$$ Which brings me to
$$A_{m,n}=\frac{8qL^2}{\rho c_pD\pi^4}\frac{1}{mn(m^2+n^2)}$$ Therefore the solution should be $$ T(x,y)=\sum_{\text{odd} m,n}\frac{8qL^2}{\rho c_pD\pi^4}\frac{1}{mn(m^2+n^2)}sin(\frac{n\pi}{l}x)sin(\frac{n\pi}{l}y)$$
Now to calculate the temperature I think I should do this, but I am really not so sure about this part: $$P=\lambda S\frac{\Delta T}{L/4}$$ where ##P=\int qdV## so $$\int qdV=\lambda S\frac{T_c-T_x}{L/4}$$ where I used notation ##T_x## for the temperature that I hope the problem is asking after...
Hmm. What do you think?