Solve the heat equation having Dirichlet boundary conditions

In summary, if you are solving for a non-eigen function initial condition, you will need to use the principle of superposition and use the appropriate boundary conditions.
  • #1
chwala
Gold Member
2,746
387
TL;DR Summary
See attached.
I am going through these notes...they are pretty easy to follow. I would like more insight on the initial condition. In this problem, (attachment below), i guess the choice of initial condition is convenient as its easier to plug in the values of ##n=2## and ##b=3## (highlighted on the attachment) to realize the required solution. My question is? supposing we had the initial condition being;

##f(x)=\dfrac{1}{2}\sin\left[\dfrac{4πx}{l}\right]## where then our ##n=4## and ##b_4=\dfrac{1}{2}##

then our solution would be;

##U(x,t)##= ##\dfrac{1}{2} e^{\left[\dfrac{-16π^2kt}{l^2}\right]} ×\sin\left[ \dfrac{4πx}{l}\right]## correct?

I think the other highlighted part is now clear to me i.e when the initial condition is not an eigen function, then we have to make use of the principle of superposition. (one thing that is clear to me here is that for us to have the eigen function then the initial condition has to be either a function of sine or cosine).

Can we have the initial condition being a non-trigonometric function which would give us an eigenfunction?

Cheers!

Find the problem and solution here from textbook;
1664017264211.png


1664017427007.png

1664017475893.png
 
Last edited:
Physics news on Phys.org
  • #2
The operator [itex]X'' + \lambda X[/itex] is self-adjoint with respect to the inner product [tex]
\langle f, g \rangle = \int_0^\ell f(x)g(x)\,dx,[/tex] so eigenfunctions corresponding to different eigenvalues will be orthogonal with respect to this inner product, and you can recover the [itex]b_n[/itex] from [tex]
\int_0^{\ell} f(x)X_m(x)\,dx = b_m \int_0^{\ell} X_m(x)^2\,dx[/tex] as long as you can do the integrals.
 
  • Like
Likes chwala
  • #3
pasmith said:
The operator [itex]X'' + \lambda X[/itex] is self-adjoint with respect to the inner product [tex]
\langle f, g \rangle = \int_0^\ell f(x)g(x)\,dx,[/tex] so eigenfunctions corresponding to different eigenvalues will be orthogonal with respect to this inner product, and you can recover the [itex]b_n[/itex] from [tex]
\int_0^{\ell} f(x)X_m(x)\,dx = b_m \int_0^{\ell} X_m(x)^2\,dx[/tex] as long as you can do the integrals.
Thanks @pasmith . Let me read more on this.
 
  • #4
chwala said:
Summary: See attached.

I am going through these notes...they are pretty easy to follow. I would like more insight on the initial condition. In this problem, (attachment below), i guess the choice of initial condition is convenient as its easier to plug in the values of ##n=2## and ##b=3## (highlighted on the attachment) to realize the required solution. My question is? supposing we had the initial condition being;

##f(x)=\dfrac{1}{2}\sin\left[\dfrac{4πx}{l}\right]## where then our ##n=4## and ##b_4=\dfrac{1}{2}##

then our solution would be;

##U(x,t)##= ##\dfrac{1}{2} e^{\left[\dfrac{-16π^2kt}{l^2}\right]} ×\sin\left[ \dfrac{4πx}{l}\right]## correct?

I think the other highlighted part is now clear to me i.e when the initial condition is not an eigen function, then we have to make use of the principle of superposition. (one thing that is clear to me here is that for us to have the eigen function then the initial condition has to be either a function of sine or cosine).

Can we have the initial condition being a non-trigonometric function which would give us an eigenfunction?

Cheers!

Find the problem and solution here from textbook;
View attachment 314544

View attachment 314545
View attachment 314546
I have been doing further reading on this by looking at heat diffusivity on insulated wire. I have noted that in general if the initial condition ##f(x)## is a sine series then the homogenous Dirichlet boundary conditions will apply for insulated wire (with non-insulated ends) resulting into a Fourier series of the odd periodic extension ##f(x)##. I have also noted that for an insulated wire (with insulated ends), the homogenous Neumann boundary conditions will apply and the ensuing solution would be Fourier series of the even peridic extension ##f(x)##...
 

FAQ: Solve the heat equation having Dirichlet boundary conditions

What is the heat equation?

The heat equation is a mathematical model that describes how heat is distributed and transferred in a given system. It is a partial differential equation that relates the rate of change of temperature with respect to time and space.

What are Dirichlet boundary conditions?

Dirichlet boundary conditions are a set of constraints that specify the values of a function at the boundaries of a given system. In the context of the heat equation, these conditions define the temperature at the boundaries of the system.

How do you solve the heat equation with Dirichlet boundary conditions?

To solve the heat equation with Dirichlet boundary conditions, you must first define the initial temperature distribution within the system and the boundary conditions. Then, you can use numerical methods such as finite difference or finite element methods to approximate the temperature distribution at different points in time.

What are some applications of the heat equation with Dirichlet boundary conditions?

The heat equation with Dirichlet boundary conditions has many applications in various fields, including physics, engineering, and environmental science. It is commonly used to model heat transfer in materials, heat flow in buildings, and temperature changes in the Earth's crust.

What are the limitations of the heat equation with Dirichlet boundary conditions?

The heat equation with Dirichlet boundary conditions assumes that the material being studied has homogeneous and isotropic properties, and that the boundary conditions remain constant over time. These assumptions may not hold true in real-world scenarios, leading to inaccuracies in the model's predictions.

Back
Top