Heat equation with non homogeneous BCs

jackkk_gatz
Messages
45
Reaction score
1
Homework Statement
$$\frac{1}{r}\frac{\partial }{\partial r}(r\frac{\partial T}{\partial r})+\frac{\partial^2T}{\partial z^2}=0$$

$$\left.-k\frac{\partial T}{\partial r}\right\rvert_{r=R}=h[T(R,z)-T_{\infty}]$$

$$\left.k\frac{\partial T}{\partial z}\right\rvert_{z=H}+h[T(r,H)-T_{\infty}]=q_s$$

$$\left.-k\frac{\partial T}{\partial z}\right\rvert_{z=0}=0$$

$$\left.-k\frac{\partial T}{\partial r}\right\rvert_{r=0}=0$$

$where \ T_0 \ ,T_{\infty} \ and \ q_s \ are \ constants$
Relevant Equations
-
I did a change of variable $$\theta(r,z) = T(r,z)-T_{\infty}$$ which resulted in

$$\frac{1}{r}\frac{\partial }{\partial r}(r\frac{\partial \theta}{\partial r})+\frac{\partial^2 \theta}{\partial z^2}=0$$

$$\left.-k\frac{\partial \theta}{\partial r}\right\rvert_{r=R}=h\theta$$

$$\left.k\frac{\partial \theta}{\partial z}\right\rvert_{z=H}+h\theta=q_s$$

$$\left.-k\frac{\partial \theta}{\partial z}\right\rvert_{z=0}=0$$

$$\left.-k\frac{\partial \theta}{\partial r}\right\rvert_{r=0}=0$$

After that I proposed $$\theta(r,z)=v(r,z)+w(r,z)$$ where w(r,z) should must satisfy

$$\left.-k\frac{\partial w}{\partial r}\right\rvert_{r=R}=hw(R,z)$$

$$\left.k\frac{\partial w}{\partial z}\right\rvert_{z=H}+hw(r,H)=q_s$$

$$\left.-k\frac{\partial w}{\partial z}\right\rvert_{z=0}=0$$

$$\left.-k\frac{\partial w}{\partial r}\right\rvert_{r=0}=0$$

I already tried interpolation, doesn't work. I don't know how w(r,z) should look like in order to satisfy the above equations. Is there an easier method?
In short, everything I have been trying has failed and I don't know what to do anymore, I have looked for books on PDEs, all the ones I have found deal with very simple cases, which are of no use to me. I have almost no experience solving this kind of equations to know what to do or to guess how w might look like
 
Last edited:
Physics news on Phys.org
I just noticed I asked my question in the wrong section 💀
 
Have you tried obtaining the general solution by the method of separation of variables? It looks like it will work in this case but I have not solved the problem.

What is ##x=0## doing in the boundary condition ##~\left.-k\dfrac{\partial T}{\partial r}\right\rvert_{x=0}=0~##. Is it a typo?

I reported this thread and it should be moved to the Advanced Homework forum by a mentor at some point in time.
 
kuruman said:
Have you tried obtaining the general solution by the method of separation of variables? It looks like it will work in this case but I have not solved the problem.

What is ##x=0## doing in the boundary condition ##~\left.-k\dfrac{\partial T}{\partial r}\right\rvert_{x=0}=0~##. Is it a typo?

I reported this thread and it should be moved to the Advanced Homework forum by a mentor at some point in time.
Yes it was a typo, fixed it already. And yes I have tried to get the the general solution by the method of separation of variables, the thing is I know how to apply it but with homogeneous BC where I do some things with Sturm-Liouville. The thing is Sturm-Liouville only works with homogeneous BCs, I know a method to transform the non homogeneous BCs to homogeneous, which is the one I already wrote where I have to guess the form of w(r,z)

And thanks for helping to move my question to the right section
 
jackkk_gatz said:
I have tried to get the the general solution by the method of separation of variables,
Wolfram gives me Bessel functions.
 
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top