I Solution of the 1D heat equation

3
0
$$\frac{\partial T}{\partial t}=\alpha\frac{\partial^2 T}{\partial^2 t}$$

with an initial condition and boundary conditions
$$T(x,0)=T_0$$
$$T(L,t)=T_0$$
$$-k\left.\frac{\partial T}{\partial x}\right|_{x=0}=2A\cos^2\left(\frac{\omega t}{2}\right)=A(\cos\omega t+1)$$

where $A=V_0^2/(8RhL)$, $V_0$ is the voltage applied to the heater, R the electrical resistance of the heater, h the thickness of the thin film, $\alpha$ the thermal diffusivity of the thin film, and $\omega/2$ the heating frequency. The solution for this problem is

$$T(x,t)-T_0=\frac Ak\sqrt{\frac\alpha\omega}\exp\left (-\sqrt{\frac{\omega}{2\alpha}}x\right)\\ \times\cos\left(\omega t-\sqrt{\frac{\omega}{2\alpha}}x-\frac\pi4\right)-\frac Ak(x-L)$$


I got this from a paper. I'm trying to derive how the author came to the solution from the boundary conditions. There is no derivation in the paper, and I searched books and the internet thoroughly but couldn't find anything. I attached the solution I'm getting. But I don't think it's correct! Any help regarding where I'm doing wrong will be greatly appreciated.
 

Attachments

Last edited:

Orodruin

Staff Emeritus
Science Advisor
Homework Helper
Insights Author
Gold Member
2018 Award
15,850
5,847
That solution does not satisfy the boundary condition at x=L.
 
3
0
That solution does not satisfy the boundary condition at x=L.
But it was published in a good journal by a very renowned professor (now). They must have some logic behind it which is unfortunately not mentioned in the paper.
 

berkeman

Mentor
55,557
5,659
But it was published in a good journal by a very renowned professor (now).
Can you provide a link to the paper in the journal? Your PDF attachment does not include the beginning of the paper with that information. Thanks. :smile:
 
19,286
3,822
Like I told you in the other forum that you and I have been interacting in, the authors of this article implicitly assume that the oscillatory part of the solution damps out before it reaches x = L. The condition for this to happen is ##\sqrt{\frac{\omega}{2\alpha}}L>>0##.

In your analysis, Eqn. 17 is incorrect. It should read (based on the above assumption), $$a(x)=C_1\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\cos{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}+C_2\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\sin{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}$$
You are missing factors of 2 in the denominators of the square root terms in your expression.

To obtain b(x), you substitute this corrected form of Eqn. 17 into your Eqn. 14. What do you obtain?
 
  • Like
Reactions: BvU
19,286
3,822
Maybe it would help also if I told you that the four roots of of ##x^4=-1## are
$$x=\frac{1+i}{\sqrt{2}}$$
$$x=\frac{-1+i}{\sqrt{2}}$$
$$x=\frac{-1-i}{\sqrt{2}}$$
and $$x=\frac{1-i}{\sqrt{2}}$$
This is where the ##\sqrt{2}## comes from in the correct form of Eqn. 17.
 
19,286
3,822
In a previous post, I indicated that the solution for a(x) is $$a(x)=C_1\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\cos{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}+C_2\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\sin{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}$$
If I substitute this into your Eqn. 14, I obtain: $$b(x)=C_1\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\sin{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}-C_2\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\cos{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}$$But, applying the boundary condition on b(x) at x = 0, we have that ##C_2=0##. This means that $$a(x)=C_1\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\cos{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}$$ and $$b(x)=C_1\exp{\left(-\sqrt{\frac{\omega}{2\alpha}}x\right)}\sin{\left(\sqrt{\frac{\omega}{2\alpha}}x\right)}$$
In your Eqn. 11, you have the wrong sign on A/k. You gotta do the math correctly.
 

Want to reply to this thread?

"Solution of the 1D heat equation" You must log in or register to reply here.

Related Threads for: Solution of the 1D heat equation

Replies
2
Views
5K
Replies
20
Views
2K
Replies
6
Views
5K
Replies
7
Views
3K
  • Posted
Replies
1
Views
2K
  • Posted
Replies
4
Views
3K
Replies
1
Views
1K

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving

Hot Threads

Top