Steady state heat equation in a rectangle with a punkt heat source

fluidistic
Gold Member
Messages
3,928
Reaction score
272
Homework Statement
Solve the steady state heat equation in a rectangle whose bottom surface is kept at a fixed temperature, left and right sides are insulated and top side too, except for a point in a corner where heat is generated constantly through time.
Relevant Equations
##\kappa \nabla ^2 T + g =0##
I have checked several textbooks about the heat equation in a rectangle and I have found none that deals with my exact problem. I have though to use separation of variables first (to no avail), then Green's function (to no avail), then simplifying the problem for example by defining a new function in terms of ##T(x,y)## such that it would satisfy a homogeneous problem instead, but to no avail. (is a problem even called homogeneous when ##dT/dx|_{x_0} = 0## rather than ##T(x=x_0)=0##? I guess not.)

Out of memory, when I went with separation of variables to tackle ##\kappa \left( \frac{\partial ^2 T}{\partial x^2}+ \frac{\partial ^2 T}{\partial y^2}\right) = 0##, I obtained solutions of the form ##X(x)Y(y)## with ##X(x)=A\cosh(\alpha x)+B\sinh(\alpha x)## and ##Y(y)=C\cos(\alpha y)+D\sin(\alpha y)## where ##\alpha## is the separation constant. The boundary conditions are of the type Dirichlet for the bottom surface: ##T(x,y=0)=T_0##. And Neumann elsewhere: ##\frac{\partial T}{\partial x}|_{x=0, y=0}## for ##y\in [0,b)##, ##\frac{\partial T}{\partial x}|_{x=a, y=0}## for ##y\in [0,b]## and ##\frac{\partial T}{\partial y}|_{x, y=b}## for ##x\in (0,a]##. The power generated translates as the Neumann boundary condition ##\nabla T \cdot \hat n## and so ##\frac{\partial T}{\partial x}|_{x=0, y=b}+ \frac{\partial T}{\partial y}|_{x=0, y=b}=p## where ##p## is the power density of the heat source.

I have been stuck there, I could not get to apply and know the constants ##A##, ##B##, ##C## and ##D##, nor ##\alpha##. All of these constants are in fact depending on ##n##, natural numbers, because the separable solutions are eigenfunctions, etc.

Any pointer would be appreciated. Thank you!
 
Physics news on Phys.org
I’m thinking of an array of sources and sinks.
 
Prove $$\int\limits_0^{\sqrt2/4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx = \frac{\pi^2}{8}.$$ Let $$I = \int\limits_0^{\sqrt 2 / 4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx. \tag{1}$$ The representation integral of ##\arcsin## is $$\arcsin u = \int\limits_{0}^{1} \frac{\mathrm dt}{\sqrt{1-t^2}}, \qquad 0 \leqslant u \leqslant 1.$$ Plugging identity above into ##(1)## with ##u...
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top